Sifting data in the real world

Abstract

In the real world, experimental data are rarely, if ever, distributed as a normal (Gaussian) distribution. As an example, a large set of data--such as the cross sections for particle scattering as a function of energy contained in the archives of the Particle Data Group--is a compendium of all published data, and hence, unscreened. Inspection of similar data sets quickly shows that, for many reasons, these data sets have many outliers--points well beyond what is expected from a normal distribution--thus ruling out the use of conventional ${\chi }^2$ techniques. This note suggests an adaptive algorithm that allows a phenomenologist to apply to the data sample a sieve whose mesh is coarse enough to let the background fall through, but fine enough to retain the preponderance of the signal, thus sifting the data. A prescription is given for finding a robust estimate of the best-fit model parameters in the presence of a noisy background, together with a robust estimate of the model parameter errors, as well as a determination of the goodness-of-fit of the data to the theoretical hypothesis. Extensive computer simulations are carried out to test the algorithm for both its accuracy and stability under varying background conditions.

29 pages, 13 figures. Version to appear in Nucl. Instr. & Meth. A

Cited by in corpus (21)

• Hadronic forward scattering: Predictions for the Large Hadron Collider and cosmic rays
• Connection of the virtual ${\gamma }^{\ast }p$ cross section of $ep$ deep inelastic scattering to real $\gamma p$ scattering, and the implications for $\nu N$ and $ep$ total cross sections
• String-Gauge Dual Description of Deep Inelastic Scattering at Small-$x$
• Analytic Expression for the Joint x and Q^2 Dependences of the Structure Functions of Deep Inelastic Scattering
• Applications of the leading-order Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations to the combined HERA data on deep inelastic scattering
• Ultrahigh energy neutrino scattering: an update
• A new numerical method for obtaining gluon distribution functions $G(x,Q^2)=xg(x,Q^2)$, from the proton structure function $F_2^{\gamma p}(x,Q^2)$
• Implications of a Froissart bound saturation of ${\gamma }^{\ast }$-$p$ deep inelastic scattering. Part I. Quark distributions at ultra small $x$
• Analyticity as a Robust Constraint on the LHC Cross Section
• Experimental consistency in parton distribution fitting
• Comprehensive fits to high energy data for $\sigma$, $\rho$, and $B$ and the asymptotic black-disk limit
• Ultra-high Energy Predictions of Proton-Air Cross Sections from Accelerator Data: an Update
• Ultra-high Energy Predictions of proton-air Cross Sections from Accelerator Data
• Test of Universal Rise of Hadronic Total Cross Sections based on pi p, Kp and pbar p,pp Scatterings
• Scope and limitations of a string theory dual description of the proton structure
• The slope, curvature, and higher parameters in $pp$ and $\bar{p}p$ scattering, and the extrapolation of measurements of $d\sigma (s,t)/dt$ to $t=0$
• Unitarisation dependence of diffractive scattering in light of high-energy collider data
• New limits on "odderon" amplitudes from analyticity constraints
• Eikonal and asymptotic fits to high energy data for $\sigma$, $\rho$, and $B$: An update with curvature corrections
• The proton inelastic cross section at ultrahigh energies
• Implications from analyticity constraints used in a Landshoff-Donnachie fit