• ### Simulation study of the correlation ($X_{max}^{\mu}$, $N^{\mu}$) in view of obtaining information on primary mass of the UHECRs(1607.03244)

July 12, 2016 hep-ex, astro-ph.IM, astro-ph.HE
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 $X_{max}^{\mu}$ of the muon production rate longitudinal profile of Extensive Air Showers (EAS) and the number of muons, $N^{\mu}$, which hit an array of detectors located in the horizontal plane. We investigate the sensitivity of the 2D distribution $X_{max}^{\mu}$ versus $N^{\mu}$ 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 $10^{19}$eV and $10^{20}$eV, at five angles of incidence, $0^{\circ}$, $37^{\circ}$, $48^{\circ}$, $55^{\circ}$ and $60^{\circ}$. Using the simulations we obtain the 2D Probability Functions $Prob(X_{max}^{\mu},N^{\mu} \ | \ p)$ and $Prob(X_{max}^{\mu},N^{\mu} \ | \ Fe)$ which give the probability that a shower induced by a proton or iron nucleus contributes to a specific point on the plane ($X_{max}^{\mu}$, $N^{\mu}$). Then we construct the probability functions $Prob(p\ | \ X_{max}^{\mu},N^{\mu})$ and $Prob(Fe \ | \ X_{max}^{\mu},N^{\mu})$ which give the probability that a certain point on the plane ($X_{max}^{\mu}$, $N^{\mu}$) 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 $X_{max}^{\mu}$ distributions.
• We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
• Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
• A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
• ### $X_{max}^{\mu}$ vs. $N^{\mu}$ from Extensive Air Showers as estimator for the mass of primary UHECR's. Application for the Pierre Auger Observatory(1409.2972)
We study the possibility of primary mass estimation for Ultra High Energy Cosmic Rays (UHECR's) using the $X_{max}^{\mu}$ (the height where the number of muons produced on the core of Extensive Air Showers (EAS) is maximum) and the number $N^{\mu}$ of muons detected on ground. We use the 2D distribution - $X_{max}^{\mu}$ against $N^{\mu}$ in order to find its sensitivity to the mass of the primary particle. For that, we construct a 2D Probability Function $Prob(p,Fe \ | \ X_{max}^{\mu},N^{\mu})$ which estimates the probability that a certain point from the plane $(X_{max}^{\mu}$, $N^{\mu})$ corresponds to a shower induced by a proton, respectively an iron nucleus. To test the procedure, we analyze a set of simulated EAS induced by protons and iron nuclei at energies of $10^{19} eV$ and $20^{\circ}$ zenith angle with CORSIKA. Using the Bayesian approach and taking into account the geometry of the infill detectors from the Pierre Auger Observatory, we observe an improvement in the accuracy of the primary mass reconstruction in comparison with the results obtained using only the $X_{max}^{\mu}$ distributions.