Difference between revisions of "GRB221009A src-ind analysis moon"

From my_wiki
Jump to: navigation, search
m (AAguascaCabot moved page GRB221009A src-ind analysis to GRB221009A src-ind analysis moon: more accurate title)
Line 289: Line 289:
  
 
=== Spectral results ===
 
=== Spectral results ===
−
 
−
=Dark/mild moon analysis=
 

Revision as of 10:02, 14 February 2023

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

  • We have two different analysis for GRB221009A. First, the analysis with strong moon conditions, which occurred during October 10th and 12th (see section LST Observations in GRB221009A). Second, the analysis with dark/mild moon conditions, which data was taken between October 15th to October 28th. However, notice that data taking in the position of GRB2210009A was extended up to November. Thus, these observations contains fantastic background (probably).
  • The analysis is done using lstchain v0.9.9 or lstchain v0.9.10. The version release between them does not affect dramatically the offsite analysis.


Strong moon analysis

Since the observations were done in strong moon conditions, the HV gain was reduced to 70%/50% of the nominal value (see Elog for more information). By the time of the data taking was done and beginning of 2023, this analysis is not standard.

Problemes

Due to the no standard observation conditions, several aspects are affected. Are the following (probably more problems are missing): (Some of the things listed here come from an internal communication with Abelardo, thanks!)

  • Pedestal pixels: several Pixels can be tagged as unusable using the standard calibration config. Thus will produce several pixels are not calibrated.
  • F-factor: Using lower HV, probably the F-factor will change since for a signal in the PMT, there will be more fluctuations due to the moon conditions. Do we have the characterisation of the PMTs in those conditions? Meaning, we do not have the gain fluctuations curve for this conditions. Then, due to the wrong F-factor, we will get wrong conversion factors for data. This will act as an additional noise with a F-factor-style noise (gain fluctuations) that we should consider in the MC simulations.
  • MC-real data match
  • To match the correct gain response of the PMTs for simulation, the gain fluctuations curve for these conditions must be known. So, the simulations will not follow the real response if we do not have this curve. This will induce an addition noise with a "F-factor-style noise" (gain fluctuations) between MC and real data.
  • The high noise conditions will produce that we have high fluctuations in the PMTs. As the the standard cleaning method is the tail cut cleaning with pedestal std condition, several pixels will have increased pixel picture/boundary threshold because of the high std of the pedestal pixels. Then, this will produce an inhomogeneous camera response with changing image cleaning values along the camera. On the other hand, the MC cleaning will not have this kind of response because the pedestal std conditions is not used in the MC data because this condition is intended for stars (which are not considered in the current MC). Then, there will be a mismatch between data and MC.
  • NSB tuning: The standard analysis procedure matches the NSB tuning of the data (Moon, stars in the FoV, ...) at the DL1 stage. The problem of doing so is that the added NSB will not be as similar as the actual NSB added in the PMT waveform. Also, if we do not assess and correct the additional F-factor-style noise mentioned above, it will contribute as a noise between MC and data to match by the script, but due to the different behaviour, it will not be correctly filtered out using the current NSB tuning. (And maybe even the tuned MC with that NSB parameters will not be as good as the output if we would have only photon noise)


  • How to face these problems?
  • Pedestal pixels: I think Franca solved this issue increasing the nominal range of pedestal pixels to be tagged as usable. Not sure...
  • F-factor: One could consider applying an empirical correction for such high-NSB data.
  • MC-real data match (1): We need the gain fluctuations curve.
  • MC-real data match (2): the levels of tail cut cleaning should be high enough such that only a handful of pixels get raised thresholds, and therefore data and MC are more or less similar.
  • MC-real data match (NSB): use the NSB tuning at the waveform level, or make specific MC simulations matching directly the NSB of the runs.


General information

Based on the problems we face in this analysis. We following procedure is used (detailed information en each section):

  • 1. Calibrations: First, Franca reprocessed the DL1 data (taken on Oct. 10th) to have a better calibration and more usable pixels.
  • 2. Image cleaning: The runs at DL1 stage where analysed to find the best image cleaning for each run that satisfied that only about 4% of the pedestal pixels will raise their picture thresholds due to the pedestal std dev condition, and not more than 5% pedestals survive the cleaning. The former condition is the restrictive one in work. The std value for the image cleaning with pedestal std is set to 2.5 (standard one). Next, we apply the best cleaning to DL1 real data.
  • 3. NSB tuning: We found the noise required to apply to the MC data in order to match the NSB conditions of the real data. Both the NSB tuning at the waveform level and at DL1 stage were used (using the tail cut cleaning values). However, for the former, the tuned MC gammas have a strange pixel charge distribution and the computation time is really high (between 1-2 days or 12h-1day if we consider lstchain PR 1066. Nevertheless, if we consider the latter code, all MC tuned files show the strange pixel charge distribution). We decided to proceed with the NSB tuning at DL1 stage
  • 4. MC productions: based on the steps 2 and 3, we decided the produce different MC production to analyse the data.

Calibrations

You can see here the difference between the DL1 files produced with the standard calibration configuration (standard when the data was taken, v0.9) and the DL1 file recalibrated by Franca (specific paths can be found later) for a run taken in Oct. 10th.

  • Pedestal pixels tagged as usable for run /fefs/aswg/data/real/DL1/20221010/v0.9/tailcut84/dl1_LST-1.Run09602.0005.h5. 0 = unusable, 1 = usable.
  • Pedestal pixels tagged as usable for run /fefs/aswg/workspace/franca.cassol/data/real/DL1/20221010_test1/tailcut84/dl1_LST-1.Run09602.0005.h5. 0 = unusable, 1 = usable.

Due to high number of unusable pixels in the standard files (v0.9), the resulting picture threshold is really high because the picture threshold set to the unusable pixels is the maximum pedestal threshold value of the pedestal pixel in the camera.

On the other hand, few unusable pixels are tagged in the standard files for runs recorded in Oct. 12th.

  • Pedestal pixels tagged as usable for run /fefs/aswg/data/real/DL1/20221012/v0.9/tailcut84/dl1_LST-1.Run09613.0000.h5. 0 = unusable, 1 = usable.
  • Pedestal pixels tagged as usable for run /fefs/aswg/data/real/DL1/20221012/v0.9/tailcut84/dl1_LST-1.Run09617.0000.h5. 0 = unusable, 1 = usable.


Image cleaning

You can see how the picture threshold value evolve for the different runs in Oct. 10th and 12th. The picture threshold increase as we increase the run_id for a given day because the moon was rising. Notice that the picture threshold for run 9613 is higher than the rest of the runs in Oct. 12th because the HV was reduced to 50%, not to 70% of the nominal value as in the rest of the runs.

  • Run-wise picture threshold distribution for runs taken in Oct. 10th.
  • Run-wise picture threshold distribution for runs taken in Oct. 12th.

NSB tuning

Here you can see the pixel charge distribution between real data and tuned MC data using NSB tuning at the waveform (right and wrong tuning?). The NSB tuning at the DL1 stage is added to the right plot.

  • Pixel charge dist for real data (orange) and tuned MC data for two runs, one with strange distribution (green) and another closer to the real data distribution (blue).
  • Pixel charge dist for real data (orange) and tuned MC data using the tuning at the waveform level (green) and at DL1 stage (blue).

It is clear that MC data and real data do not fully match in the pixel charge distribution.

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

Based on the image cleaning and NSB tuning. We decided to produce the following MC productions.

  • lstchain v0.9.10
  • lstMCpipe v0.10
  • Date: 20221010
  • MC production 1: 20230205_GRB221009A_moon_strong_ph1_NSBDL1
  • Runs to analyse (run_id): 9602-9604
  • Tail cut cleaning (picture thd., boundary thd.): 30, 15
  • NSB tuning (from mean parameters values using run 9603 subruns):
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 30.367,
   "extra_bias_in_dim_pixels": 6.881,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 78.926,
  • Configuration files
fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph1_NSBDL1/lstchain_config_MC_Franca_test1_tailcut3015.json
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph1_NSBDL1/lstmcpipe_config_2023-01-20_PathConfigAllSkyFullDL1ab.yaml
  • MC production 2: 20230205_GRB221009A_moon_strong_ph2_NSBDL1
  • Runs to analyse (run_id): 9605-9607
  • Tail cut cleaning (picture thd., boundary thd.): 34, 17
  • NSB tuning (from mean parameters values using run 9605 subruns):
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 36.21,
   "extra_bias_in_dim_pixels": 7.828,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 100.854,
  • Configuration files
fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph2_NSBDL1/lstchain_config_MC_Franca_test1_tailcut3417.json
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph2_NSBDL1/lstmcpipe_config_2023-01-20_PathConfigAllSkyFullDL1ab.yaml


  • Date: 20221012
  • MC production 1: 20230205_GRB221009A_moon_strong_ph3_NSBDL1
  • Runs to analyse (run_id): 9613
  • Tail cut cleaning (picture thd., boundary thd.): 32, 16
  • NSB tuning (from mean parameters values using run 9613 subruns):
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 29.754,
   "extra_bias_in_dim_pixels": 6.461,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 23.584,
  • Configuration files
fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph3_NSBDL1/lstchain_config_2023-01-20.json
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph3_NSBDL1/lstmcpipe_config_2023-01-20_PathConfigAllSkyFullDL1ab.yaml
  • MC production 2: 20230205_GRB221009A_moon_strong_ph4_NSBDL1
  • Runs to analyse (run_id): 9614-9615
  • Tail cut cleaning (picture thd., boundary thd.): 22, 11
  • NSB tuning (from mean parameters values using run 9614 and 9615 subruns):
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 17.049,
   "extra_bias_in_dim_pixels": 4.316,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 34.3895,
  • Configuration files
fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph4_NSBDL1/lstchain_config_2023-01-20.json
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph4_NSBDL1/lstmcpipe_config_2023-01-20_PathConfigAllSkyFullDL1ab.yaml
  • MC production 3: 20230205_GRB221009A_moon_strong_ph5_NSBDL1
  • Runs to analyse (run_id): 9616-9617
  • Tail cut cleaning (picture thd., boundary thd.): 25, 12.5
  • NSB tuning (from mean parameters values using run 9616 and 9617 subruns):
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 21.134,
   "extra_bias_in_dim_pixels": 5.073,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 48.234,
  • Configuration files
fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph5_NSBDL1/lstchain_config_2023-01-20.json
/fefs/aswg/workspace/arnau.aguasca/scripts/lstMCpipe/v0.10.X/20230205_GRB221009A_moon_strong_ph5_NSBDL1/lstmcpipe_config_2023-01-20_PathConfigAllSkyFullDL1ab.yaml

DL1 data

  • original DL1a files
/fefs/aswg/workspace/franca.cassol/data/real/DL1/20221010_test1/dl1_LST-1.Run*
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/*
  • Date: 20221010
  • MC production 1
  • Produced DL1b data: tailcut30-15 (with cleaning based on pedestal RMS), dynamic cleaning
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL1/20221010/srcind/v0.9.9/tailcut3015_dynclg_picthd9602-9604/dl1_LST-1.Run0960*
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230205_GRB221009A_moon_strong_ph1_NSBDL1/
  • Date: 20221010
  • MC production 2
  • Produced DL1b data: tailcut34-17 (with cleaning based on pedestal RMS), dynamic cleaning
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL1/20221010/srcind/v0.9.9/tailcut3417_dynclg_picthd9605-9607/dl1_LST-1.Run0960*
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230205_GRB221009A_moon_strong_ph2_NSBDL1/


  • original DL1a files
/fefs/aswg/data/real/DL1/20221012/v0.9/tailcut84/dl1_LST-1.Run0961*
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_base_prod/*
  • Date: 20221012
  • MC production 1
  • Produced DL1b data: tailcut32-16 (with cleaning based on pedestal RMS), dynamic cleaning
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL1/20221012/srcind/v0.9.9/tailcut3216_dynclg/dl1_LST-1.Run0961*
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230205_GRB221009A_moon_strong_ph3_NSBDL1/
  • Date: 20221012
  • MC production 2
  • Produced DL1b data: tailcut22-11 (with cleaning based on pedestal RMS), dynamic cleaning
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL1/20221012/srcind/v0.9.9/tailcut2211_dynclg/dl1_LST-1.Run0961*
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230205_GRB221009A_moon_strong_ph4_NSBDL1/
  • Date: 20221012
  • MC production 3
  • Produced DL1b data: tailcut25-12.5 (with cleaning based on pedestal RMS), dynamic cleaning
/fefs/aswg/workspace/arnau.aguasca/Analysis/results/real/DL1/20221012/srcind/v0.9.9/tailcut2512.5_dynclg/dl1_LST-1.Run0961*
/fefs/aswg/workspace/MC_data_simlink/DL1/AllSky/20230205_GRB221009A_moon_strong_ph5_NSBDL1/

Random forest

  • lstchain-0.9.10
  • source-indep (diffuse gamma, src_r<3deg, disp_norm)
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230205_GRB221009A_moon_strong_ph1_NSBDL1
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230205_GRB221009A_moon_strong_ph2_NSBDL1
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230205_GRB221009A_moon_strong_ph3_NSBDL1
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230205_GRB221009A_moon_strong_ph4_NSBDL1
/fefs/aswg/workspace/MC_data_simlink/models/AllSky/20230205_GRB221009A_moon_strong_ph5_NSBDL1

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

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=GRB221009A_src-ind_analysis_moon&oldid=7088"