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

From my_wiki
Revision as of 20:24, 3 March 2022 by Chaitanya (talk | contribs) (Data-taking information)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

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

General information[edit]

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

  • 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
    • 20200811 (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[edit]

  • Link to MC DL2 files used:

/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): [40]
    • AZ (deg): [180]

DL1 data[edit]

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

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

Similar to DL1 files


DL3 data selection[edit]

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

  • 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

Gammaness cuts distribution 80% efficiency different MC pointing.png Theta cuts distribution 80% efficiency different MC pointing.png

Effective area 80% efficiency different MC pointing.png Energy dispersion 80% efficiency different MC pointing.png

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

Theta2 plot[edit]

Significance map[edit]

Excess map[edit]

Spectral results[edit]

Preliminary LC and SED, for the full dataset-

LC Run 2203 5542 lp model e ref 371GeV fit 30GeV 100TeV.png SED Run 2203 5542 lp model e ref 371GeV fit range w 5 bin p decade 10 40.png