Difference between revisions of "M87 TTHA"
Tora.arnesen (talk | contribs) (→M87 analysis (TTHA)) |
Tora.arnesen (talk | contribs) (→High-level analysis) |
||
(16 intermediate revisions by the same user not shown) | |||
Line 13: | Line 13: | ||
* Analysis by Tora T. H. Arnesen (Instituto de Astrofísica de Canarias - <tora.arnesen@iac.es>) | * Analysis by Tora T. H. Arnesen (Instituto de Astrofísica de Canarias - <tora.arnesen@iac.es>) | ||
− | ==Data | + | ==Data Selection== |
Data selected with data quality notebook | Data selected with data quality notebook | ||
Line 41: | Line 41: | ||
** min_fraction_around_mode = 0.8 | ** min_fraction_around_mode = 0.8 | ||
** max_intensity_at_half_peak_rate = 70 | ** max_intensity_at_half_peak_rate = 70 | ||
+ | |||
+ | |||
+ | |||
+ | [[File:Parameters_M87_240609_selection Large.png|720px]] | ||
+ | |||
+ | [[File:RelLightYield_M87_240609_selection.png|720px]] | ||
+ | |||
+ | [[File:Nsubruns_M87_240609_selection.png|720px]] | ||
+ | |||
+ | [[File:CosmicsTriggerRate_M87_240609_selection.png|720px]] | ||
+ | |||
+ | [[File:M87_240609_SelectionRuns.png|720px]] | ||
+ | |||
===Selected Data=== | ===Selected Data=== | ||
Line 87: | Line 100: | ||
20240604 : [17641 17642 17643 17644] | 20240604 : [17641 17642 17643 17644] | ||
20240607 : [17676 17677 17678 17679 17680 17681] | 20240607 : [17676 17677 17678 17679 17680 17681] | ||
− | <nowiki> | + | </nowiki> |
− | |||
==Monte Carlo information== | ==Monte Carlo information== | ||
Line 146: | Line 158: | ||
* Gammapy 1.2 (1.1 for Theta2 and Skymaps) | * Gammapy 1.2 (1.1 for Theta2 and Skymaps) | ||
* 1D analysis (point-like IRF) | * 1D analysis (point-like IRF) | ||
+ | |||
+ | |||
+ | === Offs === | ||
+ | |||
+ | [[File:M87_240609_OFFs.png|thumb|centre|700px| All OFFs]] | ||
+ | |||
+ | |||
+ | [[File:M87_2021_OFFs.png|thumb|centre|700px| 2021 OFFs]] | ||
+ | |||
+ | [[File:M87_2022_OFFs.png|thumb|centre|700px| 2022 OFFs]] | ||
+ | |||
+ | [[File:M87_2023_OFFs.png|thumb|centre|700px| 2023 OFFs]] | ||
+ | |||
+ | [[File:M87_2024_OFFs.png|thumb|centre|700px| 2024 OFFs]] | ||
==Analysis Results== | ==Analysis Results== | ||
Line 151: | Line 177: | ||
Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here. | Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here. | ||
− | === Theta2 | + | === Theta2 === |
+ | |||
+ | [[File:theta2_single_total_M87_240609.png|thumb|centre|720px| Total Theta2]] | ||
+ | |||
+ | [[File:theta2_total_M87_240609 Medium.png|thumb|centre|720px| Theta2 in different energy bins]] | ||
+ | |||
+ | Only the energy bins of 10-40 GeV and 1-5 TeV has over 3 sigma | ||
+ | |||
+ | No individual years provided over 3 sigma | ||
+ | |||
+ | |||
+ | ==== Theta2 2021 ==== | ||
+ | |||
+ | [[File:theta2_single_total_M87_2021.png|thumb|centre|720px| 2021 total Theta2]] | ||
+ | |||
+ | [[File:theta2_total_M87_2021.png|thumb|centre|720px| 2021 theta2 in different energy bins]] | ||
+ | |||
+ | |||
+ | ==== Theta2 2022 ==== | ||
+ | |||
+ | [[File:theta2_single_total_M87_2022.png|thumb|centre|720px| 2022 total Theta2]] | ||
+ | |||
+ | [[File:theta2_total_M87_2022.png|thumb|centre|720px| 2022 theta2 in different energy bins]] | ||
+ | |||
+ | |||
+ | ==== Theta2 2023 ==== | ||
+ | |||
+ | [[File:theta2_single_total_M87_2023.png|thumb|centre|720px| 2023 total Theta2]] | ||
+ | |||
+ | [[File:theta2_total_M87_2023.png|thumb|centre|720px| 2023 theta2 in different energy bins]] | ||
+ | |||
+ | |||
+ | ==== Theta2 2024 ==== | ||
+ | |||
+ | [[File:theta2_single_total_M87_2024.png|thumb|centre|720px| 2024 total Theta2]] | ||
+ | |||
+ | [[File:theta2_total_M87_2024.png|thumb|centre|720px| 2024 theta2 in different energy bins]] | ||
+ | |||
+ | === Skymaps === | ||
+ | |||
+ | Skymaps between 0.2 TeV and 7 TeV | ||
+ | |||
+ | [[File:Significance_Off_and_Excess_maps_M87_240609.png|thumb|centre|700px|| Significance, off significance, and excess maps all observations]] | ||
+ | |||
+ | [[File:Background_residuals_map_M87_240609.png|thumb|centre|700px| Background residuals map all observations]] | ||
+ | |||
+ | |||
+ | ==== Skymaps 2021 ==== | ||
+ | |||
+ | [[File:Significance_Off_and_Excess_maps_M87_2021.png|thumb|centre|700px|| Significance, off significance, and excess maps 2021]] | ||
+ | |||
+ | [[File:Background_residuals_map_M87_2021.png|thumb|centre|700px| Background residuals map 2021]] | ||
+ | |||
+ | |||
+ | |||
+ | ==== Skymaps 2022 ==== | ||
+ | |||
+ | [[File:Significance_Off_and_Excess_maps_M87_2022.png|thumb|centre|700px|| Significance, off significance, and excess maps 2022]] | ||
+ | |||
+ | [[File:Background_residuals_map_M87_2022.png|thumb|centre|700px| Background residuals map 2022]] | ||
+ | |||
+ | |||
+ | |||
+ | ==== Skymaps 2023 ==== | ||
+ | |||
+ | [[File:Significance_Off_and_Excess_maps_M87_2023.png|thumb|centre|700px|| Significance, off significance, and excess maps 2023]] | ||
+ | |||
+ | [[File:Background_residuals_map_M87_2023.png|thumb|centre|700px| Background residuals map 2023]] | ||
+ | |||
+ | |||
+ | |||
+ | ==== Skymaps 2024 ==== | ||
+ | |||
+ | [[File:Significance_Off_and_Excess_maps_M87_2024.png|thumb|centre|700px|| Significance, off significance, and excess maps 2024]] | ||
+ | |||
+ | [[File:Background_residuals_map_M87_2024.png|thumb|centre|700px| Background residuals map 2024]] | ||
+ | |||
+ | === Light-Curve === | ||
+ | |||
+ | [[File:M87_240609_LCs_LST1.png|thumb|centre|720px| Run-wise and night-wise LCs for PWL and LP fluxes with average fluxes, compared to highest and lowest previous M87 PWL fluxes in literature (MAGIC, 2020).]] | ||
+ | |||
+ | [[File:M87_240609_LC_BB_LST1.png|thumb|centre|720px| Run-wise, night-wise, and BB-wise LP fluxes.]] | ||
+ | |||
+ | |||
+ | ==== Light-Curve 2021 ==== | ||
+ | |||
+ | [[File:M87_2021_LCs_LST1.png|thumb|centre|720px| 2021 run-wise and night-wise LCs for PWL and LP fluxes with average fluxes, compared to highest and lowest previous M87 PWL fluxes in literature (MAGIC, 2020).]] | ||
+ | |||
+ | |||
+ | ==== Light-Curve 2022 ==== | ||
+ | |||
+ | [[File:M87_2022_LCs_LST1.png|thumb|centre|720px| 2022 run-wise and night-wise LCs for PWL and LP fluxes with average fluxes, compared to highest and lowest previous M87 PWL fluxes in literature (MAGIC, 2020).]] | ||
+ | |||
+ | |||
+ | ==== Light-Curve 2023 ==== | ||
+ | |||
+ | [[File:M87_2023_LCs_LST1.png|thumb|centre|720px| 2023 run-wise and night-wise LCs for PWL and LP fluxes with average fluxes, compared to highest and lowest previous M87 PWL fluxes in literature (MAGIC, 2020).]] | ||
+ | |||
+ | |||
+ | ==== Light-Curve 2024 ==== | ||
+ | |||
+ | [[File:M87_2024_LCs_LST1.png|thumb|centre|720px| 2024 run-wise and night-wise LCs for PWL and LP fluxes with average fluxes, compared to highest and lowest previous M87 PWL fluxes in literature (MAGIC, 2020).]] | ||
+ | |||
+ | |||
+ | === SED === | ||
+ | |||
+ | [[File:M87_240609_SED_hatch_LST1.png|thumb|centre|720px| SED PWL and LP models compared to previous M87 PWL models in literature (MAGIC, 2020). The vertical shaded area marks the energy range of the literature models. All models use E_0 = 784 GeV as reference.]] | ||
+ | |||
+ | |||
+ | ==== Power law ==== | ||
+ | |||
+ | [[File:M87_240609_comparison_PWL_SED_LST1.png|thumb|centre|720px| SED PWL models per year compared to previous M87 PWL models in literature (MAGIC, 2020). The vertical shaded area marks the energy range of the literature models. All models use E_0 = 784 GeV as reference.]] | ||
+ | |||
+ | {| class="wikitable" style="margin:auto" | ||
+ | |+ M87 PWL models | ||
+ | |- | ||
+ | ! Instrument !! Year !! 𝚪 !! f_0 [1e-12 TeV-1 cm-2 s-1] | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |- | ||
+ | | LST-1 || 2021 || ± || ± | ||
+ | |} | ||
− | |||
− | === | + | ==== Log-Parabola (LP) ==== |
− | = | + | [[File:M87_240609_comparison_LP_SED_LST1.png|thumb|centre|720px| SED LP models per year compared to previous M87 PWL models in literature (MAGIC, 2020). The vertical shaded area marks the energy range of the literature models. All models use E_0 = 784 GeV as reference.]] |
Latest revision as of 15:26, 13 June 2024
Contents
M87 analysis (TTHA)[edit]
General information[edit]
- Name of the source: M87
- Brief description of the source:
- Analysis by Tora T. H. Arnesen (Instituto de Astrofísica de Canarias - <tora.arnesen@iac.es>)
Data Selection[edit]
Data selected with data quality notebook
/fefs/aswg/workspace/tora.arnesen/M87/M87_data_quality.ipynb
- Path to datacheck files
- files = glob.glob("/fefs/aswg/data/real/OSA/DL1DataCheck_LongTerm/night_wise/all/DL1_datacheck_20*")
- Dates of data-selection:
- first_date = 20191123
- last_date = 20240609
- ZD range (deg):
- min_zenith = 0 * u.deg
- max_zenith = 60 * u.deg
- maximum diffuse NSB
- max_diffuse_nsb_std = 2.3
- Pointing stability
- max_pointing_dec_std = 0.011
- Cherenkov Transparency parameters
- min_mean_fit_p = -3.
- max_LS_periodogram_maxamplitude = 1e-2
- min_drdi_index = -2.35
- max_drdi_index = -2.1
- min_drdi_at_422pe = 1.25
- min_fraction_around_mode = 0.8
- max_intensity_at_half_peak_rate = 70
Selected Data[edit]
Total number of runs: 130 out of 234
Total number of nights: 36 out of 49
Total observation time (hours): 33.70 out of 60.15
20210413 : [4396 4397 4398] 20210417 : [4460 4461 4462 4463 4464 4465] 20210419 : [4494 4495] 20210420 : [4508 4509] 20210504 : [4594 4595] 20210507 : [4636 4637] 20220301 : [7185 7186 7187] 20220323 : [7444] 20220407 : [7712 7713 7714 7715 7716] 20220408 : [7746] 20220409 : [7773] 20220503 : [8092] 20220506 : [8209 8211 8212 8213 8214] 20230412 : [12662] 20230418 : [12744 12745 12746 12751 12752 12753 12754] 20230419 : [12775 12776 12778] 20230420 : [12796 12797 12798 12799] 20230614 : [13300 13301 13302 13303] 20240110 : [16315 16316 16317 16318 16319] 20240210 : [16692 16693 16694 16695] 20240212 : [16751 16752 16753 16754 16755] 20240304 : [16960 16961 16962 16963] 20240310 : [17034 17035 17037 17045 17047 17048] 20240311 : [17072 17075 17076 17083 17084 17085 17086] 20240313 : [17103 17104 17105] 20240315 : [17160 17161] 20240502 : [17244 17245 17246] 20240503 : [17253 17254 17255 17256 17258 17259 17261 17262 17263 17264 17265 17266] 20240504 : [17284 17289 17292] 20240507 : [17359] 20240514 : [17403 17404] 20240527 : [17449 17450 17451 17452 17453] 20240601 : [17570] 20240603 : [17612 17613 17614 17615] 20240604 : [17641 17642 17643 17644] 20240607 : [17676 17677 17678 17679 17680 17681]
Monte Carlo information[edit]
- All-sky MC production used: 20240131_allsky_v0.10.5_all_dec_base
- Declination line (training): dec_931
Random forest[edit]
Source declination is 9.31 deg. The two closest MC declination tracks are 9.31 deg and 18.02 deg. Chose dec 9.31 deg as it is closer to the dec of the source at 12.32 deg.
/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931
DL2 Data[edit]
Create DL2 files by running
./M87_launch_dl1_dl2.sh /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931 /fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/dec_931/lstchain_config_2024-01-31.json
DL2 files stored in
/fefs/aswg/workspace/tora.arnesen/M87/DL2
IRF Files[edit]
Use copy of school config file
- 70 % efficiency for both gammaness and theta gamma-ray selection cuts
- Point-like IRF
/fefs/aswg/workspace/tora.arnesen/M87/irf_dl3_tool_config.json
Create IRF files by running
./M87_launch_irf_creation.sh /fefs/aswg/data/mc/DL2/AllSky/20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_931 ./irf_dl3_tool_config.json 931_gh_eff_0.7_th_cont_0.7_20240609_TT
IRF files stored in
/fefs/aswg/workspace/tora.arnesen/M87/IRF
DL3 Data[edit]
Create DL3 files by running
./M87_launch_dl2_dl3.sh irf_dl3_tool_config.json ./DL2 ./IRF/20240131_allsky_v0.10.5_all_dec_base/dec_931/931_gh_eff_0.7_th_cont_0.7_20240609_TT M87
DL3 files stored in
/fefs/aswg/workspace/tora.arnesen/M87/DL3
High-level analysis[edit]
- lstchain v0.10.9 (v0.10.7 for Theta2 and Skymaps)
- Gammapy 1.2 (1.1 for Theta2 and Skymaps)
- 1D analysis (point-like IRF)
Offs[edit]
Analysis Results[edit]
Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here.
Theta2[edit]
Only the energy bins of 10-40 GeV and 1-5 TeV has over 3 sigma
No individual years provided over 3 sigma
Theta2 2021[edit]
Theta2 2022[edit]
Theta2 2023[edit]
Theta2 2024[edit]
Skymaps[edit]
Skymaps between 0.2 TeV and 7 TeV
Skymaps 2021[edit]
Skymaps 2022[edit]
Skymaps 2023[edit]
Skymaps 2024[edit]
Light-Curve[edit]
Light-Curve 2021[edit]
Light-Curve 2022[edit]
Light-Curve 2023[edit]
Light-Curve 2024[edit]
SED[edit]
Power law[edit]
Instrument | Year | 𝚪 | f_0 [1e-12 TeV-1 cm-2 s-1] |
---|---|---|---|
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |
LST-1 | 2021 | ± | ± |