Fitting the Higgs to natural SUSY

We present a fit to the 2012 LHC Higgs data in different supersymmetric frameworks using naturalness as a guiding principle. We consider the MSSM and its D-term and F -term extensions that can raise the tree-level Higgs mass. When adding an extra chiral superfield to the MSSM, three parameters are needed determine the tree-level couplings of the lightest Higgs. Two more parameters cover the most relevant loop corrections, that affect the hγγ and hgg vertexes. Motivated by this consideration, we present the results of a five parameters fit encompassing a vast class of complete supersymmetric theories. We find meaningful bounds on singlet mixing and on the mass of the pseudoscalar Higgs mA as a function of tan β in the MSSM. We show that in the (mA, tan β) plane, Higgs couplings measurements are probing areas of parameter space currently inaccessible to direct searches. We also consider separately the two cases in which only loop effects or only tree-level effects are sizable. In the former case we study in detail stops’ and charginos’ contributions to Higgs couplings, while in the latter we show that the data point to the decoupling limit of the Higgs sector. In a particular realization of the decoupling limit, with an approximate PQ symmetry, we obtain constraints on the heavy scalar Higgs mass in a general type-II Two Higgs Doublet Model.


Introduction
With the 2012 run reaching its conclusion our knowledge of Higgs couplings is considerably improving. Some purists might argue that the new resonance [1,2] is not yet experimentally proven to be the Higgs boson, but all evidence goes in that direction and even more so after the recent updates by the LHC experiments [3,4]. In the previous round of LHC results, namely 4.7 fb −1 at √ s = 7 TeV and 5.9 fb −1 at √ s = 8 TeV, the errors in the Higgs' couplings were dominated by statistical uncertainties that are expected to scale-down with increased luminosity. Nonetheless, these measurements were already precise enough to make meaningful statements about broad classes of models [5][6][7][8][9]. After the addition of approximately 7 fb −1 at √ s = 8 TeV, we find it interesting to explore the consequences of Higgs measurements in the predictive and theoretically motivated framework of natural supersymmetry. Several groups have already fitted Higgs' rates in different contexts, focusing mainly on simplified settings in which all the couplings were determined by one or two free parameters [10][11][12][13][14], concentrating also on the cases of light stops or the type-II two Higgs doublet model (2HDM) [5,10]. Here we adopt a different perspective, discussing the implications of Higgs measurements in complete natural supersymmetric theories.
In section 2 we briefly review our definition of natural supersymmetry. In this context, it is necessary to extend the Higgs sector of the MSSM to obtain the observed Higgs mass m h ≈ 126 GeV without incurring excessive fine-tuning. Therefore we consider D-term and F -term models that modify the 2HDM structure of the MSSM potential.

JHEP03(2013)043
We study the new tree-level effects, distinguishing between the cases in which an approximately type-II 2HDM structure is preserved and those in which the mixing with extra light states in the Higgs sector is sizeable. We find that in general three couplings are enough to parameterize all the tree-level deviations in Higgs rates. Loop corrections to hgg and hγγ vertices can also be sizeable, since naturalness considerations require stops and Higgsinos to be fairly light. Two extra parameters are needed to include these effects. Therefore we perform a five dimensional fit to the data. Although naturalness has driven us to consider a certain spectrum and a certain subset of supersymmetric effects on Higgs couplings, the fine-tuning constraints are not imposed in fits. The general five dimensional framework parameterizes a vast class of complete theories, that may or may not be finetuned. We also present fit results in several simplified settings, considering the cases in which: (1) only the tree-level mixing with the second Higgs doublet plays a role, (2) only tree-level effects are important, but both the second Higgs doublet and an extra singlet are responsible for deviations in the couplings and (3) only MSSM loop effects from stops and charginos are sizeable.
The fitting procedure is described in section 3 together with the data taken as input. In section 4 we present the results of the fit and discuss their relevance in light of direct searches performed at the LHC.

The natural SUSY framework
Deviations in the Higgs couplings are probed by the measurements of the Higgs rates, which are determined by partial widths. We denote the modifications to the partial widths as with i = t, V, G, γ, b, τ standing for top, massive vector gauge boson, gluon, photon, bottom and tau, respectively. In general it is hard to write the r i 's in terms of few supersymmetry breaking parameters. However if we restrict to natural theories it is possible to identify a limited number of relevant effects, as was shown in [15]. As stated above, we do not impose fine-tuning bounds in the fits, but we still find naturalness as a good guiding principle to select models and simplified scenarios in which fit results can be interpreted. For concreteness, consider at most a ∆ −1 = 10% tuning, where This restricts, in the MSSM, an unmixed stop to be lighter than about 600 GeV. Taking into account also the relation between the Higgs mass and the Higgsino mass: −m 2 h /2 = |µ| 2 + . . ., limits one chargino to be lighter than 300 GeV. More precisely, the upper bounds obtained from fine-tuning take the form [16] µ 290 GeV m h 125 GeV 10% ∆ −1 , and Λ is the scale at which supersymmetry breaking effects are mediated to the MSSM. At tree level, the bound on µ translates into a bound on the lightest chargino Other Standard Model (SM) superpartners can be much heavier, with masses ranging from a few TeVs to tens of TeVs.
However, raising the Higgs mass to 126 GeV in the MSSM through stop radiative corrections, implies a tuning much worse than 10% [17,18]. To accommodate a 126 GeV Higgs boson, natural supersymmetry must be non-minimal, allowing for additional interactions that modify the Higgs sector. The new interactions deform the Higgs quartic potential at tree level and change the Higgs couplings beyond the significant deviations (relative to the SM) that can already occur in the MSSM. In the following we will consider both extensions that leave the type-II 2HDM structure of the potential unaltered and shift quartic couplings relative to the MSSM, and models in which the 2HDM relations between light Higgs couplings are relaxed.
Nonetheless, naturalness not only complicates the problem, as seen above, it also limits the size of A t and µ. This greatly simplifies the problem of determining the r i 's as a function of soft supersymmetry breaking parameters. First and foremost, the type-II 2HDM structure of the MSSM Higgs potential is not strongly affected by non-holomorphic loop corrections in a natural theory, regardless of our definition of tuning. As an example we can take the definition at the beginning of this section and ask for a 10% tuning, then the corrections to the bottom Yukawa for tan β ≤ 40 (and r b = O(1) at tree-level) are negligible [15]. So in the natural MSSM and its D-term extensions, Higgs couplings can be expressed in terms of only two parameters at tree-level, as in a general type-II 2HDM. Adding a singlet chiral superfield, as for example in the NMSSM, increases the number of parameters to three, if we assume all couplings to be CP conserving. The additional degree-of-freedom parameterizes the mixing between the new CP even state and the lightest scalar Higgs: c φ ≡ S|h .
We have seen that naturalness also prefers stops and charginos to be fairly light, and for this reason we consider their loop contributions to the Higgs to digluon/diphoton partial widths, introducing two new parameters. They can be chosen as δrt G , characterizing the stop contribution to the dimension-5 Higgs-gluon-gluon coupling and δrχ ± γ , characterizing the chargino contribution to the dimension-5 Higgs-photon-photon-coupling. It was suggested that light staus can be responsible for the enhancement of the h → γγ rate observed at the LHC [19][20][21], but we do not consider this possibility explicitly in what follows. We also neglect sbottom loop corrections and charged Higgs loop corrections. In all three cases we expect contributions larger than a few % only in extreme corners of the parameter space that imply a tuning much worse than 10% in our definition [15]. In this regions charge breaking minima can also be generated [21]. Nonetheless, in the 5 dimensional fits, r G and r γ are treated as independent parameters. Thus we capture any loop effect on these vertexes.

JHEP03(2013)043
In terms of the five parameters discussed above, the partial width modifiers are given by where we have chosen r b and r t to parametrize the 2HDM tree-level effects. Note that r τ = r b assumes that non-holomorphic corrections to the Higgs interactions with down quarks and leptons are small relative to the tree-level value, as discussed above. Here, and in what follows, we always neglect Higgs decays to non-SM particles, in particular to neutralinos. In section 4.1 we review all tree-level predictions in more detail, starting with a generic 2HDM with c φ = 0, and later introducing a superpotential for the singlet. In sections 4.4 and 4.5 we study the bounds on loop corrections from naturalness and direct searches and compare them with the results of the fit.

Data and fitting procedure
In order to constrain the supersymmetric parameters, we take into account all the available Higgs channel rates [3,4]. The rates depend on the product of the overall production crosssection and branching ratio for the particular channel. The results are typically reported as a confidence interval on the event rate relative to the SM prediction, denoted byμ. We take all theμ's at the two LHC best fit values for the mass: m h = 125.8 GeV for CMS and m h = 126.0 GeV for ATLAS. For the Tevatron h → bb rate we take the value ofμ at m h = 126.0 GeV. In table 1 we list the analyses relevant for our fit, regrouped by channel, and the corresponding values of theμ's. In any given model, the signal strengths are determined by the r i 's. For the LHC, the four relevant production modes and their respective theoretical dependance are • Vector boson fusion (qq → hqq): • Vector boson associated production (qq → hV ): • Top associated production (gg → htt): σ htt /σ SM htt = |r t | 2 .
In some cases, a channel can include events from several production modes; for instance, the dijet tagged signature in table 1 is dominantly produced via vector boson fusion, but contains a non-negligible contribution from gluon fusion. The collaborations have made public all the numbers needed to assess the composition in terms of the physical production modes in the references listed in table 1. In the fit we use the numbers provided channel by channel, but we find a good uniformity between different final states and the two experiments, giving roughly 75 − 80% of VBF and 20 − 25% of GF in the dijet tagged categories. The relative fraction of the untagged mode are close to be the relative crosssection fractions of the different channels. Therefore, where efficiencies are not publicly JHEP03(2013)043 available (e.g. ATLAS untagged WW), we use the ratio between cross-sections to determine the composition of the sample. The V H and ttH categories can be taken as pure. The fits to the data are performed minimizing the χ 2 where µ i (r j ) are predicted rates and σ i 's are two-sided errors (depending on the sign of µ i −μ i ). By performing the χ 2 test we are assuming that the likelihood functions forμ follow an approximate two-sided gaussian distribution and the correlations can be neglected, both of which have been shown to be valid approximations [5,13]. Furthermore we have compared our results with the two dimensional fits released by the collaborations, obtaining a good agreement.
In the following we refer to "preferred" and "allowed" regions. Preferred regions are obtained by varying an N -parameters subset of the r i 's, while fixing the other parameters to their SM-values. Strictly speaking, allowed regions are found including all the r i 's, the possibility of flavor non-universality (for instance r b = r τ ), vertex structures different from the SM and the presence of an invisible width. At this stage, with limited precision in the measurements, we find a theoretically inspired (see sections 2 and 4) five dimensional fit a reasonable approximation of the most generic setting. In this framework we obtain lower dimensional confidence intervals by treating the remaining independent parameters as nuisances (given a lower dimensional point the other parameters are varied to give the best possible fit).

Predictions and results
In this section we discuss in more detail the predictions of natural supersymmetry and compare them with the results of the fit to Higgs data. We consider four simplified scenarios: a type-II 2HDM , a type-II 2HDM with a SM singlet, light stops and light charginos both in the Higgs decoupling limit. At the end of the section we discuss the more general five parameters fit.

The type-II two Higgs doublet model
If the superpartners are relatively heavy, i.e. m t 1 500 GeV and m χ ± 1 200 GeV, loops corrections to Higgs couplings are at most of the order of 10% in a broad region of parameter space, provided that the stop mixing is limited by naturalness [15]. Therefore it is easy to decouple these effects without incurring excessive fine-tuning and without the need to single out narrow corners of parameter space.
With this in mind, we begin by exploring the tree-level corrections in the MSSM and its extensions, using the generic type-II 2HDM scalar potential and Yukawa couplings [34,35], We have not included the non-VBF τ τ categories for ATLAS. The information on the cross correlation is not provided making difficult to asses the contamination of the combined result. To extract the VBF only rate from the CLs plot we follow the procedure described in [14]. where A · B ≡ ij A i B j , and we adopt the usual definition We also take all couplings in the Higgs potential to be CP-conserving. In this setting the Higgs couplings to the SM particles depend on two parameters, the ratio of Higgs vevs [34] and the mixing angle of the CP-even Higgs states, defined by where h is the observed 126 GeV state. At tree-level, r b = r τ , and in the MSSM (4.5) Figure 1. Left: χ 2 contours corresponding to the 95% confidence level in the (sin α, tan β) plane (for a type-II 2HDM). Note that in the MSSM sin α < 0. The black dashed line corresponds to the decoupling limit of the MSSM α = β − π/2. Right: ∆χ 2 vs. ξ(≡ α − β + π/2). The red lines mark the 68% and 95% confidence levels heights of the ∆χ 2 . In black the offset between the SM value of the χ 2 and the 2HDM value at the minimum. The discontinuity in the χ 2 is physical and it is illustrated by eq. (4.6). Note that in the tree-level MSSM ξ ≥ 0.
Without loss of generality, we can take tan β > 0 [36], with eq. (4.5) implying sin α < 0. We require that the top Yukawa does not blow up above the electroweak scale, which imposes the lower bound tan β ≥ 1 [37]. As discussed in section 2, demanding loop corrections to the type-II 2HDM structure to be negligible provides the additional constraint tan β ≤ 40 (that can be relaxed or made stronger depending on the level of fine-tuning that we allow in the theory). Since the measured rates broadly agree with the SM, the results of the fit, which are shown in figure 1, point to the decoupling limit of the model, ξ ≡ α − β + π/2 ≈ 0. This translates into the bound |ξ| 0.1 at 95% C.L. Note that the discontinuity in the χ 2 is physical and can be easily understood by expanding Higgs couplings for small ξ Even for extremely small values of ξ the correction to r b can be significant for large tan β and depends on the sign of ξ. The one dimensional χ 2 was obtained by treating tan β as a nuisance, so the profiling is selecting a tan β that can contribute an observable deviation on the side where r b is depleted, while for ξ < 0 a tan β giving the smallest possible deviation is singled out. For this reason the offset between the two sides is the same as the offset between the deepest minimum and the SM. The result, albeit disappointing from the perspective of finding new physics, allows us to consider a particular realization of the decoupling limit in which more interesting statements can be made on the mass of the second Higgs doublet. In general it is difficult to translate the fit results for α and β into a direct mass exclusion or into an exclusion on JHEP03(2013)043 the Lagrangian parameters. However, if the 2HDM potential respects an approximate PQ symmetry, or a Z 2 : Assuming that the second Higgs doublet is parametrically heavier than H u , we can expand the light Higgs couplings in the Z 2 breaking spurion B/M 2 1 , where 1 and λ 35 = λ 3 + λ 5 . Up to leading order B/M 1 = 1/ tan β, thus the expansion converges rather fast in 1/ tan β, e.g., only a O(10%) error for tan β = 3. This setting is extremely predictive, giving (in the MSSM and many of its extensions where λ 7 = 0) [38] (4.8) In this limit, the corrections to r b depend on a single combination of quartics, λ 35 , and the mass of the heavy scalar Higgs m H . Note that if we assume the approximate PQ and some degree of custodial SU(2) symmetry, the differences in the masses of the heavy Higgses are small [39]: m H m A m H + . Thus, to first approximation m H can be considered the mass of the full heavy doublet.
The result of the fit is depicted in figure 2. In the MSSM, λ 3 = −(g 2 +(g ) 2 )/4 ≈ −0.14 and λ 5 = 0 at tree-level, while the loop corrections are small. This leads to a preferred region at 95% C.L. with m H 370 GeV. However the tan β expansion clearly introduces an error that in the MSSM goes in the direction of relaxing our bound 2 In the following section we obtain a limit on m A as a function of tan β using the full treelevel expression in the MSSM together with the most important loop corrections. The effect of taking into account the O(1/ tan 2 β) variations of the other couplings is much smaller and will be commented upon while repeating this exercise in the five dimensional fit case and in the MSSM with the full tan β dependence. This expansion is also relevant to a vast class of models in which the supersymmetric Higgs sector is not necessarily minimal, but it is still well described by a type-II 2HDM at low energy. This is the case for D-term models raising the Higgs mass via extra gauge interactions and also for some areas of the parameter space of F -term models.
We define D-term models as the class of theories in which the two Higgs doublets are charged under new gauge interactions. We will restrict to gauge theories that do not forbid the µ term, otherwise additional F -terms should be present to generate it. In D-term JHEP03(2013)043 with some positive definite coefficients in front. ξ D contains both gauge symmetry breaking terms that can be absorbed by the soft masses m 2 1 and m 2 2 and new fields charged under the gauge groups. Thus the new effect will be a contribution to λ 3 that is always negative, resulting in stronger bounds on m H . This is illustrated in figure 2, by the yellow band corresponding to variation of the D-term contribution to the Higgs mass from 0 to m h −m Z . Note that to obtain an observable effect from the extra gauge sector, it is necessary to have a large supersymmetry breaking mass for the new gauge bosons. If the heavy gauge bosons are integrated out supersymmetrically, there is a shift in the quartic , where µ is limited by naturalness while M V 3 TeV due to EWPTs [41]. For F -term models the applicability of the approximation in eq. (4.8) is limited, but we can consider a superpotential of the form where S is either a singlet under the SM gauge groups or an SU(2) triplet with Y = 0.
Integrating out S and neglecting v 2 /m 2 s effects (where m s is a supersymmetry breaking mass), the leading observable modification to the Higgs potential is

JHEP03(2013)043
which again shifts λ 3 by a definite sign, but in the opposite direction with respect to Dterm models. Note that in this case small tan β may be preferred to lift the tree-level Higgs mass, m 2 h ≤ m 2 Z cos 2 2β + 2|λ| 2 g 2 + (g ) 2 sin 2 2β (4.13) and the expansion given in eq. (4.8) may not be theoretically motivated. Furthermore to leading order, in gauge and gaugino mediation, the singlet supersymmetry breaking mass m s in the NMSSM vanishes at the mediation scale [42], so it is necessary to introduce some deformations of the simplest scenarios to lift the mass above the other soft terms at the electroweak scale. There are nonetheless large areas of parameter space, in which the expansion is reliable, that give the correct Higgs mass. For instance, including stop corrections with mt 1 ≈ mt 2 = 380 GeV and X t = 0, we get m h = 126 GeV for all values of λ between 1 and 2, varying tan β between 3 and 8 [15]. Considering this range for λ would correspond to taking λ 35 between 1 and 4, thus giving a favored region, that for λ = 1, is already m H 460 GeV at 95% C.L.

Heavy Higgses and m h = 126 GeV in the MSSM
As discussed above only two parameters are needed in the MSSM, at tree-level, to specify all Higgs couplings. Therefore we can express the fit constraints in the (m A , tan β) plane and compare our results with direct searches. To do so, the requirement of naturalness needs to be abandoned to allow for the stops' radiative corrections to raise the Higgs mass to 126 GeV. Considering only the leading loop correction (which is the only one that survives in the limit µ → 0) only the H u − H u element of the CP even mass matrix receives a correction. The tree-level relations in the previous section are correspondingly shifted to [43] tan 2α tan 2β = where δm 2 can be approximated by the leading stop correction but in the following we always fix δm 2 to the m A and tan β dependent value giving the correct Higgs mass Before describing the results of the fit we review current LHC searches for MSSM Higgses. Direct searches for the neutral MSSM Higgses are currently performed in the φ → τ + τ − [44,48], bbφ → bbµµ [49] and bbφ → bbbb [50] channels. The strongest bound is set by the recently updated CMS measurement [44] and ranges from m H 250 GeV for tan β = 5 to m H 700 GeV for tan β = 40. The sensitivity vanishes below tan β = 5. These searches benefit from large tan β, both thanks to new production mechanisms that become important (b and bb associated production for instance) and from the increase in the branching ratio to τ + τ − . We also include the LEP bound [45] in our comparison with Higgs rates, but we do not consider the implications of searches for the charged Higgs that are currently not as sensitive as the φ → τ + τ − one plus the LEP constraint. The experimental collaborations use the m max h scenario [46] to set their limits, and varying the soft supersymmetry breaking parameters might lead to modifications of the observed bounds [47]. The effects are stronger for large tan β mainly due to loop corrections to y b and can not produce any significant gain in sensitivity in the region 3 tan β 12, where production cross sections become too small and our analysis starts to be competitive. Additionally, the collaborations do not impose m h ≈ 126 GeV. Ideally, the bounds from [44] should be reinterpreted as a bound on σ(pp → φ)Br(φ → τ + τ − ), but it is beyond the scope of this work. Figure 3 shows the overlay of these bounds with the results of the fit to the Higgs data. In the slice 2 tan β 5 Higgs rates are probing regions of the parameter space not directly accessible to CMS and LEP. The fit results in figure 3 are an anticipation of the five dimensional fit discussed in section 4.6 since the excluded region was obtained profiling JHEP03(2013)043 the five dimensional χ 2 , treating r γ and r G as nuisances. This is reflected in the fact that for large tan β there is a lower bound m A 250 GeV, in agreement with the result of the PQ expansion in figure 8.

Singlet mixing
Finally, we consider the effects of the singlet mixing with the two Higgs doublets. We take the superpotential (4.11) and do not specify the detailed form of the soft supersymmetry breaking terms or extra superpotential interactions. For simplicity we take all new couplings to be CP conserving and assume the SM-like Higgs to be CP even and the lightest Higgs state in the theory. With these assumptions we do not need to consider exotic decays that could be triggered by approximate U(1)'s, including the U(1) R , generating light pseudo-goldstone bosons in the Higgs sector. This considerably simplifies the problem, leaving the relevant parameters to be contained in the 3 × 3 mixing matrix for the scalar CP even states, where c x = cos x, s x = sin x. In this notation, the tree-level Higgs couplings depend on (α, φ, β), and eq. (4.4) is modified to where cos φ ≡ S|h measures the amount of singlet in the lightest Higgs. The relations (4.5) are only valid in the MSSM, so that m A is no longer determined solely (at tree level) from the mixing parameters α and β. Thus, we do not impose any constraint on α and φ, but still restrict to tan β ≥ 1. In figure 4, we plot ∆χ 2 vs cos 2 φ. Since the presence of singlet mixing will only decrease the observed rates, there is no improvement in the fits over the SM and we obtain an upper bound | cos φ| 0.7 at 95% C.L. It is worth pointing out that even though this bound was derived assuming only tree-level effects are present, it is similar to the bound obtained in the full 5D space, as will be discussed in section 4.6.

Stops
LHC searches, together with naturalness, have stimulated a vast offspring of models of mediation that produce light stops, decoupling all other squarks [51][52][53][54][55][56]. In addition, bottom-up approaches, driven mainly by FCNC constraints, still allow stops to be fairly light [52,57,58]. It is then worth exploring the stop plane in view of Higgs data. Using the Higgs low energy theorem [59][60][61] it is straightforward to obtain the change in the gluon fusion rate from integrating out the stops (neglecting D-terms)  Figure 4. ∆χ 2 as a function of the singlet mixing parameter cos 2 φ. The 68% and 95% confidence levels are marked in red. In black we show the offset between the χ 2 at the minimum and the SM minimum. Figure 5. Right: χ 2 contours at the 68% and 95% confidence levels in the (r G , r γ ) plane. All other parameters are fixed to their SM value. Left: ∆χ 2 vs. r G for the stops fit, with 95% and 99.73% heights marked in red. In gray we show the offset between the χ 2 at the minimum and the SM minimum.
where X t = A t − µ cot β and m 2 are the two eigenvalues of the stop mass matrix . (4.20) The approximation in (4.19) is valid up to order m 2 h /4m 2 t 1 , but is useful to give a qualitative picture of the effect. Nevertheless all the plots and the numbers that are quoted here are obtained with the full MSSM one loop result [43].
The hγγ vertex correction can be computed rescaling the contribution to the gluongluon amplitude with the appropriate factors of N c , Q 2 and t c (the Dynkin index of the color JHEP03(2013)043 representation). Since the SM h → γγ loop is dominated by the W boson contribution, the two corrections are opposite in sign, with roughly a factor of 4 difference in magnitude The current data point to |r G | ≈ 0.9 and |r γ | ≈ 1.4 as can be seen in the left panel of figure 5. This can be realized only in a small corridor of the stop plane depicted in the left panel of figure 6, that is mostly outside of the 10% tuning region and requires the presence of a very light stop mt 1 200 GeV. 3 In this corner of parameter space, the stops' contribution to the gluon fusion rate is opposite in sign and roughly double in magnitude with respect to the top quark. This region is not only tuned by more than 1% if we demand that the gluino be heavier than 1 TeV, but can also generate color breaking minima [63].
The preferred area with |r γ | > 1.1 is rather narrow since the lightest stop mass varies fast for values of X t and mt ≡ mQ 3 = mũ 3 that give an r γ enhancement. This is reflected in a correspondingly rapid variation of δrt G as a function of the soft supersymmetry breaking parameters.
The results of the fit are presented in the right panel of figure 6. The data currently prefer the narrow corridor described above. The 3σ contours instead lie both in the deeper minimum around the large X t region and in the high mt region where stop effects decouple. This second area corresponds to the shallower SM-like minimum of the fit shown in the right panel of figure 5.
Overall the results of the fit are in tension with the requirement of 10% tuning, but point to a region directly accessible at the LHC. However a good fraction of the preferred masses falls where the sensitivity of current searches vanishes: 160 GeV mt 220 GeV [64], possibly motivating additional efforts.

Charginos
At least one chargino should be light m χ ± 1 300 GeV in view of naturalness considerations (see section 2), possibly giving an observable deviation in the h → γγ rate. However, the effects can be decoupled by raising M 2 to O(1 TeV), without introducing more than a 20% tuning [16]. Nonetheless, other than naturalness, there are further motivations to study chargino effects on Higgs couplings.
It possible that these deviations are the only observable effect of the MSSM at the LHC. Many models have the gauginos and higgsinos as the only light particles, with the possible addition of gluinos, as is the case of split supersymmetry [65,66] and models that address the moduli problem and dark matter [67][68][69]. The most studied scenarios typically decouple either gauginos or higgsinos [70,71], which would make the effects on the Higgs couplings vanishingly small, but it is possible to obtain a viable dark matter candidate together with unification keeping both species light [72]. Furthermore they are a reasonable proxy for a member of a SU(2) L multiplet providing a dark matter candidate. The mass splittings in the multiplets are expected to be of the order of a few hundreds of MeV [73], thus leaving JHEP03(2013)043  open only monojet and monophoton searches at 14 TeV. These analyses may be the only viable path to direct exclusion/discovery also in the case of compressed SUSY spectra [74]. Their sensitivity is currently limited only to colored particles [75][76][77] and it was estimated that to probe the electroweak production of particles with masses up to 200 GeV, 300 fb −1 at 14 TeV will be necessary [78]. Therefore it is worth studying the chargino effects on the hγγ coupling in the MSSM, that were previously considered also in [72,[79][80][81] In the MSSM, the chargino contribution is bounded by: −0.3 ≤ δrχ ± γ ≤ 0.13. It might be larger in the NMSSM [82] or in D-term extensions [83], but we do not consider this possibility here. Both limits come from direct chargino searches and can easily be derived by scanning (M 2 , µ, tan β) space with the LEP constraint mχ± > 94 GeV [84]. The bounds are saturated when tan β = 1, restricting tan β ≥ 2 gives −0.2 ≤ δrχ ± γ ≤ 0.1 [15]. It is thus clear that charginos alone can hardly explain the observed h → γγ excess. The results of the fit in figure 7 roughly point to µM 2 > sin 2βm 2 w > 0 where the hγγ coupling is enhanced However the allowed region of the (µ, M 2 ) plane, even for tan β = 1, is at least 2σ away from the best fit to the data, while it disappears as tan β grows and the effect decouples. Note that the leading term in a m 2 h /m 2 χ ± 1 expansion (equation (4.22)) ceases to be a good approximation in most of the region preferred by the fit, where one of the charginos can be extremely light 20 − 50 GeV and already ruled out by LEP. As in the stops case, all the plots and the numbers that are quoted were obtained with the full MSSM one loop result [43].

Five parameters fit
In the previous sections we fit the data to simplified scenarios in which all the Higgs rates could be expressed in terms of one or two parameters. This gave an idea of the level of agreement between the data and some physically motivated corners of the parameter space of natural supersymmetry. In spite of the fact that only three channels per experiments have errors below the 50% level (namely γγ, ZZ and W W untagged) we still find it interesting to explore the more general case in which r γ , r G , r b , r t and r V all play a role. This is a small modification of the four parameter natural MSSM, inspired by the possibility of adding a new singlet, and comes closer to approximating a fit with all couplings left to float. The only difference with respect to a four parameter fit resides in the fact that r t is virtually unconstrained and not artificially limited by its relation with r V . These theory inspired exclusions can be applied also outside of the framework of natural supersymmetry and to facilitate possible attempts we show profiles of the χ 2 for the five couplings in appendix A. In table 2 we show the corresponding one dimensional 95% C.L. intervals, obtained by treating the other parameters as nuisances. The errors on the single couplings vary from 30% to 100%, which alone is not enough to lose all hope of constraining the parameter space of natural supersymmetry. However the impact of leaving all couplings to float is strong for charginos and the whole plane becomes accessible. Similarly for the case of stops, where the contribution in a large part of parameter space can be compensated by a shift in r t .

JHEP03(2013)043
On the other hand not all sensitivity is lost on tree-level couplings. For instance λ 35 = λ MSSM 35 still gives a preferred region with m H 250 GeV and λ = 1 implies m H 320 GeV, as it is shown in figure 8. Similarly the tree-level statements about singlet mixing and the decoupling limit are meaningful. Repeating the exercise in section 4.3, but profiling the full five dimensional χ 2 we obtain | cos φ| 0.8 at 95% C.L. very close to the result of the purely tree-level fit. This is an indication that both the tree-level fit and the five dimensional one are dominated by the lower bound on r V that for α ∼ β − π/2, 4 becomes r V ≈ sin φ. The bound on cos φ becomes much stronger away from the decoupling limit as can be seen in the right panel of figure 9.
We can also extract information on the level of decoupling of the Higgs sector as a function of tan β. In figure 9 the 95% C.L. contour is plotted in the (ξ, tan β) plane. The key message is that large tan β is allowed only for values of ξ close to decoupling, where corrections to r b ≈ 1 − ξ tan β and r t ≈ 1 + ξ/ tan β are small. This is reflected in the (m A , tan β) exclusion discussed in section 4.2 that is even competitive with direct searches. Many of these bounds apply to a vast class of complete theories and indicate that in the MSSM and many of its motivated extensions, tree-level effects in the Higgs sector are already strongly constrained, mainly by the measurements of the h → W W and h → ZZ rates that in our setting are always below their SM value unless ξ = 0 and cos φ = 0.

Conclusion
So far most of the effort has been devoted to simplified settings or complete models in which only one or two parameters enter Higgs rates. We set upon the more ambitious task of constraining a realistic theory using only naturalness to limit the number of parameters entering the game. Some of these ideas were treated already in [15,38], but an application to the data was still missing.
We fit Higgs measurements, first treating independently the tree-level and loop effects, and then performing a five dimensional fit that parameterized complete natural theories such as the NMSSM, the MSSM and its D-term extensions. Considering a generic type-II 2HDM we found a bound on the heavy CP even Higgs mass m H 370 GeV in theories with an approximate PQ symmetry. We also obtained a strong preference for the decoupling limit: |ξ| 0.1, and in theories with an extra singlet a robust constraint on its mixing. When taking into account only loop effects we showed that the data prefer a very light stop ( 200 GeV) with large mixing, while charginos in the MSSM alone can not explain the current h → γγ enhancement.
With the full five parameter fit we found that loop-level statements cease to be valid. However we were still able to draw interesting conclusions on tree-level mixings in the Higgs sector. We found that the mass bound in a type-II 2HDM with moderate to large tan β is relaxed to m H 250 GeV and so is the upper bound on singlet mixing, | cos φ| 0.8. However there is still a strong correlation between the size of tan β and the vicinity of the theory to the decoupling limit. This was translated into an exclusion in the (m A , tan β) plane of the MSSM, that is competitive with direct searches. These bounds hold in a large class of complete theories, where loop corrections to the Yukawa couplings are small.
The comparison between the constraints discussed above and LHC searches indicates that there is still no tension between the two sets of measurements in natural supersym-JHEP03(2013)043 metry and that direct searches are still the best chance of finding superpartners, barring extreme configurations with little visible and invisible energy in the event. The only exception are extra states in the MSSM Higgs sector, that at moderate tan β may be seen first as a deviation in Higgs couplings, rather than produced directly at the LHC.

A Results of the five dimensional fit
This appendix collects the one dimensional profiles of the χ 2 obtained from the five dimensional fit described in section 4.6. In figure 10 we show them with the 68% and 95% intervals marked in orange. The details of the profiling are discussed in section 3. r t , entering only at loop level in well measured rates, is essentially unconstrained due to the compensating effect of r G and r γ .  Figure 10. ∆χ 2 profiles for the five parameters fit (r G , r γ , r V , r b , r t ). Heights corresponding to the 68% and 95% confidence levels are marked in red.

JHEP03(2013)043
Open Access. This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.