Statistical unfolding of elementary particle spectra: Empirical Bayes estimation and bias-corrected uncertainty quantification

Abstract

We consider the high energy physics unfolding problem where the goal is to estimate the spectrum of elementary particles given observations distorted by the limited resolution of a particle detector. This important statistical inverse problem arising in data analysis at the Large Hadron Collider at CERN consists in estimating the intensity function of an indirectly observed Poisson point process. Unfolding typically proceeds in two steps: one first produces a regularized point estimate of the unknown intensity and then uses the variability of this estimator to form frequentist confidence intervals that quantify the uncertainty of the solution. In this paper, we propose forming the point estimate using empirical Bayes estimation which enables a data-driven choice of the regularization strength through marginal maximum likelihood estimation. Observing that neither Bayesian credible intervals nor standard bootstrap confidence intervals succeed in achieving good frequentist coverage in this problem due to the inherent bias of the regularized point estimate, we introduce an iteratively bias-corrected bootstrap technique for constructing improved confidence intervals. We show using simulations that this enables us to achieve nearly nominal frequentist coverage with only a modest increase in interval length. The proposed methodology is applied to unfolding the $Z$ boson invariant mass spectrum as measured in the CMS experiment at the Large Hadron Collider.

Published at http://dx.doi.org/10.1214/15-AOAS857 in the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org). arXiv admin note: substantial text overlap with arXiv:1401.8274

References in corpus (3)

• Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC
• Measurement of the $Z/{\gamma }^{\ast }$ boson transverse momentum distribution in $pp$ collisions at $\sqrt{s}$ = 7 TeV with the ATLAS detector
• Statistical unfolding of elementary particle spectra: Empirical Bayes estimation and bias-corrected uncertainty quantification

Cited by in corpus (12)

• Simultaneous measurement of the muon neutrino charged-current cross section on oxygen and carbon without pions in the final state at T2K
• Statistical unfolding of elementary particle spectra: Empirical Bayes estimation and bias-corrected uncertainty quantification
• A Hierarchical Bayesian Approach to Neutron Spectrum Unfolding with Organic Scintillators
• Uncertainty quantification for wide-bin unfolding: one-at-a-time strict bounds and prior-optimized confidence intervals
• A new method for obtaining a Born cross section using visible cross section data from $e^{+}{e}^{-}$ colliders
• Model independent measurements of Standard Model cross sections with Domain Adaptation
• Convergence and error propagation results on a linear iterative unfolding method
• Uniform Deconvolution for Poisson Point Processes
• Objective frequentist uncertainty quantification for atmospheric CO${}_2$ retrievals
• Deconvolution density estimation with penalised MLE
• Frequentist-Bayes Hybrid Covariance Estimationfor Unfolding Problems
• A data-driven convergence criterion for iterative unfolding of smeared spectra