• ### The Multi-Mission Maximum Likelihood framework (3ML)(1507.08343)

July 29, 2015 astro-ph.IM, astro-ph.HE
Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### Sensitivity of the correlation between the depth of shower maximum and the muon shower size to the cosmic ray composition(1203.3732)

March 16, 2012 astro-ph.HE
The composition of ultra-high energy cosmic rays is an important issue in astroparticle physics research, and additional experimental results are required for further progress. Here we investigate what can be learned from the statistical correlation factor r between the depth of shower maximum and the muon shower size, when these observables are measured simultaneously for a set of air showers. The correlation factor r contains the lowest-order moment of a two-dimensional distribution taking both observables into account, and it is independent of systematic uncertainties of the absolute scales of the two observables. We find that, assuming realistic measurement uncertainties, the value of r can provide a measure of the spread of masses in the primary beam. Particularly, one can differentiate between a well-mixed composition (i.e., a beam that contains large fractions of both light and heavy primaries) and a relatively pure composition (i.e., a beam that contains species all of a similar mass). The number of events required for a statistically significant differentiation is ~ 200. This differentiation, though diluted, is maintained to a significant extent in the presence of uncertainties in the phenomenology of high energy hadronic interactions. Testing whether the beam is pure or well-mixed is well motivated by recent measurements of the depth of shower maximum.
• ### A method to constrain the characteristic angular size of the brightest cosmic-ray sources observed above 57 $\times$ 10^{18} eV(0909.3483)

Sept. 26, 2011 astro-ph.CO, astro-ph.HE
We introduce a method to constrain the characteristic angular size of the brightest cosmic-ray sources observed above 57 \times 1018 eV. By angular size of a source, we mean the effective angular extent over which cosmic-rays from that source arrive at earth. The method is based on the small-scale (< 10\circ) self-clustering of cosmic-ray arrival directions. The method is applicable to sparse data sets in which strong localizations of CR* directions are not yet observed. We show that useful constraints on the source size can be made in the near future and that these constraints are not strongly dependent on the assumed spatial distribution and luminosity function of the cosmic-ray sources. We suggest that an indication of the source size is quite telling. For example, an indication of the source size can be used to infer limits on the particle charge and intervening magnetic fields (not independently), both of which are not well constrained so far. This is possible because the source size is similar in scale to the magnetic deflection.
• ### On Estimating the Flux of the Brightest Cosmic Ray Source above 57x10^18 eV(0903.3981)

March 23, 2009 astro-ph.HE
The sources of ultra-high energy cosmic rays are not yet known. However, the discovery of anisotropic cosmic rays above 57x10^18 eV by the Pierre Auger Observatory suggests that a direct source detection may soon be possible. The near-future prospects for such a measurement are heavily dependent on the flux of the brightest source. In this work, we show that the flux of the brightest source above 57x10^18 eV is expected to comprise 10% or more of the total flux if two general conditions are true. The conditions are: 1.) the source objects are associated with galaxies other than the Milky Way and its closest neighbors, and 2.) the cosmic ray particles are protons or heavy nuclei such as iron and the Greisen-Zatsepin-Kuz'min effect is occurring. The Pierre Auger Observatory collects approximately 23 events above 57x10^18 eV per year. Therefore, it is plausible that, over the course of several years, tens of cosmic rays from a single source will be detected.