Analysis BLLac Aug 2021 Aguasca-Cabot

Basis

The basis of this algorithm is explained in the following LST-Reco calls:

• Content: Description of the algorithm and proof of concept with Crab data and simulations of different temporal profiles.

- https://indico.cta-observatory.org/event/3571/contributions/29877/attachments/19821/27495/Aguasca_Cabot_12-07-21.pdf

• Content: Upgrade of the algorithm and application to simulations of an orbital period of HESS J0632+057.

- https://indico.cta-observatory.org/event/3659/contributions/30756/attachments/20114/27954/20211104_Aguasca_Cabot.pdf

(Gabriel Emery developed a similar algorithm in his PhD, see https://tel.archives-ouvertes.fr/tel-02983041/document)

BL Lac light curve

- Use the same DL3 data of BL Lac obtained by the BL Lac team in order to use their light curve as the reference one. Apply the algorithm to BL Lac to obtain the minimum time interval for each flux.

Fast analysis

This fast analysis aims to show preliminary results of the algorithm to BL Lac data. MC gammas are optimized to BL Lac (see data reduction specifications). The results should improve with a refined data reduction.

Parameters values of the algorithm:

• Range of energy (TeV): 0.1 to 1
• Threshold in Li & Ma TS: 9
• Threshold for UL (sigma): 2
• Increase of events in each iteration: 1000
• Minimum threshold to consider two consecutive observations (min): 10

The best spectral model obtained fitting all data:

This model is used to compute the fluxes using forward folding

Data reduction specifications

• 

  DL1 to DL3 specifications
• 

  High-level specifications

Light curve

We can see that the flux points produced by the algorithm are well distributed around the run-wise fluxes and fluxes for a time interval of more or less three minutes. Also, we can see that when the 3 min fluxes increase or decrease, the 3 sigma fluxes also show an increase or decrease. For example, see fluxes at t=5.9435e4+0.05 MJD.

TS detection and time interval of the fluxes

• File:TS-20210808-fast-TRETS-loosecuts.png

  TS (Li & Ma) for each time interval where a flux is plotted.
• 

  Histogram of flux time intervals.

In the left figure, some TS values are much higher than the threshold value (9) because I used an increase of 1000 events in each loop to reduce the compute time.

In the right figure, there are high values of time intervals because events from two different runs are used to compute the flux.

The p-value against the null hypothesis of a constant source is 1.7350606859473474e-101.