Plastic scintillators are widely used as particle detectors in many fields, mainly, medicine, particle physics and astrophysics. Traditionally, they are coupled to a photo-multplier (PMT) but now silicon photo-multipliers (SiPM) are evolving as a promising robust alternative, specially in space born experiments since plastic scintillators may be a light option for low Earth orbit missions. Therefore it is timely to make a new analysis of the optimal design for experiments based on plastic scintillators in realistic conditions in such a configuration. We analyze here their response to an isotropic flux of electron and proton primaries in the energy range from 1 MeV to 1 GeV, a typical scenario for cosmic ray or space weather experiments, through detailed GEANT4 simulations. First, we focus on the effect of increasing the ratio between the plastic volume and the area of the photo-detector itself and, second, on the benefits of using a reflective coating around the plastic, the most common technique to increase light collection efficiency. In order to achieve a general approach, it is necessary to consider several detector setups. Therefore, we have performed a full set of simulations using the highly tested GEANT4 simulation tool: several parameters have been analyzed such as the energy lost in the coating, the deposited energy in the scintillator, the optical absorption, the fraction of scintillation photons that are not detected, the light collection at the photo-detector, the pulse shape and its time parameters and finally, other design parameters as the surface roughness, the coating reflectivity and the case of a scintillator with two decay components. This work could serve as a guide on the design of future experiments based on the use of plastic scintillators.

The current methods to determine the primary energy of ultra-high energy cosmic rays (UHECRs) are different when dealing with hadron or photon primaries. The current experiments combine two different techniques, an array of surface detectors and fluorescence telescopes. The latter allow an almost calorimetric measurement of the primary energy. Thus, hadron-initiated showers detected by both type of detectors are used to calibrate the energy estimator from the surface array (usually the interpolated signal at a certain distance from the shower core S(r0)) with the primary energy. On the other hand, this calibration is not feasible when searching for photon primaries since no high energy photon has been unambiguously detected so far. Therefore, pure Monte Carlo parametrizations are used instead. In this work, we present a new method to determine the primary energy of hadron-induced showers in a hybrid experiment based on a technique previously developed for photon primaries. It consists on a set of calibration curves that relate the surface energy estimator, S(r0), and the depth of maximum development of the shower, Xmax, obtained from the fluorescence telescopes. Then, the primary energy can be determined from pure surface information since S(r0) and the zenith angle of the incoming shower are only needed. Considering a mixed sample of ultra-high energy proton and iron primaries and taking into account the reconstruction uncertainties and shower to shower fluctuations, we demonstrate that the primary energy may be determined with a systematic uncertainty below 1% and resolution around 16% in the energy range from 10^{18.5} to 10^{19.6} eV. Several array geometries, the shape of the energy error distributions and the uncertainties due to the unknown composition of the primary flux have been analyzed as well.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

The search for photons at EeV energies and beyond has considerable astrophysical interest and will remain one of the key challenges for ultra-high energy cosmic ray (UHECR) observatories in the near future. Several upper limits to the photon flux have been established since no photon has been unambiguously observed up to now. An improvement in the reconstruction efficiency of the photon showers and/or better discrimination tools are needed to improve these limits apart from an increase in statistics. Following this direction, we analyze in this work the ability of the surface parameter Sb, originally proposed for hadron discrimination, for photon search. Semi-analytical and numerical studies are performed in order to optimize Sb for the discrimination of photons from a proton background in the energy range from 10^18.5 to 10^19.6 eV. Although not shown explicitly, the same analysis has been performed for Fe nuclei and the corresponding results are discussed when appropriate. The effects of different array geometries and the underestimation of the muon component in the shower simulations are analyzed, as well as the Sb dependence on primary energy and zenith angle.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

The current methods to determine the primary energy in surface arrays are different when dealing with hadron or photon initiated showers. In this work, we adapt a method previously developed for photon-initiated showers to hadron primaries. We determine the Monte Carlo parametrizations that relate the surface energy estimator and the maximum of shower development for both, proton and Iron primaries. Using for each primary their own set of calibration curves, which is of course impossible in practice, we show that the energy could be inferred with a negligible bias and 12% resolution. However, we show that a mixed calibration could also be performed, including both type of primaries, such that the bias still remains low and the achieved resolution is around 15%. In addition, the method allows the simultaneous determination of Xmax in pure surface arrays with resolution better than 7%.