RS Oph Analysis Mab Bernardos v09X

From my_wiki
Jump to: navigation, search

Overview[edit]

  • 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[edit]

  • 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[edit]

Event selection[edit]

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[edit]

Theta2 plots[edit]

Full energy range[edit]

Theta2 standard cuts.png Theta2 excess standard cuts.png

Per energy bin[edit]

Theta2 perbin standard cuts.png

SED[edit]

SED standard cuts.png

Best fit spectral parameters (fit range 0.01-40 TeV):[edit]

Spectral parameters standard cuts.png

Lightcurve[edit]

Lightcurve standard cuts.png

RS Oph cuts analysis[edit]

Event selection[edit]

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[edit]

Theta2 plots[edit]

Full energy range[edit]

Theta2 rsoph cuts.png Theta2 excess rsoph cuts.png

Per energy bin[edit]

Theta2 perbin rsoph cuts.png

SED[edit]

SED rsoph cuts.png

Best fit spectral parameters (fit range 0.05-40 TeV):[edit]

Spectral parameters rsoph cuts.jpg

Lightcurve[edit]

Lightcurve rsoph cuts.png

After moon break[edit]

20210829 Theta2 plot[edit]

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

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

20210830 Theta2 plot[edit]

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

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

20210901 Theta2 plot[edit]

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

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

20210902 Theta2 plot[edit]

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

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