Update of the global electroweak fit and constraints on two-Higgs-doublet models

We present an update of the global fit of the Standard Model electroweak sector to latest experimental results. We include new kinematic top quark and $W$ boson mass measurements from the LHC, a $\sin^2\theta^{\ell}_{\mathrm{eff}}$ result from the Tevatron, and a new evaluation of the hadronic contribution to $\alpha(M_Z^2)$. We present tests of the internal consistency of the electroweak Standard Model and updated numerical predictions of key observables. The electroweak data combined with measurements of the Higgs boson coupling strengths and flavour physics observables are used to constrain parameters of two-Higgs-doublet models.


Introduction
Since the 1990'ies, electroweak precision data from LEP and SLD [1] were used together with accurate Standard Model (SM) calculations to predict parameters of the theory. A first impressive confirmation of the predictive power of global fits in high-energy physics (HEP) was the discovery of the top quark at the Tevatron [2,3] in 1995, with a mass in agreement with the predictions from global fits. Knowledge of the top quark mass (m t ) made it possible to constrain the mass of the Higgs boson (M H ). Increasing experimental and theoretical precision and the inclusion of constraints from direct Higgs boson searches from LEP and Tevatron narrowed the allowed mass range over time [4][5][6][7][8]. The discovery of the Higgs boson at the Large Hadron Collider (LHC) [9, 10] with a mass around 125 GeV impressively confirmed the SM at the quantum level. The historical development of the constraints is illustrated in Figs. 1 and 2, where the predictions of, respectively, m t and M H as derived from various global fits and direct measurements [2, 3,[9][10][11][12][13][14][15][16][17][18][19][20]] are shown versus time.
With the measurement of M H the electroweak sector of the SM is overconstrained and the strength of global fits can be exploited to predict key observables such as the W boson mass and the effective electroweak mixing angle, with a precision exceeding that of the direct measurements [21]. Since the last update of our fit [22] improved experimental results have become available that allow for more accurate tests of the internal consistency of the SM. Among these are the first determination of the W boson mass at the LHC by the ATLAS collaboration [23], new combined results of the top quark mass by the LHC experiments [14,20], a new combination of measurements of the effective leptonic electroweak mixing angles from the Tevatron experiments [24], a Higgs boson mass combination released by the ATLAS and CMS collaborations [17], and an updated value of the hadronic contribution to the running of the electromagnetic coupling strength at the Z boson mass [25]. In the first part of this paper we present an update of the electroweak fit including these new experimental results and up-to-date theoretical predictions.
While the Higgs boson measurements so far agree with a minimal scalar sector as implemented in the SM, the question remains whether a more complex scalar sector may be realised in nature, possibly featuring a variety of Higgs boson states. Two-Higgs-doublet models (2HDM) [26] are a popular SM extension in which an additional SU (2) L ×U (1) Y scalar doublet field with hypercharge Y = 1 is added to the SM leading to the existence of five physical Higgs boson states, h, H, A, H + , and H − , where the neutral h may be identified with the discovered 125 GeV Higgs boson as is assumed in this paper. The scalar H boson has CP-even quantum number, A is a CPodd pseudo-scalar, and H + and H − carry opposite electric charge but have identical mass. No experimental hint for additional scalar states has been observed so far in direct searches [27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42][43]. In this situation global 2HDM fits, exploiting observables sensitive to these additional Higgs boson states via quantum corrections, can be used to constrain the allowed mass ranges and 2HDM mixing parameters. In the second part of this article such constraints are derived from a global fit using a combination of electroweak precision data, flavour physics observables, the anomalous magnetic moment of the muon, and measurements of the Higgs boson coupling strength to SM particles.  Figure 1: Prediction of the top quark mass versus year as obtained by various analysis groups using electroweak precision data (grey [5], light blue [4], green [6]). The bands indicate the 68% confidence level. The direct m t measurements after the top quark discovery are displayed by the data points (orange [2, 3,11,15,16], red [13,14,20], black [12]).  Figure 2: Prediction of the Higgs boson mass versus year as obtained by various analysis groups using electroweak precision data (grey [5], light blue [4], dark blue [6]) and including direct search results (green [6]).
The bands indicate the 68% confidence level. The direct M H measurements after the Higgs boson discovery are displayed by the red data points [9, 10, 17-19].

Update of the global electroweak fit
The updated global electroweak fit presented in this section uses the Gfitter framework. For a detailed discussion of the experimental data, the implementation of the theoretical predictions, and the statistical procedure employed by Gfitter we refer the reader to our previous publications [8,21,22,44]. A detailed list of all the observables, their values and uncertainties used in the fit, is given in the first two columns of Table 1. The description below discusses recent changes in the input quantities and calculations.

Input measurements and theoretical predictions
The electroweak precision data measured at the Z pole and their correlations [1]

Top quark mass
For lack of a recent m t world average, we attempt here for the purpose of the fit a conservative combination of the most precise kinematic m t measurements obtained at the LHC. We combine the m t averages from ATLAS (172.51 ± 0.27 stat ± 0.42 syst GeV) [20] and CMS (172.47 ± 0.13 stat ± 0.47 syst GeV) [14], which are based on 7 and 8 TeV data. These averages include results from the dilepton [51-53], lepton+jets [13, 54] and fully hadronic [55] channels. Assuming the overlapping fraction of the systematic uncertainties to be fully correlated (which corresponds to a correlation coefficient of 72% between the two measurements) we obtain the combined value m t = 172.47 ± 0.46 GeV (p-value of 0.84), which we use as input in the fit.
The latest average from the D0 collaboration m t = 174.95 ± 0.40 stat ± 0.64 syst GeV [16] is barely compatible with the aforementioned average of the LHC measurements. A combination of the D0 average with the LHC average would result in p-values between 5·10 −3 and 3·10 −5 , depending on the assumed correlation between the systematic uncertainties. The result from the CDF collaboration, m t = 173.16 ± 0.57 stat ± 0.74 syst GeV [56], agrees with the LHC average, with p-values between 0.40 and 0.51 depending on the correlation.
As in our previous work [22] we assign an additional theoretical uncertainty of 0.5 GeV to the value of m t from hadron collider measurements due to the ambiguity in the kinematic top quark mass definition [57][58][59][60][61], the colour structure of the fragmentation process [62,63], and the perturbative relation between pole and MS mass currently known to three-loop order [64][65][66].

Theoretical calculations
The theoretical higher-order calculations used in Gfitter have not changed since our last update [22], except for new bosonic two-loop corrections to the Zbb vertex [67].
For the effective weak mixing angle sin 2 θ f eff we use the parametrisations provided in [67][68][69], which include full two-loop electroweak [68,69] and partial three-loop and four-loop QCD corrections [70][71][72][73][74][75][76][77]. For bottom quarks, the calculations from Refs. [67,78] are used. The new bosonic two-loop corrections are numerically small. They shift the prediction of the forward-backward asymmetry for b quarks A 0,b FB by 1.3 · 10 −5 , which is two orders of magnitude smaller than the experimental uncertainty and thus does not alter the fit results. We use the parametrisation of the full two-loop result [79] for predicting the mass of the W boson, where we also include four-loop QCD corrections [75][76][77]. Full fermionic two-loop corrections for the partial widths and branching ratios of the Z boson and the hadronic peak cross section σ 0 had are used [80][81][82]. The dominant contributions from final-state QED and QCD radiation are included in the calculations [83][84][85][86][87][88]. The width of the W boson is known up to one electroweak loop order, where we use the parametrisation given in Ref. [89].
The size and treatment of theoretical uncertainties are unchanged with respect to our last analysis [22].

Results
The fit uses as input observables the quantities and values given in the left rows of Table 1. The fit parameters are M H , M Z , m c , m b , m t , ∆α (5) had (M 2 Z ), α S , as well as ten theoretical uncertainty (nuisance) parameters constrained by Gaussian functions (see Ref. [22] for more details).
The fit results in a minimum χ 2 value of 18.6 for 15 degrees of freedom, corresponding to a p-value of 0.23. The results of the full fit for each observable are given in the fourth column of Table 1, together with the uncertainties estimated from their ∆χ 2 = 1 profiles. The fifth column in Table 1 gives the results obtained without using the experimental measurement corresponding to that row in the fit (indirect determination of the observable). The last column in Table 1 corresponds to the fits of the previous column but ignoring all theoretical uncertainties [22].
The left-hand panel of Fig. 3 displays the pulls each given by the difference of the global fit result of an observable (fourth column of Table 1) and the corresponding input measurement (second column of Table 1) in units of the measurement uncertainty. The right-hand panel of Fig. 3 shows the difference between the global fit result (fourth column of Table 1) as well as the input measurements (first column of Table 1) with the indirect determination (fifth column of Table 1) for each observable in units of the total uncertainty obtained by adding in quadrature the uncertainties of the indirect determination and the input measurement. The analog result using the value of the indirect determination, trivially centered around zero, are shown to illustrate the impact of its uncertainty on the total uncertainty. As in our previous fits, a tension is observed in the leptonic and hadronic asymmetry observables, which is largest in the forward-backward asymmetry of the b quarks, A 0,b FB . The impact of the new Tevatron sin 2 θ eff measurement on the fit result is small due to yet insufficient precision.
which is in agreement with the direct measurement within 1.7 standard deviations. The value is lower by 3 GeV than in our previous result (93 +25 −21 GeV) [22] due to the lower value of m t used here. The reduced uncertainty of +21 −18 GeV compared to +25 −21 GeV previously, is due to the smaller uncertainty in m t . When assuming perfect knowledge of m t , ∆α (5) had (M 2 Z ) and α S (M 2 Z ), the uncertainty is reduced by +4.5 −3.5 , +5 −4 and ±2 GeV, respectively. The predictions of M H using A , A 0,b FB and M W (LEP and Tevatron) concur with earlier findings [8]. The predictions derived   Table 1: Input values and fit results for the observables used in the global electroweak fit. The first and second columns list respectively the observables/parameters used in the fit, and their experimental values or phenomenological estimates (see text for references). The third column indicates whether a parameter is floating in the fit. The fourth column gives the results of the fit including all experimental data. In the fifth column, the fit results are given without using the corresponding experimental or phenomenological estimate in the given row (indirect determination). The last column shows for illustration the result using the same fit setup as in the fifth column, but ignoring all theoretical uncertainties.   4 The results are shown for fits including (blue) and excluding (grey) the direct M H measurement highlighting the strong impact of the M H measurement on the fit constraints. The direct measurement of each observable with its 1σ uncertainty are indicated by the data points at ∆χ 2 = 1. The detailed predictions of the fit are given in Table 1.
The fit indirectly determines the W mass to be and the effective leptonic weak mixing angle as  When evaluating sin 2 θ eff through the parametric formula from Ref. [69], an upward shift of 2 · 10 −5 with respect to the fit result is observed, mostly due to the inclusion of M W in the fit. Using the parametric formula the total uncertainty is larger by 0.6 · 10 −5 , as the global fit exploits the additional constraint from M W . The fit also constrains the nuisance parameter associated with the theoretical uncertainty in the calculation of sin 2 θ eff , resulting in a reduced theoretical uncertainty of 4.0 · 10 −5 compared to the 4.7 · 10 −5 input uncertainty.
The mass of the top quark is indirectly determined to be with a theoretical uncertainty of 0.6 GeV induced by the theoretical uncertainty on the prediction of M W . The largest potential to improve the precision of the indirect determination of m t is through a more precise measurement of M W . Perfect knowledge of M W would result in an uncertainty on m t of 0.9 GeV.
The strong coupling strength at the Z-boson mass scale is determined to be which corresponds to a determination at full next-to-next-to leading order (NNLO) for electroweak and strong contributions, and partial strong next-to-NNLO (NNNLO) corrections. The theory uncertainty of this result is 0.0009, which is shared in equal parts between missing higher orders in the calculations of the radiator functions and the partial widths of the Z boson. The most important constraints on α S (M 2 Z ) come from the measurements of R 0 , Γ Z and σ 0 had , also shown in

Global fits in the two-Higgs-doublet model
Combining information from the electroweak precision data, Higgs boson coupling measurements, flavour observables and the anomalous magnetic moment of the muon we derive in this section constraints on parameters of various 2HDM scenarios.
Besides the four mass parameters for the scalars, M h , M H , M A , and M H ± , the 2HDM introduces the angle α, which describes the mixing of the two neutral Higgs fields h and H, and the angle β that fixes the ratio of the vacuum expectation values of the two Higgs doublets, tanβ = v 2 /v 1 . We only consider 2HDM scenarios with a Z 2 symmetric potential with a dimension-two softly broken term proportional to the scale parameter M 2 12 .  100 < M H ± < 1000 GeV, 0 ≤ β − α ≤ π, 0.001 < tanβ < 50, and −8 · 10 5 < M 2 12 < 8 · 10 5 GeV 2 . No contribution from new physics other than the 2HDM is assumed.
Direct searches for additional Higgs bosons in collider experiments can be interpreted in the context of the 2HDM (see, for example, Ref. [93]). However, due to the large freedom in the choice of the 2HDM parameters, these search results provide only weak absolute exclusion limits on the masses of the scalars. From searches for a charged Higgs boson by the LEP experiments [27] a lower limit of M H ± > 72.5 GeV was reported for the Type-I scenario, while a limit of M H ± 150 GeV can be derived from searches at the LHC for the Type-II scenario [93]. Stronger mass limits mainly on M H ± can be obtained for specific regions of tanβ.

Constraints from Higgs boson coupling measurements
A second Higgs doublet modifies the coupling strengths of the lightest neutral Higgs boson h to SM particles compared to those of the SM Higgs boson. The modifications depend on the 2HDM scenario and parameters in particular the angles α and β. Constraints on h are derived from the joint ATLAS and CMS Higgs boson coupling analysis [94] in which measurements sensitive to five Higgs boson production modes (ggF, VBF, W H, ZH, ttH) and five decay modes (γγ, W W , ZZ, τ τ , bb) were combined. We make use of the relative signal strengths µ ij defined as the ratio of measured over predicted cross section times branching ratio, µ ij = (σ i · B j )/(σ SM i · B SM j ). We include the 20 (out of the 25 possible) µ ij parameters determined by ATLAS and CMS together with their uncertainties and correlations. A validation of our results is discussed in the Appendix on page 23.
The corresponding SM predictions and uncertainties are taken from Ref. [95]. The signal strength measurements are compared with the theory predictions for the 2HDM calculated with the program 2HDMC [96]. 5 In the calculation of the µ ij for the 2HDM also the denominator σ SM i · B SM j is determined using 2HDMC for consistency. Since more precise theory predictions for the SM cross sections and branching ratios exist and are used for the normalisation of the results in [94], theory uncertainties in the SM prediction from [95] are taken into account as additional scaling (nuisance) parameters of the µ ij .
The constraints from the Higgs boson signal strength measurements on the four 2HDM scenarios are shown as 68% and 95% CL allowed regions in the tanβ versus cos(β − α) plane in Fig. 8. 6 The angles α and β are highly constrained in all 2HDM scenarios except for Type-I. The allowed parameter regions are concentrated in two bands corresponding to solutions with β ± α = π/2. For β − α = π/2, the Yukawa structure of the SM is reproduced (alignment limit). The case β + α = π/2 differs from the SM-like Yukawa couplings by a sign flip that is still allowed by the combined coupling strengths measurements. These constraints are differently pronounced in the four 2HDM scenarios as they depend on the Yukawa coupling strengths. In the Type-I scenario (top left panel in Fig. 8) the Yukawa couplings of h to all fermions are proportional to cos α/ sin β. The constraints are stronger in the other three scenarios as the Yukawa coupling for at least one fermion type is proportional to − sin α/ cos β. In the flipped scenario (bottom right panel) only the Yukawa coupling to down-type quarks is given by − sin α/ cos β, which is constrained by the measurements of H → bb. Measurements of H → τ + τ − give stronger bounds in the Type-II (top right panel) and lepton specific (bottom left panel) scenarios where the Yukawa couplings to leptons is given by − sin α/ cos β. In all scenarios, the measurements of Higgs boson decays to W and Z boson pairs disfavour large values of cos(β − α). Similar constraints have been obtained by the ATLAS collaboration [97].

Constraints from flavour observables
Because tree-level FCNC transitions are forbidden by construction in the four 2HDM scenarios considered, flavour violation only arises at loop level by the exchange of a charged Higgs boson with observable strength depending on the parameters M H ± and tanβ. 5 2HDMC computes the couplings of all five Higgs boson states to SM particles for a given set of parameters in a CP conserving 2HDM with general Yukawa structure. From these couplings, production and decay rates of the Higgs boson states can be derived. Most decay widths are calculated at leading QCD order in 2HDMC. Higher order QCD corrections are included for couplings to fermion and gluon pairs. 6 Theoretical bounds from positivity of the Higgs potential, tree-level unitarity, and perturbativity of the quartic Higgs boson couplings as implemented in 2HDMC were found to give no additional constraints in these figures.

Experimental input data and theory calculation
The flavour physics observables taken into account in our analysis are listed in Table 2 and briefly described below.
For the branching fraction of the radiative decay B(B → X s γ) with E γ > 1.6 GeV we use the value of the Heavy Flavour Averaging Group (HFLAV) [98] which combines measurements from the BABAR [105][106][107], Belle [108][109][110], and CLEO [111] experiments. The prediction for B(B → X s γ) has been adopted from Ref. [100] and includes QCD corrections up to NNLO [112]. We make use of the code implementation kindly provided by M. Misiak.
For the branching ratio B(B → τ ν) we use the HFLAV average [98] of measurements from BABAR [120] and Belle [121]. For the prediction of B(B → τ ν) in the 2HDM we use the calculation from Ref. [122], which contains tree-level contributions of a charged Higgs boson where the leading tanβ corrections are resummed to all orders [122]. The theoretical uncertainties in |V ub | and f B d (see below) are included.
The latest measurements of B(B s → µµ) and B(B d → µµ) from LHCb [104] are combined in our fits with the CMS result [103], assuming them uncorrelated. Their theoretical predictions in the 2HDM include NLO corrections given in Refs. [123,124]. The SM contribution to these observables are known up to three-loop level in QCD and include NLO electroweak corrections [125][126][127]. The predictions depend on the CKM matrix elements |V tb | and |V ts | or |V td |, respectively, and on the respective hadronic parameters f Bs and f B d . Uncertainties in these parameters are taken into account in the fit.
The charged Higgs boson of the 2HDM contributes to the leptonic decays of D s mesons. For the observables B(D s → µν) and B(D s → τ ν) we use the HFLAV averages [98] of measurements from BABAR [128], Belle [129], and CLEO [130][131][132]. For the 2HDM predictions we use the analytic expression for the 2HDM tree-level contribution to B(D s → ν) from Ref. [133] that allows us to vary the dependencies on |V cs | and f Ds in the fit.   [138]. For the CKM matrix element |V ub | we use the average of inclusive and exclusive measurements [140], while all other CKM matrix elements are calculated from the Wolfenstein parameters.
The charged Higgs boson also contributes via box diagrams to the mixing of the neutral B d and B s mesons altering the mixing frequencies ∆m d and/or ∆m s . We use again the HFLAV [98] experimental averages for these quantities. Their predictions in the 2HDM are obtained from analytic expressions of the full one-loop calculation of Refs. [119,134] neglecting small terms proportional to m 2 b /M 2 W . The predictions depend on the CKM matrix elements |V td | and |V ts |, the bag parametersB d andB s , and the decay constants f B d and f Bs , respectively, and the correction factor η B .
Finally, the 2HDM contributes at leading order to the ratio B(K → µν)/B(π → µν) for which we use a value adopted from Ref.
[45], based on the measurement of the kaon decay rates [135], and the 2HDM prediction from Ref. [119]. The ratio involves the CKM matrix elements |V us | and |V ud |, the decay constants f K and f π , and an electromagnetic correction δ K/π EM .
As input values for the unitarity CKM matrix we use the latest available results for the allorders Wolfenstein parameters A, λ, ρ, η from Ref. [136][137][138], taking them uncorrelated. A fully consistent analysis would require a combined fit of the Wolfenstein and 2HDM parameters within the 2HDM [139], which is however beyond the scope of this paper. Studies in Ref. [119] and by ourselves have shown that the numerical impact of the 2HDM on the CKM parameters is modest. For the CKM element |V ub |, occurring mainly in the prediction of the leptonic B ± branching fraction, we take the average of inclusive and exclusive measurements [140] instead of the CKM fit prediction to allow for a more conservative uncertainty in view of the tension between the inclusive and exclusive results.
The input parameters used in the fit are summarised in Table 3.

Results
Since most flavour observables are only sensitive to M H ± and tanβ, separate scans of these parameters are performed for each observable. The other 2HDM parameters are ignored in these scans, with the exception of B(B s/d → µµ), where in addition M H , M A , and M 2 12 are allowed to float freely within the bounds defined in the introduction of Section 3 as these two observables depend at NLO level on these parameters. In all fits the CKM matrix elements and the other parameters given in Table 3 are allowed to vary within their uncertainties. Figure 9 shows for the four 2HDM scenarios the one-sided 95% CL excluded regions in the tanβ versus M H ± plane as obtained from fits using the most sensitive individual flavour observables. The CLs are derived assuming a Gaussian behaviour of the test statistic with one degree of freedom. The Type-I (top left) and lepton specific (bottom left) scenarios are only weakly constrained allowing to exclude tanβ < 1. Stronger constraints are obtained for the Type-II (top right) and flipped (bottom right) scenarios in which in particular B(B → X s γ) allows to exclude M H ± < 590 GeV. 7 The measurements of R(D) and R(D * ) differ from their SM predictions [101,102,144]. In the 2HDM only the Type-II scenario features a compatible parameter region (at large tanβ and relatively small M H ± , not shown in the upper right plot of Fig. 9), which is, however, excluded by several other observables. Similar results have been reported in Ref. [119]. Because of this incompatibility R(D) and R(D ( * ) ) are excluded from our analysis in the following.

Constraints from the anomalous magnetic moment of the muon
The measured value of the anomalous magnetic moment of the muon a µ = (g µ − 2)/2 shows a long-standing tension with the SM prediction of ∆a µ = (268 ± 63 ± 43) · 10 −11 [25,145], where the first uncertainty is due the the measurement and the second the prediction (see also the recent reanalysis in Ref. [146]). Loops involving 2HDM bosons can modify the coupling between photons and muons. We have adopted the two-loop 2HDM prediction of ∆a µ from Ref. [147], which depends on all 2HDM parameters. We make use of the code implementation kindly provided by H. Stöckinger-Kim. Figure 10 shows the 68% and 95% CL allowed regions in the tanβ versus M H ± plane for the four 2HDM scenarios using only ∆a µ as input. All other parameters of the 2HDM are left free to vary within their respective bounds. Compatibility is found in a narrow band with tanβ 1 and M H ± below about 600 GeV (depending on the scenario), as well as for a region with larger tanβ that broadens with decreasing M H ± . When combined with the constraints from the other flavour observables (cf. Fig. 9), values of tanβ above about 5∼10 remain allowed.

Constraints from electroweak precision data
The electroweak precision data can be used to constrain the 2HDM via the oblique parameters determined in Eq. (6). We use the predictions from Refs. [148][149][150] similar to our previous analysis [44]. The oblique corrections to electroweak observables in the 2HDM are independent of the Yukawa interactions and their impact is identical in the four 2HDM scenarios considered. Figure 11 shows the 68% and 95% CL allowed parameter regions in the neutral Higgs-boson mass plane M A versus M H for fixed charged Higgs-boson masses of 250, 500, and 750 GeV as obtained from fits using only the oblique parameters as input. All other parameters of the 2HDM (including β − α) are free to vary in these scans. While no information on the absolute mass scale of the 2HDM bosons is obtained from the electroweak data, relative masses are constrained. In our previous analysis [44] we showed that the oblique parameters constrain the values of M H and M A to be close to M H ± for fixed β − α = π/2. Removing this restriction (cf. Fig. 11) relaxes the constraint to having either M A close to M H ± , or M H larger than M H ± .

Combined fit
We combine in this section the 2HDM constraints from the Higgs-boson coupling strength measurements, flavour observables, muon anomalous magnetic moment, and electroweak precision data. Figure 12 shows for the four 2HDM scenarios considered the resulting 68% and 95 % CL allowed regions in the M A versus M H plane for fixed (benchmark) charged Higgs-boson masses of 250, 500, and 750 GeV. All other 2HDM parameters are allowed to vary freely within their bounds. Depending on the 2HDM scenario and M H ± , the minimum χ 2 values found lie between 48 and 59 for N dof = 53 (corresponding to p-values between 25% and 68%).
The combined fit leads in all four 2HDM scenarios to a strong alignment of either the H or the A boson mass with that of the H ± boson, owing to the constraint on β − α from the Higgs coupling strength measurements (cf. Fig. 8) in addition to those from the electroweak precision data. In this sense, the fit resembles the result from our previous analysis [44], but replacing the fixed restriction of β − α = π/2 by the Higgs couplings strengths measurements.
The absolute mass limits on M H ± obtained from the flavour observables in the Type-II and flipped scenarios (cf. Fig. 9) exclude the low-M H ± benchmarks, as indicated by the hatched regions in the two right-hand panels of Fig. 12 (where in addition different statistical assumptions are compared: one-sided versus two-sided test statistic and one versus two degrees of freedom 8 ). For these two scenarios pairs of (H, A) masses below ∼400 GeV are excluded for any set of values of the other 2HDM parameters. For the Type-I and lepton specific scenarios no absolute limits on the Higgs boson masses can be derived.

Conclusion
We have presented results for an updated global fit of the electroweak sector of the Standard Model using latest experimental and theoretical input. We include new precise kinematic top quark and W boson mass measurements from the LHC, a sin 2 θ eff measurement from the Tevatron, and a new evaluation of the hadronic contribution to α(M 2 Z ). The fit confirms the consistency of the Standard Model and slightly improves the precision of the indirect determination of key observables.
Using constraints from Higgs-boson coupling strength measurements, flavour observables, the muon anomalous magnetic moment, and electroweak precision data, we studied allowed and excluded parameter regions of four CP conserving two-Higgs-doublet models. Strong constraints on the extended Higgs boson masses are found for the so-called Type-II and flipped scenarios.   Figure 13: Validation of our implementation of the combined ATLAS and CMS Higgs boson coupling measurements: preferred regions from a two-dimensional scan of the coupling strength modifiers κ V and κ F for individual Higgs boson decay channels and their combination. and κ F = 0.92 ± 0.11 with a correlation coefficient of −0.37. Decent agreement with Ref. [94] is seen.
[9] ATLAS Collaboration, Observation of a new particle in the search for the