Tuning Monte Carlo Models to Reproduce Cosmic Radiation Interacting with the Earth’s Atmosphere

In this work, we performed a comparative study between HIJING, Sibyll, and QGSJET model-based event generators. Such Monte Carlo (MC) models are used to simulate the interaction and propagation of high-energy cosmic radiation (e.g., coming from the sun) with the Earth’s atmosphere. The global event observables selected for the study were the transverse momentum ((Formula presented.)) spectra and rapidity density distributions of strange particles ((Formula presented.), (Formula presented.), and (Formula presented.)). This study was performed in the STAR and CMS fiducial phase spaces by simulating the strange particles in (Formula presented.) collisions at (Formula presented.) = 200 GeV, 900 GeV, and 7 TeV, and the simulations were then compared to the experimental measurements. It was observed that none of the discussed model-based event generators ultimately predicted the experimental results, except QGSJET, which generally agrees reasonably with the data. However, QGSJET does not produce (Formula presented.) particles; therefore, it does not provide any predictions for (Formula presented.). The other two models reproduced the data only in a limited rapidity or transverse momentum region while mainly underpredicting the data in the rest of the areas. These cosmic radiation simulation models are capable of covering the mid-rapidity regions of density distributions. Utilizing model-based observations, some fundamental parameters can be re-tuned and extrapolations to the highest energies can be investigated. Furthermore, these observations can provide valuable insights that could potentially constrain and improve perturbative- and non-perturbative-based QCD event generators, thereby facilitating a better understanding of the underlying physics.