Mrk 501 Preliminary Analysis with lstchain v0.7.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, 20200713, 20200820, 20200909, 20210310, 20210317, 20210319, 20210320, 20210321, 20210416, 20210417, 20210418, 20210419, 20210503, 20210504, 20210506, 20210507, 20210508, 20210509, 20210605, 20210606, 20210608, 20210612, 20210613, 20210702, 20210708, 20210710, 20210716, 20210801, 20210805, 20210808, 20210906

• • 20200710 (0.97 h)
    • ZD range (deg): [11.19, 16.25]
    • AZ range (deg): [-45.85, -9.97]
    • Runs: 2177 - 2180
    • Joint observations with MAGIC? : No

• • 20200909 (0.73 h)
    • ZD range (deg): [38.48, 61.18]
    • AZ range (deg): [-61.24, -57.49]
    • Runs: 2606, 2607, 2608, 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.01 h)
    • ZD range (deg): [12.52, 19.98]
    • AZ range (deg): [26.89, 51.3]
    • Runs: 4191 - 4194
    • 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.77 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.56 h)
    • ZD range (deg): [27.25, 35.9]
    • AZ range (deg): [58.21, 61.65]
    • Runs: 4480, 4483
    • 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 (1.22 h)
    • ZD range (deg): [23.48, 32.04]
    • AZ range (deg): [-60.92, 59.2]
    • Runs: 4921, 4922, 4933, 4934
    • Joint observations with MAGIC? : Not for 4921 and 4922

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

• • 20210608 (0.43 h)
    • ZD range (deg): [29, 35.09]
    • AZ range (deg): [-60.86, -59.07]
    • Runs: 4977, 4979
    • Joint observations with MAGIC? : Yes

• • 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.62 h)
    • ZD range (deg): [35.39, 43.85]
    • AZ range (deg): [-61.82, -60.71]
    • Runs: 5094 - 5097
    • Joint observations with MAGIC? : No

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

• • 20210716 (0.09 h)
    • ZD range (deg): [13.42]
    • AZ range (deg): [-32.76]
    • Runs: 5277
    • Joint observations with MAGIC? : No

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

Monte Carlo information

• Link to MC DL2 files used:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/ /fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/

• • Particle types: point-gamma (offset=0.4 deg), proton, electron
  • ZD (deg): [20, 40]
  • AZ (deg): [180, 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.7.1 and v0.7.3)

• original real DL1a data

/fefs/aswg/data/real/DL1/{Date}/v0.7.{1,3}/tailcut84/

• lstchain v0.7.1 tailcut8-4 (with cleaning based on pedestal RMS)

 /fefs/aswg/data/real/DL1/20210319/v0.7.1/log/lstchain_standard_config_tailcut84.json

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

 /fefs/aswg/data/real/DL1/20210320/v0.7.3/tailcut84/log/lstchain_standard_config_tailcut84.json

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

 /fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/lstchain_standard_config_v075_Tel1_tail84_dynamic_clean.json

• MC (zen 40 deg): tailcut8-4 (with cleaning based on pedestal RMS)

/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/lstchain_standard_config_v073_Tel1_tail_8_4_MC_config.json

Random forest

• lstchain v0.7.5 source-indep, zenith 20 deg, azimuth 180 deg, dynamic cleaning

 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/lstchain_standard_config_v075_Tel1_tail84_dynamic_clean.json

• lstchain v0.7.3 source-indep, zenith 40 deg, azimuth 180 deg

 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/lstchain_standard_config_v073_Tel1_tail_8_4_MC_config.json

DL2 data

Similar to DL1 files

DL3 data selection

Information about your DL3 data selection.

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

• source-indep, energy-dependent cuts
  • 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 24.5 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
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~490 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

Preliminary LC and SED, for the full dataset-