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

From my_wiki
Jump to: navigation, search

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

Gammaness cuts distribution 80% efficiency different MC pointing.png Theta cuts distribution 80% efficiency different MC pointing.png

Effective area 80% efficiency different MC pointing.png Energy dispersion 80% efficiency different MC pointing.png

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

Theta2 plot

Significance map

Excess map

Spectral results

Preliminary LC and SED, for the full dataset-

LC Run 2177 6078 lp model e ref 487GeV fit 300GeV 100TeV.png LC Run 2177 6078 lp model e ref 487GeV fit 30GeV 100TeV.png

SED Run 2177 6078 lp model e ref 487GeV fit range w 5 bin p decade 10 40.png

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Mrk_501_Preliminary_Analysis_with_lstchain_v0.7.x_Priyadarshi&oldid=4641"