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)

DL2 data

The DL2 files of the real data are stored here - /fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/Mrk501/v09x/tailcut84/DL2/real/src_indep/20220531_dec3476/{Date}

Link to MC DL2 files used: /fefs/aswg/data/mc/DL2/AllSky/20220531_dec3476_std/TestingDataset/dec_3476/

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

Runs included:

[2606 2607 2608 2610 4167 4168 4169 4452 4453 4454 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 5094 5095
5134 5202 5203 5204 5205 5206 5438 5439 5489 5490 5550 5551 6076 6077
6078 6863 6864 6865 6880 6881 6882 6928 6930 6931 6959 6978 6979 6980
6981 7155 7156 7158 7353 7354 7356 7357 7358 7359 7360 7361 7362 7363
7364 7381 7382 7383 7384 7385 7401 7402 7403 7482 7483 7485 7564 7565
7635 7636 7637 7718 7719 7720 7721 7722 7723 7724 7725 7775 7776 7777
7778 7860 7861 7862 7863 7875 7876 7877 8274 8275 8276 8325 8326 8381
8382 8383 8416 8417 8418 8419]

• 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 -
  • Log Parabola + EBL (Dominguez) (6.5 sigmas preference over Power Law + EBL)
• Reference energy used by calculating the decorrelation energy with a Power Law fit, 561.5 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

Spectral model: LogParabola

amplitude  4.44e-11 +/- 1.05e-12, alpha=2.27 +\- 0.02, beta=0.097 +\- 0.018, ref = 0.5615