The CMSSM and NUHM1 after LHC Run 1

We analyze the impact of data from the full Run 1 of the LHC at 7 and 8 TeV on the CMSSM with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu > 0$$\end{document}μ>0 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0$$\end{document}<0 and the NUHM1 with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu > 0$$\end{document}μ>0, incorporating the constraints imposed by other experiments such as precision electroweak measurements, flavour measurements, the cosmological density of cold dark matter and the direct search for the scattering of dark matter particles in the LUX experiment. We use the following results from the LHC experiments: ATLAS searches for events with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E\!\!/}_{T}$$\end{document}E/T accompanied by jets with the full 7 and 8 TeV data, the ATLAS and CMS measurements of the mass of the Higgs boson, the CMS searches for heavy neutral Higgs bosons and a combination of the LHCb and CMS measurements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{BR}(B_s \rightarrow \mu ^+\mu ^-)$$\end{document}BR(Bs→μ+μ-) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{BR}(B_d \rightarrow \mu ^+\mu ^-)$$\end{document}BR(Bd→μ+μ-). Our results are based on samplings of the parameter spaces of the CMSSM for both \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu >0$$\end{document}μ>0 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu <0$$\end{document}μ<0 and of the NUHM1 for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu > 0$$\end{document}μ>0 with 6.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^6$$\end{document}×106, 6.2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^6$$\end{document}×106 and 1.6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^7$$\end{document}×107 points, respectively, obtained using the MultiNest tool. The impact of the Higgs-mass constraint is assessed using FeynHiggs 2.10.0, which provides an improved prediction for the masses of the MSSM Higgs bosons in the region of heavy squark masses. It yields in general larger values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_h$$\end{document}Mh than previous versions of FeynHiggs, reducing the pressure on the CMSSM and NUHM1. We find that the global \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2$$\end{document}χ2 functions for the supersymmetric models vary slowly over most of the parameter spaces allowed by the Higgs-mass and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E\!\!/}_{T}$$\end{document}E/T searches, with best-fit values that are comparable to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/\mathrm{dof}$$\end{document}χ2/dof for the best Standard Model fit. We provide 95 % CL lower limits on the masses of various sparticles and assess the prospects for observing them during Run 2 of the LHC.


Introduction
In addition to establishing the mechanism for electroweak symmetry breaking, one of the primary objectives of experiments at the LHC has been to search for possible physics beyond the Standard Model (SM), such as new particles that might alleviate the naturalness problem and/or be associated with cosmological dark matter. In contrast with the triumphant discovery at the LHC of a particle that resembles the Higgs boson of the SM [1,2], and the observation of BR(B s → μ + μ − ) decay at a rate close to the SM prediction [3,4], the first run of the LHC has not revealed any convincing evidence for physics beyond the SM. In particular, the LHC searches for jets + E / T events [5,6] and for heavy Higgs bosons H ± /H/A [7] have drawn blanks so far. In parallel, neither direct nor indirect searches for astrophysical dark matter have found any convincing signals [8,9], posing questions regarding the implications of those results for supersymmetric models.
Prior to the LHC start-up, the discrepancy between the experimental measurement of (g − 2) μ [34,35] and theoretical calculations (see [36,37] and references therein), favoured relatively light sparticle masses, but these have not appeared in ATLAS and CMS E / T searches, disfavouring small values of the CMSSM or NUHM1 SUSY-breaking mass parameters [5]. On the other hand, the discovery of a SM-like Higgs boson by ATLAS and CMS [1,2] with a mass consistent with the predictions of SUSY models has provided an important indirect constraint on SUSY model parameters such as m 0 , m 1/2 , A 0 and tan β. A significant rôle is also played by the observation by CMS and LHCb of BR(B s → μ + μ − ) decay [3,4], which imposes a complementary constraint on the CMSSM and NUHM1 parameter spaces. Our most recent analyses [33] of these models were based on the 7-TeV ATLAS 5/fb jets + E / T data set [38,39]. In this paper we update our analyses to include the 8-TeV ATLAS 20/fb jets + E / T data set [5], providing a complete study of the implications of the LHC Run 1 for the CMSSM and NUHM1 scenarios. On the basis of this study, we also discuss the prospects for discovering sparticles in the LHC Run 2.
As described below, these constraints are analyzed in a frequentist approach using an overhauled version of the MasterCode [40] framework to calculate the global χ 2 function. For other recent post-LHC analyses of the CMSSM and NUHM, see . In this paper we sample the CMSSM and NUHM1 parameter spaces using the MultiNest tool [67,68], which is more efficient than the Markov Chain Monte Carlo technique we used previously. We implement the ATLAS 20/fb jets + E / T constraint using scaling laws to extrapolate the sensitivity to regions of the parameter space where documentation is not available [33]. In our implementation of the M h constraint we use a new version of FeynHiggs, 2.10.0 [69], which incorporates a resummation of leading and subleading logarithmic corrections from the top/stop sector to provide improved results for larger stop masses. Since in the SUSY models we study BR(B s → μ + μ − ) and BR(B d → μ + μ − ) are expected to have the same ratio as in the SM, we combine these con-straints by assuming this SM ratio and incorporating the experimental correlations between BR(B s → μ + μ − ) and BR(B d → μ + μ − ) reported by the LHCb and CMS Collaborations [3,4]. Recent cosmological observations, including those by the Planck satellite [70], have refined the estimate of the cold dark matter density, but this does not have a relevant impact on our study. Concerning direct searches for dark matter, the only constraint we apply is that from LUX [9] on spin-independent dark matter scattering, which we incorporate taking due account of the uncertainties in the hadronic scattering matrix element, as discussed later.
We find that the global χ 2 function varies relatively little across most of the regions of the (m 0 , m 1/2 ) planes that are allowed by the E / T , Higgs and dark matter density constraints on the CMSSM and NUHM1 parameter spaces, with a global minimum at large m 0 and m 1/2 , which is similar to the χ 2 /dof for the SM. Within the CMSSM, there are four principal mechanisms for bringing the SUSY relic density χ h 2 into the range favoured by Planck and other measurements [70], namely coannihilation with the lighter stauτ 1 and other sleptons, coannihilation with the lighter stop t 1 , rapid annihilation through the heavy Higgs bosons H, A in the direct channel and annihilation in the focus-point region where the lightest neutralinoχ 0 1 has an enhanced Higgsino component. In the following, we comment on the respective rôles of these mechanisms. 1 Within the range of the CMSSM parameter space examined in this paper, the χ h 2 constraint sets an upper bound on m 1/2 but not on m 0 . In the case of the NUHM1, more annihilation mechanisms may come into play, and we find no upper bound on either m 0 or m 1/2 .
One of our key findings is that the LHC measurement of M h is not in tension with other constraints on the CMSSM and NUHM1 parameter spaces except for (g − 2) μ . The M h constraint does not impact them as strongly as had previously been thought [62,63], since the improved prediction incorporated in FeynHiggs 2.10.0 [69] yields in general a higher value of M h than previous versions of FeynHiggs (as well as SoftSUSY) for the same values of the model parameters [64], as will be discussed in detail in Sect. 2.4. The best-fit point in the CMSSM with μ > 0 (< 0) has tan β ∼ 51 (36), and tan β ∼ 39 is preferred in the NUHM1 with μ > 0. All these points have relatively large values of m 0 and m 1/2 , but the likelihood functions of these models are quite flat, and each of the models also has a local minimum of the likelihood function at low mass, with smaller tan β and small χ 2 ≤ 1 relative to the global minimum. We present 95 % CL lower limits on mg (the gluino mass), mq R (the average over the right-handed squark masses of the first two generations), mt 1 (the light scalar top mass) and mτ 1 (the light scalar tau mass) in each of these models. In each case, we find that the lighter stopt 1 may be significantly lighter than the other strongly interacting sparticles.
The structure of this paper is as follows. In Sect. 2 we discuss the updated MasterCode framework and the more important changes in our implementations of the experimental constraints. There are no significant changes in the ways we treat the constraints not discussed explicitly. In Sect. 3 we describe the results of our fits within the CMSSM and NUHM1. Finally, in Sect. 4 we summarise our conclusions and discuss the prospects for future studies of these and other SUSY models, in particular during the LHC Run 2.  [85][86][87] and SSARD [88] codes, which are interfaced using the SUSY Les Houches Accord [89,90]. The MasterCode is used to construct a global likelihood function that includes contributions from electroweak precision observables, flavour measurements, the cosmological dark matter density and direct searches for dark matter, as well as the LHC Higgs-mass measurement and E / T searches.

Implementation of MultiNest
There has been a major overhaul of the MasterCode since [33], with the aim of simplifying its use and facilitating its application to different SUSY models. The most important change in its implementation has been to use the MultiNest algorithm [67,68]  tion approaches, it can be used to sample well multidimensional parameter spaces, and thereby estimate efficiently and robustly frequentist confidence intervals. The main requirements for our purposes are that no nodes of high likelihood are missed, and that the regions with low χ 2 are well sampled. For the scans used in this paper we use the ranges 0 < m 0 < 7000 GeV, 0 < m 1/2 < 4000 GeV, 2 < tan β < 68 and −5000 GeV < A 0 < 5000 GeV 3 in the CMSSM, for both signs of μ, thereby extending significantly the m 0 range compared to [33]. In the case of the NUHM1, we use the same ranges for m 1/2 , tan β and A 0 , sample 0 < m 0 < 4000 and study the range −5 × 10 7 GeV 2 < m 2 H < 5 × 10 7 GeV 2 , restricting our attention to μ > 0. The total numbers of points sampled in the CMSSM with μ > 0 and μ < 0 and the NUHM1 are 6.8×10 6 , 5.3×10 6 and 1.6×10 7 , respectively. In all cases, the best-fit points were checked by running Minuit [91] on the parameter space, and the differences in total χ 2 between MultiNest and Minuit were 1 %. In this analysis we make several changes in our implementations of the constraints, of which the most important are described in the following subsections.

The ATLAS 20/fb jets + E / T constraint
The ATLAS Collaboration has made public preliminary updates of their SUSY searches using the entire available 8 TeV data set, including the results of many different searches targeting different E / T final states and topologies. Here we follow the same prescription as in [33], restricting ourselves to using the 0-lepton + 2 to 6 jets + E / T search [5]. This is done in order to ensure that the limits presented by ATLAS in the CMSSM m 0 , m 1/2 plane for tan β = 30 and A 0 = 2m 0 can be extrapolated to other values of tan β and A 0 in the ranges used in our scan. As in [33], we have performed a dedicated validation to check that the 0-lepton E / T limit reported in [5] is quite independent of tan β and A 0 . As was to be to be anticipated given the similarity of the search methodologies between the ATLAS 0-lepton analyses at 7 and 8 TeV, we find very similar results to [33]. Therefore, we assume that the 95 % CL exclusion contour in the (m 0 , m 1/2 ) plane presented in [5] may be used irrespective of tan β and A 0 , and we apply a penalty term to points in our scan according to their distance from the stated 95 % CL limit, using the same scaling function as in [33].
We do not use other LHC SUSY searches (such as those relying on lepton and/or b-jet signatures) in our analysis, for two reasons. One is that in the CMSSM and associated models the jet + E / T searches give the strongest constraints in the most relevant region of parameter space, as we discuss in more detail below. The other is that these are more sensitive to the model parameters, whereas the 0-lepton + jets + E / T search sensitivity is known [33] to be independent of the details of the models we study.

The Higgs mass constraint
In view of the relatively large value of the Higgs mass [92,93], M h = 125.7 ± 0.4 GeV (where the quoted uncertainty is purely experimental) and the stronger lower limits on sparticle masses from direct LHC searches [5] within the CMSSM and NUHM1, the calculation of the Higgs boson masses using FeynHiggs has been improved [69] to achieve a higher accuracy for large stop-mass scales. The calculations implemented in FeynHiggs 2.8.7, which we used previously 4 included the full one-loop contributions and the leading and subleading two-loop corrections. The calculations included in the new version FeynHiggs 2.10.0 used here [69] include a resummation to all orders of the leading and next-to-leading logarithms of the type log(m t /m t ) (where m t denotes the geometric average of the two scalar top masses), based on the relevant two-loop Renormalisation-Group Equations (RGEs) [94]; see [95] and references therein for details. The effects of this new correction start at the three-loop order. It has been ensured that the resummed logarithms, which are obtained in the MS scheme, are correctly matched onto the one-and two-loop corrections in the on-shell scheme that were already included previously [69]. The main effect is an upward shift of M h for stop masses in the multi-TeV range, as well as the possibility of a refined estimate of the theoretical uncertainty that is incorporated in our global fits. This shift in M h relaxes substantially the constraints from the Higgs mass on the CMSSM and NUHM1 and related models [64].
A numerical analysis in the CMSSM including leading three-loop corrections to M h using the code H3m [96][97][98]) was presented in [99]. It was shown that the leading threeloop terms can have a strong impact on the interpretation of the measured Higgs-mass value in the CMSSM. Here, with the new version of FeynHiggs, we go beyond this analysis by including (formally) subleading three-loop corrections as well as a resummation to all orders of the logarithmic contributions to M h ; see above.
The new version of FeynHiggs also includes an updated estimate of the theoretical uncertainty, M h | FH , due to missing higher-ordercontributions to M h [69], which is typically 4 This version was an extension of the publicly available FeynHiggs 2.8.6, which differed in the conversion of the trilinear coupling A b from the DR scheme to the on-shell (OS) scheme. This issue was treated in an improved way in FeynHiggs 2.9.5. The implementation of FeynHiggs used here slightly differs from the public FeynHiggs 2.10.0, with a small difference of 0.5 GeV in the M h calculation. in the range 1.0 to 1.5 GeV in the favoured regions of the parameter spaces we sample. The theoretical uncertainty is to be incorporated in the global χ 2 function via a contribution of the form (1) Conservatively, in this paper we assume a fixed value M h | FH = 1.5 GeV in our evaluation of (1), pending a more complete evaluation of M h | FH in a future version of FeynHiggs.
We do not include in our analysis the current measurements of Higgs production strengths. The measurements in individual channels currently have accuracies measured in the tens of %, whereas the best-fit models we discuss later make predictions that differ from the SM predictions by only a few %, typically considerably less than the current theoretical uncertainties. For this reason, at the present stage the Higgs measurements do not have important impact on the specific models studied models studied in this paper, though they may be relevant now for other models and in the future also for the CMSSM and related models as the precisions of the Higgs measurements improve.

The BR(B
To date, the most precise measurements of BR(B s → μ + μ − ) and BR(B d → μ + μ − ) have been provided by the CMS Collaboration [3]: and the LHCb Collaboration [4]: These numbers correspond to time averaged (TA) branching fractions, 5 and are in good agreement with the SM TA expectations [103][104][105] (see also [106,107]): An official combination of the CMS and LHCb results can be found in the conference note [108]: In all new physics (NP) models with minimal flavour violation (MFV) [109,110], including the CMSSM and the NUHM1, BR(B s → μ + μ − ) and BR(B d → μ + μ − ) can deviate from their corresponding SM predictions, but their ratio remains fixed at the SM value [111,112]: 6 We exploit this property to combine BR(B s → μ + μ − ) and BR(B d → μ + μ − ) measurements into a single constraint in the CMSSM (NUHM1) parameter space. In particular, for each of the four measurements in (2) and (3) we determine the ratio that is independent of q in the context of MFV models. The four constraints are then combined into a single weighted mean (hereafter denoted R exp μμ ), taking into account the correlations between the different measurements. It should also be noted that SuFla computes directly a theoretical prediction for R μμ , allowing one to separate the theory uncertainties into three sources: SUSY theory uncertainties (which are negligible), the uncertainty from (6), and those affecting the SM prediction of the branching fractions.
The CMS Collaboration has provided an estimate of R CMS μμ = 1.01 +0. 31 −0.26 [113] by combining its BR(B s → μ + μ − ) and BR(B d → μ + μ − ) measurements (and using the SM values for the branching ratios in [106,107]). Here we construct a joint likelihood for the four measurements (2)-(3) using correlation coefficients between BR(B s → μ + μ − ) and BR(B d → μ + μ − ) of −50 % in CMS and +3 % in LHCb [114]. In order to do that, we use Eqs. (3) and (4) of [108], where the uncertainty coming from the ratio of hadronisation fractions of the b quark (the ratio of probabilities of the b quark to hadronise into a B 0 or a B s meson, and only common systematic between LHCb and CMS), f d / f s , is factorised out.
Thus the full likelihood becomes where B R s and B R d stand for BR(B s → μ + μ − )and BR(B d → μ + μ − ), respectively. The log-likelihoods of quantities with asymmetric errors are approximated using a treatment equivalent to formula (4) in [115]. Equation 8 assumes the same value of f d / f s for CMS and LHCb, which is the same assumption done in [108]. We then reparameterise the joint likelihood as a function of the single parameter of  (2) and (3), as calculated using the prescription described in the text (blue solid line) compared with the contribution calculated previously in [33] (dashed red line) interest, R ex p μμ , imposing the constraint in (6) and using the SM prediction of BR(B s → μ + μ − ) from [103][104][105].
Our final estimate after profiling on the theory uncertain- We have checked that our approach reproduces with good accuracy both the results in (5) when we use (8) and the R CMS μμ value when we drop LHCb input from the final likelihood. Of course we performed both checks using the SM values for the branching ratios in [106,107]. The contribution this function makes to the global χ 2 function is shown as the solid blue line in Fig. 1, where it is compared with the contribution calculated previously in [33] (dashed red line).

The dark matter constraints
There are two important dark matter constraints on the CMSSM and NUHM1 parameter spaces. One is the cosmological dark matter density χ h 2 = 0.1198 ± 0.0026 estimated from Planck data [70], and the other is the upper limit on the spin-independent elastic cold dark matter scattering cross section σ SI p from the LUX experiment [9], which is stronger by a factor ∼2 than that from the XENON100 experiment [8] in the range of neutralino masses relevant to this study. Upper limits on the spin-dependent cross section do not impinge on the parameter spaces of the models we study [30].
Previously, we used Micromegas 2.4.5 to calculate χ h 2 , which we checked gave results similar to the independent SSARD code in the regions of interest. Here we use Micromegas 3.2 [85][86][87]. The recent results from the Planck satellite [70] refine the previous observational estimate of χ h 2 , but this does not alter significantly the implications for other observables.
We compute the elastic scattering cross section, σ SI p , using [88]. There are important uncertainties in the calculation of σ SI p , which are now incorporated in the present analysis, which are also computed using [88]. There are two major sources for these uncertainties that we review here briefly. The first is the uncertainty related to the shift in the nucleon mass due to finite quark masses, σ 0 = 36 ± 7 MeV, and the second is due to the uncertainty in the π -nucleon sigma term, πN , which we take here as 50 ± 7 MeV. The spin-independent matrix element forχ 0 1 -nucleon scattering is proportional to a parameter f N that can be written as where the parameters f with [116,117] An expression for α 3q in terms of supersymmetric model parameters is given in [118]: it does not contribute significantly to the uncertainty in the calculation of the cross section, which is dominated by uncertainties in hadronic parameters. These matrix elements are all directly proportional to πN . It is well known that the elastic cross section is very sensitive to the strange scalar density in the nucleon, Indeed, f T s is proportional to πN y, and hence the uncertainties in both πN and σ 0 enter. Our calculation of the uncertainty in the elastic cross section propagates the independent uncertainties in πN and σ 0 , as well as uncertainties in the quark-mass ratios m d /m u and m s /m d , though the latter two are much smaller than the former two. For a more complete discussion of these uncertainties, see [119]. While the uncertainty in σ SI p is often attributed to the uncertainty in πN , there is an almost equally large contribution to the uncertainty in σ SI p coming from σ 0 , particularly in the determination of the important strangeness contribution, f T s .
We display in Fig. 2 the contribution to the global χ 2 function that we calculate on the basis of the LUX 90 % CL upper limit on the spin-independent cross section σ SI p [9], without (red points) and with (blue points) taking into account the Fig. 2 The contribution to the global χ 2 function that we calculate on the basis of the LUX 90 % CL upper limit on the spin-independent cross section σ SI p [9], without (red points) and with (blue points) taking into account the theoretical uncertainty in the calculation of σ SI p . The horizontal blue bar exhibits the effect of this uncertainty in the hadronic matrix element on the calculation of σ SI p for a specific CMSSM point uncertainty in the calculation of σ SI p . The horizontal blue bar is a representative example of the effect of the theoretical uncertainty in the hadronic matrix element on the calculation of σ SI p for one of the CMSSM points.

CMSSM fits
We now present the results of our new CMSSM fit using the above new inputs, considering first the case of μ > 0 and then the case of μ < 0. In each case, we present first some illustrative parameter planes, and then some one-dimensional likelihood functions.

Parameter planes for μ > 0
In each of the parameter planes in Fig. 3, the best-fit point is indicated by a green star, the χ 2 = 2.30 contour that corresponds approximately to the 68 % CL is shown as a red line, and the χ 2 = 5.99 contour that corresponds approximately to the 95 % CL is shown as blue line. The results of the current fit are indicated by solid lines and solid stars. The results of fits to the same set of data constraints as used in our previous paper [33], but using current theoretical codes including  We see in the upper left panel of Fig. 3 that the current CMSSM fit has two disjoint χ 2 = 2.30 contours in the (m 0 , m 1/2 ) plane, one enclosing an 'island' at relatively low masses, centred around (m 0 , m 1/2 ) ∼ (500, 1000) GeV, and a larger 'continent' extending from (m 0 , m 1/2 ) ∼ (500, 1500) GeV to larger mass values, beyond the range m 0 < 6000 GeV studied here. As we discuss below, the lowmass 'island' lies in the stau-coannihilation region, where the (g − 2) μ contribution to the global χ 2 function is reduced, whereas in the high-mass 'continent' the relic density is brought within the cosmological range by rapidχ 0 1 annihilations due to direct-channel heavy Higgs boson resonances. Our current best-fit point lies in the outer 68 % CL region and has (m 0 , m 1/2 ) ∼ (5650, 2100) GeV. There is a single χ 2 = 5.99 contour enclosing both the inner and the outer 68 % CL regions. We note that the global χ 2 function is quite flat in the outer region and very similar to the χ 2 value for the SM. 7 The lower limit on m 1/2 at small m 0 ∼ 500 GeV is provided mainly by the ATLAS 20/fb search for events with E / T and 2 to 6 jets, whereas at large m 0 3000 GeV there are several relevant ATLAS limits using different event topologies with jets, leptons, b quarks and E / T . These are quite sensitive to tan β and A 0 , and they have little impact on the preferred regions of the CMSSM parameter space, so we do not include them in our analysis. The lower limit on m 0 and the low-mass 'island' corresponds to the stau LSP boundary and the nearby coannihilation strip. The region at large m 0 and m 1/2 containing the best-fit point is in the rapid-annihilation funnel region, with the upper bound on m 1/2 being provided by the cosmological constraint on χ h 2 . The region at small m 1/2 and large m 0 is in the focus-point region.
Looking now at the (m 0 , tan β) plane in the upper right panel of Fig. 3, we see that the low-mass coannihilation 'island' corresponds to values of tan β 30, whereas the Footnote 7 continued and A 0 = 100 GeV, discarding the contributions from M h and χ h 2 and evaluating the σ SI p contribution as 0.14, which yields χ 2 = 36.5. See the text for a detailed comparison of the current fit to that in [33] lower-χ 2 part of the high-mass continent corresponds to a band with larger values of tan β ∼ 50 in the rapidannihilation funnel region, connected to a 'continental shelf' in the focus-point region extending to lower tan β when m 0 3500 GeV. The best-fit point has tan β ∼ 50 and lies in the funnel region.
The lower left panel of Fig. 3 shows the (tan β, m 1/2 ) plane. As already mentioned, the only ATLAS 20/fb jets + E / T limit we use is that on E / T + 2 to 6 jets: other limits using topologies with leptons and/or b jets could have an impact when m 1/2 500 GeV, depending in particular on the value of tan β and/or A 0 . We have not attempted to model these limits, but note that they would not affect the 68 % CL region displayed.
Finally, the lower right panel of Fig. 3 displays the (M A , tan β) plane of the CMSSM. We see that in the lowmass coannihilation 'island' typical values of M A ∼ 1500 to 2500 GeV. The best-fit point has a similar value of M A , but with a much larger value of tan β. The band at large tan β corresponds to the rapid-annihilation funnel region. It is clear that the larger values of m 0 seen in the other panels correspond to large values of M A ∼ 2500 GeV and more.
We compare in Fig. 4 the results of the current analysis (solid lines and filled stars) with the results that were shown in [33] using the previous data set and the previous implementations of the constraints (dashed lines and open stars). As already mentioned, the strengthened ATLAS E / T constraints with 20/fb of data at 8 TeV have had little impact except to strengthen the lower limit on m 1/2 at low m 0 ∼ 500 GeV. At larger m 0 , the range of m 1/2 values is broader than that shown in [33], largely because of changes in the MicrOMEGAs code between versions 2.4.5 and 3.2, which allow points with the correct dark matter relic density to be found in this region of parameter space. This is the case, in particular, for the best-fit point we now find, which lies in the outer 68 % CL region and has (m 0 , m 1/2 ) ∼ (5650, 2100) GeV. Our new treatment of the uncertainty in σ SI p discussed in Sect. 2.6, combined with larger Higgs mass found in FeynHiggs 2.10.0, has the effect of disfavouring the focus-point region [121] less than in [33], leading to an expansion in the region allowed at the 95 % CL at large m 0 and small m 1/2 . The extension of the CMSSM 95 % CL region to larger m 0 in the left panel of Fig. 4 is due to the extended sampling range we use here: the MultiNest technique used here does not have a big impact beyond improving the density of sampling.
The global likelihood function calculated in [33] had two local minima with almost equal values of χ 2 . The new ATLAS constraint and the new implementations of the M h and σ SI p constraints combine to slightly disfavour the local minimum in the low-mass 'island' by only ∼ 0.7 in χ 2 compared to the global minimum at (m 0 , m 1/2 ) ∼ (5650, 2100) GeV, where the main contribution comes from the ATLAS E / T constraint. As in the case of our previous analysis [33], the BR(B s,d → μ + μ − ) constraint does not play a large rôle in the current fit. Its main importance is at large tan β and small m 0 , m 1/2 and M A , but the low-mass 'island' has small tan β and the current best-fit point has large values of m 0 , m 1/2 and M A . Table 1 summarises the values of χ 2 and the locations of the best-fit points found in the current analysis in the lowmass ('island') and high-mass ('continent') regions of the CMSSM parameter space with μ > 0. 8 We see that although the minimum value of χ 2 in the 'continent' is smaller than in the 'island', the difference is less than unity and hence is not significant. As already mentioned, the 'island' best- Table 1 The best-fit points found in global CMSSM fits for both signs of μ and an NUHM1 fit with μ > 0, using the ATLAS 20/fb jets + E / T constraint [5], and the combination of the CMS [3] and LHCb [4] constraints on BR(B s,d → μ + μ − ) [108], as well as an update of the FeynHiggs calculation of M h and a more conservative treatment of the hadronic matrix element uncertainties in σ SI p , as discussed in the text. The results for the CMSSM with μ > 0 and the NUHM1 are compared with those found previously in global fits based on the ATLAS 7-TeV E / T data and the previous experimental constraint on BR(B s,d → μ + μ − ), and with a current SM fit made using the procedure discussed in the text. We list the parameters of the best-fit points in both the low-and the high-mass 'islands' in Figs. 3, 8 and 12. The last column shows the fine-tuning for each point; see the text for more details. We note that the overall likelihood function is quite flat in both the CMSSM and the NUHM1, so that the precise locations of the best-fit points are not very significant, and we do not quote uncertainties. For completeness, we note that in the best NUHM1 fits m 2 H = −2.54 × 10 7 GeV 2 at the low-mass point and m Comparing with the best fit found previously in the CMSSM using the ATLAS 7-TeV E / T constraint and the previous BR(B s,d → μ + μ − ) measurement, we see that the best-fit χ 2 has increased by about 2.1 and the 'island' χ 2 by about 3. Thus the pressure exerted by the ATLAS 20/fb jets + E / T and BR(B s,d → μ + μ − ) constraints does not change significantly the overall picture for the CMSSM. Specifically, the values of m 0 and m 1/2 at the new best-fit 'island' point are not very different from those at the previous CMSSM best-fit point, though the values of A 0 and tan β have changed substantially. Also shown for comparison is the value of χ 2 for the "SM", as calculated using the MasterCode by setting m 0 = m 1/2 = 15 TeV. We see that the CMSSM is unable to reduce χ 2 much below the "SM" value, with a similar fit probability.

Characteristics of the best-fit points for μ > 0
The last column in Table 1 lists the amounts of fine-tuning of SUSY parameters required at each of the best-fit points, as calculated using the measure suggested in [127,128]. The CMSSM and NUHM1 low-mass best-fit points do not require significantly more fine-tuning than those found in our previous analysis [33], whereas the high-mass points do. However, we caution against using this as a strong argument in favour of the low-mass points, in view of the subjective nature of the fine-tuning measure and the ambibuities in its interpretation. Table 2 gives more details of the contributions to the global χ 2 function from different observables in the CMSSM at the high-and low-mass best-fit points, compared with our implementation of the SM. At both the high-and the lowmass points, the M h measurement makes a small contribution to the global χ 2 function. We see that the low-mass point has less tension with (g − 2) μ , and it is favoured by both M W and BR(B s,d → μ + μ − ), in particular, whereas the highmass point is preferred by BR(b → sγ ) and ATLAS 20/fb jets + E / T , in particular. The "SM" fit is noticeably worse for (g − 2) μ and M W .

One-dimensional likelihood functions for μ > 0
We now present the one-dimensional χ 2 likelihood functions for various particle masses and other observables when μ > 0, which are shown as continuous lines in Fig. 5 (the dotted lines are discussed below). The upper left panel displays the χ 2 function for mg. We see that it falls essentially monotonically for mg 1000 GeV, a feature that masks the structures seen in the upper left panel of Fig. 3. The onedimensional projection merges the low-mass 'island' and the high-mass 'continent' that are separated in the (m 0 , m 1/2 ) plane of the CMSSM. It is to be expected that the χ 2 function continues close to zero also at larger values of mg.
The χ 2 function for mq R seen in the upper right panel of Fig. 5 exhibits more structure, with a local minimum at mq R ∼ 2200 GeV, a local maximum at mq R ∼ 3000 GeV, and then an essentially monotonic fall at larger mq R . The appearance of the local minimum can be understood by remembering that m 2 q R ∼ m 2 0 + 5m 2 1/2 , so that the value of mq R is fixed along elliptical contours in the (m 0 , m 1/2 ) plane. The local minimum in the χ 2 function for mq R corresponds to an ellipse passing through the red 'island' in the upper left panel of Fig. 3, and the local maximum corresponds to an ellipse passing between the 'island' and the 'continent'. However, we should emphasise that neither the local mini- Table 2 Summary of the contributions of the most relevant observables to the global χ 2 function at the best-fit high-and low-mass points in the CMSSM (with both signs of μ) and NUHM1 (with μ > 0), including the recently updated observables ATLAS 20/fb jets + E / T , BR(B s,d → μ + μ − ) and the LUX upper limit on dark matter scattering. As noted in parentheses, within the SM, χ 2 ∼ 1.5 is found in [73] due to the (small) tension between the measured value of M h and the precision electroweak data Observable χ 2 CMSSM μ > 0 (high)  mum nor the local maximum is very significant, since they have χ 2 ∼ 1, 4 relative to the minimum value of χ 2 . Similar features are seen in the χ 2 function for the mass of the lighter stop squark, mt 1 , as seen in the lower left panel of Fig. 5. However, in this case the local minimum appears at a lower mass mt 1 ∼ 1000 GeV, and the local maximum is also at a lower mass mt 1 ∼ 2000 GeV, reflecting the fact that the isomass contours for mt 1 and mq R are different. As in many other models, we find that thet 1 is likely to be considerably lighter than the other strongly interacting sparticles. This is due to a large mixing in the scalar top sector, driven by the relatively large value of M h .
Similar local structures can also be seen in the χ 2 function for the lighter stau, mτ 1 , as seen in the lower right panel of Fig. 5. In this case, the local minimum is at mτ 1 ∼ 450 GeV, nearly degenerate with the lightest neutralino, and placing theτ 1 and other sleptons beyond the reach of an e + e − collider with E CM 900 GeV. We also find that χ 2 > 9 for mτ 1 < 300 GeV. We also see a second local minimum of χ 2 at mτ 1 ∼ 1000 GeV, which arises from the lobe at (m 0 , m 1/2 ) ∼ (1500, 2200) in the upper left panel of Fig. 3. However, we emphasise that these observations are very model-dependent.
We now comment briefly on the differences between the one-dimensional likelihood functions found in our analysis of the current data, and those found using the same implementations of the M h and σ SI p constraints for the data set used in [33], shown in Fig. 5 as dotted lines. The current likelihood functions for mg, mq R , mt 1 and mτ 1 are generally higher at small masses, where the ATLAS E / T search has the most impact, but are similar at high masses. Figure 6 displays the χ 2 functions for the mass of the lightest SUSY Higgs boson, M h , shown in the left panel and the mass of the pseudoscalar Higgs boson, M A , shown in the right panel. We see that the likelihood for M h is well Table 3 The 95 % CL lower limits (in GeV) on various sparticle masses in the CMSSM with both signs of μ and the NUHM1 with μ > 0. We emphasise that these limits are specific to the models studied. In the case of the CMSSM with μ < 0 and the NUHM1, the parentheses indicate the approximate locations of small mass ranges where the χ 2 function dips briefly below the 95 % CL maximised close to the measured Higgs mass. The likelihood for M A is very flat for M A 1000 GeV, with χ 2 rising rapidly to reach > 9 for M A < 500 GeV, and it is very similar to the likelihood found using the same data set as in [33].
On the basis of these one-dimensional likelihood functions we can establish 95 % CL lower limits on mg, mq R , mt 1 and mτ 1 for the CMSSM with μ > 0, which are listed in the second column of Table 3. Reflecting the relatively large values of m 0 favoured in this analysis, we see that the lower limit on mq R is considerably stronger than that on mg. On the other hand, thet 1 could be substantially lighter than the other strongly interacting sparticles.
The left panel of Fig. 7 displays the likelihood function for BR(B s,d → μ + μ − ), which is seen to be minimised close to the SM value. The rise at larger BR(B s,d → μ + μ − ) is largely due to the direct experimental constraint on this quantity, but the steep rise at lower BR(B s,d → μ + μ − ) is due to the other constraints on the CMSSM, which are hard to reconcile with R μμ < 1. The rise at large BR(B s,d → μ + μ − ) Fig. 7 The one-dimensional χ 2 likelihood function in the CMSSM for μ > 0 for BR(B s,d → μ + μ − ) (left) and the (mχ0 1 , σ SI p ) plane (right). In both panels, the solid lines are derived from a global analysis of the present data and the dotted lines are derived from a reanalysis of the data used in [33], using the implementations of the M h and σ SI p con-straints discussed in Sect. 2. In the right panel, the red lines denote the χ 2 = 2.30 contours, the blue lines denote the χ 2 = 5.99 contours in each case and the filled (open) green star denotes the corresponding best-fit point found from the data set used in [33] is less steep, reflecting the evolution in the measurement of BR(B s,d → μ + μ − ). The right panel of Fig. 7 displays the (mχ0 1 , σ SI p ) plane, again with solid (dashed) lines representing the current analysis (the constraints of [33]), respectively, with the filled (open) green star denoting the corresponding best-fit point whereas the red (blue) lines representing 68 (95) % CL contours, respectively. We see that a range 10 −47 cm 2 σ SI p 10 −43 cm 2 is allowed at the 95 % CL, and the best-fit point yields a value in the middle part of this range ∼ 10 −45 cm 2 . The mass of mχ0 1 at the best-fit point is 935 GeV. Since the favoured range of mχ0 1 is high, in this and the other models discussed later, and the predicted values of σ SI p correspondingly small, the search for spin-independent dark matter scattering does not have a strong impact on the global fits.

Comparisons between analyses
We restrict our attention here to the only other analysis that incorporates the latest ATLAS 20/fb jets + E / T constraint. Preliminary results from a new global frequentist analysis of the CMSSM with μ > 0 within the FITTINO framework have recently been presented [47]. The best-fit point found in [47] is very similar to the best-fit point we find in the low-mass region of the CMSSM with μ > 0. However, the regions of the parameter space favoured at the 68 and 95 % CL in the FITTINO analysis do not extend to values of (m 0 , m 1/2 ) as large as those we find in the present analysis. In addition to ATLAS 20/fb jets + E / T , this analysis also uses HIGGSSIGNALS [120] to derive constraints from the Higgsmass and signal-strength measurements. The latter do not change substantially the results, since the Higgs rate predictions in the favoured regions of the CMSSM parameter space, which are in the in the decoupling regime, 9 are quite similar to those in the SM and do not vary significantly. 10 3.2 CMSSM with μ < 0 The case μ < 0 has been studied less than μ > 0 (but see, e.g., [65,66,124,125]), for various reasons: It worsens the discrepancy between the experimental value of (g − 2) μ and the SM calculation, it is in general more restricted by BR(b → sγ ) and it yields a smaller value of M h for fixed values of the other CMSSM parameters. However, since the ATLAS 20/fb jets + E / T and other constraints require relatively large values of m 0 and m 1/2 where the SUSY contribution to (g − 2) μ and BR(b → sγ ) are small, it is appropriate to reconsider the μ < 0 case.

Parameter planes with μ < 0
We see in the upper left panel of Fig. 8 that there are three regions of the (m 0 , m 1/2 ) plane that are allowed at the 95 % level, two small 'reefs' at relatively low masses (m 0 , m 1/2 ) ∼ (300, 1000) and (600, 2000) GeV and a more extensive 'continent' at larger masses m 0 4000 GeV. The lower-mass 'reef' is in the stau-connihilation region, as in the μ > 0 case, but the higher-mass 'reef' is in the stop-coannihilation region. Compared to the high-mass 'continent' in the rapid-annihilation funnel and focus-point regions, the 'reef' has smaller contributions to the global χ 2 Fig. 8 As in Fig. 3, but for μ < 0 in the CMSSM. For the reason discussed in the text, only the ranges tan β ≤ 40 are displayed. The yellow star in the lower left panel marks the best-fit point in the CMSSM with μ < 0, which is out of the ranges of the other panels function for some electroweak and flavour observables, but is disfavoured by ATLAS 20/fb jets + E / T . The best-fit point in the CMSSM for μ < 0 is shown as a yellow star: it is located in the high-mass 'continent', in the focus-point region.
The (m 0 , tan β) plane for μ < 0 is shown in the upper right panel of Fig. 8. 11 Here we see that the low-mass 'reefs' are restricted to 5 tan β 15, whereas the 'continent' extends over all tan β 8. In the lower left panel of Fig. 8, we see in the (tan β, m 1/2 ) plane that the 'reefs' and 'continent' merge in this projection of the CMSSM parameter space with μ < 0. We also see that 500 GeV m 1/2 2500 GeV is allowed at the 95 % CL for the range m 0 < 6000 GeV studied here. 12 The small region within the red 68 % contour does not appear in the other panels, because it corresponds to values of m 0 > 6000 GeV and M A > 4000 GeV, which 11 Here and in subsequent panels, we restrict attention to tan β ≤ 40. The parameters required to solve the electroweak vacuum conditions for larger values of tan β are relatively unstable, often with values of m 0 and A 0 outside the ranges studied here, so the limits on the ranges of these parameters that we sample imply that our sampling is incomplete for tan β > 40 and μ < 0. 12 As in the μ > 0 case, we have not studied in detail the sensitivity to tan β of the lower bound on m 1/2 due to ATLAS 20/fb jets + E / T searches with leptons and/or b jets, which are not important near the 'reef' or the best-fit point for μ < 0.
are not displayed in the other panels of Fig. 8. Finally, in the lower right panel of Fig. 8 we see in the (M A , tan β) plane that only in the 'reefs' are values of M A 3000 GeV are allowed at the 95 % CL when tan β ≤ 40. 13 The 'reef's are again clearly separated at relatively small values of tan β, with a restricted range of M A ∈ (2000, 3000) GeV.

Characteristics of the best-fit points for μ < 0
We display in Table 1 the characteristics of the best-fit points in the CMSSM with μ < 0 in the low-mass 'reef' region and the high-mass 'continent'. Unlike the case of the CMSSM with μ > 0,t 1 coannihilation is important at the best-fit point in the 'reef' region, andχ ± 1 ,χ 0 2 ,χ 0 3 coannihilation at the best-fit 'continental' point. In both cases, the global χ 2 function is somewhat higher than in the corresponding regions for μ > 0, by ∼ 3.1 in the low-mass region and by ∼ 1.5 in the high-mass region. The main origins of the differences can be seen in Table 2. The high-mass model with μ < 0 receives larger contributions from (g − 2) μ , M W and σ SI p , whereas there are larger contributions from (g − 2) μ and M h in the low-mass case, compensated only partially by smaller χ 2 contributions from BR(b → sγ ) and BR(B s,d → μ + μ − ). As a result, the best-fit CMSSM points with μ < 0 have higher χ 2 and lower fit probabilities than the SM.

One-dimensional likelihood functions for μ < 0
We display in Fig. 9 the one-dimensional χ 2 functions for various sparticle masses in the CMSSM with μ < 0. We see in the upper left panel that the χ 2 function for mg falls essentially monotonically as mg → 5000 GeV towards χ 2 ∼ 2.5 relative to the global minimum. The best fit for μ < 0 has χ 2 ∼ 1.8 at mg ∼ 5300 GeV, and hence it is not seen in this plot.
On the other hand, the one-dimensional χ 2 function for mq R , shown in the upper right panel of Fig. 9 has a very different form. After falling initially to χ 2 ∼ 4, there is a local maximum at mq R ∼ 2000 GeV with χ 2 ∼ 8. This is followed by a region where χ 2 falls again to ∼ 4, followed by a sharp rise to χ 2 > 9. Finally, the χ 2 function falls again below χ 2 = 9 when mq R > 3800 GeV and contin-ues falling with increasing mq R . The low-mass structures are in the 'reef' regions and the high-mass fall is in the 'continental' region. Similar features are seen in the χ 2 function for mt 1 , but at lower masses, in the lower left panel of Fig. 9. The χ 2 function for mτ 1 shown in the lower right panel of Fig. 9 exhibits sharp local minima at mτ 1 ∼ 400 (associated with the dip in the gluino χ 2 at 2000 GeV) and 2000 GeV (associated with the high gluino-mass region), followed again by a decrease across the 'continent' at large masses.
We display in Fig. 10  122 GeV do not acquire a heavier χ 2 penalty is due to the theoretical uncertainty in the calculation of M h , which we take to be 1.5 GeV. The χ 2 function for M A has a local minimum at M A ∼ 2000 GeV followed by a rise to a local maximum at M A ∼ 2300 GeV and then a decrease towards χ 2 ∼ 4 when M A ∼ 4000 GeV. These one-dimensional likelihood functions can be used to set 95 % lower limits on various sparticle masses by requiring , σ SI p ) plane (right). In both panels, the solid lines are derived from a global analysis of the present data using the implementations of the M h and σ SI p con-straints discussed in Sect. 2. In the right panel, the red lines denote the χ 2 = 2.30 contours, the blue lines denote the χ 2 = 5.99 contours in each case and the filled yellow star denotes the corresponding best-fit point χ 2 < 4 relative to the global minimum for the CMSSM, which occurs for μ > 0 as discussed earlier. These lower limits are tabulated in the third column of Table 3. We indicate in parentheses the approximate locations of limited ranges of masses where the χ 2 function dips briefly below the 95 % CL. Figure 11 shows the one-dimensional χ 2 functions for BR(B s,d → μ + μ − ) (left panel) and σ SI p (right panel) for μ < 0. We see that BR(B s,d → μ + μ − ) is expected to be very similar to the SM value, reflecting the previous observation that the lowest χ 2 values for μ > 0 are attained in the 'continent' at large sparticle masses and large M A , and the secondary minima in the 'reefs' at low masses has small values of tan β. We also see that the preferred values of σ SI p for μ < 0 are ∼ 10 −44 to 10 −45 cm 2 at large mχ0 1 , whereas σ SI p is 10 −48 cm 2 in the 'reef' region.

The NUHM1 with μ > 0
We now turn our attention to the NUHM1, concentrating on the case μ > 0, since our study of the CMSSM indicates that this sign is still preferred by the data, albeit less strongly than in [33]. Figure 12 displays our selection of NUHM1 parameter planes, with the same conventions for solid/dashed lines as in Fig. 3. We see in the upper left panel that the likelihood function is relatively flat for m 1/2 2000 GeV, and that there is a low-mass 'peninsula' extending down to (m 0 , m 1/2 ) ∼ (500, 1200) GeV, which is analogous to the 'island' in the CMSSM. The 68 % CL region extends to val- See the text for a detailed comparison of the current fit to that in [33]. The red lines denote χ 2 = 2.30 contours (corresponding approximately to the 68 % CL), and the red lines denote χ 2 = 5.99 (95 % CL) contours ues of m 1/2 > 4000 GeV, which was not the case in the CMSSM. This is because the NUHM1 is able to satisfy the χ h 2 constraint for larger values of m 1/2 than are possible in the CMSSM, thanks to the extra degree of freedom associated with the soft SUSY-breaking contribution to the Higgs masses. This permits values of μ or M A that allow χ h 2 to fall within the astrophysical range even if m 1/2 is large. We also note that the NUHM1 can satisfy the electroweak vacuum conditions in regions of the parameter space with m 2 0 < 0, though we have not studied this possibility in any detail.

NUHM1 parameter planes
The differences in Fig. 12 between the results of the current analysis (solid lines and filled stars) with our current implementations of the data constraints used in [33] are relatively minor. On the other hand, looking back at the right panel of Fig. 4 where our current NUHM1 results are compared with those shown previously in [33], cf, the dashed lines and open star, we see that both the 68 % and the 95 % CL regions now extend to much larger m 1/2 . This is largely the result of sampling an extended range in m 2 H , as well as using FeynHiggs 2.10.0 to calculate M h . As in the CMSSM case shown in the left panel of Fig. 4, the extension of the 95 % CL region to lower m 1/2 at large m 0 is due to the new implementation of the dark matter scattering constraint discussed in Sect. 2.6.
The upper right panel of Fig. 12 displays the (m 0 , tan β) plane in the NUHM1. We see a general trend for the preferred range of tan β to increase with the value of m 0 1000 GeV. Values of tan β as low as ∼ 5 are allowed in the 'peninsula' region. In the (tan β, m 1/2 ) plane shown in the lower left panel of Fig. 12, we see that values of tan β ∼ 5 to 30 are preferred when m 1/2 2000 GeV, whereas larger values of m 1/2 are associated with tan β 15. Finally, we see in the lower right panel of Fig. 12 that values of M A 500 GeV are generally preferred, with most of the favoured region appearing in a lobe with M A 2000 GeV.

Characteristics of the best-fit points in the NUHM1
The best-fit point in the 'continental' region has nearly degenerateχ 0 1 ,χ 0 2 andχ ± 1 , since μ m 1/2 and the LSP is nearly a pure higgsino, and theτ 1 is ∼ 20 GeV heavier in this case. Fig. 13 The one-dimensional χ 2 likelihood functions in the NUHM1 for μ > 0 for mg (upper left), mq R (upper right), mt 1 (lower left) and mτ 1 (lower right). In each panel, the solid line is derived from a global analysis of the present data, and the dotted line is obtained from a reanalysis of the data used in [33], using the implementations of the M h and σ SI p constraints discussed in Sect. 2 Thusχ ± 1 ,χ 0 1,2 coannihilation is important in fixing χ h 2 , butτ 1 coannihilation is not negligible. As could be expected from the shape of the 68 % CL region in the lower right panel of Fig. 12, whilstχ ± 1 ,χ 0 2 coannihilation is important in most of the 'continental' region, different dynamical processes are important in different regions of the NUHM1 parameter space. For example,τ 1 coannihilation and rapid annihilation via direct-channel poles are both important in the lobe where M A ∼ 1000 GeV and tan β ∼ 10, which includes the best-fit point to the previous data set (open star). On the other hand, only rapid annihilation via direct-channel poles is important in the lobe where M A ∼ 1000 GeV and tan β ∼ 30, and onlyχ ± 1 ,χ 0 2 coannihilation is important in the narrow strip where M A ∼ 1000 GeV and tan β ∼ 55. Finally, both χ ± 1 coannihilation and rapid annihilation via direct-channel poles are important in the lobe where M A 2000 GeV and tan β 60.
We see in Table 2, comparing the contributions to the global χ 2 functions for the high-mass points in the NUHM1 and the CMSSM with μ > 0, that the NUHM1 point has a noticeably smaller χ 2 contribution from M W . Comparing the low-mass points in the NUHM1 and the CMSSM with μ > 0, we see that the NUHM1 point has smaller χ 2 contributions from BR(b → sγ ) and ATLAS 20/fb jets + E / T , in particular. The M h constraint does not make an important contribution to χ 2 at either of the NUHM points.
3.5 One-dimensional likelihood functions in the NUHM1 Figure 13 displays the one-dimensional χ 2 functions for various sparticle masses. The likelihood function for mg (upper left panel) decreases essentially monotonically until mg ∼ 2600 GeV, which is followed by a local maximum at mg ∼ 3500 GeV. The global minimum is at 6800 GeV and hence not visible on this plot. The χ 2 function for mq R shown in the upper right panel of Fig. 13 has similar behaviour. On the other hand, the χ 2 function for mt 1 , shown in the lower left panel of Fig. 13, manifests an important local minimum at mt 1 ∼ 700 GeV followed by a local maximum at mt 1 ∼ 1000 GeV before exhibiting a second local minimum and local maximum at mt 1 ∼ 2000 and 2700 GeV, respectively. Finally, the χ 2 function for mτ 1 ,  Fig. 14 The one-dimensional χ 2 likelihood functions in the NUHM1 for μ > 0 for M h (left) and M A (right). In each panel, the solid line is derived from a global analysis of the present data, and the dotted line is derived from a reanalysis of the data used in [33], using the implementations of the M h and σ SI p constraints discussed in Sect. 2

Fig. 15
The one-dimensional χ 2 likelihood function in the NUHM1 for μ > 0 for BR(B s,d → μ + μ − ) (left) and the (mχ0 1 , σ SI p ) plane (right). In both panels, the solid lines are derived from a global analysis of the present data, and the dotted lines are derived from a reanalysis of the data used in [33], using the implementations of the M h and σ SI p con-straints discussed in Sect. 2. In the right panel, the red lines denote the χ 2 = 2.30 contours, the blue lines denote the χ 2 = 5.99 contours in each case and the filled (open) green star denotes the corresponding best-fit point seen in the lower right panel of Fig. 13, exhibits a lowmass local minimum at mτ 1 ∼ 500 GeV associated with the above-mentioned 'peninsula' followed by a local maximum at mτ 1 ∼ 700 GeV, and then it falls to a shallow minimum at mτ 1 ∼ 1000 GeV, eventually rising slowly at larger masses.
Turning now to the one-dimensional χ 2 functions for the SUSY Higgs bosons shown in Fig. 14, we see in the left panel that the likelihood function for the mass of the lightest supersymmetric Higgs boson M h is maximised very close to the experimental value, though with tail extending to lower and higher masses reflecting the theoretical uncertainty in the calculation. As for M A , we see in the right panel of Fig. 14 that the likelihood function is rather flat for M A 1000 GeV. The 95 % CL lower bounds on mg, mq R , mt 1 , mτ 1 and M A inferred from the one-dimensional χ 2 functions in Figs. 13 and 14 are tabulated in Table 3. As in the CMSSM cases stud-ied, thet 1 may be significantly lighter than the other strongly interacting sparticles.
We see in Fig. 15 that the one-dimensional χ 2 function for BR(B s → μ + μ − ) is minimised close to the SM value. The NUHM1 offers very little scope for values of BR(B s → μ + μ − ) below this, but values larger than in the SM are not so strongly disfavoured. The right plot of Fig. 15 shows the NUHM1 results in the (mχ0 1 , σ SI p ) plane. Similar ranges of mχ0 1 and σ SI p are favoured as for the CMSSM with μ > 0. Finally we show some in Fig. 16 a comparison of the onedimensional likelihoods for the mass of the lightest charginõ χ ± 1 for the CMSSM with μ > 0 and the NUHM1. (The corresponding likelihood function for theχ 0 2 would be similar.) In the CMSSM theχ ± 1 likelihood function is quite flat between around 800 GeV and 2100 GeV, before rising steeply after that. The NUHM1 likelihood is similar to the CMSSM one In each panel, the solid line is derived from a global analysis of the present data, and the dotted line is obtained from a reanalysis of the data used in [33], using the implementations of the M h and σ SI p constraints discussed in Sect. 2 below 1000 GeV, but maintains its flatness out to 2500 GeV. This reflects the flatness of the NUHM1 likelihood at out to large values of m 1/2 .

Summary and prospects
We have presented in this paper analyses of the CMSSM with both signs of μ and the NUHM1 with μ > 0 that take into account all the relevant constraints from the first run of the LHC with ∼ 5/fb of luminosity at 7 TeV and ∼ 20/fb of luminosity at 8 TeV, as well as flavour and precision electroweak observables and the first results from the LUX search for spin-independent dark matter scattering [9]. We have sampled the model parameter spaces using the MultiNest technique, made SUSY model calculations of M h using version 2.10.0 of the FeynHiggs code, and we have taken account of uncertainties in these calculations and in the estimation of hadronic matrix elements for dark matter scattering.
It is a general feature of our analysis that we find larger values of m 0 and m 1/2 to be allowed than were found in our previous analyses, largely because of our updated interpretation of the experimental M h constraint using FeynHiggs 2.10.0 and the newer version of MicrOMEGAs that we use. The parameters of the best fits we find in the CMSSM and NUMH1 are displayed in Table 1: we note that they also have larger values of m 0 and m 1/2 than were favoured previously. Also shown for comparison are the model parameters for local minima of the global χ 2 functions at low masses, which are disfavoured by the ATLAS 20/fb jets + E / T constraint, in particular. We note that all the favoured CMSSM and NUHM1 model points can accommodate the measured value of M h . None of the SUSY models studied has a global χ 2 value that is much lower than the SM. This is because none of the SUSY models discussed reduces significantly the contributions to the global χ 2 functions from the observables that make the largest contributions to the global χ 2 functions in the SM fit, namely (g − 2) μ , A fb (b), A (SLD), σ 0 had , K and BR(B u → τ ν τ ), as seen in Table 2. The 95 % CL lower limits on sparticle masses found in our CMSSM and NUHM1 analysis are displayed in Table 3. We see that gluino masses above ∼ 1300 GeV are preferred in the models analyzed. The right-handed squark mass is restricted to even higher values, because of the preferred values of m 0 , whereas the lighter stop squark may be significantly lighter. The lighter stau slepton may also be relatively light in the CMSSM and NUHM1 with μ > 0. On the other hand, the heavier Higgs bosons A, H and H ± are all expected to have masses above 500 GeV in these models.
Estimates of the discovery reach of the LHC at 14 TeV have been provided in [126]. With 300/fb of luminosity, the 5-σ discovery reach for squarks and gluinos should extend to (mq R , mg) ∼ (3500, 2000) GeV in the CMSSM with μ > 0, would include the low-mass 'reef' in the CMSSM with μ < 0, and would reach the first local minimum of the χ 2 function in the NUHM1 with μ > 0, at (mq R , mg) ∼ (2500, 3000) GeV. The discovery range with 3000/fb of luminosity would extend a few hundred GeV further, and it would be very similar to the 95 % CL exclusion reach with 300/fb. The reach for 95 % CL exclusion with 3000/fb would extend several hundred GeV further still, e.g., to (mq R , mg) ∼ (4000, 2700) GeV in the CMSSM with μ > 0.
We conclude that large parts of the preferred parameter regions of the CMSSM and NUHM1 are accessible in future runs of the LHC, although the strongly interacting sparticle masses might be so high as to escape the searches at the LHC. That said, we re-emphasise that all the likelihood estimates made in this paper and the estimates of the LHC physics reach are specific to the models studied, and are quite model depen-dent. The approach we have followed here for constructing the global likelihood function can easily be extended to other models, a subject to which we will turn in future work.