Speaker
Description
We present a new 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_\text{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 $\Delta X_\text{max}$ of the predicted shower maximum, which is the same for every mass and energy, and a rescaling $R_\text{Had}$ of the hadronic component at the ground level which depends on 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$^\circ$. Given the modeling assumptions made in this analysis, the best fit reaches its optimum value when shifting the $X_\text{max}$ predictions of hadronic interaction models to deeper values and increasing the hadronic signal at both extreme zenith angles. This change in the predicted $X_\text{max}$ scale alleviates the previously identified model deficit in the hadronic signal (commonly called the muon puzzle) but does not remove it. Because of the size of the adjustments $\Delta X_\text{max}$ and $R_\text{Had}$ and the large number of events in the sample, the statistical significance of these assumed adjustments is large, greater than 5$\sigma_\text{stat}$, even for the combination of the systematic experimental shifts within 1$\sigma_\text{sys}$ that are the most favorable for the models.
