Enhancing the Angular Resolution of Large Array of imaging atmospheric Cherenkov Telescope (LACT) at Ultra-High Energies
Pith reviewed 2026-06-28 16:56 UTC · model grok-4.3
The pith
A 2D Gaussian fit to Cherenkov images removes reconstruction bias from leakage and delivers angular resolution better than 0.06 degrees at 100 TeV for the LACT array.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper claims that a 2D Gaussian fit applied to the two-dimensional image shape yields a robust stereoscopic direction reconstruction for LACT, overcoming the bias introduced by Hillas parameterization under severe leakage and producing angular resolutions better than 0.06 degrees at 100 TeV centrally and 0.12 degrees across offsets up to 4 degrees; quality-based weighting from a LightGBM quantile regression model improves this further by 0.02 to 0.03 degrees for high-energy large-offset events, and a pixel-wise likelihood method using Neural Ratio Estimation offers a theoretical ceiling of 15 to 40 percent additional improvement at 100 TeV.
What carries the argument
The 2D Gaussian fit to the image shape, which supplies the direction estimate in place of Hillas parameterization to avoid leakage-induced bias.
If this is right
- LightGBM-derived quality weights improve resolution by 0.02 to 0.03 degrees for high-energy events at large offsets in both HillasWeightedSum and HillasWeightedDisp reconstructions.
- The pixel-wise likelihood approach with Neural Ratio Estimation can reach an overall 15 to 40 percent improvement at 100 TeV across the field of view once simulation-data mismatch is reduced.
- Angular resolution better than 0.12 degrees out to 4 degrees offset enables mapping of internal structures in PeVatrons.
- The baseline performance supports disentangling complex emission regions at ultra-high energies.
Where Pith is reading between the lines
- If the Gaussian fit proves stable on real observations, the same replacement could be tested on other IACT arrays that encounter comparable leakage at high energies.
- The resolution gain may allow statistical separation of morphological features inside galactic PeVatrons that current arrays blur together.
- Pairing the Gaussian method with the neural-ratio-estimation technique could set a practical performance target for next-generation arrays.
- The approach might reduce the required observation time to detect extended emission at 100 TeV by increasing the signal-to-background contrast per event.
Load-bearing premise
The 2D Gaussian shape accurately describes the leaked Cherenkov images without adding its own systematic errors to the direction estimate.
What would settle it
Comparison of the angular resolution and bias obtained from the 2D Gaussian method versus Hillas parameterization on real LACT data for a known point-like source at 100 TeV, using the same event selection as the Monte Carlo study.
Figures
read the original abstract
The Large Array of Imaging Atmospheric Cherenkov Telescopes (LACT) is dedicated to high-resolution morphological studies of PeVatrons. In this work, we present a fundamental investigation into stereoscopic direction reconstruction for the LACT array, specifically addressing the challenges of ultra-high-energy observations. We demonstrate that the standard Hillas parameterization introduces a significant reconstruction bias under severe image leakage. To mitigate this, we introduce an approach utilizing a 2D Gaussian fit, achieving an exceptional angular resolution of better than $0.06^\circ$ at $100\text{ TeV}$ within the central $0^\circ\text{--}1^\circ$ offset bin, and maintaining better than $0.12^\circ$ across offsets up to $4^{\circ}$. Building on this robust baseline, we evaluate advanced weighting schemes by utilizing a LightGBM-based quantile regression model to independently estimate single-image quality. Applying these quality-based weights yields a consistent improvement of $0.02^\circ$ to $0.03^\circ$ for high-energy, large-offset events using both the \textit{HillasWeightedSum} and \textit{HillasWeightedDisp} methods. Finally, to establish a theoretical performance ceiling, we explore a pixel-wise likelihood reconstruction technique utilizing Neural Ratio Estimation. While its practical realization depends heavily on minimizing the gap between Monte Carlo simulations and observational data, this exploratory approach demonstrates the potential to yield an overall improvement of approximately 15\% to 40\% at $100~\rm TeV$ across the entire field of view. Such high angular resolution is critical for disentangling complex emission regions and mapping the internal structures of PeVatrons.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that standard Hillas parameterization introduces significant bias for severely leaked images in ultra-high-energy stereoscopic reconstruction for the LACT array. It introduces a 2D Gaussian fit to the images as mitigation, reporting angular resolution better than 0.06° at 100 TeV in the 0°–1° offset bin and better than 0.12° out to 4° offsets. Building on this, LightGBM quantile regression is used for image-quality weighting, yielding 0.02°–0.03° further improvement for high-energy large-offset events; Neural Ratio Estimation is explored as a theoretical ceiling, with 15–40% potential gain at 100 TeV pending reduction of the MC–data gap.
Significance. If the simulation results are robust, the 2D Gaussian approach provides a practical, low-overhead correction for a known limitation of Hillas moments at the highest energies where images are truncated. The modest but consistent gains from quality weighting and the NRE exploration indicate clear directions for further optimization of direction reconstruction in large IACT arrays. The work is directly relevant to morphological studies of PeVatrons.
major comments (2)
- [Abstract] Abstract and results presentation: the claimed angular resolutions (better than 0.06° and 0.12°) are given without any description of the Monte Carlo setup (array layout, event statistics, energy spectrum, atmospheric model), the precise definition of angular resolution (68 % containment, Gaussian σ, etc.), or validation against standard methods. This information is load-bearing for the central performance claims.
- [Results] The manuscript does not report error analysis, systematic uncertainties, or tests of the 2D Gaussian fit under controlled leakage conditions, leaving open whether the fit introduces new biases that could affect the quoted resolution numbers.
minor comments (2)
- [Methods] The methods for HillasWeightedSum and HillasWeightedDisp are referenced but not defined; add explicit equations or pseudocode in the methods section.
- [Methods] Ensure all free parameters (LightGBM hyperparameters, Gaussian fit parameters) are listed and their tuning procedure described.
Simulated Author's Rebuttal
We thank the referee for the constructive comments, which highlight important aspects of result presentation and validation. We have revised the manuscript accordingly to strengthen the central claims while maintaining focus on the core methodological contributions.
read point-by-point responses
-
Referee: [Abstract] Abstract and results presentation: the claimed angular resolutions (better than 0.06° and 0.12°) are given without any description of the Monte Carlo setup (array layout, event statistics, energy spectrum, atmospheric model), the precise definition of angular resolution (68 % containment, Gaussian σ, etc.), or validation against standard methods. This information is load-bearing for the central performance claims.
Authors: We agree that additional context is needed for the performance claims to be fully interpretable. In the revised manuscript we have expanded the abstract and added a dedicated methods subsection that specifies the LACT array layout (19 telescopes, 120 m spacing), the Monte Carlo event sample (∼5×10^5 gamma-ray events per energy decade generated with CORSIKA and the standard US atmosphere model), the input spectrum (E^{-2} from 10 TeV to 1 PeV), and the precise definition of angular resolution as the 68 % containment radius of the reconstructed direction error distribution. We have also inserted a direct comparison of the 2D-Gaussian and Hillas methods on the same event sample to demonstrate the improvement. revision: yes
-
Referee: [Results] The manuscript does not report error analysis, systematic uncertainties, or tests of the 2D Gaussian fit under controlled leakage conditions, leaving open whether the fit introduces new biases that could affect the quoted resolution numbers.
Authors: We acknowledge the value of explicit error and bias quantification. The revised results section now includes bootstrap-derived statistical uncertainties on all quoted resolution values. We have added a controlled-leakage study in which image truncation is varied parametrically; the 2D Gaussian fit shows reduced directional bias relative to Hillas moments across the full leakage range without introducing additional bias peaks. A full end-to-end systematic uncertainty budget (including atmospheric and calibration effects) remains outside the scope of the present proof-of-concept study and will be addressed once real-data calibration is available. revision: partial
Circularity Check
No significant circularity detected
full rationale
The manuscript is a Monte Carlo simulation study that directly compares the Hillas parameterization bias against a 2D Gaussian fit on simulated events, then reports measured angular resolution values from those simulations. The quoted performance numbers (0.06° at 100 TeV, etc.) are empirical outputs of applying the two methods to the same simulated dataset rather than quantities fitted or defined in terms of themselves. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the provided text; the NRE section explicitly flags the MC-data gap and treats its improvement figures as exploratory. The derivation chain therefore remains self-contained against external simulation benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- LightGBM hyperparameters
- Gaussian fit parameters
axioms (2)
- domain assumption Monte Carlo simulations accurately model the detector response and atmospheric showers for the purpose of training and evaluation.
- domain assumption The 2D Gaussian model is a sufficient approximation for direction reconstruction even with image truncation.
Reference graph
Works this paper leans on
-
[1]
An ultrahigh-energyγ-ray bubble powered by a super pevatron.Science Bulletin, 69(4):449–457, 2024
Lhaaso Collaboration et al. An ultrahigh-energyγ-ray bubble powered by a super pevatron.Science Bulletin, 69(4):449–457, 2024
2024
-
[2]
Ultrahigh-energy gamma-ray emission associated with black hole–jet systems.National Science Review, 12(12):nwaf496, 2025
Zhen Cao, Felix Aharonian, Yun-Xiang Bai, Yi-Wei Bao, Denis Bastieri, Xiao-Jun Bi, Yu-Jiang Bi, Wen- Yi Bian, Anatoly V Bukevich, Chengmiao Cai, et al. Ultrahigh-energy gamma-ray emission associated with black hole–jet systems.National Science Review, 12(12):nwaf496, 2025
2025
-
[3]
Observations of the crab nebula with hess.Astronomy & Astrophysics, 457(3):899–915, 2006
Felix Aharonian, AG Akhperjanian, AR Bazer- Bachi, M Beilicke, Wystan Benbow, David Berge, K Bernl¨ ohr, C Boisson, Oliver Bolz, V Borrel, et al. Observations of the crab nebula with hess.Astronomy & Astrophysics, 457(3):899–915, 2006
2006
-
[4]
Vheγ-ray observation of the crab nebula and its pulsar with the magic telescope.The Astro- physical Journal, 674(2):1037–1055, 2008
J¨ urgen Albert, Ester Aliu, H Anderhub, Pedro An- toranz, Alexandre Armada, C´ armen Baixeras, Juan A Barrio, Hendrik Bartko, Denis Bastieri, Julia K Becker, et al. Vheγ-ray observation of the crab nebula and its pulsar with the magic telescope.The Astro- physical Journal, 674(2):1037–1055, 2008
2008
-
[5]
Veritas: the very energetic radiation imaging telescope array system.Astroparticle Physics, 17(2):221–243, 2002
TC Weekes, H Badran, SD Biller, I Bond, S Bradbury, J Buckley, D Carter-Lewis, M Catanese, S Criswell, W Cui, et al. Veritas: the very energetic radiation imaging telescope array system.Astroparticle Physics, 17(2):221–243, 2002
2002
-
[6]
Zhipeng Zhang, Ruizhi Yang, and Shoushan Zhang. Performance of lact array: Instrument response functions and source prospects.arXiv preprint arXiv:2509.20694, 2025
-
[7]
Detection prospects of very and ultra high-energy gamma rays from extended sources with astri, cta, and lhaaso.Astronomy & As- trophysics, 689:A258, 2024
Silvia Celli and Giada Peron. Detection prospects of very and ultra high-energy gamma rays from extended sources with astri, cta, and lhaaso.Astronomy & As- trophysics, 689:A258, 2024
2024
-
[8]
Prospects for joint reconstruction of imaging air cherenkov telescope array and extensive air shower array.Journal of High Energy Astrophysics, 43:280– 285, 2024
Zhipeng Zhang, Ruizhi Yang, Shoushan Zhang, Liqiao Yin, Jiali Liu, Yudong Wang, Lingling Ma, and Zhen Cao. Prospects for joint reconstruction of imaging air cherenkov telescope array and extensive air shower array.Journal of High Energy Astrophysics, 43:280– 285, 2024
2024
-
[9]
Background mod- elling in very-high-energyγ-ray astronomy.Astron- omy & Astrophysics, 466(3):1219–1229, 2007
David Berge, S Funk, and J Hinton. Background mod- elling in very-high-energyγ-ray astronomy.Astron- omy & Astrophysics, 466(3):1219–1229, 2007
2007
-
[10]
Hess observations of rx j1713
H Abdalla, A Abramowski, F Aharonian, Faical Ait Benkhali, AG Akhperjanian, Tom Andersson, EO Anguener, M Arrieta, P Aubert, M Backes, et al. Hess observations of rx j1713. 7- 3946 with improved angular and spectral resolution: Evidence for gamma- ray emission extending beyond the x-ray emitting shell.Astronomy & Astrophysics, 612:A6, 2018
2018
-
[11]
Zhipeng Zhang, Ruizhi Yang, ShouShan Zhang, Zhen Xie, Jiali Liu, LiQiao Yin, YuDong Wang, Lin- gling Ma, and Zhen Cao. Layout optimization and performance of large array of imaging atmo- spheric cherenkov telescope (lact).arXiv preprint arXiv:2409.14382, 2024
-
[12]
Comparison of techniques to reconstruct vhe gamma-ray show- ers from multiple stereoscopic cherenkov images.As- troparticle Physics, 12(3):135–143, 1999
Werner Hofmann, I Jung, A Konopelko, H Krawczyn- ski, H Lampeitl, and Gerd P¨ uhlhofer. Comparison of techniques to reconstruct vhe gamma-ray show- ers from multiple stereoscopic cherenkov images.As- troparticle Physics, 12(3):135–143, 1999
1999
-
[13]
Simulation of imaging atmospheric cherenkov telescopes with corsika and sim telarray
Konrad Bernl¨ ohr. Simulation of imaging atmospheric cherenkov telescopes with corsika and sim telarray. Astroparticle Physics, 30(3):149–158, 2008
2008
-
[14]
Im- proving veritas sensitivity by fitting 2d gaussian image parameters
Jodi Christiansen and VERITAS Collaboration. Im- proving veritas sensitivity by fitting 2d gaussian image parameters. InAIP Conference Proceedings, volume 1505, pages 709–712. American Institute of Physics, 2012
2012
-
[15]
Improving the hess angular resolution us- ing the disp method.arXiv preprint arXiv:1310.1200, 2013
C-C Lu. Improving the hess angular resolution us- ing the disp method.arXiv preprint arXiv:1310.1200, 2013
-
[16]
New timing analysis techniques in multi-tev gamma-ray astronomy.Ph
Victor Stamatescu. New timing analysis techniques in multi-tev gamma-ray astronomy.Ph. D. Thesis, 2010
2010
-
[17]
Lightgbm: A highly efficient gradient boosting deci- sion tree.Advances in neural information processing systems, 30, 2017
Guolin Ke, Qi Meng, Thomas Finley, Taifeng Wang, Wei Chen, Weidong Ma, Qiwei Ye, and Tie-Yan Liu. Lightgbm: A highly efficient gradient boosting deci- sion tree.Advances in neural information processing systems, 30, 2017
2017
-
[18]
ctapipe-prototype open event reconstruc- tion pipeline for the cherenkov telescope array
Maximilian Linhoff, Lukas Beiske, Noah Biederbeck, Stefan Fr¨ ose, Karl Kosack, Lukas Nickel, K Abe, Jann Aschersleben, Andres Delgado Giler, CTA Consor- tium, et al. ctapipe-prototype open event reconstruc- tion pipeline for the cherenkov telescope array. In38th International Cosmic Ray Conference (ICRC2023)- Gamma-ray Astronomy (GA), page 703. Sissa Med...
2024
-
[19]
New methods of atmo- spheric cherenkov imaging for gamma-ray astronomy
VP Fomin, AA Stepanian, RC Lamb, c DA Lewis, M Punch, and TC Weekes. New methods of atmo- spheric cherenkov imaging for gamma-ray astronomy. i. the false source method.Astroparticle Physics, 2(2):137–150, 1994
1994
-
[20]
A new analysis method for reconstructing the arrival direction of tev gamma rays using a single imaging atmospheric cherenkov telescope.Astroparti- cle Physics, 15(1):1–18, 2001
RW Lessard, JH Buckley, V Connaughton, and S Le Bohec. A new analysis method for reconstructing the arrival direction of tev gamma rays using a single imaging atmospheric cherenkov telescope.Astroparti- cle Physics, 15(1):1–18, 2001
2001
-
[21]
Observations of the crab nebula and pulsar with the large-sized telescope prototype of the cherenkov telescope array.The As- trophysical Journal, 956(2):80, 2023
H Abe, K Abe, S Abe, Arnau Aguasca-Cabot, I Agudo, N Alvarez Crespo, LA Antonelli, C Aramo, A Arbet-Engels, C Arcaro, et al. Observations of the crab nebula and pulsar with the large-sized telescope prototype of the cherenkov telescope array.The As- trophysical Journal, 956(2):80, 2023
2023
- [22]
-
[23]
The disp method for analysing large zenith angle gamma-ray data.arXiv preprint arXiv:1109.6044, 2011
GD S ¸ent¨ urk. The disp method for analysing large zenith angle gamma-ray data.arXiv preprint arXiv:1109.6044, 2011
-
[24]
Domingo-Santamar´ ıa, J
E. Domingo-Santamar´ ıa, J. Flix, V. Scalzotto, W. Wittek, and J. Rico. The DISP analysis method for point-like or extendedγsource searches/studies with the MAGIC telescope.Part 3: MAGIC Detector and Analysis Details, page 23, 2005
2005
-
[25]
A new analysis method for very high definition imaging atmospheric cherenkov telescopes as applied to the cat telescope
St´ ephane Le Bohec, B Degrange, M Punch, A Barrau, R Bazer-Bachi, H Cabot, LM Chounet, G Debiais, JP Dezalay, A Djannati-Ata¨ ı, et al. A new analysis method for very high definition imaging atmospheric cherenkov telescopes as applied to the cat telescope. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detecto...
1998
-
[26]
A high per- formance likelihood reconstruction ofγ-rays for imag- ing atmospheric cherenkov telescopes.Astroparticle Physics, 32(5):231–252, 2009
Mathieu De Naurois and Lo¨ ıc Rolland. A high per- formance likelihood reconstruction ofγ-rays for imag- ing atmospheric cherenkov telescopes.Astroparticle Physics, 32(5):231–252, 2009
2009
-
[27]
A monte carlo template based analysis for air-cherenkov arrays.Astroparticle Physics, 56:26–34, 2014
RD Parsons and JA Hinton. A monte carlo template based analysis for air-cherenkov arrays.Astroparticle Physics, 56:26–34, 2014
2014
-
[28]
Acceler- ation and transport of relativistic electrons in the jets of the microquasar ss 433.Science, 383(6681):402– 406, 2024
HESS Collaboration*†, F Aharonian, F Ait Benkhali, J Aschersleben, H Ashkar, M Backes, V Barbosa Mar- tins, R Batzofin, Y Becherini, D Berge, et al. Acceler- ation and transport of relativistic electrons in the jets of the microquasar ss 433.Science, 383(6681):402– 406, 2024
2024
-
[29]
A hybrid approach to event reconstruction for at- mospheric cherenkov telescopes combining machine learning and likelihood fitting.Astroparticle Physics, 163:103008, 2024
Georg Schwefer, Robert Parsons, and Jim Hinton. A hybrid approach to event reconstruction for at- mospheric cherenkov telescopes combining machine learning and likelihood fitting.Astroparticle Physics, 163:103008, 2024
2024
-
[30]
Likelihood-free mcmc with amortized approximate ra- tio estimators
Joeri Hermans, Volodimir Begy, and Gilles Louppe. Likelihood-free mcmc with amortized approximate ra- tio estimators. InInternational conference on ma- chine learning, pages 4239–4248. PMLR, 2020
2020
-
[31]
Ctlearn: Deep learning for gamma-ray astronomy.arXiv preprint arXiv:1912.09877, 2019
Daniel Nieto, A Brill, Q Feng, TB Humensky, B Kim, T Miener, R Mukherjee, and J Sevilla. Ctlearn: Deep learning for gamma-ray astronomy.arXiv preprint arXiv:1912.09877, 2019
-
[32]
Lact: Status and future plans.Nu- clear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 1086:171373, 2026
Shoushan Zhang. Lact: Status and future plans.Nu- clear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 1086:171373, 2026
2026
-
[33]
Masked autoencoders are scalable vision learners
Kaiming He, Xinlei Chen, Saining Xie, Yanghao Li, Piotr Doll´ ar, and Ross Girshick. Masked autoencoders are scalable vision learners. InProceedings of the IEEE/CVF conference on computer vision and pat- tern recognition, pages 16000–16009, 2022
2022
-
[34]
Fermi large area telescope observations of the crab pulsar and nebula.The Astrophysical Journal, 708(2):1254–1267, 2010
AA Abdo, Markus Ackermann, Marco Ajello, WB At- wood, Magnus Axelsson, Luca Baldini, Jean Bal- let, Guido Barbiellini, MG Baring, Denis Bastieri, et al. Fermi large area telescope observations of the crab pulsar and nebula.The Astrophysical Journal, 708(2):1254–1267, 2010
2010
-
[35]
Tarek Hassan, O Gueta, G Maier, Maximilian N¨ othe, M Peresano, and I Vovk. Performance of a proposed event-type based analysis for the cherenkov telescope array.arXiv preprint arXiv:2108.10135, 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.