Crab Nebula Performance

Crab Nebula Performance

Draft versions

Version 1, 20220923: [[1]]

Version 2, 20230124: [[2]]

Version 3, 20230126: [[3]]

Version 4, 20230321: [[4]]

Version 5, 20230329: [[5]]

Format was changes after A&A considered the paper topic was not suitable for their journal

Version 6, sent to ApJ, 20230421: [[6]]

Version 9, addressing the ApJ referee's comments, 20230623: [[7]]

Version 10, minor changes (acknowledgements, author list, fig 16 typo, zenodo doi for repo containing the paper's plots data), 20230719: [[8]]

General information

• Name of the source: Crab Nebula and pulsar
• Brief description of the source:
  • Object type: PWN
  • Other relevant information
  • RA: 05 34 31.9 (hh mm ss), Dec: +20 00 52.2 (dd mm ss)
  • RA, Dec in deg (ICRS): 83.6331, 22.0145

• Analysis by:
  • Abelardo Moralejo (IFAE, moralejo@ifae.es)
  • Rubén López-Coto (IAA-CSIC, rlopezcoto@iaa.es)
  • Thomas Vuillaume (LAPP, thomas.vuillaume@lapp.in2p3.fr)
  • Seiya Nozaki (MPP, nozaki@mpp.mpg.de)
  • Daniel Morcuende (UCM, dmorcuen@ucm.es)

Data-taking information

• Data selection based on the notebook: https://github.com/cta-observatory/cta-lstchain/blob/master/notebooks/data_selection.ipynb
  • Excluded runs with issues in interleaved pedestals
  • Run selection also based on "healthy intensity spectra", i.e., rate of cosmics with intensity between 80 and 120 p.e. larger than 800 events/s.
  • Runs taken in wobble mode. Selected wobble pointings at ~0.4 deg from the source [0.35, 0.45] deg

• Observation time: 35.94 hours

• Dates of data-taking:

1 : 2020-11-17 : [2914]
2 : 2020-11-18 : [2929, 2930, 2931, 2932, 2933, 2934]
3 : 2020-11-19 : [2949, 2950]
4 : 2020-11-20 : [2967, 2968, 2969, 2970, 2971, 2972, 2973, 2974, 2975, 2976, 2977]
5 : 2020-11-21 : [2988, 2989, 2990, 2991, 2992]
6 : 2020-11-22 : [3004, 3005, 3006, 3007, 3008]
7 : 2020-12-07 : [3093, 3094, 3095, 3096]
8 : 2020-12-13 : [3231, 3232]
9 : 2020-12-14 : [3243]
10 : 2020-12-15 : [3270, 3271, 3272, 3273, 3274, 3275, 3276, 3277, 3278, 3279]
11 : 2020-12-17 : [3318, 3319, 3320, 3321]
12 : 2020-12-18 : [3328, 3329, 3330]
13 : 2020-12-19 : [3338, 3339, 3340]
14 : 2020-12-20 : [3355, 3356]
15 : 2020-12-21 : [3373]
16 : 2021-02-11 : [3598, 3599, 3600, 3601]
17 : 2021-02-12 : [3615]
18 : 2021-02-13 : [3632, 3633, 3634, 3635]
19 : 2021-02-15 : [3672, 3673, 3674, 3675, 3676, 3677]
20 : 2021-02-16 : [3706, 3707, 3708]
21 : 2021-03-15 : [4067, 4068]
22 : 2021-03-16 : [4086, 4087]
23 : 2021-09-04 : [6045]
24 : 2021-09-05 : [6073]
25 : 2021-09-14 : [6304]
26 : 2022-02-04 : [6872, 6873, 6874, 6875]
27 : 2022-02-05 : [6892, 6893, 6894, 6895]
28 : 2022-02-23 : [7097, 7098, 7099]
29 : 2022-02-28 : [7133, 7136]
30 : 2022-03-01 : [7161]
31 : 2022-03-02 : [7195, 7196, 7197, 7199, 7200]
32 : 2022-03-03 : [7227, 7228, 7231, 7232, 7233]
33 : 2022-03-04 : [7253, 7254, 7255, 7256]
34 : 2022-03-05 : [7274, 7275, 7276, 7277]

• List of all selected runs:

[2914, 2929, 2930, 2931, 2932, 2933, 2934, 2949, 2950, 2967, 2968, 2969,
 2970, 2971, 2972, 2973, 2974, 2975, 2976, 2977, 2988, 2989, 2990, 2991,
 2992, 3004, 3005, 3006, 3007, 3008, 3093, 3094, 3095, 3096, 3231, 3232,
 3243, 3270, 3271, 3272, 3273, 3274, 3275, 3276, 3277, 3278, 3279, 3318,
 3319, 3320, 3321, 3328, 3329, 3330, 3338, 3339, 3340, 3355, 3356, 3373,
 3598, 3599, 3600, 3601, 3615, 3632, 3633, 3634, 3635, 3672, 3673, 3674,
 3675, 3676, 3677, 3706, 3707, 3708, 4067, 4068, 4086, 4087, 6045, 6073,
 6304, 6872, 6873, 6874, 6875, 6892, 6893, 6894, 6895, 7097, 7098, 7099,
 7133, 7136, 7161, 7195, 7196, 7197, 7199, 7200, 7227, 7228, 7231, 7232,
 7233, 7253, 7254, 7255, 7256, 7274, 7275, 7276, 7277]

• ZD range (deg): 0 - 35 deg
• AZ range (deg): 90 - 270 deg

• Transmission at 9 km: See corresponding values in the ELOG entries. We do not use it for data selection. Instead we apply selection cuts based on rates (see above).
• Joint observations with MAGIC? : No

Monte Carlo information

• Link to MC files used:
  • Particle types:
  • ZD (deg): Zenith Range
  • AZ (deg): Azimuth Range

DL1 data

Used standard DL1b files produced by lstosa (using lstchain v0.9.x):

• real data: tailcut8-4 (plus cleaning based on pixel-wise pedestal std dev & dynamic cleaning, see paper):

/fefs/aswg/data/real/DL1/YYYYMMDD/v0.9/tailcut84/dl1_LST-1.Run?????.h5

• MC: tailcut8-4 (plus cleaning based on pixel-wise pedestal std dev, dynamic cleaning). NSB tuning to Crab:

/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/TrainingDataset/dec_2276/
/fefs/aswg/data/mc/DL1/AllSky/20221027_v0.9.9_crab_tuned/TestingDataset/ (see paper, closest pointing node used for each data run)

• 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.

• lstchain-0.9.9
• source-indep:
  • config: /fefs/aswg/data/models/AllSky/20221027_v0.9.9_crab_tuned/dec_2276/lstchain_config_2022-10-27.json

• source-dep:
  • config: /fefs/aswg/workspace/seiya.nozaki/data/MC/v0.7.5/tailcut84_dynamic_bllac/srcdep/RF/lstchain_src_dep_config.json
  • models: /fefs/aswg/data/models/AllSky/20221027_v0.9.9_crab_tuned/dec_2276/

DL2 data

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

• lstchain-0.9.9
• source-indep
• source-dep

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.