-
Simple and statistically sound recommendations for analysing physical theories
Authors:
Shehu S. AbdusSalam,
Fruzsina J. Agocs,
Benjamin C. Allanach,
Peter Athron,
Csaba Balázs,
Emanuele Bagnaschi,
Philip Bechtle,
Oliver Buchmueller,
Ankit Beniwal,
Jihyun Bhom,
Sanjay Bloor,
Torsten Bringmann,
Andy Buckley,
Anja Butter,
José Eliel Camargo-Molina,
Marcin Chrzaszcz,
Jan Conrad,
Jonathan M. Cornell,
Matthias Danninger,
Jorge de Blas,
Albert De Roeck,
Klaus Desch,
Matthew Dolan,
Herbert Dreiner,
Otto Eberhardt
, et al. (50 additional authors not shown)
Abstract:
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 mul…
▽ More
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 https://doi.org/10.5281/zenodo.4322283.
△ Less
Submitted 11 April, 2022; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Global Analysis of Dark Matter Simplified Models with Leptophobic Spin-One Mediators using MasterCode
Authors:
E. Bagnaschi,
J. C. Costa,
K. Sakurai,
M. Borsato,
O. Buchmueller,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flächer,
K. Hahn,
S. Heinemeyer,
M. Lucio,
D. Martínez Santos,
K. A. Olive,
S. Trifa,
G. Weiglein
Abstract:
We report the results of a global analysis of dark matter simplified models (DMSMs) with leptophobic mediator particles of spin one, considering the cases of both vector and axial-vector interactions with dark matter (DM) particles and quarks. We require the DMSMs to provide all the cosmological DM density indicated by Planck and other observations, and we impose the upper limits on spin-independe…
▽ More
We report the results of a global analysis of dark matter simplified models (DMSMs) with leptophobic mediator particles of spin one, considering the cases of both vector and axial-vector interactions with dark matter (DM) particles and quarks. We require the DMSMs to provide all the cosmological DM density indicated by Planck and other observations, and we impose the upper limits on spin-independent and -dependent scattering from direct DM search experiments. We also impose all relevant LHC constraints from searches for monojet events and measurements of the dijet mass spectrum. We model the likelihood functions for all the constraints and combine them within the MasterCode framework, and probe the full DMSM parameter spaces by scanning over the mediator and DM masses and couplings, not fixing any of the model parameters. We find, in general, two allowed regions of the parameter spaces: one in which the mediator couplings to Standard Model (SM) and DM particles may be comparable to those in the SM and the cosmological DM density is reached via resonant annihilation, and one in which the mediator couplings to quarks are $\lesssim 10^{-3}$ and DM annihilation is non-resonant. We find that the DM and mediator masses may well lie within the ranges accessible to LHC experiments. We also present predictions for spin-independent and -dependent DM scattering, and present specific results for ranges of the DM couplings that may be favoured in ultraviolet completions of the DMSMs.
△ Less
Submitted 2 May, 2019;
originally announced May 2019.
-
Likelihood Analysis of the Sub-GUT MSSM in Light of LHC 13-TeV Data
Authors:
J. C. Costa,
E. Bagnaschi,
K. Sakurai,
M. Borsato,
O. Buchmueller,
M. Citron,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flächer,
S. Heinemeyer,
M. Lucio,
D. Martínez Santos,
K. A. Olive,
A. Richards,
G. Weiglein
Abstract:
We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gau…
▽ More
We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gaugino mass $m_{1/2}$, a common soft supersymmetry-breaking scalar mass $m_0$, a common trilinear mixing parameter $A$ and the ratio of MSSM Higgs vevs $\tanβ$, assuming that the Higgs mixing parameter $μ> 0$. We take into account constraints on strongly- and electroweakly-interacting sparticles from $\sim 36$/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for $M_{in} \sim 10^5$ to $10^9$ GeV, with $M_{in} \sim M_{GUT}$ disfavoured by $Δχ^2 \sim 3$ due to the ${\rm BR}(B_{s, d} \to μ^+μ^-)$ constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with $\tilde{χ^0_1} \sim 1$ TeV, with annihilation via heavy Higgs bosons $H/A$ and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the $(g-2)_μ$ constraint, finding similar results whether or not it is included.
△ Less
Submitted 15 March, 2018; v1 submitted 1 November, 2017;
originally announced November 2017.
-
Likelihood Analysis of the pMSSM11 in Light of LHC 13-TeV Data
Authors:
E. Bagnaschi,
K. Sakurai,
M. Borsato,
O. Buchmueller,
M. Citron,
J. C. Costa,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flächer,
S. Heinemeyer,
M. Lucio,
D. Martínez Santos,
K. A. Olive,
A. Richards,
V. C. Spanos,
I. Suárez Fernández,
G. Weiglein
Abstract:
We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from ~ 36/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the $(g-2)_μ$ constraint. The pMSSM…
▽ More
We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from ~ 36/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the $(g-2)_μ$ constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses $M_{1,2,3}$, a common mass for the first-and second-generation squarks $m_{\tilde{q}}$ and a distinct third-generation squark mass $m_{\tilde{q}_3}$, a common mass for the first-and second-generation sleptons $m_{\tilde l}$ and a distinct third-generation slepton mass $m_{\tilde τ}$, a common trilinear mixing parameter $A$, the Higgs mixing parameter $μ$, the pseudoscalar Higgs mass $M_A$ and $\tanβ$. In the fit including $(g-2)_μ$, a Bino-like $\tildeχ^0_1$ is preferred, whereas a Higgsino-like $\tilde χ^0_1$ is favoured when the $(g-2)_μ$ constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, $\tildeχ^0_1$, into the range indicated by cosmological data. In the fit including $(g-2)_μ$, coannihilations with $\tilde χ^0_2$ and the Wino-like $\tildeχ^{\pm}_1$ or with nearly-degenerate first- and second-generation sleptons are favoured, whereas coannihilations with the $\tilde χ^0_2$ and the Higgsino-like $\tildeχ^{\pm}_1$ or with first- and second-generation squarks may be important when the $(g-2)_μ$ constraint is dropped. Prospects remain for discovering strongly-interacting sparticles at the LHC as well as for discovering electroweakly-interacting sparticles at a future linear $e^+ e^-$ collider such as the ILC or CLIC.
△ Less
Submitted 1 May, 2018; v1 submitted 30 October, 2017;
originally announced October 2017.
-
Likelihood Analysis of the Minimal AMSB Model
Authors:
E. Bagnaschi,
M. Borsato,
K. Sakurai,
O. Buchmueller,
R. Cavanaugh,
V. Chobanova,
M. Citron,
J. C. Costa,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flächer,
S. Heinemeyer,
G. Isidori,
M. Lucio,
F. Luo,
D. Martínez Santos,
K. A. Olive,
A. Richards,
G. Weiglein
Abstract:
We perform a likelihood analysis of the minimal Anomaly-Mediated Supersymmetry Breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that a wino-like or a Higgsino-like neutralino LSP, $m_{\tilde χ^0_{1}}$, may provide the cold dark matter (DM) with similar likelihood. The upper limit on the DM density from Planck and other experiments enforces…
▽ More
We perform a likelihood analysis of the minimal Anomaly-Mediated Supersymmetry Breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that a wino-like or a Higgsino-like neutralino LSP, $m_{\tilde χ^0_{1}}$, may provide the cold dark matter (DM) with similar likelihood. The upper limit on the DM density from Planck and other experiments enforces $m_{\tilde χ^0_{1}} \lesssim 3~TeV$ after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the $\tilde χ_0^1$, the measured value of the Higgs mass favours a limited range of $\tan β\sim 5$ (or for $μ> 0$, $\tan β\sim 45$) but the scalar mass $m_0$ is poorly constrained. In the wino-LSP case, $m_{3/2}$ is constrained to about $900~TeV$ and ${m_{\tilde χ^0_{1}}}$ to $2.9\pm0.1~TeV$, whereas in the Higgsino-LSP case $m_{3/2}$ has just a lower limit $\gtrsim 650TeV$ ($\gtrsim 480TeV$) and $m_{\tilde χ^0_{1}}$ is constrained to $1.12 ~(1.13) \pm0.02~TeV$ in the $μ>0$ ($μ<0$) scenario. In neither case can the anomalous magnetic moment of the muon, ${(g-2)_μ}$, be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, {though there} are some prospects for direct DM detection. On the other hand, if the ${m_{\tilde χ^0_{1}}}$ contributes only a fraction of the cold DM density, {future LHC $E_T$-based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant}, and interference effects enable ${\rm BR}(B_{s, d} \to μ^+μ^-)$ to agree with the data better than in the SM in the case of wino-like DM with $μ> 0$.
△ Less
Submitted 15 December, 2016;
originally announced December 2016.
-
Likelihood Analysis of Supersymmetric SU(5) GUTs
Authors:
E. Bagnaschi,
J. C. Costa,
K. Sakurai,
M. Borsato,
O. Buchmueller,
R. Cavanaugh,
V. Chobanova,
M. Citron,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flächer,
S. Heinemeyer,
G. Isidori,
M. Lucio,
D. Martínez Santos,
K. A. Olive,
A. Richards,
K. J. de Vries,
G. Weiglein
Abstract:
We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has 7 parameters: a universal gaugino mass $m_{1/2}$, distinct masses for the scalar partners of matter fermions in five- and ten-di…
▽ More
We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has 7 parameters: a universal gaugino mass $m_{1/2}$, distinct masses for the scalar partners of matter fermions in five- and ten-dimensional representations of SU(5), $m_5$ and $m_{10}$, and for the $\mathbf{5}$ and $\mathbf{\bar 5}$ Higgs representations $m_{H_u}$ and $m_{H_d}$, a universal trilinear soft SUSY-breaking parameter $A_0$, and the ratio of Higgs vevs $\tan β$. In addition to previous constraints from direct sparticle searches, low-energy and flavour observables, we incorporate constraints based on preliminary results from 13 TeV LHC searches for jets + MET events and long-lived particles, as well as the latest PandaX-II and LUX searches for direct Dark Matter detection. In addition to previously-identified mechanisms for bringing the supersymmetric relic density into the range allowed by cosmology, we identify a novel ${\tilde u_R}/{\tilde c_R} - \tildeχ^0_1$ coannihilation mechanism that appears in the supersymmetric SU(5) GUT model and discuss the role of ${\tilde ν_τ}$ coannihilation. We find complementarity between the prospects for direct Dark Matter detection and SUSY searches at the LHC.
△ Less
Submitted 26 April, 2017; v1 submitted 31 October, 2016;
originally announced October 2016.
-
Supersymmetric Dark Matter after LHC Run 1
Authors:
E. A. Bagnaschi,
O. Buchmueller,
R. Cavanaugh,
M. Citron,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flaecher,
S. Heinemeyer,
G. Isidori,
S. Malik,
D. Martinez Santos,
K. A. Olive,
K. Sakurai,
K. J. de Vries,
G. Weiglein
Abstract:
Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, neutralino_1, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate next-to-lightest supersymmetric particle (NLSP) such as the lighte…
▽ More
Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, neutralino_1, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate next-to-lightest supersymmetric particle (NLSP) such as the lighter stau (stau_1), stop (stop_1) or chargino (chargino_1), resonant annihilation via direct-channel heavy Higgs bosons H/A, the light Higgs boson h or the Z boson, and enhanced annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2 and pMSSM10. We analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM mechanisms within these scenarios. We find that the stau_1 coannihilation regions of the CMSSM, NUHM1, NUHM2 can largely be explored at the LHC via searches for missing E_T events and long-lived charged particles, whereas their H/A funnel, focus-point and chargino_1 coannihilation regions can largely be explored by the LZ and Darwin DM direct detection experiments. We find that the dominant DM mechanism in our pMSSM10 analysis is chargino_1 coannihilation: {parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches.
△ Less
Submitted 5 August, 2015;
originally announced August 2015.
-
The pMSSM10 after LHC Run 1
Authors:
K. J. de Vries,
E. A. Bagnaschi,
O. Buchmueller,
R. Cavanaugh,
M. Citron,
A. De Roeck,
M. J. Dolan,
J. R. Ellis,
H. Flaecher,
S. Heinemeyer,
G. Isidori,
S. Malik,
J. Marrouche,
D. Martinez Santos,
K. A. Olive,
K. Sakurai,
G. Weiglein
Abstract:
We present a frequentist analysis of the parameter space of the pMSSM10, in which the following 10 soft SUSY-breaking parameters are specified independently at the mean scalar top mass scale Msusy = Sqrt[M_stop1 M_stop2]: the gaugino masses M_{1,2,3}, the 1st-and 2nd-generation squark masses M_squ1 = M_squ2, the third-generation squark mass M_squ3, a common slepton mass M_slep and a common triline…
▽ More
We present a frequentist analysis of the parameter space of the pMSSM10, in which the following 10 soft SUSY-breaking parameters are specified independently at the mean scalar top mass scale Msusy = Sqrt[M_stop1 M_stop2]: the gaugino masses M_{1,2,3}, the 1st-and 2nd-generation squark masses M_squ1 = M_squ2, the third-generation squark mass M_squ3, a common slepton mass M_slep and a common trilinear mixing parameter A, the Higgs mixing parameter mu, the pseudoscalar Higgs mass M_A and tan beta. We use the MultiNest sampling algorithm with 1.2 x 10^9 points to sample the pMSSM10 parameter space. A dedicated study shows that the sensitivities to strongly-interacting SUSY masses of ATLAS and CMS searches for jets, leptons + MET signals depend only weakly on many of the other pMSSM10 parameters. With the aid of the Atom and Scorpion codes, we also implement the LHC searches for EW-interacting sparticles and light stops, so as to confront the pMSSM10 parameter space with all relevant SUSY searches. In addition, our analysis includes Higgs mass and rate measurements using the HiggsSignals code, SUSY Higgs exclusion bounds, the measurements B-physics observables, EW precision observables, the CDM density and searches for spin-independent DM scattering. We show that the pMSSM10 is able to provide a SUSY interpretation of (g-2)_mu, unlike the CMSSM, NUHM1 and NUHM2. As a result, we find (omitting Higgs rates) that the minimum chi^2/dof = 20.5/18 in the pMSSM10, corresponding to a chi^2 probability of 30.8 %, to be compared with chi^2/dof = 32.8/24 (31.1/23) (30.3/22) in the CMSSM (NUHM1) (NUHM2). We display 1-dimensional likelihood functions for SUSY masses, and show that they may be significantly lighter in the pMSSM10 than in the CMSSM, NUHM1 and NUHM2. We discuss the discovery potential of future LHC runs, e+e- colliders and direct detection experiments.
△ Less
Submitted 13 April, 2015;
originally announced April 2015.
