6. Conclusion
In this paper, we have analyzed the muon content of air showers, proposed a parametrization of the first two moments of the number of muons with energy and primary particle mass and showed how the measured mean and σ of log10Nμ can be used to discriminate between composition scenarios. We proposed a model to describe log10 Nmeas μ and σ[log10 Nmeas μ ] as a function of energy and primary particle mass (A). This model was conceived to keep the most relevant hadronic interaction uncertainties concentrated in only two parameters (a and b). We have validated the model with Monte Carlo simulation of the air shower and its capability to describe the log10 Nmeas μ moments was proven. Six composition scenarios were considered. The particle flux predicted by these scenarios was transformed into the corresponding log10 Nmeas μ and σ[log10 Nmeas μ ] evolution with energy. The log10 Nmeas μ and σ[log10 Nmeas μ ] evolution with energy was fitted using the proposed parametrization. A comparison of the log10 Nmeas μ and σ[log10 Nmeas μ ] model using a simple χ2 test allows the discrimination between the scenarios. The discrimination is effective even considering the systematic uncertainties on the Nμ prediction and on energy scale uncertainty