RS Oph Analysis v0.9.x AllSky

From my_wiki
Jump to: navigation, search

Overview

  • Source-independent analysis using LSTCHAIN v0.9.10 (DL1->DL2) and LSTCHAIN v0.9.14.dev79+g6a03d48.d20230508 (DL2 to DL3 stage (irf interpolation)). The MC production processed with lstMCpipe v0.9.0 using LSTCHAIN v0.9.9.
  • Source-dependent analysis using LSTCHAIN v0.9.14.dev139+g90d0450 (DL1->DL2) and LSTCHAIN v0.9.14.dev79+g6a03d48.d20230508 (DL2 to DL3 stage (irf interpolation)). The MC production processed with lstMCpipe v0.9.1.dev14+gfb9c7c using LSTCHAIN v0.9.14.dev139+g90d0450.


  • Cleaning method tail cut 84 and dynamic cleaning. Also, the MC files are tuned to match the NSB in the FoV.
  • Observations used in the analysis: 5580->5584, 5630->5639, 5696->5704, 5890->5893, 5912->5917, 5942->5945, 5970->5973.


Monte Carlo information

  • AllSky MC prodution.
  • Link to MC files used (srcind): /fefs/aswg/data/mc/DL1/AllSky/20221108_galsource_min_413_tuned_nsb/...
- Particle types: gamma diffuse and protons
- MC prod (deg): 20221108_galsource_min_413_tuned_nsb
- Dec band (deg): dec_min_431
- Other information: MC production through PR in lstMCpipe Github webpage (https://github.com/cta-observatory/lstmcpipe/tree/master/production_configs/20221108_galsource_min_413_tuned_nsb)
- MC files are tuned to match the NSB in the FoV.
  • Link to MC files used (srcdep): /fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230301_srcdep_galsource_min_413_tuned_nsb/...
- Particle types: gamma diffuse and protons
- MC prod (deg): 20230301_srcdep_galsource_min_413_tuned_nsb
- Dec band (deg): dec_min_431
- Other information: MC processed by Arnau (config files in path: /fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230301_srcdep_galsource_min_413_tuned_nsb/)
- MC files are tuned to match the NSB in the FoV.


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

  • Tuned NSB
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 0.84,
   "extra_bias_in_dim_pixels": 0.326,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 1.064,


  • Original DL1a files (srcind)
/fefs/aswg/data/mc/DL1/AllSky/20221108_galsource_min_413_tuned_nsb/{TrainingDataset,TestingDataset}/dec_min_413/{particle}/node...
  • Original DL1a files (srcdep)
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230301_srcdep_galsource_min_413_tuned_nsb/{TrainingDataset,TestingDataset}/dec_min_413/{particle}/node...


Random forest

  • Disp_norm parameterisation
  • Source-independent
- Processed with lstMCpipe using LSTCHAIN v0.9.9.
- Path:
/fefs/aswg/data/mc/models/AllSky/20221108_galsource_min_413_tuned_nsb/dec_min_413/
  • Source-dependent
- LSTCHAIN v0.9.14.dev139+g90d0450
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230301_srcdep_galsource_min_413_tuned_nsb/lstchain_config_2023-03-01.json
- path:
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230301_srcdep_galsource_min_413_tuned_nsb/dec_min_413/


DL1 to DL2 data

Real data

  • Source-independent
- LSTCHAIN v0.9.10
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/standard_real_config_lstchain_belong_v0.10.0.json
  • Source-dependent
- LSTCHAIN v0.9.14.dev139+g90d0450
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/RSOph/analysis_v0.9.13/standard_real_srcdep_config_lstchain_belong_extra_srcdep_params.json


MC data

  • Source-independent
- Processed with lstMCpipe using LSTCHAIN v0.9.9.
- Config file:
https://github.com/cta-observatory/lstmcpipe/blob/master/production_configs/20221108_galsource_min_413_tuned_nsb/lstchain_config_2022-11-08.json
  • Source-dependent
- LSTCHAIN v0.9.14.dev139+g90d0450
- Config file:
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230301_srcdep_galsource_min_413_tuned_nsb/lstchain_config_2023-03-01.json

Quality cuts optimization for detection (srcind)

Since the AllSky MC is a tailored production than the Fixed MC production, we cannot use high-zenith Crab observations to find the best cuts for detection. To do so, I found the best cuts (gammaness and theta) that minimise the sensitivity of MC gammas re-weightening the slope the the spectrum reported by MAGIC (Gamma = -3.7) and using OFF data from real OFF data (3 reflected OFF regions from runs 5890,5891,5892,5893, 5912,5913,5914,5915,5916,5917, 5942, 5943, 5944, 5945, 5970, 5971, 5972, 5973). I used a grid of gammaness and theta2 values. The best cuts are gammaness = 0.9 and theta2 = 0.01 deg2.

Theta2 plot

  • Right after the outburst (run id <= 5704)

File:Theta2 plot RSOph.pdf

  • After moon break (run id >= 5890)

File:Theta2 plot RSOph afterMoonBreak.pdf

DL3 data selection

  • - Quality cuts
 "EventSelector": {
   "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]
   }
 },
 "DL3FixedCuts": {
   "gh_efficiency": 0.7,
   "theta_containment": 0.7,
   "alpha_containment": 0.7,
   "allowed_tels": [1]
 }

High-level analysis

  • Source-independent
    - IRF: point-like, single-offset, node interpolation.
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/DL3IRF-generated-configs/config_thetaChange_int50_leak0.2_gh0.9_geff0.7_th0.2_tc0.7_a15_ac0.68.json
- Produced DL3 files
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/LST/real/DL3/RSOph/srcind/v0.9.14.dev79+g6a03d48.d20230508/tailcut84_dynclg_tNSBtoRSOph/AllSky_20221108_galsource_min_413_tuned_nsb_dec_min_413/int_50_lk_0.2_gh_dyn_0.7_th_dyn_0.7/interp_node/
  • Source-dependent
    - IRF: point-like, single-offset, nearest node.
/fefs/aswg/workspace/arnau.aguasca/scripts/_configs/DL3IRF-generated-configs/config_thetaChange_int50_leak0.2_gh0.9_geff0.5_th0.2_tc0.7_a15_ac0.7.json
- Produced DL3 files
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/LST/real/DL3/RSOph/nearest_node/srcdep/v0.9.14.dev79+g6a03d48.d20230508/tailcut84_dynclg_tNSBtoRSOph/AllSky_20230301_srcdep_galsource_min_413_tuned_nsb_dec_min_413/int50_lk0.2_geff0.5_aeff0.7/


Analysis results src ind

  • The srcind high-level results will be used in the paper publication.
  • High-level analysis performed with gammapy-v1.0. Fitting process using the joint method in gammapy for the srcind results.
  • Spectral Energy distribution
  • 20210809
File:LST-1MAGICHESS SED RSOph 20210809 interp node srcind v0.9.14dev AllSkyMC.pdf
  • 20210810
File:LST-1MAGICHESS SED RSOph 20210810 interp node srcind v0.9.14dev AllSkyMC.pdf
  • 20210812
File:LST-1MAGICHESS SED RSOph 20210812 interp node srcind v0.9.14dev AllSkyMC.pdf
  • All data from Aug 9, 10 and 12
File:LST-1MAGIC SED RSOph 09to12 interp node srcind v0.9.14dev AllSkyMC.pdf
  • Light curve
  • Above E>100 GeV (MAGIC comparison)
File:LST-1MAGIC lightcurve RSOph interp node srcind v0.9.14dev AllSkyMC.pdf
  • Above E>250 GeV (HESS comparison)
File:LST-1HESS lightcurve RSOph interp node srcind v0.9.14dev AllSkyMC.pdf


Analysis results src dep

  • Comparison with srcind and srcdep
  • 20210809
File:Plot LST-1MAGICHESS SED RSOph final srcind&srcdep comparison 20210809.pdf
  • 20210810
File:Plot LST-1MAGICHESS SED RSOph final srcind&srcdep comparison 20210810.pdf
  • 20210812
File:Plot LST-1MAGICHESS SED RSOph final srcind&srcdep comparison 20210812.pdf
  • All data from Aug 9, 10 and 12
File:Plot LST-1MAGICHESS SED RSOph final srcind&srcdep comparison 09to12.pdf
  • Light curve
  • Above E>100 GeV (MAGIC comparison)
File:LST-1MAGIC lightcurve RSOph interp node srcind v0.9.14dev AllSkyMC.pdf
  • Above E>250 GeV (HESS comparison)
File:LST-1HESS lightcurve RSOph interp node srcind v0.9.14dev AllSkyMC.pdf