-
Ultrahigh energy cosmic ray air showers probe particle physics at energies
beyond the reach of accelerators. Here we introduce a new method to test
hadronic interaction models without relying on the absolute energy calibration,
and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose
longitudinal development and lateral distribution were simultaneously measured
by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16
(1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models
EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.
-
Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre
Auger Observatory to extend its range of detection and to directly measure the
muon content of the particle showers. It consists of an infill of surface
water-Cherenkov detectors accompanied by buried scintillator detectors used for
muon counting. The main objectives of the AMIGA engineering array, referred to
as the Unitary Cell, are to identify and resolve all engineering issues as well
as to understand the muon-number counting uncertainties related to the design
of the detector. The mechanical design, fabrication and deployment processes of
the muon counters of the Unitary Cell are described in this document. These
muon counters modules comprise sealed PVC casings containing plastic
scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier
tubes, and acquisition electronics. The modules are buried approximately 2.25 m
below ground level in order to minimize contamination from electromagnetic
shower particles. The mechanical setup, which allows access to the electronics
for maintenance, is also described in addition to tests of the modules'
response and integrity. The completed Unitary Cell has measured a number of air
showers of which a first analysis of a sample event is included here.
-
The Pierre Auger Observatory has begun a major Upgrade of its already
impressive capabilities, with an emphasis on improved mass composition
determination using the surface detectors of the Observatory. Known as
AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors
on top of all 1660 water-Cherenkov detectors, updated and more flexible surface
detector electronics, a large array of buried muon detectors, and an extended
duty cycle for operations of the fluorescence detectors. This Preliminary
Design Report was produced by the Collaboration in April 2015 as an internal
document and information for funding agencies. It outlines the scientific and
technical case for AugerPrime. We now release it to the public via the arXiv
server. We invite you to review the large number of fundamental results already
achieved by the Observatory and our plans for the future.
-
The azimuthal asymmetry in the risetime of signals in Auger surface detector
stations is a source of information on shower development. The azimuthal
asymmetry is due to a combination of the longitudinal evolution of the shower
and geometrical effects related to the angles of incidence of the particles
into the detectors. The magnitude of the effect depends upon the zenith angle
and state of development of the shower and thus provides a novel observable,
$(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays
above $3 \times 10^{18}$ eV. By comparing measurements with predictions from
shower simulations, we find for both of our adopted models of hadronic physics
(QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass
increases slowly with energy, as has been inferred from other studies. However,
the mass estimates are dependent on the shower model and on the range of
distance from the shower core selected. Thus the method has uncovered further
deficiencies in our understanding of shower modelling that must be resolved
before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.
-
To exploit the full potential of radio measurements of cosmic-ray air showers
at MHz frequencies, a detector timing synchronization within 1 ns is needed.
Large distributed radio detector arrays such as the Auger Engineering Radio
Array (AERA) rely on timing via the Global Positioning System (GPS) for the
synchronization of individual detector station clocks. Unfortunately, GPS
timing is expected to have an accuracy no better than about 5 ns. In practice,
in particular in AERA, the GPS clocks exhibit drifts on the order of tens of
ns. We developed a technique to correct for the GPS drifts, and an independent
method is used for cross-checks that indeed we reach nanosecond-scale timing
accuracy by this correction. First, we operate a "beacon transmitter" which
emits defined sine waves detected by AERA antennas recorded within the physics
data. The relative phasing of these sine waves can be used to correct for GPS
clock drifts. In addition to this, we observe radio pulses emitted by
commercial airplanes, the position of which we determine in real time from
Automatic Dependent Surveillance Broadcasts intercepted with a software-defined
radio. From the known source location and the measured arrival times of the
pulses we determine relative timing offsets between radio detector stations. We
demonstrate with a combined analysis that the two methods give a consistent
timing calibration with an accuracy of 2 ns or better. Consequently, the beacon
method alone can be used in the future to continuously determine and correct
for GPS clock drifts in each individual event measured by AERA.
-
This paper presents the results of different searches for correlations
between very high-energy neutrino candidates detected by IceCube and the
highest-energy cosmic rays measured by the Pierre Auger Observatory and the
Telescope Array. We first consider samples of cascade neutrino events and of
high-energy neutrino-induced muon tracks, which provided evidence for a
neutrino flux of astrophysical origin, and study their cross-correlation with
the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular
separation. We also study their possible directional correlations using a
likelihood method stacking the neutrino arrival directions and adopting
different assumptions on the size of the UHECR magnetic deflections. Finally,
we perform another likelihood analysis stacking the UHECR directions and using
a sample of through-going muon tracks optimized for neutrino point-source
searches with sub-degree angular resolution. No indications of correlations at
discovery level are obtained for any of the searches performed. The smallest of
the p-values comes from the search for correlation between UHECRs with IceCube
high-energy cascades, a result that should continue to be monitored.
-
Joint contributions of the Pierre Auger Collaboration and the Telescope Array
Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6
August 2015, The Hague, The Netherlands.
-
We have conducted three searches for correlations between ultra-high energy
cosmic rays detected by the Telescope Array and the Pierre Auger Observatory,
and high-energy neutrino candidate events from IceCube. Two cross-correlation
analyses with UHECRs are done: one with 39 cascades from the IceCube
`high-energy starting events' sample and the other with 16 high-energy `track
events'. The angular separation between the arrival directions of neutrinos and
UHECRs is scanned over. The same events are also used in a separate search
using a maximum likelihood approach, after the neutrino arrival directions are
stacked. To estimate the significance we assume UHECR magnetic deflections to
be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and
$9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field
strength and UHECR charge. A similar analysis is performed on stacked UHECR
arrival directions and the IceCube sample of through-going muon track events
which were optimized for neutrino point-source searches.
-
Contributions of the Pierre Auger Collaboration to the 34th International
Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands
-
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.
-
We present experimental results on inclusive spectra and mean multiplicities
of negatively charged pions produced in inelastic p+p interactions at incident
projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7,
8.8, 12.3 and 17.3 GeV, respectively). The measurements were performed using
the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton
Synchrotron.
Two-dimensional spectra are determined in terms of rapidity and transverse
momentum. Their properties such as the width of rapidity distributions and the
inverse slope parameter of transverse mass spectra are extracted and their
collision energy dependences are presented. The results on inelastic p+p
interactions are compared with the corresponding data on central Pb+Pb
collisions measured by the NA49 experiment at the CERN SPS.
The results presented in this paper are part of the NA61/SHINE ion program
devoted to the study of the properties of the onset of deconfinement and search
for the critical point of strongly interacting matter. They are required for
interpretation of results on nucleus-nucleus and proton-nucleus collisions.
-
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.
-
Spectra of K0S mesons and Lambda hyperons were measured in p+C interactions
at 31 GeV/c with the large acceptance NA61/SHINE spectrometer at the CERN SPS.
The data were collected with an isotropic graphite target with a thickness of
4% of a nuclear interaction length. Interaction cross sections, charged pion
spectra, and charged kaon spectra were previously measured using the same data
set. Results on K0S and Lambda production in p+C interactions serve as
reference for the understanding of the enhancement of strangeness production in
nucleus-nucleus collisions. Moreover, they provide important input for the
improvement of neutrino flux predictions for the T2K long baseline neutrino
oscillation experiment in Japan. Inclusive production cross sections for K0S
and Lambda are presented as a function of laboratory momentum in intervals of
the laboratory polar angle covering the range from 0 up to 240 mrad. The
results are compared with predictions of several hadron production models. The
K0S mean multiplicity in production processes <n_K0S> and the inclusive cross
section for K0S production were measured and amount to 0.127 +- 0.005 (stat) +-
0.022 (sys) and 29.0 +- 1.6 (stat) +- 5.0 (sys) mb, respectively.
-
The T2K long-baseline neutrino oscillation experiment in Japan needs precise
predictions of the initial neutrino flux. The highest precision can be reached
based on detailed measurements of hadron emission from the same target as used
by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The
corresponding data were recorded in 2007-2010 by the NA61/SHINE experiment at
the CERN SPS using a replica of the T2K graphite target. In this paper details
of the experiment, data taking, data analysis method and results from the 2007
pilot run are presented. Furthermore, the application of the NA61/SHINE
measurements to the predictions of the T2K initial neutrino flux is described
and discussed.
-
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 modular hadron calorimeter with micro-pixel avalanche photodiodes readout
for the NA61/SHINE experiment at the CERN SPS is presented. The calorimeter
consists of 44 independent modules with lead-scintillator sandwich structure.
The light from the scintillator tiles is captured by and transported with
WLS-fibers embedded in scintillator grooves. The construction provides a
longitudinal segmentation of the module in 10 sections with independent MAPD
readout. MAPDs with pixel density of $~10^{4}$/mm$^2$ ensure good linearity of
calorimeter response in a wide dynamical range. The performance of the
calorimeter prototype in a beam test is reported.
-
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.
-
Spectra of positively charged kaons in p+C interactions at 31 GeV/c were
measured with the NA61/SHINE spectrometer at the CERN SPS. The analysis is
based on the full set of data collected in 2007 with a graphite target with a
thickness of 4% of a nuclear interaction length. Interaction cross sections and
charged pion spectra were already measured using the same set of data. These
new measurements in combination with the published ones are required to improve
predictions of the neutrino flux for the T2K long baseline neutrino oscillation
experiment in Japan. In particular, the knowledge of kaon production is crucial
for precisely predicting the intrinsic electron neutrino component and the high
energy tail of the T2K beam. The results are presented as a function of
laboratory momentum in 2 intervals of the laboratory polar angle covering the
range from 20 up to 240 mrad. The kaon spectra are compared with predictions of
several hadron production models. Using the published pion results and the new
kaon data, the K+/\pi+ ratios are computed.
-
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.
-
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.
Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.
The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.
The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.
To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used for cross-checks that indeed we reach nanosecond-scale timing accuracy by this correction. First, we operate a "beacon transmitter" which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.
This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.
Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.
We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.
Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands
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.
We present experimental results on inclusive spectra and mean multiplicities of negatively charged pions produced in inelastic p+p interactions at incident projectile momenta of 20, 31, 40, 80 and 158 GeV/c ($\sqrt{s} = $ 6.3, 7.7, 8.8, 12.3 and 17.3 GeV, respectively). The measurements were performed using the large acceptance NA61/SHINE hadron spectrometer at the CERN Super Proton Synchrotron. Two-dimensional spectra are determined in terms of rapidity and transverse momentum. Their properties such as the width of rapidity distributions and the inverse slope parameter of transverse mass spectra are extracted and their collision energy dependences are presented. The results on inelastic p+p interactions are compared with the corresponding data on central Pb+Pb collisions measured by the NA49 experiment at the CERN SPS. The results presented in this paper are part of the NA61/SHINE ion program devoted to the study of the properties of the onset of deconfinement and search for the critical point of strongly interacting matter. They are required for interpretation of results on nucleus-nucleus and proton-nucleus collisions.
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.
Spectra of K0S mesons and Lambda hyperons were measured in p+C interactions at 31 GeV/c with the large acceptance NA61/SHINE spectrometer at the CERN SPS. The data were collected with an isotropic graphite target with a thickness of 4% of a nuclear interaction length. Interaction cross sections, charged pion spectra, and charged kaon spectra were previously measured using the same data set. Results on K0S and Lambda production in p+C interactions serve as reference for the understanding of the enhancement of strangeness production in nucleus-nucleus collisions. Moreover, they provide important input for the improvement of neutrino flux predictions for the T2K long baseline neutrino oscillation experiment in Japan. Inclusive production cross sections for K0S and Lambda are presented as a function of laboratory momentum in intervals of the laboratory polar angle covering the range from 0 up to 240 mrad. The results are compared with predictions of several hadron production models. The K0S mean multiplicity in production processes <n_K0S> and the inclusive cross section for K0S production were measured and amount to 0.127 +- 0.005 (stat) +- 0.022 (sys) and 29.0 +- 1.6 (stat) +- 5.0 (sys) mb, respectively.
The T2K long-baseline neutrino oscillation experiment in Japan needs precise predictions of the initial neutrino flux. The highest precision can be reached based on detailed measurements of hadron emission from the same target as used by T2K exposed to a proton beam of the same kinetic energy of 30 GeV. The corresponding data were recorded in 2007-2010 by the NA61/SHINE experiment at the CERN SPS using a replica of the T2K graphite target. In this paper details of the experiment, data taking, data analysis method and results from the 2007 pilot run are presented. Furthermore, the application of the NA61/SHINE measurements to the predictions of the T2K initial neutrino flux is described and discussed.
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 modular hadron calorimeter with micro-pixel avalanche photodiodes readout for the NA61/SHINE experiment at the CERN SPS is presented. The calorimeter consists of 44 independent modules with lead-scintillator sandwich structure. The light from the scintillator tiles is captured by and transported with WLS-fibers embedded in scintillator grooves. The construction provides a longitudinal segmentation of the module in 10 sections with independent MAPD readout. MAPDs with pixel density of $~10^{4}$/mm$^2$ ensure good linearity of calorimeter response in a wide dynamical range. The performance of the calorimeter prototype in a beam test is reported.
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.
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.