Analysis of cosmic ray data

K.J.Orford

DOI: https://doi.org/10.1088/0954-3899/26/4/201

2000-02-07

Abstract:Reviews statistical methods for 1- and 2-dimensional analysis of cosmic ray data including recent Bayesian analyses.

Astrophysics

What problem does this paper attempt to address?

Statistical Methods for Investigating the Cosmic Ray Energy Spectrum

J. D. Hague,B. R. Becker,M. S. Gold,J. A. J. Matthews,J. Urbář

DOI: https://doi.org/10.48550/arXiv.0710.3600

2007-10-19

Abstract:Two separate statistical tests are described and developed in order to test un-binned data sets for adherence to the power-law form. The first test employs the TP-statistic, a function defined to deviate from zero when the sample deviates from the power-law form, regardless of the value of the power index. The second test employs a likelihood ratio test to reject a power-law background in favor of a model signal distribution with a cut-off.

Astrophysics,Statistics Theory
Statistical Tools for Analyzing the Cosmic Ray Energy Spectrum

J. D. Hague,B. R. Becker,M. S. Gold,J. A. J. Matthews

DOI: https://doi.org/10.48550/arXiv.0808.2685

2008-08-20

Abstract:In this paper un-binned statistical tools for analyzing the cosmic ray energy spectrum are developed and illustrated with a simulated data set. The methods are designed to extract accurate and precise model parameter estimators in the presence of statistical and systematic energy errors. Two robust methods are used to test for the presence of flux suppression at the highest energies: the Tail-Power statistic and a likelihood ratio test. Both tests give evidence of flux suppression in the simulated data. The tools presented can be generalized for use on any astrophysical data set where the power-law assumption is relevant and can be used to aid observational design.

Astrophysics
Techniques for Searching for Anisotropy in the Arrival Directions of Ultra-High Energy Cosmic Rays

dos Santos, Diego C.,Lago, Bruno L.

DOI: https://doi.org/10.1007/s13538-024-01640-6

2024-11-09

Brazilian Journal of Physics

Abstract:Understanding the origin and distribution of ultra-high-energy cosmic rays (UHECRs) remains a cornerstone challenge in astrophysics. These cosmic rays, with energies exceeding 8 EeV, provide insights into the most energetic phenomena in the universe. In this paper, we review and apply various statistical and analytical techniques to a simulated dataset to investigate the anisotropy in the arrival directions of UHECRs with energies exceeding 8 EeV. We simulate cosmic-ray events based on the distribution of local matter from the 2MASS Redshift Survey, supplemented with an isotropic background component, to create a robust dataset for analysis. Our methods include blind searches for excesses over the sky, Rayleigh analysis in right ascension, reconstruction of dipole and quadrupole moments using spherical harmonics, and the angular power spectrum technique — a method extensively utilized in cosmic microwave background studies — to assess anisotropies across multiple angular scales. The results underscore the potential of these methodologies to reveal significant anisotropies and contribute to our understanding of the possible extragalactic sources and propagation mechanisms of UHECRs.

physics, multidisciplinary
II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements

Nelson Christensen,Renate Meyer,Lloyd Knox,Ben Luey

DOI: https://doi.org/10.1088/0264-9381/18/14/306

2001-03-08

Abstract:We present a strategy for a statistically rigorous Bayesian approach to the problem of determining cosmological parameters from the results of observations of anisotropies in the cosmic microwave background. Our strategy relies on Markov chain Monte Carlo methods, specifically the Metropolis-Hastings algorithm, to perform the necessary high-dimensional integrals. We describe the Metropolis-Hastings algorithm in detail and discuss the results of our test on simulated data.

Astrophysics
Data analysis methods for the cosmic microwave background

Matthieu Tristram,Ken Ganga

DOI: https://doi.org/10.1088/0034-4885/70/6/R02

2008-06-10

Abstract:In this review, we give an overview of some of the major aspects of data reduction and analysis for the cosmic microwave background (CMB). Since its prediction and discovery in the last century, the CMB radiation has proven itself to be one of our most valuable tools for precision cosmology. Recently, and especially when combined with complementary cosmological data, measurements of the CMB anisotropies have provided us with a wealth of quantitive information about the birth, evolution and structure of our Universe. We begin with a simple, general introduction to the physics of the CMB, including a basic overview of the experiments which record CMB data. The focus, however, will be the data analysis treatment of CMB data sets.

Astrophysics
Novel Simulation Framework for Analyzing Cosmic Ray Particle Distributions at a Global Scale

Olesya Sarajlic,Xiaochun He

2024-11-05

Abstract:Cosmic ray measurements have inspired numerous interesting applications over several decades worldwide. These applications encompass non-invasive cosmic ray muon tomography, which enables the imaging of concealed dense objects or structures, the monitoring of area-averaged soil moisture with cosmic ray neutrons in agriculture and climate studies, real-time monitoring of the dynamical changes of the space and earth weather, etc. The demand for a quantitative characterization of cosmic ray shower particles near the Earth's surface is substantial, as it provides realistic particle spectra and rates for these diverse applications. In this study, we introduce Earth Cosmic Ray Shower (ECRS), a GEANT4-based software designed to simulate cosmic ray particle interactions in the atmosphere. ECRS incorporates the U.S. Standard Atmospheric Model and integrates a time-dependent geomagnetic field based on the Tsyganenko and IGRF models. Additionally, we present two case studies illustrating variations in the location-dependent average particle energy for muons, electrons, neutrons, and gammas at sea level. An outlook of this project is provided toward the conclusion.

High Energy Astrophysical Phenomena,Earth and Planetary Astrophysics,Instrumentation and Methods for Astrophysics,Atmospheric and Oceanic Physics,Space Physics
Analysis Methods for Gamma-ray Astronomy

Denys Malyshev,Lars Mohrmann

DOI: https://doi.org/10.1007/978-981-16-4544-0_177-1

2023-09-06

Abstract:The launch of the Fermi satellite in 2008, with its Large Area Telescope (LAT) on board, has opened a new era for the study of gamma-ray sources at GeV ($10^9$ eV) energies. Similarly, the commissioning of the third generation of imaging atmospheric Cherenkov telescopes (IACTs) - H.E.S.S., MAGIC, and VERITAS - in the mid-2000's has firmly established the field of TeV ($10^{12}$ eV) gamma-ray astronomy. Together, these instruments have revolutionised our understanding of the high-energy gamma-ray sky, and they continue to provide access to it over more than six decades in energy. In recent years, the ground-level particle detector arrays HAWC, Tibet, and LHAASO have opened a new window to gamma rays of the highest energies, beyond 100 TeV. Soon, next-generation facilities such as CTA and SWGO will provide even better sensitivity, thus promising a bright future for the field. In this chapter, we provide a brief overview of methods commonly employed for the analysis of gamma-ray data, focusing on those used for Fermi-LAT and IACT observations. We describe the standard data formats, explain event reconstruction and selection algorithms, and cover in detail high-level analysis approaches for imaging and extraction of spectra, including aperture photometry as well as advanced likelihood techniques.

Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena
Reanalysis of the Systematic Uncertainties in Cosmic-Ray Antiproton Flux

Xing-Jian Lv,Xiao-Jun Bi,Kun Fang,Peng-Fei Yin,Meng-Jie Zhao

DOI: https://doi.org/10.1103/PhysRevD.109.043006

2024-02-12

Abstract:Recent studies on cosmic rays (CRs) have reported the possibility of an excess in the antiproton flux around $10-20$ GeV. However, the associated systematic uncertainties have impeded the interpretation of these findings. In this study, we conduct a global Bayesian analysis to constrain the propagation parameters and evaluate the CR antiproton spectrum, while comprehensively accounting for uncertainties associated with interstellar CR propagation, production cross sections for antiprotons and other secondaries, and the charge and energy dependent effects of solar modulation. We establish that the most recent AMS-02 $\bar{p}$ spectrum is in agreement with a pure secondary origin. Based on this, we establish upper limits on dark matter (DM) annihilation. We also determine that the AMS-02 data favors the empirical hadronic interaction models over phenomenological ones. Finally, we find that the latest AMS-02 antiproton data from 2011 to 2018 disfavors the antiproton excess at $\mathcal{O}$(10) GeV and the corresponding DM interpretation that can simultaneously account for the Galactic Center excess in the gamma-ray observation.

High Energy Astrophysical Phenomena,High Energy Physics - Phenomenology
Reanalysis of the systematic uncertainties in the cosmic-ray antiproton flux

Xing-Jian Lv,Xiao-Jun Bi,Kun Fang,Peng-Fei Yin,Meng-Jie Zhao

DOI: https://doi.org/10.1103/physrevd.109.043006

IF: 5.407

2024-02-07

Physical Review D

Abstract:Recent studies on cosmic rays (CRs) have reported the possibility of an excess in the antiproton flux around 10–20 GeV. However, the associated systematic uncertainties have impeded the interpretation of these findings. In this study, we conduct a global Bayesian analysis to constrain the propagation parameters and evaluate the CR antiproton spectrum, while comprehensively accounting for uncertainties associated with interstellar CR propagation, production cross sections for antiprotons and other secondaries, and the charge- and energy-dependent effects of solar modulation. We establish that the most recent Alpha Magnetic Spectrometer (AMS-02) p ̄ spectrum is in agreement with a pure secondary origin. Based on this, we establish upper limits on dark matter (DM) annihilation. We also determine that the AMS-02 data favor the empirical hadronic interaction models over phenomenological ones. Finally, we find that the latest AMS-02 antiproton data from 2011 to 2018 disfavors the antiproton excess at O(10) GeV and the corresponding DM interpretation that can simultaneously account for the Galactic Center excess in the gamma-ray observation. https://doi.org/10.1103/PhysRevD.109.043006 © 2024 American Physical Society

astronomy & astrophysics,physics, particles & fields
Analysis and packaging of radiochemical solar neutrino data. 1. Bayesian approach

P.A. Sturrock,M.S. Wheatland

DOI: https://doi.org/10.1007/s11207-007-9112-5

2007-10-26

Abstract:According to current practice, the results of each run of a radiochemical solar neutrino experiment comprise an estimate of the flux and upper and lower error estimates. These estimates are derived by a maximum-likelihood procedure from the times of decay events in the analysis chamber. This procedure has the following shortcomings: (a) Published results sometimes include negative flux estimates. (b) Even if the flux estimate is non-negative, the probability distribution function implied by the flux and error estimates will extend into negative territory; and (c) The overall flux estimate derived from the results of a sequence of runs may differ substantially from an estimate made by a global analysis of all of the timing data taken together. These defects indicate that the usual packaging of data in radiochemical solar neutrino experiments provides an inadequate summary of the data, which implies a loss of information. This article reviews this problem from a Bayesian perspective, and suggests an alternative scheme for the packaging of radiochemical solar neutrino data, which is we believe free from the above objections.

High Energy Physics - Phenomenology,Astrophysics,General Relativity and Quantum Cosmology
Statistical Issues in Neutrino Physics Analyses

Louis Lyons

DOI: https://doi.org/10.48550/arXiv.1607.03549

2016-07-13

Abstract:Various statistical issues relevant to searches for new physics or to parameter determination in analyses of data in neutrino experiments are briefly discussed.

High Energy Physics - Experiment,Data Analysis, Statistics and Probability
Data-driven analysis for understanding ultrahigh energy cosmic ray source spectra

Marco Stein Muzio,Luis A. Anchordoqui,Michael Unger

DOI: https://doi.org/10.22323/1.444.0293

2023-08-14

Abstract:One of the most challenging open questions regarding the origin of ultrahigh energy cosmic rays (UHECRs) deals with the shape of the source emission spectra. A commonly-used simplifying assumption is that the source spectra of the highest energy cosmic rays trace a Peters cycle, in which the maximum cosmic-ray energy scales linearly with $Z$, i.e., with the charge of the UHECR in units of the proton charge. However, this would only be a natural assumption for models in which UHECRs escape the acceleration region without suffering significant energy losses. In most cases, however, UHECRs interact in the acceleration region and/or in the source environment changing the shape of the source emission spectra. Energy losses are typically parameterized in terms of $Z$ and the UHECR baryon number $A$, and therefore one would expect the source emission spectra to be a function of both $Z$ and $A$. Taking a pragmatic approach, we investigate whether existing data favor any region of the $(Z,A)$ parameter space. Using data from the Pierre Auger Observatory, we carry out a maximum likelihood analysis of the observed spectrum and nuclear composition to shape the source emission spectra for the various particle species. We also study the impact of possible systematic uncertainties driven by hadronic models describing interactions in the atmosphere.

High Energy Astrophysical Phenomena
Stochastic modelling of cosmic ray sources for diffuse high-energy gamma-rays and neutrinos

Anton Stall,Leonard Kaiser,Philipp Mertsch

DOI: https://doi.org/10.22323/1.444.0687

2023-09-06

Abstract:Cosmic rays of energies up to a few PeV are believed to be of galactic origin, yet individual sources have still not been firmly identified. Due to inelastic collisions with the interstellar gas, cosmic-ray nuclei produce a diffuse flux of high-energy gamma-rays and neutrinos. Fermi-LAT has provided maps of galactic gamma-rays at GeV energies which can be produced by both hadronic and leptonic processes. Neutrinos, on the other hand, are exclusively produced by the sought-after hadronic processes, yet they can be detected above backgrounds only at hundreds of TeV. Oftentimes, diffuse emission maps are extrapolated from GeV to PeV energies, but the sources contributing at either energies likely differ. We have modelled the production of diffuse emission from GeV through PeV energies in a Monte Carlo approach, taking into consideration the discrete nature of sources. We can generate realisations of the diffuse sky in a matter of seconds, thus allowing for characterising correlations in direction and energy. At hundreds of TeV, relevant for observations with LHAASO, Tibet AS-gamma, IceCube and the upcoming SWGO, variations between different realisations are sizeable. Specifically, we show that extrapolations of diffuse emission from GeV to PeV energies must fail and apply our results on the recent experimental findings.

High Energy Astrophysical Phenomena,Astrophysics of Galaxies
Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy

Sanjib Sharma

DOI: https://doi.org/10.1146/annurev-astro-082214-122339

IF: 37.226

2017-08-18

Annual Review of Astronomy and Astrophysics

Abstract:Markov chain Monte Carlo–based Bayesian data analysis has now become the method of choice for analyzing and interpreting data in almost all disciplines of science. In astronomy, over the past decade, we have also seen a steady increase in the number of papers that employ Monte Carlo–based Bayesian analysis. New, efficient Monte Carlo–based methods are continuously being developed and explored. In this review, we first explain the basics of Bayesian theory and discuss how to set up data analysis problems within this framework. Next, we provide an overview of various Monte Carlo–based methods for performing Bayesian data analysis. Finally, we discuss advanced ideas that enable us to tackle complex problems and thus hold great promise for the future. We also distribute downloadable computer software ( https://github.com/sanjibs/bmcmc ) Python that implements some of the algorithms and examples discussed here.

astronomy & astrophysics
A novel method of data analysis for hadronic physics

C. N. Papanicolas,E. Stiliaris

DOI: https://doi.org/10.48550/arXiv.1205.6505

2012-05-29

High Energy Physics - Phenomenology

Abstract:A novel method for extracting physical parameters from experimental and simulation data is presented. The method is based on statistical concepts and it relies on Monte Carlo simulation techniques. It identifies and determines with maximal precision parameters that are sensitive to the data. The method has been extensively studied and it is shown to produce unbiased results. It is applicable to a wide range of scientific and engineering problems. It has been successfully applied in the analysis of experimental data in hadronic physics and of lattice QCD correlators.
What Do We Expect to Learn from Very High Energy Cosmic Ray Observations?

H. Rebel

DOI: https://doi.org/10.48550/arXiv.astro-ph/0109520

2001-09-27

Astrophysics

Abstract:By a short tour through the exciting field of very high-energy and ultra high-energy cosmic rays studies, a brief review is given about the current questions approached, in particular by the KASCADE experiment and the Pierre Auger project. The present status of the investigations of the knee region of the cosmic ray spectrum by KASCADE is presented and open problems are discussed.
Two approaches to the cosmic ray "knee" problem and the experimental data

Yuri V. Stenkin

DOI: https://doi.org/10.48550/arXiv.astro-ph/0507018

2005-07-01

Astrophysics

Abstract:Two approaches to the old problem of the so-called "knee" in cosmic ray spectrum are compared with the existing experimental data. It is shown that "standard astrophysical model" of the knee is not supported by some experimental data, while the alternative model agrees better with the data.
Statistics for Searches at the LHC

Glen Cowan

DOI: https://doi.org/10.48550/arXiv.1307.2487

2013-07-09

Abstract:These lectures describe several topics in statistical data analysis as used in High Energy Physics. They focus on areas most relevant to analyses at the LHC that search for new physical phenomena, including statistical tests for discovery and exclusion limits. Particular attention is payed to the treatment of systematic uncertainties through nuisance parameters.

High Energy Physics - Experiment,Data Analysis, Statistics and Probability
Cosmic Ray Transport in the Galaxy: a Review

Elena Amato,Pasquale Blasi

DOI: https://doi.org/10.1016/j.asr.2017.04.019

2017-04-19

Abstract:The physics of energetic particle propagation in magnetised environments plays a crucial role in both the processes of acceleration and transport of cosmic rays. Recent theoretical developments in the field of cosmic ray research have been mainly in the direction of exploring non-linear aspects of the processes in which these particles are involved, namely the action of cosmic rays on the environment in which the transport and/or acceleration take place. When cosmic rays propagate outside of the acceleration region, such action is mainly in two forms: 1) they generate hydromagnetic waves, through streaming instabilities, leading to a dependence of the scattering properties of the medium on the spectrum and spatial distribution of the energetic particles, and 2) they exert a dynamical action on the plasma, which may cause the launching of cosmic ray driven Galactic winds. In this article we discuss these and other recent developments and how they compare with the bulk of new observations on the spectra of primary nuclei (mainly H and He) and secondary to primary ratios, such as the B/C ratio and the $\bar p$/ratio, and the positrons ratio $e^+/(e^-+e^+)$. We also comment on some radically new models of the origin of CRs, in which the physical meaning of the secondary to primary ratios is not the same as in the standard model.

High Energy Astrophysical Phenomena
Bayesian Computation in Astronomy: Novel methods for parallel and gradient-free inference

Minas Karamanis

DOI: https://doi.org/10.7488/era/2932

2023-03-29

Abstract:The goal of this thesis is twofold; introduce the fundamentals of Bayesian inference and computation focusing on astronomical and cosmological applications, and present recent advances in probabilistic computational methods developed by the author that aim to facilitate Bayesian data analysis for the next generation of astronomical observations and theoretical models. The first part of this thesis familiarises the reader with the notion of probability and its relevance for science through the prism of Bayesian reasoning, by introducing the key constituents of the theory and discussing its best practices. The second part includes a pedagogical introduction to the principles of Bayesian computation motivated by the geometric characteristics of probability distributions and followed by a detailed exposition of various methods including Markov chain Monte Carlo (MCMC), Sequential Monte Carlo (SMC), and Nested Sampling (NS). Finally, the third part presents two novel computational methods (Ensemble Slice Sampling and Preconditioned Monte Carlo) and their respective software implementations (zeus and pocoMC). [abridged]

Instrumentation and Methods for Astrophysics,Cosmology and Nongalactic Astrophysics,Earth and Planetary Astrophysics

Analysis of cosmic ray data

Statistical Methods for Investigating the Cosmic Ray Energy Spectrum

Statistical Tools for Analyzing the Cosmic Ray Energy Spectrum

Techniques for Searching for Anisotropy in the Arrival Directions of Ultra-High Energy Cosmic Rays

II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements

Data analysis methods for the cosmic microwave background

Novel Simulation Framework for Analyzing Cosmic Ray Particle Distributions at a Global Scale

Analysis Methods for Gamma-ray Astronomy

Reanalysis of the Systematic Uncertainties in Cosmic-Ray Antiproton Flux

Reanalysis of the systematic uncertainties in the cosmic-ray antiproton flux

Analysis and packaging of radiochemical solar neutrino data. 1. Bayesian approach

Statistical Issues in Neutrino Physics Analyses

Data-driven analysis for understanding ultrahigh energy cosmic ray source spectra

Stochastic modelling of cosmic ray sources for diffuse high-energy gamma-rays and neutrinos

Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy

A novel method of data analysis for hadronic physics

What Do We Expect to Learn from Very High Energy Cosmic Ray Observations?

Two approaches to the cosmic ray "knee" problem and the experimental data

Statistics for Searches at the LHC

Cosmic Ray Transport in the Galaxy: a Review

Bayesian Computation in Astronomy: Novel methods for parallel and gradient-free inference