PKS1424+240 (LST) analysis - neutrino candidate

Analysis of LST-1 PKS 1424+240 - neutrino candidate

General information

Name of the source: PKS 1424+240

Brief description of the source:

• It is the farthest BL Lac detected at TeV
• It is an IBL/HBL
• z ~ 0.6
• RA, Dec (deg) : 216.75163247, 23.80001042

Analysis by Name (Institute - email)

• Davide Cerasole
• Gloria Maria Cicciari (gloriamaria.cicciari@unipa.it)
• Sara Fogliaceo
• Serena Loporchio
• Manuela Mallamaci
• Salvatore Mangano
• Giovanni Marsella
• Elisa Prandini
• Ilaria Viale
• Elisa Visentin

Data-taking information

Cycle I observations 2023-01-23 to 2023-07-17

• Dates of data-taking: 2023/01/23, 2023/01/27, 2023/03/10, 2023/04/18, 2023/05/12, 2023/05/15, 2023/05/22, 2023/07/17
• Observation runs selected are: [11805 11806 11807 11808 11933 11934 11935 11936 11937 12747 12748 12749 12750 12957 12958 12959 12960]
• Total livetime of all observations (hrs): 19.7
• Total livetime of selected observations (hrs): 3.29
• ZD range (deg): 0 - 60

Analysis steps

• R0 : uncalibrated RAW waveforms from the camera
• R1 : calibrated waveforms from the camera
• DL1a : integrated charge and peak position of the waveform
• DL1b : image parameters (width, length, intensity…)
• DL2 : event parameters (energy, direction, time…)
• DL3 : lists of reconstructed events after event selection with corresponding IRFs

Analysis information

Data quality

The list with the runs passing data quality

Monte Carlo simulations

/fefs/aswg/data/mc/DL2/AllSky 20240131_allsky_v0.10.5_all_dec_base/TestingDataset/dec_2276/

DL1 data

/fefs/aswg/data/real/DL1/YYYYMMDD/v0.10/tailcut84/dl1_LST-1.RunXXXXX.h5

Random forest

/fefs/aswg/data/models/AllSky/20240131_allsky_v0.10.5_all_dec_base/

DL2 data

 /fefs/aswg/workspace/gloriamaria.cicciari/DL3/PKS1424/

DL3 data and IRF

 /fefs/aswg/workspace/gloriamaria.cicciari/DL3/PKS1424/

Analysis results

Preliminary Theta2 plots

Theta2 plots obtained from DL2 data using the explore_DL2.ipynb code

Energy range: 0.0 - 0.2 TeV
Excess: 690.000; Off: 155056
Gamma rate: 0.675 events / minute
Off rate: 151.779 events / minute
alpha (backg normalization): 1.000
Li & Ma Significance: 1.24 standard deviations

Energy range: 0.2 - 1000000.0 TeV
Excess: -6.333; Off: 415
Gamma rate: -0.006 events / minute
Off rate: 0.406 events / minute
alpha (backg normalization): 0.333
Li & Ma Significance: 0.00 standard deviations

Inspect the DL3 file content

• Spatial distribution of events (typically on a plane of celestial coordinates such as RA/Dec or local altitude/azimuth coordinates)
• Energy distribution (spectrum of events)
• Temporal distribution (histogram of events over time)
• Angular offset distribution (angular distance from the centre of the field of view)

Selection filters for the observations

Filters applied to select just a subset of the observation list based on source name, zenith angle or livetime

max_zen = 70 # in deg for a maximum limit on zenith pointing of observations min_time = 300 # in seconds for minimum livetime of each observation

Theta2 plots

Theta2 plots obtained from DL3 data using the post_DL3_analysis.ipynb code

Excess plots

Excess plots obtained from DL3 data using the post_DL3_analysis.ipynb code

SED and LC

SED and LC obtained from DL3 data using the post_DL3_analysis.ipynb code

SED and LC preliminaries that only give us the upper limits. (We reserve the use of a more polished code for analysis).

The analysis will need to be fine-tuned because of the lack of detection.