Editing 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 [https://pos.sissa.it/395/792]
** 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

[[File:Gammaness_cuts_distribution_80%_efficiency_different_MC_pointing.png]] [[File:Theta_cuts_distribution_80%_efficiency_different_MC_pointing.png]]

[[File:Effective_area_80%_efficiency_different_MC_pointing.png]] [[File:Energy_dispersion__80%_efficiency_different_MC_pointing.png]]

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==

=== Theta2 plot ===

=== Significance map ===

=== Excess map ===

=== Spectral results ===

Preliminary LC and SED, for the full dataset-

[[File:LC_Run_3281_3364_lp_model_e_ref_56GeV_fit_30GeV_100TeV.png]]
[[File:SED_Run_3281_3364_lp_model_e_ref_56GeV_fit_range_w_5_bin_p_decade_10_40.png]]