Difference between revisions of "Analysis A. Aguasca-Cabot lstchain v0.9.X"

From my_wiki
Jump to: navigation, search
 
(20 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
= Overview =
 
= Overview =
−
*Source-dependent analysis using LSTCHAIN v0.9.1 and for stage DL0 to DL1 using lstMCpipe v0.5.1 commit a2b7cd26b94c8a91915140d8b1c0dce04f48a4ce.  
+
*Source-dependent/independent analysis using LSTCHAIN v0.9.4 (MC production processed with LSTCHAIN v0.9.3, R0 to DL1 process by LSTOSA using LSTCHAIN v0.9.2)
−
*Cleaning method tail cut 84 and dynamic cleaning. Also, tuned NSB and PSF.
+
*Cleaning method tail cut 84 and dynamic cleaning. Also, the MC files are tuned to match the NSB in the FoV and the PSF.
 
*Observations used in the analysis:
 
*Observations used in the analysis:
  
Line 9: Line 9:
 
:- ZD (deg): 40
 
:- ZD (deg): 40
 
:- AZ (deg): 180
 
:- AZ (deg): 180
−
:- Other information: DL1 MC files are produced with lstMCpipe v0.5.1 last commit a2b7cd26b94c8a91915140d8b1c0dce04f48a4ce (using LSTCHAIN v0.9.1).
+
:- Other information: MC production through PR in lstMCpipe Github webpage (https://github.com/cta-observatory/lstmcpipe/pull/185)
−
:- Directory where config files used in lstMCpipe are:
+
:- MC files are tuned to match the NSB in the FoV and the PSF.
−
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.9.1/
 
−
:- Link to DL0 MC files used:
 
−
/fefs/aswg/data/mc/DL0/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing
 
  
 
== DL1 data ==  
 
== DL1 data ==  
−
* LSTCHAIN v0.9.1
+
* Processed with LSTOSA using LSTCHAIN v0.9.2.
−
**Source-dependent parameters added using script `lstchain_add_source_dependent_parameters.py`.
+
:- Dynamic cleaning and tail cut cleaning with pedestal threshold applied by the the standard parameters using pipeline LSTOSA.
−
**Dynamic cleaning applied using script `lstchain_dl1ab.py`.
 
  
 
   "dynamic_cleaning": {
 
   "dynamic_cleaning": {
Line 25: Line 21:
 
     "fraction_cleaning_intensity": 0.03
 
     "fraction_cleaning_intensity": 0.03
 
   }
 
   }
 +
 
=== Real data ===
 
=== Real data ===
  
−
* Original DL1a files
+
* Parameters for the tail cut cleaning with pedestal threshold
−
  /fefs/aswg/data/real/DL1/{}/v0.9.1
+
 
 +
  "tailcuts_clean_with_pedestal_threshold": {
 +
    "picture_thresh":8,
 +
    "boundary_thresh":4,
 +
    "sigma":2.5,
 +
    "keep_isolated_pixels":false,
 +
    "min_number_picture_neighbors":2,
 +
    "use_only_main_island":false,
 +
    "delta_time": 2
 +
  }
 +
 
 +
* Original DL1a files (processed up to DL1 by LSTOSA using LSTCHAIN v0.9.2)
 +
  /fefs/aswg/data/real/DL1/{}/v0.9/tailcut84/dl1_LST-1.RunXXXXX.XXXX.h5
  
 
=== MC data ===
 
=== MC data ===
  
−
* Tunned NSB:
+
* Tunned NSB
−
** MC:
 
−
** DL1a run:
 
  
 
     "increase_nsb": true,
 
     "increase_nsb": true,
−
     "extra_noise_in_dim_pixels":,
+
     "extra_noise_in_dim_pixels": 0.6,
−
     "extra_bias_in_dim_pixels": ,
+
     "extra_bias_in_dim_pixels": 0.269,
 
     "transition_charge": 8,
 
     "transition_charge": 8,
−
     "extra_noise_in_bright_pixels": ,
+
     "extra_noise_in_bright_pixels": 0.566,
  
 
* Tunned PSF
 
* Tunned PSF
Line 48: Line 55:
 
   
 
   
 
* Original DL1a files
 
* Original DL1a files
−
  /fefs/aswg/workspace/arnau.aguasca/Analysis/results/MC/lstMCpipe/DL1/20200629_prod5_trans_80/{particle}/{zenith}/{azimuth}/
+
  /fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/{}/
−
/{if gamma: off_0.4deg}/
 
  
 
== Random forest ==
 
== Random forest ==
−
* LSTCHAIN v0.9.1
+
 
−
* Dispnorm
+
* Dispnorm parameterisation
 +
 
 
*; Source-dependent  
 
*; Source-dependent  
−
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/XXXXXXXXX
+
:- LSTCHAIN v0.9.4
 +
:- Config file:
 +
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
 +
:- path:
 +
/fefs/aswg/workspace/arnau.aguasca/Analysis/models/20200609_XXXX/srcdep/v0.9.4/
  
 
*; Source-independent  
 
*; Source-independent  
−
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/XXXXXXXXX
+
:- Processed with lstMCpipe using LSTCHAIN v0.9.3.
 +
:- Path:
 +
  /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/
  
 
== DL1 to DL2 data ==
 
== DL1 to DL2 data ==
−
* LSTCHAIN v0.9.1
 
−
* Config file
 
  
 +
=== Real data ===
 
*; Source-dependent  
 
*; Source-dependent  
−
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/XXXXXXXXX
+
:- LSTCHAIN v0.9.4
 +
:- Config file:
 +
  //fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
  
 
*; Source-independent
 
*; Source-independent
−
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/XXXXXXXXX
+
:- LSTCHAIN v0.9.4
 +
:- Config file:
 +
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_standard_config_tailcut84_dispnorm_v0.9.3.json
 +
 
 +
=== MC data ===
 +
 
 +
*; Source-dependent
 +
:- LSTCHAIN v0.9.4
 +
:- Config file:
 +
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
 +
 
 +
*; Source-independent
 +
:- Processed with lstMCpipe using LSTCHAIN v0.9.3.
 +
:- Config file:
 +
/fefs/aswg/data/mc/analysis_logs/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/off0.4deg/lstchain_config_NSBtuning.json
 +
 
  
 
== DL3 data selection ==
 
== DL3 data selection ==
Line 75: Line 104:
 
   "EventSelector": {
 
   "EventSelector": {
 
     "filters": {
 
     "filters": {
−
       "intensity": [XXX, Infinity],
+
       "intensity": [50, Infinity],
 
       "width": [0, Infinity],
 
       "width": [0, Infinity],
 
       "length": [0, Infinity],
 
       "length": [0, Infinity],
 
       "r": [0, 1],
 
       "r": [0, 1],
 
       "wl": [0.1, 1],
 
       "wl": [0.1, 1],
−
       "leakage_intensity_width_2": [0, XXX]
+
       "leakage_intensity_width_2": [0, 0.2]
 
     }
 
     }
 
   },
 
   },
 
   "DL3FixedCuts": {
 
   "DL3FixedCuts": {
−
     "fixed_gh_cut": ,
+
     "fixed_gh_cut": 0.8,
−
     "fixed_alpha_cut": ,
+
     "fixed_alpha_cut": 15,
 
     "allowed_tels": [1]
 
     "allowed_tels": [1]
 
   }
 
   }
Line 95: Line 124:
 
*:- IRF: point-like, single-offset
 
*:- IRF: point-like, single-offset
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/XXXX
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/XXXX
−
*:- Produced DL3 files
+
::- Produced DL3 files
  
 
  /fefs/aswg/workspace/arnau.aguasca/XXXX
 
  /fefs/aswg/workspace/arnau.aguasca/XXXX
Line 106: Line 135:
 
=== Alpha plot ===
 
=== Alpha plot ===
 
* Filters in EventSelector applied:
 
* Filters in EventSelector applied:
 +
  "EventSelector": {
 +
    "filters": {
 +
      "intensity": [50, Infinity],
 +
      "width": [0, Infinity],
 +
      "length": [0, Infinity],
 +
      "r": [0, 1],
 +
      "wl": [0.1, 1],
 +
      "leakage_intensity_width_2": [0, 0.2]
 +
    }
 +
  }
 +
See section [[#Crab Check|Crab Check]] to know more about the justification of the applied cuts.
 +
Stacked alpha plots:
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Alpha plot 09to12.png|thumb|none|400px|Alpha plot using runs in 20210809, 20210810 and 20210812.]] </li>
 +
<li style="display: inline-block;"> [[File:Cumulative SNr sig gc08 ac15.png|thumb|none|450px|Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812]] </li>
 +
</ul></div>
 +
Daily alpha plots:
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210809 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210809.]] </li>
 +
<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210810 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210810.]] </li>
 +
<li style="display: inline-block;"> [[File:Alphaplot RSOph 1off 20210812 gc08 ac15.png|thumb|none|400px|Alpha plot using runs in 20210812.]] </li>
 +
</ul></div>
 +
 +
* Energy threshold for the analysis
 +
[[File:Energy thd srcdep RSOph int50 gcut08 acut15.jpg|400px]]
  
 
=== Theta plot ===
 
=== Theta plot ===
 +
Stacked theta plots:
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Theta2 plot 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210809, 20210810 and 20210812.]] </li>
 +
<li style="display: inline-block;"> [[File:Cumulative SNr sqrt-ts 3off gc06 tc02.png|thumb|none|450px|Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812]] </li>
 +
</ul></div>
 +
Daily theta plots:
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210809 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210809.]] </li>
 +
<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210810 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210810.]] </li>
 +
<li style="display: inline-block;"> [[File:Thetaplot RSOph 20210812 3off gc06 tc02.png|thumb|none|400px|Theta2 plot using runs in 20210812.]] </li>
 +
 +
</ul></div>
 +
 +
 +
* Energy threshold for the analysis
 +
[[File:Energy thd srcind RSOph int50 gcut06 tcut02.jpg|400px]]
  
 
=== Spectral results ===
 
=== Spectral results ===
  
 
*; Source-dependent  
 
*; Source-dependent  
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:RSOph SED 09to12 int50 gc08 ac15.png|thumb|none|400px|Stacked SED (Aug 9-12)]] </li>
 +
</ul></div>
 +
:Daily SEDs
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Srcdep spec 20210809 gc08 ac15.png|thumb|none|400px|20210809 SED]] </li>
 +
<li style="display: inline-block;"> [[File:Srcdep spec 20210810 gc08 ac15.png|thumb|none|400px|20210810 SED]] </li>
 +
<li style="display: inline-block;"> [[File:Srcdep spec 20210812 gc08 ac15.png|thumb|none|400px|20210812 SED]]</li>
 +
</ul></div>
 +
 +
:*; HESS and MAGIC comparison
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:LST-1HESSLAT SED RSOph srcdep 09to12 int50 gc08 ac15.png|thumb|none|400px|HESS/Fermi and LST-1 SED.]] </li>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 09to12 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9-12).]] </li>
 +
</ul></div>
 +
 +
::Daily MAGIC/LST-1 SEDs:
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210809 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9).]] </li>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210810 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 10).]] </li>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcdep 20210812 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 12).]] </li>
 +
</ul></div>
 +
 +
 
*; Source-independent  
 
*; Source-independent  
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Srcind spec 20210809to12 gc06 tc02.png|thumb|none|400px|Stacked SED (Aug 9-12)]] </li>
 +
</ul></div>
 +
:Daily SEDs
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Srcind spec 20210809 gc06 tc02.png|thumb|none|400px|20210809 SED]] </li>
 +
<li style="display: inline-block;"> [[File:Srcind spec 20210810 gc06 tc02.png|thumb|none|400px|20210810 SED]] </li>
 +
<li style="display: inline-block;"> [[File:Srcind spec 20210812 gc06 tc02.png|thumb|none|400px|20210812 SED]]</li>
 +
</ul></div>
 +
 +
:*; HESS and MAGIC comparison
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:LST-1HESSLAT SED RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|400px|HESS/Fermi and LST-1 SED.]] </li>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC SED RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|450px|MAGIC and LST-1 SED (Aug 9-12).]] </li>
 +
</ul></div>
 +
 +
No daily SED can be plotted because the SED from LST-1 using gc=0.6 and tc=0.2 deg are not reliable.
  
 
=== Light curve results ===
 
=== Light curve results ===
  
 
*; Source-dependent  
 
*; Source-dependent  
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:LST-1HESSLAT lightcurve RSOph srcind 09to12 int50 gc08 ac15.png|thumb|none|500px|HESS/Fermi and LST-1 light curve.]] </li>
 +
<li style="display: inline-block;"> [[File:LST-1MAGIC lightcurve RSOph srcind 09to12 int50 gc08 ac15.png|thumb|none|450px|MAGIC and LST-1 light curve.]] </li>
 +
</ul></div>
 
*; Source-independent  
 
*; Source-independent  
 +
:;- Notice that since the daily SEDs with quality cuts (gc=0.6 and tc=0.2 deg) are not reliable, the daily flux are computed with the stacked SED using data between 20210809-12!!!
 +
:;- No comparison with MAGIC light curve is done because the energy threshold using (gc=0.6 and tc=0.2 deg) is higher than 100 GeV.
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:LST-1HESSLAT lightcurve RSOph srcind 09to12 int50 gc06 tc02.png|thumb|none|500px|HESS/Fermi and LST-1 light curve.]] </li>
 +
</ul></div>
 +
 +
  
 
= Crab Check =
 
= Crab Check =
 +
Run selection using Crab_data_selection_notebook: https://indico.cta-observatory.org/event/3984/contributions/32878/attachments/20998/29572/Crab_data_selection_notebook.pdf
 +
 +
Threshold for the rate30_mask is set to 4.3. I exclude run 6016 because I obtain a strange SED at low energies.
  
 
== Data reduction Specifications ==
 
== Data reduction Specifications ==
−
*; Source-dependent  
+
*; Source-dependent
 +
 
 +
[[File:SNr sig diff ghc ac.png|600px]]
 +
 
 +
Cuts at alpha cut = 15 deg are studied.
 +
 
 +
* Energy threshold for the analysis
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Energy thd srcdep Crab int50 gcut06 acut15.jpg|thumb|none|800px|gammaness cut>0.6 and alpha cut<15 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcdep Crab int50 gcut07 acut15.jpg|thumb|none|800px|gammaness cut>0.7 and alpha cut<15 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcdep Crab int50 gcut08 acut15.jpg|thumb|none|800px|gammaness cut>0.8 and alpha cut<15 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcdep Crab int100 gcut09 acut15.jpg|thumb|none|800px|gammaness cut>0.9 and alpha cut<15 deg, int>100 p.e.]]</li>
 +
</ul></div>
 +
 
 
*; Source-independent
 
*; Source-independent
 +
 +
[[File:SNr sig diff ghc tc.png|600px]]
 +
 +
Cuts at theta cut = 0.2 deg are studied.
 +
 +
* Energy threshold for the analysis
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Energy thd srcind Crab int50 gcut06 tcut02.jpg|thumb|none|800px|gammaness cut>0.6 and theta cut< 0.2 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcind Crab int50 gcut07 tcut02.jpg|thumb|none|800px|gammaness cut>0.7 and theta cut< 0.2 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcind Crab int50 gcut08 tcut02.jpg|thumb|none|800px|gammaness cut>0.8 and theta cut< 0.2 deg, int>50 p.e.]] </li>
 +
<li style="display: inline-block;"> [[File:Energy thd srcind Crab int100 gcut09 tcut02.jpg|thumb|none|800px|gammaness cut>0.9 and theta cut<0.2 deg, int>100 p.e.]]</li>
 +
</ul></div>
 +
 +
I decide to use only the loose cuts in order to have a low energy threshold.
  
 
== Spectral results ==
 
== Spectral results ==
 
*; Source-dependent  
 
*; Source-dependent  
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab 20201120 loose.jpg|thumb|none|400px|Loose cuts, 20 November 2020. 2 bins per decade]] </li>
 +
<li style="display: inline-block;"> [[File:Crab 20201120 hard.jpg|thumb|none|400px|Hard cuts, 20 November 2020. 2 bins per decade]] </li>
 +
</ul></div>
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab 202109 loose.jpg|thumb|none|400px|Loose cuts, September 2021. 2 bins per decade]] </li>
 +
<li style="display: inline-block;"> [[File:Crab 202109 hard.jpg|thumb|none|400px|Hard cuts, September 2021. 2 bins per decade]] </li>
 +
</ul></div>
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab 202203 loose.jpg|thumb|none|400px|Loose cuts, March 2022. 2 bins per decade]] </li>
 +
<li style="display: inline-block;"> [[File:Crab 202203 hard.jpg|thumb|none|400px|Hard cuts, March 2022. 2 bins per decade]] </li>
 +
</ul></div>
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab all loose.jpg|thumb|none|600px|All runs. 4 bins per decade]] </li>
 +
</ul></div>
 +
 +
Conclusions: I decide to apply the following cuts: gammaness cut = 0.8, alpha cut = 15 deg and intensity = 50 p.e.
 +
 
*; Source-independent
 
*; Source-independent
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab srcind 20201120 loose.jpg|thumb|none|400px|Loose cuts, 20 November 2020]] </li>
 +
<li style="display: inline-block;"> [[File:Crab srcind 202109 loose.jpg|thumb|none|400px|Loose cuts, September 2021]] </li>
 +
<li style="display: inline-block;"> [[File:Crab srcind 202203 loose.jpg|thumb|none|400px|Loose cuts, March 2022]] </li>
 +
</ul></div>
 +
 +
<div><ul>
 +
<li style="display: inline-block;"> [[File:Crab srcind all loose.jpg|thumb|none|600px|All runs]] </li>
 +
</ul></div>

Latest revision as of 11:28, 29 April 2022

Overview[edit]

  • Source-dependent/independent analysis using LSTCHAIN v0.9.4 (MC production processed with LSTCHAIN v0.9.3, R0 to DL1 process by LSTOSA using LSTCHAIN v0.9.2)
  • Cleaning method tail cut 84 and dynamic cleaning. Also, the MC files are tuned to match the NSB in the FoV and the PSF.
  • Observations used in the analysis:

Monte Carlo information[edit]

  • Link to MC files used:
- Particle types: gamma off0.4deg, proton
- ZD (deg): 40
- AZ (deg): 180
- Other information: MC production through PR in lstMCpipe Github webpage (https://github.com/cta-observatory/lstmcpipe/pull/185)
- MC files are tuned to match the NSB in the FoV and the PSF.

DL1 data[edit]

  • Processed with LSTOSA using LSTCHAIN v0.9.2.
- Dynamic cleaning and tail cut cleaning with pedestal threshold applied by the the standard parameters using pipeline LSTOSA.
 "dynamic_cleaning": {
   "apply": true,
   "threshold": 267,
   "fraction_cleaning_intensity": 0.03
 }

Real data[edit]

  • Parameters for the tail cut cleaning with pedestal threshold
 "tailcuts_clean_with_pedestal_threshold": {
   "picture_thresh":8,
   "boundary_thresh":4,
   "sigma":2.5,
   "keep_isolated_pixels":false,
   "min_number_picture_neighbors":2,
   "use_only_main_island":false,
   "delta_time": 2
 }
  • Original DL1a files (processed up to DL1 by LSTOSA using LSTCHAIN v0.9.2)
/fefs/aswg/data/real/DL1/{}/v0.9/tailcut84/dl1_LST-1.RunXXXXX.XXXX.h5

MC data[edit]

  • Tunned NSB
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 0.6,
   "extra_bias_in_dim_pixels": 0.269,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 0.566,
  • Tunned PSF
   "increase_psf": true,
   "smeared_light_fraction": 0.175
  • Original DL1a files
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/{}/

Random forest[edit]

  • Dispnorm parameterisation
  • Source-dependent
- LSTCHAIN v0.9.4
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
- path:
/fefs/aswg/workspace/arnau.aguasca/Analysis/models/20200609_XXXX/srcdep/v0.9.4/
  • Source-independent
- Processed with lstMCpipe using LSTCHAIN v0.9.3.
- Path:
/fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/

DL1 to DL2 data[edit]

Real data[edit]

  • Source-dependent
- LSTCHAIN v0.9.4
- Config file:
//fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
  • Source-independent
- LSTCHAIN v0.9.4
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_standard_config_tailcut84_dispnorm_v0.9.3.json

MC data[edit]

  • Source-dependent
- LSTCHAIN v0.9.4
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
  • Source-independent
- Processed with lstMCpipe using LSTCHAIN v0.9.3.
- Config file:
/fefs/aswg/data/mc/analysis_logs/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/off0.4deg/lstchain_config_NSBtuning.json


DL3 data selection[edit]

  • Source-dependent
    - Quality cuts
 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2]
   }
 },
 "DL3FixedCuts": {
   "fixed_gh_cut": 0.8,
   "fixed_alpha_cut": 15,
   "allowed_tels": [1]
 }
  • Source-independent
    - Quality cuts

High-level analysis[edit]

  • Source-dependent
    - IRF: point-like, single-offset
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/XXXX
- Produced DL3 files
/fefs/aswg/workspace/arnau.aguasca/XXXX
  • Source-independent

Analysis results[edit]

  • High level analysis performed with gammapy-v0.19

Alpha plot[edit]

  • Filters in EventSelector applied:
 "EventSelector": {
   "filters": {
     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2]
   }
 }

See section Crab Check to know more about the justification of the applied cuts. Stacked alpha plots:

  • Alpha plot using runs in 20210809, 20210810 and 20210812.
  • Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812

Daily alpha plots:

  • Alpha plot using runs in 20210809.
  • Alpha plot using runs in 20210810.
  • Alpha plot using runs in 20210812.
  • Energy threshold for the analysis

Energy thd srcdep RSOph int50 gcut08 acut15.jpg

Theta plot[edit]

Stacked theta plots:

  • Theta2 plot using runs in 20210809, 20210810 and 20210812.
  • Cumulative signal to noise ratio and sqrt_ts. Black vertical line marks the the last run in 20210812

Daily theta plots:

  • Theta2 plot using runs in 20210809.
  • Theta2 plot using runs in 20210810.
  • Theta2 plot using runs in 20210812.


  • Energy threshold for the analysis

Energy thd srcind RSOph int50 gcut06 tcut02.jpg

Spectral results[edit]

  • Source-dependent
  • Stacked SED (Aug 9-12)
Daily SEDs
  • 20210809 SED
  • 20210810 SED
  • 20210812 SED
  • HESS and MAGIC comparison
  • HESS/Fermi and LST-1 SED.
  • MAGIC and LST-1 SED (Aug 9-12).
Daily MAGIC/LST-1 SEDs:
  • MAGIC and LST-1 SED (Aug 9).
  • MAGIC and LST-1 SED (Aug 10).
  • MAGIC and LST-1 SED (Aug 12).


  • Source-independent
  • Stacked SED (Aug 9-12)
Daily SEDs
  • 20210809 SED
  • 20210810 SED
  • 20210812 SED
  • HESS and MAGIC comparison
  • HESS/Fermi and LST-1 SED.
  • MAGIC and LST-1 SED (Aug 9-12).

No daily SED can be plotted because the SED from LST-1 using gc=0.6 and tc=0.2 deg are not reliable.

Light curve results[edit]

  • Source-dependent
  • HESS/Fermi and LST-1 light curve.
  • MAGIC and LST-1 light curve.
  • Source-independent
- Notice that since the daily SEDs with quality cuts (gc=0.6 and tc=0.2 deg) are not reliable, the daily flux are computed with the stacked SED using data between 20210809-12!!!
- No comparison with MAGIC light curve is done because the energy threshold using (gc=0.6 and tc=0.2 deg) is higher than 100 GeV.
  • HESS/Fermi and LST-1 light curve.


Crab Check[edit]

Run selection using Crab_data_selection_notebook: https://indico.cta-observatory.org/event/3984/contributions/32878/attachments/20998/29572/Crab_data_selection_notebook.pdf

Threshold for the rate30_mask is set to 4.3. I exclude run 6016 because I obtain a strange SED at low energies.

Data reduction Specifications[edit]

  • Source-dependent

SNr sig diff ghc ac.png

Cuts at alpha cut = 15 deg are studied.

  • Energy threshold for the analysis
  • gammaness cut>0.6 and alpha cut<15 deg, int>50 p.e.
  • gammaness cut>0.7 and alpha cut<15 deg, int>50 p.e.
  • gammaness cut>0.8 and alpha cut<15 deg, int>50 p.e.
  • gammaness cut>0.9 and alpha cut<15 deg, int>100 p.e.
  • Source-independent

SNr sig diff ghc tc.png

Cuts at theta cut = 0.2 deg are studied.

  • Energy threshold for the analysis
  • gammaness cut>0.6 and theta cut< 0.2 deg, int>50 p.e.
  • gammaness cut>0.7 and theta cut< 0.2 deg, int>50 p.e.
  • gammaness cut>0.8 and theta cut< 0.2 deg, int>50 p.e.
  • gammaness cut>0.9 and theta cut<0.2 deg, int>100 p.e.

I decide to use only the loose cuts in order to have a low energy threshold.

Spectral results[edit]

  • Source-dependent
  • Loose cuts, 20 November 2020. 2 bins per decade
  • Hard cuts, 20 November 2020. 2 bins per decade
  • Loose cuts, September 2021. 2 bins per decade
  • Hard cuts, September 2021. 2 bins per decade
  • Loose cuts, March 2022. 2 bins per decade
  • Hard cuts, March 2022. 2 bins per decade
  • All runs. 4 bins per decade

Conclusions: I decide to apply the following cuts: gammaness cut = 0.8, alpha cut = 15 deg and intensity = 50 p.e.

  • Source-independent
Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Analysis_A._Aguasca-Cabot_lstchain_v0.9.X&oldid=5629"