Testing Model Predictions of Depth of Air-Shower Maximum and Signals in Surface Detectors using Hybrid Data of the Pierre Auger Observatory

Institute of Physics of the Czech Academy of Sciences

^* Corresponding author: vicha@fzu.cz
^** Corresponding author: spokespersons@auger.org

Published online: 28 April 2023

Abstract

We present a method for testing the predictions of hadronic interaction models and improving their consistency with observed two-dimensional distributions of the depth of shower maximum, X _max, and signal at the ground level as a function of zenith angle. The method relies on the assumption that the mass composition is the same at all zenith angles, while the atmospheric shower development and attenuation depend on composition in a correlated way. In the present work, for each of the three leading LHC-tuned hadronic interaction models, we allow a global shift ∆ X _max of the predicted shower maximum, which is the same for every mass and energy, and a rescaling R _Had of the hadronic component at the ground level which is constant with the zenith angle. We apply the analysis to 2297 events reconstructed with both the fluorescence and surface detectors of the Pierre Auger Observatory with energies 10^18.5−19.0 eV and zenith angles below 60°. Given the modeling assumptions made in this analysis, the best fit reaches its optimum value when shifting the X _max predictions of hadronic interaction models to deeper values and increasing the hadronic signal. This change in the predicted X _max scale alleviates the previously identified model deficit in the hadronic signal (commonly called the muon puzzle) but does not fully remove it. Because of the size of the adjustments ∆ X _max and R _Had and the large number of events in the sample, the statistical significance of need for these adjustments is large, greater than 5σ_stat, even for the combination of the systematic experimental shifts within 1σ_sys that are the most favorable for the models.