We compare the predictions of aafrag for the spectra of secondary photons, neutrinos, electrons, and positrons produced in proton-proton collisions to those of the parametrizations of Kamae et al., Kelner et al. and Kafexhiu et al. We find that the differences in the normalization of the photon energy spectra reach 20%–50% at intermediate values of the transferred energy fraction $x$, growing up to a factor of two for $x\to 1$, while the differences in the neutrino spectra are even larger. We argue that LHCf results on the forward production of photons and neutral pions favor the use of the qgsjet-II-04m model on which aafrag is based. The differences in the normalization have important implications in the context of multimessenger astronomy, in particular, for the prediction of neutrino fluxes, based on gamma-ray flux measurements, or regarding the inference of the cosmic ray spectrum, based on gamma-ray data. We note also that the positron-electron ratio from hadronic interactions increases with energy toward the cutoff, an effect which is missed using the average electron-positron spectrum from Kelner et al. Finally, we describe the publicly available python package aafragpy, which provides the secondary spectra of photons, neutrinos, electrons, and positrons. This package complements the aafrag results for protons with energies above 4 GeV with previous analytical parameterizations of particle spectra for lower energy protons.

• Received 4 October 2021
• Accepted 24 November 2021

DOI:https://doi.org/10.1103/PhysRevD.104.123027