-
Silicon photomultipliers present saturation effects as they have a limited
number of pixels and work in Geiger mode. Their response to light pulses in the
nonlinear region is very complex for two reasons: pixel recharging after an
avalanche affects the trigger probability and charge multiplication of
subsequent avalanches, and non-trivial effects due to crosstalk and
afterpulsing.
A parametrization of the nonlinear response of silicon photomultipliers was
developed where the above effects were readily accounted for. The model was
tested on a setup of gamma-ray spectrometry using different combinations of
scintillation crystals and detectors. The model parameters were interpreted in
terms of fundamental characteristics of the setup (e.g., lifetime of the
scintillation crystal and pixel recovery time). The proper conversion from
signal resolution to energy resolution was provided.
-
It is generally assumed that fluorescence radiation does not play a
significant role in the performance of Cherenkov telescopes. However, this
assumption is required to be verified using detailed Monte Carlo simulations.
In order to do this, we have implemented the production and tracking of
fluorescence radiation inside the CORSIKA code, and simulated gamma-ray induced
showers in the very high energy range. The most accurate fluorescence-yield
data available so far was used for this purpose.
The distribution of both light components on the ground has been studied as a
function of various parameters affecting the detection and reconstruction of
gamma-ray showers such as the angular aperture. From these distributions, we
determined the conditions under which fluorescence radiation becomes
significant. These results will also be useful to estimate the corresponding
systematic errors in Cherenkov telescope observations. The full simulation
results have been cross-checked, on a small sample of events, against numerical
calculations based on a one-dimension shower profile and found to be compatible
with each other. Both tools can be used for further investigations, like
studying the possibility to modify Cherenkov telescopes for the measurement of
fluorescence induced by extensive air showers.
-
In gamma-ray astronomy, a variety of data formats and proprietary software
have been traditionally used, often developed for one specific mission or
experiment. Especially for ground-based imaging atmospheric Cherenkov
telescopes (IACTs), data and software are mostly private to the collaborations
operating the telescopes. However, there is a general movement in science
towards the use of open data and software. In addition, the next-generation
IACT instrument, the Cherenkov Telescope Array (CTA), will be operated as an
open observatory.
We have created a Github organisation at
https://github.com/open-gamma-ray-astro where we are developing high-level data
format specifications. A public mailing list was set up at
https://lists.nasa.gov/mailman/listinfo/open-gamma-ray-astro and a first
face-to-face meeting on the IACT high-level data model and formats took place
in April 2016 in Meudon (France). This open multi-mission effort will help to
accelerate the development of open data formats and open-source software for
gamma-ray astronomy, leading to synergies in the development of analysis codes
and eventually better scientific results (reproducible, multi-mission).
This write-up presents this effort for the first time, explaining the
motivation and context, the available resources and process we use, as well as
the status and planned next steps for the data format specifications. We hope
that it will stimulate feedback and future contributions from the gamma-ray
astronomy community.
-
We measure the energy emitted by extensive air showers in the form of radio
emission in the frequency range from 30 to 80 MHz. Exploiting the accurate
energy scale of the Pierre Auger Observatory, we obtain a radiation energy of
15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV
arriving perpendicularly to a geomagnetic field of 0.24 G, scaling
quadratically with the cosmic-ray energy. A comparison with predictions from
state-of-the-art first-principle calculations shows agreement with our
measurement. The radiation energy provides direct access to the calorimetric
energy in the electromagnetic cascade of extensive air showers. Comparison with
our result thus allows the direct calibration of any cosmic-ray radio detector
against the well-established energy scale of the Pierre Auger Observatory.
-
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger
Observatory and is used to detect the radio emission of cosmic-ray air showers.
These observations are compared to the data of the surface detector stations of
the Observatory, which provide well-calibrated information on the cosmic-ray
energies and arrival directions. The response of the radio stations in the 30
to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of
the incoming electric field. For the latter, the energy deposit per area is
determined from the radio pulses at each observer position and is interpolated
using a two-dimensional function that takes into account signal asymmetries due
to interference between the geomagnetic and charge-excess emission components.
The spatial integral over the signal distribution gives a direct measurement of
the energy transferred from the primary cosmic ray into radio emission in the
AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air
shower arriving perpendicularly to the geomagnetic field. This radiation energy
-- corrected for geometrical effects -- is used as a cosmic-ray energy
estimator. Performing an absolute energy calibration against the
surface-detector information, we observe that this radio-energy estimator
scales quadratically with the cosmic-ray energy as expected for coherent
emission. We find an energy resolution of the radio reconstruction of 22% for
the data set and 17% for a high-quality subset containing only events with at
least five radio stations with signal.
-
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are
detectable with the Surface Detector array of the Pierre Auger Observatory. We
report here on searches through Auger data from 1 January 2004 until 20 June
2013. No neutrino candidates were found, yielding a limit to the diffuse flux
of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound
predictions. Neutrino identification is attempted using the broad
time-structure of the signals expected in the SD stations, and is efficiently
done for neutrinos of all flavors interacting in the atmosphere at large zenith
angles, as well as for "Earth-skimming" neutrino interactions in the case of
tau neutrinos. In this paper the searches for downward-going neutrinos in the
zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for
upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L.
single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an
$E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times
10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~
s^{-1}~ sr^{-1}}$.
-
A measurement of the cosmic-ray spectrum for energies exceeding
$4{\times}10^{18}$ eV is presented, which is based on the analysis of showers
with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger
Observatory between 1 January 2004 and 31 December 2013. The measured spectrum
confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$
eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with
index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by
a smooth suppression region. For the energy ($E_\text{s}$) at which the
spectral flux has fallen to one-half of its extrapolated value in the absence
of suppression, we find
$E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$
eV.
-
The crosstalk and afterpulsing in Hamamatsu silicon photomultipliers, called
Multi-Pixel Photon Counters (MPPCs), have been studied in depth. Several
components of the correlated noise have been identified according to their
different possible causes and their effects on the signal. In particular, we
have distinguished between prompt and delayed crosstalk as well as between
trap-assisted and hole-induced afterpulsing. The prompt crosstalk has been
characterized through the pulse amplitude spectrum measured at dark conditions.
The newest MPPC series, which incorporate isolating trenches between pixels,
exhibit a very low prompt crosstalk, but a small component remains likely due
to secondary photons reflected on the top surface of the device and
photon-generated minority carriers diffusing in the silicon substrate. We
present a meticulous procedure to characterize the afterpulsing and delayed
crosstalk through the amplitude and delay time distributions of secondary
pulses. Our results indicate that both noise components are due to minority
carriers diffusing in the substrate and that this effect is drastically reduced
in the new MPPC series as a consequence of an increase of one order of
magnitude in the doping density of the substrate. Finally, we have developed a
Monte Carlo simulation to study the different components of the afterpulsing
and crosstalk. The simulation results support our interpretation of the
experimental data. They also demonstrate that trenches longer than those
employed in the Hamamatsu MPPCs would reduce the crosstalk to a much greater
extent.
-
An experimental method to characterize the crosstalk and afterpulsing in
silicon photomultipliers has been developed and applied to two detectors
fabricated by Hamamatsu. An analytical model of optical crosstalk that we
presented in a previous publication has been compared with new measurements,
confirming our results. Progresses on a statistical model to describe
afterpulsing and delayed crosstalk are also shown and compared with preliminary
experimental data.
-
Contributions of the Pierre Auger Collaboration to the 33rd International
Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
-
Pulsed UV lasers can provide useful "testbeams" for observatories that use
optical detectors, especially fluorescence detectors, to measure high energy
cosmic-rays. The light observed by the detector is proportional to the energy
of the laser pulse. Since the absolute laser energy can be measured locally, a
well-calibrated laser offers a practical way to test the photometric
calibration of the cosmic-ray detector including atmospheric corrections. This
poster will describe a robotic system for laser polarization and energy
calibration. Laboratory measurements of laser energies and polarizations by
energy probes from different manufactures will be presented
Silicon photomultipliers present saturation effects as they have a limited number of pixels and work in Geiger mode. Their response to light pulses in the nonlinear region is very complex for two reasons: pixel recharging after an avalanche affects the trigger probability and charge multiplication of subsequent avalanches, and non-trivial effects due to crosstalk and afterpulsing. A parametrization of the nonlinear response of silicon photomultipliers was developed where the above effects were readily accounted for. The model was tested on a setup of gamma-ray spectrometry using different combinations of scintillation crystals and detectors. The model parameters were interpreted in terms of fundamental characteristics of the setup (e.g., lifetime of the scintillation crystal and pixel recovery time). The proper conversion from signal resolution to energy resolution was provided.
It is generally assumed that fluorescence radiation does not play a significant role in the performance of Cherenkov telescopes. However, this assumption is required to be verified using detailed Monte Carlo simulations. In order to do this, we have implemented the production and tracking of fluorescence radiation inside the CORSIKA code, and simulated gamma-ray induced showers in the very high energy range. The most accurate fluorescence-yield data available so far was used for this purpose. The distribution of both light components on the ground has been studied as a function of various parameters affecting the detection and reconstruction of gamma-ray showers such as the angular aperture. From these distributions, we determined the conditions under which fluorescence radiation becomes significant. These results will also be useful to estimate the corresponding systematic errors in Cherenkov telescope observations. The full simulation results have been cross-checked, on a small sample of events, against numerical calculations based on a one-dimension shower profile and found to be compatible with each other. Both tools can be used for further investigations, like studying the possibility to modify Cherenkov telescopes for the measurement of fluorescence induced by extensive air showers.
In gamma-ray astronomy, a variety of data formats and proprietary software have been traditionally used, often developed for one specific mission or experiment. Especially for ground-based imaging atmospheric Cherenkov telescopes (IACTs), data and software are mostly private to the collaborations operating the telescopes. However, there is a general movement in science towards the use of open data and software. In addition, the next-generation IACT instrument, the Cherenkov Telescope Array (CTA), will be operated as an open observatory. We have created a Github organisation at https://github.com/open-gamma-ray-astro where we are developing high-level data format specifications. A public mailing list was set up at https://lists.nasa.gov/mailman/listinfo/open-gamma-ray-astro and a first face-to-face meeting on the IACT high-level data model and formats took place in April 2016 in Meudon (France). This open multi-mission effort will help to accelerate the development of open data formats and open-source software for gamma-ray astronomy, leading to synergies in the development of analysis codes and eventually better scientific results (reproducible, multi-mission). This write-up presents this effort for the first time, explaining the motivation and context, the available resources and process we use, as well as the status and planned next steps for the data format specifications. We hope that it will stimulate feedback and future contributions from the gamma-ray astronomy community.
We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
The crosstalk and afterpulsing in Hamamatsu silicon photomultipliers, called Multi-Pixel Photon Counters (MPPCs), have been studied in depth. Several components of the correlated noise have been identified according to their different possible causes and their effects on the signal. In particular, we have distinguished between prompt and delayed crosstalk as well as between trap-assisted and hole-induced afterpulsing. The prompt crosstalk has been characterized through the pulse amplitude spectrum measured at dark conditions. The newest MPPC series, which incorporate isolating trenches between pixels, exhibit a very low prompt crosstalk, but a small component remains likely due to secondary photons reflected on the top surface of the device and photon-generated minority carriers diffusing in the silicon substrate. We present a meticulous procedure to characterize the afterpulsing and delayed crosstalk through the amplitude and delay time distributions of secondary pulses. Our results indicate that both noise components are due to minority carriers diffusing in the substrate and that this effect is drastically reduced in the new MPPC series as a consequence of an increase of one order of magnitude in the doping density of the substrate. Finally, we have developed a Monte Carlo simulation to study the different components of the afterpulsing and crosstalk. The simulation results support our interpretation of the experimental data. They also demonstrate that trenches longer than those employed in the Hamamatsu MPPCs would reduce the crosstalk to a much greater extent.
An experimental method to characterize the crosstalk and afterpulsing in silicon photomultipliers has been developed and applied to two detectors fabricated by Hamamatsu. An analytical model of optical crosstalk that we presented in a previous publication has been compared with new measurements, confirming our results. Progresses on a statistical model to describe afterpulsing and delayed crosstalk are also shown and compared with preliminary experimental data.
Pulsed UV lasers can provide useful "testbeams" for observatories that use optical detectors, especially fluorescence detectors, to measure high energy cosmic-rays. The light observed by the detector is proportional to the energy of the laser pulse. Since the absolute laser energy can be measured locally, a well-calibrated laser offers a practical way to test the photometric calibration of the cosmic-ray detector including atmospheric corrections. This poster will describe a robotic system for laser polarization and energy calibration. Laboratory measurements of laser energies and polarizations by energy probes from different manufactures will be presented