GRB221009A src-ind analysis dark conditions

From my_wiki
Revision as of 12:27, 3 March 2023 by AAguascaCabot (talk | contribs)
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 A. Aguasca-Cabot (Universitat de Barcelona - arnau.aguasca@fqa.ub.edu)


General information

This wiki web page shows the analysis of GRB221009A data recorded in dark conditions (run_id: 9625-9960) and in mild moon conditions (run_id: 9629 and 9630) using the source-independent analysis.

For the analysis of GRB221009A for runs in moon conditions go to the following page: GRB221009A_src-ind_analysis_moon

The wiki page with the summary of the current results of the src-ind moon analysis (as of February 8th) for the discussion with LST-Reco conveners is the following: [1]

  • Slides presented at the GRB221009A meetings about the source-independent dark/mild-moon analysis:
Start and Dark analysis I
Dark analysis II
Dark analysis III
Dark analysis IV and mild moon
Dark analysis V and mild moon


Dark analysis

Monte Carlo information

  • Link to MC files used: /fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/*
    • Particle types: The standard ones for the AllSky MC production.
    • DEC Band (deg): dec_2276

The MC production for the Crab is used because the NSB tuning parameters at the DL1 level are very similar between them:

  • Pixel charge distribution for real cosmics (blue, run 9672), standard MC gammas (orange), tuned MC gammas to Crab (green) and to GRB221009A (red) using the DL1 NSB parameters. Plot with y-axis in log scale.

DL1 parameters obtained from GRB221009A runs (the NSB tuning parameters for the Crab can be found in the following subsection):

Runs: 9672, 9802, 9862 and 9890
{
    "number_of_files_used": 46,
    "std": [
    0.181,
    0.077,
    0.209
]
}{
    "increase_nsb": true,
    "extra_noise_in_dim_pixels": 1.782,
    "extra_bias_in_dim_pixels": 0.601,
    "transition_charge": 8,
    "extra_noise_in_bright_pixels": 2.25
}

DL1 data

No reprocessing of the real DL1a data was required, only NSB tuning was added to the MC files. DL1a files produced by LSTOSA (lstchain v0.9.6?) and lstMCpipe (lstMCpipe v0.9.0 using lstchain v0.9.9).

  • original DL1a files
/fefs/aswg/data/real/DL1/20221010/v0.9/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/
  • lstchain v0.9.9
  • real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning.
 /fefs/aswg/data/real/DL1/20221015/v0.9/tailcut84/log/lstchain_config_v0.9.4_from2022onwards.json
  • lstMCpipe-v0.9.0
  • lstchain v0.9.9
  • MC: tailcut8-4, dynamic cleaning, NSB tuning
   "image_modifier": {
       "increase_nsb": true,
       "extra_noise_in_dim_pixels": 1.62,
       "extra_bias_in_dim_pixels": 0.655,
       "transition_charge": 8,
       "extra_noise_in_bright_pixels": 2.08,
       "increase_psf": false,
       "smeared_light_fraction": 0
   },
Config file: https://github.com/cta-observatory/lstmcpipe/blob/master/production_configs/20221027_v0.9.9_crab_tuned/lstchain_config_2022-10-27.json

Random forest

RFs were produced with lstMCpipe

  • lstchain-v0.9.9
  • lstMCpipe-v0.9.0
  • source-indep (the standard one for the AllSky MC production)
/fefs/aswg/data/models/AllSky/20221027_v0.9.9_crab_tuned/dec_2276/

DL2 data

  • lstchain-v0.9.9
  • source-indep
/fefs/aswg/data/mc/DL2/AllSky/20221027_v0.9.9_crab_tuned/TestingDataset/dec_2276/node_theta_*
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL2/202210*/srcind/v0.9.9/tailcut84_dynclg_tNSBtoCrabPerfPaper/AllSky_20221027_v0.9.9_crab_tuned_dec_2276/dl2_*
  • config file
 /fefs/aswg/data/real/DL1/20221015/v0.9/tailcut84/log/lstchain_config_v0.9.4_from2022onwards.json

DL2 data selection to produce theta2 plots

Information about the quality cuts to produce the theta2 plots.

Example

  • intensity > 50
  • r: [0, 1 ]
  • wl: [0.01, 1 ]
  • leakage_intensity_width_2: [0, 0.2 ]
  • event_type: [32, 32]
  • source-indep
    • fixed_gh_cut: 0.7 (standard cut)
    • fixed_theta_cut: 0.2 (standard cut)


Crab check

Mild moon analysis

Monte Carlo information

  • Link to MC files used: /fefs/aswg/data/mc/DL1/AllSky/20221215_v0.9.12_base_prod/*
    • Particle types: The standard ones for the AllSky MC production.
    • DEC Band (deg): dec_2276

DL1 data

No reprocessing of the real DL1a data was required, only NSB tuning was added to the MC files. DL1a files produced by LSTOSA (lstchain v0.9.6?) and lstMCpipe (lstMCpipe v0.9.0 using lstchain v0.9.9).

  • original DL1a files
/fefs/aswg/data/real/DL1/20221010/v0.9/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/
  • lstchain v0.9.9
  • real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning.
 /fefs/aswg/data/real/DL1/20221015/v0.9/tailcut84/log/lstchain_config_v0.9.4_from2022onwards.json
  • lstMCpipe-v0.9.0
  • lstchain v0.9.9
  • MC: tailcut8-4, dynamic cleaning, NSB tuning
   "image_modifier": {
       "increase_nsb": true,
       "extra_noise_in_dim_pixels": 1.62,
       "extra_bias_in_dim_pixels": 0.655,
       "transition_charge": 8,
       "extra_noise_in_bright_pixels": 2.08,
       "increase_psf": false,
       "smeared_light_fraction": 0
   },
Config file: https://github.com/cta-observatory/lstmcpipe/blob/master/production_configs/20221027_v0.9.9_crab_tuned/lstchain_config_2022-10-27.json

Random forest

RFs were produced with lstMCpipe

  • lstchain-v0.9.9
  • lstMCpipe-v0.9.0
  • source-indep (the standard one for the AllSky MC production)
/fefs/aswg/data/models/AllSky/20221027_v0.9.9_crab_tuned/dec_2276/

DL2 data

  • lstchain-v0.9.9
  • source-indep
/fefs/aswg/data/mc/DL2/AllSky/20221027_v0.9.9_crab_tuned/TestingDataset/dec_2276/node_theta_*
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL2/202210*/srcind/v0.9.9/tailcut84_dynclg_tNSBtoCrabPerfPaper/AllSky_20221027_v0.9.9_crab_tuned_dec_2276/dl2_*
  • config file
 /fefs/aswg/data/real/DL1/20221015/v0.9/tailcut84/log/lstchain_config_v0.9.4_from2022onwards.json

DL2 data selection to produce theta2 plots

Information about the quality cuts to produce the theta2 plots.

Example

  • intensity > 50
  • r: [0, 1 ]
  • wl: [0.01, 1 ]
  • leakage_intensity_width_2: [0, 0.2 ]
  • event_type: [32, 32]
  • source-indep
    • fixed_gh_cut: 0.7 (standard cut)
    • fixed_theta_cut: 0.2 (standard cut)

Analysis Results

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

Theta2 plot

Significance map

Excess map

Spectral results