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 probe of concept in Crab data and simulations.

- 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 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. Neither the cuts nor the MC simulations are optimized to BL Lac. The results should improve with a refined data reduction and proper parameter values for the algorithm, for example, the number of events added in each iteration.

Events

• Range of energy: 0.1 to 1 TeV
• Threshold in Li & Ma significance: 3
• Threshold for UL: 2 sigma
• Increase of events in each iteration: 50
• Minimum threshold to consider two consecutive observations: 10 min

Data reduction specifications

• 

  DL1 to DL3 specifications
• 

  High-level specifications

Light curve

TS detection and time interval of the fluxes

• 

  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 because I did not use soft cuts and the number of events added for each iteration was a bit high to reduce computation time.

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