Determination of the muonic branching ratio of the W boson and its total width via cross-section measurements at the Tevatron and LHC

The total W-boson decay width ΓW\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma _W$$\end{document} is an important observable which allows testing of the standard model. The current world average value is based on direct measurements of final state kinematic properties of W-boson decays, and has a relative uncertainty of 2%. The indirect determination of ΓW\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma _W$$\end{document} via the cross-section measurements of vector-boson production can lead to a similar accuracy. The same methodology leads also to a determination of the leptonic branching ratio. This approach has been successfully pursued by the CDF and D0 experiments at the Tevatron collider, as well as by the CMS collaboration at the LHC. In this paper we present for the first time a combination of the available measurements at hadron colliders, accounting for the correlations of the associated systematic uncertainties. Our combination leads to values of BR(W→μν)=(10.72±0.16)%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {BR}(W\rightarrow \mu \nu )=(10.72 \pm 0.16)\%$$\end{document} and ΓW=2113±31\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma _W = 2113 \pm 31$$\end{document} MeV, respectively, both compatible with the current world averages.


Introduction
Precise measurements of the W -boson properties, such as its mass m W and its decay width W , allows testing of the standard model of particle physics. As a matter of fact, the relation between the W -boson mass, m W , the top-quark mass, m t , and the Higgs-boson mass, m H , via loop corrections, allowed a prediction of the mass of the Higgs boson with an uncertainty smaller than 25 GeV. Models beyond the standard model could alter the relation between m W and m H , since new particles can appear in virtual loops. Similarly, the total decay width of the W boson can be altered by new particles.
Within the standard model, the total decay width of the W boson is predicted to be equal to the sum of the partial widths a e-mail: stefano.camarda@cern.ch over three generations of lepton doublets and two generations of quark doublets. The partial widths are expressed as where M ff = 1, N C = 1 for leptonic decays, M ff corresponds to the CKM matrix elements, and N C = 3 · (1 + α s (m W )/π + · · · ) is the colour factor for the quarkdecay modes [1]. Radiative corrections are represented by δ rad f . They depend, among other parameters, on the topquark mass, m t , and the Higgs-boson mass, m H , and are small in the standard model (SM) since a large part of the corrections is absorbed in the measured values of G F = 1.1663787(8)·10 −5 GeV −2 and m W = 80.385±0.015 GeV. The radiative corrections correspond to δ rad ≈ 0.34% for leptons and of δ rad q ≈ 0.40% for quarks [2]. New particle candidates that couple to the W boson and are lighter than m W , would open a new decay channel and alter W . One very prominent example is supersymmetric models in which the W boson can decay to the lightest super-partner of the charged gauge bosons and the lightest super-partner of the neutral gauge bosons. Hence a precise measurement of W might reveal physics beyond the standard model. In addition, assuming standard model relations, the dependence of the partial and total widths of the W boson on the strongcoupling constant allows to determine the value of α s from the hadronic and leptonic branching ratios of the W boson [3].
The total width of the W boson can be measured directly by kinematic fits to the measured decay lepton spectra, such as the transverse momentum of the charged lepton decay p T or the high-mass tail of the transverse mass m T as was performed at CDF and D0 [4][5][6], or via fits to the invariant mass distributions in the qqqq and qqlν final states as was done at the LEP experiments [7][8][9][10][11]. A combination of these direct results, based on kinematic measurements, leads to W = 2085 ± 42 MeV, which is currently used as world average value [12].
An independent determination of the width of the W boson is based on the measurement of the ratio of cross sections of W -and Z -boson production in hadron collisions, defined as where BR(V → ) = V → / V denotes the leptonic branching ratio of the vector-boson (V = W, Z ) decays. The ratio R can be written as where the total cross-section ratio σ W /σ Z is known theoretically to high accuracy [13]. The ratio σ Z → /σ Z was precisely measured by the LEP experiments and therefore the leptonic branching ratio of the W boson, BR(W → ν) = W → ν / W , can be inferred from the measurement of R. The advantage of extracting BR(W → ν) from the crosssection ratio R lies in the fact that many experimental systematic uncertainties of each vector-boson cross-section measurement, such as the uncertainty on the integrated luminosity, are highly correlated and cancel in the ratio. The leptonic width of the W boson in the SM can be predicted by Eq. (1) and is (W → ν) = 226.5 ± 0.1 MeV. 1 The dominant uncertainty is due to the accuracy of m W . Using this value, the total width of the W boson can be extracted by a measurement of the leptonic branching ratio. This approach for the determination of the W -boson width was already pursued by the CDF [14], D0 [15], and CMS [16,17] experiments, leading to measurements of W which have an accuracy comparable to the current world average.
In this paper we present a procedure for a first combination of the individual measurements of the muonic branching ratio of the W boson and of W , accounting for the correlations of the individual systematic uncertainties. We have chosen to focus on the muon decay channel, as it has smaller experimental uncertainties.
The paper is structured as follows: we introduce the basic methodology in Sect. 2 and discuss the selected measurements for the combination in Sect. 3, where we also derive the corresponding fiducial cross-section ratios. The theoretical predictions of the cross-section ratios are discussed in Sect. 4 and the final extraction and combination of W for the different experiments is presented in Sect. 5. The paper concludes with a brief summary and a discussion of the consistency of the results with the direct measurements and with the global electroweak fit in Sect. 6. 1 Taken from Ref. [2], with updated values of m W and G F , and α s (m W ).

Methodology
The production cross section of W and Z bosons in hadron collisions is described by the Drell-Yan process [18] and can be experimentally defined as where N Cand and N Bkg are the number of vector-boson candidates and the expected background events, respectively, and Ldt is the integrated luminosity of the corresponding data sample. The factor is the efficiency of the signal events passing the signal selection criteria, which is typically estimated with simulated samples of the signal process, and corrected for differences in the detector response between data and MC simulation. The efficiency correction can be decomposed as the product of a fiducial acceptance, A, and a detectorinduced correction factor, C, i.e. = A · C. The fiducial acceptance is the ratio of the number of events that pass the geometrical and kinematic requirements in the analysis at generator level over the total number of generated events in a simulated sample of signal process. The advantage of this decomposition is the separation to a large extent of detector and analysis related uncertainties, which enter the factor C, while all model and theoretical uncertainties, such as QCD scales and parton density function (PDF) uncertainties, enter A. The fiducial production cross section σ fid within the detector acceptance volume defined by A, is barely affected by model uncertainties, and is related to the fully inclusive cross section by σ fid = σ incl · A.
The strategy for the combination of several indirect BR(W → μν) and W measurements from various experiments is therefore based on the measured fiducial crosssection ratio which has only negligible model uncertainties and uncorrelated experimental uncertainties between the different experiments. The fiducial ratio R fid can be related to the inclusive ratio R, by where the A W and A Z are the acceptance correction factors for the Z -and W -boson analyses, respectively. Some published results only present a value for the inclusive crosssection ratio R, but do not publish a value for R fid . In these cases, we have used the fiducial volume definition, the PDF set, and the MC generator of the corresponding analysis that were used to extract the acceptance ratio A W /A Z , in order to reconstruct the value of R fid . The uncertainty on the extrapolated values of R fid is estimated by subtracting the published model and PDF uncertainties from the total uncertainty on R.
Once the fiducial ratios are determined for each measurement, we can coherently predict the acceptance correction ratios A W /A Z and the inclusive fiducial cross-section ratios σ W /σ Z , and extract the corresponding branching ratio BR(W → μν) and decay width W from the measurements of R fid . Each model variation, e.g. one particular eigenvector variation of a given PDF set, leads to new predictions of A W /A Z and σ W /σ Z , thus also to new determined values of BR(W → μν) and W for each experiment. The measurements are combined treating the experimental uncertainties as uncorrelated, and the PDF and model uncertainties with a correlation model based on a common baseline for the theoretical predictions.

Measurements used for the combination
One of the first precise measurements of the σ W /σ Z crosssection ratio was published by the D0 collaboration in proton anti-proton collisions at a centre-of-mass energy of √ s = 1.8 TeV [15]. However, this measurement was performed only in the electron decay channel and hence is not used for this combination. The most precise measurement at the Tevatron collider was performed by the CDF collaboration at √ s = 1.96 TeV [14], using the electron and muon decay channels. Only the inclusive cross-section ratio was published (Table 1), but the clear definition of the fiducial volume, as reported in the paper, allows the extrapolation of the value of R fid . The extrapolation factor A W /A Z is estimated using the Pythia 6.2 [19] generator with the CTEQ5L PDF set [20].
Several measurements of R were performed at the LHC by the CMS and ATLAS collaborations in proton-proton collisions at √ s = 7 TeV and √ s = 8 TeV [16,17,21], which are all used for the combination. In contrast to the Tevatron experiments, fiducial ratios together with a fiducial volume definition have also been published by the ATLAS and CMS experiments, as summarised in Table 1. Hence no additional extrapolation to R fid is performed for these measurements.

Theoretical predictions and systematic uncertainties
The total W -and Z -boson production cross sections and their ratio, corresponding to the fiducial volume definitions of Table 1, are calculated at next-to-next-to-leading order in the perturbative expansion of the strong-coupling constant with FEWZ [22] using the MMHT2014 PDF set [23]. The calculations are based on the G μ electroweak parameter scheme and the strong-coupling constant at the Z -boson mass is set to α s (m Z ) = 0.118, as used in the MMHT2014 PDF determination. 2 The uncertainties of the PDF set are estimated by a reevaluation of the predicted cross-section ratio for each error eigenvector within the MMHT2014 PDF set, as well as the comparison to the central prediction using a second PDF set, which is chosen to be the CT10 [24] in this study. The uncertainties, at 68% CL, include contributions from the strong-coupling constant α s as well as variations of the renormalisation scale, μ R , and factorisation scale, μ F .
The correct description of the vector-boson transverse momentum, p T (V = W, Z ), is essential for the estimation of A W and A Z . Since fixed order perturbative QCD predictions do not provide a sufficiently good description of the low p T (V = W, Z ) spectrum, we use the Powheg MC generator interfaced to Pythia8, henceforth referred to as Powheg+Pythia8, to estimate the central values for A W and A Z , using the MMHT2014 PDF set.
The uncertainties due to missing higher order QCD corrections are estimated by varying the renormalisation and factorisation scales, μ R and μ F , by a factor of two up and down, as well as by reevaluating the acceptance factors with hdamp set to m V (V = W, Z ), instead of the default value hdamp = ∞ [25], in the Powheg generator. The correlation of the μ R and μ F variations on the W -and Z -boson cross sections and acceptances can be treated according to various prescriptions. In the most conservative approach they are considered as fully uncorrelated, leading to an uncertainty of 0.5% on the predicted inclusive cross-section ratio. The uncertainty reduces by more than a factor of two when assuming a fully correlated behaviour. In the following we adopt an intermediate approach, and assume a correlation of 50%. In addition to these uncertainties, the acceptance factor ratio A W /A Z is also affected by other effects, which change the kinematic distribution of the final states, but has little effect on the inclusive cross sections. In particular, the uncertainties due to soft non-perturbative effects and initial-state radiation (ISR), which vary the transverse momentum spectrum of the vector boson, p T (V ), have to be estimated. To perform a conservative estimation, we reweight the predicted p T (V ) from Powheg+Pythia8 to corresponding predictions of the Resbos generator [26][27][28]. Resbos is based on a resummed calculation, which is performed at next-to-next-to-leading logarithmic order and matched to approximate NNLO perturbative QCD calculations at large boson momenta. The difference between the nominal Powheg+Pythia8 predictions of A W and A Z and the Resbos reweighted samples, is considered as an ISR and resummation uncertainty. The corresponding uncertainties on A W /A Z vary between 0.1% (ATLAS) and 0.4% (CMS). This difference can be explained by the larger effect on the W -boson selection for CMS, as it Table 1 The collider beams, the corresponding centre-of-mass energy, the fiducial volume definitions, the published inclusive cross-section ratio R, as well as the fiducial cross-section ratio R fid are given for each analysis used for the combination.  requires only a minimum threshold on the p T of the decay muons. Furthermore, NLO electroweak corrections can be comparable in size to NNLO QCD corrections. We distinguish between the corrections due to QED final-state radiation (FSR) and loop-induced electroweak corrections (EWK). The QED FSR related uncertainties are estimated by comparing Sherpa [29] and Pythia8 [30], where the acceptances are derived for both generators using dressed and bare leptons. The resulting differences in the predicted acceptance ratios are taken as the QED FSR model uncertainty and amount to 0.1%. The uncertainties due to loop-induced electroweak corrections are taken from literature [31] and are accounted for by 0.1% variations on A W /A Z .
We obtain 57 predictions for the cross-section ratios σ W /σ Z and the acceptance ratios A W /A Z , accounting for 50 MMHT PDF eigenvector variations, the central prediction of the CT10 PDF set, μ R and μ F scale variations, and variations of α s . In addition, we have further uncertainties on A W /A Z due to ISR/resummation effects, QED FSR and NLO EWK effects. A summary of the cross-section ratios σ W /σ Z and acceptance ratios A W /A Z for each measurement, including the relevant model uncertainties, is given in Table 2. The PDF uncertainties are evaluated with the Hessian method [32].
The uncertainties due to ISR and resummation, QED FSR, and electroweak corrections, as well as the variations of μ F and μ R , are symmetrised by taking the average of the positive and negative variations.

Extraction of the W -boson width and combination
The total inclusive cross-section ratio for each experiment is estimated by combining the central values of A W /A Z reported in Table 2 and the fiducial cross-section ratios R fid from Table 1. It should be noted that these derived values for R will differ from the original published values, as our baseline prediction for the estimation of A W /A Z differs from the approach followed by each experiment. Clearly, the advantage of having a common model for the theoretical predictions lies in the traceability of correlated systematic uncertainties. The published values of R are compared to the values obtained using the newly derived acceptance ratios as a first sanity check of our extrapolation. The derived values agree with the published values of the experiments within their associated model uncertainties.
In a second step, the expected leptonic branching ratios can be rederived for each experiment individually, using the predicted cross-section ratios, the measured fiducial ratios of the experiments reported in Table 2, and the relation where a value of the leptonic Z -boson branching ratio of Z →μμ / Z = 0.033658 ± 0.000023 [12] is used. The results are presented in Table 3. Assuming the validity of the SM, the partial leptonic W -boson width is predicted by Eq. (1), leading to W →μν = 226.5 MeV, where the (1+δ rad ) corrections are taken from Ref. [2]. Finally, the total W -boson width can be derived from the measured leptonic branching ratios. The resulting values for W = BR· W →μν for each experimental measurement are illustrated in Fig. 2 and reported in Table 4, together with the associated statistical, experimental systematic, and combined model uncertainties.
For the combination, we use the measured inclusive crosssection ratio R and the corresponding inclusive cross-section prediction for each model systematic variation, thus leading to new values of BR(W → μν) and W for each measurement, respectively. In a second step, we combine the individual measurements following the BLUE method [33], again, separately for all model variations. For the combination of the four experimental values of BR(W → μν) and W , we treat the statistical and experimental systematic uncertainties as fully uncorrelated. In a last step, we calculate the difference between the combined values of BR(W → ν) and W for each model variation and their central combined values, and evaluate the theoretical and model systematic uncertainties from these differences. The PDF uncertainties are evaluated with the MMHT2014 PDF set using the Hessian method [32]. The other model systematic uncertainties are added in quadrature. The

Summary and interpretation
In this paper we have presented a combination of measurements of the muonic branching ratio of the W boson and its total decay width, extracted from the cross-section ratios of W -and Z -boson production from the ATLAS, CMS, and CDF experiments at various centre-of-mass energies. Special emphasis was drawn to the correct treatment of the correlations between systematic uncertainties, in particular uncertainties due to the limited knowledge of the parton distribution functions and variations of the renormalisation    [6,11], which are based solely on direct measurements.
The total width of the W boson is potentially sensitive to new physics scenarios in the context of the global electroweak fit [38]. The indirect determination via the electroweak fit yields a value of EW W = 2091 ± 1 MeV [39], which is in good agreement with our combined value. The relation expressed in Eq. (1) shows that W depends, among other SM parameters, on m W , α s , and m H . However, the small uncertainties on the determination of EW W indicates that the sensitivity of W to these parameters of the SM is rather weak.
It should be noted that loop corrections arising from contributions of new physics to the W -boson width would alter the term δ rad in Eq. (1) independently from the decay channel. As a consequence, the branching ratio is insensitive to effects that could appear in the corresponding loop correction terms. Only new physics effects that directly alter the leptonic branching ratio can be tested with our combined value BR(W → μν) R . We find a perfect agreement with the SM prediction BR(W → μν) SM = (10.83±0.01)% [1,2]. 3