Confidence Level Computation for Combining Searches with Small Statistics
Thomas Junk
TL;DR
The paper presents an efficient, approximate method to compute confidence levels for excluding new particles when multiple small-statistics searches are combined. It uses a likelihood-ratio-based test statistic that factorizes across channels, with $CL_{s+b}$, $CL_b$, and $CL_s$ defined in the modified frequentist framework. By constructing and binning the PDF of the test statistic and convolving across channels, the approach avoids exhaustive enumeration and can incorporate uncorrelated systematic uncertainties via Gaussian averaging. The method is validated against exact calculations and Monte Carlo benchmarks, and is demonstrated on Higgs/MSSM LEP analyses and mock experiments, offering a practical tool for rapid model-scanning and multi-channel combinations.
Abstract
This article describes an efficient procedure for computing approximate confidence levels for searches for new particles where the expected signal and background levels are small enough to require the use of Poisson statistics. The results of many independent searches for the same particle may be combined easily, regardless of the discriminating variables which may be measured for the candidate events. The effects of systematic uncertainty in the signal and background models are incorporated in the confidence levels. The procedure described allows efficient computation of expected confidence levels.
