-
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
-
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.
-
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.
-
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.
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.
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.
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.
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.