Difference between revisions of "December January paper"
Jorge.otero (talk | contribs) (→Conference plots) |
m (→LST analysis) |
||
| (12 intermediate revisions by 2 users not shown) | |||
| Line 9: | Line 9: | ||
=== Parameters === | === Parameters === | ||
| − | RA = 197.6324 degrees | + | RA = 197.6324 degrees |
| − | + | Dec = 32.3547 degrees | |
| − | Dec = 32.3547 degrees | + | Start = 715617241.0 MET (s) |
| − | + | Stop = 727281241.0 MET (s) | |
| − | Start = 715617241.0 MET (s) | + | ROI = 10.0 degrees |
| − | + | E min = 100.0 MeV | |
| − | Stop = 727281241.0 MET (s) | + | E max = 300000.0 MeV |
| − | + | E min ext = 50.11872336272722 MeV | |
| − | ROI = 10.0 degrees | + | E max ext = 598578.6944906638 MeV |
| − | + | IRFs = P8R3_SOURCE_V3 | |
| − | E min = 100.0 MeV | + | evclass = 128 |
| − | + | evtype = 3 | |
| − | E max = 300000.0 MeV | ||
| − | |||
| − | E min ext = 50.11872336272722 MeV | ||
| − | |||
| − | E max ext = 598578.6944906638 MeV | ||
| − | |||
| − | IRFs = P8R3_SOURCE_V3 | ||
| − | |||
| − | evclass = 128 | ||
| − | |||
| − | evtype = 3 | ||
| Line 53: | Line 42: | ||
== LST analysis == | == LST analysis == | ||
| + | Find all the analysis details, including crosscheck plots at [[OP313_Dec2023_Jan2024_LST1_analysis_details]] | ||
| + | |||
| + | === Data-taking information === | ||
| + | |||
| + | * NSB level: all data was taken in dark NSB conditions | ||
| + | * Zenith angles: 15-54 deg (December 2023 data have zenith angle > 30 deg) | ||
| + | |||
| + | === Good-quality selected dataset === | ||
| + | |||
| + | The list of data runs below was obtained using the notebook https://github.com/cta-observatory/cta-lstchain/blob/main/notebooks/data_quality.ipynb. | ||
| − | = | + | We used default settings, but <code>min_drdi_at_422pe = 1.2</code> instead of the default 1.5 to recover 1.7 h from January 2024 where low statistics limit us. Still, the data selected have a relative light yield above ~80%. |
| − | |||
| + | Also, short data runs (< 5 min) were excluded, because usually they are a symptom of problems during data acquisition. | ||
| − | |||
'2023-12-09': [15942, 15943, 15944, 15945], | '2023-12-09': [15942, 15943, 15944, 15945], | ||
| Line 75: | Line 73: | ||
'2024-01-19': [16412, 16413, 16414, 16415], | '2024-01-19': [16412, 16413, 16414, 16415], | ||
'2024-01-21': [16418, 16419, 16420] | '2024-01-21': [16418, 16419, 16420] | ||
| + | |||
| + | === lstchain settings === | ||
| + | * lstchain version: v0.10.7 | ||
| + | * MC production: 20230901_v0.10.4_allsky_base_prod | ||
| + | * RF models: /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_3476/ | ||
| + | * DL2 MC for IRF calculation: /fefs/aswg/data/mc/DL2/AllSky/20230901_v0.10.4_allsky_base_prod/TestingDataset/dec_3476 | ||
| + | * Used IRF interpolation | ||
| + | * Energy threshold: ~ 40 GeV | ||
| + | * Selection cuts IRF and DL3: 70% efficiency in both gammaness and theta gamma-ray selection cuts. Intensity > 50 p.e. | ||
| + | |||
| + | === Gammapy settings === | ||
| + | * Gammapy version: v1.1 | ||
| + | * Energy binning | ||
| + | * Safe mask: 5% of max. Aeff | ||
| + | * Energy range for spectral fitting: 70 GeV to 1 TeV | ||
| + | * EBL model: Saldaña-Lopez 2021 | ||
| + | |||
| + | === Cross-checks === | ||
| + | LST-1 data analysis performed by Seiya Nozaki, Jorge Otero Santos, Daniel Morcuende and Mireia Nievas Rosillo | ||
| + | |||
| + | Find all the analysis details, including crosscheck plots at [[OP313_Dec2023_Jan2024_LST1_analysis_details]] | ||
== Summed fit of Fermi-LAT and LST data == | == Summed fit of Fermi-LAT and LST data == | ||
| Line 83: | Line 102: | ||
For 2024 conferences, we decided to have a common set of plots to be shown, including theta2 plots, high and low state lightcures and SEDs showing the LST-1 spectrum, LAT spectrum, and joint fit butterfly shape (but no parameters or spectral points from the joint fit). | For 2024 conferences, we decided to have a common set of plots to be shown, including theta2 plots, high and low state lightcures and SEDs showing the LST-1 spectrum, LAT spectrum, and joint fit butterfly shape (but no parameters or spectral points from the joint fit). | ||
| − | |||
| − | |||
| − | Fermi-LAT and LST results together with the joint fit on top. Not corrected for EBL. | + | '''December Data set |
| + | |||
| + | [[File:OPGammaSED_LSTTime.png|thumb|centre|400px|Fermi-LAT and LST results. Not corrected for EBL.]] | ||
| + | |||
| + | [[File:OPGammaSED_LSTTime_summed.png|thumb|centre|400px|Fermi-LAT and LST results together with the joint fit on top. Not corrected for EBL.]] | ||
| + | |||
| + | |||
| + | |||
| + | '''January Data set | ||
| − | + | [[File:OPGammaSED_LSTTime_summed_Jan.png|thumb|centre|400px|Fermi-LAT and LST results together with the joint fit on top. Not corrected for EBL.]] | |
Latest revision as of 15:20, 22 February 2024
This page is devoted to the first paper on OP313 following the LST detection and the ATel
The aim is to focus on the December 2023 and January 2024 data set and publish this data set in a quick paper
Contents
Fermi analysis[edit]
The Fermi LAT analysis has been performed on a large time scale to derive a long-term light curve and run a Bayesian block analysis to defined bocks of observation
Parameters[edit]
RA = 197.6324 degrees Dec = 32.3547 degrees Start = 715617241.0 MET (s) Stop = 727281241.0 MET (s) ROI = 10.0 degrees E min = 100.0 MeV E max = 300000.0 MeV E min ext = 50.11872336272722 MeV E max ext = 598578.6944906638 MeV IRFs = P8R3_SOURCE_V3 evclass = 128 evtype = 3
Long term Fermi-LAT Light curve[edit]
For more information, see https://www.lst1.iac.es/wiki/index.php/Fermi_LAT_Main_OP313
This plot is the Fermi-LAT LC with the final bayesian block description.
The SEDs corresponding to each block is given here :
Focus on the December / January data[edit]
The Fermi analysis has been run on the LST GTI and lead to a detection at TS=485.58. Power-law and logparabola model were tested and the later is prefered at a TS 0f 8.
File:SED OP313 LogParabola LSTTime Dec Jan.png
LST analysis[edit]
Find all the analysis details, including crosscheck plots at OP313_Dec2023_Jan2024_LST1_analysis_details
Data-taking information[edit]
- NSB level: all data was taken in dark NSB conditions
- Zenith angles: 15-54 deg (December 2023 data have zenith angle > 30 deg)
Good-quality selected dataset[edit]
The list of data runs below was obtained using the notebook https://github.com/cta-observatory/cta-lstchain/blob/main/notebooks/data_quality.ipynb.
We used default settings, but min_drdi_at_422pe = 1.2 instead of the default 1.5 to recover 1.7 h from January 2024 where low statistics limit us. Still, the data selected have a relative light yield above ~80%.
Also, short data runs (< 5 min) were excluded, because usually they are a symptom of problems during data acquisition.
'2023-12-09': [15942, 15943, 15944, 15945], '2023-12-10': [15975, 15976, 15977, 15978], '2023-12-11': [16007, 16008, 16009, 16010], '2023-12-12': [16042, 16043, 16044, 16045], '2023-12-13': [16079, 16080, 16081, 16082, 16083], '2023-12-14': [16114, 16116, 16117, 16118, 16119], '2023-12-15': [16150, 16151, 16152, 16153, 16154, 16155, 16156], '2023-12-16': [16188, 16189, 16191, 16192, 16193], '2023-12-17': [16220, 16221, 16222, 16223, 16224, 16225], '2023-12-18': [16240, 16241, 16242, 16243, 16244, 16246], '2024-01-09': [16293, 16294, 16295], '2024-01-13': [16350, 16351, 16352, 16353, 16354, 16356, 16357, 16358, 16359], '2024-01-17': [16391, 16392], '2024-01-19': [16412, 16413, 16414, 16415], '2024-01-21': [16418, 16419, 16420]
lstchain settings[edit]
- lstchain version: v0.10.7
- MC production: 20230901_v0.10.4_allsky_base_prod
- RF models: /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_3476/
- DL2 MC for IRF calculation: /fefs/aswg/data/mc/DL2/AllSky/20230901_v0.10.4_allsky_base_prod/TestingDataset/dec_3476
- Used IRF interpolation
- Energy threshold: ~ 40 GeV
- Selection cuts IRF and DL3: 70% efficiency in both gammaness and theta gamma-ray selection cuts. Intensity > 50 p.e.
Gammapy settings[edit]
- Gammapy version: v1.1
- Energy binning
- Safe mask: 5% of max. Aeff
- Energy range for spectral fitting: 70 GeV to 1 TeV
- EBL model: Saldaña-Lopez 2021
Cross-checks[edit]
LST-1 data analysis performed by Seiya Nozaki, Jorge Otero Santos, Daniel Morcuende and Mireia Nievas Rosillo
Find all the analysis details, including crosscheck plots at OP313_Dec2023_Jan2024_LST1_analysis_details
Summed fit of Fermi-LAT and LST data[edit]
Conference plots[edit]
For 2024 conferences, we decided to have a common set of plots to be shown, including theta2 plots, high and low state lightcures and SEDs showing the LST-1 spectrum, LAT spectrum, and joint fit butterfly shape (but no parameters or spectral points from the joint fit).
December Data set
January Data set
