Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present a new approach for the identification of ultra-high energy cosmic rays from sources using dynamic graph convolutional neural networks. These networks are designed to handle sparsely arranged objects and to exploit their short- and long-range correlations. Our method searches for patterns in the arrival directions of cosmic rays, which are expected to result from coherent deflections in cosmic magnetic fields. The network discriminates astrophysical scenarios with source signatures from those with only isotropically distributed cosmic rays and allows for the identification of cosmic rays that belong to a deflection pattern. We use simulated astrophysical scenarios where the source density is the only free parameter to show how density limits can be derived. We apply this method to a public data set from the AGASA Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

The arrival directions of ultra-high-energy cosmic rays appear to be approximately isotropically distributed over the whole sky, but the last-generation UHECR detector arrays, the Pierre Auger Observatory (Auger) and the Telescope Array (TA), have detected thousands of events, allowing us to study small deviations from isotropy. So far, Auger has detected a large-scale anisotropy consistent with a $\sim 6.5\%$ dipole moment in the distribution of cosmic rays with $E > 8$~EeV, and both collaborations have reported indications for smaller-scale anisotropies at higher energies. On the other hand, neither array has full-sky coverage, the Auger field of view being limited to declinations $\delta < +45^\circ$ and the TA one to $\delta > -16^\circ$. Searches for anisotropies with full-sky coverage thus require combining data from both arrays. A working group with members from both collaborations has been established for this task. Since even a minor systematic error in the energy determinations at either of the arrays could result in a sizeable spurious anisotropy along the north--south axis, we devised a method to cross-calibrate the energy scales of the two experiments with respect to each other by using events in a declination band within the intersection of their fields of view. In this contribution, we report on both updates on blind searches for anisotropies and the first full-sky studies of flux models based on possible local extragalactic sources.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

Despite deflections by Galactic and extragalactic magnetic fields, the distribution of ultra-high energy cosmic rays (UHECRs) over the celestial sphere remains a most promising observable for the identification of their sources. Thanks to a large number of detected events over the past years, a large-scale anisotropy at energies above 8 EeV has been identified, and there are also indications from the Telescope Array and Pierre Auger Collaborations of deviations from isotropy at intermediate angular scales (about 20 degrees) at the highest energies. In this contribution, we map the flux of UHECRs over the full sky at energies beyond each of two major features in the UHECR spectrum -- the ankle and the flux suppression --, and we derive limits for anisotropy on different angular scales in the two energy regimes. In particular, full-sky coverage enables constraints on low-order multipole moments without assumptions about the strength of higher-order multipoles. Following previous efforts from the two Collaborations, we build full-sky maps accounting for the relative exposure of the arrays and differences in the energy normalizations. The procedure relies on cross-calibrating the UHECR fluxes reconstructed in the declination band around the celestial equator covered by both observatories. We present full-sky maps at energies above 10 EeV and 50 EeV, using the largest datasets shared across UHECR collaborations to date. We report on anisotropy searches exploiting full-sky coverage and discuss possible constraints on the distribution of UHECR sources.