All sky mc bllac

From my_wiki
Revision as of 14:09, 28 October 2022 by Nozaki (talk | contribs) (How much do we need to tune?)
Jump to: navigation, search

All Sky MC simulation data for BL Lac analysis

Do we need to apply NSB tuning?

Probably yes. It is because NSB level is 1.3-1.4 x dark NSB in BL Lac FoV due to the low galactic latitude (-10.44 degree)

Bllac nsb level.png

How much do we need to tune?

(lst_v0.9.6) [seiya.nozaki@cp02 DL1]$ lstchain_tune_nsb --input-mc 
/fefs/aswg/workspace/georgios.voutsinas/AllSky/TrainingDataset/GammaDiffuse/dec_4822/sim_telarray/node_corsika_theta_19.456_az_0.0_/output_v1.4/simtel_corsika_theta_19.456_az_0.0_run10.simtel.gz  --input-data 
/fefs/aswg/data/real/DL1/20210808/v0.9/tailcut84/dl1_LST-1.Run05555.0000.h5 --config lstchain_standard_config.json
Real data: median across camera of good pixels' pedestal std 2.101 p.e.
Number of pixels beyond 3 std dev of median: 31, (above 2.80 p.e.)
Good and not too bright pixels: 1818
{
  "increase_nsb": true,
  "extra_noise_in_dim_pixels": 0.977,
  "extra_bias_in_dim_pixels": 0.388,
  "transition_charge": 8,
  "extra_noise_in_bright_pixels": 1.197
}
(lst_v0.9.9) [seiya.nozaki@cp02 v0.9.9]$ lstchain_tune_nsb --input-mc /fefs/aswg/workspace/georgios.voutsinas/AllSky/TrainingDataset/GammaDiffuse/dec_4822/sim_telarray/node_corsika_theta_19.456_az_0.0_/output_v1.4/simtel_corsika_theta_19.456_az_0.0_run10.simtel.gz   --input-data /fefs/aswg/data/real/DL1/20210808/v0.9/tailcut84/dl1_LST-1.Run05555.0000.h5 --config 
lstchain_standard_config.json 
Real data: median across camera of good pixels' pedestal std 2.101 p.e.
Number of pixels beyond 3 std dev of median: 31, (above 2.80 p.e.)
Good and not too bright pixels: 1818


{
  "increase_nsb": true,
  "extra_noise_in_dim_pixels": 0.955, 
  "extra_bias_in_dim_pixels": 0.388,
  "transition_charge": 8,
  "extra_noise_in_bright_pixels": 1.197
}

Which MC is proper for BL Lac analysis?

dec_3476 and dec_4822 is the closest nodes to the BL Lac path.

training MC node

How much does NSB tuning affect the analysis?

You can see some differences below a few hundreds of GeV on intensity distribution. A ratio of surviving events after image cleaning is also different (~5% below 100 GeV).

quick check