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

General information

Name of the source: Mrk 421
Brief description of the source:
- Object type : AGN
- Redshift (z) : 0.0308
- HBL
- RA: 11 04 27.314 (hh mm ss), Dec: +38 12 31.8 (dd mm ss)
- RA, Dec in deg (ICRS): 166.113808, 38.20883287

Analysis by Chaitanya Priyadarshi (IFAE - cpriyadarshi@ifae.es)

Data-taking information

Dates of data-taking: 20201212, 202201213, 20201214, 20210212, 20210213, 20210214, 20210215, 20210216, 20210217, 20210219, 20210303, 20210310, 20210313, 20210316, 20210317, 20210318, 20210320, 20210321, 20210412, 20210416, 20210417, 20210418, 20210430, 20210505, 20210508, 20210509, 20210511, 20210531, 20210601, 20210606, 20210608

- 20201212 (3.07 h)
  - ZD range (deg): [11.2, 52.1]
  - AZ range (deg): [34.68, 64.37]
  - Runs: 3216 - 3227
  - Joint observations with MAGIC? : Not for 3219 and 3220

- 20201213 (0.7 h)
  - ZD range (deg): [10.16, 13.4]
  - AZ range (deg): [20.97, 42.6]
  - Runs: 3238, 3239
  - Joint observations with MAGIC? : Yes

- 20201214 (1.25 h)
  - ZD range (deg): [11.25, 22.41]
  - AZ range (deg): [31.59, 57.87]
  - Runs: 3247 - 3250
  - Joint observations with MAGIC? : Yes

- 20210212 (2.15 h)
  - ZD range (deg): [19.4, 45.78]
  - AZ range (deg): [-64.84, -55.85]
  - Runs: 3617 - 3624
  - Joint observations with MAGIC? : No

- 20210213 (1.07 h)
  - ZD range (deg): [27.45, 45.24]
  - AZ range (deg): [-64.58, -61.95]
  - Runs: 3641, 3642, 3644, 3645, 3646
  - Joint observations with MAGIC? : No

- 20210214 (2.42 h)
  - ZD range (deg): [14.38, 42.66]
  - AZ range (deg): [-64.31, -44.05]
  - Runs: 3661 - 3668
  - Joint observations with MAGIC? : No

- 20210215 (0.62 h)
  - ZD range (deg): [9.44, 10.42]
  - AZ range (deg): [-5.43, 24.23]
  - Runs: 3687, 3688
  - Joint observations with MAGIC? : No

- 20210216 (1.81 h)
  - ZD range (deg): [21.56, 48.6]
  - AZ range (deg): [-64.41, -58.23]
  - Runs: 3715 - 3721
  - Joint observations with MAGIC? : No

- 20210217 (1.59 h)
  - ZD range (deg): [9.11, 14.3]
  - AZ range (deg): [-45.62, 40.88]
  - Runs: 3732 - 3737
  - Joint observations with MAGIC? : No

- 20210219 (0.07 h)
  - ZD range (deg): [51.35]
  - AZ range (deg): [-62.46]
  - Runs: 3800
  - Joint observations with MAGIC? : No

- 20210313 (0.77 h)
  - ZD range (deg): [27.71, 37.26]
  - AZ range (deg): [-64.51, -62.07]
  - Runs: 4016 - 4019
  - Joint observations with MAGIC? : No

- 20210316 (1.26 h)
  - ZD range (deg): [12.55, 24.53]
  - AZ range (deg): [-61.46, -39.1]
  - Runs: 4097 - 4100
  - Joint observations with MAGIC? : No

- 20210317 (2.56 h)
  - ZD range (deg): [12.82, 44.27]
  - AZ range (deg): [-64.14, -40.23]
  - Runs: 4131 - 4138
  - Joint observations with MAGIC? : No

- 20210318 (1.28 h)
  - ZD range (deg): [33.49, 47.24]
  - AZ range (deg): [-63.86, -63.18]
  - Runs: 4153 - 4156
  - Joint observations with MAGIC? : No

- 20210320 (1.44 h)
  - ZD range (deg): [24.68, 43.9]
  - AZ range (deg): [-64.41, -60.53]
  - Runs: 4184 - 4188
  - Joint observations with MAGIC? : No

- 20210321 (0.31 h)
  - ZD range (deg): [40.09]
  - AZ range (deg): [-63.82]
  - Runs: 4208
  - Joint observations with MAGIC? : No

- 20210412 (0.55 h)
  - ZD range (deg): [52.74, 62.03]
  - AZ range (deg): [-62.18, -59.17]
  - Runs: 4441 - 4443
  - Joint observations with MAGIC? : Yes

- 20210417 (0.96 h)
  - ZD range (deg): [9.78, 12.51]
  - AZ range (deg): [-9.47, 38.73]
  - Runs: 4457 - 4459
  - Joint observations with MAGIC? : No

- 20210418 (0.31 h)
  - ZD range (deg): [32.75]
  - AZ range (deg): [-63.38]
  - Runs: 4479
  - Joint observations with MAGIC? : No

- 20210430 (1.59 h)
  - ZD range (deg): [11.94, 32.75]
  - AZ range (deg): [-62.11, -35.77]
  - Runs: 4568 - 4572
  - Joint observations with MAGIC? : No

- 20210505 (0.71 h)
  - ZD range (deg): [16.13, 23.12]
  - AZ range (deg): [-58.44, -50.43]
  - Runs: 4612 - 4614
  - Joint observations with MAGIC? : Yes

- 20210508 (0.27 h)
  - ZD range (deg): [14.12]
  - AZ range (deg): [-45.14]
  - Runs: 4649
  - Joint observations with MAGIC? : Yes

- 20210509 (1.59 h)
  - ZD range (deg): [15.27, 34.35]
  - AZ range (deg): [-63.6, -48.35]
  - Runs: 4667 - 4671
  - Joint observations with MAGIC? : No

- 20210531 (0.3 h)
  - ZD range (deg): [41.6, 45.69]
  - AZ range (deg): [-63.77, -63.39]
  - Runs: 4852, 4853
  - Joint observations with MAGIC? : No

- 20210601 (0.64 h)
  - ZD range (deg): [46.22, 53.68]
  - AZ range (deg): [-63.32, -61.95]
  - Runs: 4868, 4869
  - Joint observations with MAGIC? : No

- 20210606 (0.32 h)
  - ZD range (deg): [37.16]
  - AZ range (deg): [-64.47]
  - Runs: 4939
  - Joint observations with MAGIC? : No

- 20210608 (0.25 h)
  - ZD range (deg): [38.21]
  - AZ range (deg): [-63.78]
  - Runs: 4962
  - Joint observations with MAGIC? : Yes

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/20201212/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 27 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