1ES0647+250 Preliminary Analysis with lstchain v0.7.x Priyadarshi

1ES0647+250 Preliminary Analysis with lstchain v0.7.x Priyadarshi

General information

• Name of the source: 1ES 0647+250
• Brief description of the source:
  • Object type : AGN
  • Redshift (z) : 0.41 [1]
  • HBL
  • RA: 6 50 46.49 (hh mm ss), Dec: +25 2 59.62 (dd mm ss)
  • RA, Dec in deg (ICRS): 102.69370933, 25.04989555

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

Data-taking information

• Dates of data-taking: 20201215, 20201218, 20201219, 20201220

• • 20201215 (2.42 h)
    • ZD range (deg): [5.57, 39.61]
    • AZ range (deg): [233.4, 275.37]
    • Runs: 3281 - 3288
    • Joint observations with MAGIC? : Not for 3285 - 3288

• • 20201218 (1.57 h)
    • ZD range (deg): [9.75, 28.38]
    • AZ range (deg): [249.43, 269.94]
    • Runs: 3331 - 3335
    • Joint observations with MAGIC? : Yes

• • 20201219 (2.21 h)
    • ZD range (deg): [12.53, 46.32]
    • AZ range (deg): [253.63, 278.23]
    • Runs: 3341 - 3344, 3346, 3347, 3348
    • Joint observations with MAGIC? : Yes for 3341 - 3343

• • 20201220 (1.95 h)
    • ZD range (deg): [4.48, 24.04]
    • AZ range (deg): [151.1, 268.37]
    • Runs: 3357 - 3364
    • Joint observations with MAGIC? : Not for 3357 - 3359

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)

• original real DL1a data

/fefs/aswg/data/real/DL1/{Date}/v0.7.1/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

• 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.16 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, ~56 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-