• ### Description of longitudinal profiles of showers dominated by Cherenkov light(1709.01458)

Sept. 5, 2017 astro-ph.HE
With the aim to describe the longitudinal development of Cherenkov dominated showers we investigate the energy deposit and the number of charged particles in air showers induced by energetic cosmic rays. Based on the Monte Carlo simulations, discrepancies between different estimates of calorimetric energies are documented. We focus on the energy deposit profiles of air showers deducible from the fluorescence and Cherenkov light generated along CONEX and CORSIKA cascades.
• ### Mass Composition of Cosmic Rays with Combined Surface Detector Arrays(1708.06164)

Aug. 21, 2017 astro-ph.HE
Our study exploits the Constant Intensity Cut principles applied simultaneously to muonic and electromagnetic detectors of cosmic rays. We use the fact that the ordering of events according to their signal sizes induced in different types of surface detectors provides information about the mass composition of primary cosmic-ray beam, with low sensitivity to details of hadronic interactions. Composition analysis at knee energies is performed using Monte Carlo simulations for extensive air showers having maxima located far away from a hypothetical observatory. Another type of a hypothetical observatory is adopted to examine composition of ultra-high energy primaries which initiate vertical air showers with maxima observed near surface detectors.
• ### A Bayesian on-off analysis of cosmic ray data(1707.03155)

July 11, 2017 astro-ph.IM
We deal with the analysis of on-off measurements designed for the confirmation of a weak source of events whose presence is hypothesized, based on former observations. The problem of a small number of source events that are masked by an imprecisely known background is addressed from a Bayesian point of view. We examine three closely related variables, the posterior distributions of which carry relevant information about various aspects of the investigated phenomena. This information is utilized for predictions of further observations, given actual data. Backed by details of detection, we propose how to quantify disparities between different measurements. The usefulness of the Bayesian inference is demonstrated on examples taken from cosmic ray physics.
• 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.
• ### On Bayesian analysis of on-off measurements(1603.03386)

March 10, 2016 astro-ph.IM
We propose an analytical solution to the on-off problem within the framework of Bayesian statistics. Both the statistical significance for the discovery of new phenomena and credible intervals on model parameters are presented in a consistent way. We use a large enough family of prior distributions of relevant parameters. The proposed analysis is designed to provide Bayesian solutions that can be used for any number of observed on-off events, including zero. The procedure is checked using Monte Carlo simulations. The usefulness of the method is demonstrated on examples from gamma-ray astronomy.
• ### Maximum entropy analysis of cosmic ray composition(1512.09248)

Jan. 1, 2016 astro-ph.HE
We focus on the primary composition of cosmic rays with the highest energies that cause extensive air showers in the Earth's atmosphere. A way of examining the two lowest order moments of the sample distribution of the depth of shower maximum is presented. The aim is to show that useful information about the composition of the primary beam can be inferred with limited knowledge we have about processes underlying these observations. In order to describe how the moments of the depth of shower maximum depend on the type of primary particles and their energies, we utilize a superposition model. Using the principle of maximum entropy, we are able to determine what trends in the primary composition are consistent with the input data, while relying on a limited amount of information from shower physics. Some capabilities and limitations of the proposed method are discussed. In order to achieve a realistic description of the primary mass composition, we pay special attention to the choice of the parameters of the superposition model. We present two examples that demonstrate what consequences can be drawn for energy dependent changes in the primary composition.
• ### Variability of VHE $\gamma$-ray emission from the binary PSR B1259-63/LS 2883(1512.04849)

Dec. 15, 2015 astro-ph.HE
We examine changes of the $\gamma$-ray intensity observed from the direction of the binary system PSR B1259-63/LS 2883 during campaigns around its three periastron passages. A simple and straightforward method is applied to the published data obtained with the Imaging Atmospheric Cherenkov Technique. Regardless of many issues of the detection process, the method works only with numbers of very high energetic photons registered in the specified regions. Within the realm of this scheme, we recognized changes attributable to the variations of the intrinsic source activity at high levels of significance.
• 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.
• ### Significance for signal changes in gamma-ray astronomy(1509.00353)

Sept. 1, 2015 astro-ph.IM, astro-ph.HE
We describe a straightforward modification of frequently invoked methods for the determination of the statistical significance of a gamma-ray signal observed in a counting process. A simple criterion is proposed to decide whether a set of measurements of the numbers of photons registered in the source and background regions is consistent with the assumption of a constant source activity. This method is particularly suitable for immediate evaluation of the stability of the observed gamma-ray signal. It is independent of the exposure estimates, reducing thus the impact of systematic inaccuracies, and properly accounts for the fluctuations in the number of detected photons. The usefulness of the method is demonstrated on several examples. We discuss intensity changes for gamma-ray emitters detected at very high energies by the current gamma-ray telescopes (e.g. 1ES 0229+200, 1ES 1959+650 and PG 1553+113). Some of the measurements are quantified to be exceptional with large statistical significances.
• ### A branching model for hadronic air showers(1509.00364)

Sept. 1, 2015 astro-ph.HE
We introduce a simple branching model for the development of hadronic showers in the Earth's atmosphere. Based on this model, we show how the size of the pionic component followed by muons can be estimated. Several aspects of the subsequent muonic component are also discussed. We focus on the energy evolution of the muon production depth. We also estimate the impact of the primary particle mass on the size of the hadronic component. Even though a precise calculation of the development of air showers must be left to complex Monte Carlo simulations, the proposed model can reveal qualitative insight into the air shower physics.
• ### Unexpected gamma-ray signal in the vicinity of 1ES 0229+200(1509.00333)

Sept. 1, 2015 astro-ph.HE
We report on an unidentified gamma-ray signal found in the region around the BL Lac object 1ES 0229+200. It was recognized serendipitously in our analysis of 6.2 years of Fermi-LAT data at a distance less than 3{\deg} away from the blazar. The observed excess of counts manifests itself as an unexpected local maximum in the test statistic map. Although several Fermi-LAT sources have been identified in this area we were not able to link them to the position of this residual signal. A clear association with sources visible in other wavebands was not successful either. We briefly discuss characteristics of this unresolved phenomenon. Our results suggest a steep energy spectrum and a point-like nature of this candidate gamma-ray emitter.
• ### Study of Dispersion of Mass Distribution of Ultra-High Energy Cosmic Rays using a Surface Array of Muon and Electromagnetic Detectors(1503.07734)

March 26, 2015 astro-ph.IM, astro-ph.HE
We consider a hypothetical observatory of ultra-high energy cosmic rays consisting of two surface detector arrays that measure independently electromagnetic and muon signals induced by air showers. Using the constant intensity cut method, sets of events ordered according to each of both signal sizes are compared giving the number of matched events. Based on its dependence on the zenith angle, a parameter sensitive to the dispersion of the distribution of the logarithmic mass of cosmic rays is introduced. The results obtained using two post-LHC models of hadronic interactions are very similar and indicate a weak dependence on details of these interactions.
• ### Variability of VHE $\gamma$-ray sources(1412.2050)

Dec. 5, 2014 astro-ph.HE
We study changes in the $\gamma$-ray intensity at very high energies observed from selected active galactic nuclei. Publicly available data collected by Cherenkov telescopes were examined by means of a simple method utilizing solely the number of source and background events. Our results point to some degree of time variability in signal observed from the investigated sources. Several measurements were found to be excessive or deficient in the number of source events when compared to the source intensity deduced from other observations.
• ### Time variability of the $\gamma$-ray binary HESS J0632+057(1409.5395)

Sept. 18, 2014 astro-ph.HE
We study changes in the $\gamma$--ray intensity at very high energies observed from the $\gamma$--ray binary HESS J0632+057. Publicly available data collected by Cherenkov telescopes were examined by means of a simple method utilizing solely the number of source and background events. Our results point to time variability in signal from the selected object consistent with periodic modulation of the source intensity.
• ### Testing time variability of gamma-ray flux(1309.6476)

Sept. 25, 2013 astro-ph.IM, astro-ph.HE
A way of examining a hypothetical non--zero $\gamma$--ray signal for the time changes is presented. The time variability of the recently observed $\gamma$--ray source PKS 2155--304 is discussed. Several measurements were found to be excessive or deficient with large significances on time scales of months and days.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### New method for atmospheric calibration at the Pierre Auger Observatory using FRAM, a robotic astronomical telescope(0706.1710)

June 12, 2007 astro-ph
FRAM - F/(Ph)otometric Robotic Atmospheric Monitor is the latest addition to the atmospheric monitoring instruments of the Pierre Auger Observatory. An optical telescope equipped with CCD camera and photometer, it automatically observes a set of selected standard stars and a calibrated terrestrial source. Primarily, the wavelength dependence of the attenuation is derived and the comparison between its vertical values (for stars) and horizontal values (for the terrestrial source) is made. Further, the integral vertical aerosol optical depth can be obtained. A secondary program of the instrument, the detection of optical counterparts of gamma-ray bursts, has already proven successful. The hardware setup, software system, data taking procedures, and first analysis results are described in this paper.
• ### The bright optical flash from GRB 060117(astro-ph/0606004)

June 2, 2006 astro-ph
We present a discovery and observation of an extraordinarily bright prompt optical emission of the GRB 060117 obtained by a wide-field camera atop the robotic telescope FRAM of the Pierre Auger Observatory from 2 to 10 minutes after the GRB. We found rapid average temporal flux decay of alpha = -1.7 +- 0.1 and a peak brightness R = 10.1 mag. Later observations by other instruments set a strong limit on the optical and radio transient fluxes, unveiling an unexpectedly rapid further decay. We present an interpretation featuring a relatively steep electron-distribution parameter p ~ 3.0 and providing a straightforward solution for the overall fast decay of this optical transient as a transition between reverse and forward shock.