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

From my_wiki
Jump to: navigation, search
Line 1: Line 1:
 
= Overview =
 
= Overview =
 
*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)
 
*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 10: Line 10:
 
:- AZ (deg): 180
 
:- AZ (deg): 180
 
:- Other information: MC production through PR in lstMCpipe Github webpage (https://github.com/cta-observatory/lstmcpipe/pull/185)
 
:- 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 ==  
 
== DL1 data ==  
Line 58: Line 59:
 
== Random forest ==
 
== Random forest ==
  
* Dispnorm
+
* Dispnorm parameterisation
 +
 
 
*; Source-dependent  
 
*; Source-dependent  
* LSTCHAIN v0.9.4
+
:- LSTCHAIN v0.9.4
 
:- Config file:
 
:- Config file:
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
Line 67: Line 69:
  
 
*; Source-independent  
 
*; Source-independent  
* Processed with lstMCpipe using LSTCHAIN v0.9.3.
+
:- Processed with lstMCpipe using LSTCHAIN v0.9.3.
 
:- Path:
 
:- 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/
 
  /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/
Line 75: Line 77:
 
=== Real data ===
 
=== Real data ===
 
*; Source-dependent  
 
*; Source-dependent  
* LSTCHAIN v0.9.4
+
:- LSTCHAIN v0.9.4
::- config file:
+
:- Config file:
 
  //fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
 
  //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
* LSTCHAIN v0.9.4
+
:- LSTCHAIN v0.9.4
* Config file:
+
:- Config file:
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_standard_config_tailcut84_dispnorm_v0.9.3.json
 
  /fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_standard_config_tailcut84_dispnorm_v0.9.3.json
  
Line 87: Line 89:
  
 
*; Source-dependent  
 
*; Source-dependent  
* LSTCHAIN v0.9.4
+
:- LSTCHAIN v0.9.4
* Config file:
+
:- Config file:
  //fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.4/lstchain_src_dep_tailcut84_config_dispnorm_v0.9.4.json
+
  /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
* Processed with lstMCpipe using LSTCHAIN v0.9.3.
+
:- Processed with lstMCpipe using LSTCHAIN v0.9.3.
* Config file:
+
:- 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
 
  /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
  
Line 102: 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 122: 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

Revision as of 07:33, 24 March 2022

Overview

  • 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

  • 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

  • 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

  • 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

  • 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

  • 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

Real 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
- 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

  • 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

  • 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

  • High level analysis performed with gammapy-v0.19

Alpha plot

  • Filters in EventSelector applied:

Theta plot

Spectral results

  • Source-dependent
    Source-independent

Light curve results

  • Source-dependent
    Source-independent

Crab Check

Data reduction Specifications

  • Source-dependent
    Source-independent

Spectral results

  • Source-dependent
    Source-independent