Mrk 501 Preliminary Analysis with lstchain v0.9.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, 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

Theta2 plot

For energy-dependent cuts

Significance map

Excess map

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:
Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Mrk_501_Preliminary_Analysis_with_lstchain_v0.9.x_Priyadarshi&oldid=4852"