Analysis A. Aguasca-Cabot lstchain v0.9.X

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:

 "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]
   }
 }

Spectral results

• Source-dependent

  Source-independent

Light curve results

• Source-dependent

  Source-independent

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

• Source-dependent

Cuts at alpha cut = 15 deg are studied.

• Energy threshold for the analysis

• 

  gammaness cut>0.6 and alpha cut<15 deg, int>50
• 

  gammaness cut>0.7 and alpha cut<15 deg, int>50
• 

  gammaness cut>0.8 and alpha cut<15 deg, int>50
• 

  gammaness cut>0.9 and alpha cut<15 deg, int>100

• Source-independent

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
• 

  gammaness cut>0.7 and theta cut< 0.2 deg, int>50
• 

  gammaness cut>0.8 and theta cut< 0.2 deg, int>50
• 

  gammaness cut>0.9 and theta cut<0.2 deg, int>100

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

Spectral results

• Source-dependent

• 

  Loose cuts, 20 November 2020. 2 bins per decade
• 

  Loose 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
• 

  All runs. 4 bins per decade

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

• Source-independent

• 

  Loose cuts, 20 November 2020
• 

  Loose cuts, September 2021
• 

  Loose cuts, March 2022

• 

  All runs
• 

  All runs