Abstract:ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015--2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons ($N_\mu>4$) and in the zenith angle range $0^{\circ}<\theta<50^{\circ}$. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about $ 4 \times 10^{15}<E_\mathrm{prim}< 6 \times 10^{16}$~eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than $20\%$ and $30\%$, respectively. The rate of high muon multiplicity events ($N_\mu>100$) obtained with QGSJET and SIBYLL is compatible with the data, while EPOS-LHC produces a significantly lower rate ($55\%$ of the measured rate). For both QGSJET and SIBYLL, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy $E_\mathrm{prim} \sim 10^{17}$~eV of these events. This result places significant constraints on more exotic production mechanisms.

Muons in EASs with $E_0 = 10^{19}$ eV according to data of the Yakutsk Array

Muon puzzle in ultra-high energy EASs according to Yakutsk array and Auger experiment data

Zenith-Angular Characteristics of Particles in EASs with $E_0 \simeq 10^{18}$ eV According to the Yakutsk Array Data

Estimation of the composition of ultra-high energy cosmic rays using the muon correlation method based on Yakutsk EAS array data

The rise of muon content in extensive air showers after the 3 PeV knee of the primary cosmic ray spectrum according to data of the Tien Shan mountain installation

A Direct Comparison of Muon Measurements at the Yakutsk Array and the Pierre Auger Observatory

Study of muons from ultra-high energy cosmic ray air showers measured with the Telescope Array experiment

Muon deficit in simulations of air showers inferred from AGASA data

UHECR primary identification using the lateral profile of muons in EAS

Azimuthal fluctuations and number of muons at the ground in muon-depleted proton air showers at PeV energies

Explaining muon excess in cosmic rays using the gluon condensation model

Mass composition of cosmic rays above 0.1 EeV by the Yakutsk array data

The forward particle production in the energy range of 1 PeV as seen with the Tibet hybrid experiment

On the model uncertainties for the predicted muon content of extensive air showers

Precise determination of muon shower content from shower universality property

On the calibration of ultra-high energy EASs at the Yakutsk array and Telescope Array

Reconstructing the Energy Spectrum of Cosmic Muons Measured with the Underground Detector of the Tien Shan High-Altitude Station

Air showers, hadronic models, and muon production

Multimuons in cosmic-ray events as seen in ALICE at the LHC