Difference between revisions of "Fast Offline Analysis Template"

From my_wiki
Jump to: navigation, search
(Fast offline analysis template)
(Fast offline analysis template)
Line 27: Line 27:
 
** Date (time of observation h)
 
** Date (time of observation h)
 
***ZD range (deg): Zenith Range
 
***ZD range (deg): Zenith Range
***AZ range (deg): Azimuth Range
+
***Observations: (see example below)
***Transmission9km:
 
***Runs:
 
***Observations:
 
  
  '2024-02-28': [16895, 16896, 16897, 16898]
+
  2024-02-28: [16895, 16896, 16897, 16898]
  '2024-03-03': [16912, 16913, 16917, 16918, 16919, 16920, 16921, 16922, 16923, 16924, 16925, 16926, 16927, 16928, 16929, 16930, 16931, 16932]
+
  2024-03-03: [16912, 16913, 16917, 16918, 16919, 16920, 16921, 16922, 16923, 16924, 16925, 16926, 16927, 16928, 16929, 16930, 16931, 16932]
'2024-03-04': [16953, 16954, 16955, 16956, 16957, 16958]
 
'2024-03-06': [16985, 16986, 16987, 16988, 16989, 16990, 16991]
 
'2024-03-11': [17077, 17078, 17079, 17080, 17081, 17082]
 
'2024-03-13': [17107, 17108, 17109, 17110, 17111, 17112, 17113, 17114, 17115]
 
'2024-03-14': [17124, 17125, 17126, 17127, 17128, 17129, 17130, 17131, 17132, 17133, 17134, 17135, 17136, 17137, 17138]
 
'2024-03-15': [17162, 17163, 17164, 17165, 17166, 17167, 17168, 17169]
 
'2024-03-16': [17176, 17177, 17178, 17179, 17180, 17181, 17182]
 
'2024-03-17': [17189, 17191, 17192, 17193]
 
  
*Data quality selection not yet applied on above data sample
+
*Data quality selection: it is not mandatory to apply data selection for these preliminary analyses. But in case you do it and select good-quality runs, please also add the plot in which the
*Joint observations with MAGIC? : Simultaneous observations but MAGIC GPS not working so timestamp matching between LST-1 and MAGIC cannot be properly done.
+
*Joint observations with MAGIC?:  
*Information on data-taking (observations summaries): https://www.lst1.iac.es/datacheck/observation_overview/2024/03/run_summary_2024_03_17.pdf
+
*Information on data-taking (observations summaries): Was there any problem with data-taking (check ELOG and e.g. add link https://www.lst1.iac.es/datacheck/observation_overview/2024/03/run_summary_2024_03_17.pdf)
  
 
==General analysis information==
 
==General analysis information==
Line 56: Line 45:
  
 
==Monte Carlo information==
 
==Monte Carlo information==
Write below the MC production id among the existing MC
+
Write below the MC production ID among the existing MC
 
* All-sky MC production used: (e.g. 20230901_v0.10.4_allsky_base_prod)
 
* All-sky MC production used: (e.g. 20230901_v0.10.4_allsky_base_prod)
 
** Declination line (training): dec_XXXX
 
** Declination line (training): dec_XXXX
Line 62: Line 51:
  
 
===Random forest===
 
===Random forest===
Path to the models used:
+
Path to the models used: (example below)
  
 
  /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_3476
 
  /fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_3476

Revision as of 14:27, 2 May 2024

Fast offline analysis template

Name your fast offline analysis wiki page following the scheme:

Fast_Offline_Analysis_SourceName_Month(s)_Year_ObsPeriod

where:

  • ObsPeriod should be e.g. P063

For example: https://www.lst1.iac.es/wiki/index.php/Fast_Analysis_OP313_ToO_Feb_Mar_2024_P062

General information

  • Name of the source:
  • Brief description of the source:
    • Object type :
    • Redshift (z) :
    • Other relevant information
    • RA: XX XX XX.XXX (hh mm ss), Dec: YY YY YY.YYY (dd mm ss)
    • RA, Dec in deg (ICRS): XX.XX, YY.YY
  • Analysis by:
    • Analyzer 1 name (Institute - email)
    • Analyzer 2 name (Institute - email)

Data-taking information

  • Dates of data-taking:
    • Date (time of observation h)
      • ZD range (deg): Zenith Range
      • Observations: (see example below)
2024-02-28: [16895, 16896, 16897, 16898]
2024-03-03: [16912, 16913, 16917, 16918, 16919, 16920, 16921, 16922, 16923, 16924, 16925, 16926, 16927, 16928, 16929, 16930, 16931, 16932]

General analysis information

  • source-independent analysis
  • lstchain version: (e.g. v0.10.10)
  • DL1b standard production (v0.10/tailcut84), run-wise merged files in:
/fefs/aswg/data/real/DL1/YYYYMMDD/v0.10/tailcut84/dl1_LST-1.Run?????.h5


Monte Carlo information

Write below the MC production ID among the existing MC

  • All-sky MC production used: (e.g. 20230901_v0.10.4_allsky_base_prod)
    • Declination line (training): dec_XXXX


Random forest

Path to the models used: (example below)

/fefs/aswg/data/models/AllSky/20230901_v0.10.4_allsky_base_prod/dec_3476

DL2 data

Config file:


Files in:

/fefs/aswg/workspace/daniel.morcuende/data/real/DL2/OP313_ToO_Dec2023/dec_3476

DL3 data & IRF

 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.01, 1],
     "leakage_intensity_width_2": [0, 1],
     "event_type": [32, 32]
   }
 },
 "DL3Cuts": {
   "min_event_p_en_bin": 100,
   "min_gh_cut": 0.1,
   "max_gh_cut": 0.95,
   "min_theta_cut": 0.0,
   "max_theta_cut": 0.2,
   "fill_theta_cut": 0.2,
   "allowed_tels": [1]
 },
  • 70 % efficiency for both gammaness and theta gamma-ray selection cuts
  • Point-like IRF

High-level analysis

  • Gammapy 1.1
  • 1D analysis (point-like IRF)

Analysis Results

Theta2 plot

Settings:

  • Theta2 cut for signal extraction: X.XX deg2
  • Normalization range: X.X - Y.Y deg2
  • Number of off positions:
  • Energy range:

Daily results

  • YYYY-MM-DD

600 px

Stacked observations

Significance and Excess sky maps

Spectral results

Include the spectral parameters (see example below) and the SED plot (best-fit model + spectral flux points)

For example:

  type      name      value         unit        error      min    max frozen is_norm link
-------- ---------- ---------- -------------- --------- --------- --- ------ ------- ----
spectral      index 2.8255e+00                3.298e-01 1.000e+00 nan  False   False     
spectral  amplitude 1.7730e-09 TeV-1 s-1 cm-2 4.988e-10       nan nan  False    True     
spectral  reference 1.0000e+02            GeV 0.000e+00       nan nan   True   False       


Light curve