Analysis BLLac Aug 2021 Priyadarshi

From my_wiki
Jump to: navigation, search

BL Lacertae flare August 8, 2021[edit]

General information[edit]

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

  • Dates of data-taking:
    • 2021-08-08 (1.93 h)
      • ZD: 13 - 32 deg
      • Transmission9km: >80%
      • Runs: 5552-5559
      • Joint observations with MAGIC

Monte Carlo information[edit]

For energy-dependent cuts[edit]

  • Link to MC DL2 files used:

/fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf /fefs/aswg/data/mc/DL2/20200629_prod5_trans_80/{particle}/zenith_40deg/south_pointing/20220304_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_psf

    • Particle types: point-gamma (offset=0.4 deg), proton, electron
    • ZD (deg): [20, 40]
    • AZ (deg): [180, 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.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/20210808/v0.9/tailcut84/log/lstchain_config_tailcut84_v092.json
  • MC (zen 20 deg): tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF tuning
 /fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf/lstchain_config.json
  • MC (zen 40 deg): tailcut8-4 (with cleaning based on pedestal RMS), dynamic cleaning, PSF tuning
/fefs/aswg/data/mc/DL1/20200629_prod5_trans_80/proton/zenith_40deg/south_pointing/20220304_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_psf/lstchain_config.json

Random forest[edit]

  • lstchain v0.9 source-indep, zenith 20 deg, azimuth 180 deg, dynamic cleaning, PSF tuning
 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_20deg/south_pointing/20220303_v0.9.3_prod5_trans_80_zen20az180_dl1ab_tuned_psf/lstchain_config.json
  • lstchain v0.9 source-indep, zenith 40 deg, azimuth 180 deg, dynamic cleaning, PSF tuning
 /fefs/aswg/data/models/20200629_prod5_trans_80/zenith_40deg/south_pointing/20220304_v0.9.3_prod5_trans_80_zen40az180_dl1ab_tuned_psf/lstchain_config.json

DL2 data[edit]

/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/BLLac/v09x/tailcut84/DL2/real/src_indep/{Date} The DL2 files were trained using RF model, with zenith 20 deg, for runs with zenith pointing less than 30 deg, and with RF model of zenith 40 deg, for otherwise

DL3 data selection[edit]

For energy-dependent cuts[edit]

Information about your DL3 data selection.

  • intensity > 0, 50
  • r: [0, 1]
  • wl: [0.01, 1]
  • leakage_intensity_width_2: [0, 1]
  • source-indep
    • gh_efficiency: 0.8
    • theta_containment: 0.8

Off regions Run 5552 5559 wo safe en mask 1 off.png

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

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 10 TeV with 5 bins per decade

For spectral analysis,

  • Model used - Log Parabola, Smooth Broken Power Law
  • Reference energy used by calculating the decorrelation energy with a Power Law fit, 149 GeV (without intensity cut) and 142 GeV (with intensity cut)
  • Spectral Fit energy range - 10 GeV to 10 TeV with 5 bins per decade
  • Spectral fit with gammapy default optimization - minuit

For LC,

  • E_min = 39.81, 63.1, 100 Gev
  • Also plotted with time bins of 1/3 of each GTI

Analysis Results[edit]

Theta2 plot[edit]

Significance map[edit]

Excess map[edit]

Spectral results[edit]

For energy-dependent cuts[edit]

  • Without intensity cut, Log Parabola

Preliminary LC and SED

Model parameters Run 5552 5559 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png

LC Run 5552 5559 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night 1.png LC Run 5552 5559 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night 2.png LC Run 5552 5559 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png

SED Run 5552 5559 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png


LC with 1/3 GTI time bins

Model parameters Run 5552 0 5559 2 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare short time scale 2.png

LC Run 5552 0 5559 2 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare short time scale 3.png LC Run 5552 0 5559 2 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare short time scale 2.png LC Run 5552 0 5559 2 lp e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare short time scale 1.png

  • Without intensity cut, Smooth Broken Power Law

Preliminary LC and SED

Model parameters Run 5552 5559 sbpl e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night.png

LC Run 5552 5559 sbpl e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night.png SED Run 5552 5559 sbpl e ref 149GeV fit range w 5 bin p decade 10GeV 10TeV flare night.png


  • With intensity cut, Log Parabola

Model parameters Run 5552 5559 lp e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png

LC Run 5552 5559 lp e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 1.png LC Run 5552 5559 lp e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 2.png LC Run 5552 5559 lp e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png

SED Run 5552 5559 lp e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 3.png

  • With intensity cut, Smooth Broken Power Law

Model parameters Run 5552 5559 sbpl e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 1.png

LC Run 5552 5559 sbpl e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 1.png SED Run 5552 5559 sbpl e ref 142GeV fit range w 5 bin p decade 10GeV 10TeV flare night 1.png

Retrieved from "https://www.lst1.iac.es/wiki/index.php?title=Analysis_BLLac_Aug_2021_Priyadarshi&oldid=5344"