# Fits of the Electroweak Standard Model and Beyond using Gfitter

###### Abstract

The global fit of the Standard Model to electroweak precision data, routinely performed by the LEP electroweak working groups and others, has been revisited in view of (i) the development of the new generic fitting package, Gfitter, (ii) the insertion of constraints from direct Higgs searches at LEP and Tevatron, and (iii) a more thorough statistical interpretation of the results. This paper describes the Gfitter project, and presents state-of-the-art results for the global electroweak fit in the Standard Model, and for a model with an extended Higgs sector. Example results are an estimation of the mass of the Higgs boson () and a forth-order result for the strong coupling strength (). Using toy Monte Carlo techniques the -value of the SM has been determined (). As an example of a New Physics model constraints are derived for the Two Higgs Doublet Model of Type-II using observables from the and physics sectors.

## I Introduction

Precision measurements allow us to probe physics at much higher energy scales than the masses of the particles directly involved in experimental reactions by exploiting contributions from quantum loops. A prominent example is the global fit of the Standard Model (SM) to electroweak precision data, routinely performed by the LEP electroweak working group and others (for latest results see renton ()), which demonstrated impressively the predictive power of electroweak unification and quantum loop corrections. Several theoretical libraries within and beyond the SM have been developed in the past, which allowed to constrain the unbound parameters of the SM and models of New Physics. However, most of these programs are relatively old, were implemented in outdated programming languages, and are difficult to maintain in line with the theoretical and experimental progress expected during the forthcoming era of the LHC. These considerations led to the development of the generic fitting package Gfitter gfitterweb (); gfitter (), designed to provide a modular framework for complex fitting tasks in high-energy physics, like model testing and parameter estimation problems. Gfitter is implemented in C++ and relies on ROOT functionality. It consists of a core package providing the tools for data handling, fitting and statistical analyses and allows a consistent treatment of statistical, systematic and theoretical errors, possible correlations, and inter-parameter dependencies. Tools provided for statistical analyses include e.g. parameter scans, contours, Monte Carlo (MC) toy analyses and goodness-of-fit -value evaluation. More details on the framework can be found at gfitterweb (). Gfitter performs the minimisation of a test statistics quantifying the deviation of the experimental data from the predictions in a certain physics model. The theoretical calculations are implemented via plug-in libraries for the Gfitter framework. In this paper we report results which are obtained using the first libraries implemented in the Gfitter package: SM predictions of the electroweak precision observables and predictions of and physics observables in a model with two Higgs doublets (2HDM).

## Ii The Global Electroweak Fit

In the global electroweak fit with Gfitter state-of-the-art calculations are compared with the most recent experimental data to constraint the free parameters of the fit and to test the goodness-of-fit. The SM parameters relevant for the global electroweak analysis are the coupling constants of the electromagnetic, weak and strong interactions, and the masses of the elementary bosons and fermions. Electroweak unification and simplifications arising from fixing parameters with insignificant uncertainties compared to the sensitivity of the fit allow to reduce the number of free fit parameters. The remaining free parameters are the coupling parameters and , the masses , , , and . In addition, four free parameters enter to include the theoretical uncertainties of , and the electroweak form factors and .

For the prediction of the electroweak precision observables as measured by the LEP, SLC and Tevatron experiments the most up-to-date calculations are implemented in the Gfitter SM library using the OMS scheme. Wherever possible the results have been cross-checked against the ZFITTER package zfitter (). The full two-loop and leading beyond-two-loop correction are available for the computation of and awramik1 (); awramik2 (). The partial and total widths of the are known to leading order, while for the second order only the leading corrections are available. Among the new developments included is the NNNLO perturbative calculation of the massless QCD Adler function adler (), contributing to the vector and axial-vector radiator functions in the prediction of the hadronic width (and other observables). It allows to fit the strong coupling constant with unique theoretical accuracy. More details on the theoretical computations in Gfitter can be found in gfitter ().

The following experimental measurements are used: The mass and width of the boson, the hadronic pole cross section , the partial widths ratio , and the forward-backward asymmetries for leptons , have been determined by fits to the line-shape measured precisely at LEP (see zsummary () and references therein). Measurements of the polarisation at LEP zsummary () and the left-right asymmetry at SLC zsummary () have been used to determine the lepton asymmetry parameter . The corresponding and -quark asymmetry parameters , the forward-backward asymmetries , and the widths ratios and , have been measured at LEP and SLC zsummary (). In addition, the forward-backward charge asymmetry measurement in inclusive hadronic events at LEP was used to directly determine zsummary (). For the running quark masses and the world average values are used. For we use the most recent phenomenological result hagiwara (). Results presented in this paper are obtained using the combined LEP and Tevatron results on the mass and the width of the boson newwmass (), , , and the latest combined result on the top mass newtopmass (), , presented at this conference.

The direct searches for the SM Higgs Boson at LEP higgsLEP () and the most recent results from the Tevatron higgstev1 (); higgstev2 (), leading to a 95% confidence level (CL) exclusion for and at respectively, are included using a Gaussian approach that quantifies the difference between the observed test statistics (the log-likelihood ratios) and the expected values for the s+b hypothesis using the published values of the respective confidence level (). A contribution to the estimator of the fit is derived for each Higgs mass. We perform global fits in two versions: the standard (“blue-band”) fit makes use of all the available information except for the direct Higgs searches; the complete fit uses also the constraints from the direct Higgs searches.

Due to the restricted space available only example results of the electroweak fit are reported in the following. More fit results and thorough studies of its statistical properties are given in gfitter () where also the perspectives for the LHC, ILC and GigaZ data are discussed. The standard (complete) fit converges at the global minimum value () for 13 (14) degrees of freedom. The estimation for from the standard fit without the direct Higgs searches is and the and intervals are respectively and . The complete fit represents the most accurate estimation of considering all available data. We find . The resulting curve versus is shown in Fig. 1 (left). The shaded band indicates the influence of theoretical uncertainties. The one, two and three standard deviations from the minimum are indicated by the crossings with the corresponding horizontal lines. The and allowed regions of , including all errors, are and , respectively. The inclusion of the direct Higgs search results from LEP leads to a strong rise of the curve below . The data points from the searches at the Tevatron, available in the range increases the estimator for Higgs masses above beyond that obtained from the standard fit.

The strong coupling at the -mass scale is determined by the complete fit with where the first error is experimental and the second due to the truncation of the perturbative QCD series.

Figure 1 (right) compares the direct measurements of and , shown by the shaded/green bands, with the 68%, 95% and 99% CL obtained for three sets of fits: the largest/blue (narrower/yellow) allowed regions are derived from the standard fit (complete fit) excluding the measured values in the fits. The inclusion of the LEP and Tevatron Higgs searches significantly impacts the constraints obtained. Good agreement between direct measurements and indirect fit results is observed. The third set of fits (narrowest/green) results from the complete fit including the measured values. Hence it uses all available information and leads to the narrowest allowed region.

The -value of the global SM fit, quantifying the probability of wrongly rejecting the SM hypothesis, has been evaluated by means of toy MC experiments. For each MC experiment, the complete fit is performed yielding the distribution shown by the light shaded histogram in Fig. 2 (left). The monotonously decreasing curves give the -value of the SM fit as a function of obtained by integrating the distribution between and infinity. The value of the global SM fit is given by , where the first error is statistical and the second accounts for the shift resulting from theoretical uncertainties.

## Iii Constraints in the 2hdm

As an example for a study beyond the SM we investigate models with an extended Higgs sector of two doublets. In the Type-II 2HDM, we constrain the mass of the charged Higgs and the ratio of the vacuum expectation values of the two Higgs doublets using current measurements of observables from the and physics sectors and their most recent theoretical 2HDM predictions, namely zsummary (); haber (), the branching ratio (BR) of hfag (); misiak (), the BR of leptonic decays of charged pseudoscalar mesons ( chang (); hou (), babar (); hou () and flavianet ()) and the BR of the semileptonic decay dtaunu (); kamenik ().

For each observable, individual constraints have been derived in the () plane. Figure 2 (right) displays the resulting 95% excluded regions derived assuming Gaussian behaviour of the test statistics, and one degree of freedom. The figure shows that is mainly sensitive to excluding small values. BR() is only sensitive to for values below 1. For larger values it provides an almost constant exclusion of a charged Higgs lighter than . For all leptonic observables the 2HDM contribution can be either positive or negative since signed terms enter the prediction of the BRs resulting in a two-fold ambiguity in the space.

In addition, we have performed a global fit combining the information from all observables. For the CL calculation in the 2-dim plane we performed toy MC tests in each scan point which allows to avoid the problem of ambiguities in the effective number of degrees of freedom. The 95% CL excluded region obtained are indicated in Fig. 2 (right) by the area below the single solid line. We can exclude a charged Higgs mass below independently of . This limit increases towards larger , e.g., are excluded for .

## References

- (1) P. Renton, “Global Electroweak Fits and the Higgs Boson Mass”, these proceedings.
- (2) Gfitter Group, Gfitter web page, http://cern.ch/Gfitter.
- (3) H. Flächer, M. Goebel, J. Haller, A. Höcker, K. Mönig, J. Stelzer, “Gfitter - Revisiting the Global Electroweak Fit of the Standard Model and Beyond”, to be submitted to Eur. Phys. J. C, 2008.
- (4) A. B. Arbuzov et al., Comput. Phys. Commun. 174, 728 (2006), [hep-ph/0507146].
- (5) M. Awramik et al., Phys. Rev. D69, 053006 (2004), [hep-ph/0311148], and references therein.
- (6) M. Awramik et al., JHEP 11, 048 (2006), [hep-ph/0608099] and references therein.
- (7) P. A. Baikov, K. G. Chetyrkin and J. H. Kühn, arXiv:0801.1821, SFB-CPP-0804, TTP08-01.
- (8) LEP and SLD Elektroweak and Heavy Flavour Working Groups, Phys.Rept. 427,257 (2006), [hep-ex/0509008].
- (9) K. Hagiwara et al., Phys. Lett. B649, 173 (2007), [hep-ph/0611102].
- (10) CDF Collaboration and others, arXiv:0808.0147, FERMILAB-TM-2415.
- (11) TeVatron Electroweak Working Group and others, arXiv:0808.1089, FERMILAB-TM-2413-E.
- (12) ADLO Collaborations and LEP Higgs Working Group, Phys. Lett. B565, 61 (2003), [hep-ex/0306033].
- (13) The TEVNPH Working Group, arXiv:0804.3423, FERMILAB-PUB-08-069-E.
- (14) The TEVNPH Working Group, arXiv:0808.0534, FERMILAB-PUB-08-270-E.
- (15) H. E. Haber and H. E. Logan, Phys. Rev. D62, 015011 (2000), [hep-ph/9909335].
- (16) Heavy Flavour Averaging Group (HFAG), arXiv:0704.3575.
- (17) M. Misiak et al., Phys. Rev. Lett. 98, 022002 (2007), [hep-ph/0609232].
- (18) P. Chang, “Rare decays and new physics”, these proceedings.
- (19) W. S. Hou, Phys. Rev. D48, 2342 (1993).
- (20) BaBar Collaboration, arXiv:0807.4187.
- (21) FlaviaNet Working Group on Kaon Decays, arXiv:0801.1817, FERMILAB-PUB-08-101-T.
- (22) BaBar Collaboration, Phys. Rev. Lett. 100, 021801 (2008,[arXiv:0709.1698].
- (23) J. F. Kamenik and F. Mescia, Phys. Rev. D78, 014003 (2008), [arXiv:0802.3790].