Simulation study of the correlation (, ) in view of obtaining information on primary mass of the UHECRs
In this paper we study, using Monte Carlo simulations, the possibility to discriminate the mass of the Ultra High Energy Cosmic Rays (UHECRs) by combining information obtained from the maximum of the muon production rate longitudinal profile of Extensive Air Showers (EAS) and the number of muons, , which hit an array of detectors located in the horizontal plane. We investigate the sensitivity of the 2D distribution versus to the mass of the primary particle generating the air shower. To this purpose we analyze a set of CORSIKA showers induced by protons and iron nuclei at energies of eV and eV, at five angles of incidence, , , , and . Using the simulations we obtain the 2D Probability Functions and which give the probability that a shower induced by a proton or iron nucleus contributes to a specific point on the plane (, ). Then we construct the probability functions and which give the probability that a certain point on the plane (, ) corresponds to a shower initiated by a proton or an iron nucleus, respectively. Finally, a test of this procedure using a Bayesian approach, confirms an improved accuracy of the primary mass estimation in comparison with the results obtained using only the distributions.