Simulation study of the correlation (Xmaxμ,Nμ) 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 Xmaxμ of the muon production rate longitudinal profile of Extensive Air Showers (EAS) and the number of muons, Nμ  , which hit an array of detectors located in the horizontal plane. We investigate the sensitivity of the 2D distribution Xmaxμ versus Nμ 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 1019 eV and 1020 eV, at five angles of incidence, 0°, 37°, 48°, 55° and 60°. Using the simulations we obtain the 2D Probability Functions Prob(Xmaxμ,Nμ|p) and Prob(Xmaxμ,Nμ|Fe) which give the probability that a shower induced by a proton or iron nucleus contributes to a specific point on the plane (Xmaxμ,Nμ  ). Then we construct the probability functions Prob(p|Xmaxμ,Nμ) and Prob(Fe|Xmaxμ,Nμ) which give the probability that a certain point on the plane (Xmaxμ,Nμ  ) 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 Xmaxμ distributions.