Difference between revisions of "BLLac2022flare"
Jorge.otero (talk | contribs) (→Lightcurve) |
Jorge.otero (talk | contribs) (→Lightcurve) |
||
| Line 264: | Line 264: | ||
Flux value for Emin = 150 GeV | Flux value for Emin = 150 GeV | ||
| − | *Average flux: f = (4. | + | *Average flux: f = (4.5 +/- 0.2)e-10 cm-2 s-1 = 1.5 C.U. |
*Maximum flux: f = (9.1 +/- 0.8)e-10 cm-2 s-1 = 3.1 C.U. | *Maximum flux: f = (9.1 +/- 0.8)e-10 cm-2 s-1 = 3.1 C.U. | ||
Revision as of 08:58, 1 July 2024
Contents
- 1 BL Lacertae 2022 flare project
BL Lacertae 2022 flare project
Major flare from BL Lacertae between September and November 2022, one of the (if not the) longest VHE flare observed from the source, and the brightest together with the 2021 flare.
We aim to characterize the VHE and MWL variability during the flare, both in flux and spectral variability, with focus on the intranight minute timescale variations in the VHE band, and to perform a sort of time dependent modelling of the emission to interpret the evolution of the flare in terms of physical mechanisms and paramters.
General information
- Name of the source: BL Lacertae
- Brief description of the source:
- - Object type : Blazar
- - RA: 22 02 43.3 (hh mm ss), Dec: +42 16 40 (dd mm ss)
- Analysis by:
- Jorge Otero Santos (IAA-CSIC, joteros@iaa.es)
- Daniel Morcuende (IAA-CSIC, dmorcuende@iaa.es)
Data-taking information
- Dates of data-taking:
- Date: 16/09/2022-25/11/2022 (53 h before cuts)
- ZD range (deg): 15-60 deg
- AZ range (deg): 0-60 and 300-360 deg
- Joint observations with MAGIC? : No
- Runs:
- Date: 16/09/2022-25/11/2022 (53 h before cuts)
1 : 2022-09-16 : 9228 to 9232 2 : 2022-09-21 : 9262 to 9268 3 : 2022-09-22 : 9347 to 9350 4: 2022-10-01: 9453 to 9461 5: 2022-10-03: 9544 to 9555 6: 2022-10-15: 9631 to 9635 7: 2022-10-16: 9648 8: 2022-10-17: 9678 to 9685 9: 2022-10-18: 9697 to 9706 10: 2022-10-19: 9729 to 9739 11: 2022-10-20: 9753 to 9761 12: 2022-10-21: 9808 to 9815 13: 2022-10-23: 9836 to 9848 14: 2022-10-24: 9864 to 9873 15: 2022-10-25: 9892 to 9900 16: 2022-10-26: 9921 to 9924 17: 2022-10-27: 9940 to 9944 18: 2022-10-28: 9969 to 9973 19: 2022-10-29: 9993 to 9995 20: 2022-10-31: 10022 to 10029 21: 2022-11-01: 10069 to 10076 22: 2022-11-02: 10193 to 10204 23: 2022-11-13: 10508 to 10512 24: 2022-11-14: 10521 to 10523 25: 2022-11-15: 10586 to 10588 26: 2022-11-16: 10616 to 10620 27: 2022-11-17: 10652 to 10655 28: 2022-11-18: 10827 to 10830 29: 2022-11-19: 10870 to 10872 30: 2022-11-25: 11033 to 11037
Data-quality selection
Used latest version of the data quality selection notebook.
Standard configuration + dR/dI cosmics rate @ 422 p.e. = 1.2
Cut of max_diffuse_nsb_std = 2.3 for moon data selection
- Selected dark runs: 104 of 193 runs (29 h)
1 : 2022-09-21 : 9262 to 9268 2: 2022-10-01: 9459 to 9461 3: 2022-10-03: 9544 4: 2022-10-16: 9648 5: 2022-10-17: 9678 to 9683 6: 2022-10-18: 9698 to 9705 7: 2022-10-19: 9730 to 9739 8: 2022-10-20: 9753 to 9761 9: 2022-10-21: 9809 to 9814 10: 2022-10-23: 9843, 9844, 9845, 9846, 9848 11: 2022-10-24: 9866 12: 2022-10-25: 9892 to 9900 13: 2022-10-26: 9921 to 9924 14: 2022-10-27: 9940 to 9944 15: 2022-10-29: 9993 to 9995 16: 2022-10-31: 10025, 10026, 10028, 10029 17: 2022-11-13: 10508 18: 2022-11-15: 10586 to 10587 19: 2022-11-16: 10616 to 10620 20: 2022-11-17: 10652 to 10655 21: 2022-11-18: 10827 to 10830 22: 2022-11-19: 10871 to 10872 23: 2022-11-25: 11033 to 11037
- Selected moon runs: 30 of 193 runs (8 h)
1: 2022-10-01: 9453 to 9457 2: 2022-10-03: 9544, 9545, 9546, 9548, 9549, 9550, 9551, 9552 3: 2022-10-15: 9631, 9633, 9634, 9635 4: 2022-10-17: 9684 5: 2022-10-31: 10022 to 10023 6: 2022-11-01: 10071, 10072, 10073, 10074, 10076 7: 2022-11-13: 10509 to 10512 8: 2022-11-14: 10523
Monte Carlo information
- Link to MC files used:
- Dark data: /fefs/aswg/data/mc/DL2/AllSky/20230925_v0.10.4_src3_dec3476_4822_tuned/TestingDataset/dec_3476/
- Moon data: /fefs/aswg/data/mc/DL2/AllSky/20240426_v0.10.9_src3_dec3476_4822_tunedNSB/TestingDataset/dec_3476/ (tuned NSBs and cleaning to match that needed for the real data)
DL1 data
Dark data :
Standard data from /fefs/aswg/data/real/DL1/ with version v0.9 and cleaning tailcut84
Moon data :
Moon data reprocessed to produce the DL1b files with the optimal cleaning for moon analysis as calculated with lstchain_dvr_pixselector -f "/fefs/aswg/data/real/DL1/2022XXXX/v0.9/tailcut84/dl1_LST-1.RunXXXXX'.*.h5" Reprocessed using DL1a files produced by LSTOSA (lstchain v0.10.XX) and dl1ab script (v0.7.5)
- original DL1a files
/fefs/aswg/data/real/DL1/2022XXXX/v0.9/tailcut84/dl1_LST-1.Run0XXXX.XXXX.h5
- Produced DL1b files
/fefs/aswg/workspace/daniel.morcuende/data/real/DL1b/BLLac_2022_NSB/tailcut_*_*/merged
tailcut105: 09454, 09455, 09456, 09457, 09631, 09684 tailcut126: 09453, 09551, 09552, 10022, 10023, 10509, 10523 tailcut147: 09546, 09548, 09549, 09550, 09633, 09634, 09635, 10071, 10072, 10073, 10074, 10076, 10510, 10511 tailcut168: 09544, 09545, 10512
Low-level analysis
- Interpolated IRFs
- Custom MC and models (as explained above for dark/moon data)
- Standard 70% efficiency energy dependent cuts
- Max theta cut = 0.2
- Intensity dark = 50
- Intensity moon = 100
High-level analysis
- lstchain v0.10.11 to generate source-indep interpolated IRF, energy dependent cuts DL3
- Science Tool: gammapy 1.1
- point-like IRF, 1D analysis
- Aeff>5%
- 3OFF
Energy threshold
Calculated from the dark MCs up to the Zd of BL Lac observations (evaluated from the MCs between 30 and 60 deg): ~80 GeV
Calculated from the moon NSB-tuned MCs up to the Zd of BL Lac observations (evaluated from the MCs between 30 and 60 deg , to be re-evaluated with tuned-NSB+tuned-cleaning MCs): ~100 GeV
NOTES ON THE ENERGY THRESHOLD CALCULATION:
- For the dark data, the distribution is calculated in a rather conservative way, with MCs >32deg Zd.
- For moon data:
- The distribution is calculated above 10deg (~minimum of BL Lacs observations)
- Above 30deg the Eth is ~130 GeV
- Only 9 runs with moon above 40 deg zenith, most of the data (except ~3 hours) described by an Eth<100-130 GeV
- 21 runs remaining <30deg Zd, where the Eth is 130 GeV or less (~100 GeV for the lowest zeniths)
Jorge's Analysis Results
We produce the stacked and daily theta2 plots and average SED.
With the average SED, a preliminary LC is also computed and a Bayesian Block analysis. We then calculate the SED of each BB to look for spectral variability (see below the BB definition).
For a more precise LC, we then use the SED of each block to calculate the LC block by block to build the long-term final LC. This is to avoid miscalculations in the integral flux caused by using the wrong spectral shape if there is strong spectral variability.
We take as reference the energy threshold estimated from the moon MCs for the minimum energy in which we integrate the flux for the LC. This energy threshold is specified above.
Configuration of the analysis:
- E_reco=[20,20000] GeV
- E_true=[10,50000] GeV
- n_bins_reco=5-10 (depending on the flux level)
- n_bins_true=10-20 (depending on the flux level)
- E_fit=[20,20000] GeV
- Spectral shape: Log parabola+EBL (Saldaña-Lopez21, z=0.069)
- Analysis type: stack
- LC energy: 3 different Emin tested to account for possible threshold increase due to moon
- E=[100,20000] GeV
- E=[150,20000] GeV
- E=[200,20000] GeV
The rest of the details of the high-level analysis are as specified above.
Theta2 plot
All data stacked (>20 GeV left, >100 GeV right)
20 October >100 GeV (left), 13 November >100 GeV (right)
Lightcurve
Complete LC calculated as detailed above (using the SED of each BB, shown below, for the flux integration). We test different energy mins as specified above.
Emin 100 GeV (left), 150 GeV (center) and 200 GeV (right)
Flux value for Emin = 100 GeV (average flux calculated only from significant flux points, i.e. TS>4)
- Average flux:
- Maximum flux (nightwise): (9.9 +/- 0.2)e-10 cm-2 s-1 = 2.1 C.U.
Flux value for Emin = 150 GeV (average flux calculated only from significant flux points, i.e. TS>4)
- Average flux:
- Maximum flux (nightwise): (4.7 +/- 0.2)e-10 cm-2 s-1 = 1.6 C.U.
Flux value for Emin = 200 GeV (average flux calculated only from significant flux points, i.e. TS>4)
- Average flux:
- Maximum flux (nightwise): (2.7 +/- 0.1)e-10 cm-2 s-1 = 1.4 C.U.
We also calculate the intranight lightcurves for the nights of Oct. 20 and Nov. 13, 2022.
For the former, given that it is all dark data, we use an energy threshold of 100 GeV, and we compare a 5-min (left) and a 10-min binning (right).
- Average flux: f = (9.9 +/- 0.2)e-10 cm-2 s-1 = 2.1 C.U.
- Maximum flux: f = (2.0 +/- 0.2)e-9 cm-2 s-1 = 4.3 C.U.
For the latter we test as Emin 100 GeV, 150 GeV and 200 GeV due to the moon, with a 5-min binning.
Emin 100 GeV (left), 150 GeV (center) and 200 GeV (right)
Flux value for Emin = 100 GeV
- Average flux: f = (8.8 +/- 0.3)e-10 cm-2 s-1 = 1.9 C.U.
- Maximum flux: f = (1.8 +/- 0.1)e-9 cm-2 s-1 = 3.8 C.U.
Flux value for Emin = 150 GeV
- Average flux: f = (4.5 +/- 0.2)e-10 cm-2 s-1 = 1.5 C.U.
- Maximum flux: f = (9.1 +/- 0.8)e-10 cm-2 s-1 = 3.1 C.U.
Flux value for Emin = 200 GeV
- Average flux: f = (2.7 +/- 0.1)e-10 cm-2 s-1 = 1.4 C.U.
- Maximum flux: f = (5.1 +/- 0.5)e-10 cm-2 s-1 = 2.5 C.U.
Spectral results
Average SED (parameters correspond to the EBL-corrected spectrum using the EBL model from Saldaña-Lopez 2021)
- f0 = (4.28 +/- 0.17)e-13 GeV-1 cm-2 s-1
- E0 = 223 GeV
- alpha = 2.97 +/- 0.04
- beta = 0.13 +/- 0.04
- alpha_norm = 1
- redshift = 0.069
We calculate also the SED of each Bayesian Block (also considering a log-parabola plus EBL). BB4 and BB8 correspond to the nights of Oct. 20 and Nov. 13, 2022, respectively.
BB1 (left) and BB2 (right)
BB3 (left) and BB4 (right)
BB5 (left) and BB6 (right)
BB7 (left) and BB8 (right)
BB9 (left) and BB10 (right)
Daniel's Analysis Results
Relevant details
Theta2 plot
Lightcurve
Spectral results
Crosschecked plots
Lightcurve
Complete LC calculated as detailed above (using the SED of each BB, shown below, for the flux integration). We test different energy mins as specified above.
Emin 100 GeV (left), 150 GeV (center) and 200 GeV (right)
We also calculate the intranight lightcurves for the nights of Oct. 20 and Nov. 13, 2022.
For the former, given that it is all dark data, we use an energy threshold of 100 GeV, and we compare a 5-min (left) and a 10-min binning (right).
For the latter we test as Emin 100 GeV, 150 GeV and 200 GeV due to the moon, with a 5-min binning.
Emin 100 GeV (left), 150 GeV (center) and 200 GeV (right)
Spectral results
Average SED (parameters correspond to the EBL-corrected spectrum using the EBL model from Saldaña-Lopez 2021)
We calculate also the SED of each Bayesian Block (also considering a log-parabola plus EBL). BB4 and BB8 correspond to the nights of Oct. 20 and Nov. 13, 2022, respectively.
BB1 (left) and BB2 (right)
BB3 (left) and BB4 (right)
BB5 (left) and BB6 (right)
BB7 (left) and BB8 (right)
BB9 (left) and BB10 (right)
Bayesian blocks
1. 21/09/2022 22:50:00 to 04/10/2022 00:57:00 2. 15/10/2022 23:25:00 to 18/10/2022 01:37:00 3. 18/10/2022 22:09:00 to 20/10/2022 01:57:00 4. 20/10/2022 21:56:00 to 21/10/2022 00:39:00 5. 21/10/2022 22:57:00 to 26/10/2022 23:38:00 6. 27/10/2022 23:37:00 to 28/10/2022 01:00:00 7. 29/10/2022 23:44:00 to 03/11/2022 00:11:00 8. 13/11/2022 22:26:00 to 13/11/2022 23:47:00 9. 14/11/2022 23:18:00 to 20/11/2022 00:04:00 10. 25/11/2022 21:45:00 to 25/11/2022 22:30:00
