Difference between revisions of "December January paper"

From my_wiki
Jump to: navigation, search
m (→‎LST analysis)
 
(23 intermediate revisions by 3 users not shown)
Line 5: Line 5:
  
 
== Fermi analysis ==
 
== Fermi analysis ==
 +
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 ===
 +
 +
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 ==
 +
 +
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.
 +
 +
[[File:Fermi_Final_LC.png]]
 +
 +
The SEDs corresponding to each block is given here :
 +
 +
[[File:OPGammaSED_BB.png]]
 +
 +
=== Focus on the December / January data ===
 +
 +
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 ==
 
== 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-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 ===
 +
* 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 ==
 +
 +
 +
== Conference plots ==
 +
 +
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
 +
 +
[[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.]]
 +
 +
 +
  
−
===Data-taking information===
+
'''January Data set
−
* observation summary
 
−
{| class="wikitable"
 
−
| Date || Runs || time (h)|| Zd || Az || Transmission || Comments
 
  
−
|-
+
[[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.]]
−
| 2023-12-09 || 15942, 15943, 15944, 15945 || 1.31 || 35-53  || - || - || 50-sec shifted timestamp -> solved!
 
−
|-
 
−
| 2023-12-10 || 15975, 15976, 15977, 15978 || 1.61 || 30-51 || - || - || 
 
−
|-
 
−
| 2023-12-11 || 16007, 16008, 16009, 16010 || 1.39 || 30-54  || - || - || 
 
−
|-
 
−
| 2023-12-12 || 16042, 16043, 16044, 16045 || 1.63 || 30-54  || - || - || 
 
−
|-
 
−
| 2023-12-13 || 16079, 16080, 16081, 16082, 16083 || 1.69 || 30-53  || - || - || 
 
−
|-
 
−
| 2023-12-14 || 16114, 16116, 16117, 16118, 16119|| 1.60 || 30-52  || - || - || 
 
−
|-
 
−
| 2023-12-15 || 16150, 16151, 16152, 16153, 16154, 16155, 16156 || 1.75 || 30-52  || - || - ||
 
−
|-
 
−
| 2023-12-16 || 16188, 16189, 16191, 16192, 16193 || 1.67 ||  25-52 || - || - || Drive issue during run 16190 (no tracking)
 
−
|-
 
−
| 2023-12-17 || 16220, 16221, 16222, 16223, 16224, 16225 || 1.72 ||  25-50 || - || - || 
 
−
|-
 
−
| 2023-12-18 || 16240, 16241, 16242, 16243, 16244, 16246 ||  ||  25-50 || - || - || online gain-selected data
 
−
|-
 
−
| '''2023-12-19 || '''16258, 16259, 16260, 16261, 16262, 16263 ||  ||  25-50 || - || - ||  '''online gain-selected data.''' Removed ?
 
−
|-
 
−
| '''2023-12-20 ||''' 16269, 16270 ||  || ''' 40-43 || - || - ||''' online gain-selected data, humidity is rising during run 16270'''  Removed ?
 
−
|-
 
−
| 2024-01-09 || 16293, 16294, 16295 || 0.63 ||  20-26 || - || - || 
 
−
|-
 
−
| 2024-01-13 || 16350, 16351, 16352, 16353, 16354, 16356, 16357, 16358, 16359 || 2.08 ||  15-43 || - || - || 
 
−
|-
 
−
| 2024-01-17 || 16391, 16392 || 0.46 ||  25-30 || - || - || Timestamp is shifted?? -> solved!
 
−
|-
 
−
| ''' 2024-01-18 ||''' 16407 ||0.10 || 4-6 || - || - || '''  Removed ?
 
−
|-
 
−
| 2024-01-19 || 16412, 16413, 16414, 16415 || 1.34 || 5-30 || - || - ||
 
−
|-
 
−
| 2024-01-21 || 16418, 16419, 16420 || 0.70 || 4-8 || - || - ||
 
−
|}
 

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


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.

File:Fermi Final LC.png

The SEDs corresponding to each block is given here :

OPGammaSED BB.png

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

Fermi-LAT and LST results. Not corrected for EBL.
Fermi-LAT and LST results together with the joint fit on top. Not corrected for EBL.



January Data set

Fermi-LAT and LST results together with the joint fit on top. Not corrected for EBL.
Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=December_January_paper&oldid=9059"