Analysis April 2022 with lstchain v0.9.X Pirola (cross-check)

Relevant Analysis Name

General information

• Name of the source: LHAASO J2108+5157
• Brief description of the source:

- Object type : Unidentified Galactic PeVatron candidate

- Distance (pc) : Unknown

- RA, Dec in deg (ICRS): 317.22, 51.95

• Analysis by Giorgio Pirola (MPP, gpirola@mpp.mgp.de)

Data-taking information (Run selection)

• Run selection with the use of the notebook from Abelardo Moralejo [1]
• Atmospheric transmission extracted from ELOG
• The analysis has been performed on 2 samples, with different limits on the Zenith range:

- Zenith < 45 deg

- Zenith < 55 deg (in agreement with the primary analysis by Jakub Jurysek)

• Summary of selection cuts:

- Wobble in (0.45, 0.55)deg

- elapsed_time > 5 min

- transmission_cut > 0.65

- pedestal charge std dev < 1.8 p.e.

- Cosmic rate > 3000 ev/s

- Cosmic rate( > 10 p.e.) > 20 ev/s

- Cosmic rate( > 30 p.e.) > 3 ev/s

- muon ring width std dev < 0.023

• Run selection after cuts:

- Zenith < 45 deg:

 Total wobble runs: 162
 Observation time: 47.16 hours
 Selected Runs:
 1 : 2021-06-04 : [4913, 4914, 4915, 4916, 4917]
 2 : 2021-06-05 : [4935, 4936]
 3 : 2021-06-12 : [5028, 5029, 5030, 5031]
 4 : 2021-06-30 : [5072]
 5 : 2021-07-01 : [5082, 5083, 5084]
 6 : 2021-07-02 : [5091, 5092, 5093]
 7 : 2021-07-03 : [5101, 5102, 5103, 5104, 5105, 5106, 5107, 5108]
 8 : 2021-07-04 : [5118, 5119, 5120, 5121]
 9 : 2021-07-05 : [5138, 5139, 5140, 5141, 5142]
 10 : 2021-07-15 : [5270, 5272]
 11 : 2021-08-01 : [5440, 5441, 5442]
 12 : 2021-08-03 : [5461, 5462, 5463, 5464, 5465]
 13 : 2021-08-04 : [5473, 5474, 5475, 5476, 5477, 5478, 5479, 5480]
 14 : 2021-08-05 : [5491, 5492, 5493, 5494, 5497, 5498, 5499, 5500]
 15 : 2021-08-06 : [5505, 5506, 5507, 5508, 5509, 5510, 5511, 5512, 5513, 5514, 5515, 5516, 5517]
 16 : 2021-08-08 : [5576]
 17 : 2021-08-09 : [5590, 5591]
 18 : 2021-08-10 : [5641, 5642, 5643]
 19 : 2021-08-11 : [5681, 5682, 5683, 5684, 5685, 5686, 5687]
 20 : 2021-08-12 : [5707, 5708, 5709, 5710, 5711, 5712, 5713]
 21 : 2021-08-13 : [5727]
 22 : 2021-09-01 : [5947, 5948, 5949, 5950, 5952]
 23 : 2021-09-02 : [5980, 5981, 5982, 5983, 5984, 5985, 5986, 5987, 5988]
 24 : 2021-09-03 : [5999, 6000, 6001, 6002, 6003, 6004, 6005, 6006, 6007, 6008, 6009]
 25 : 2021-09-04 : [6023, 6024, 6034, 6035, 6036, 6037, 6038]
 26 : 2021-09-05 : [6057, 6058, 6059, 6060, 6061, 6062, 6063, 6064, 6065, 6066]
 27 : 2021-09-06 : [6079, 6080, 6082, 6083, 6084, 6085]
 28 : 2021-09-07 : [6130, 6131, 6132, 6133, 6134]
 29 : 2021-09-09 : [6175, 6176, 6177, 6178, 6179, 6180, 6181, 6182, 6183]
 30 : 2021-09-11 : [6230, 6231]
 31 : 2021-09-12 : [6254, 6255, 6256]

- Zenith < 55 deg:

 Total wobble runs: 177
 Observation time: 51.40 hours
 Selected Runs:
 1 : 2021-06-04 : [4913, 4914, 4915, 4916, 4917]
 2 : 2021-06-05 : [4935, 4936]
 3 : 2021-06-12 : [5028, 5029, 5030, 5031]
 4 : 2021-06-30 : [5071, 5072]
 5 : 2021-07-01 : [5080, 5081, 5082, 5083, 5084]
 6 : 2021-07-02 : [5091, 5092, 5093]
 7 : 2021-07-03 : [5101, 5102, 5103, 5104, 5105, 5106, 5107, 5108]
 8 : 2021-07-04 : [5115, 5116, 5117, 5118, 5119, 5120, 5121]
 9 : 2021-07-05 : [5135, 5136, 5137, 5138, 5139, 5140, 5141, 5142]
 10 : 2021-07-15 : [5270, 5272]
 11 : 2021-08-01 : [5440, 5441, 5442]
 12 : 2021-08-03 : [5461, 5462, 5463, 5464, 5465]
 13 : 2021-08-04 : [5473, 5474, 5475, 5476, 5477, 5478, 5479, 5480]
 14 : 2021-08-05 : [5491, 5492, 5493, 5494, 5497, 5498, 5499, 5500]
 15 : 2021-08-06 : [5505, 5506, 5507, 5508, 5509, 5510, 5511, 5512, 5513, 5514, 5515, 5516, 5517]
 16 : 2021-08-08 : [5576]
 17 : 2021-08-09 : [5590, 5591]
 18 : 2021-08-10 : [5641, 5642, 5643]
 19 : 2021-08-11 : [5681, 5682, 5683, 5684, 5685, 5686, 5687]
 20 : 2021-08-12 : [5707, 5708, 5709, 5710, 5711, 5712, 5713]
 21 : 2021-08-13 : [5727]
 22 : 2021-09-01 : [5947, 5948, 5949, 5950, 5952]
 23 : 2021-09-02 : [5980, 5981, 5982, 5983, 5984, 5985, 5986, 5987, 5988, 5989, 5990, 5991]
 24 : 2021-09-03 : [5999, 6000, 6001, 6002, 6003, 6004, 6005, 6006, 6007, 6008, 6009, 6010]
 25 : 2021-09-04 : [6023, 6024, 6034, 6035, 6036, 6037, 6038]
 26 : 2021-09-05 : [6057, 6058, 6059, 6060, 6061, 6062, 6063, 6064, 6065, 6066]
 27 : 2021-09-06 : [6079, 6080, 6082, 6083, 6084, 6085]
 28 : 2021-09-07 : [6130, 6131, 6132, 6133, 6134]
 29 : 2021-09-09 : [6175, 6176, 6177, 6178, 6179, 6180, 6181, 6182, 6183]
 30 : 2021-09-11 : [6230, 6231, 6233]
 31 : 2021-09-12 : [6254, 6255, 6256, 6257]

Monte Carlo information

- 20200629_prod5_trans_80

- 20, 40 deg zenith angle, south-pointing

- Note: 20 deg zenith angle MC was tuned manually using lstchain scripts, while for 40 deg MC we used lstmcpipe production [2] and thus the directory structure for the two is different

• PSF and NSB tuned:

"image_modifier": {
      "increase_nsb": true,
      "extra_noise_in_dim_pixels": 0.705,
      "extra_bias_in_dim_pixels": 0.269,
      "transition_charge": 8,
      "extra_noise_in_bright_pixels": 0.673,
      "increase_psf": true,
      "smeared_light_fraction": 0.125
      }

• Tuned DL1b MC:

- 20 deg zenith angle:

/fefs/aswg/workspace/jakub.jurysek/simulations/dl1/20200629_prod5_trans_80/gamma/zenith_20deg/south_pointing/20220225_v0.9.1_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/off0.4deg/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl1/20200629_prod5_trans_80/gamma-diffuse/zenith_20deg/south_pointing/20220225_v0.9.1_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl1/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20220225_v0.9.1_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl1/20200629_prod5_trans_80/electron/zenith_20deg/south_pointing/20220225_v0.9.1_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/

- 40 deg zenith angle:

/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/off0.4deg/
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/gamma-diffuse/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/electron/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/

• Random Forests - source independent:

- 20 deg zenith angle:

/fefs/aswg/workspace/jakub.jurysek/simulations/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20220225_v0.9.2_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/

- 40 deg zenith angle:

/fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/

• DL2 MC (testing only):

- 20 deg zenith angle:

/fefs/aswg/workspace/jakub.jurysek/simulations/dl2/20200629_prod5_trans_80/gamma/zenith_20deg/south_pointing/20220225_v0.9.2_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/off0.4deg/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl2/20200629_prod5_trans_80/gamma-diffuse/zenith_20deg/south_pointing/20220225_v0.9.2_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl2/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20220225_v0.9.2_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/
/fefs/aswg/workspace/jakub.jurysek/simulations/dl2/20200629_prod5_trans_80/electron/zenith_20deg/south_pointing/20220225_v0.9.2_prod5_trans_80_dl1ab_tune_MC_to_LHAASO/

- 40 deg zenith angle:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/proton/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/
/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/gamma/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/off0.4deg/
/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/gamma-diffuse/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/
/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/electron/zenith_40deg/south_pointing/20220308_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_nsb_psf_src2/

DL1 data

Please include your settings to produce your DL1a and DL1b files. This includes whether you use LSTOSA, specific versions of lstchain, dllab scripts, cleaning levels, and calibration information.

Example :

Reprocessed using DL1a files produced by LSTOSA (lstchain v0.7.3) and dl1ab script (v0.7.5)

• original DL1a files

/fefs/aswg/data/real/DL1/20210808/v0.7.3/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/XXX/zenith_20deg/south_pointing/20210416_v0.7.3_prod5_trans_80_local_taicut_8_4/

• lstchain v0.7.5
• real data: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning

 /fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/lstchain_standard_config.json

• MC: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF+NSB tuning

 "image_modifier": {
   "increase_nsb": true,
   "extra_noise_in_dim_pixels": 0.74,
   "extra_bias_in_dim_pixels": 0.38,
   "transition_charge": 8,
   "extra_noise_in_bright_pixels": 0.052,
   "increase_psf": true,
   "smeared_light_fraction": 0.125
 },

 /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/lstchain_dl1ab_tune_MC_to_BLLac_config.json

• Produced DL1b files

/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/DL1_raw/data/
/fefs/aswg/workspace/seiya.nozaki/data/BLLac/v0.7.3/tailcut84_dynamic_v075/DL1_raw/Run0XXXX

Random forest

Please include your specific .json files used in producing your RFs and for which source analysis (source dependent or independent), and variabules used in the RF.

Example

• lstchain-0.7.6.dev242+g2086cb9

• source-indep (diffuse gamma, src_r<3deg, disp_norm)

/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcindep/RF/lstchain_standard_config.json

• source-dep

/fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcdep/RF/lstchain_src_dep_config.json

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.