Mrk 501 Xcheck lstchain v0.9.x Heckmann

Xcheck following the analysis in: Mrk 501 Preliminary Analysis with lstchain v0.9.x Priyadarshi w IRF interpolation

General information

• Name of the source: Mrk 501
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.034
  • 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 Heckmann Lea (MPP - heckmann@mpp.mpg.de)

DL1 data and models

The DL1 files of the real data are stored here - /remote/ceph/user/h/heckmann/Mrk501_LST/data/DL1/v0.9.3/tailcut84/

Models used: /fefs/aswg/data/models/AllSky/20230127_v0.9.12_base_prod_az_tel/dec_3476/

DL2 data

The DL2 files of the real data are stored here - /fefs/aswg/workspace/lea.heckmann/Data/Mrk501/v09x/tailcut84/DL2/real/src_indep/20230127_dec3476/

Link to MC DL2 files used: /fefs/aswg/data/mc/DL2/AllSky/20230127_v0.9.12_base_prod/TestingDataset/dec_3476/

IRFs: /fefs/aswg/workspace/lst1-agn-paper/LST1/Mrk501/Analysis2/source-independent/irfs/irf_20230127_v0.9.12_base_prod_az_tel_node_theta*

DL3 data selection

For energy-dependent cuts

Information about your DL3 data selection.

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

• source-independent
  • gh_efficiency: 0.8
  • theta_containment: 0.8

High-level analysis

Data selection:

max_zen = 52 
min_livetime = 60s 
wobble_offset_diff = 0.1 (from 0.4)
max_ped_std = 1.8 phe
pixel rate of >10 and >30 phe pulses in cosmics

153 Runs included:

[2177 2178 2179 2606 2607 2608 2610 4139 4140 4167 4168 4169 4170 4191
4192 4193 4194 4195 4212 4213 4214 4215 4452 4455 4473 4474 4475 4476
4480 4483 4484 4496 4497 4498 4586 4587 4588 4597 4598 4599 4600 4601
4602 4625 4626 4627 4628 4642 4643 4644 4645 4656 4657 4672 4673 4674
4675 4676 4677 4922 4933 4934 4944 4945 4946 4947 5011 5037 5038 5039
5040 5094 5095 5134 5202 5203 5204 5205 5206 5227 5228 5438 5439 5489
5490 5550 5551 6076 6077 6078 6862 6863 6864 6865 6879 6880 6881 6882
6928 6929 6930 6931 6958 6959 6978 6979 6980 6981 7155 7156 7158 7351
7353 7354 7355 7356 7357 7358 7359 7360 7361 7362 7363 7364 7381 7382
7383 7384 7385 7401 7402 7403 7482 7483 7564 7565 7635 7636 7637 7718
7719 7720 7721 7722 7723 7724 7775 7777 7778 7875 7876 8276 8416] 

• lstchain Tools to generate source-independent IRF and DL3
• Science Tool: gammapy 1.0
• point-like IRF, 1D analysis

For DL3 data reduction to DL4 dataset,

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

For spectral analysis,

• Model used -
  • Power law or Log Parabola + EBL (Dominguez) - depending on if log parabola is preferred at >3sigma
• Reference energy used by calculating the decorrelation energy with a Power Law fit, 574 GeV (Full dataset)
• Spectral Fit energy range - 10 GeV to 100 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

For LC,

• E_min = 100 Gev

Analysis Results

Preliminary LC and SEDs in comparison

• Full Dataset