|
|
Line 1: |
Line 1: |
− | = 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
| |
− |
| |
− | * 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
| |
− |
| |
− | [[File:Gammaness_cuts_distribution_80%_efficiency_different_MC_pointing.png]] [[File:Theta_cuts_distribution_80%_efficiency_different_MC_pointing.png]]
| |
− |
| |
− | [[File:Effective_area_80%_efficiency_different_MC_pointing.png]] [[File: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-
| |
− |
| |
− | [[File:LC_Run_2177_6078_lp_model_e_ref_487GeV_fit_300GeV_100TeV.png]]
| |
− |
| |
− | [[File:SED_Run_2177_6078_lp_model_e_ref_487GeV_fit_range_w_5_bin_p_decade_10_40.png]]
| |