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

From my_wiki
Revision as of 15:18, 15 June 2023 by Epons (talk | contribs) (gammapy analysis parameters)
Jump to: navigation, search

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
      • 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

1ES0647 202304 v0 countmap.png

Theta2 plot

1ES0647 202304 v0 theta2.png

SED

  • power-law + EBL = best fit

1ES0647 202304 v0 spectra power law+ebl 2.png 1ES0647 202304 v0 spectra power law+ebl model.png

  • log-parabola + EBL

1ES0647 202304 v0 spectra log parabola+ebl 2.png 1ES0647 202304 v0 spectra log parabola+ebl model.png