1ES0647+250 Analysis with lstchain v0.9.12 E. Pons

1ES 0647+250 Analysis with lstchain v0.9.12 E. Pons

General information

• Name of the source: 1ES 0647+250
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.41
  • HBL
  • RA: 6 50 46.49 (hh mm ss), Dec: +25 2 59.62 (dd mm ss)
  • RA, Dec in deg (ICRS): 102.69370933, 25.04989555

• Analysis by Estelle Pons (LAPP - estelle.pons@lapp.in2p3.fr)

Dataset (< 2022-05-24)

• Number of selected observations : 28
• Runs: 3281 - 3364
• Date range: 2020-12-16 to 2020-12-21

obsids:

[3281 3282 3283 3284 3285 3287 3288 3331 3332 3333 3334 3335 3341 3342
3343 3344 3345 3346 3347 3348 3357 3358 3359 3360 3361 3362 3363 3364]

DL1 data

DL1 files produced by LSTOSA (lstchain v0.9)

• original real DL1a data

/fefs/aswg/data/real/DL1/{Date}/v0.9/tailcut84/

DL1 to DL2

• Random Forest

Standard ‘un-tuned’ production for processing standard AGN data, with lstchain v0.9.12 standard settings (az instead of sin(az) in training features).

/fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_2276

• DL2 data

/fefs/aswg/workspace/estelle.pons/1ES0647/DL2/{Date}/v0.9.12/tailcut84

NB: The lstchain config file from the RF production is specified for lstchain_dl1_to_dl2

Example:

lstchain_dl1_to_dl2 
  -f /fefs/aswg/data/real//DL1/20201219/v0.9/tailcut84/dl1_LST-1.Run03348.0012.h5                 
  --output-dir /fefs/aswg/workspace/estelle.pons/1ES0647//DL2/20201219/v0.9.12/tailcut84/                 
  --path-models /fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_2276/ 
  --config /fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_2276/lstchain_config.json

DL2 to DL3

• DL3 cuts: source-independent analysis & global cuts:
  • gh_cut = 0.7
  • theta_cut = 0.2

• MC DL2 files used to produce IRF

/fefs/aswg/data/mc/DL2/AllSky/20230414_v0.9.13_dec_931_az_tel/TestingDataset/dec_931

The closest node compared to the telescope pointing position is then selected for each run to create the corresponding IRF.

Example:

lstchain_create_irf_files  
  -g /fefs/aswg/data/mc/DL2/AllSky/20230414_v0.9.13_dec_931_az_tel/TestingDataset/dec_931/node_theta_32.059_az_248.099_/dl2_20230414_v0.9.13_dec_931_az_tel_node_theta_32.059_az_248.099__merged.h5                 
  -o /fefs/aswg/workspace/estelle.pons/PG1553/IRF/allsky/irf_20230414_v0.9.13_dec_931_az_tel_dec_931_node_theta_32.059_az_248.099_ghcut0.7_run8787.fits.gz
  --point-like  --global-theta-cut 0.2 --global-gh-cut 0.7

• DL3 data

/fefs/aswg/workspace/estelle.pons/PG1553/DL3/lstchain/v0.9.13/point_like/offset0.4/ghcut0.7/{Date}

Example:

lstchain_create_dl3_file 
  -d /fefs/aswg/workspace/estelle.pons/PG1553//DL2/20220617/v0.9.13/tailcut84/dl2_LST-1.Run08787.h5     
  -o /fefs/aswg/workspace/estelle.pons/PG1553/DL3/lstchain/v0.9.13/point_like/offset0.4/ghcut0.7/20220617     
  --input-irf /fefs/aswg/workspace/estelle.pons/PG1553/IRF/allsky/irf_20230414_v0.9.13_dec_931_az_tel_dec_931_node_theta_32.059_az_248.099_ghcut0.7_run8787.fits.gz     
  --source-name PG1553+113     --source-ra " 238.92935002022deg"     --source-dec " 11.19010156943deg"  
  --global-gh-cut 0.7 --gzip

gammapy analysis parameters

• gammapy version v1.0.1
• DL4 datasets
  • on_radius = 0.2 deg / exclusion_radius = 0.35 deg
  • Reflected Regions Background method
  • Safe Energy mask:
    • MapDataset: methods=["offset-max"], offset_max=4 deg
    • SpectrumDataset: NO safe mask
  • Energy axes:
    • True: 10 GeV to 10 TeV with 5 bins per decade
    • Reco: 10 GeV to 40 TeV with 10 bins per decade
• DL5 datasets
  • Spectral analysis
    • Joint fit
    • Energy axes Flux Points: 10 GeV to 10 TeV with 5 bins per decade
    • Reference energy: pivot energy from a power-law fit = 0.12 TeV
    • SED fitting: best fit between power-law + EBL (Dominguez), log-parabol + EBL, exponantial cutoff power-law + EBL

Analysis Results

Livetime : 8.17 h

Significance : 7.33 sigma

Excess : 1711.1

S/B : 0.040

Rate : 3.49 per min

Count Map

Theta2 plot

SED

• power-law + EBL = best fit

• log-parabola + EBL