Testing model energy spectra of charged particles produced in hadron interactions on the basis of atmospheric muons

Abstract

An original method for calculating the spectrum of atmospheric muons with the aid of the CORSIKA 7.4 code package and numerical integration is proposed. The first step consists in calculating the energy distribution of muons for various fixed energies of primary-cosmic-ray particles and within several chosen hadron-interaction models included in the CORSIKA 7.4 code package. After that, the spectrum of atmospheric muons is calculated via integrating the resulting distribution densities with the chosen spectrum of primary-cosmic-ray particles. The atmospheric-muon fluxes that were calculated on the basis of the SIBYLL 2.1, QGSJET01, and QGSJET II-04 models exceed the predictions of the wellknown Gaisser approximation of this spectrum by a factor of 1.5 to 1.8 in the range of muon energies between about 103 and 104 GeV.Under the assumption that, in the region of extremely highmuon energies, a dominant contribution to the muon flux comes from one to two generations of charged π± and K ± mesons, the production rate calculated for these mesons is overestimated by a factor of 1.3 to 1.5. This conclusion is confirmed by the results of the LHCf and TOTEM experiments.