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

Mrk 421 Preliminary Analysis with lstchain v0.9.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, 20210214, 20210215, 20210216, 20210217, 20210313, 20210316, 20210317, 20210318, 20210320, 20210321, 20210416, 20210417, 20210418, 20210430, 20210505, 20210508, 20210509, 20210531, 20210601, 20210606, 20220205, 20220206, 20220208, 20220301, 20220305, 20220309

• • 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

• • 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

• • 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 (0.31 h)
    • ZD range (deg): [33.49]
    • AZ range (deg): [-63.52]
    • Runs: 4153
    • 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

• • 20210416 (0.32 h)
    • ZD range (deg): [11.24]
    • AZ range (deg): [-34.66]
    • Runs: 4441
    • Joint observations with MAGIC? : Yes

• • 20210417 (0.32 h)
    • ZD range (deg): [9.96]
    • AZ range (deg): [-9.47]
    • Runs: 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

• • 20220205 (1.26 h)
    • ZD range (deg): [9.23, 20.14]
    • AZ range (deg): [11.92, 56.75]
    • Runs: 6900-6903
    • Joint observations with MAGIC? : NA

• • 20220206 (2.13 h)
    • ZD range (deg): [9.01, 43.2]
    • AZ range (deg): [-28.74, 63.86]
    • Runs: 6940, 6941, 6947-6952
    • Joint observations with MAGIC? : NA

• • 20220208 (1.0 h)
    • ZD range (deg): [9.56, 19.05]
    • AZ range (deg): [-9.97, 55.47]
    • Runs: 6994, 6997-6999
    • Joint observations with MAGIC? : NA

• • 20220301 (0.69 h)
    • ZD range (deg): [29.78, 36.38]
    • AZ range (deg): [-63.56, -62.75]
    • Runs: 7182-7184
    • Joint observations with MAGIC? : NA

• • 20220305 (0.2 h)
    • ZD range (deg): [32.65, 34.99]
    • AZ range (deg): [-63.71, -63.40]
    • Runs: 7289, 7290
    • Joint observations with MAGIC? : NA

• • 20220309 (0.11 h)
    • ZD range (deg): [9.14]
    • AZ range (deg): [-11.78]
    • Runs: 7370
    • 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/20220208/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/Mrk421/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 29.9 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, 293 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Analysis Results

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 3216 - 3250:

• Preliminary SED, for Runs 3661 - 3737:

• Preliminary SED, for Runs 4016 - 4441:

• Preliminary SED, for Runs 4459 - 4671:

• Preliminary SED, for Runs 4852 - 7370: