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

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:

File:Alpha plot 20210809.png
Alpha plot using runs in 20210809.
File:Alpha plot 20210810.png
Alpha plot using runs in 20210810.
File:Alpha plot 20210812.png
Alpha plot using runs in 20210812.

Energy threshold for the analysis

Theta plot

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:

File:Theta2 plot 3off gc06 tc02 20210809.png
Theta2 plot using runs in 20210809.
File:Theta2 plot 3off gc06 tc02 20210810.png
Theta2 plot using runs in 20210810.
File:Theta2 plot 3off gc06 tc02 20210812.png
Theta2 plot using runs in 20210812.

Energy threshold for the analysis

Spectral results

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

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.