Muon deficit in air shower simulations estimated from AGASA muon measurements
In this work, direct measurements of the muon density at $1000\,\textrm{m}$ from the shower axis obtained by the Akeno Giant Air Shower Array (AGASA) are analysed. The selected events have zenith angles $θ\leq 36^{\textrm{o}}$ and reconstructed energies in the range $18.83\,\leq\,\log_{10}(E_{R}/\textrm{eV})\,\leq\,19.46$. These are compared to the predictions corresponding to proton, iron, and mixed composition scenarios obtained by using the high-energy hadronic interaction models EPOS-LHC, QGSJetII-04, and Sibyll2.3c. The mass fractions of the mixed composition scenarios are taken from the fits to the depth of the shower maximum distributions performed by the Pierre Auger Collaboration. The cross-calibrated energy scale from the Spectrum Working Group [D. Ivanov, for the Pierre Auger Collaboration and the Telescope Array Collaboration, PoS(ICRC2017) 498 (2017)] is used to combine results from different experiments. The analysis shows that the AGASA data are compatible with a heavier composition with respect to the one predicted by the mixed composition scenarios. Interpreting this as a muon deficit in air shower simulations, the incompatibility is quantified. The muon density obtained from AGASA data is greater than that of the mixed composition scenarios by a factor of $1.49\pm0.11\,\textrm{(stat)}\pm0.18\,\textrm{(syst)}$, $1.54\pm0.12\,\textrm{(stat)}\pm0.18\,\textrm{(syst)}$, and $1.66\pm0.13\,\textrm{(stat)}\pm0.20\,\textrm{(syst)}$ for EPOS-LHC, Sibyll2.3c, and QGSJetII-04, respectively.