Mrk 501 Preliminary Analysis with lstchain v0.9.x Priyadarshi

Mrk 501 Preliminary Analysis with lstchain v0.7.x Priyadarshi

General information

• Name of the source: Mrk 501
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.03364
  • EHBL
  • RA: 16 53 52.216 (hh mm ss), Dec: +39 45 36.613 (dd mm ss)
  • RA, Dec in deg (ICRS): 253.46756699, 39.76017034

• Analysis by Chaitanya Priyadarshi (IFAE - cpriyadarshi@ifae.es)

Data-taking information

• Dates of data-taking: 20200710, 20200909, 20210317, 20210319, 20210320, 20210321, 20210416, 20210417, 20210418, 20210419, 20210503, 20210504, 20210506, 20210507, 20210508, 20210509, 20210605, 20210606, 20210612, 20210613, 20210702, 20210705, 20210708, 20210710, 20210801, 20210805, 20210808, 20210906, 20220203, 20220204, 20220205, 20220206, 20220207, 20220228, 20220308, 20220309, 20220310

• • 20200710 (0.47 h)
    • ZD range (deg): [11.19, 14.53]
    • AZ range (deg): [-37.81, -9.97]
    • Runs: 2177, 2179
    • Joint observations with MAGIC? : No

• • 20200909 (0.82 h)
    • ZD range (deg): [38.48, 61.18]
    • AZ range (deg): [-61.24, -57.49]
    • Runs: 2606, 2607, 2608, 2610, 2612, 2613
    • Joint observations with MAGIC? : No

• • 20210317 (0.64 h)
    • ZD range (deg): [18.7, 23.5]
    • AZ range (deg): [49.17, 55.47]
    • Runs: 4139, 4140
    • Joint observations with MAGIC? : No

• • 20210319 (0.91 h)
    • ZD range (deg): [32.56, 42.27]
    • AZ range (deg): [60.24, 61.3]
    • Runs: 4167 - 4170
    • Joint observations with MAGIC? : Yes

• • 20210320 (1.05 h)
    • ZD range (deg): [11.42, 19.98]
    • AZ range (deg): [15.32, 51.3]
    • Runs: 4191 - 4195
    • Joint observations with MAGIC? : No

• • 20210321 (0.95 h)
    • ZD range (deg): [13.28, 21.8]
    • AZ range (deg): [34.58, 53.74]
    • Runs: 4212 - 4215
    • Joint observations with MAGIC? : No

• • 20210416 (0.47 h)
    • ZD range (deg): [10.8, 15.46]
    • AZ range (deg): [-39.5, -10.35]
    • Runs: 4452, 4455
    • Joint observations with MAGIC? : No

• • 20210417 (0.95 h)
    • ZD range (deg): [10.98, 17.1]
    • AZ range (deg): [-44.23, -1.53]
    • Runs: 4473 - 4476
    • Joint observations with MAGIC? : Yes

• • 20210418 (0.87 h)
    • ZD range (deg): [24.1, 35.9]
    • AZ range (deg): [54.96, 61.65]
    • Runs: 4480, 4483, 4484
    • Joint observations with MAGIC? : Yes

• • 20210419 (0.85 h)
    • ZD range (deg): [16.05, 22.83]
    • AZ range (deg): [43.04, 53.73]
    • Runs: 4496 - 4498
    • Joint observations with MAGIC? : No

• • 20210503 (1.13 h)
    • ZD range (deg): [10.62, 19.4]
    • AZ range (deg): [3.84, 50.41]
    • Runs: 4586 - 4590
    • Joint observations with MAGIC? : Yes

• • 20210504 (1.76 h)
    • ZD range (deg): [20.28, 41.63]
    • AZ range (deg): [51.78, 61.31]
    • Runs: 4597 - 4602
    • Joint observations with MAGIC? : No

• • 20210506 (1.26 h)
    • ZD range (deg): [13.52, 25.26]
    • AZ range (deg): [35.58, 56.91]
    • Runs: 4625 - 4628
    • Joint observations with MAGIC? : No

• • 20210507 (0.81 h)
    • ZD range (deg): [17.19, 25.34]
    • AZ range (deg): [-57.98, -46.16]
    • Runs: 4642 - 4645
    • Joint observations with MAGIC? : Yes

• • 20210508 (0.47 h)
    • ZD range (deg): [20.82, 24.77]
    • AZ range (deg): [52.53, 56.58]
    • Runs: 4656, 4657
    • Joint observations with MAGIC? : No

• • 20210509 (1.9 h)
    • ZD range (deg): [13.64, 34.37]
    • AZ range (deg): [36.21, 60.13]
    • Runs: 4672 - 4677
    • Joint observations with MAGIC? : No

• • 20210605 (0.58 h)
    • ZD range (deg): [28.46, 32.04]
    • AZ range (deg): [-60.92, -57.97]
    • Runs: 4933, 4934
    • Joint observations with MAGIC? : Yes

• • 20210606 (1.28 h)
    • ZD range (deg): [11.03, 14.28]
    • AZ range (deg): [-39, 26.44]
    • Runs: 4944 - 4947
    • Joint observations with MAGIC? : No

• • 20210612 (0.32 h)
    • ZD range (deg): [26.47]
    • AZ range (deg): [57.72]
    • Runs: 5011
    • Joint observations with MAGIC? : No

• • 20210613 (1.08 h)
    • ZD range (deg): [23.78, 37.05]
    • AZ range (deg): [55.73, 61.12]
    • Runs: 5036 - 5040
    • Joint observations with MAGIC? : Not for 5036

• • 20210702 (0.51 h)
    • ZD range (deg): [35.39, 41.29]
    • AZ range (deg): [-61.26, -60.71]
    • Runs: 5094 - 5096
    • Joint observations with MAGIC? : No

• • 20210705 (0.08 h)
    • ZD range (deg): [11.49]
    • AZ range (deg): [16.3]
    • Runs: 5134
    • Joint observations with MAGIC? : NA

• • 20210708 (1.38 h)
    • ZD range (deg): [11.13, 16.05]
    • AZ range (deg): [-43.17, 35.63]
    • Runs: 5202 - 5206
    • Joint observations with MAGIC? : No

• • 20210710 (0.64 h)
    • ZD range (deg): [19.64, 23.63]
    • AZ range (deg): [-56.65, -49.36]
    • Runs: 5227, 5228
    • Joint observations with MAGIC? : Yes

• • 20210801 (0.8 h)
    • ZD range (deg): [29.09, 34.2]
    • AZ range (deg): [-60.67, -59.89]
    • Runs: 5438, 5439
    • Joint observations with MAGIC? : Not for 5438

• • 20210805 (0.57 h)
    • ZD range (deg): [26, 29.64]
    • AZ range (deg): [-58.5, -57.43]
    • Runs: 5489, 5490
    • Joint observations with MAGIC? : Yes

• • 20210808 (0.53 h)
    • ZD range (deg): [29.9, 33.46]
    • AZ range (deg): [-60.52, -60.26]
    • Runs: 5550, 5551
    • Joint observations with MAGIC? : Yes

• • 20210906 (0.64 h)
    • ZD range (deg): [27.02, 33.13]
    • AZ range (deg): [-60.45, -57.15]
    • Runs: 6076 - 6078
    • Joint observations with MAGIC? : Not for 6078

• • 20220203 (2.21 h)
    • ZD range (deg): [34.39, 60.0]
    • AZ range (deg): [58.33, 61.82]
    • Runs: 6858-6865
    • Joint observations with MAGIC? : NA

• • 20220204 (2.06 h)
    • ZD range (deg): [39.12, 63.33]
    • AZ range (deg): [57.3, 61.32]
    • Runs: 6876-6882
    • Joint observations with MAGIC? : NA

• • 20220205 (1.28 h)
    • ZD range (deg): [32.78, 48.04]
    • AZ range (deg): [60.38, 61.69]
    • Runs: 6927-6931
    • Joint observations with MAGIC? : NA

• • 20220206 (1.94 h)
    • ZD range (deg): [34.81, 62.09]
    • AZ range (deg): [57.7, 61.79]
    • Runs: 6953-6959
    • Joint observations with MAGIC? : NA

• • 20220207 (1.96 h)
    • ZD range (deg): [30.79, 54.04]
    • AZ range (deg): [59.01, 61.92]
    • Runs: 6975-6981
    • Joint observations with MAGIC? : NA

• • 20220228 (0.49 h)
    • ZD range (deg): [19.45, 22.8]
    • AZ range (deg): [51.92, 53.74]
    • Runs: 7157, 7158
    • Joint observations with MAGIC? : NA

• • 20220308 (1.88 h)
    • ZD range (deg): [14.62, 51.92]
    • AZ range (deg): [40.53, 61.73]
    • Runs: 7351, 7352, 7355, 7356, 7359, 7360, 7363, 7364
    • Joint observations with MAGIC? : NA

• • 20220309 (0.42 h)
    • ZD range (deg): [26.06, 28.97]
    • AZ range (deg): [58.23, 58.45]
    • Runs: 7382, 7383
    • Joint observations with MAGIC? : NA

• • 20220310 (0.59 h)
    • ZD range (deg): [14.14, 17.58]
    • AZ range (deg): [38.38, 45.44]
    • Runs: 7402, 7403
    • Joint observations with MAGIC? : NA

Monte Carlo information

For energy-dependent cuts

• Link to MC DL2 files used:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf

• • Particle types: point-gamma (offset=0.4 deg), proton, electron
  • ZD (deg): [20]
  • AZ (deg): [180]

DL1 data

This includes whether you use LSTOSA, specific versions of lstchain, dllab scripts, cleaning levels, and calibration information.

DL1a files produced by LSTOSA (lstchain v0.9)

• original real DL1a data

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

• lstchain v0.9 tailcut8-4 (with cleaning based on pedestal RMS, dynamical cleaning)

 /fefs/aswg/data/real/DL1/20220228/v0.9/tailcut84/log/lstchain_config_tailcut84_v092.json

• MC (zen 20 deg): tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF tuning

 /fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf/lstchain_config.json

Random forest

• lstchain v0.9 source-indep, zenith 20 deg, azimuth 180 deg, dynamic cleaning, PSF tuning

 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf/lstchain_config.json

DL2 data

/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/Mrk501/v09x/tailcut84/DL2/real/src_indep/{Date} The DL2 files were trained using RF model, with zenith 20 deg

DL3 data selection

For energy-dependent cuts

Information about your DL3 data selection.

• intensity > 0
• r: [0, 1]
• wl: [0.01, 1]
• leakage_intensity_width_2: [0, 1]

• source-indep
  • gh_efficiency: 0.8
  • theta_containment: 0.8

High-level analysis

• lstchain Tools to generate source-independent IRF and DL3
• Science Tool: gammapy 0.20.dev309+gb11005a (Developer build to extract signal from energy-dependent regions)
• point-like IRF, 1D analysis

A total of 36.58 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

• OFF regions selected - 1 Wobble
• Safe Energy mask for systematics - No
• Energy axes -
  • True: 10 GeV to 100 TeV with 10 bins per decade
  • Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

• Model used - Log Parabola, Power Law
• Reference energy used by calculating the decorrelation energy with a Power Law fit, 450 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Analysis Results

Spectral results

For energy-dependent cuts

Preliminary LC and SED, for the full dataset-

Checking SEDs for different time ranges based on the LC,

• Preliminary SED, for Runs 4586 - 4947:

• Preliminary SED, for Runs 7157 - 7403: