Abstract
In this paper we study the depth at which a cosmic ray shower reaches its maximum (Xmax) as predicted by Monte Carlo simulation. The use of Xmax in the determination of the primary particle mass can only be done by comparing the measured values with simulation predictions. For this reason it is important to study the differences between the available simulation models. We have done a study of the first and second moments of the Xmax distribution using the Corsika and Conex programs. The study was done with high statistics in the energy range from 1017 to 1020.4eV. We focus our analysis in the different implementations of the hadronic interaction models Sibyll2.1 and QGSJetII in Corsika and Conex. We show that the predictions of the 〈Xmax〉 and RMS(Xmax) depend slightly on the combination of simulation program and hadronic interaction model. Although these differences are small, they are not negligible in some cases (up to 5g/cm2 for the worse case) and they should be considered as a systematic uncertainty of the model predictions for 〈Xmax〉 and RMS(Xmax). We have included a table with the suggested systematic uncertainties for the model predictions. Finally, we present a parametrization of the Xmax distribution as a function of mass and energy according to the models Sibyll2.1 and QGSJetII, and showed an example of its application to obtain the predicted Xmax distributions from cosmic ray propagation models.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.