Open Access
September 2015 Statistical unfolding of elementary particle spectra: Empirical Bayes estimation and bias-corrected uncertainty quantification
Mikael Kuusela, Victor M. Panaretos
Ann. Appl. Stat. 9(3): 1671-1705 (September 2015). DOI: 10.1214/15-AOAS857

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.

Citation

Download Citation

Mikael Kuusela. Victor M. Panaretos. "Statistical unfolding of elementary particle spectra: Empirical Bayes estimation and bias-corrected uncertainty quantification." Ann. Appl. Stat. 9 (3) 1671 - 1705, September 2015. https://doi.org/10.1214/15-AOAS857

Information

Received: 1 January 2014; Revised: 1 July 2015; Published: September 2015
First available in Project Euclid: 2 November 2015

zbMATH: 06526003
MathSciNet: MR3418740
Digital Object Identifier: 10.1214/15-AOAS857

Keywords: bootstrap , high energy physics , Large Hadron Collider , Monte Carlo EM Algorithm , Poisson inverse problem , Poisson process , regularization

Rights: Copyright © 2015 Institute of Mathematical Statistics

Vol.9 • No. 3 • September 2015
Back to Top