PG1553+113 Preliminary Analysis with lstchain v0.7.x Priyadarshi

General information

Name of the source: PG 1553+113
Brief description of the source:
- Object type : AGN
- Redshift (z) : 0.433 Dorigo Jones et al
- BL Lac
- RA: 15 5543.044 (hh mm ss), Dec: +11 11 24.365 (dd mm ss)
- RA, Dec in deg (ICRS): 238.92934976, 11.19010155

Analysis by Chaitanya Priyadarshi (IFAE - cpriyadarshi@ifae.es)

Data-taking information

Dates of data-taking: 20210321, 20210407, 20210408, 20210410, 20210420, 20210508, 20210604, 20210605, 20210606, 20210609

- 20210321 (0.31 h)
  - ZD range (deg): [39.76]
  - AZ range (deg): [108.39]
  - Runs: 4207
  - Joint observations with MAGIC? : Yes

- 20210407 (1.11 h)
  - ZD range (deg): [18.42, 26.43]
  - AZ range (deg): [197.5, 231.14]
  - Runs: 4333 - 4336
  - Joint observations with MAGIC? : No

- 20210408 (0.43 h)
  - ZD range (deg): [17.71, 18.46]
  - AZ range (deg): [183.62, 198.43]
  - Runs: 4343, 4344
  - Joint observations with MAGIC? : No

- 20210410 (1.96 h)
  - ZD range (deg): [17.48, 48.84]
  - AZ range (deg): [100.79, 190.59]
  - Runs: 4363 - 4370
  - Joint observations with MAGIC? : No

- 20210420 (0.29 h)
  - ZD range (deg): [20.15, 22.06]
  - AZ range (deg): [211.01, 219.81]
  - Runs: 4514, 4515
  - Joint observations with MAGIC? : No

- 20210508 (0.83 h)
  - ZD range (deg): [34, 42.45]
  - AZ range (deg): [105.21, 114.84]
  - Runs: 4652 - 4654
  - Joint observations with MAGIC? : Yes

- 20210604 (0.88 h)
  - ZD range (deg): [18.8, 24.67]
  - AZ range (deg): [131.9, 162.59]
  - Runs: 4902 - 4905
  - Joint observations with MAGIC? : No

- 20210605 (0.79 h)
  - ZD range (deg): [33.95, 41.33]
  - AZ range (deg): [245.96, 253.11]
  - Runs: 4930 - 4932
  - Joint observations with MAGIC? : Yes

- 20210606 (1.28 h)
  - ZD range (deg): [18.21, 27.27]
  - AZ range (deg): [125.81, 172.56]
  - Runs: 4940 - 4943
  - Joint observations with MAGIC? : No

- 20210609 (0.89 h)
  - ZD range (deg): [19.11, 28.35]
  - AZ range (deg): [123.68, 139.46]
  - Runs: 4987 - 4989
  - Joint observations with MAGIC? : No

Monte Carlo information

For energy-dependent cuts

Link to MC DL2 files used:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/ /fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/

- Particle types: point-gamma (offset=0.4 deg), proton, electron
- ZD (deg): [20, 40]
- AZ (deg): [180, 180]

For global cuts

Link to MC DL2 files used:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_20deg/south_pointing/20210416_v0.7.3_prod5_trans_80_local_taicut_8_4/

- Particle types: point-gamma (offset=0.4 deg), proton, electron
- ZD (deg): [20]
- AZ (deg): [180]

DL1 data

This includes whether you use LSTOSA, specific versions of lstchain, dllab scripts, cleaning levels, and calibration information.

DL1a files produced by LSTOSA (lstchain v0.7.1 and v0.7.3)

original real DL1a data

/fefs/aswg/data/real/DL1/{Date}/v0.7.{1,3}/tailcut84/

lstchain v0.7.1 tailcut8-4 (with cleaning based on pedestal RMS)

 /fefs/aswg/data/real/DL1/20210319/v0.7.1/log/lstchain_standard_config_tailcut84.json

lstchain v0.7.3 tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning

 /fefs/aswg/data/real/DL1/20210320/v0.7.3/tailcut84/log/lstchain_standard_config_tailcut84.json

MC (zen 20 deg): tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning

 /fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/lstchain_standard_config_v075_Tel1_tail84_dynamic_clean.json

MC (zen 40 deg): tailcut8-4 (with cleaning based on pedestal RMS)

/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/lstchain_standard_config_v073_Tel1_tail_8_4_MC_config.json

Random forest

lstchain v0.7.5 source-indep, zenith 20 deg, azimuth 180 deg, dynamic cleaning

 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20210923_v0.7.5_prod5_trans_80_dynamic_cleaning/lstchain_standard_config_v075_Tel1_tail84_dynamic_clean.json

lstchain v0.7.3 source-indep, zenith 40 deg, azimuth 180 deg

 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20210506_v0.7.3_prod5_trans_80_zen40deg_local_tailcut_8_4/lstchain_standard_config_v073_Tel1_tail_8_4_MC_config.json

DL2 data

Similar to DL1 files

DL3 data selection

For energy-dependent cuts

Information about your DL3 data selection.

intensity > 0
r: [0, 1]
wl: [0.01, 1]
leakage_intensity_width_2: [0, 1]

source-indep
- gh_efficiency: 0.8
- theta_containment: 0.8

For global cuts

Information about your DL3 data selection.

intensity > 0
r: [0, 1]
wl: [0.01, 1]
leakage_intensity_width_2: [0, 1]

source-indep
- global_gh_cut: 0.7
- global_theta_cut: 0.2

High-level analysis

lstchain Tools to generate source-independent IRF and DL3
Science Tool: gammapy 0.20.dev309+gb11005a (Developer build to extract signal from energy-dependent regions)
point-like IRF, 1D analysis

A total of 8.69 hrs of data was analyzed

For DL3 data reduction to DL4 dataset,

OFF regions selected - 1 Wobble
Safe Energy mask for systematics - No
Energy axes -
- True: 10 GeV to 100 TeV with 10 bins per decade
- Reco: 10 GeV to 40 TeV with 5 bins per decade

For spectral analysis,

Model used - Log Parabola, Power Law
Reference energy used by calculating the decorrelation energy with a Power Law fit, ~86 GeV
Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
Spectral fit with gammapy default optimization - minuit