RS Oph Analysis Mab Bernardos v09X

From my_wiki
Revision as of 15:15, 29 March 2022 by Misabelber (talk | contribs)
Jump to: navigation, search

Overview

  • Source independent analysis using lstchain v0.9.4
  • Cleaning method: tailcut84 + pedestal cleaning + dynamic cleaning.
  • Observations: 20210809(Run5580-5584), 20210810(Run5630-5639), 20210812(Run 5696-5704)

Monte Carlo information

  • Particle types: point gamma, proton
  • Zenith 40deg, Azimuth 180 deg
  • Path to MC DL1 data:
   "/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/[gamma, gamma-diffuse, proton, electron]/zenith_40deg/south_pointing/0220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1"
  • Path to trained models:
  "/fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1"
  • Path to IRFs:
   "/fefs/aswg/data/mc/IRF/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/off0.4deg"
  • Path to config file:
  "/fefs/aswg/data/mc/IRF/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220307_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_src1/off0.4deg/lstchain_config_NSBtuning.json"
  • NSB is tuned to RS Oph observations:
    • Tuned 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,
    • Tuned PSF
   "increase_psf": true,
   "smeared_light_fraction": 0.175

Standard Cuts analysis

Event selection

Standard lstchain cuts were used:

     "intensity": [100, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2],
     "event_type": [32, 32]
   "global_gh_cut": 0.7,
   "gh_efficiency": 0.9,
   "global_alpha_cut": 10,
   "global_theta_cut": 0.2,
   "theta_containment": 0.68,
   "allowed_tels": [1]

First three nights (20210809, 20210810, 20210812) combined

  • Theta2 plots
    • Full energy range

Theta2 standard cuts.png Theta2 excess standard cuts.png

    • Per energy bin

Theta2 perbin standard cuts.png

  • SED

SED standard cuts.png

  • Best fit spectral parameters:

Spectral parameters standard cuts.png

  • Lightcurve

Lightcurve standard cuts.png

RS Oph cuts analysis

Event selection

Standard lstchain cuts were used:

     "intensity": [50, Infinity],
     "width": [0, Infinity],
     "length": [0, Infinity],
     "r": [0, 1],
     "wl": [0.1, 1],
     "leakage_intensity_width_2": [0, 0.2],
     "event_type": [32, 32]
   "global_gh_cut": 0.9,
   "gh_efficiency": 0.9,
   "global_alpha_cut": 10,
   "global_theta_cut": 0.2,
   "theta_containment": 0.68,
   "allowed_tels": [1]

First three nights (20210809, 20210810, 20210812) combined

  • Theta2 plots
    • Full energy range

Theta2 rsoph cuts.png Theta2 excess rsoph cuts.png

    • Per energy bin

Theta2 perbin rsoph cuts.png

  • SED

SED rsoph cuts.png

  • Best fit spectral parameters (fit range 0.05-40 TeV):

Spectral parameters rsoph cuts.jpg

  • Lightcurve

Lightcurve rsoph cuts.png

After moon break

  • 20210829 Theta2 plot

(Typo in number of wobbles [2] --> 4)

20210829 theta2 rsoph cuts.png 20210829 theta2 excess rsoph cuts.png

  • 20210830 Theta2 plot

(Typo in number of wobbles [2] --> 4)

20210830 theta2 rsoph cuts.png 20210830 theta2 excess rsoph cuts.png

  • 20210901 Theta2 plot

(Typo in number of wobbles [2] --> 4)

20210901 theta2 rsoph cuts.png 20210901 theta2 excess rsoph cuts.png

  • 20210902 Theta2 plot

(Typo in number of wobbles [2] --> 4)

20210902 theta2 rsoph cuts.png 20210902 theta2 excess rsoph cuts.png