The Light and Heavy Higgs Interpretation of the MSSM

We perform a parameter scan of the phenomenological Minimal Supersymmetric Standard Model (pMSSM) with eight parameters taking into account the experimental Higgs boson results from Run I of the LHC and further low-energy observables. We investigate various MSSM interpretations of the Higgs signal at 125 GeV. First, we consider the case where the light CP-even Higgs boson of the MSSM is identified with the discovered Higgs boson. In this case it can impersonate the SM Higgs-like signal either in the decoupling limit, or in the limit of alignment without decoupling. In the latter case, the other states in the Higgs sector can also be light, offering good prospects for upcoming LHC searches and for searches at future colliders. Second, we demonstrate that the heavy CP-even Higgs boson is still a viable candidate to explain the Higgs signal - albeit only in a highly constrained parameter region, that will be probed by LHC searches for the CP-odd Higgs boson and the charged Higgs boson in the near future. As a guidance for such searches we provide new benchmark scenarios that can be employed to maximize the sensitivity of the experimental analysis to this interpretation.


Introduction
The discovery of a Higgs-like scalar boson in Run I of the large hadron collider (LHC) [1,2] marks a milestone in the exploration of electroweak symmetry breaking (EWSB). Within experimental and theoretical uncertainties, the properties of the new particle are compatible with the Higgs boson of the Standard Model (SM) [3]. However, a variety of other interpretations of the Higgs signal are possible, corresponding to very different underlying physics. Here, a prime candidate for the observed scalar boson is a CPeven Higgs boson of the Minimal Supersymmetric Standard Model (MSSM) [4][5][6], as it possesses SM Higgs-like properties over a significant part of the model parameter space with only small deviations from the SM in the Higgs production and decay rates [7].
One of the main tasks of the LHC Run II will be to determine whether the observed scalar boson forms part of the Higgs sector of an extended model. In contrast to the SM, two Higgs doublets are needed in the MSSM to give mass to upand down-type fermions. The extended Higgs sector entails the existence of five scalar bosons, namely a light and heavy CP-even Higgs bosons, h and H , a CP-odd Higgs boson, A, and a pair of charged Higgs bosons, H ± . Mixing between the neutral CP-even and CP-odd states are possible in the CPviolating case [8][9][10][11][12][13], which we will not considered here. At lowest order, the Higgs sector of the MSSM can be fully specified in terms of the W and Z boson masses, M W and M Z , the CP-odd Higgs boson mass, M A , and tan β ≡ v U /v D , the ratio of the two neutral Higgs vacuum expectation values. However, higher-order corrections are crucial for a pre-cise prediction of the MSSM Higgs boson properties and introduce dependences on other model parameters, see e.g. Refs. [14][15][16] for reviews.
Many fits for the Higgs rates in various models and within the effective field theory approach have been performed over the last years, see e.g. Refs. [17,18]. Focusing on the MSSM, recent fits have shown that the interpretation of the observed scalar as the light CP-even MSSM Higgs boson ("light Higgs case") is a viable possibility, providing a very good description of all data [19][20][21][22][23][24]; see also Refs. [25][26][27] for global fits including also astrophysical data. In the light-Higgs case, decoupling of the heavy Higgs bosons (M A M Z ) [28][29][30][31] naturally explains the SM-like couplings of the light MSSM Higgs boson h [7]. Another interesting possibility to explain the SM-like behavior of h without decoupling in the MSSM -the so-called limit of alignment without decoupling -has been outlined in Refs. [32,33], relying on an (accidental) cancellation of tree-level and loop contributions in the CP-even Higgs boson mass matrix. This led to the definition of a specific benchmark scenario [33], which has since been ruled out in the interesting low M A region via pp → H/A → τ + τ − searches [23].
Alternatively, it was demonstrated that the heavy CP-even Higgs boson can also be identified with the observed signal [7,19,[34][35][36][37] ("heavy Higgs case"). 1 In this scenario all five MSSM Higgs bosons are relatively light, and in particular the lightest CP-even Higgs boson has a mass (substantially) smaller than 125 GeV with suppressed couplings to gauge bosons. This led to the development of the low-M H benchmark scenario [43]. This particular scenario has meanwhile been ruled out by ATLAS and CMS searches for a light charged Higgs boson [44,45]. However the heavy Higgs interpretation in the MSSM remains viable, as we will discuss in this paper.
The questions arise, whether, and if so by how much, the MSSM can improve the theoretical description of the experimental data compared to the SM, and which parts of the MSSM parameter space are favored. In a previous analysis [19] we analyzed these questions within the MSSM. We performed a scan over the seven most relevant parameters for MSSM Higgs boson phenomenology, taking into account the data up to July 2012, which showed in particular some enhancement in the measured rate for h → γ γ . We found that both the light and the heavy Higgs case provided a good fit to the data. In particular, the MSSM light Higgs case gave a better fit than the SM when the data in the γ γ channel and low-energy data was included. 1 Such a situation is more common in extensions of the MSSM. In particular, in the NMSSM it occurs generically if the singlet-like CP-even state is lighter than the doublet-like Higgs bosons, see e.g. Refs. [38][39][40][41][42].
The situation has changed in several respects with the release of additional Higgs data by the ATLAS and CMS Collaborations [46]. In particular, the final data obtained in the LHC Run I does not show a significant enhancement over the SM prediction in the γ γ channel anymore, and the heavy Higgs case is much more restricted due to light charged Higgs boson searches. The main aim of the present paper is to study the MSSM Higgs sector in full detail taking into account the current experimental data and in particular the final LHC Run I results, and to propose paths towards a complete exploration of the heavy Higgs case at the LHC in the ongoing Run II. We incorporate the available measurements of the Higgs boson mass and signal strengths, as well as measurements of the relevant low-energy observables. Furthermore we take into account all relevant constraints from direct Higgs and supersymmetric (SUSY) particle searches. We investigate whether the MSSM can still provide a good theoretical description of the current experimental data, and which parts of the parameter space of the MSSM are favored. Within the light Higgs case we analyze the situation with very large M A (decoupling), as well as for small/moderate M A (alignment without decoupling). We also investigate the feasibility of the heavy Higgs case and define new benchmark scenarios in which this possibility is realized, in agreement with all current Higgs constraints.
The paper is organized as follows. We employ the phenomenological MSSM with 8 parameters ( pM SS M 8), which is introduced in detail in Sect. 2. In this section, we also expand upon the theoretical background of the two possible limits that lead to alignment in the CP even Higgs sector, i.e. when one of the CP-even neutral MSSM Higgs bosons behaves like the SM Higgs boson. In particular, we outline how leading two-loop effects on the conditions for alignment can be assessed and present a brief quantitative discussion of these effects. 2 The parameter scan with O(10 7 ) sampling points, the techniques to achieve good coverage, as well as the considered experimental observables and constraints are described in Sect. 3. In Sect. 4 we present our results for the best-fit points and the preferred parameter regions for the light Higgs and the heavy Higgs interpretation. The effects of the Higgs mass and Higgs rates measurements, precision observables, and direct Higgs and SUSY searches are discussed, and the phenomenology of the other MSSM Higgs states is outlined. In particular, in Sect. 4.4 we propose new benchmark scenarios for the study of the heavy Higgs case, which can be probed at the LHC Run II. We conclude in Sect. 5. In Appendix A, we discuss the extent of the tuning associated with the regions of the MSSM parameter space that exhibit approximate Higgs alignment without decoupling. Finally in Appendix B, we provide tables listing the signal strength measurements from ATLAS, CMS and the Tevatron (DØ and CDF) that are included in our analysis.

The MSSM Higgs sector
In this section we briefly review the most important features of the MSSM Higgs sector and motivate the choice of the eight free pMSSM parameters in our scan. We provide a detailed description of the relevant MSSM parameter sectors and our notations, which remain unchanged compared to [19].
In the supersymmetric extension of the SM, an even number of Higgs multiplets consisting of pairs of Higgs doublets with opposite hypercharge is required to avoid anomalies due to the supersymmetric Higgsino partners. Consequently the MSSM employs two Higgs doublets, denoted by H D and H U , with hypercharges −1 and +1, respectively. After minimizing the scalar potential, the neutral components of H D and H U acquire vacuum expectation values (vevs), v D and v U . Without loss of generality, we assume that the vevs are real and non-negative (this can be achieved by appropriately rephasing the Higgs doublet fields). The vevs are normalized such that In addition, we define Without loss of generality, we may assume that 0 ≤ β ≤ 1 2 π (i.e., tan β is non-negative). This can always be achieved by a rephasing of one of the two Higgs doublet fields.
The two-doublet Higgs sector gives rise to five physical Higgs states. The mass eigenstates correspond to the neutral Higgs bosons h, H (with M h < M H ) and A, and the charged Higgs pair H ± . Neglecting possible CP-violating contributions of the soft-supersymmetry-breaking terms (which can modify the neutral Higgs properties at the loop level), h and H are the light and heavy CP-even Higgs bosons, and A is CP-odd.
At lowest order, the MSSM Higgs sector is fully described by M Z and two MSSM parameters, often chosen as the CPodd Higgs boson mass, M A , and tan β. In the MSSM at the tree-level the mass of the light CP-even Higgs boson does not exceed M Z . However, higher order corrections to the Higgs masses are known to be sizable and must be included, in order to be consistent with the observed Higgs signal at 125 GeV [3]. Particularly important are the oneand two-loop contributions from top quarks and their scalar top ("stop") partners. In order to shift the mass of h up to 125 GeV, large radiative corrections are necessary, which require a large splitting in the stop sector and/or heavy stops. For large values of tan β, the sbottom contributions to the radiative corrections also become sizable. The stop (sbottom) sector is governed by the soft SUSY-breaking mass parameter Mt L and Mt R (Mb L and Mb R ), where SU(2) gauge invariance requires Mt L = Mb L , the trilinear coupling A t (A b ) and the Higgsino mass parameter μ.
To achieve a good sampling of the full MSSM parameter space with O(10 7 ) points, we restrict ourselves to the eight MSSM parameters most relevant for phenomenology of the Higgs sector (the scan ranges will be given in Sect. 3.1), under the assumption that the third generation squark and slepton parameters are universal. That is, we take Mq 3 := Mt L (= Mb L ) = Mt R = Mb R , M˜ 3 := Mτ L = Mτ R = Mν τ and A f := A t = A b = A τ . Note that the soft SUSY-breaking mass parameter in the stau sector, M˜ 3 , can significantly impact the Higgs decays as light staus can modify the loop-induced diphoton decay. M˜ 3 is therefore taken as an independent parameter in our scans. Even though the other slepton and gaugino parameters are generally of less importance for the Higgs phenomenology, we scan over the SU(2) gaugino mass parameter M 2 as well as over the mass parameter of the first two generation sleptons, M˜ 1,2 , (assumed to be equal) as these parameters are important for the low-energy observables included in our analysis. The remaining MSSM parameters are fixed, We choose relatively high values for the squark and gluino mass parameters, which have a minor impact on the Higgs sector, in order to be in agreement with the limits from direct SUSY searches. Finally, the U(1) Y gaugino mass parameter, M 1 , is fixed via the GUT relation with s W = 1 − c 2 W and c W = M W /M Z . For more details on the definition of the MSSM parameters, we refer to [19].

The Higgs alignment limit
In light of the Higgs data, which indicates that the properties of the observed Higgs boson are SM-like, we seek to explore the region of the MSSM parameter space that yields a SMlike Higgs boson. In general, a SM-like Higgs boson arises if one of the neutral Higgs mass eigenstates is approximately aligned with the direction of the Higgs vev in field space. Thus, the limit of a SM Higgs boson is called the alignment limit.
To analyze the alignment limit, it is convenient to define where 12 = − 21 = 1 and 11 = 22 = 0, and there is an implicit sum over the repeated SU(2) index j = 1, 2. For consistency of the notation, we denote the corresponding neutral Higgs vevs by v 1 ≡ v D and v 2 ≡ v U . We now define the following linear combinations of Higgs doublet fields, such that H 0 1 = v/ √ 2 and H 0 2 = 0, which defines the socalled Higgs basis [48][49][50]. 3 It is straightforward to express the scalar Higgs potential in terms of the Higgs basis fields H 1 and H 2 , where the most important terms of the scalar potential are highlighted above. The quartic couplings Z 1 , Z 5 and Z 6 are linear combinations of the quartic couplings that appear in the MSSM Higgs potential expressed in terms of H D and H U . In particular, at tree-level, where g and g are the SU(2) and U(1) Y gauge couplings, respectively, c 2β ≡ cos 2β and s 2β ≡ sin 2β. Hence, the Z i are O(1) parameters. One can then evaluate the squared-mass matrix of the neutral CP-even Higgs bosons, with respect to the neutral Higgs states, 3 Since the tree-level MSSM Higgs sector is CP-conserving, the Higgs basis is defined up to an overall sign ambiguity, where H 2 → −H 2 . However, since we have adopted the convention in which tan β is nonnegative [cf. the comment below Eq. (2)], the overall sign of the Higgs basis field H 2 is now fixed.
If √ 2 Re H 0 1 − v were a Higgs mass eigenstate, then its treelevel couplings to SM particles would be precisely those of the SM Higgs boson. This would correspond to the exact alignment limit. To achieve a SM-like neutral Higgs state, it is sufficient for one of the neutral Higgs mass eigenstates to be approximately given by In light of the form of the squared-mass matrix given in (11), we see that a SM-like neutral Higgs boson can arise in two different ways: In particular, the CP-even mass eigenstates are: where c β−α ≡ cos(β −α) and s β−α ≡ sin(β −α) are defined in terms of the mixing angle α that diagonalizes the CP-even Higgs squared-mass matrix when expressed in the original basis of scalar fields, 1 [52]. The case of a SM-like H necessarily corresponds to alignment without decoupling.
In the case of exact alignment without decoupling, Z 6 = 0, the tree-level couplings of the SM-like Higgs boson are precisely those of the Higgs boson of the Standard Model. Nevertheless, deviations from SM Higgs boson properties can arise due to two possible effects. First, there might exist new particles that enter in loops and modify the loop-induced Higgs couplings to gg, γ γ and Z γ . For example, if H is the SM-like Higgs boson, then the charged Higgs boson mass is not significantly larger than the observed Higgs mass, in which case the charged Higgs loop can shift the one-loop induced couplings of the observed Higgs boson to γ γ and Z γ [52]. Similarly, SUSY particles can give a contribution at the loop-level to other, at the tree-level SM-like, couplings. Second, there might exist new particles with mass less that half the Higgs mass, allowing for new decay modes of the SM-like Higgs boson. An example of this possibility arises if H is the SM-like Higgs boson and M h < M H /2, in which case the decay mode H → hh is allowed. Indeed, in the exact alignment limit where s β−α = 0, the tree-level H hh coupling in the MSSM is given by [53] The possibility of alignment without decoupling has been analyzed in detail in Refs. [31][32][33]51,52,[54][55][56] (see also the "τ -phobic" benchmark scenario in Ref. [57]). It was pointed out that exact alignment via Z 6 = 0 can only happen through an accidental cancellation of the tree-level terms with contributions arising at the one-loop level (or higher). In this case the Higgs alignment is independent of M 2 A , Z 1 and Z 5 . This has two phenomenological consequences. First, the remaining Higgs states can be light, which would imply good prospects for LHC searches. Second, either the light or the heavy neutral Higgs mass eigenstate can be aligned with the SM Higgs vev and thus be interpreted as the SM-like Higgs boson observed at 125 GeV.
The leading one-loop contributions to Z 1 , Z 5 and Z 6 pro- is the top quark Yukawa coupling, have been obtained in Ref. [33] in the limit M Z , M A M S (using results from Ref. [58]): where s β ≡ sin β, M S ≡ √ mt 1 mt 2 denotes the SUSY mass scale, given by the geometric mean of the light and heavy stop masses, and In Eqs. (15)-(18), we have assumed for simplicity that μ and A t (as well as the gaugino mass parameters that contribute subdominantly at one-loop to the Z i ) are real parameters. That is, we are neglecting CP-violating effects that can enter the MSSM Higgs sector via radiative corrections.
The approximate expression for Z 6 v 2 given in (17) depends only on the unknown parameters μ, A t , tan β and M S . Exact alignment arises when Z 6 = 0. Note that Z 6 = 0 trivially occurs when β = 0 or 1 2 π (corresponding to the vanishing of either v 1 or v 2 ). But, this choice of parameters is not relevant for phenomenology as it leads to a massless b quark or t quark, respectively, at tree-level. Henceforth, we assume that tan β is non-zero and finite. In our convention, tan β is positive with 0 < β < 1 2 π . We can simplify the analysis of the condition Z 6 = 0 by solving Eq. (15) for ln(M 2 S /m 2 t ) and inserting this result back into Eq. (17). The resulting expression for Z 6 now depends on Z 1 , tan β, and the ratios, Using Eq. (18) to rewrite the final expression in terms of A t and μ, we obtain, Setting Z 6 = 0, we can identify Z 1 v 2 with the mass of the observed (SM-like) Higgs boson (which may be either h or H depending on whether s β−α is close to 1 or 0, respectively). We can then numerically solve for tan β for given values of A t and μ. The values of the real positive tan β solutions of Z 6 = 0 obtained by using the one-loop approximate formula given in Eq. (20) are illustrated by the contour plots shown in the three left panels of Fig. 1, where each panel corresponds to a different solution of Z 6 = 0. Note that at every point in the ( μ, A t ) plane, the value of M S has been adjusted according to Eq. (15) such that the squared-mass of the SM-like Higgs boson in the alignment limit is given by Taking the three left panels together, one can immediately discern the regions of zero, one, two and three positive tan β solutions of Eq. (20), and their corresponding values. A more detailed discussion of these solutions will be presented in a separate paper [47]. It is instructive to obtain an approximate analytic expression for the value of the largest real positive tan β solution. Assuming | μ A t | tan β 1 the following approximate alignment condition, first written in Ref. [33], is obtained, t < 2 and | μ| is sufficiently large such that the numerator of Eq. (21) is negative. Keeping in mind that Eq. (21) was Taking the three panels on each side together, one can immediately discern the regions of zero, one, two and three values of tan β in which exact alignment is realized. In the overlaid blue regions we have (unstable) values of |X t /M S | ≥ 3 derived under the assumption that μ A t tan β 1, one easily verifies that the largest of the three roots of Eq. (20) shown in Fig. 1 always satisfies the stated conditions above. Another consequence of Eq. (21) is that by increasing the value of | μ A t | (in the region where 2 < A 2 t < 6), it is possible to lower the tan β value at which alignment occurs.
If | A t | 1, then Eq. (21) is no longer a good approximation. Returning to Eq. (20), we set A t = 0 and again assume that tan β 1. We can then solve approximately for tan β, For example, in the parameter regime where A t 0 and | μ| 1, we obtain tan β | μ|/ √ 6. The question of whether the light or the heavy CP-even Higgs boson possesses SM-like Higgs couplings in the alignment without decoupling regime depends on the relative size of Z 1 v 2 and Z 5 v 2 + M 2 A . Combining Eqs. (16) and (17), it follows that in the limit of exact alignment where Z 6 = 0, we identify Z 1 v 2 as the squared mass of the observed SM-like Higgs boson and We define a critical value of M 2 A , where Z 1 v 2 = (125 GeV) 2 and Z 5 v 2 is given by Eq. (23). Note further that the squared-mass of the non-SM-like CPeven Higgs boson in the exact alignment limit, M 2 A + Z 5 v 2 , is positive, which implies that the minimum value possible for the squared-mass of the CP-odd Higgs boson is That is, if Z 5 is sufficiently large and negative, then the minimal allowed value of M 2 A is non-zero and positive. We focus again on the parameter region in the ( μ, A t ) plane, and compute Z 5 from Eq. (23) using the value of tan β obtained from setting Z 6 = 0 in Eq. (20) Fig. 2, which are in one-to-one correspondence with the three left panels of Fig. 1.
As previously noted, the analysis above was based on approximate one-loop formulae given in Eqs. (15)- (17), where only the leading terms proportional to m 2 t h 2 t are included. In the exact alignment limit, we identify Z 1 v 2 given by Eq. (15) as the squared-mass of the observed SMlike Higgs boson. However, it is well known that Eq. (15) overestimates the value of the radiatively corrected Higgs mass. Remarkably, one can obtain a significantly more accurate result simply by including the leading two-loop radiative corrections proportional to α s m 2 t h 2 t . In Ref. [59], it was shown that the dominant part of these two-loop corrections can be obtained from the corresponding one-loop formulae with the following very simple two step prescription. First, we replace where m t (m t ) 165.6 GeV is the MS top quark mass [60], and the running top quark mass in the one-loop approximation is given by In our numerical analysis, we take α s = α s (m t (m t )) 0.10826. Second, when m 4 t multiplies that threshold corrections due to stop mixing (i.e., the one-loop terms proportional to X t and Y t ), then we make the replacement, where Note that the running top-quark mass evaluated at M S includes a threshold correction proportional to X t that enters at the scale of supersymmetry breaking. Here, we only keep the leading contribution to the threshold correction under the assumption that m t M S (a more precise formula can be found in Appendix B of Ref. [59]). The above two step prescription can now be applied to Eqs. (15)- (17), which yields a more accurate expression for the radiatively corrected Higgs mass and the condition for exact alignment without decoupling. Details of this analysis will be presented in a forthcoming work [47].
The end results are summarized below. We have derived analogous expressions to Eqs. (20) and (23)   introduce the following notation, where m t ≡ m t (m t ) is the MS top quark mass, and Then, the two-loop corrected condition for the exact alignment limit corresponding to Z 6 = 0 is given by, which supersedes Eq. (20), and the O(α s ) correction to Eq. (23) is given by, One can now define two-loop improved versions of M 2 A,c and M 2 A,m [cf. Eqs. (24) and (25)]. In the right panels of Figs. 1 and 2, we plot the twoloop improved versions of the corresponding one-loop results shown in the left panels. There are a few notable changes, which we now discuss. First, in our scan of the ( μ, A t ) plane, we have observed numerically that there is a new solution to the alignment condition [cf. Eq. (32)] that is unrelated to the solutions found in the one-loop analysis. However, this solution always corresponds to a value of |X t | > 3M S , which lies outside our region of interest. Henceforth, we simply discard this possibility. What remains are solutions that can be identified as the two-loop corrected versions of the oneloop results obtained above. The right panels of Fig. 1 exhibit the remaining real positive tan β solutions of Eq. (32).
We can now see the effects of including the leading O(α s h 2 t ) corrections. The regions where positive solutions to Eq. (32) exist shown in the right panels of Fig. 1 have shrunk somewhat as compared to the corresponding positive solutions to Eq. (20) shown in the left panels of Fig. 1. For example, only one positive solution for tan β exists for large μ and A t in the two-loop approximation, whereas three positive solutions exist in the one-loop approximation. Using the values of tan β found in the right panels of Fig. 1, one can now produce the corresponding two-loop corrected plots shown in the right panels of Fig. 2. The qualitative features of the one-loop and two-loop results are similar, after taking note of the slightly smaller regions in which positive solutions for tan β exist in the two-loop approximation.
One new feature of the two-loop approximation not yet emphasized is that we must now carefully define the input parameters μ and A t . In the above formulae and plots we interpret these parameters as MS parameters. However, it is often more convenient to re-express these parameters in terms of on-shell parameters. In Ref. [59], the following expression was obtained for the on-shell squark mixing parameter X OS t in terms of the MS squark mixing parameter X t , where only the leading O(α s ) corrections are kept, Since the on-shell and MS versions of μ are equal at this level of approximation, we also have A more detailed examination of the above results, when expressed in terms of the on-shell parameters, will be treated in Ref. [47]. The approximations employed in the section capture some of the most important radiative corrections relevant for analyzing the alignment limit of the MSSM. However, it is important to appreciate what has been left out. First, higherorder corrections beyond O(α s h 2 t ) are known to be relevant (see, e.g., Ref. [61]). In particular the O(h 4 t ) corrections are in magnitude roughly 20% of the O(α s h 2 t ) corrections, and enter with a different sign, thus leading to potentially nonnegligible corrections to the approximate two-loop results obtained above. On more general grounds, the analysis of this section ultimately corresponds to a renormalization of cos(β − α), which governs the tree-level couplings of the Higgs boson and its departure from the alignment limit. However, radiative corrections also contribute other effects that modify Higgs production cross sections and branching ratios. It is well-known that for M A M S , the effective low-energy theory below the scale M S is a general two Higgs doublet model with the most general Higgs-fermion Yukawa couplings. These include the so-called wrong-Higgs couplings of the MSSM [62,63], which ultimately are responsible for the b and τ corrections that can significantly modify the coupling of the Higgs boson to bottom quarks and tau leptons. 4 The implication of these couplings will be briefly reviewed in Sect. 2.3. In addition, integrating out heavy SUSY particles at the scale M S can generate higher dimensional operators that can also modify Higgs production cross sections and branching ratios [65]. None of these effects are accounted for in the analysis presented in this section.

Implications of the wrong-Higgs couplings
At tree-level, the Higgs-fermion Yukawa couplings follow the Type-II pattern [53,66] of the two-Higgs doublet model (2HDM), in which the hypercharge −1 Higgs doublet field H D couples exclusively to right-handed down-type fermions and the hypercharge +1 Higgs doublet field H U couples exclusively to right-handed up-type fermions. When radiative corrections are included, the so-called wrong-Higgs Yukawa couplings are induced by supersymmetrybreaking effects, in which H * D couples to right-handed uptype fermions and H * U couples to right-handed down-type fermions. We shall denote by M SUSY a generic scale that characterizes the size of supersymmetric mass parameters. In the limit where M Z , M A M SUSY , the radiatively-corrected Higgs-quark Yukawa couplings can be summarized by an effective Lagrangian, 5 which yields a modification of the tree-level relations between h t , h b and m t , m b as follows [64,[67][68][69][70][71][72][73][74]: The dominant contributions to b are tan β-enhanced, is suppressed, whereas h b does not decouple. This non-decoupling can be explained by the fact that h b arises from the radiatively-generated wrong-Higgs couplings. Below the scale M SUSY , the effective low energy theory is the 2HDM which contains the most general set of Higgs-fermion Yukawa couplings allowed by gauge invariance, and is no longer restricted to be of Type-II [62]. Similarly, δh t ∼ O(M 2 Z /M 2 SUSY ) is suppressed, whereas h t does not decouple. However, t is not tan β-enhanced and thus yields only small corrections to the Higgs boson couplings to fermions in the parameter regime of interest (i.e., where tan β > ∼ 1).
where mg is the gluino mass, mb 1,2 and mt 1,2 are the bottom and top squark masses, respectively, and smaller electroweak corrections have been ignored. The loop integral I (a 2 , b 2 , c 2 ) is given by Note that and I (a, a, a) = 1/(2a 2 ). Thus, in the limit in which all supersymmetric parameters appearing in Eq.
Expressions for the Higgs couplings to the third generation quarks can be found in Ref. [64]. In particular, the charged Higgs coupling to the third generation quarks is noteworthy. It is convenient to write the approximate one-loop corrected H − tb coupling in the following form, where P R,L ≡ 1 2 (1 ± γ 5 ). One of the important constraints on the MSSM Higgs sector is derived from the decay rate for b → sγ due to the presence of one loop diagrams involving a charged Higgs boson. At large tan β, it is important to incorporate SUSY corrections to the charged Higgs couplings to quarks 6 in the computation of BR(b → sγ ) [76]. Including the radiatively corrected H +t b and H − ts couplings using Eq. (43), suitably generalized to include intergenerational quark mixing, and taking tan β 1, after comparing the result obtained from the contribution of the charged Higgs loop in the MSSM, including the leading SUSY radiative corrections to the charged Higgs-fermion couplings, to the corresponding results of the 2HDM with Type-II Yukawa couplings. In Eq. (44), b is given by Eq. (39) and h t is given by [76] h t h t 2α s 3π μmg cos 2 θt I (ms L , mt 2 , mg) where m s L is the mass of the SUSY partner of the left-handed strange quark, θt is the t Lt R mixing angle [77], and I is defined in Eq. (40). Once again, the non-decoupling behavior of h t is evident in the limit in which all supersymmetric parameters appearing in Eq. (45) are of O(M SUSY ). As previously emphasized, the non-decoupling properties b and h t arise due to the wrong-Higgs Yukawa couplings, and are responsible for the significance of the deviation from Type-II behavior of the two-Higgs doublet sector of the MSSM.

Sampling of the parameter space
We sample the pM SS M 8 parameter space with uniformly distributed random values in the eight input parameters. Scans are performed separately for the light Higgs and heavy Higgs interpretation of the observed Higgs signal (see below for details) over the parameter ranges given in Table 1. Besides the scan parameters listed in Table 1, the remaining MSSM parameters are chosen as described in Sect. 2.1.
In both cases, we start with O(10 7 ) randomly sampled points in the ranges given in Table 1 and identify interesting regions where either h or H has a mass close to the observed signal at 125 GeV (i.e. we select points with M h/H ∈ [120,130] GeV) and the global χ 2 function is low, see Sect. 3.2 for details on how the global χ 2 function is evaluated. In a second step we perform dedicated smaller scans over more restricted parameter ranges in order to obtain high sampling densities in the interesting regions of the parameter space.
The choices of the parameter ranges for the light Higgs and heavy Higgs case differ in particular for M A and tan β, where the ranges in the heavy Higgs case are quite restricted. This is because M H ∼ 125 GeV can only be obtained in a rather small region of the parameter space, and a high sampling density in this region is desired. Furthermore, while we scan the third generation squark masses, Mq 3 , up to 5 TeV in the light Higgs case, we restrict Mq 3 to be at most 1.5 TeV in the heavy Higgs case. As mentioned before, the SM-like Higgs boson mass can be lifted to the observed value of ∼125 GeV by radiative corrections from either a large stop mass scale, M S , or from a large stop mixing parameter, X t . We consider a larger Mq 3 range in the light Higgs case in order to allow for solutions with small to moderate X t , μ and A t values. In contrast, in the heavy Higgs case, the SM-like properties can only be obtained in the alignment limit (without decoupling) which already requires large values of μ/M S and/or A t /M S (see Sect. 2.2), and we restrict ourselves to Mq 3 < 1.5 TeV in this case. Lastly, the choice of the scanning range in the Higgsino mass parameter, μ, differs in the two cases. In the light Higgs case we restrict |μ| ≤ 3Mq 3 , thus allowing μ 15 TeV for very large third generation squark masses Mq 3 ∼ 5 TeV. Parameter points with more extreme values of |μ/Mq 3 | beyond ∼3 often face severe constraints from vacuum stability requirements [78-86] (for a recent analysis see also Ref. [87]). Nevertheless, in the heavy Higgs case we include such more extreme values of |μ/Mq 3 | and do not impose a specific upper limit on this ratio. As we discussed in Sect. 2.2, |μ/M S | greatly influences the tan β value where the alignment limit occurs as well as the critical M A value that indicates the crossover of the light and heavy Higgs case mass hierarchies. As we will see, a large ratio μ/Mq 3 will be crucial to obtain an acceptable fit of the heavy Higgs to the observed Higgs signal. We will comment on the fit outcome in the case where the requirement |μ/Mq 3 | ≤ 3 is imposed.
In our scans we allow both signs of the Higgsino mass parameter μ. The sign of the SUSY contributions to the anomalous magnetic moment of the muon, (g − 2) μ , is given by the sign of μ, thus the negative μ branch is significantly disfavored in the light of this observable, receiving a χ 2 penalty like the SM or higher, depending on the mass scale of the relevant SUSY particles. In this work we will therefore present fit results where (g − 2) μ is either included or excluded from the global χ 2 function, see below for details.
In addition to the eight pM SS M 8 scan parameters we sample the top quark pole mass from a Gaussian distribution with m t = 173.34 ± 0.76 GeV [88], using a cutoff at ±2σ . Effects from other parametric uncertainties of SM quantities are estimated to be small and therefore neglected in this analysis.
We calculate the SUSY particle spectrum and the MSSM Higgs masses using FeynHiggs (version 2.11.2) 7 [61,[90][91][92][93][94] and estimate the remaining theoretical uncertainty (e.g. from unknown higher-order corrections) in the Higgs mass calculation to be 3 GeV [61]. Following Refs. [19,35], we demand that all points fulfill a Z-matrix criterion, in the light (heavy) Higgs case), in order to ensure a reliable and stable perturbative behavior in the calculation of propagator-type contributions in the MSSM Higgs sector. 8 In the Feynmandiagrammatic approach of FeynHiggs all model parameters (except for tan β, which is a DR parameter defined at the scale m t ) are defined in the on-shell (OS) renormalization scheme, which we adopt for the definition of our fit parameters [cf. (3)].

Observables
In our scan we take into account the following experimental measurements (we denote all experimental measurements with a hat, while unhatted quantities correspond to the model predictions of the respective quantity):

• Higgs boson mass:
We use the combined result from the ATLAS and CMS Higgs mass measurements [3], where we linearly combine the uncertainties. In our pM SS M 8 scans, where the measured Higgs mass corresponds to either the light Higgs or the heavy Higgs mass, we linearly add the theoretical Higgs mass uncertainty of 3 GeV. Thus, the total mass uncertainty in the MSSM case is σM H = 3.32 GeV.

• Higgs signal rates:
We employ the public code HiggsSignals-1.4.0 [17,95,96] to evaluate a χ 2 value, χ 2 HS , for the compatibility of the pM SS M 8 predictions with rate measurements in 85 different Higgs signal channels from the LHC experiments ATLAS and CMS, as well as the Tevatron experiments CDF and DØ. A detailed list of all Higgs rate observables is given in Appendix B. HiggsSignals takes into account the correlations among major systematic uncertainties, including the uncertainties of the integrated luminosity and the theoretical uncertainties for the cross section and branching ratio predictions for a SM Higgs boson. It furthermore takes into account a potential overlap of signals from nearby Higgs bosons by simply adding the signal rates if the mass difference of the Higgs bosons is less than the experimental mass resolution of the search channel. 9 This feature is of relevance for the heavy Higgs interpretation where all three neutral Higgs bosons can be within the mass range ∼(100-150) GeV. For instance, the mass resolution of the H SM → τ + τ − analyses is typically assumed to be ∼25 GeV, thus the τ + τ − signal rates of Higgs bosons within the above mass range will potentially be added by HiggsSignals. The Higgs production cross sections are evaluated, both in the MSSM and the SM, with the code FeynHiggs (version 2.11.2) [61,[90][91][92]. This includes an implementation of the SM cross sections of the LHC Higgs cross section working group (LHCHXSWG) [98-100] (using the gg → H cross section prediction from Refs. [101,102]). The MSSM Higgs production cross sections are calculated in the effective coupling approximation [103]. More details on the calculation of the production cross section in the various channels can be found in Refs. [103,104]. The Higgs decay widths are also calculated with FeynHiggs, including the full one-loop corrections for the Higgs decay to fermions and leading higher-order contributions [105,106].

• Low energy observables (LEOs):
We include the rare B meson decays B → X s γ (X s represents any hadronic system containing a strange quark), B s → μ + μ − and B + → τ + ν τ , which have small branching ratios in the SM, being either loop or helicity suppressed. In SUSY, however, they can be mediated through SUSY particles and/or charged Higgs bosons, which can 9 Interference effects can be incorporated in this context using the method developed in Ref. [97]. Since in our analysis non-negligible interference contributions only occur between the CP-even states h and H in the parameter regions where they are nearly mass-degenerate, we neglect these effects here.
give sizable contributions. Thus these observables feature a high sensitivity to physics beyond the Standard Model (BSM). 10 We evaluate the SM prediction -using a top mass value of m t = 173.34 GeV and a Higgs mass value of M SM H = 125.09 GeV as input -and MSSM predictions for these flavor observables with the public code SuperIso (version 3.5) [111][112][113]. These predictions are listed besides the latest (combinations of) experimental measurements in Table 2. For BR(B → X s γ ) we use the current world average of Ref. [114]. Here we assign an additional uncertainty on the MSSM prediction for BR(B → X s γ ) of ±0.15 × 10 −4 [115,116]. The process B s → μ + μ − was observed for the first time by LHCb and CMS [117] and recently also by ATLAS [118]. 11 For B + → τ + ν τ we use a Belle combination of measurements using hadronic and semi-leptonic tagging methods with a combined significance of 4.6σ [119,120]. For all observables in Table 2, theoretical and experimental uncertainties are combined linearly. The anomalous magnetic moment of the muon, a μ = 1 2 (g − 2) μ , comprises another very sensitive low energy probe of BSM physics [121][122][123]. The experimentally observed value exhibits a very persistent deviation from the Standard Model prediction at the level of 3-4σ [124][125][126]. We obtain the MSSM contribution to the anomalous magnetic moment of the muon from SuperIso, which includes the one-loop result [127] as well as leading two-loop contributions [128][129][130][131]. We cross-checked the SuperIso result with results from FeynHiggs and found good agreement. Besides the flavor observables and a μ , we also include the MSSM prediction of the W boson mass into our fit. The SM value for M W shows a 1.8σ deviation [132] from the latest experimental value [133,134]. Our MSSM evalua- 10 We do not include the B → D ( * ) τ −ν τ measurements [107][108][109], which show some tension with respect to the SM prediction. For an explanation within the MSSM (requiring a mass degeneracy between the lightest chargino and neutralino), see [110]. 11 The value used in this work does not include the ATLAS measurement yet. tion of M W follows Refs. [132,135] and includes, besides the most advanced SM calculation, the full SUSY oneloop contributions as well as leading SUSY two-loop contributions. The uncertainties from unknown higher-order corrections have been estimated to be around 4 MeV in the SM [136] and somewhat larger (∼ (4 − 9) MeV) in the MSSM [137,138], depending on the SUSY mass scale. The main parametric uncertainty on M W stems from the top quark mass and does not need to be included in our χ 2 evaluation because we vary m t within its 2σ uncertainty in the scan. The remaining parametric uncertainties from M Z and α had are ∼ 3 MeV. Combining these two sources of theoretical uncertainties linearly we estimate 10 MeV for the MSSM uncertainty and 7 MeV for the SM uncertainty.
From these observables and their predictions we evaluate for every parameter point in the scan the global χ 2 function As mentioned above, we denote all experimental measurements with a hat. Unhatted quantities correspond to the model predictions of the respective quantity. The sum over the low energy observables O i runs over the five observables mentioned above. The last term, −2 ln L limits , denotes the contribution from Higgs search limits at LEP and LHC, for which the likelihood information about the level of exclusion is available. Details will be given below in Sect. 3.3. The total number of degrees of freedom, ν, is given by the number of observables, n obs , minus the number of scan parameters, n para . We count every observable and constraint that contributes to the global χ 2 function, Eq. (47), to n obs , thus we have in total n obs = 93 if all observables are included in the fit. In the SM we have only one free parameter (n para = 1), namely the Higgs mass, whereas in both MSSM cases we have eight fit parameters (n para = 8).

• Exclusion limits from Higgs collider searches:
For every scan point we test the neutral and charged Higgs bosons against the exclusion limits from Higgs searches at the LEP, Tevatron and LHC experiments by employing the public computer code HiggsBounds-4.2.1 [23,[139][140][141][142]. HiggsBounds determines for each model parameter point the most sensitive exclusion limit, based on the expected exclusion limit given by the experiments. It then judges whether the parameter point is excluded at the 95% C.L. by comparing the signal prediction against the observed exclusion limit from the most sensitive analysis. In this way, the quoted C.L. of the limit is preserved even though many different Higgs analyses are considered at the same time.
Besides the hard cut imposed by testing the parameter points at the standard 95% C.L. limit, HiggsBounds enables us to obtain a likelihood value for the model exclusion by LEP Higgs searches [13], as well as by the CMS search for non-standard Higgs bosons decaying into τ lepton pairs [143,144] (see Ref. [23] for details).
While the LEP Higgs searches are only relevant in the heavy Higgs case, i.e. where the heavier Higgs state is the SM-like Higgs boson at 125 GeV, the CMS search yields important constraints in either case, and in particular at larger values of tan β. Each of these likelihoods, which we commonly denote as −2 ln L limits , approximately resembles a χ 2 contribution corresponding to one degree of freedom, and can therefore simply be added to the global χ 2 function, Eq. (47).

• Exclusion limits from SUSY collider searches:
Lower limits on sfermion and chargino masses from mostly model-independent direct searches at LEP are typically at the level of ∼100 GeV (summarized in the PDG review [145]) and are applied in our scan. We furthermore require the lightest supersymmetric particle (LSP) to be the lightest neutralino, however, we do not apply any dark matter relic density constraints. Exclusions from SUSY searches at Run 1 of the LHC are tested by employing the public computer code CheckMATE-1.2.2 [146], which includes all relevant 8 TeV SUSY analyses from ATLAS and CMS. However, due to the large computational effort and the large scan samples it is neither feasible nor relevant to test all scan points with CheckMATE. Therefore, in a postprocessing step, we select the most interesting parameter points, i.e. points with a χ 2 difference to the minimal χ 2 of less than ∼10, and only test the LHC SUSY search constraints on these points. 12 For each of these 12 We explicitly check that the point with minimal χ 2 is not excluded by CheckMATE, or, in case it is excluded, we select more points for points we evaluate the sparticle decay spectrum with SUSY-HIT-1.5 [147] and feed these into Herwig++ (version 2.7.1) [148,149] for Monte-Carlo generation of inclusive sparticle pair production and the evaluation of the leading-order production cross section. It is not computationally feasible to evaluate the NLO corrections to the leading-order cross section for each parameter point. Instead, we multiply the leading-order cross section by an estimated global k-factor of 1.5 to approximately account for these corrections. CheckMATE processes the MC events through the implemented ATLAS and CMS analyses and follows a similar statistical procedure as HiggsBounds: It first determines which analysis is the most sensitive one, based on the expected exclusion limit, and then applies the observed exclusion limit from only this search in order to judge whether the parameter point is excluded at the 95% C.L. or not.

Results
In this section we discuss the results of our numerical analysis. We first discuss the results for the best fit points in both the light and heavy Higgs case in order to give an impression on the overall fit quality. Recall that approximate alignment without decoupling relies on an approximate cancellation between tree-level and loop level contributions to the effective Higgs basis parameter Z 6 as discussed in Sect. 2.2. The extent of the tuning associated with the regions in the pM SS M 8 scan that exhibit approximate Higgs alignment without decoupling is discussed in Appendix A.

Best-fit points and fit quality
The minimal χ 2 over the number of degrees of freedom, ν, indicates the best achievable level of agreement with the observations within a specific model and defines our best-fit (BF) points. These are summarized in Table 3 for the SM and the two MSSM fits that are separately performed for the light Higgs interpretation (h) and heavy Higgs interpretation (H ). The results are given for fits to three different selections of observables (cf. Sect. 3.2): (i) only Higgs data (i.e. Higgs the CheckMATE test in order to retain the maximal χ 2 difference of ∼10 to the minimum. We furthermore provide the reduced χ 2 value, χ 2 ν ≡ χ 2 /ν, as well as the corresponding p value (assuming an idealized χ 2 probability distribution) for each scenario in Table 3.
In total, we have 92 (93) observables in the SM and light Higgs MSSM interpretation (heavy Higgs MSSM interpretation) contributing to the global χ 2 value (cf. Sect. 3.2): 85 Higgs signal rate measurements, one Higgs mass measurement, and one (two) Higgs exclusion observable(s) for the SM and light Higgs MSSM interpretation (heavy Higgs MSSM interpretation -here also the LEP exclusion bounds apply), as well as five low energy observables.
We treat the SM as a one-parameter model, where the free parameter is the Higgs mass, M H . Its best-fit value is mainly set by the Higgs mass measurement. The χ 2 contribution from the Higgs mass measurement is therefore negligible in the SM. In both MSSM cases we have eight free model parameters.
Taking into account only the Higgs data, the minimal χ 2 values found in all three cases are very similar, with the lowest value being found in the light Higgs case of the MSSM. However, accounting for the additional degrees of freedom in the two MSSM cases, the overall fit quality is slightly better in the SM. Nevertheless, all three scenarios give very high p values, indicating excellent agreement with the observations in Higgs searches in each case. 13 13 It should be noted that the Higgs signal rates of certain search channels are based on the same physical degrees of freedom of the model The picture does not change much when the three flavor observables BR(B → X s γ ), BR(B s → μ + μ − ) and BR(B u → τ ν τ ), as well as the W boson mass are included in the fit. In the SM the largest χ 2 contributions from these additional observables come from the W boson mass (χ 2 ∼ 1.5) and BR(B s → μ + μ − ) (χ 2 ∼ 0.6). Both MSSM best-fit points yield slightly better agreement with these observables in comparison to the SM.
Taking into account also the anomalous magnetic moment of the muon, a μ , as observable in the fit, the SM receives a large χ 2 penalty (χ 2 ∼ 11.3) and thus becomes disfavored with respect to the two MSSM interpretations. Both the light and heavy Higgs case of the MSSM are capable to accommodate the a μ measurement, receiving χ 2 contributions of only ∼ 0.3 and ∼ 1.8 from this observable at the best fit point, respectively. The various Higgs mass and LEO predictions and the respective pull, defined as (Ô i − O i )/σ i , for the BF points of the two MSSM interpretations are summarized in Table 4.
The parameters for the best-fit points in the light and heavy Higgs case are shown in Table 5. Naturally, the M A values differ significantly, with a value in the decoupling regime for the light Higgs case and a low value ∼ 170 GeV in the heavy Higgs case. Large mixing in the stop sector is required and thus their predictions cannot be varied independently. For example, in the MSSM, the predicted rates in h → Z Z * and h → W W * searches are directly related. The effective number of degrees of freedom is thus lower, see e.g. Ref. [150] for a detailed discussion and analysis. However, since the naive p value found in the study presented here is of O(> 50 %), no significant change of the conclusion on the validity of the model can be expected from a pseudo-data based study as performed in Ref. [150].  in the light Higgs case to yield M h ∼ 125 GeV. As explained above, large higher-order corrections are also required in the heavy Higgs case. For the best fit value in the latter case the trilinear coupling A t is small, while X t is still sizable because of the large contribution from the term μ/ tan β. We find large and positive values of μ for the best fit points in both scenarios. In the light Higgs case we find parameter points providing a very good fit in the entire positive μ range (as we will discuss below), whereas in the heavy Higgs case large values for μ/M S are crucial to achieve the approximate alignment limit. The scalar leptons of the first and second generation are relatively light to accommodate a μ , and also the preferred value for M 2 is relatively low.
A comparison of the Higgs signal rates -given in terms of signal strength modifiers μ, which are defined by the measured value of σ ×BR in the respective channel normalized to the SM prediction (see e.g. Ref. [95]) -between the predictions of the BF points in the light and heavy Higgs case and the Tevatron and LHC measurements is displayed in Fig. 3. The BF point predictions are closely centered around the SM prediction (μ = 1), especially in the light Higgs case, where deviations from the SM prediction are O(2 − 3%). For the heavy Higgs case BF point the deviations are slightly larger O(10%), however, these mostly appear in less accurately measured channels.

The light Higgs interpretation
We start our presentation of the fit results in the light Higgs interpretations of the MSSM with the predicted Higgs signal rates and their correlations for the preferred parameter points. We then give an overview of the preferred MSSM parameter regions. Here we separate the discussion between the full parameter space and a region with low CP-odd Higgs mass, M A ≤ 350 GeV. The latter selects most of the preferred parameter points that feature the limit of alignment without decoupling. The last two subsections provide dedicated discussions of the impact of the low energy observables and direct LHC SUSY searches on the fit.

Higgs signal rates
In Fig. 4 we show the χ 2 h = χ 2 h − χ 2 h,min distributions (the subscript 'h' refers to the light Higgs (h) interpretation), based on all observables, for four different Higgs signal rates, defined by Here where the upper and lower values are determined by χ 2 h ≤ 1, which approximately corresponds to the one-dimensional 68% C.L. interval. All central values lie very close to the SM prediction (R = 1).
The narrowest range is found for the Higgs signal rate R V h bb , i.e. for V h production with the Higgs boson decaying into bb. This is mainly because a direct variation of BR(h → bb) through a modification of the partial width (h → bb) also leads to a substantial change of the total Higgs decay width, h,tot , and thus affects significantly the branching ratios of all other Higgs decay modes. This global rescaling of the branching ratios of all decay modes except h → bb could only be compensated by an inverse rescaling in the rates of all Higgs production modes. However, this cannot be accomplished within the MSSM, and hence the modification of BR(h → bb) is severely constrained.
In Fig. 4 we can furthermore observe a spread of parameter points deviating substantially from the SM value, R = 1, within the χ 2 ≤ 4 interval (approximately corresponding to the 95% C.L. interval), e.g. R h γ γ 1.35 and R h τ τ 0.65. Such modifications can easily appear in the MSSM: An enhancement of the h → γ γ partial decay width can appear through the loop contribution of light charged SUSY particles such as light scalar tau leptons (staus) or charginos. 14 The observed reduction of R h τ τ for some points in Fig. 4 originates from a simultaneous (but small) suppression of the gluon fusion production cross section, σ (gg → h), and the decay rate BR(h → τ + τ − ) with respect to their SM predictions.
The two-dimensional correlations among the Higgs signal rates are shown in Fig. 5. Compared to Fig. 4 we introduced two new colors in order to indicate regions close to the minimum χ 2 . Points with χ 2 h < 2.30 (5.99) are highlighted in red (yellow), corresponding to points in a two-dimensional 68% (95%) C.L. region in the Gaussian limit. We shall denote these regions simply by favored/preferred (yellow) and most favored/preferred regions (red). The best fit point is indicated by a black star in the figures.
As already noticed in our previous analysis [19] the diphoton rate R h γ γ exhibits a strong correlation with R h V V and a strong anti-correlation with R V h bb . The latter arises through The color coding is the same as in Fig. 5 the strong influence of the h → bb partial width on the total Higgs decay width (see the discussion above). In contrast, the rate R h τ τ only shows mild correlations with R h V V , R h γ γ (not shown here) and R V h bb . The latter is easily understood from the fact that the same Higgs doublet couples to down-type quarks and leptons in the MSSM, thus, at tree-level, the light Higgs coupling to τ leptons and b quarks is affected in the same way (and enhanced at large tan β). However, we also find favoured points that feature a significant suppression of R h τ τ while R h bb shows no (or only a small) suppression. Differences between R h τ τ and R h bb arise from loop contributions, which modify the hτ + τ − and hbb couplings in a different manner, see e.g the discussion of b corrections in Ref. [19].

Parameter space
The distribution of preferred parameter points in the plane of the parameters M A and tan β, which determine the Higgs sector at lowest order, is shown in Fig. 6 (left). The bulk of the favored points is found at large CP-odd Higgs mass values, M A 350 GeV, i.e. in a region where the decoupling limit is already approximately realized. Large tan β values at moderate values of M A are disfavored by the non-observation of a signal in LHC H/A → τ + τ − searches, which is incorporated in our study by adding the exclusion likelihood from CMS (provided by HiggsBounds) to the global χ 2 function instead of applying the hard cut at the 95% C.L. (see Sect. 3.3) [23]. Thus, points excluded (e.g. at the 95% C.L.) by the CMS search alone can still appear as blue points here, but are unlikely to show up as yellow or red points.
While all most favored (red) points are found for M A 350 GeV, some preferred parameter points (yellow) are found at low CP-odd Higgs masses down to M A 170 GeV in a narrow range of tan β values ∼4-10. These points are far away from the decoupling limit, however, they turn out to feature an (approximate) realization of the limit of alignment without decoupling. As discussed in Sect. 2.2 these points must have either large values of μA t /M 2 S or, in the case where |A t /M S | is small, even larger values of |μ/M S | (beyond ∼3) in order to achieve the Higgs alignment limit at reasonably small values of tan β which are unexcluded by LHC H/A → τ + τ − searches. This is indeed the case, as can be seen in Fig. 6 (right), where we show the preferred parameter points in the (M A , μA t /M 2 S ) plane. All preferred parameter points with M A 350 GeV have μA t /M 2 S values of at least 3, but typically between 5 and 9. It should be noted that in the light Higgs scan we constrain both scan parameters μ and A t to be ≤ 3Mq 3 (see Table 1), limiting the possible values of μA t /M 2 S to 9. In the following, in order to study the parameter points close to the limit of alignment without decoupling, we apply a cut M A ≤ 350 GeV (denoted as low-M A selection) whenever relevant to isolate these points from the other parameter points.
The MSSM parameters from the stop sector, namely the stop mixing parameter, X t /M S , and the light stop mass, mt 1 , are shown in Fig. 7 for the full scan (left) and the low-M A selection (right). In the full scan we find preferred parameter points in both the positive and negative X t /M S branches near the value where the contribution to the Higgs mass from stop mixing is maximized, |X t /M S | ∼ 2. 15 Light stop masses, mt 1 , down to values 300 (400) GeV are possible in the positive (negative) X t /M S branch. 16 In the low-M A selection preferred parameter points are found only in the positive X t /M S branch for X t /M S 2. Here, the lowest light stop mass value in the preferred parameter region is found at around mt 1 ∼ 580 GeV.
In order to understand why the favored parameter points near the limit of alignment without decoupling are found only at X t /M S = A t /M S − (μ/M S ) cot β 2, we show the cor- 15 The highest M h values are reached for |X t /M S | ∼ 2 due to the inclusion of the higher-order corrections in the on-shell renormalization scheme, see Ref. [151] for details. 16 Note that we assumed universality of the left-and right-handed softbreaking stop mass parameter here. Lower light stop masses even below the top quark mass can be obtained while being consistent with the Higgs rates in the presence of a large mass splitting in the stop sector [152].  A t /M S 3, while the range in μ/M S is larger (roughly between 1.4 and 3). In the full scan (left) it can be seen that even values with A t /M S and μ/M S close to zero can yield a very good fit. Consequently, the quite large best fit value of μ, see also Table 5, should be regarded as accidental.
In the following we will analyze in detail which observables and constraints lead to this particular favored region of parameter space in the alignment limit. We recall the parametric dependence of the approximate one-loop alignment condition in the limit |μA t | tan β M 2 S , (21), In order to find viable solutions of this condition (while restricting ourselves to |A t |/M S ≤ 3 and |μ|/M S ≤ 3) we need to have 17 We illustrate the impact of the various constraints and observables on the possible solutions for the limit of alignment without decoupling in the (μ/M S , A t /M S ) plane in Fig. 9.
Here we plot only the parameter points in the low-M A selection with χ 2 ≤ 5.99 (approximately corresponding to the 95% C.L. region) based on different sets of observables: only Higgs mass and signal rates (top left), Higgs mass, signal rates and h/H/A → τ + τ − exclusion likelihood (upper right), all observables except a μ (lower left) and all observables (lower right). In color we indicate the tan β value of the parameter points (see legend). The upper left plot of Fig. 9 shows all possible regions where the alignment condition, Z 6 v 2 = 0, (cf. (17)), is approximately fulfilled. . It can clearly be seen that |A t /M S | < √ 6 is required in the (+, −) and (−, +) quadrants, consistent with (51), and that these parameter points tend to have larger tan β values than those in the (+, +) and (−, −) quadrants. Once the constraints from LHC H/A → τ + τ − searches are taken into account, as shown in the upper right plot of Fig. 9, the large tan β pointsand thus the (+, −) and (−, +) quadrants -become strongly disfavored. Small |A t | values in the (+, +) and (−, −) quadrants also require larger tan β and are thus equally disfavored. Adding also the flavor observables BR(B → X s γ ), BR(B s → μ + μ − ) and BR(B + → τ + ν τ ) as well as the W boson mass observable to the fit, the negative μ region (as well as the regions with μA t < 0) become mostly disfavored, as shown in the lower left plot of Fig. 9. Interestingly, this feature emerges already without including the observable a μ in the fit. The negative μ region is mostly disfavored by BR(B → X s γ ), while the negative μA t region is disfavored by BR(B s → μ + μ − ), as we will discuss in more detail in the next section. In addition, the negative μ region becomes strongly disfavored after adding a μ to the fit, as the sign of the SUSY contribution to a μ depends on the sign of μ. Thus μ needs to be positive in order to account for the currently observed discrepancy between measurement and theory prediction. This is shown in the lower right plot of Fig. 9 where only points with positive μ and A t at low tan β values remain, reproducing the distribution of favored points in the right plot of Fig. 8.

Impact of low energy observables
In this section we discuss the interplay of the Higgs observables and limits from Higgs searches with the low-energy observables, in particular the rare B decays, in the global fit. The light Higgs case features a very good fit to all lowenergy observables as we have already seen in Tables 3 and  4. For the most part we concentrate on the low-M A selection and study the low-energy observables for the parameter points close to the limit of alignment without decoupling. This is particularly interesting since the low M A value implies that all MSSM Higgs bosons, and in particular the charged Higgs boson, are relatively light, which can lead to large contributions to the B decays. In contrast, the MSSM (loop-)contributions to the anomalous magnetic moment of The leading contribution to the FCNC process b → sγ occurs in the SM via a W ± -t loop, allowing new-physics contributions to be of similar size. The branching ratio BR(B → X s γ ) can receive sizable positive contributions from an H ± -t loop if the charged Higgs boson mass is not too large. It was recently pointed out [153] that b → sγ excludes charged Higgs bosons with M H ± < 480 GeV at the 95% C.L. in a 2HDM with Type-II Yukawa couplings [53,66,154]. However, as discussed in Sect. 2.3, wrong-Higgs Yukawa couplings are induced radiatively by SUSY-breaking effects in the MSSM. These lead to important modifications of the charged Higgs couplings to up-and down-type quarks [cf. Eq. (43)], which in turn change the BR(B → X s γ ) prediction in the MSSM from the corresponding prediction in the Type-II 2HDM according to Eq. (44). In particular, if b [given by Eq. (39)] is positive, then the MSSM prediction for BR(B → X s γ ) is smaller than in the Type-II 2HDM. Furthermore, supersymmetric particles in the loop can also contribute to the b → sγ amplitude. For example, charginostop (χ ± -t) loops can contribute with either sign, depending on sign and magnitude of the parameters μ, M 2 and A t , and thus may partially cancel the effects of the H ± -t loop.
The impact of the LEOs on the global fit is illustrated best by studying the LEO predictions of the parameter points that are preferred before the LEOs are included in the fit. In the following we therefore focus on the (approximate) 95% C.L. preferred parameter points in the low-M A selection after the Higgs signal rates, Higgs mass and h/H/A → τ + τ − exclusion likelihood are included in the fit, i.e. the points in the upper right plot of Fig. 9. The left plot in Fig. 10 shows the charged Higgs mass dependence of BR(B → X s γ ). Generally we observe that for light charged Higgs values corre-   Fig. 10. The color coding of the displayed points is the same as in Fig. 9 sponding to the low-M A selection (slightly) too large predictions for BR(B → X s γ ) are obtained. For M H ± < ∼ 350 GeV we do not find any preferred points within the 1σ region of the experimental measurement, which is indicated by the green band. One can see that the BR(B → X s γ ) prediction tends to increase when going to smaller charged Higgs masses, M H ± . The two branches in this parameter region visible at low M A values correspond to μ > 0 (lower branch) and μ < 0 (upper branch), which can also be seen in the right plot in Fig. 10, where BR(B → X s γ ) is shown as a function of μ/M S (for M A < 350 GeV). The lower branch lies mostly within the 2σ region of the experimental measurement whereas most of the upper branch is inconsistent with the measurement at the 2σ level. This confirms our previous statement that negative μ is disfavored by B → X s γ .
The μ dependence of the BR(B → X s γ ) prediction is easily understood from the discussion in Sect. 2.3. In the  44)]. This b dependence is also directly shown in Fig. 11 for the parameter region with M A ≤ 350 GeV. The color coding furthermore illustrates that b ∝ tan β, i.e. the largest suppression of BR(B → X s γ ) is obtained for large, positive μ/M S and large tan β. Note that the chargino-stop contributions to BR(B → X s γ ) are found to be relatively small in this parameter region and thus play only a minor role in this discussion.
The decay B s → μ + μ − can be mediated in the MSSM by a neutral Higgs boson and receives loop corrections involving squarks, sleptons and electroweakinos. As can be seen in Fig. 12, in the limit of alignment without decoupling we find a very good fit to BR(B s → μ + μ − ), while for larger values of M H ± larger deviations are possible. In the right plot of Fig. 12, where BR(B s → μ + μ − ) is shown as a function of μA t /M 2 S (for M A < 350 GeV), one can see that the few points with negative μA t yield a value of BR(B s → μ + μ − ) close to the SM result, whereas the many points with positive μA t predict a smaller BR(B s → μ + μ − ) than in the SM, which are in better agreement with the current experimental central value. These negative corrections to BR(B s → μ + μ − ) become more sizable for larger values of tan β 5. It is interesting to note that predictions for BR(B s → μ + μ − ) that precisely match the experimental central value are possible over the whole range of M H ± values displayed in the left plot of Fig. 12.
The decay B + → τ + ν is helicity suppressed in the SM. The tree-level exchange of a charged Higgs boson constitutes the dominant MSSM contribution. These contributions can be sizable for large values of tan β and small M H ± . In the low-M A selection we have small tan β after including the h/H/A → τ + τ − exclusion likelihood and thus the contributions to BR(B + → τ + ν) are small. We have checked that the MSSM predictions are close to the SM value over the entire range of M H ± covered in our scan.
Finally, Fig. 13 shows the prediction for the SUSY contribution to the anomalous magnetic moment of the muon, δa μ , against the first and second generation slepton mass, mẽ ,μ for all scan points (left) and in the low-M A selection (right). 18 The MSSM contributions to a μ consist of smuon-neutralino (μ ± −χ 0 ) and sneutrino-chargino (ν−χ ± ) loops. Clearly δa μ is strongly correlated with the slepton mass, whereas the correlation with the chargino and neutralino masses (not shown) is less pronounced. In the full scan the favored region extends to large tan β values, implying that sizable contributions to δa μ are possible also for moderately large slepton masses. Parameter points within the 1σ region around the desired δa μ value are found all the way up to mẽ ,μ ∼ 1000 GeV for large tan β. In the low-M A selection, tan β is small ( 10) and therefore the sleptons need to be very light, 300 GeV, in order to give a large enough contribution to δa μ . This has interesting consequences for the direct SUSY searches, discussed below, where we find that a larger fraction of the preferred points in the alignment region is excluded compared to the decoupling region. This is because SUSY searches with multilepton final states targeting direct stop, gaugino or slepton production yield stronger exclusion if the 1st/2nd generation sleptons are light. 18 Note that in order to be able to separately analyze the corresponding effects in the fit to the Higgs rates and the fit to a μ , we have chosen to treat M˜ 1,2 and M˜ 3 as independent fit parameters in this work. Indeed, light staus can significantly influence the Higgs rates, in particular the γ γ rate [19,152,155].

Fig. 13
The SUSY contribution to the anomalous magnetic moment of the muon, δa μ , as a function of the (1st and 2nd generation) slepton mass, mẽ ,μ , for the favored points in the fit without taking into account the LEOs, for all scan points (left) and for the low-M A selection, M A < 350 GeV (right). The green line indicates the desired new physics contribution needed to achieve agreement with the observed deviation from the SM, and the hatched region corresponds to the 1σ experimental uncertainty. The color coding of the displayed points is the same as in Fig. 9

Impact of direct LHC SUSY searches
As discussed in Sect. 3.3, we have included all relevant 8 TeV SUSY searches in our analysis. 19 Overall, we find that the LHC SUSY searches impact the validity of the scan points in the light Higgs interpretation considerably. In total 38% of the 2.4×10 5 tested pM SS M 8 parameter points are excluded by the limits from LHC SUSY searches. 20 Restricting ourselves to the points within the (approximate) 95% C.L. favored region, LHC SUSY searches exclude 43% of these points. The impact of the direct LHC SUSY searches is mostly "orthogonal" to the impact of the Higgs constraints in the sense that the direct SUSY searches have no impact on the χ 2 profile. This observation can already be inferred from the light gray points of Fig. 4.
This point is further exemplified in the left plots of Fig. 14 in two-dimensional mass planes of the lightest stop, stau, selectron/smuon and neutralino, which are the most important SUSY particles in this context. The impact of the LHC searches for SUSY particles has been tested with CheckMATE for all displayed points. The points excluded by direct SUSY search limits are plotted in pale colors above the points allowed by the direct SUSY searches, which are shown in bright colors. The excluded points are found with a similar distribution as the allowed points in all projections of these mass parameters. The reduced density of excluded points for low values of the lightest neutralino mass, mχ0 1 , is explained by a reduction of the signal acceptance in most of the searches applied (thett → 1 + N jets + E miss T search [164] being a notable exception). This is caused by 19 CheckMATE-1.2.2 does not yet include the first 13 TeV analyses of ATLAS and CMS. However so far the 13 TeV exclusion limits are comparable (or still weaker) than the 8 TeV limits. 20 Recall that only points with a χ 2 ≤ 10 with respect to the BF point are fed into the CheckMATE analysis. the fact that often mχ0 due to the assumption of (6). For small neutralino masses the signal acceptance ofχ ± /χ 0 → (W ± /Z )χ 0 1 final states reduces as the phase space of the decay decreases, mχ0 /χ ± − mχ0 Fig. 11 of Ref. [168]), yielding reduced final state activity. Furthermore, the missing energy is smaller for small neutralino masses, affecting in particular the signal acceptance of searches for hadronic final states. Nevertheless, we do not find any specific effect of this feature on the allowed Higgs sector phenomenology.
Although the impact of the Higgs data seems orthogonal to the impact of the LHC SUSY searches, the impact of some of the included LEOs, in particular a μ , is not. The latter favors light first/second generation sleptons and light gauginos, which can be probed by searches at the LHC. As a result, we find a significant amount of favored points being excluded by multilepton searches for direct slepton, electroweak gaugino or stop pair production. As already mentioned above, this effect is particularly prominent in the low-M A region, where the requirement of light 1st/2nd generation sleptons to accommodate δa μ is much stricter than in the decoupling region (cf. Fig. 13).
For the SUSY exclusion using CheckMATE, only 8 TeV searches and LO cross sections are used with a global kfactor of 1.5 (see Sect. 3.3), due to their availability in the applied computer codes (as mentioned above, the version of CheckMATE used for our analysis does not yet include the first limits from ATLAS and CMS searches at 13 TeV). The results shown in the left plots of Fig. 14 imply that possible limits from upcoming searches at 13 TeV are not expected to significantly alter the Higgs signal interpretation in the pM SS M 8: Due to the described "orthogonality" of the SUSY searches to the Higgs rate constraints, a refinement of the k-factor calculation or a further strengthening of the limits, both in terms of rate constraints and in terms of mass reach, are not expected to change the parameter ranges pre- In the three left panels, we compare the CheckMATE allowed (bright colored) and excluded (pale colored) points for all scan points (blue), the (approximate) 95% C.L. (yellow) and 68% C.L. (red) preferred points (plotted in this order, with excluded points on top of the allowed points in each step); In the three right panels, we exhibit LHC analyses that yield the exclusion (in the order of plotting): (i) N jets + E miss T searches (blue circles) [156,157], (ii) hadronictt searches with b-jets (red squares) [158][159][160][161], (iii) χ ± 1 χ 0 2 → 3 + E miss T searches (green up-triangles) [162], (iv)tt → 2 + E miss T searches (yellow downtriangle) [163], (v)tt → 1 + N jets + E miss T searches (orange stars) [164], (vi) χ ± 1 χ 0 2 ,˜ ˜ → + − + E miss T searches (magenta diamonds) [165], (vii) 2 + N jets + E miss T (gray plus) [166,167] ferred by the Higgs rate observables significantly. Instead, it can be expected that strengthened SUSY limits will only decrease the point density, thus indicating that parameters need an increasingly refined tuning in order to still find a good fit for the Higgs mass and the Higgs rates.
The sensitivity of the various SUSY searches in the different kinematical regions of the parameter space is shown in the right plots of Fig. 14 Fig. 14 it can be seen that the sensitivity of LHC searches for direct electroweak gaugino and slepton production with dilepton final states (magenta diamonds) is centered at low mass parameters for sleptons, staus and neutralinos. Dilepton searches for stop pairs (yellow down-triangles) and multilepton searches for electroweak gauginos (green up-triangles) also provide important constraints at low stau and slepton masses, and in particular remain sensitive at larger neutralino masses.

The heavy Higgs interpretation
We now consider the more exotic MSSM interpretation where the heavy CP-even Higgs boson is identified as the observed Higgs boson at 125 GeV. We recall that this interpretation has a similarly good fit quality as the light Higgs interpretation (see Table 3). As in the light Higgs case we first show the predictions for the Higgs signal rates and their correlations and then discuss the preferred parameter region. After commenting on the impact of the low energy observables and LHC SUSY searches on the fit we conclude this section by discussing the discovery prospects for the other neutral and charged Higgs states in the heavy Higgs interpretation.

Higgs signal rates
We show the χ 2 profiles in the most important Higgs signal rates (cf. (48)) in Fig. 15. From the sparseness of the allowed points (blue) in these distributions it is evident that the heavy Higgs interpretation is quite constrained. Moreover, the actual χ 2 minimum of all scan points lies signif-icantly deeper than the determined best-fit point, however, all points with a lower total χ 2 than the best-fit point are excluded by either direct Higgs search limits, SUSY search limits, or other more technical requirements such as the Zmatrix criterion (see Sect. 3.1).
We find the preferred ranges for the Higgs signal rates to be The Higgs rates R P(H ) X X are defined as in the light Higgs case, see (48), but with h ↔ H . Due to the sparseness of points in the vicinity of the χ 2 minimum, these ranges should be taken only as indicative results for the actual 68% C.L. range (and therefore we also refrain from giving the central values here). Nevertheless, these distributions indicate on the one hand rather good agreement with the SM prediction (R H X X = 1) for the H → V V and V H → V bb channels, and on the other hand a modest suppression of the H → γ γ and H → τ τ signal rates with respect to the SM prediction. These tendencies appear not only for the most favored points (with χ 2 ≤ 1) but also for the bulk of allowed points with larger χ 2 values. They may thus be tested in the current and upcoming LHC runs. In Ref. [19] we analyzed two mechanisms for modifying R γ γ . It can be enhanced (suppressed) either via a suppression (enhancement) of the total width, which is dominated by H → bb, or via a directly enhancement (suppression) of the H → γ γ decay width through light charged particles, e.g. the lightest stau. Since we observe R V H bb ∼ 1 here, the suppression of R H γ γ is due to a direct suppression of the H → γ γ width induced by SUSY loops in the H → γ γ decay amplitude. Indeed we find for the most favored points (H → γ γ )/ (H → γ γ ) SM ∼ 0.80 − 0.95, from loops of a light charged Higgs boson and a moderately light stop, mt 1 ∼ (350 − 650) GeV, where the stop mixing parameter is close to X t ∼ −1.5M S (see Sect. 4.3.2) [152].
We show some of the correlations of the four Higgs signal rates in Fig. 16 Fig. 5). We note that the (approximate) 95% C.L. region extends over smaller values for the H → τ + τ − rate than in the light Higgs interpretation, down to values of R H τ τ 0.5.

Parameter space
We show the fit results for the heavy Higgs interpretation in Fig. 17 in the (M A , tan β) plane (left) and the (μ/M S , A t /M S ) plane (right). The preferred parameter points expand over only a narrow range in the parameters determining the Higgs sector at lowest order, M A ∼ (140, 185) GeV and tan β ∼ 6 − 11. Compared to our previous results [19], where we found smaller values M A ∼(110-140) GeV being preferred, the favored parameter region has shifted towards larger M A values, caused by several reasons. Firstly, at small values M A 150 GeV the CP-odd Higgs boson A potentially contributes to the predicted signal rate at 125 GeV in the τ + τ − channel. 21 In that case, 21 HiggsSignals automatically adds the signal rates of Higgs bosons that overlap within the combined experimental and theoretical  Fig. 16 the predicted signal rate would tend to be higher than the total observed τ + τ − rate, resulting in a larger χ 2 from HiggsSignals. In Ref. [19] we also took a possible signal overlap of H and A in the τ + τ − channel into account; the measurements at that time, however, were not accurate enough to notably affect the fit outcome. Secondly, parameter points with charged Higgs masses M H + below 160 GeV are strongly constrained by exclusion limits from LHC searches for a charged Higgs boson in top quark decays, t → H + b, with successive decay to τ leptons, H + → τ + ν τ [44,45]. At tree-level, the CP-odd and charged Higgs masses are related as M 2 H ± ,tree = M 2 A + M 2 W , thus, these constraints apply in particular at low values M A 140 GeV. In Ref. [19] we found good discovery prospects for the heavy Higgs case in t → H + b → (τ ν τ )b searches. Based on the most recent limits from such searches performed by ATLAS and CMS [44,45] the favored parameter regions of Ref. [19] are now excluded and the new preferred parameter space has moved towards larger M A values in the light of the updated limits. Thirdly, another reason for disfavoring M A values below ∼ 150 GeV is the prediction of somewhat too large values of BR(B → X s γ ), as will be discussed below.
The distribution of preferred parameter points in the (μ/M S , A t /M S ) plane, Fig. 17 (right), singles out values of μ ∼ (6 − 9)M S and A t from −M S to zero. This is due to an interplay of various observables and constraints, as we will outline in the following. Agreement of the heavy Higgs boson H with the LHC Higgs rate measurements requires firstly that M H ∼ 125 GeV and m h < M H (corresponding to M A < M A,c , see Sect. 2.2) and secondly that the alignment condition, Z 6 v 2 = 0, is approximately fulfilled. Furthermore, the LHC H/A → τ + τ − constraints impose Footnote 21 Continued mass uncertainties. For most Higgs channels with τ + τ − final states, the experimental mass resolution is assumed to be 20% · m H ≈ 25 GeV, thus the signals of a 125 GeV heavy Higgs H and a 150 GeV CP-odd Higgs A would be added.

Fig. 18
Preferred parameter regions in the (X t /M S , mt 1 ) plane in the heavy Higgs case. The color coding is the same as in Fig. 16 that this alignment must occur at not too large tan β values. As we analyzed in detail in Sect. 2.2, these requirements single out small regions of the (μ, A t ) plane, namely either μ/M S ∼ A t /M S ≈ ±(3 − 4) or |μ|/M S ≈ 4 − 9 (or even larger) and negative A t with |A t |/M S 1. Quite generically, in the heavy Higgs interpretation we have a slightly too high prediction of BR(B → X s γ ) due to the light charged Higgs boson. This discrepancy, however, decreases for very large positive values of μ/M S (similar to the alignment limit in the light Higgs case, cf. Fig. 10). Combining these arguments we find that large positive μ values, together with rather small negative A t are favored.
In order to avoid problems with vacuum instability, we applied a cut of |μ|/M S < 3 in the light Higgs case. Obviously, all our favored points in the heavy Higgs case would be cut away by such a constraint. 22   The color coding of the displayed points is the same as in Fig. 9 here that this cut only provides an approximate limit. Testing if any of our favored points in the heavy Higgs case are still allowed by vacuum stability would require a more thorough analysis, which is beyond the scope of this paper. The fact that A t and μ are predicted to be confined to a narrow range in the heavy Higgs interpretation also results in a definite prediction for X t and a quite small range for the light stop mass, as shown in Fig. 18. The favored region has X t ∼ −1.5 M S , and the light stop mass is found between 350 and 650 GeV -a prediction that can be tested at the upcoming LHC stop searches. However, note that the upper limit on the light stop mass is a consequence of our restricted scan range, μ ≤ 5 TeV, and the fact that large μ/M S values are favored.

Impact of low energy observables
The discussion of the low energy observables in the heavy Higgs case closely follows the corresponding discussion in the alignment limit of the light Higgs case, see Sect. 4.2.3, as the parameters most relevant for the low energy observables are similar, i.e. we again have small M A , M H ± and small tan β. Figure 19 shows BR(B → X s γ ) as a function of M H ± (left plot) and μ/M S (right plot) for the favored points of the fit to the Higgs signal rates and Higgs mass (including the exclusion likelihood from LEP and LHC searches, but without LEOs). As discussed above the parameter points with μ/M S ≈ 3 − 4 feature A t ∼ μ, whereas the points with μ/M S 6 correspond to small and negative A t values. By comparing the two plots in Fig. 19 we can also see that the parameter points with μ/M S ∼ A t /M S ≈ 3 − 4 have relatively light charged Higgs masses, M H ± < m t . Again we find the expected increase of BR(B → X s γ ) when going to small charged Higgs masses. Overall, we observe surprisingly good agreement with the experimental measurement.  Fig. 19. The color coding of the displayed points is the same as in Fig. 9 Even parameter points with a charged Higgs mass below the top quark mass lie mostly within the 2σ region of the experimental measurement. As in the light Higgs case (cf. Fig. 10) we again find that large positive μ values give the best agreement with the measurement (within the 1σ region or just above). However, whereas it takes μ/M S > ∼ 3 in order that the BR(B → X s γ ) prediction be within the 1σ region of the measurement in the light Higgs case, here we need much larger values μ/M S 6 to achieve this. Again, this can be understood from the sizable b corrections to the charged Higgs coupling to top and bottom quarks, as discussed in Sect. 2.3, and shown in Fig. 20, which impressively resembles the relation in Eq. (44).
The BR(B s → μ + μ − ) prediction is shown in the left plot of Fig. 21 in dependence of μA t /M 2 S . We find that the few points with positive μA t are in excellent agreement with BR(B s → μ + μ − ), whereas for negative μA t the predictions are slightly too high, in particular for μA t /M 2 S −5. Again this can be compared to Fig. 12 where we found a   Fig. 22 The SUSY contribution to the anomalous magnetic moment of the muon, δa μ , as a function of the (1st and 2nd generation) slepton mass, mẽ ,μ , for the favored points in the heavy Higgs case in the fit without taking into account the LEOs. The green line and hatched region show the desired new physics contribution needed to achieve agreement with the observed deviation from the SM and the 1σ experimental uncertainty. The colors of the displayed points is the same as in Fig. 9 similar result (even though there the typical μA t values of the favored points were quite different). Again we find only very small SUSY corrections to BR(B u → τ ν τ ) (despite the light charged Higgs) due to the small tan β values, as one can see in the right plot of Fig. 21.
The prediction for the SUSY contribution to the anomalous magnetic moment of the muon, δa μ , as a function of the slepton mass are shown in Fig. 22. Again we see a very strong correlation, i.e. the value for δa μ is almost entirely determined by the slepton mass. Light slepton masses 300 GeV are favored in the heavy Higgs case.

Impact of direct LHC SUSY searches
The comparison of the impact of the direct SUSY searches with the impact of the other constraints qualitatively agrees between the light Higgs interpretation (see Sect. 4.3.4) and the heavy Higgs interpretation. Also for the heavy Higgs interpretation, the SUSY search limits are thinning the parameter space in an "orthogonal" way to the Higgs observables and limits from the Higgs searches. We again show the distribution of CheckMATE allowed and excluded points in terms of the most relevant mass parameters in Fig. 23  (left), and indicate the corresponding experimental search that yields the exclusion in Fig. 23 (right). Similar to the light Higgs case, the CheckMATE exclusion is correlated with low-energy constraints, particularly the constraints from the anomalous magnetic moment of the muon. As discussed in the previous section, a μ strongly favors low slepton masses. Hence we observe a large fraction of favored parameter points excluded by multilepton searches for gaugino pair production (green up-triangles) as long as M 2 − mẽ ,μ 20 GeV is fulfilled. As we have relatively light stops in the preferred region of parameter space, constraints from dilepton searches (yellow down-triangles) as well as hadronic searches with b-jets (red squares) targeting stop pair production are also relevant. Overall, LHC SUSY searches exclude only ∼5% of all tested parameter points (i.e. points with χ 2 ≤ 10), but ∼65% of the favored ( χ 2 < 5.99) parameter points. This illustrates again the complementarity between the impact of a μ and the constraints from SUSY searches sensitive to light sleptons.

Phenomenology of the other Higgs states
Here we discuss the prospects for the discovery of the other Higgs bosons at the LHC in the heavy Higgs interpretation. We start with the phenomenology of the charged Higgs boson, which is a crucial test of this scenario. In fact, the previous benchmark scenario for the heavy Higgs interpretation [19] has been excluded with limits that have meanwhile been obtained from searches for charged Higgs bosons [44,45].  [156,157], (ii) hadronictt searches with b-jets (red squares) [158][159][160][161], (iii) χ ± 1 χ 0 2 → 3 + E miss T searches (green up-triangles) [162], (iv)tt → 2 + E miss T searches (yellow down-triangle) [163], (v) tt → 1 + N jets + E miss T searches (orange stars) [164], (vi) [165], (vii) 2 + N jets + E miss T (gray plus) [166,167] In Fig. 24 we show the rate for the main production and decay channel of a light charged Higgs boson with M H ± < m t , BR(t → H ± b) × BR(H ± → τ ν τ ), as a function of M H ± . The impact of the current limits from charged Higgs searches in this channel [44,45] can be seen by the gray area in Fig. 24, cutting out the region with M H ± < 160 GeV and Only a few favored points (and none of the most favored points) have M H ± < m t and are therefore displayed in Fig. 24. As one can see, these points have charged Higgs masses close to the top quark mass and thus BR(t → H ± b) × BR(H ± → τ ν τ ) is strongly suppressed due to the limited phase space of the top quark decay. Also the decay BR(H ± → τ ν τ ) is suppressed by the competing decay H ± → hW ± , which is open for most of these favored points (we will discuss the m h range of the favored points below). Consequently it is very difficult to detect charged Higgs bosons in this mass range in the t → H + b → (τ + ν τ )b channel at the LHC.
Many of our favored and most favored points have M H ± > m t (and are thus not visible in Fig. 24). Charged Higgs Fig. 24 Branching ratio of the top quark decay into a charged Higgs boson and a bottom quark, with the successive decay of the charged Higgs boson into a tau lepton and neutrino, in the heavy Higgs case. The color coding is the same as in Fig. 16 bosons with masses above the top quark mass are searched for in the pp → t H ± production channel with H ± → τ ν τ [44,45,169] or H ± → tb [45,170]. These searches, although concentrating on the charged Higgs mass region that is relevant for the heavy Higgs interpretation, are not yet sensitive to constrain the favored parameter space. However, they will become more sensitive with increasing integrated luminosity. Furthermore, we emphasize again that the decay H ± → hW ± is possible and unsuppressed in large parts of the parameter space, but currently not directly searched for at the LHC. In Sect. 4.4 we will present specific benchmark scenarios, inspired by our best-fit point in the heavy Higgs case, that can be employed to study the sensitivity of these searches.
We will now turn to the discussion of the phenomenology of the light CP-even Higgs boson, h, in the preferred parameter region in the heavy Higgs case. The light CP-even Higgs boson has a mass in the range  GeV and a strongly reduced coupling to vector bosons. This is shown in the top left plot of Fig. 25, where the squared coupling g 2 hV V is displayed, normalized to the corresponding coupling in the SM with the same value of the Higgs boson mass. One can see that the squared coupling is reduced by a factor of 10 3 or more with respect to the SM, as the heavy CP even Higgs boson H in this scenario acquires the coupling to vector bosons with approximately SM Higgs strength. This results in a strongly reduced cross section for the LEP Higgs-Strahlung process, e + e − → Zh. Consequently, the light Higgs boson in this case would have escaped detection in corresponding LEP Higgs searches. The limits from the Higgs searches at LEP occur for higher values of the relative squared coupling g 2

hV V
and are not visible in this plot.
The reduced light Higgs coupling to vector bosons furthermore leads to a reduced rate of the h → γ γ decay, which happens through a W -boson loop (amongst other contribut-ing diagrams). In contrast, the light Higgs coupling to gluons is up to ten times stronger than the SM Higgs boson coupling at very low light Higgs masses, M h , as shown in the center left plot of Fig. 25, where the (SM normalized) squared light Higgs-gluon coupling, g 2 hgg , is shown in dependence of M h . This results in an abundant production of the light CP-even Higgs boson via gluon fusion. The resulting LHC cross section for gg → h (at 8 TeV) with subsequent decay h → γ γ is shown in the bottom left plot of Fig. 25. Limits from LHC searches in this channel [171] have been taken into account in our analysis (using HiggsBounds). Their effect can be seen in this plot as the gray excluded region above ∼ 0.1 pb. Clearly these searches are currently very far from being sensitive to detect the light Higgs boson in this scenario.
The light Higgs boson predominantly decays to bottom quarks (∼(70-80)% of the time) or tau leptons (∼(15-30)% of the time), as shown in the top and middle right plots in Fig. 25, respectively. Given the mass range 20 ≤ M h ≤ 90 GeV we expect direct LHC searches in both channels to be rather challenging, given the huge SM background and the difficulty to trigger the events. The majority of the favored points have M h > M H /2, however we also observe some favored points with lower M h (down to 20 GeV) for which the Higgs-to-Higgs decay channel H → hh is kinematically open. If this decay rate is sizable the rates of heavy Higgs H decays to SM particle final states can be significantly affected. We find that for our preferred points the branching ratio of this decay is at most ∼20%, but for most points 10%, as shown in the bottom right plot of Fig. 25. Branching fractions of this size lead to only moderate changes of the Higgs decay rates to SM particles. Consequently we still find a good fit in this scenario even for the case M h < M H /2 (the allowed points for M h > M H /2 with vanishing branching ratio of H → hh are not visible in the plot).
The CP-odd Higgs boson has a mass between 140 GeV and 185 GeV, as shown already in Fig. 17, for tan β in the range between 6 and 11. Consequently, in this scenario it should signify itself with higher luminosity at the LHC in A → τ τ searches. At the lower end of its mass range it might be visible possibly as an enhanced decay rate of the SM-like Higgs boson due to the limited mass resolution. At the higher end of its mass range it would appear as a new resonance decaying to τ + τ − .

Updated low−M H benchmark scenarios
In the previous section we demonstrated that the heavy Higgs interpretation in which the heavy CP-even Higgs boson is identified with the observed signal [7,[34][35][36][37] is still a viable scenario. While the original low-M H benchmark scenario [43] has meanwhile been ruled out by ATLAS and CMS via the search for a light charged Higgs boson in top quark decays [44,45], we showed that viable realizations of the heavy Higgs interpretation exist outside the parameter region where these searches are sensitive. In this section we define new versions of the low−M H benchmark scenario that are valid after taking into account all current experimental constraints. They are inspired by the best-fit point found in our global analysis of the heavy Higgs case, see Table 5, however, we slightly increased the stop mass scale (while roughly retaining the preferred μ to M S ratio) in order to evade potential exclusion limits on the stop mass from the upcoming 13 TeV LHC results. These new scenarios could provide a useful benchmark for the ongoing light charged Higgs boson searches in the MSSM. In particular, they exhibit the not-yet-sought-for MSSM decay signature H + → W + h and, in some parameter regions, even the decay H + → W + H .
We define three different scenarios, which we call low−M alt H in order to distinguish them from the previous low-M H benchmark scenario [43]. The first two low−M alt H scenarios follow the original idea in Ref. [43] and fix a "heavy" Higgs boson mass, in this case the charged Higgs boson mass, M H ± , to a certain value, whereas μ and tan β are taken as free parameters. We suggest two variants given by different choices of M H ± below (low−M alt− H ) and above (low−M alt+ H ) the top quark mass, m t . The second scenario fixes μ and explores the (M H ± , tan β) plane. This scenario may be utilized for charged Higgs boson searches in the relevant mass range. The parameters of these benchmark scenarios are given in Table 6.
In the following we discuss the compatibility of these benchmark planes with the current experimental constraints.  unphysical  [13], LHC H/A → τ + τ − searches [143,144] and LHC t → H + b → (τ ν)b searches [44,45], respectively. The green area indicates parameter regions that are compatible with the Higgs signal (at ∼ 95% C.L., see text for details), unphysical regions are displayed in gray (see text). In the two left panels, contour lines indicate the Higgs masses M h and M H (in GeV). In the two right panels, contours lines indicate the charged Higgs branching ratios, as well as the branching ratio for the top quark decay t → H + b (upper row) or the 13 TeV LHC cross section for charged Higgs production in association with a top quark, σ LHC13 (gb → t H + ) (in fb) [172,173] [44,45], LEP light Higgs searches [13] (mostly from the e + e − → Zh → Z (bb) channel) and LHC H/A → τ + τ − searches [143,144], respectively. For the latter two we again employ the χ 2 implementation of these results in HiggsBounds, see Sect. 3.3, and define the 95% C.L. excluded region by χ 2 ≥ 6.0 (given the two free parameters of the model). Moreover, we indicate the regions compatible with the Higgs signal based on the total χ 2 constructed from Higgs signal rates, Higgs mass and the two exclusion likelihoods from LEP light Higgs searches and LHC H/A → τ + τ − searches. The green area indicates where the p value -estimated from this total χ 2 value under the assumption of independent and Gaussian observables -is above 5%.  In the low−M alt− H scenario the branching fraction for the top quark decay into a charged Higgs boson, t → H + b, ranges between ∼0.1 and 0.25%. The charged Higgs H + successively decays either to τ + ν τ or W + h, where the branching ratios are highly dependent on the kinematical phase space of the H + → W + h decay, and thus on the light Higgs mass, M h . Either decay can be completely dominating while the other is suppressed, depending on the parameter space. The charged Higgs boson decay H + → cs is negligible.
In the low−M alt+ H scenario, the charged Higgs boson is predominantly produced in association with a top quark via gb → t H ± . The charged Higgs boson branching fractions in the low−M alt+ H scenario are very similar to the low−M alt− H scenario, with the exception that the decay H + → tb is present. However, its decay rate amounts to at most ∼20% (at the low end of the μ range). An interesting situation in this benchmark scenario occurs in the region around μ ∼ (5.0-5.5) TeV and tan β ∼ 6-7, where both the light and heavy Higgs boson have masses between 120 and 130 GeV. In our analysis, part of this region is even compatible with the Higgs signal and at the same time not directly excluded by limits from Higgs searches, i.e. this part of the parameter space gives rise to two CP-even Higgs bosons close to 125 GeV. A dedicated experimental analysis of such a sce-nario, taking into account also interference effects between the two nearly mass-degenerate Higgs bosons, see Ref. [97], would be desirable.
The low−M alt v H benchmark scenario is illustrated in Fig. 27. Here it can nicely be seen that the limit of alignment without decoupling occurs roughly at tan β ∼ 7, as the green area is centered around this value. The charged Higgs phenomenology is quite rich: at lower charged Higgs masses, M H ± 180 GeV, the decay modes H + → τ + ν τ and H + → W + h completely dominate. At larger M H ± , the decay mode H + → tb and even H + → W + H become non-negligible, albeit they remain small for most of the unexcluded parameter space.
The best prospects for exploring these scenarios may still be via LHC searches for H/A → τ + τ − . However, these benchmark scenarios will hopefully also provide useful guidance for upcoming charged Higgs searches and in particular motivate searches for new charged Higgs signatures such as H ± → W ± h, with the light Higgs h decaying into bottom quark or tau lepton pairs, h → bb, τ + τ − . Excluding these scenarios, which appear to be cornered from all sides, with the upcoming searches at the LHC would strongly restrict the heavy Higgs interpretation in the MSSM, approaching an exclusion of this interesting possibility, whose phenomenology drastically differs from the most commonly considered light Higgs scenario.

Conclusions and outlook
We have analyzed the compatibility of the phenomenological Minimal Supersymmetric Standard Model (pMSSM) with the SM-like Higgs boson with a mass ∼125 GeV as measured by ATLAS and CMS. We performed a parameter scan of the pMSSM with the eight most relevant parameters var-ied freely (pM SS M 8): the CP-odd Higgs boson mass, M A , the ratio of the two neutral Higgs vacuum expectation values, tan β, a common soft SUSY-breaking parameter for the scalar top-and bottom quarks, Mq 3 , a soft SUSY-breaking parameter for the scalar tau and neutrino sector, M˜ 3 , and similarly for the first and second generation of sleptons, M˜ 1,2 , a common trilinear coupling for the third generation, A f , the Higgsino mass parameter, μ, as well as the SU(2) gaugino mass parameter, M 2 . The U(1) gaugino mass parameter M 1 was fixed from the value of M 2 using the GUT relation. The other parameters have been set to fixed values that are generically in agreement with recent SUSY searches at the LHC.
A random parameter scan with O(10 7 ) scan points has been performed. For each scan point a χ 2 function was evaluated, taking into account the combined Higgs boson mass measurement of the LHC experiments and the measured rates in 85 individual Higgs search channels from ATLAS, CMS, and the Tevatron experiments (via the code HiggsSignals), the exclusion bounds from the search for additional Higgs bosons (via the code HiggsBounds, both relying on the evaluations done by FeynHiggs), exclusion bounds from the direct search for SUSY particles (via the code CheckMate), as well as the following low-energy observables BR(B → X s γ ), BR(B s → μ + μ − ), BR(B u → τ ν τ ), (g − 2) μ (as obtained from SuperIso) and M W (via the prediction of Refs. [132,135]).
Taking into account only the Higgs measurements and direct searches, we find that the SM, the MSSM scenario where the light CP-even Higgs boson corresponds to the observed signal ("light Higgs case"), as as well as the MSSM scenario where the heavy CP-even Higgs boson corresponds to the observed signal ("heavy Higgs case") provide similarly good fits to the data with a χ 2 /d.o.f. of 70.2/86, 67.9/79 and 70.0/80, respectively. In a naive evaluation of p values that neglects the correlations between different Higgs observables, this translates into p values of 89%, 81% and 78%, respectively. Including also the low-energy observables we find (via the same evaluation) p values of 69%, 89% and 80%, where the SM suffers in particular from the inclusion of (g − 2) μ . Thus, the "light Higgs case" of the MSSM and even the rather exotic "heavy Higgs case" of the MSSM provide a slightly better descriptions of the data in a global fit than the SM.
Within the MSSM, a SM-like Higgs boson at ∼125 GeV can be realized in three ways and with similarly good fit qualities. For M A M Z the light CP-even Higgs boson is SM-like in the decoupling limit. For M A < ∼ 350 GeV (the "low-M A case") the light CP-even Higgs boson can be SM-like in the limit of alignment without decoupling. This limit also offers additionally the unique possibility that the heavy CP-even Higgs boson can have a mass of ∼ 125 GeV with SM-like couplings. We have analyzed analytically in which parts of the MSSM parameter space the limit of align-ment without decoupling can be realized. In this latter scenario all MSSM Higgs bosons are relatively light, offering good prospects for the searches for additional Higgs bosons at the LHC and future colliders. Our analytic expressions contain the leading two-loop contributions of O(α s h 2 t ) [47], which somewhat modify the leading one-loop contributions of O(h 2 t ) that had been considered previously [33]. For the light Higgs case in the decoupling limit and in the limit of alignment without decoupling, as well as in the heavy Higgs case we have analyzed the predictions for the various Higgs boson production and decay rates. In the light Higgs case the various rates are predicted to be close to the SM Higgs boson rates, where the largest allowed deviation to smaller values is found in the h → τ + τ − channel. The light scalar top is found to have masses down to ∼ 300 GeV, while the preferred region in the fit provides no upper limit on the scalar top quarks. For the low-M A case we find lighter stop masses down to ∼500 GeV and X t /M S ∼ +2. While in the decoupling limit the parameters μ and A t can vary from very small to very large values, the low-M A case (and thus the limit of alignment without decoupling in the light Higgs interpretation) can be realized only for μ/M S ∼ 1.4 to 3 and A t /M S ∼ 2.4 to 3, where the upper values correspond to the upper scan ranges imposed in our study. For larger M A values also relatively large values of tan β are still allowed, leading to a sizable contribution to (g −2) μ even for relatively large chargino/neutralino and slepton masses. Concerning the impact of the low-energy observable, it is interesting to note that a clear preference for positive μ (and also positive A t ) already emerges when combining the Higgs-and B-physics observables in the fit, i.e. already before including (g − 2) μ . We have furthermore found that the preferred region in the fit without taking into account the low-energy observables includes predictions for BR(B s → μ + μ − ) that are close to the experimental central value (i.e., below the SM prediction) both for the decoupling and the alignment without decoupling region of the light Higgs case.
We included the limits from direct LHC searches for SUSY particles from the 8 TeV run via CheckMATE. These searches constrain the parameter space of the pM SS M 8 in an orthogonal way to the Higgs mass and signal rate constraints and therefore do not directly alter the Higgs phenomenology, neither in the light nor the heavy Higgs case. Furthermore, due to this orthogonality we also expect that future stronger SUSY limits from the 13 TeV run would not substantially alter the conclusions found in this paper.
In the heavy Higgs case the preferred rates are also SMlike, however with a possibly larger suppression of H → γ γ and/or H → τ + τ − . We find that in the heavy Higgs case the CP-odd Higgs boson A is restricted to have a mass of 140 GeV < ∼ M A < ∼ 185 GeV, with the charged Higgs boson being the heaviest Higgs boson with M H ± < ∼ 210 GeV, and 6 < ∼ tan β < ∼ 11. We furthermore find μ/M S ∼ 6 to 9 and −1 < ∼ A t /M S < ∼ 0. The light scalar top is predicted to have a mass of 350 GeV < ∼ mt 1 < ∼ 650 GeV. Due to the relatively small tan β values the scalar leptons have to be relatively light with masses below ∼ 450 GeV to bring the prediction into agreement with the observed discrepancy of the experimental measurement of (g − 2) μ with the SM prediction. In particular, we have checked that the charged Higgs corrections to the B-physics observables, B → X s γ , B u → τ ν τ and B s → μ + μ − , are consistent with the experimental results. The preferred region in the fit has mostly light Higgs boson masses above M H /2, in which case H → hh decays are kinematically closed. However, also smaller M h values are possible, where BR(H → hh) does not exceed 20% (for the most favored parameter region we find BR(H → hh) < 2%.) The coupling of the light Higgs boson to W ± and Z bosons is strongly suppressed, much below the existing bounds from LEP Higgs searches. Note that such a light Higgs boson re-opens the possibility of light neutralino dark matter in the sub-GeV to 65 GeV range by acting as s-channel (near-)resonance in dark matter pair-annihilation [174].
As a guidance for the Higgs boson searches in the heavy Higgs interpretation we provide a new set of benchmark scenarios that can be employed to maximize the sensitivity of the experimental analysis to this interpretation. In the (μ, tan β) plane we define the low−M alt− H and the low−M alt+ H scenario with M H ± = 155 GeV < m t and M H ± = 185 GeV > m t , respectively. In the (M H ± , tan β) plane we define the low−M alt v H scenario with μ = 6 TeV. We have shown that in all three scenarios a parameter regime exists, where the mass and rates of the heavy CP-even Higgs boson are in agreement with all available measurements, and which is also not excluded by searches for additional Higgs bosons. In the low−M alt+ H scenario we find a very restricted part of the parameter space in which M h ∼ M H ∼ 125 GeV, with both CP-even Higgs bosons contributing to the Higgs boson rates. The proposed benchmark scenarios will be of interest for upcoming charged Higgs boson searches and will provide motivation for searches for new charged Higgs signatures such as H ± → W ± h, with the light Higgs h decaying into bottom quark or tau lepton pairs, h → bb, τ + τ − .
New data from the ATLAS and CMS Higgs measurements and the search for new Higgs boson states are now rapidly emerging in the current run of the LHC. It is critical to improve the precision of the measurements of the properties of the SM-like Higgs boson, while improving the sensitivity of the searches for new Higgs bosons with masses either above or below the observed Higgs boson mass. In particular, these searches will yield new constraints on the parameter space of the MSSM. The observation of the SMlike Higgs boson already implies that the MSSM Higgs sector lies close to the alignment limit. Indeed, the regions of the MSSM parameter space in which the approximate alignment limit is realized provide a description of the data that is as good (or in some cases slightly better) than the SM. In addition to the possibility that the observed Higgs boson corresponds to the lighter CP-even Higgs boson of the MSSM, the more exotic possibility in which the heavier CP-even Higgs boson is identified as the observed Higgs boson cannot yet be ruled out. Higgs studies at Run 2 of the LHC may prove decisive in determining whether the cracks in the Standard Model facade finally shatter, and whether a supersymmetric interpretation of Higgs phenomena is ultimately viable. der how difficult it is to achieve regions of the MSSM parameter space with very small values of Z 6 . That is, how tuned is approximate Higgs alignment without decoupling in the MSSM?
Consider the case of exact alignment in the two Higgs doublet model (2HDM), which corresponds to Z 6 = 0. If exact alignment is a consequence of a global (discrete or continuous) symmetry, then it is natural to consider 2HDM Higgs sectors that exhibit approximate alignment without decoupling. The possibility of achieving exact alignment by a symmetry of the scalar potential was considered in Refs. [56,175]. However, it is not clear whether these symmetries survive once a realistic Higgs-fermion Yukawa sector is considered [176]. The one known example of a realistic 2HDM with exact alignment is the inert doublet model [177], where Z 6 = 0 is a consequence of an exact Z 2 discrete symmetry, under which the inert doublet field is odd and all other fields are even.
In the MSSM, there is no symmetry associated with the parameter regime corresponding to exact alignment. Indeed, as explained in Sect. 2.2, exact alignment in the MSSM is a result of an accidental cancellation between tree-level and radiatively corrected loop effects that contribute to the effective Z 6 parameter. In an approximate one-loop expression, this cancellation can be explicitly seen in Eq. (20), where the expression inside the braces vanishes for a particular choice of tan β. Of course, given the limited statistics of the present day Higgs data, a region of the MSSM space that exhibits approximate alignment without decoupling can still be consistent with all known experimental constraints. Thus, in any comprehensive scan of the MSSM parameter space, one must necessarily encounter regions of approximate alignment without decoupling. In the absence of a fundamental underlying theory of supersymmetry breaking, the correct measure that governs the MSSM parameter space is unknown. Thus, in order to get a sense of the extent of the tuning associated with the parameter regime of alignment without decoupling, the best one can do is to examine the frequency that such parameter points occur in a comprehensive parameter scan (with uniform priors).
In this appendix, we address this last point from a practical point of view. We examine the values of Z 6 as a function of M A that arise in our pM SS M 8 parameter scans. Note that in the approximate alignment limit, m 2 h Z 1 v 2 = (125 GeV) 2 , i.e., Z 1 0.26, which sets the "natural" size for the other Higgs basis quartic self-couplings (including Z 6 ). In Fig. 28 we show |Z 6 /Z 1 | as a function of M A for the light Higgs (left) and heavy Higgs (right) interpretation. The parameters Z 1 and Z 6 are calculated using the approximate two-loop formulas as described in Sect. 2.2. In the case where h is SM-like, one can see that |Z 6 /Z 1 | ∼ 0.2 in the region of approximate alignment without decoupling. That is, Z 6 is suppressed by less than one order of magnitude relative to Z 1 . 24 In the case where H is SM-like, values between 0 and 0.2 are found for |Z 6 /Z 1 |, with the best-fit value of |Z 6 /Z 1 | 0.1. In both the light and heavy Higgs case, the values of |Z 6 | in the regions of approximate alignment without decoupling are not unnaturally small. Indeed, even in the decoupling regime of large M A , the typical values of |Z 6 | in the preferred MSSM parameter regime are not significantly different in magnitude. We conclude that given the present precision of the Higgs data, approximate alignment without decoupling can be achieved without resorting to an excessive fine tuning of the MSSM parameters. Tables 7 and 8 list the 85 signal strength measurements from ATLAS, CMS and the Tevatron (DØ and CDF), which we