1ES1959+650 Preliminary Analysis with lstchain v0.7.x Priyadarshi

1ES1959+650 Preliminary Analysis with lstchain v0.7.x Priyadarshi

General information

• Name of the source: 1ES 1959+650
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.048
  • nearly HBL
  • RA: 19 59 59.85 (hh mm ss), Dec: +65 8 54.65 (dd mm ss)
  • RA, Dec in deg (ICRS): 299.99938562, 65.14851421

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

Data-taking information

• Dates of data-taking: 20200711, 20200713, 20200715, 20200726, 20200809, 20200811, 20200812, 20200820, 20210319, 20210530, 20210807

• • 20200711 (1.01 h)
    • ZD range (deg): [37.63, 41.1]
    • AZ range (deg): [12.27, 18.9]
    • Runs: 2193 - 2196
    • Joint observations with MAGIC? : No

• • 20200713 (1.08 h)
    • ZD range (deg): [37.48, 43.12]
    • AZ range (deg): [11.77, 21.73]
    • Runs: 2202, 2203, 2205, 2206, 2207
    • Joint observations with MAGIC? : No

• • 20200726 (0.79 h)
    • ZD range (deg): [46.6, 51.9]
    • AZ range (deg): [-28.33, -25.02]
    • Runs: 2348, 2350, 2351, 2352
    • Joint observations with MAGIC? : No

• • 20200809 (1.28 h)
    • ZD range (deg): [36.27, 37.04]
    • AZ range (deg): [-7.86, 6.38]
    • Runs: 2370 - 2375
    • Joint observations with MAGIC? : No

• • 2020811 (1.52 h)
    • ZD range (deg): [36.39, 40.11]
    • AZ range (deg): [-14.98, 2.75]
    • Runs: 2408, 2411 - 2415
    • Joint observations with MAGIC? : No

• • 20200812 (0.73 h)
    • ZD range (deg): [43.1, 45.65]
    • AZ range (deg): [-24.57, -21.72]
    • Runs: 2420 - 2422
    • Joint observations with MAGIC? : No

• • 20210319 (1.14 h)
    • ZD range (deg): [48.05, 60.3]
    • AZ range (deg): [25.98, 28.77]
    • Runs: 4171 - 4177
    • Joint observations with MAGIC? : Not for 4175 - 4177

• • 20210530 (1.04 h)
    • ZD range (deg): [40.3, 44.54]
    • AZ range (deg): [17.47, 23.28]
    • Runs: 4840, 4841, 4843, 4844
    • Joint observations with MAGIC? : Yes

• • 20210807 (4.08 h)
    • ZD range (deg): [36.08, 49.04]
    • AZ range (deg): [-27.06, 21.7]
    • Runs: 5520 - 5533, 5535, 5540, 5542
    • Joint observations with MAGIC? : Yes

Monte Carlo information

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

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

Information about your DL3 data selection.

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

• source-indep, energy-dependent cuts
  • gh_efficiency: 0.8
  • theta_containment: 0.8

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
• Reference energy used by calculating the decorrelation energy with a Power Law fit, ~370 GeV
• Spectral Fit energy range - 10 GeV to 40 TeV with 5 bins per decade
• Spectral fit with gammapy default optimization - minuit

Analysis Results

Spectral results

Preliminary LC and SED, for the full dataset-