Analysis BLLac Aug 2021 w IRF interpolation Priyadarshi

From my_wiki
Jump to: navigation, search

BL Lacertae August 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]

See https://www.lst1.iac.es/wiki/index.php/BLLac#Data-taking_information

Monte Carlo information[edit]

For getting a closer data-MC comparison, we produce MC data along various sky positions (nodes) as a preliminary method to facilitate 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[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 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[edit]

  • Old RF

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/


  • New RF

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 and similarly for "4822". We have combined the Training nodes of the 2 declination lines to have a better RF model for the declination path of the source NGC 1275. See All_sky_mc_bllac for more information.

In this new RF model production, we include zenith and azimuth pointing parameters as well.

  • lstchain v0.9.13dev source-independent, dynamic cleaning, NSB tuning on BL Lac,
 /fefs/aswg/workspace/seiya.nozaki/BLLac_paper/20221104_src3_dec3476_4822_tuned/merged/srcindep/RF/trained_models/

For IRF interpolation[edit]

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:
    • Old files

/fefs/aswg/data/mc/DL2/AllSky/20220601_dec3476_tuned/TestingDataset/

    • New files

/fefs/aswg/workspace/seiya.nozaki/BLLac_paper/20221104_src3_dec3476_4822_tuned/merged/srcindep/DL2_MC/data/

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

Source pointings Run 5443 5706.png

DL2 data[edit]

The DL2 files of the real data are stored here -

  • Old files

/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/BLLac/v09x/tailcut84/DL2/real/src_indep/20220601_dec3476/{Date}

  • New files

/fefs/aswg/workspace/chaitanya.priyadarshi/real_data/AGN/BLLac/v09x/tailcut84/DL2/real/src_indep/20221104_dec3476_4822_tuned/{Date}

DL3 data selection[edit]

For energy-dependent cuts[edit]

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

Off regions Run 5443 5706 wo safe en mask 1 off.png Off regions Run 5443 5706 w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

Gammaness cuts all Run 5443 5706.png Rad max all.png

High-level analysis[edit]

  • 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 and 2 Wobbles
  • Safe Energy mask for systematics - Yes, 25.1 GeV < E < 40 TeV
  • Energy axes -
    • True: 10 GeV to 100 TeV with 10 bins per decade
    • Reco: 10 GeV to 10 TeV with 5 (10 for Aug 3, 8) bins per decade
  • DL4 dataset plots for full August (1 and 2 OFF regions)

DL4 significance excess plots Run 5443 5706 w safe en mask min 25.1 GeV 1 off 10 energy bin p dec.png DL4 significance excess plots Run 5443 5706 w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

  • DL4 dataset plots for Aug 3 (1 and 2 OFF regions)

DL4 significance excess plots Run 5454 5466 w safe en mask min 25.1 GeV 1 off 10 energy bin p dec.png DL4 significance excess plots Run 5454 5466 w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

  • DL4 dataset plots for Aug 8 (1 and 2 OFF regions)

DL4 significance excess plots Run 5552 5559 w safe en mask min 25.1 GeV 1 off 10 energy bin p dec.png DL4 significance excess plots Run 5552 5559 w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

For spectral analysis,

  • Model used - Log Parabola + EBL (Dominguez)
  • Reference energy used by calculating the decorrelation energy with a Power Law fit, 124 (82 for Aug 3, 125 for Aug 8) GeV
  • Spectral Fit energy range - 10 GeV to 10 TeV with 5 (10 for Aug 3, 8) bins per decade
  • Spectral fit with gammapy default optimization - minuit

For LC,

  • E_min = 63.1, 100 Gev

Analysis Results[edit]

Theta2 plot[edit]

Significance map[edit]

Excess map[edit]

Spectral results[edit]

For energy-dependent cuts[edit]

  • Full Aug dataset
    • 1 OFF region

Preliminary LC and SED

Model parameters correlation Run 5443 5706 lp ebl dominguez e ref 125GeV fit range w 5 bin p decade 10GeV 10TeV.png

Model parameters Run 5443 5706 lp ebl dominguez e ref 125GeV fit range w 5 bin p decade 10GeV 10TeV.png

LC min 63 Run 5443 5706 lp ebl dominguez e ref 125GeV fit range w 5 bin p decade 10GeV 10TeV.png

LC min 100 Run 5443 5706 lp ebl dominguez e ref 125GeV fit range w 5 bin p decade 10GeV 10TeV.png

SED Run 5443 5706 lp ebl dominguez e ref 125GeV fit range w 5 bin p decade 10GeV 10TeV.png

    • 2 OFF regions (10 energy bins per decade)

Model parameters correlation Run 5443 5706 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

Model parameters Run 5443 5706 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

LC min 100 max 10000 Run 5443 5706 w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

SED 2 Run 5443 5706 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png


  • Aug 3 night
    • 2 OFF regions

Model parameters correlation Run 5454 5466 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

Model parameters Run 5454 5466 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

LC min 63 max 10000 Run 5454 5466 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

LC min 100 max 10000 Run 5454 5466 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

LC min 100 max 100000 Run 5454 0 5466 2 lp ebl dominguez w safe en mask min 25.1 GeV w short gti obs w 3 gti bins 2 off 10 energy bin p dec.png

SED Run 5454 5466 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

  • Aug 8 night
    • 1 OFF regions

Preliminary LC and SED

Model parameters correlation Run 5552 5559 lp ebl dominguez e ref 124GeV fit range w 10 bin p decade 10GeV 10TeV w safe mask.png

Model parameters Run 5552 5559 lp ebl dominguez e ref 124GeV fit range w 10 bin p decade 10GeV 10TeV w safe mask.png

LC min 63 Run 5552 5559 lp ebl dominguez e ref 124GeV fit range w 10 bin p decade 10GeV 10TeV w safe mask.png SED Run 5552 5559 lp ebl dominguez e ref 124GeV fit range w 10 bin p decade 10GeV 10TeV w safe mask.png

    • 2 OFF regions

Preliminary LC and SED

Model parameters correlation Run 5552 5559 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

Model parameters Run 5552 5559 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

LC here is built with standard GTI time intervals and with 1/3rd of each GTI time interval

LC min 63 max 10000 Run 5552 5559 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png LC min 63 max 10000 Run 5552 0 5559 2 lp ebl dominguez w safe en mask min 25.1 GeV w short gti obs w 3 gti bins 2 off 10 energy bin p dec.png LC min 100 Run 5552 5559 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png LC min 100 max 10000 Run 5552 0 5559 2 lp ebl dominguez w safe en mask min 25.1 GeV w short gti obs w 3 gti bins 2 off 10 energy bin p dec.png

SED Run 5552 5559 lp ebl dominguez w safe en mask min 25.1 GeV 2 off 10 energy bin p dec.png

  • Comparison of different OFF regions for Aug 8 dataset

LC comparison flare 8Aug 2021 3 gti bins diff off regions emin 100GeV.png

SED comparison flare 8Aug 2021 diff off regions emin 100GeV.png

  • Comparison of different OFF regions for Aug 3 dataset

SED comparison flare 3Aug 2021 diff off regions emin 100GeV.png

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