On the challenge of estimating diphoton backgrounds at large invariant mass

We examine, using the analyses of the 750 GeV diphoton resonance as a case study, the methodology for estimating the dominant backgrounds to diphoton resonance searches. We show that close to the high energy tails of the distributions, where background estimates rely on functional extrapolations or Monte Carlo predictions, large uncertainties are introduced, in particular by the challenging photonjet background. Analyses with loose photon and low photon pT cuts and those susceptible to high photon rapidity regions are especially affected. Given that diphoton-based searches beyond 1TeV are highly motivated as discovery modes, these considerations are relevant for future analyses. We first consider a physics-driven deformation of the photonjet spectrum by next-to-leading order effects and a phase space dependent fake rate and show that this reduces the local significance of the excess. Using a simple but more general ansatz, we demonstrate that the originally reported local significances of the 750 GeV excess could have been overestimated by more than one standard deviation. We furthermore cross-check our analysis by comparing fit results based on the 2015 and 2016 LHC data sets. Finally we employ our methodology on the available 13 TeV LHC data set assessing the systematics involved in the current diphoton searches beyond the TeV region.