Introduction
The global fit to electroweak precision data routinely performed by the LEP electroweak working group and others, demonstrates impressively the predictive power of electroweak unification and quantum loop corrections. We have performed the fit using the most recent experimental measurements and stateoftheart SM predictions.Note: We interpret the new boson discovered at the LHC as the SM Higgs boson and use the measurements from ATLAS and CMS for the Higgs boson mass M_{H}. More details can be found in our latest publication The Electroweak Fit of the Standard Model after the Discovery of a New Boson at the LHC.
New since the last publication:
 Changed R_{b} calculation following the publication JHEP 1208 (2012) 050 and Erratumibid. 1305 (2013) 074, see also [arXiv:1205.0299].
 New top mass measurement by the Tevatron experiments [arXiv:1305.3929]
Theoretical predictions
For the following results the latest theoretical predictions of electroweak observables are used. In particular:
 The mass of the W boson (M_{W}) is calculated with the full twoloop corrections and known beyondtwoloop corrections from
(M. Awramik et al., Phys. Rev. D69, 053006 (2004), hepph/0311148).  The effective weak mixing angle (sin^{2}θ_{eff}^{l}) is calculated with the full twoloop corrections and known beyondtwoloop corrections from
(M. Awramik et al., JHEP 0611, 048 (2006), hepph/0608099)
(M. Awramik et al., Nucl.Phys.B813:174187 (2009), arXiv:0811.1364).  The partial and total widths of the Z and of the total width of the W boson make use of the parametrizations of (Cho et. al, arXiv:1104.1769, see also older papers 1, 2, 3).
 The determination of the strong coupling makes use of the complete fourthorder (3NLO) calculation of the hadronic Z width (P. A. Baikov et al., arXiv:1201.5804) .
 Electroweak twoloop corrections to R_{b} (Freitas and Huang, arXiv:1205.0299, v3).
Fit results of the current global fit
In the following tables and figures the experimental input used in the fit and the fit results are given. All fits discussed here minimise the test statistics χ^{2} which accounts for the deviations between the observables given in the table below and their SM predictions. The fit converges at the global minimum value χ^{2}_{min} =18.1 for 14 degrees of freedom, giving the naive pvalue Prob(χ^{2}_{min},14)=0.20. Excluding the experimental information on M_{H} the fit converges at the global minimum value χ^{2}_{min} =16.7 for 13 degrees of freedom, giving a naive pvalue of Prob(χ^{2}_{min},13)=0.21. See the Section Probing the SM for a more accurate toyMCbased determination of the pvalue.With the measurement of M_{H} it is for the first time possible to fully predict all Standard model parameters / observables with only a minimal set of input parameters. The minimal set of parameters needed are M_{H} together with all fermion masses, α_{s}(M_{Z}^{2}) and three parameters defining the electroweak sector and its radiative corrections, here chosen to be M_{Z}, G_{F} and Δα_{had}(M_{Z}^{2}). The prediction based on these sets of parameters is shown in some figures below where it is called "SM fit with minimal input".
Table: Input values and fit results for the observables and parameters of 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. The subscript ``theo'' labels theoretical error ranges. The third column indicates whether a parameter is floating in the fit. The fourth column quotes the results of the complete fit including all experimental data. The fifth column gives the fit results for each parameter without using the M_{H} measurement in the fit. In the last column the fit results are given without using the corresponding experimental or phenomenological estimate in the given row. 

Comparing fit results with direct measurements: pull values for the SM fit, i.e. deviations between experimental measurements and theoretical calculations in units of the experimental uncertainty. 

Comparing fit results with direct measurements: pull values for the SM fit with and without inclusion of M_{H} in the fit. The pull values are defined as deviations between experimental measurements and theoretical calculations in units of the experimental uncertainty. 

Comparing fit results (orange bars) with indirect determinations (blue bars) and direct measurements (data points): pull values for the SM fit defined as deviations to the indirect determinations. The total error is taken to be the error of the direct measurement added in quadrature with the error from the indirect determination. This plot is a graphical representation of the numbers presented in the above table. 

Determination of M_{H} excluding all the sensitive observables from the fit,
except for the one given.


Δχ^{2} as a function of Higgs boson mass M_{H}, shown as blue band. Also shown is the result of the fit without the M_{H} measurements (grey band). The solid and dashed lines give the results when including and ignoring theoretical errors, respectively. 

Δχ^{2} as function of the top mass m_{t}, shown as blue band. Also shown is the result of the fit without the M_{H} measurements (grey band). In both cases the direct determinations of m_{t} were excluded from the fit. Direct determinations of m_{t} are indicated by dots with 1σ error bars. 

Δχ^{2} versus M_{W}, shown as blue band. Also shown is the result of the fit without the M_{H} measurements (grey band). In both cases the direct measurements of M_{W} were excluded from the fit. The Standard Model fit with minimal input (see above) is shown as a black line. The experimental world average of M_{W} is indicated by a dot with 1σ error bars. 

Δχ^{2} versus the effective weak mixing, shown as blue band. Also shown is the result of the fit without the M_{H} measurement (grey band). In both cases all precision observables sensitive to sin^{2}θ_{eff} are excluded from the fit. The Standard Model fit with minimal input (see above) is shown as a black line. The average of the LEP and SLD measurements of sin^{2}θ_{eff} is indicated by a dot with 1σ error bars. 

Δχ^{2} versus the running strong coupling at M_{Z}, shown as blue band. The green band shows the Standard Model fit with minimal input, where the value of α_{S}(M_{Z}^{2}) is left free in the fit and instead measurements of σ^{0}_{had} and R^{0}_{l} are used as input parameters. The value of α_{S}(M_{Z}^{2}) determined from τ decays at N^{3}LO is indicated by a dot with 1σ error bars. 

Contours of 68% and 95% confidence level obtained from scans of fits with fixed variable pairs M_{W} vs. m_{t}. The narrower blue and larger grey allowed regions are the results of the fit including and excluding the M_{H} measurements, respectively. The horizontal bands indicate the 1σ regions of the M_{W} and m_{t} measurements (world averages). 

Probing the Standard Model
We evaluate the pvalue of the global SM fit using a toy MC simulation with 20000 experiments. These are generated using as true values for the SM parameters the outcomes of the global fit. For each toy simulation, the central values of all of the observables used in the fit are generated according to Gaussian distributions around their expected SM values (given the parameter settings) with standard deviations equal to the full experimental errors taking into account all correlations. Fair agreement is observed between the empirical toy MC distribution and the χ^{2} function expected for Gaussian observables.
Result of the MC toy analysis of the SM fit. Shown are the χ^{2}_{min} distribution of a toy MC simulation (open histogram), the corresponding distribution for a fit ignoring theoretical errors (shaded/green histogram), an ideal χ^{2} distribution assuming a Gaussian case with n_{dof}=14 (black line) and the pvalue as a function of the χ^{2}_{min} of the global fit. 
