Pith. sign in

REVIEW 3 major objections 6 minor 32 references

Reconstructing Air-Shower Observables using a Universality-Based Model at the Pierre Auger Observatory

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A universality-based shower model reconstructs the depth of the shower maximum and the relative muon number from Pierre Auger surface-detector time traces, and the shower-depth estimates agree with fluorescence measurements.

desk verdict Preliminary but honest application of the universality model to Auger Phase-I data: event-level Xmax correlation is real, but the mean scale is FD-calibrated and the constant-bias assumption is the main soft spot. read the letter →

arxiv 2507.13209 v1 pith:VGICUYVJ submitted 2025-07-17 astro-ph.HE hep-ph

Maximilian Stadelmaier (for the Pierre Auger Collaboration) This is my paper
classification astro-ph.HEhep-ph
keywords air-showeruniversalityultra-high-energycosmicraysPierreAugerObservatoryshowermaximumdepthsurfacedetectortimetracesmuonnumbermasscompositionelongationrate
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that air-shower universality, the principle that showers sharing the same energy, arrival direction, depth of maximum, and muon content have identical particle distributions, is strong enough to serve as a reconstruction tool rather than only a simulation concept. The authors fit a universality-based model directly to the 25-ns time traces measured by the water-Cherenkov stations of the Pierre Auger surface detector and show that the depth of the shower maximum $X_{\rm max}$ can be inferred from the timing of the arriving particles, in particular from the time quantile $t_{40}$ at which 40% of a station's signal is deposited. When an independent fluorescence energy estimate is available, the same model also returns the relative muon number $R_\mu$ per event. Since the surface detector operates around the clock while the fluorescence detector is limited to clear moonless nights, the method would extend mass-composition measurements to many more events and to the highest energies, which matters because the cosmic-ray flux is too low for fluorescence-only statistics there.

What carries the argument

The central object is the universality shower model, in which the shower's particle content is split into components whose longitudinal, lateral, and time-at-ground distributions are parametrized as functions of $E_0$, $X_{\rm max}$, $R_\mu$, and the event geometry. For the $X_{\rm max}$ reconstruction, the decisive mechanism is the relation between the time quantile $t_{40}$, the time at which 40% of the total signal of a station has arrived, and the depth difference between that station and the shower maximum, obtained with a quasi-spherical description of the shower front. This mechanism converts the shape of each 25-ns water-Cherenkov time trace into a geometric constraint on the shower development, and it is what allows the fit to work without fluorescence information; normalizing the traces to probability densities removes most of the sensitivity to energy and muon number.

What would settle it

Compare universality-model $X_{\rm max}$ with fluorescence-measured $X_{\rm max}$ for the same Golden Hybrid events in narrow energy and zenith-angle bins; if the mean residual is not constant across bins, or if the residual width grows with energy, the constant-offset calibration fails and the reported elongation rate and $\sigma(X_{\rm max})$ would be biased.

Watch

Extended reading notes

Core claim

On its own terms, the paper's central claim is that a parametric model built on air-shower universality, describing every shower through its primary energy $E_0$, zenith angle, depth of maximum $X_{\rm max}$, and relative muon number $R_\mu$, can be fitted to the surface-detector time traces to recover these observables at the event level. For $X_{\rm max}$, the load-bearing step is the modeled relation between the time quantile $t_{40}$ of each normalized trace and the distance of the detector from the shower maximum, expressed through a quasi-spherical shower model; the paper validates this relation against Epos-LHC simulations and against Golden Hybrid events where the same showers are also seen by the fluorescence detector, finding a significant correlation between the reconstructed and directly measured $X_{\rm max}$. The paper then applies the model to Phase-I surface-detector data above 4 EeV and reports a reconstructed mean $X_{\rm max}$ that rises with energy, an elongation rate that is not constant but is best described by a broken line, and fluctuations $\sigma(X_{\rm max})$ that agree with fluorescence results above about 10 EeV; at lower energies the reconstructed fluctuations are larger than expected, which the authors attribute to performance differences on data or to outliers. A constant-offset calibration against fluorescence is used to correct the overall mean, and the reported systematic uncertainties are dominated by the energy-scale and $X_{\rm max}$-calibration uncertainties inherited from the fluorescence detector.

Load-bearing premise

The argument assumes that the simulation-derived mapping from the 40% arrival-time quantile of each detector trace to the depth of the shower maximum, built with the quasi-spherical shower model and Epos-LHC, is accurate for real air showers, and that any bias is a single constant offset that an energy-independent calibration can remove.

Editorial extensions

If this is right

  • Mass-composition studies can be pushed to energies and event numbers where fluorescence data alone are statistically too thin, because the surface detector runs continuously.
  • The same event-level fit yields $X_{\rm max}$ and $R_\mu$ inside one physically motivated framework, so correlations between the two mass-sensitive observables can be studied event by event.
  • The reconstructed $X_{\rm max}$ can separate lighter and heavier primary events, which the paper notes can reduce background in arrival-direction analyses.
  • Combining the muon-number estimate with an independent energy estimator gives a handle on the known muon deficit between air-shower data and hadronic-interaction predictions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct test of the calibration assumption would bin the residuals between universality-model and fluorescence $X_{\rm max}$ by energy and zenith angle; if the mean residual drifts across bins, the constant-offset calibration would distort the reported elongation rate and $\sigma(X_{\rm max})$.
  • Because the $t_{40}$--$X_{\rm max}$ mapping is derived from Epos-LHC simulations, regenerating the model with a different hadronic interaction model, such as QGSJetII-04, would show how much of the event-level scatter is physics-model dependent.
  • The same universality machinery could reconstruct other profile-related observables, such as the muon production depth, from the same time traces without new detector hardware.
  • The larger-than-expected $\sigma(X_{\rm max})$ below 10 EeV could be investigated by comparing per-event reconstructions against independent machine-learning methods on the same data set, isolating whether the excess comes from the model or from outlier events.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript presents a universality-based model of air-shower particle densities and uses it to reconstruct observables from data of the Pierre Auger Observatory. The model is applied in two modes: a Golden Hybrid analysis in which the FD energy is used as input to estimate the relative muon number, and an SD-only analysis in which the time structure of the detector traces, specifically the t40 time quantile, is used to estimate Xmax. The method is validated on simulations and on Golden Hybrid events, and preliminary Phase-I SD results are shown for the energy evolution of the mean Xmax and its fluctuations. The mean Xmax is calibrated to FD measurements with an energy-independent offset, and the paper reports a broken-line elongation rate and a low-energy excess in sigma(Xmax).

Significance. If the reconstruction method is validated, it would be valuable because the SD array operates with a much higher duty cycle than the FD, allowing Xmax measurements at the highest energies where FD statistics are limited. The paper's strengths are its physically motivated model, the event-level validation against Golden Hybrid data shown in Fig. 2, and the explicit, honest discussion of systematic limitations. The muon-number result in Fig. 1 also provides a useful cross-check of the known muon deficit. However, the central physical conclusions about the mean Xmax evolution and its fluctuations depend on an energy-independent calibration assumption and on an unexplained low-energy sigma discrepancy; these issues are load-bearing and need to be resolved before the results can be considered quantitatively reliable.

major comments (3)
  1. [Section 7 and Fig. 3] The absolute scale of the reported mean Xmax is calibrated to FD measurements with a single energy-independent offset, and the same FD data are then used as the comparison benchmark in Fig. 3. This makes the mean-level agreement partly circular and does not validate the energy dependence of any residual reconstruction bias. The first three energy bins in Fig. 3 were omitted from the calibration fit and have no assigned systematics, but the text does not explain how the constant-offset assumption is justified for those bins. Please show the calibration residuals as a function of energy and zenith angle, and provide evidence from simulations (for example with a different hadronic interaction model) that the bias is indeed constant.
  2. [Section 6, sigma(Xmax) at low energies] The paper reports that below about 10 EeV the reconstructed sigma(Xmax) is consistently larger than expected and attributes this to possible discrepancies between data and simulation performance or to over-represented outliers. This unexplained excess appears in exactly the energy range where the subtraction of the expected intrinsic resolution is most critical, and the sigma result is presented as a physics output with composition implications. The current treatment is not sufficient to rule out a reconstruction bias as the cause. Please quantify the resolution estimate used for the subtraction, report closure tests in simulations as a function of energy, and examine the sensitivity of the low-energy sigma to outlier removal and to the calibration assumption.
  3. [Section 6, t40-to-Xmax mapping and broken-line fit] The mapping from the t40 time quantile to Xmax is derived from Epos-LHC simulations through the quasi-spherical shower model in Ref. [20] and is corrected only by an overall constant offset. Any energy-, zenith-, or mass-dependent bias in this mapping would distort the reported elongation rate, the broken-line behavior, and sigma(Xmax). The one-break fit is quoted as chi2/ndf = 45.2/11 with p = 3e-6, but no fit parameters or systematic uncertainties are given, and this poor chi2 could equally be a signature of a non-constant reconstruction bias. Please provide a cross-check with another hadronic model (e.g., QGSJetII-04) and show residuals of the Golden Hybrid validation as a function of energy and zenith angle.
minor comments (6)
  1. [Section 2] The phrase 'approximately 109 particles per EeV' should render the exponent correctly as 10^9 particles per EeV.
  2. [Section 6] The definition of t40 and the trace-normalization procedure are only referenced to Ref. [20]; a brief explicit definition in the text would make the paper easier to evaluate on its own.
  3. [Fig. 2] The Pearson correlation coefficients for both simulations and Golden Hybrid data appear only in the figure legend; they should be quoted in the text as well.
  4. [Fig. 3] The caption does not describe what the black caps ('systematic uncertainties') include; please specify the components and how they are propagated.
  5. [Section 6] The phrase 'removing the expected intrinsic precision, as estimated using Monte-Carlo simulations' should state explicitly how this resolution is obtained and whether it is the same resolution shown in the simulation panel of Fig. 2.
  6. [References] Reference [24] appears to duplicate Reference [16] (both are J. Matthews, Astropart. Phys. 22 (2005) 387); please consolidate the duplicate citation.

Circularity Check

1 steps flagged · score 6.0 of 10

Partial circularity: the absolute mean Xmax is calibrated to FD data and then reported as agreeing with FD, while the event-level correlation and the shape of the Xmax evolution remain independent.

  1. fitted input called prediction [Section 7 (Systematic Uncertainties and Calibration) and Section 6 (Surface Detector Data: Depth of the Shower Maximum), Figs. 2-3]
    "The latter is necessary to correct for an overall difference of the mean Xmax provided by raw results of the universality fit and direct measurements; a similar calibration was performed in Refs. [7, 19, 31]. ... Where they are available, the results of the universality-model fit (after calibration of the mean) are in reasonable agreement with the data from direct FD measurements in terms of the evolution with primary energy."

    The calibration step fits a constant offset in the mean Xmax to the FD measurements. Since the same FD mean is then used as the reference, the reported agreement of the calibrated mean with FD is enforced by that fit, not demonstrated by the data. The offset is energy-independent, so the elongation-rate shape and sigma(Xmax) are not fixed by the calibration and retain independent content; the event-level correlation in Fig. 2 also remains an independent check. But the absolute Xmax scale and any statement about agreement of the mean are inputs, not predictions.

full rationale

The main claim, that the universality model can estimate Xmax from SD time traces at the event level, is supported by the event-by-event correlation with Golden Hybrid FD data (Fig. 2, right), which is not calibrated and is genuinely informative. The muon-number analysis explicitly uses FD energy as input and is therefore an external-reference measurement, not circular. The one concrete reduction-by-construction is the absolute normalization of the mean Xmax: Section 7 calibrates the overall mean difference between the universality fit and FD direct measurements, and Section 6 then reports that the calibrated mean Xmax agrees with FD. The constant offset cannot affect the elongation-rate slope or the width sigma(Xmax), so those results remain non-circular. The t40-to-Xmax relation is imported from the same authors' Ref. [20], but because it is tested against independent FD data and is not itself the target claim, I do not count it as a separate circular step. The paper is transparent about the calibration and about the unexplained low-energy sigma excess, which keeps this at partial rather than full circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central measurement rests on a simulation-based universality model (from Ref. [20]), an assumed constant calibration offset to FD data, and the reliability of the Epos-LHC interaction model. No new physical entities are introduced. The free parameters are the calibration offset, the inherited model parameters, and the MC-derived resolution correction.

free parameters (3)
  • Mean Xmax calibration offset = not provided
    Applied in Section 7 to bring the universality-fit mean Xmax into agreement with FD measurements; the offset is fitted to data and sets the absolute scale of Fig. 3.
  • Universality model parameters (lateral, longitudinal, temporal profiles per component) = not provided in this paper
    Parametrized in Ref. [20] using Corsika+Epos-LHC detector-response simulations (Section 4); inherited here without re-derivation.
  • Expected intrinsic Xmax resolution (used for sigma estimate) = from Monte Carlo simulations
    Used in Section 6 to subtract the reconstruction resolution from the observed variance; the MC resolution values depend on the simulation model and energy.
assumptions (6)
  • domain assumption Air-shower universality: showers with same energy, zenith, Xmax, and muon number have the same electromagnetic particle distributions.
    Invoked in Section 4 as the basis of the model.
  • domain assumption Epos-LHC is an adequate hadronic-interaction model for parametrizing the universality model.
    Section 4: the model is parametrized with Epos-LHC showers.
  • domain assumption The time quantile t40 of detector traces is directly related to the depth difference to Xmax via a quasi-spherical shower model.
    Section 6, based on Ref. [20]; central to the Xmax fit.
  • ad hoc to paper The calibration of the mean Xmax to FD is a constant offset, independent of energy and zenith angle.
    Section 7 states a calibration of the mean is applied; the paper does not test whether the offset is energy-dependent, which would bias the elongation rate.
  • domain assumption The FD energy scale (with ~14% systematic uncertainty) is reliable as input for the muon fit.
    Sections 5 and 7 use EFD as input for the muon-number reconstruction and inherit its systematics.
  • domain assumption The decadal elongation rate D=56 g/cm2 used to compute X19max is correct.
    Section 6 uses this constant to remove the energy trend when comparing across energies.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Reconstructing Air-Shower Observables using a Universality-Based Model at the Pierre Auger Observatory." pith.science (2026). https://pith.science/paper/VGICUYVJ

@misc{pith2026250713209,
  author       = {Pith},
  title        = {Pith review of: Reconstructing Air-Shower Observables using a Universality-Based Model at the Pierre Auger Observatory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VGICUYVJ}},
  note         = {Machine review of arXiv:2507.13209}
}
read the original abstract

Based on solutions of the cascade equations, the air-shower universality is a framework that for all air showers with the same energy, zenith angle, depth of shower maximum, and muon number predicts the same longitudinal, lateral, and energy distributions of electromagnetic shower particles. We employ a universality-based model of shower development that incorporates hadronic particle components to reconstruct observables from extensive air showers produced by ultra-high-energy cosmic rays. The model can estimate key parameters, such as the depth of the shower maximum and the number of muons at the event level. We discuss the performance of the reconstruction algorithm using both air-shower simulations, and preliminary results obtained from the Phase-I data of the Pierre Auger Observatory.

Figures

Figures reproduced from arXiv: 2507.13209 by the authors.

Figure 1
Figure 1. Left: The relative muon number 𝑅μ as a function of the primary energy 𝐸 FD 0 measured by the fluorescence detector. The expectations from different hadronic interaction models and primary particles are given as blue and red reference lines. Expectations according to the 𝑋max measurements of the showers are shown in orange. The systematic uncertainties of the estimated number of muons are shown as black brackets arou… view at source ↗
Figure 2
Figure 2. Left: True and reconstructed values of 𝑋max for simulated showers above a primary energy of 3 EeV. The number of data points as well as overall mean and width of the residuals Δ𝑋 rec max are given alongside the Pearson correlation coefficient in the upper left corner. Dot colors correspond to proton, helium, oxygen, and iron. Right: Correlation of the estimated (univ) and directly measured (FD) values of the depth o… view at source ↗
Figure 3
Figure 3. Left: ⟨𝑋max⟩ as a function of the (SD reconstructed) energy with a broken-line fit. Right: Fluctuations of 𝑋max given as the estimated standard deviation of the 𝑋max-distribution, 𝜎(𝑋max), as a function of the primary energy. The black markers denote the mean (or standard-deviation) values from the universality fit, orange markers show the results from direct FD measurements. Black caps stand for the systematic unce… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

32 extracted references · 19 canonical work pages

  1. [20]

    A Model of the Response of Surface Detectors to Extensive Air Showers Based on Shower Universality

    M. Stadelmaier, R. Engel, M. Roth et al.,Phys. Rev. D110 (2024) 023030 [2405.03494]

    work page Pith review arXiv 2024

  2. [1]

    Telescope ArrayCollaboration, Nucl. Phys. B Proc. Suppl.175(2008) 221

    work page 2008

  3. [2]

    Pierre AugerCollaboration, Nucl. Instrum. Meth. A798(2015) 172 [1502.01323]

    arXiv 2015

  4. [3]

    Pierre AugerCollaboration, Science 357 (2017) 1266–1270

    work page 2017

  5. [4]

    Pierre AugerCollaboration, Phys. Rev. D90(2014) 122006 [1409.5083]

    arXiv 2014

  6. [5]

    Pierre AugerCollaboration, Phys. Lett. B762(2016) 288 [1609.08567]

    arXiv 2016

  7. [6]

    Pierre AugerCollaboration, Phys. Rev. Lett.125(2020) 121106 [2008.06488]

    work page Pith review arXiv 2020

  8. [7]

    Pierre AugerCollaboration, Phys. Rev. D96(2017) 122003 [1710.07249]

    arXiv 2017

Show all 32 references
  1. [8]

    Hillas,J

    A.M. Hillas,J. Phys. G8 (1982) 1461

    1982

  2. [9]

    Lipari,Phys

    P. Lipari,Phys. Rev. D79(2009) 063001 [0809.0190]

    2009 arXiv

  3. [10]

    Nerling, J

    F. Nerling, J. Blümer, R. Engel et al.,Astropart. Phys.24 (2006) 421

    2006

  4. [11]

    Pierre AugerCollaboration, Nucl. Instrum. Meth. A620 (2010) 227 [0907.4282]

    2010 arXiv

  5. [12]

    Pierre AugerCollaboration, Nucl. Instrum. Meth. A586 (2008) 409 [0712.2832]

    2008 arXiv

  6. [13]

    Pierre AugerCollaboration, Nucl. Instrum. Meth. A568 (2006) 839 [2102.01656]

    2006 arXiv

  7. [14]

    Pierre AugerCollaboration, JCAP 08(2014) 019 [1407.3214]

    2014 arXiv

  8. [15]

    Instrum.15(2020) 10021 [2007.09035]

    Pierre AugerCollaboration, J. Instrum.15(2020) 10021 [2007.09035]

    2020 arXiv

  9. [17]

    Pierre AugerCollaboration, Phys. Rev. D91(2015) 032003 [1408.1421]

    2015 arXiv

  10. [18]

    Pierre AugerCollaboration, PoS ICRC2023(2023) 339

    2023

  11. [19]

    Pierre AugerCollaboration, Phys. Rev. D111(2025) 022003 [2406.06319]

    2025 arXiv

  12. [21]

    Argiro, S.L.C

    S. Argiro, S.L.C. Barroso, J. Gonzalez et al.,Nucl. Instrum. Meth. A580 (2007) 1485 [0707.1652]

    2007 arXiv

  13. [22]

    D. Heck, J. Knapp, J. Capdevielle et al.,Report FZKA6019 (1998)

    1998

  14. [23]

    Pierog, I

    T. Pierog, I. Karpenko, J.M. Katzy et al.,Phys. Rev. C92(2015) 034906 [1306.0121]

    2015 arXiv

  15. [24]

    Matthews,Astropart

    J. Matthews,Astropart. Phys.22(2005) 387

    2005

  16. [25]

    Ostapchenko,Phys

    S. Ostapchenko,Phys. Rev. D83(2011) 014018 [1010.1869]

    2011 arXiv

  17. [26]

    Arteaga Velazquez,PoS ICRC2023(2023) 466

    J.C. Arteaga Velazquez,PoS ICRC2023(2023) 466

    2023

  18. [27]

    Pierre AugerCollaboration, Eur. Phys. J. Plus131 (2016) 301 [1609.02498]

    2016 arXiv

  19. [28]

    Pierre AugerCollaboration, PoS UHECR2024 (2025) 091

    2025

  20. [29]

    Pierre AugerCollaboration, PoS ICRC2023(2023) 371

    2023

  21. [30]

    Pierre AugerCollaboration, PRL (2024) [2406.06315]

    2024 arXiv

  22. [31]

    Pierre AugerCollaboration, JINST 16 (2021) P07019 [2101.02946]

    2021 arXiv

  23. [32]

    Vícha, A

    J. Vícha, A. Bakalová, A.L. Müller et al.,2504.11985

    arXiv

  24. [33]

    Tor Vergata

    Pierre AugerCollaboration, PoS ICRC2025(2025) 274. 8 Shower Universality The Pierre Auger Collaboration A. Abdul Halim 13, P. Abreu 70, M. Aglietta 53,51, I. Allekotte 1, K. Almeida Cheminant 78,77, A. Almela7,12, R. Aloisio 44,45, J. Alvarez-Muñiz 76, A. Ambrosone 44, J. Amme...

    2025

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.