Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2011). The exposure of the hybrid detector of the Pierre Auger Observatory. Astropart Phys., 34(6), 368–381.
Abstract: The Pierre Auger Observatory is a detector for ultra-high energy cosmic rays. It consists of a surface array to measure secondary particles at ground level and a fluorescence detector to measure the development of air showers in the atmosphere above the array. The “hybrid” detection mode combines the information from the two subsystems. We describe the determination of the hybrid exposure for events observed by the fluorescence telescopes in coincidence with at least one water-Cherenkov detector of the surface array. A detailed knowledge of the time dependence of the detection operations is crucial for an accurate evaluation of the exposure. We discuss the relevance of monitoring data collected during operations, such as the status of the fluorescence detector, background light and atmospheric conditions, that are used in both simulation and reconstruction.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2011). Search for first harmonic modulation in the right ascension distribution of cosmic rays detected at the Pierre Auger Observatory. Astropart Phys., 34(8), 627–639.
Abstract: We present the results of searches for dipolar-type anisotropies in different energy ranges above 2.5 x 10(17) eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% C.L. for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2011). Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory. Nucl. Instrum. Methods Phys. Res. A, 635(1), 92–102.
Abstract: The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs “radio-hybrid” measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluorescence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.
|
Menjo, H. et al, Faus-Golfe, A., & Velasco, J. (2011). Monte Carlo study of forward pi(0) production spectra to be measured by the LHCf experiment for the purpose of benchmarking hadron interaction models at 10(17) eV. Astropart Phys., 34(7), 513–520.
Abstract: The LHCf experiment aims to improve knowledge of forward neutral particle production spectra at the LHC energy which is relevant for the interpretation of air shower development of high energy cosmic rays. Two detectors, each composed of a pair of sampling and imaging calorimeters, have been installed at the forward region of IP1 to measure pi(0) energy spectra above 600 GeV. In this paper, we present a Monte Carlo study of the pi(0) measurements to be performed with one of the LHCf detectors for proton-proton collisions at root s = 14 TeV. In approximately 40 min of operation at luminosity 0.8 x 10(29) cm(-2) s(-1) during the beam commissioning phase of LHC, about 1.5 x 10(4) pi(0) events are expected to be obtained at two transverse positions of the detector. The backgrounds from interactions of secondary particles with beam pipes and interactions of beam particles with residual gas in the beam pipes are expected to be less than 0.1% of the signal from pi(0)s. We also discuss the capability of LHCf measurements to discriminate between the various hadron interaction models that are used for simulation of high energy air showers, such as DPMJET3.03, QGSJETII-03, SIBYLL2.1 and EPOS1.99.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2011). The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray showers detected by the Pierre Auger Observatory. Astropart Phys., 35(5), 266–276.
Abstract: 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 degrees. 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
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2012). Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory. Astropart Phys., 35(6), 354–361.
Abstract: 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.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2012). Description of atmospheric conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS). Astropart Phys., 35(9), 591–607.
Abstract: 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 Malargue and averaged monthly models, the utility of the GDAS data is shown.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2012). A search for anisotropy in the arrival directions of ultra high energy cosmic rays recorded at the Pierre Auger Observatory. J. Cosmol. Astropart. Phys., 04(4), 040–21pp.
Abstract: Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2012). Search for Point-like Sources of Ultra-high Energy Neutrinos at the Pierre Auger Observatory and Improved Limit on the Diffuse Flux of Tau Neutrinos. Astrophys. J. Lett., 755(1), L4–7pp.
Abstract: The surface detector array of the Pierre Auger Observatory can detect neutrinos with energy E-nu between 10(17) eV and 10(20) eV from point-like sources across the sky south of +55 degrees and north of -65 degrees declinations. A search has been performed for highly inclined extensive air showers produced by the interaction of neutrinos of all flavors in the atmosphere (downward-going neutrinos), and by the decay of tau leptons originating from tau neutrino interactions in Earth's crust (Earth-skimming neutrinos). No candidate neutrinos have been found in data up to 2010 May 31. This corresponds to an equivalent exposure of similar to 3.5 years of a full surface detector array for the Earth-skimming channel and similar to 2 years for the downward-going channel. An improved upper limit on the diffuse flux of tau neutrinos has been derived. Upper limits on the neutrino flux from point-like sources have been derived as a function of the source declination. Assuming a differential neutrino flux k(PS) . E-nu(-2). from a point-like source, 90% confidence level upper limits for k(PS) at the level of approximate to 5x10(-7) and 2.5x10(-6) GeV cm(-2) s(-1) have been obtained over a broad range of declinations from the searches for Earth-skimming and downward-going neutrinos, respectively.
|
