Time-Structured Tail Probabilities for Ultra-High-Energy Gamma-Hadron Discrimination in Water-Cherenkov Arrays

TL;DR

The paper addresses the challenge of discriminating ultra-high-energy gamma rays from hadronic backgrounds using sparse water-Cherenkov detector arrays. It introduces the time-structured tail observable $P^{\alpha, T}_{\rm tail}$, extending traditional tail-probability methods to include radial and time-bin information via cumulative distributions. In SD-433-like simulations around $E \sim 10^{17}$ eV, this observable achieves a background contamination of about $2\times 10^{-2}$ at $50\%$ gamma efficiency, outperforming existing WCD-only discriminators by roughly a factor of five and approaching an idealized muon-isolating benchmark. The study demonstrates that incorporating time-domain information substantially enhances photon searches in large-scale WCD arrays and outlines practical calibration pathways and prospects for future, denser detectors to further close the gap to muon-based discrimination.

Abstract

Gamma-hadron discrimination based on shower observables is essential for identifying gamma-ray astrophysical sources at the highest energies. In this work, we introduce $P^{α, T}_{\rm tail}$, a new discrimination variable for ultra-high-energy photon searches within the framework of a water-Cherenkov detector (WCD) array. The observable extends signal-integrated methods by incorporating the time structure of WCD traces, using cumulative signal distributions. Using simulated proton- and gamma-induced air showers at energies around $10^{17}\,\mathrm{eV}$, we evaluate the performance of $P^{α, T}_{\rm tail}$ and compare it with established WCD-based observables such as $S_b$, risetime-based variables, and the SWGO-inspired, $P^α_{\rm tail}$. The new variable attains a background contamination of roughly $2 \times 10^{-2}$ at $50\%$ gamma efficiency, improving upon existing WCD-only methods by nearly a factor of five and approaching the performance of an idealized muon-isolating reference. These results demonstrate the effectiveness of exploiting time-resolved signal tails to enhance ultra-high-energy photon searches in sparse surface arrays.

Figures (10)

• Figure 1: Reconstructed energy distributions of the set of simulated proton (red) and gamma-ray (blue) events.
• Figure 2: Reconstructed zenith distributions of the set of simulated proton (red) and gamma-ray (blue) events.
• Figure 3: Diagram illustrating the differences in development and time structure of extensive air showers at different energies. The left hand side roughly represents the typical longitudinal profile of the showers, while the right hand presents the characteristics of the showers at $X_{\rm max}$ (dotted lines), and at the ground level (dashed lines). The difference in arrival time between two possible particle trajectories, $\Delta t \equiv t_2 - t_1$, is of the order of hundreds of nanoseconds for the lower energy shower ($100\,$TeV, blue) and of the order of a microsecond for the higher energy shower ($100\,$PeV, green). The number of particles in the longitudinal profile of the lower energy shower has been scaled up by a factor of 100 for visualization purposes.
• Figure 4: Cumulative signal distributions of the time trace for the ring spanning from $r_i=300$ m to $r_i=350$ m. The time of each bin is expressed in relation to the arrival time of the shower core. The two areas highlighted in the top figure correspond to the distributions shown in the bottom figure.
• Figure 5: Cumulative signal distributions of the trace time bins for the ring spanning from $r_i=500$ m to $r_i=550$ m. The time of each bin is expressed in relation to the arrival time of the shower core. The two areas highlighted in the top figure correspond to the distributions shown in the bottom figure.