Enhancing the Cherenkov Telescope Array Observatory high-level performance through an event-type-based analysis
Juan Bernete, Silvia García-Soto, Tarek Hassan, Orel Gueta, Maximilian Linhoff, Atreyee Sinha, Gernot Maier
astro-ph.IM, astro-ph.HE
Submitted: 2026-07-10
Comments: 37 pages, 12 figures
Journal ref: Astroparticle Physics, Volume 182, 2026, page 103227
DOI: 10.1016/j.astropartphys.2026.103277
Code: https://github.com/cta-observatory/iact_event_types
License: http://creativecommons.org/licenses/by-nc-sa/4.0/
The gist: The analysis traditionally employed by Imaging Atmospheric Cherenkov Telescopes involves optimizing quality cuts to select a sub-sample of high-quality events.
Terminology
Abstract
The analysis traditionally employed by Imaging Atmospheric Cherenkov Telescopes involves optimizing quality cuts to select a sub-sample of high-quality events. These events are used for the scientific interpretation of the data, employing a single set of Instrument Response Functions (IRFs). All selected events are treated equally and assumed to be well represented by these IRFs, while the rest are discarded. An alternative approach, successfully applied in experiments such as Fermi-LAT, is an event-type-based analysis. This method divides datasets into subsamples, each containing events of a given expected reconstruction quality. IRFs are computed for each subsample independently, improving the accuracy with which IRFs represent the reconstruction quality of each event. The high-level analysis of these subsamples is performed treating them as independent observations, each with their own set of IRFs, and analyzed jointly. In this work we present a proof-of-concept implementation of an event-type-based analysis for the future CTAO using simulated data. A neural network (specifically a multi-layer perceptron) is trained to predict the direction reconstruction error of each event, and the simulated dataset is divided into event types based on this predicted variable. We compute IRFs for each event type and compare them with those from the standard analysis (without event types). Finally, we simulate observations using these event-type-wise IRFs and analyze them with high-level analysis tools to test the performance of both approaches. This implementation demonstrates notable improvements: 25% to 50% boost in spatial resolving power and 25% in sensitivity. This boost in performance will have strong implications in the scientific exploitation of the CTAO data, especially in crowded regions such as the Galactic Plane or searching for spectral signatures like Dark Matter annihilation lines.
Sources
- GammaLib and ctools: A software framework for the analysis of astronomical gamma-ray data
- Pass 8: Toward the Full Realization of the Fermi-LAT Scientific Potential
- Performance update of an event-type based analysis for the Cherenkov Telescope Array
- Simulation of Imaging Atmospheric Cherenkov Telescopes with CORSIKA and sim_telarray
- A Hybrid Approach to Event Reconstruction for Atmospheric Cherenkov Telescopes Combining Machine Learning and Likelihood Fitting
- Gamma-Hadron Separation Methods for the VERITAS Array of Four Imaging Atmospheric Cherenkov Telescopes
Related papers
- A signal dedispersion algorithm for imaging-based transient searches
- AVICA: A fully automated CASA pipeline for large volume VLBI data calibration
- Spectral Map Making with SPHEREx
- Long-Integration Magnetar Burst Observatory (LIMBO): Instrument Summary and Early FRB Rate Constraints
- Towards independent event horizon imaging of the supermassive black holes in M87 and the Milky Way
- A PINK update: Improvements to the CELEBI fast radio burst data reduction and analysis pipeline