Measurement of the Mass and Width of the W Boson in e+e- Collisions at 189 GeV

The mass and width of the W boson are determined in e+e- collisions at LEP using 183 pb^-1 of data recorded at a centre-of-mass energy roots=189 GeV with the OPAL detector. The invariant mass distributions from 970 WW->qqqq and 1118 WW->qqln candidate events are used to measure the mass of the W boson, Mw = 80.451 +- 0.076(stat.) +- 0.049(syst.) GeV. A direct measurement of the width of the W boson gives Gw=2.09 +- 0.18(stat.) +- 0.09(syst.) GeV. The results are combined with previous OPAL results from 78 pb^-1 of data recorded with roots from 161 to 183 GeV, to obtain: Mw = 80.432 +- 0.066(stat.) +- 0.045(syst.) GeV, Gw = 2.04 +- 0.16(stat.) +- 0.09(syst.) GeV. The consistency of the direct measurement of Mw with that inferred from other measurements of electroweak parameters provides an important test of the Standard Model of electroweak interactions.


Introduction
The LEP e + e − collider at CERN provides an ideal environment for the study of the properties of the gauge bosons of the Standard Model (SM) of electroweak interactions [1].In the first stage of its operation, LEP produced e + e − collisions at centre-of-mass energies, √ s, within a few GeV of the Z 0 resonance, allowing precise measurements of the properties of the Z 0 boson including its mass and fermionic couplings [2,3].In the context of the SM these measurements place constraints on the mass of the Higgs boson and provide an indirect determination of the mass of the W boson, M W . Since 1996, the LEP collider has operated above the threshold for W + W − production (LEP2), allowing measurements of the trilinear gauge boson couplings [4] and a direct measurement of M W .When combined with the direct measurements of the top quark mass at the Tevatron [5], measurements of M W enable further constraints to be set on the mass of the Higgs boson via electroweak radiative corrections [6,7].Comparison between the direct measurements of the mass of the W boson and the value determined indirectly from data recorded at √ s ≈ M Z provides an important test of the self-consistency of the Standard Model.The direct measurement of Γ W further tests the consistency of the Standard Model.
The published OPAL measurements of M W are based on approximately 78 pb −1 of data.This paper describes a measurement of the mass of the W boson using a further 183 pb −1 of data recorded by OPAL during 1998 at √ s∼ 189 GeV.This result is combined with previous OPAL measurements to give a direct measurement of the mass of the W boson with a total uncertainty of 79 MeV.It is expected that the ultimate LEP precision on M W will be approximately 30 MeV when all data are included and the results of the four LEP experiments are combined [16].

Data and Monte Carlo Samples
A detailed description of the OPAL detector can be found in [17].The data sample used for this analysis corresponds to an accepted integrated luminosity, evaluated using small angle Bhabha scattering events observed in the forward calorimeters [18], of 183.1 ± 0.2(stat.)± 0.4(syst.)pb −1 .The luminosity weighted mean centre-of-mass energy for the data sample is √ s = 188.635± 0.040 GeV [19].

Monte Carlo samples
A number of Monte Carlo generators are used to simulate the physics processes relevant to the studies presented in this paper.All samples include a full simulation of the OPAL detector [20].The main physics processes at LEP2 can be broken down into three main categories: four-fermion (4f ) production, including e + e − → W + W − → 4f and e + e − → Z 0 Z 0 → 4f , but excluding contributions from multi-peripheral diagrams; two-fermion production; and multi-peripheral two-photon mediated processes.For the measurement of the W boson mass and width, only the four-fermion processes and the two-fermion background process e + e − → Z 0 /γ → qq play an important role.
The Koralw program [21], which uses matrix elements calculated with grc4f [22], is used to simulate the production of most four-fermion final states, e + e − → 4f .The main Monte Carlo samples are generated at √ s = 188.634GeV with M W = 80.33 GeV.The samples were generated using the running width scheme for the Breit-Wigner distribution.The four-fermion samples are divided into final states which have contributions from processes involving the W-boson propagator and those which do not.In the invariant mass region close to M W , the four-fermion cross-section is dominated by doubly resonant W-pair production diagrams (CC03) 1 .Additional Koralw four-fermion Monte Carlo samples are produced with different centre-of-mass energies and with different values of M W .The four-fermion background from the process e + e − → e + e − qq is simulated using the grc4f [22] generator.The most important two-fermion background process, e + e − → Z 0 /γ → qq, is simulated using Pythia [23], with Herwig [24] used to assess possible systematic uncertainties.

W + W − Event Selection
The event selections are described in [25] and references therein.The selections are sensitive to the leptonic W + W − → ℓν ℓ ℓν ℓ , semi-leptonic W + W − → qqℓν ℓ , and hadronic W + W − → qqqq final states.Due to the presence of two unobserved prompt neutrinos in the W + W − → ℓν ℓ ℓν ℓ final state there is little sensitivity to M W and the leptonic final state is therefore not used here.
Semi-leptonic W + W − → qqℓν ℓ decays comprise 44% of the total W + W − cross-section.The event selection employs three multivariate relative likelihood discriminants, one for each of the W + W − → qqeν e , W + W − → qqµν µ , and W + W − → qqτ ν τ final states.The W + W − → qqeν e and W + W − → qqµν µ channels are characterised by two well-separated hadronic jets, a high-momentum lepton and missing momentum due to the prompt neutrino from the leptonic W decay.The signature for the W + W − → qqτ ν τ channel is similar, with the exception that the τ lepton is identified as an isolated, low-multiplicity jet typically consisting of one or three tracks.W + W − → qqℓν ℓ events are selected with an efficiency of 87% and a purity of 91%.The dominant backgrounds are Z 0 /γ → qq and four-fermion processes such as e + e − → Weν e and e + e − → (Z 0 /γ) * (Z 0 /γ) * → qqℓ + ℓ − .Hadronic W + W − → qqqq decays comprise 46% of the total W + W − cross-section and are identified by requiring four energetic hadronic jets and little or no missing energy.A preselection removes approximately 98% of the dominant background process, Z 0 /γ → qq.A multivariate relative likelihood discriminant is then employed to select the W + W − → qqqq candidates with an efficiency of 87% and a purity of 77%.
After the selections are applied, 1546 W + W − → qqqq and 1246 W + W − → qqℓν ℓ candidate events remain, consistent with Standard Model expectations.Not all events are used in the measurements described here.As discussed in Section 4.1, additional cuts are applied to remove poorly reconstructed events and further reduce backgrounds.

Measurement of the Mass and Width of the W Boson
The measurement of the mass and width of the W boson proceeds in two stages.Firstly the invariant masses of the W decay products are reconstructed on an event-by-event basis.Kinematic fits are applied to each selected W + W − → qqℓν ℓ and W + W − → qqqq event to improve the mass resolution.The reconstructed invariant mass spectra are then used to determine M W and Γ W .
Fits to the invariant mass spectra to obtain M W are performed using three different techniques.The central results of this paper are obtained using a Monte Carlo reweighting technique [16] to fit the observed mass spectra to obtain M W and Γ W .The W mass is also determined using two alternative methods, which are used as cross-checks.In the first, an analytic fit to the measured mass spectrum uses an unbinned likelihood fit method to determine M W .To describe the signal shape, the fit uses a parametrisation based on a Breit-Wigner function [9].The second method uses a convolution technique [12,16].The three different fitting techniques have similar expected statistical sensitivities and similar estimated systematic uncertainties.

Invariant mass reconstruction
The three methods for extracting M W use nearly identical procedures to reconstruct the invariant mass of the W candidates.The description here applies to the reweighting method.Small variations relevant for the alternative analyses are discussed in Section 4.3.
In previous OPAL measurements of M W the tracks and clusters in selected W + W − → qqqq events were grouped into four jets using the Durham algorithm [26].For the results presented here, W + W − → qqqq events are grouped into either four or five jet topologies depending on the value of ln(y 45 ), where y 45 is the value of the Durham jet resolution parameter at which the transition from 5 to 4 jets occurs.Events with ln(y 45 ) > −5.6 are treated as five jets.This separation allows for the possibility of hard gluon radiation from one of the quarks and is found to improve the expected statistical sensitivity to M W by approximately 5% compared with treating all events as four jets.The improvement in statistical sensitivity comes from events where, under the four jet hypothesis, a gluon jet from one W boson is combined with a quark jet from a different W boson.In the semi-leptonic decay channels, W + W − → qqℓν ℓ , the lepton candidate is removed and the hadronic part of the event is reconstructed as two jets.The division of W + W − → qqℓν ℓ events into two and three jet would not improve the mass resolution.After the association of tracks and clusters into jets, corrections derived from the Monte Carlo simulation are applied to the measured jet momenta to account for double counting from particles which deposit energy in more than one sub-detector [27].
The invariant masses of the two W bosons can be determined directly from the reconstructed momenta of the observed decay products.However, the mass resolution is limited by the relatively large uncertainty on the measured energies of the jets, σ E /E ≈ 20%, rather than by uncertainties on the measured jet directions.For this reason the use of a kinematic fit which imposes the four constraints of energy and momentum conservation [28] (4-C fit) significantly improves the invariant mass resolution.The improvements mainly arise from the energy in event being constrained to the well measured centre-of-mass energy.The output of the 4-C fit consists of two reconstructed masses per event, one for each W boson in the final state.The resolution of the kinematic fit can be improved by imposing a further constraint that the masses of the two reconstructed W boson candidates are equal (5-C fit), yielding a single reconstructed mass per event.For W + W − → qqℓν ℓ events the three unmeasured variables corresponding to the neutrino momentum means that the effective number of constraints in the semi-leptonic is two, giving a 2-C fit.For the results presented in this paper, the mass reconstruction is performed using the 5-C fit for W + W − → qqqq events and the 2-C fit for W + W − → qqℓν ℓ events.
A common kinematic fitting algorithm is used for W + W − → qqeν e , W + W − → qqµν µ and W + W − → qqqq events.The fitted mass is obtained using an iterative χ 2 -minimisation procedure where the constraints in the kinematic fit are implemented using Lagrange multipliers.The presence of initial state radiation (ISR) is neglected.The corrected jet momenta, their associated errors, and the measured jet masses are input to the kinematic fit.The errors on the measured jet momenta are parameterised by expressions derived from Monte Carlo studies, which are functions of the visible energy and polar angle of the jet.The jet masses are fixed to their measured values.This is found to improve the mass resolution from the fit compared to treating jets as massless and also to reduce the bias in the fitted mass distribution.
The mass reconstruction for W + W − → qqτ ν τ events is different.In the W + W − → qqτ ν τ system, the absence of a measurement of the tau lepton energy means that the W mass information is determined entirely by the hadronic system.The absence of a measurement of the tau lepton energy reduces the effective number of constraints in the kinematic to one.Previously [10] the reconstructed mass in W + W − → qqτ ν τ events was obtained from the invariant mass of the jet-jet system, scaled by the ratio of the beam energy to the sum of the jet energies.For the results presented here an analytic formula which reproduces the results of the 1-C kinematic fit is used.
For W + W − → qqeν e and W + W − → qqµν µ events the lepton direction is taken to coincide with the direction of the track associated with the electron or muon candidate.The energy is estimated from the associated electromagnetic calorimeter cluster for electrons and from the momentum of the track for muons.Unassociated electromagnetic clusters close to the lepton track, consistent with being from final state radiation (FSR), are included in the energy calculation.For each event a single mass is determined from a 2-C kinematic fit.In addition to the fitted mass, m rec , the error on the fitted mass, σ rec , and the χ 2 fit probability, P FIT , are calculated.Events with m rec > 65 GeV and P FIT > 0.001 are retained.About half of the W + W − → qqτ ν τ events selected as either W + W − → qqeν e or W + W − → qqµν µ fail these.The mass information is recovered by treating events with m rec < 65 GeV or P FIT < 0.001 as W + W − → qqτ ν τ events.In addition, about 4% of selected W + W − → qqeν e and W + W − → qqµν µ events have the identified lepton beyond the effective tracking acceptance of the OPAL detector, within 20 • of the beam axis, in which case the lepton energy is either poorly measured or not measured at all.These events are included in the W + W − → qqτ ν τ channel.The numbers of events used in each of the W + W − → qqℓν ℓ channels are given in Table 1.
The situation for W + W − → qqqq events is complicated by the fact that there are three possible assignments of four jets to the two W bosons.For events reconstructed as five jets there are ten possible assignments of the jets to the W bosons. Incorrect combinations contain little or no information on the mass of the W boson.The incorrect jet-pairings result in a combinatorial background.For each W + W − → qqqq event, three (or ten) kinematic fits are performed, corresponding to the possible jetpairings.To eliminate poorly reconstructed events and reduce backgrounds, only combinations which give a successful 5-C kinematic fit with a resulting χ 2 fit probability P FIT > 0.01 and m rec > 65 GeV are considered.A multivariate relative likelihood discriminant similar to that described in [10] is employed to pick out a single combination for each event and reduce combinatorial background.Different variables are used for the jet-pairing likelihoods for the four and five jet cases.For four jet events, the combination used in the fit is selected on the basis of two variables: the difference between the two fitted masses from the 4-C fit and the sum of the di-jet opening angles.In the five jet sample four variables are used: the 5-C fit mass, the difference between the two fitted masses from the 4-C fit, the minimum opening angle between the jets in the system assigned as W ± → qqg and the cosine of the polar angle of the reconstructed W ± → qqg system.For events reconstructed as four (five) jets, the combination corresponding to the largest jet-pairing likelihood is retained provided it has a likelihood output exceeding 0.40(0.42).Monte Carlo studies indicate that in 89% (70%) of the surviving signal events, the selected combination corresponds to the correct jet-pairing.The number of surviving events in the W + W − → qqqq channel is given in Table 1. Figure 1 shows the reconstructed mass distributions in both four and five jet channels before and after the jet-pairing likelihood cuts.
For events without ISR the average m rec resolution (as defined in [10]) for the correct jet pairing in W + W − → qqqq events is 1.7 GeV.For W + W − → qqℓν ℓ events the average m rec resolution is 2.4 GeV, 2.8 GeV and 3.4 GeV in the W + W − → qqeν e , W + W − → qqµν µ and W + W − → qqτ ν τ channels respectively.

Extraction of the W mass and width
The Monte Carlo reweighting technique is used to provide the central results of this paper.The W boson mass and width are measured by directly comparing the reconstructed mass spectra in the data to Monte Carlo mass spectra corresponding to different values of M W and Γ W .A likelihood fit is then used to extract M W and Γ W by determining which Monte Carlo spectrum best describes the data.The Monte Carlo spectra for arbitrary (M W , Γ W ) are obtained using the reweighting technique described in [9].In previous OPAL publications, W + W − Monte Carlo samples generated using only CC03 diagrams were reweighted using the ratio of Breit-Wigner functions.For this analysis, Koralw e + e − → 4f reference samples generated with (M W = 80.33 GeV, Γ W = 2.093 GeV) are reweighted to (M ′ W , Γ ′ W ) using the ratio of Breit-Wigner functions.The reference samples used for the reweighting include only final states which have contributions from diagrams involving the W propagator.This procedure is found to give a good approximation to the more exact treatment of using the full fourfermion matrix elements, introducing a bias of less than 5 MeV in the W + W − → qqeν e channel and less than 2 MeV in the other channels.Using the ratio of Breit-Wigner functions rather than four-fermion matrix elements results in a much faster fit.
The mass spectra for background events are taken from Monte Carlo and are assumed to be independent of M W and Γ W .The main sources of background are e + e − → Z 0 /γ → qq and fourfermion processes.The background mass distributions are normalised to the expected number of background events.The reweighted signal spectra are then normalised such that the total number of signal plus background events corresponds to the observed number of events.This is done separately for the W + W − → qqqq, W + W − → qqeν e , W + W − → qqµν µ and W + W − → qqτ ν τ channels, with the W + W − → qqqq channel split into four and five jet topologies.In addition, the W + W − → qqℓν ℓ channels are divided into four subsamples according to the error on the reconstructed invariant mass, σ rec .This division gives a larger weight to events with reconstructed masses which are known with better precision (i.e.small σ rec ) and reduces the expected statistical uncertainty on the fitted W mass from the W + W − → qqℓν ℓ channels by approximately 5%.In the W + W − → qqqq channels, the width of the reconstructed mass distribution is dominated by the intrinsic width of the W. Consequently a similar subdivision into bins of σ rec does not improve the M W sensitivity.However, jet-pairings which give a large jet-pairing likelihood are more likely to be correct and have a better mass resolution as shown in Figure 2.For this reason, the W + W − → qqqq events reconstructed as four jets are subdivided into four bins of jet-pairing likelihood, resulting in a 7% improvement in the statistical sensitivity to M W .No subdivision is performed for the five jet events.
A binned log-likelihood fit to the m rec distributions of the data is performed in the range m rec > 65 GeV.The log-likelihood function is identical to that used previously [10].A log-likelihood curve is determined separately for each channel.For the W + W − → qqℓν ℓ channels, the results are obtained by adding the log-likelihood curves separately determined from each channel in each bin of σ rec .For the W + W − → qqqq channel, the results are obtained by adding the likelihood curves obtained from the five jet events to the likelihood curves obtained from the four bins of jet-pairing likelihood for events reconstructed as four jets.
Two types of fit are performed.In the one parameter fit, Γ W is constrained by the SM relation [16], and only M W is determined.The results of this fit for each channel are given in Table 2 and displayed in Figure 3.The combined result is discussed in Section 6.In the two parameter fit, both M W and Γ W are determined simultaneously.
In the reweighting method the fitted parameters are expected to be unbiased since any offsets in the reconstructed mass introduced in the analysis are implicitly accounted for in the Monte Carlo reconstructed mass spectra used in the reweighting procedure.This is verified using several Monte Carlo samples generated at various M W and Γ W .In addition, tests using a large ensemble of Monte Carlo subsamples, each corresponding to 183 pb −1 and including background contributions, are used to verify for each channel separately and for all channels combined, that the measured fit errors accurately reflect the root-mean-squared spread of the residual distribution for both the M W and Γ W fits. Since M W in the Monte Carlo corresponds to the running width definition, so does the fitted mass.The expected statistical error on the W mass from the combination of the qqℓν ℓ and qqqq channels is 76 ± 1 MeV, where the weights given to the two channels are determined by the both statistical and systematic error contributions.The quoted uncertainty on the expected statistical error is from Monte Carlo statistics.

Breit-Wigner fit
The Breit-Wigner method is analogous to that described in [10].It employs an unbinned maximumlikelihood fit to the reconstructed mass spectrum using an analytic Breit-Wigner function to describe the mass spectrum from e + e − → W + W − .Due to initial-state radiation, the reconstructed mass spectrum is asymmetric.For the W + W − → qqℓν ℓ channels a relativistic Breit-Wigner function, with different widths above and below the peak, gives a satisfactory description of the m rec lineshape.The fitting function used is where Γ +(−) is the width assumed for all m rec above (below) the peak centred at m 0 .This empirical choice of fitting function provides an adequate description for samples up to ten times the integrated luminosity of the data.The widths, Γ +(−) , are fixed to values determined from fits to W + W − signal Monte Carlo samples.Different widths are obtained for each W + W − → qqℓν ℓ decay channel.The shapes of the background distributions and the background fractions are determined from Monte Carlo.The background fractions are held constant in the fit.The fit is restricted to the range 70 < m rec < 88 GeV.
In the W + W − → qqqq channel events are divided into four and five jet samples which are fitted separately.The division is made on the basis of the 5-C kinematic fit probability.Events are reconstructed as five jets if any of the ten possible jet-pairings in the five jet assignment gives a fit probability which is greater than twice that of highest probability of the three jet-pairings for the four jet hypothesis.For both four jet and five jet samples, the fitting function S(m rec ) is multiplied by G(m rec ) = exp[−(m 0 − m rec ) 2 /2σ 2 ].This empirical choice provides a good description of the reconstructed mass spectra in the fit range for W + W − → qqqq events.The value of σ is determined from Monte Carlo.For the four jet W + W − → qqqq sample either one or two jet-pairings are used following the procedure described in [9].In the five jet sample, a jet-pairing likelihood is used.Here the variables used in the likelihood are those used for the reweighting fit described in Section 4.1, with the exception that the 5-C fit mass is not included.
In contrast to the procedure employed for the reweighting method, the W + W − → qqℓν ℓ events are not divided into subsamples according to σ rec nor are the W + W − → qqqq events divided into subsamples according to the jet-pairing likelihood.However, W + W − → qqτ ν τ events are sub-divided into fully-leptonic and semi-leptonic decays of the tau lepton.
The fitted mass, m 0 , must be corrected for offsets not accounted for in the fit, e.g. from initial-state radiation and event selection.The correction is determined using fully simulated Monte Carlo samples generated at different values of M W with the expected background contributions included and is found to depend linearly on M W .The results from the W + W − → qqqq and W + W − → qqℓν ℓ channels, after correction, are given in Table 3.The expected statistical error on the combined W mass measurement is 78 ± 2 MeV.

Convolution fit
The convolution method [12,16] attempts to exploit all available information by constructing a likelihood curve for each selected event.The likelihood is calculated using a functional form where p s is the probability of a candidate event being a true signal event and P s is the probability function for M W given the observed reconstructed mass, m rec .The function P s (M W , m rec ) is defined as where BW(M W , m, s ′ ) is the relativistic Breit-Wigner function for producing off-shell W bosons with mass m for a W mass of M W including the effects of phase space.Γ W is fixed relative to M W using Equation 1.The radiator function, ISR(s, s ′ ), is used to account for the effects of ISR (on a statistical basis).The inclusion of this term accounts for the fact that in the kinematic fit the observed energy is constrained to √ s, where in reality the effective centre-of-mass energy, √ s ′ , depends on the amount of ISR.Finally, R(m, m rec ) is the probability density function relating the experimentally reconstructed mass, m rec , to the average of the two true masses, m, of the off-shell W bosons in the event.
In the W + W − → qqℓν ℓ channels, the full error information on the fitted mass of each event is used in an unbinned maximum likelihood fit.In a previous OPAL publication [10], the event probability function, R(m, m rec ), was assumed to be Gaussian with its central value and standard deviation taken from a 2-C kinematic fit.In this analysis, the 2-C fit is modified, replacing the equal W mass constraint by one in which the masses are individually constrained to a given W boson mass, m.The χ 2 distribution from this fit is converted into the event probability density function, R(m, m rec ), which is non-Gaussian.
In the W + W − → qqqq channel, all events are forced into a five jet topology.A relative likelihood, constructed from the difference between the reconstructed masses in a 4-C kinematic fit, is used to select on average about three jet-pairings per event thereby reducing significantly the combinatorial background.In the Monte Carlo 92% of selected events the correct jet-pairing combination is used in the fit.The convolution is performed in two dimensions, namely the reconstructed masses of the two W bosons, with joint probability density functions corresponding to each jet-pairing obtained from the 4-C kinematic fit.The probability density functions for the different combinations are added to form an event probability density function [12].
In both channels the log-likelihood curves from each selected event are summed to yield a single curve from which a fitted mass is determined.The fitted mass is corrected for offsets not accounted for in the fit in the manner described for the Breit-Wigner fit.The results from the W + W − → qqqq and W + W − → qqℓν ℓ channels, after all corrections, are given in Table 3.The expected statistical error on the combined W mass measurement is 72 ± 2 MeV.The total combined expected error, statistical and systematic, from the convolution method is slightly smaller than that obtained from the reweighting fit.However, at present the reweighting fit is the only method in which a simultaneous fit for (M W , Γ W ) is implemented.For this reason, and for consistency with previous OPAL publications, the reweighting method is retained for the results given in this publication.

Systematic Uncertainties
The systematic uncertainties for the measurement of M W using the reweighting method are estimated as described below and summarised in Table 4.The contributions from each of the sources are added in quadrature to yield the total systematic uncertainty.All contributions were evaluated for both the M W and Γ W determinations.For the alternative analyses, the systematics are estimated similarly and yield comparable results.The determination of the individual systematic uncertainties is described below.

Beam energy
The average LEP beam energy is currently known with a precision of 20 MeV [19] which leads to a systematic uncertainty of 16 MeV on M W and 2 MeV on Γ W .The RMS spread in the LEP centre-ofmass energy of 237 ± 12 MeV [19] and the average difference between the electron and positron beam energies results in negligible impact on the measurement of M W .The Monte Carlo does not include the intrinsic beam energy spread.This introduces a bias of +10 MeV to the value of Γ W obtained from the reweighting fit.A −10 MeV correction is applied and the size of the correction taken as an additional systematic uncertainty on the measurement Γ W .

Initial state radiation
The main Koralw Monte Carlo samples include an O(α 3 ) treatment of initial state radiation (ISR).The systematic error associated with ISR is estimated by reweighting the generated Koralw events to correspond to an O(α 2 ) or O(α) treatment of ISR using the matrix elements calculated in Koralw.The resulting shifts from the comparison of O(α 3 ) to O(α) are less than 2 MeV and are neglected.

Full O(α) corrections
In the process e + e − → 4f γ ISR diagrams dominate.The accuracy of the treatment of ISR in the Koralw generator is sufficient for this analysis.However the Koralw generator does not include all O(α) effects.For example, interference between ISR and FSR graphs and photon radiation from the W bosons are not implemented.Recent calculations using the double-pole approximation indicate that possible mass biases of order 10 MeV could arise [29].However, these estimates were based on generator level studies which compared mass distributions from the full O(α) treatment with the corresponding Born level treatment which does not include any final state radiation.To obtain a proper estimate of the potential mass bias when applied to the experimental measurement it would be necessary to use the full OPAL detector simulation and mass reconstruction procedure.Currently this is not possible since the necessary Monte Carlo programs are not yet available.For this paper possible systematic biases from the full O(α) treatment of photon radiation are not included since the size of the effect is significantly smaller than the current statistical precision of the M W measurement.

Hadronisation in W ± → qq decays
The Koralw Monte Carlo samples taken as references in the extraction of M W use the Lund string model as implemented in Jetset for the simulation of the hadronisation2 of W ± → qq decays.The Jetset model was tuned to OPAL hadronic data recorded at the Z 0 resonance [30].The model was tuned to describe event shape variables and inclusive particle production rates.To assess a possible systematic error from the uncertainties in the tuning of the Jetset model, the Jetset parameters σ q , b, Λ QCD and Q 0 are each varied in turn by ±1σ about their tuned values at the Z 0 [30].The maximum resulting biases to the fitted values of M W from these high statistics samples are 30 MeV and 20 MeV in the qqℓν ℓ and qqqq channels respectively.In addition, an earlier OPAL tune of the Jetset model is used.This version was tuned to describe global event shape variables alone [31].The resulting differences in the fitted value of M W using the earlier OPAL tune are (−33 ± 14) MeV and (−30 ± 12) MeV in the W + W − → qqqq and W + W − → qqℓν ℓ channels respectively.
As a further test of hadronisation uncertainties, the Jetset string model is replaced by the Herwig cluster model.For this comparison, two samples were generated using the same W + W − final states and differing only in the hadronisation modelling.The resulting fitted mass values from the reweighting method from Jetset and Herwig are compared.Differences (Jetset-Herwig) of (−18 ± 14) MeV and (−11 ± 12) MeV are obtained for the W + W − → qqℓν ℓ and W + W − → qqqq channels.A similar comparison is made using the Ariadne colour dipole model [32].Differences (Jetset-Ariadne) of (+15 ± 14) MeV and (+6 ± 12) MeV are obtained for the W + W − → qqℓν ℓ and W + W − → qqqq channels.
On the basis of these comparisons systematic uncertainties of 30 MeV are assigned to the W mass measurements from both W + W − → qqℓν ℓ and W + W − → qqqq.The hadronisation uncertainty is taken to be fully correlated between the two channels.For the measurement of Γ W the corresponding systematic errors in the W + W − → qqℓν ℓ and W + W − → qqqq channels are 25 MeV and 60 MeV respectively.

Detector calibration and simulation
The effects of detector mis-calibrations and deficiencies in the Monte Carlo are investigated by varying the jet and lepton energy scales over reasonable ranges.The ranges used for the systematic variations are dependent on polar angle and are determined from detailed comparisons of data and Monte Carlo utilising both data recorded at 189 GeV and 3.1 pb −1 of data collected at √ s ≈ M Z during 1998.From e + e − → e + e − events the electro-magnetic calorimeter energy scale is known to 0.3%.Similarly from e + e − → Z 0 → µ + µ − events, the uncertainty on momentum scale from the tracking detectors is determined to be 0.3%.These scale uncertainties dominate the detector related systematics in the W + W − → qqeν e and W + W − → qqµν µ channels, and are 14 MeV and 16 MeV respectively.Uncertainties on energy (momentum) resolution and angular resolution are also evaluated but have much smaller impact.Uncertainties on the jet energy scale are determined from e + e − → Z 0 → qq events to be 0.2%-1.0%depending on the polar angle.The dominant effect of the jet energy scale uncertainty is a 6 MeV uncertainty in the W + W − → qqτ ν τ channel.In the W + W − → qqqq channel the sensitivity to the overall jet energy scale is greatly reduced due to the kinematic fit.Here the dominant effect is from the global uncertainty in the jet polar angle scale, i.e. the relative length-towidth scale of the detector, resulting in a 3 MeV uncertainty on the W mass.The linearity of the jet energy scale is checked using e + e − → Z 0 → qq events with a clear three-jet topology.For the measurement of Γ W the main systematic error arises from uncertainties in the jet energy scale and resolution and is 47 MeV in both W + W − → qqℓν ℓ and W + W − → qqqq channels.

Four-fermion modelling
Possible systematic effects associated with the modelling of the four-fermion final state, including interference between W + W − diagrams and other four-fermion processes, are investigated by comparing the grc4f and the Excalibur [33] generators.The comparison is performed using the ratio of the Excalibur to grc4f four-fermion matrix elements to reweight the Koralw Monte Carlo sample.This procedure avoids introducing an additional statistical error associated with generating additional Monte Carlo.The resulting shift in the reconstructed mass is 1 ± 1 MeV in the W + W − → qqℓν ℓ channel and 0 ± 1 MeV in the W + W − → qqqq channel.Similar results are obtained when the Excalibur Monte Carlo is reweighted using the ratio of the grc4f to Excalibur four-fermion matrix elements.Consequently no additional systematic uncertainty is assigned.The corresponding systematic uncertainties on the measurement of Γ W , of order 30 MeV, are determined by the statistical precision of the comparisons.In the W + W − → qqeν e channel an additional 5 MeV uncertainty is assigned, associated with using the Breit-Wigner functions in the reweighting fit.

Background treatment
Uncertainties in the modelling of background processes mainly arise from the simulation of the Z 0 /γ → qq background in the W + W − → qqτ ν τ and W + W − → qqqq channels.A number of systematic checks of the normalisation and shape of the Z 0 /γ → qq background are performed.The background normalisation in the reweighting fit is varied by one standard deviation, using the uncertainties evaluated in [25].The Herwig cluster model is used as an alternative to the Lund string model used in Pythia.In addition, for the W + W − → qqqq channel, data taken at √ s ≈ M Z , scaled by (189 GeV/M Z ), are also substituted for the Z 0 /γ → qq background.The sum in quadrature of the above systematic shifts leads to the assignment of a background systematic error of 6 MeV in the W + W − → qqqq channel and 8 MeV in the W + W − → qqℓν ℓ channel.The systematic error on Γ W due to uncertainties in background modelling is 25 MeV and 41 MeV in the W + W − → qqℓν ℓ and W + W − → qqqq channels respectively arising from differences between Jetset and Herwig and uncertainties on the accepted background cross section.

Monte Carlo statistics
The finite statistics of the Monte Carlo samples used in the reweighting procedure result in an additional uncertainty of 10 MeV on the W mass determined separately in the W + W − → qqqq and W + W − → qqℓν ℓ channels.The corresponding systematic uncertainties for the W width measurement are 51 MeV.

Final state interactions (FSI)
A significant bias to the apparent W mass measured in the W + W − → qqqq channel could arise if the hadronisation of the two W bosons is not independent [16,34].At LEP2 energies, the decay length for W bosons is about 0.1 fm, significantly less than the typical hadronisation scale, approximately 1 fm.The hadronisation of the two W bosons in the W + W − → qqqq channel may have space-time overlap allowing possible final state interactions between the decay products of different W bosons.Two possible sources of final state interactions, both leading to non-independent hadronisation, have been widely considered: Colour Reconnection (CR) and Bose-Einstein Correlations (BEC).Whether these effects play a significant role at LEP2 has yet to be experimentally determined.As a result, to assign a systematic error related to possible final state interactions it is necessary to rely on phenomenological models.The strategy adopted here is to consider only those models for CR and BEC which are consistent with both LEP1 and LEP2 data.
The CR systematic error refers to possible biases in the W mass measurement arising from QCD reconnection effects in the non-perturbative phase.Reconnection effects in the perturbative phase, i.e. multiple gluon exchange between quarks from different W bosons, have been shown to be small [35].To investigate the systematic biases originating from non-perturbative CR effects, the models implemented in the Pythia and Ariadne [32] Monte Carlo generators are studied.In both cases the results from the mass fit to a Monte Carlo sample including CR effects are compared to the results obtained from the corresponding sample without colour reconnection.The results are summarised in Table 5.The largest biases are seen in the Ariadne 2 and Ariadne 3 models3 .However, these models have been strongly disfavoured by studies of three jet events in LEP1 Z 0 data [37].For this reason the current Ariadne implementation of CR is not used to assign a systematic uncertainty.Of the remaining models considered, the SK I model [35] produces the largest bias.The SK Monte Carlo samples use the Jetset tune of [30], where the parton shower cut-off, which directly influences the behaviour of the CR model [35], was set to 1.9 GeV.The strength parameter k i (also called ρ) of the SK I model was set to 0.9, corresponding to 35% of events being reconnected.A systematic uncertainty of 66 MeV on M W is assigned.The corresponding uncertainty on Γ W is 68 MeV.
As is the case in Z 0 decays, BEC between like-sign charged pions are observed in W ± decays at LEP2 [38].These correlations are not implemented in the standard Monte Carlo programs used in this analysis.Correlations between pions from the same W do not affect significantly the global jet properties and are unlikely to bias the M W measurement.If BEC exist between pions from different W bosons in W + W − → qqqq events, significant biases to the reconstructed value of M W may arise.To investigate possible systematic biases from BEC, a sample of W + W − → qqqq events is produced using the Pythia Monte Carlo generator which includes BEC as described in [39].For the evaluation of the systematic error the BEC model BE 32 of [39] with λ = 1.0 and R = 0.42 fm is used.The difference in the fitted W mass between a sample generated with BEC between different W bosons and a sample generated with BEC only between particles from the same W is +67 ± 14 MeV (different-same).The resulting systematic uncertainties on M W and Γ W from BEC in the W + W − → qqqq channel are 67 MeV and 39 MeV respectively.

Fit procedure
A comparison of the results from the reweighting method with those from the Breit-Wigner and convolution fits, summarised in Table 3, is used to test for residual biases in the reweighting method.The expected RMS of the differences of the results from the three methods are determined using an ensemble of Monte Carlo subsamples including background, each corresponding to an integrated luminosity of 183 pb −1 .For each subsample the difference in the fitted M W determined using the reweighting method and that determined using each of the other methods is calculated for the W + W − → qqqq and W + W − → qqℓν ℓ channels.These differences have mean values consistent with zero and RMS values of approximately 54 MeV and 41 MeV (62 MeV and 84 MeV) when comparing the reweighting fits to the Breit-Wigner (convolution) fits in the W + W − → qqℓν ℓ and W + W − → qqqq samples, respectively.Since the alternative analyses yield results consistent with those obtained using the default reweighting analysis, no additional systematic error is assigned.
The correlation coefficient between M W and Γ W is 0.04.For this fit, the central values are determined by adding the log-likelihood curves from the W + W − → qqqq and W + W − → qqℓν ℓ channels including the effects of systematic uncertainties.The likelihood contours for this fit are displayed in Figure 4.
Due to the larger systematic errors in the W + W − → qqqq channel, mainly from FSI, it carries a reduced weight of 0.36 in the combination.The combined W + W − → qqqq and W + W − → qqℓν ℓ results from the alternative analyses, after all corrections, are for the Breit-Wigner fit, M W = 80.436 ± 0.078(stat.)± 0.050(syst.)GeV, and for the convolution fit, M W = 80.378 ± 0.073(stat.)± 0.051(syst.)GeV.
where the uncertainties from the LEP beam energy (lep) and from final state interactions (fsi) are quoted separately.In the combination the systematic uncertainties for the colour reconnection and Bose-Einstein correlations of the previous data were set to those presented in this paper.
The difference between the fitted M W in the W + W − → qqqq and W + W − → qqℓν ℓ channels is ∆M W ≡ (M qqqq W − M qqℓν ℓ W ) = −0.032± 0.127(stat.)± 0.047(syst.)GeV.A significant non-zero value for ∆M W could indicate that FSI effects are biasing the value of M W determined from W + W − → qqqq events.The systematic error on the quoted value of ∆M W does not include contributions from CR/BEC effects.For the ∆M W measurement the hadronization systematic uncertainties in the W + W − → qqqq and W + W − → qqℓν ℓ channels are taken to be uncorrelated; a conservative choice in the absence of a good estimate of the true correlation coefficient.
The M W results from direct reconstruction in the W + W − → qqℓν ℓ and W + W − → qqqq channels are combined with the OPAL measurement of M W from the W + W − production cross-section at √ s ≈ 161 GeV [8].The combination is made assuming that the mass measurements from direct reconstruction and from the threshold cross-section are uncorrelated, apart from the uncertainty associated with the LEP beam energy, which is taken to be fully correlated.The direct reconstruction measurements from the W + W − → qqℓν ℓ and W + W − → qqqq channels from three different centre-of-mass energies are combined using the full covariance matrix.The combined result is M W = 80.432 ± 0.066(stat.)± 0.032(syst.)± 0.028(fsi) ± 0.017(lep) GeV.
The measurement of Γ W presented in this paper is combined with the previous OPAL measurements from direct reconstruction [9,10] to give Γ W = 2.04 ± 0.16(stat.)± 0.09(syst.)GeV.
For the earlier data, with relatively low statistics, the measured statistical error is correlated with the measured width.Consequently, to avoid biasing the combination, the separate width measurements are weighted using the expected statistical error rather than the measured statistical error.

Summary
Using the 183 pb −1 of data recorded by the OPAL detector at a mean centre-of-mass energy of approximately 189 GeV, a total of 2088 W + W − → qqqq and W + W − → qqℓν ℓ candidate events are used in a fit constraining Γ W using Equation 1 to obtain a direct measurement of the W boson mass, M W = 80.451 ± 0.076(stat.)± 0.049(syst.)GeV, while a second fit is used to determine directly the width of the W boson, Γ W = 2.09 ± 0.18(stat.)± 0.09(syst.)GeV.

Figure 1 :Figure 2 :Figure 3 :
Figure 1: The reconstructed mass, m rec , distributions for W + W − → qqqq events divided into four and five jet samples.The top two plots show the distributions before the jet-pairing likelihood (JPLH) cut.The lower two plots show the mass distributions which are used to determine M W , i.e. after the cut on JPLH.The points correspond to the OPAL data and the histograms to the Monte Carlo predictions.The contribution from the non-WW background is shown as the cross-hatched histogram and the combinatorial background is indicated by the singly-hatched histogram.

Table 1 :
[41]13,40]t is consistent with both the other direct measurements and the indirect value inferred from fits to electroweak data[3].It should be noted that the systematic error is dominated by the uncertainties arising from hadronisation and final state interactions.For future LEP combinations these uncertainties will need to be reduced to benefit from the ultimate LEP combined statistical uncertainty of about 25 MeV.This measurement of Γ W is in good agreement with the Standard Model expectation and with other direct measurements at LEP[12,13,40]and at the Tevatron[41].Numbers of events used in the W mass and width determination for each channel and all channels combined.Only events surviving the cuts described in Section 4.1 are included.The W + W − → qqτ ν τ numbers include some events initially selected as either W + W − → qqeν e or W + W − → qqµν µ .
The result for the W boson width is combined with the previous OPAL result from data recorded at √ s ≈ 172 − 183 GeV to obtain Γ W = 2.04 ± 0.16(stat.)± 0.09(syst.)GeV.

Table 2 :
Results using the reweighting method for the fit from 183 pb −1 of data taken at √ s ≈ 189 GeV for each of the channels separately and for the combined W + W − → qqℓν ℓ channel.Equation 1 is used to constrain Γ W to M W .The expected errors are estimated using an ensemble of Monte Carlo subsamples.The errors obtained in the data are consistent with the RMS spread of the errors obtained from the Monte Carlo subsamples.W /GeV Measured M W /GeV Measured M W /GeV W + W − →qqℓν ℓ 80.48 ± 0.10 ± 0.04 80.47 ± 0.11 ± 0.04 80.38 ± 0.11 ± 0.04 W + W − → qqqq 80.40 ± 0.10 ± 0.10 80.38 ± 0.11 ± 0.10 80.37 ± 0.09 ± 0.09

Table 3 :
Comparison of M W fit results obtained using the three different fitting techniques for the W + W − → qqℓν ℓ and W + W − → qqqq channels separately.In each case the first uncertainty is statistical and the second systematic.The observed differences are compatible with the expected RMS differences between the results from the different fits.

Table 4 :
Summary of the systematic uncertainties for the fit results.For the fits to determine M W , Γ W is constrained to its Standard Model relation.The uncertainties are given separately for fits to the W + W − → qqqq and W + W − → qqℓν ℓ samples.Uncertainties from ISR give negligible contributions to both M W and Γ W .The uncertainty related to the Monte Carlo models of the four-fermion process give is negligible for the measurement of M W .

Table 5 :
The predicted shifts in the fitted mass and width from the W + W − → qqqq channel for various colour-reconnection (CR) models.The numbers correspond to the fitted mass (width) using the reweighting method in the reconnected sample minus the fitted mass (width) for the corresponding sample without CR.The errors are statistical only.