Editing HEPRO VIII GRB221009A analysis (dark)

=Dark analysis=

This page explains the analysis (source independent) done to produce the high-level results from GRB 221009A for the HEPRO VIII conference. The analysis uses the dark observations of this GRB (see below)

==General information==
* '''Name of the source''': GRB 221009A
* '''Brief description of the source''': 
** Object type : GRB
** Redshift (z) : 0.15
** Other relevant information 
** RA, Dec in deg (ICRS): 288.26452 deg, +19.77350 deg

* Analysis by Arnau Aguasca-Cabot (UB - arnau.aguasca@fqa.ub.edu)

==Data-taking information==
** '''2022-10-15'''
*** Zenith range (deg): 23.8 -> 57.9
*** Transmission9km:   n/a 
*** Observation condition: dark and moon (9629,9630)
*** Observation mode: std wobble
*** Runs: 9625 - 9628
*** Joint observations with MAGIC? : Yes, except the first 2 runs
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2441]
*** Important comments: Data were affected by timestamp issue. (Reprocessing is already done and results below are from reprocessed data)

** '''2022-10-16'''
*** Zenith range (deg): 32.1 -> 64.1
*** Transmission9km:   n/a 
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9640 - 9647
*** Joint observations with MAGIC? : No
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2443]
*** Important comments: Data were affected by timestamp issue. (Reprocessing is already done and results below are from reprocessed data)

** '''2022-10-17'''
*** Zenith range (deg): 25.9 -> 58.6
*** Transmission9km:   n/a 
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9670 - 9677
*** Joint observations with MAGIC? : Yes
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2445]
*** Important comments: Data were affected by timestamp issue. (Reprocessing is already done and results below are from reprocessed data)

** '''2022-10-18'''
*** Zenith range (deg): 23.1 -> 43.5
*** Transmission9km:   n/a 
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9691 - 9696
*** Joint observations with MAGIC? : No
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2447]
*** Important comments: unstable rate in run 9691

** '''2022-10-21'''
*** Zenith range (deg): 34.6 -> 55.6
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9802 - 9807
*** Joint observations with MAGIC? : Yes
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2451]
*** Important comments:

** '''2022-10-23'''
*** Zenith range (deg): 32.07 -> 59.21
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9829 - 9835
*** Joint observations with MAGIC? : Yes
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2454]
*** Important comments: 

** '''2022-10-24'''
*** Zenith range (deg): 49.6 -> 58.8
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9861 - 9863
*** Joint observations with MAGIC? : No
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2454]
*** Important comments: some clouds

** '''2022-10-25'''
*** Zenith range (deg): 43.6 -> 56
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9888 - 9891
*** Joint observations with MAGIC? : Yes
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2458]
*** Important comments: 

** '''2022-10-26'''
*** Zenith range (deg): 40 -> 56
*** Observation condition: dark
*** Observation mode: std wobble
*** Runs: 9917 - 9920
*** Joint observations with MAGIC? : Only run 9919
*** ELOG link: [https://www.lst1.iac.es/elog/LST+commissioning/2460]
*** Important comments: run 9918,9919 are both in W2

==Monte Carlo information==
* MC production reprocessed with lstMCpipe (v0.9.9) starting at DL1a level for AllSky production. NSB tuning added al DL1 level.
* Link to lstMCpipe production: https://github.com/cta-observatory/lstmcpipe/tree/master/production_configs/20221027_v0.9.9_crab_tuned

* Link to MC files used: /fefs/aswg/data/mc/*/AllSky/20221027_v0.9.9_crab_tuned
** Particle types: Diffuse, point-like gammas, diffuse protons
** dec band (deg): 22.67

==DL1 data==

Reprocessed using DL1a files produced by LSTOSA. The files were reprocessed solving the problem with the time stamp.

* original DL1a files (real)
* lstchain v0.9.6
* real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning

 /fefs/aswg/data/real/DL1/*/v0.9/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5

* MC: tailcut8-4, dynamic cleaning, NSB tuning
* lstchain v0.9.9

  "image_modifier": {
    "increase_nsb": true,
     "extra_noise_in_dim_pixels": 1.62,
     "extra_bias_in_dim_pixels": 0.655,
     "transition_charge": 8,
     "extra_noise_in_bright_pixels": 2.08,
     "increase_psf": false,
     "smeared_light_fraction": 0
  },

 /fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/*

==Random forest==

* lstchain-0.9.9
* source-independent analysis (diffuse gamma, disp_norm)
 https://github.com/cta-observatory/lstmcpipe/tree/master/production_configs/20221027_v0.9.9_crab_tuned

==DL2 data==

* lstchain-0.9.9
* source-independent analysis
 /fefs/aswg/workspace/arnau.aguasca/Analysis/results/LST/real/DL2/202210*/srcind/v0.9.9/tailcut84_dynclg_tNSBtoCrabPerfPaper/AllSky_20221027_v0.9.9_crab_tuned_dec_2276/

==DL3 data selection==

DL3 files produced considering the edisp normalisation correction (using the fix_edisp.py script) and the interpolation mode to compute the IRFs with energy-dependent cuts (both for the gammaness and theta cut).

* lstchain-0.9.14dev (commit: 90d0450, after IRF interpolation method merged to main branch)

  "filters": {
    "intensity": [50, Infinity],
    "width": [0, Infinity],
    "length": [0, Infinity],
    "r": [0, 1],
    "wl": [0.01, 1],
    "leakage_intensity_width_2": [0, 0.2],
    "event_type": [32, 32]
  }
  "gh_efficiency": 0.9,
  "theta_containment": 0.7,

==High-level analysis==

*Science Tool: gammapy 1.0
* point-like IRF, 1D analysis
* Energy threshold based on Crab check (200 GeV)

==Analysis Results==

=== Spectral results ===
<div><ul> 
<li style="display: inline-block;"> [[File:HEPROVIII_SED_GRB221009A_alldataDark.png|thumb|none|400px|Spectral energy distribution using all observations in dark conditions. The Crab spectrum is shown as a dashed red line.]] </li>
</ul></div>

=== Light curve results ===
<div><ul> 
<li style="display: inline-block;"> [[File:v2_HEPROVIII_eflux_GRB221009A_200GeV.png|thumb|none|400px|Energy flux for LST-1 dark observations (Oct. 15-18 and 21-28). The ULs from HAWC (green) and H.E.S.S. (orange) are shown. The dashed line marks the Fermi-GBM trigger time (t0).]] </li>
</ul></div>

== Crab check ==
The quality cuts used to produce the high-level results are based on the analysis of Crab data in similar conditions to the GRB 221009A data in dark conditions
* Crab runs (2022-03-04) 
:* Runs: 7255, 7256
* Crab runs (2022-03-05) 
:* Runs: 7277, 7278, 7279, 7280, 7281, 7282
* Crab runs (2022-11-21) 
:* Runs: 10914, 10915, 10916, 10917, 10918, 10986, 10987
* Crab runs (2022-11-23) 
:* Runs: 10988, 10989

<div><ul> 
<li style="display: inline-block;"> [[File:crabcheck_SED_dark_int50_geff0.9_teff0.7_ethd0.1_interp_node_corrected_edisp.png|thumb|none|400px|Spectral energy distribution.]] </li>
</ul></div>

<div><ul> 
<li style="display: inline-block;"> [[File:crabcheck_LC_dark_int50_geff0.9_teff0.7_ethd0.1_interp_node_corrected_edisp.png|thumb|none|400px|Light curve above 200 GeV.]] </li>
</ul></div>


== Cross-check by Kenta ==
The plots below show the cross-check by Kenta (source-dependent analysis). No edisp normalisation correction was applied.

<div><ul> 
<li style="display: inline-block;"> [[File:GRB221009A_SED_kenta_dark.png|thumb|none|400px|Spectral energy distribution.]] </li>
</ul></div>

<div><ul> 
<li style="display: inline-block;"> [[File:GRB221009A_LC_kenta_dark.png|thumb|none|400px|Light curve.]] </li>
</ul></div>