A statistical procedure for the identification of positrons in the PAMELA experiment

The PAMELA satellite experiment has measured the cosmic-ray positron fraction between 1.5 GeV and 100 GeV. The need to reliably discriminate between the positron signal and proton background has required the development of an ad hoc analysis procedure. In this paper, a method for positron identification is described and its stability and capability to yield a correct background estimate is shown. The analysis includes new experimental data, the application of three different fitting techniques for the background sample and an estimate of systematic uncertainties due to possible inaccuracies in the background selection. The new experimental results confirm both solar modulation effects on cosmic-rays with low rigidities and an anomalous positron abundance above 10 GeV.

💡 Research Summary

The paper presents a comprehensive statistical methodology developed for the PAMELA satellite experiment to isolate the cosmic‑ray positron signal from the overwhelming proton background across the energy range of 1.5 GeV to 100 GeV. PAMELA, launched in 2006, combines a magnetic spectrometer, time‑of‑flight system, silicon‑tungsten calorimeter, and anticoincidence shields to simultaneously measure electrons, positrons, protons, and light nuclei. Because protons outnumber positrons by three orders of magnitude, even a tiny mis‑identification probability can dominate the measured positron fraction. The authors therefore designed a multi‑step selection procedure that first exploits charge‑sign discrimination (via curvature in the magnetic field), energy‑loss patterns in the silicon layers, and shower topology in the calorimeter to define a clean candidate sample.

To estimate the residual proton contamination, three independent fitting techniques are applied to a dedicated proton control sample selected with inverted calorimetric cuts. (1) A conventional Gaussian fit captures the core of the proton distribution in the discriminant variable. (2) A polynomial‑based non‑linear fit accounts for energy‑dependent asymmetries and long tails. (3) Kernel density estimation (KDE) provides a non‑parametric reconstruction of the full shape, preserving subtle features that parametric models may miss. Each method is evaluated using χ² per degree of freedom, the Akaike Information Criterion, and cross‑validation residuals. The final background model is obtained by averaging the three fits and propagating their spread as a systematic component.

Systematic uncertainties are quantified in two complementary ways. First, the selection thresholds (e.g., calorimeter energy‑deposit cuts, track‑length requirements) are varied by ±5 % to gauge the sensitivity of the background estimate to the definition of the control sample. Second, a detailed GEANT4 simulation of the PAMELA instrument is used to model detection efficiencies, energy‑reconstruction biases, and track‑reconstruction errors. By comparing simulated and real data, the authors constrain the overall systematic error on the positron fraction to be between 2.5 % and 4 %, which is then combined in quadrature with the statistical uncertainties.

The resulting positron fraction shows two distinct regimes. Below ~10 GeV the fraction exhibits a clear solar‑modulation pattern, decreasing during periods of high solar activity, consistent with earlier balloon‑borne and ground‑based measurements. Above ~10 GeV the positron fraction rises sharply, exceeding the predictions of conventional secondary‑production models that attribute positrons to spallation of primary cosmic‑ray protons on the interstellar medium. This excess persists up to the highest measured energy (≈100 GeV) and mirrors the “positron anomaly” reported by other experiments such as AMS‑02 and Fermi‑LAT. The authors discuss possible interpretations, including dark‑matter particle annihilation or decay, nearby pulsar wind nebulae, and supernova‑remnant acceleration of secondary particles.

Importantly, the paper demonstrates that the adopted statistical framework yields a stable and unbiased background estimate, even when the control sample is perturbed within realistic limits. The three‑fold fitting strategy, together with rigorous systematic studies, provides a template for future high‑precision cosmic‑ray experiments where signal‑to‑background ratios are extremely low. The authors conclude that the observed high‑energy positron excess is robust against methodological variations, reinforcing the case for new astrophysical or particle‑physics phenomena contributing to the cosmic‑ray lepton spectrum.