Editing Analysis BLLac Oct 2022 w IRF interpolation Priyadarshi

=BL Lacertae October 2022=


==General information==
* '''Name of the source''': BL Lacertae
* '''Brief description of the source''': 
** IBL / LBL
** highly variable
** z = 0.069
** RA: 22 02 43.3 (hh mm ss), Dec: +42 16 40 (dd mm ss)
** RA, Dec in deg (ICRS): 330.68038, 42.277772

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

==Data-taking information==
* '''Dates of data-taking''': 20221020

** 20221020 (2.46 h)
***ZD range (deg): [16.73, 43.89]
***AZ range (deg): [-57.25, -29.91]
***Runs: 9753 - 9761
***Joint observations with MAGIC? : Yes

==Monte Carlo information==
For getting a closer data-MC comparison, we produce MC data along various sky positions (nodes) as a preliminary method to facilitate a better Random Forest model for the reconstruction and to interpolate IRFs from various sky positions to the target position of each Observation Run.

== 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.9)

* original real DL1a data
 /fefs/aswg/data/real/DL1/{Date}/v0.9/tailcut84/

* lstchain v0.9 tailcut8-4 (with cleaning based on pedestal RMS, dynamical cleaning)
  /fefs/aswg/data/real/DL1/20221020/v0.9/tailcut84/log/lstchain_config_v0.9.4_from2022onwards.json

* MC nodes: tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, NSB tuning
  /fefs/aswg/data/mc/DL1/AllSky/20220601_dec3476_tuned/{Testing, Training}Dataset/

==Random forest==

Here we use the TrainingDataset nodes, selected along the declination path (or closer to) of the selected source (BL Lac here). The "dec3476" part in the production ID refers to the declination angle of +34.76 deg.
In this new RF model production, we include zenith and azimuth pointing parameters as well.

* lstchain v0.9.6dev source-independent, dynamic cleaning, NSB tuning
  /fefs/aswg/data/models/AllSky/20220601_dec3476_tuned/

=== For IRF interpolation ===

Here we create an almost uniform grid of nodes of sky positions for TestingDataset MC, so as to facilitate IRF interpolation in the parameter space of cos zenith and the sin of angle (delta) between the geomagnetic field and the particle shower direction. IRF interpolation is performed using pyirf v0.7.

* Link to MC DL2 files used:
/fefs/aswg/data/mc/DL2/AllSky/20220601_dec3476_tuned/TestingDataset/

* Particle types: ring-wobble gamma (offset range = [0.3999, 0.4001] deg)

[[File:Source_pointings_BLLac_Oct20_2022.png|1000 px]]

==DL2 data==

The DL2 files of the real data are stored here -
/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/BLLac/v09x/tailcut84/DL2/real/src_indep/20220601_dec3476/int50/Oct_2022/{Date}

==DL3 data selection==

=== For energy-dependent cuts ===

Information about your DL3 data selection. 

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

* source-independent
**gh_efficiency: 0.8
**theta_containment: 0.8

[[File:Off_regions_Run_9753_9761_w_safe_en_mask_min_25.1_GeV_1_off_10_energy_bin_p_dec.png ]]

* Energy-dependent Gammaness  and Theta cut distribution
[[File:Gammaness_cuts_all_BLLac_Oct20_2022.png|400 px]]
[[File:Rad_max_all_BLLac_Oct20_2022.png|400 px]]

==High-level analysis==

* lstchain Tools to generate source-independent IRF and DL3
* Science Tool: gammapy 0.20.1
* point-like IRF, 1D analysis

For DL3 data reduction to DL4 dataset,
* OFF regions selected - 1 Wobble
* Safe Energy mask for systematics - Yes
* Energy axes -
** True: 10 GeV to 10 TeV with 10 bins per decade
** Reco: 10 GeV to 10 TeV with 10 bins per decade

For spectral analysis, 
* Model used - Log Parabola + EBL (Dominguez)
* Reference energy used by calculating the decorrelation energy with a Power Law fit, 116 GeV
* Spectral Fit energy range - 25.1 GeV to 10 TeV with 10 bins per decade
* Spectral fit with gammapy default optimization - minuit

For LC,
* E_min = 100 Gev

==Analysis Results==

=== Theta2 plot ===
* Single energy bin
[[File:Theta2_counts_BLLac_9_wob_obs_Run_9753_to_9761_1_en_bins_cut_04_std.png|400 px]]
[[File:Theta2_excess_BLLac_9_wob_obs_Run_9753_to_9761_1_en_bins_cut_04_std.png|400 px]]

* 10 bins per decade in energy
[[File:Theta2_counts_BLLac_9_wob_obs_Run_9753_to_9761_30_en_bins_cut_04_std.png|1000 px]]
[[File:Theta2_excess_BLLac_9_wob_obs_Run_9753_to_9761_30_en_bins_cut_04_std.png|1000 px]]

=== Significance map ===

=== Excess map ===

=== Spectral results ===

==== For energy-dependent cuts ====
* 20 Oct dataset

Preliminary LC and SED

[[File:Model_parameters_correlation_Run_9753_9761_lp_ebl_dominguez_1_off_e_ref_116GeV_fit_range_w_10_bin_p_decade_10GeV_10TeV_w_safe_en_mask.png|500 px]]

[[File:Model_parameters_Run_9753_9761_lp_ebl_dominguez_1_off_e_ref_116GeV_fit_range_w_10_bin_p_decade_10GeV_10TeV_w_safe_en_mask.png|1000 px]]

[[File:LC_min_100_Run_9753_9761_lp_ebl_dominguez_1_off_e_ref_116GeV_fit_range_w_10_bin_p_decade_10GeV_10TeV_w_safe_en_mask.png|1000 px]]

[[File:SED_Run_9753_9761_lp_ebl_dominguez_1_off_e_ref_116GeV_fit_range_w_10_bin_p_decade_10GeV_10TeV_w_safe_en_mask.png]]

* SED comparison with Flare of Aug 8, 2021
[[File:SED_comparison_flares_8Aug_2021_v_20Oct_2022.png]]