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: [[]]

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

General information

• Name of the source: Crab Nebula
• 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 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

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.

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

• original DL1 files

/fefs/aswg/data/real/DL1/YYYYMMDD/v0.9/tailcut84/dl1_LST-1.Run?????.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.