GRB221009A srcdep analysis

From my_wiki
Revision as of 10:01, 20 February 2023 by Kenta.T (talk | contribs) (Created page with " Back to the Data analysis page Go back to Transient Working Group. Go Back to [https://www.lst1.iac.es/wiki/index.php/Transient_observation_logbook G...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Back to the Data analysis page

Go back to Transient Working Group.

Go Back to Gamma-Ray Bursts (GRBs).

Go back to GRB221009A.

  • Analysis by Kenta Terauchi (Kyoto University - terauchi.kenta.74s@st.kyoto-u.ac.jp)

General information

Superplots

  • Superplots on 2022-10-10. Made by Seiya


NSB level

  • Median of pedestal std on 2022-10-10
  • Median of pedestal std on 2022-10-12


Calibrations

  • Pedestal std distribution of run9602(1st subrun)
  • Pedestal std distribution of run9602(1st subrun). Calibration was processed with new calibration file made by Franca
  • Pedestal std distribution in the camera frame (run9602;1st subrun).
  • Pedestal std distribution in the camera frame (run9602;1st subrun). Calibration was processed with new calibration file made by Franca

Image cleaning

  • Cleaning method: tailcut cleaning w/ pedestal cleaning, dynamic cleaning
  • Cleaning threshold: (picture threshold, boundary threshold)

The following threshold values are obtained so that ~5% of pedestal events survive after image cleaning. For now, the ratio of picture threshold to boundary threshold are fixed to 2:1 (same as default ratio, 8:4).

<2022_10_10>
run9602 - 9603: (24, 12)
run9604: (26, 13)
run9605 - 9607: (28, 14)
<2022_10_12>
run9613: (26, 13)
run9614: (18,9)
run9615 - 9616: (20, 10)
run9617: (22, 11)
  • Sigma value for pedestal cleaning: 2.5 (default)


NSB tuning of MC data

I produced run-wise MC taking into account the NSB level, cleaning threshold, and Zd/Az range of each run.

Monte Carlo information

  • Link to MC files used (base): /fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/*
  • Particle types: standard particle types for AllSky MC prodution
  • DEC band: dec_2276

MC productions

  • lstchain v0.9.10
  • lstMCpipe v0.10
  • NSB tuning method: Adding NSB noise to DL1 images. (Adding NSB noise to R0 waveform is VERY time consuming! It makes not so much sense to tune waveform for now. Maybe we could produce more dedicated MC in later stage before publication)


Date: 2022-10-10

  • run9602
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 28.061,
   "extra_bias_in_dim_pixels": 6.6885,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 69.279,
   "increase_psf": false,
   "smeared_light_fraction": 0
 },
  • run9603
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 29.361,
   "extra_bias_in_dim_pixels": 6.632,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 78.425,
   "increase_psf": false,
   "smeared_light_fraction": 0
 }
  • run9604
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 31.232,
   "extra_bias_in_dim_pixels": 6.977,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 85.884,
   "increase_psf": false,
   "smeared_light_fraction": 0
 }
  • run9605
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 33.304,
   "extra_bias_in_dim_pixels": 7.233,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 94.252,
   "increase_psf": false,
   "smeared_light_fraction": 0
 }
  • run9606
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 33.860,
   "extra_bias_in_dim_pixels": 7.2195,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 100.024,
   "increase_psf": false,
   "smeared_light_fraction": 0
 }
  • run9607
"image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 37.096,
   "extra_bias_in_dim_pixels": 7.972,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 107.688,
   "increase_psf": false,
   "smeared_light_fraction": 0
 }

Date: 2022-10-12

DL1 data

  • original DL1a files
/fefs/aswg/workspace/franca.cassol/data/real/DL1/20221010/tailcut84/
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/*
  • Produced DL1b data (2022-10-10)
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/source/GRB221009A/2022_10_10/DL1_new_calib_tuned/dl1_LST-1.Run0960[2-7].*.h5
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/NSB_tuning/AllSky_nsb_tuned/run960[2-7]/*


Random forest

  • lstchain-0.9.10, lstchain-0.9.14-dev (branch: extra_srcdep_params)


  • source-dep (diffuse gamma, src_r<1deg)
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/RF/AllSky_nsb_tuned/run960[2-7]/trained_models/*
  • source-dep (diffuse gamma, src_r<1deg, w/ extra srcdep parameters)
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/RF/AllSky_nsb_tuned/run960[2-7]/trained_models/*
  • source-dep (diffuse gamma, 0.2<src_r<0.6deg)
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/RF/AllSky_nsb_tuned/run960[2-7]/trained_models/*
  • source-dep (diffuse gamma, src_r<1deg)
/fefs/aswg/workspace/kenta.terauchi/Work/lst-analysis/RF/AllSky_nsb_tuned/run960[2-7]/trained_models/*

DL2 data

Information about your DL2 data and settings such as: specifc versions of lstchain, specfic .json version used.

Example:

  • lstchain-0.7.6.dev242+g2086cb9
  • source-indep
/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcindep/DL2/data/
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/srcindep/DL2/data/dl2_LST-1.Run0XXXX_merged.h5
  • source-dep
/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcdep/DL2/data/ 
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/srcdep/DL2/data/dl2_LST-1.Run0XXXX_merged.h5

DL3 data selection

Information about your DL3 data selection.

Example

  • intensity > 50
  • r: [0, 1 ]
  • wl: [0.1, 1 ]
  • leakage_intensity_width_2: [0, 0.2 ]
  • source-indep
    • fixed_gh_cut: 0.3
    • fixed_theta_cut: 0.2
  • source-dep
    • fixed_gh_cut: 0.7
    • fixed_alpha_cut: 10

High-level analysis

Please put any information about the production of higher level analysis here.

Example

  • lstchain to generate source-dep IRF, DL3
  • Science Tool: gammapy 0.18.2
  • point-like IRF, 1D analysis

Analysis Results

Please place higher-level analysis results (Spectra, SkyMaps, Lightcurves, etc) here.

Theta2 plot

Significance map

Excess map

Spectral results