Simple and statistically sound recommendations for analysing physical theories.

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at Zenodo.

Jet cross sections at the LHC and the quest for higher precision.

We perform a phenomenological study of Z plus jet, Higgs plus jet and di-jet production at the Large Hadron Collider. We investigate in particular the dependence of the leading jet cross section on the jet radius as a function of the jet transverse momentum. Theoretical predictions are obtained using perturbative QCD calculations at the next-to and next-to-next-to-leading order, using a range of renormalization and factorization scales. The fixed order predictions are compared to results obtained from matching next-to-leading order calculations to parton showers. A study of the scale dependence as a function of the jet radius is used to provide a better estimate of the scale uncertainty for small jet sizes. The non-perturbative corrections as a function of jet radius are estimated from different generators.

Status of the scalar singlet dark matter model.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a Z 2 symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above ∼ 1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.