Measurement of the Mass and Width of the W Boson in e+e- Collisions at sqrt(s) = 189 GeV

A measurement of the W mass and width has been performed by the DELPHI collaboration using the data collected during 1998. The data sample has an integrated luminosity of 155 pb^{-1} and an average centre-of-mass energy of 188.6 GeV. Results are obtained by applying the method of direct reconstruction of the mass of the W from its decay products in both the W+W- ->lvqq and W+W- ->qqqq channels. The W mass result for the 1998 data set is M_W = 80.387 +/- 0.087 (stat) +/- 0.034 (syst) +/- 0.017 (LEP) +/- 0.035 (FSI) GeV/c2, where FSI represents the uncertainty due to final state interaction effects in the qqqq channel, and LEP represents that arising from the knowledge of the beam energy of the accelerator. Combining this result with those previously published by the DELPHI collaboration gives the result M_W = 80.359 +/- 0.074 (stat) +/- 0.032 (syst) +/- 0.017 (LEP) +/- 0.033 (FSI) GeV/c2. The combined value for the W width is \Gamma_W = 2.266 +/- 0.176 (stat) +/- 0.056 (syst) +/- 0.052 (FSI) GeV/c2.

where F SI represents the uncertainty due to final state interaction effects in the qq ′ qq ′ channel, and LEP represents that arising from the knowledge of the beam energy of the accelerator.Combining this result with those previously published by the Delphi collaboration gives the result M W = 80.359 ± 0.074(stat) ± 0.032(syst) ± 0.017(LEP ) ± 0.033(F SI) GeV/c 2 .
(Accepted by Physics Letters B)

Introduction
The W mass has been measured by the Delphi collaboration using the data collected during 1998.This data sample has allowed a significant improvement in the accuracy of the collaboration's W mass determination as the integrated luminosity is more than twice that on which previous Delphi results are based [1,2,3].The W mass has also been determined by the other Lep collaborations [4] and at hadron colliders [5].Using the same reconstruction method as for the W mass, results on the direct measurement of the W width are also obtained in this paper, and can be compared with those of the other Lep collaborations and of the Cdf collaboration [6].
Section 2 of this paper describes the characteristics of the 1998 data sample and of the event generators used in this analysis.
The analysis is performed through the direct reconstruction of the mass of the W from its decay products in the W + W − → qq ′ qq ′ (fully-hadronic) and W + W − → ℓν ℓ qq ′ (semi-leptonic) decay channels.The applied methods are described in Sec. 3, and have been refined from those in previous publications.Results are now reported for τ ν τ qq ′ events, and in both the fully-hadronic and semi-leptonic channels improvements in the handling of initial-state radiation (ISR) have been included.
The systematic error evaluation described in Sec. 4 has increased in sophistication from that previously reported.New techniques have been applied, such as the use of mixed Lorentz boosted Z's, and a wider range of higher precision simulation studies have been performed.
The results of this analysis are reported in Sec. 5, and are combined with the previous Delphi results.

Data and Simulation Samples 2.1 Data
A detailed description of the Delphi apparatus and its performance can be found in [7] 1 .In 1998 the detector was used to record data at the Z peak and at a nominal centre-of-mass energy of 189 GeV.
The Z peak data were recorded before (1.8 pb −1 ) and towards the end (0.6 pb −1 ) of the high energy data, thus facilitating checks of detector stability.These data were the principal sample used for calibration and alignment of the detector and, in this analysis, assist the study of systematic uncertainties.
The average centre-of-mass collision energy for the high energy data was 188.6 GeV.The luminosity weighted r.m.s. of this value, assessed on a fill by fill basis, was 50 MeV.In the data sample considered for analysis all the detectors essential for this measurement were required to be fully efficient; the operation of the central tracking detectors was important for all decay channels, in the ℓν ℓ qq ′ analysis stricter requirements than in the qq ′ qq ′ channel were placed on the electromagnetic calorimeters.The selected samples correspond to integrated luminosities of 152.9 pb −1 for the ℓν ℓ qq ′ analysis and 157.4 pb −1 for the qq ′ qq ′ decay channel.

Simulation
The response of the detector to various physical processes was modelled using the simulation program DELSIM [7], which includes modelling of the resolution, granularity and efficiency of the detector components.In addition, detector correction factors, described in Sec. 4, were included to improve the description of jets, electrons and muons.For systematic uncertainty studies a fast simulation program, relying on a relatively simple set of smearing and efficiency parametrisations, was also used.
WW events and all other four-fermion processes were produced using the event generator EXCALIBUR [8], with initial-state radiation described using the QEDPS program [9].The W mass (M W ) and width (Γ W ) definition used throughout this paper correspond to a W propagator with an s-dependent width.The background process e + e − → q q(γ) was simulated with the PYTHIA 5.7 [10] event generator.Two photon backgrounds were studied using the TWOGAM generator [11].The fragmentation of all events was performed using JETSET 7.4 [10] tuned to the Delphi Lep1 data [12].Systematic error checks were performed using other generators and variations in the fragmentation tuning as described in Sec. 4. The systematics section also reports Delphi results for the common Lep event samples produced in the context of the Lep WW Workshops: the event generation was performed by Aleph using KORALW 1.21 [13], and these events were then passed through the DELSIM program.

Semi-Leptonic Decay Channel
The fitting procedure presented here is a development of that described in [3] for the eν e qq ′ and µν µ qq ′ channels, where the fitting function now includes a description of the ISR spectrum of WW events.In addition, we present here an analysis of the τ ν τ qq ′ channel and the event selection in all semi-leptonic channels is now based on a neural network.

Lepton Identification
Charged particles were identified as muons if they were associated with a hit in the muon chambers, or had an energy deposit in the hadron calorimeter that was consistent with a minimum ionising particle.Muon identification was performed in the polar angle range between 10 • and 170 • .
Electron identification was performed in the polar angle range between 12 • and 168 • by selecting charged particles with a characteristic energy deposition in the electromagnetic calorimeters.In the central region of the detector covered by the HPC electromagnetic calorimeter, electrons were selected using the energy over momentum (E/p) ratio of the candidate.For lower energy candidates (below 30 GeV) this was supplemented by selection criteria on the shape of the calorimetric shower and a more stringent comparison of the track extrapolation and calorimetric shower positions, while for electrons of higher energy, negligible energy deposition in the hadron calorimeter was required.In the polar angle region below 36 • and above 144 • , where the momentum resolution is poorer, tracks associated to electromagnetic energy showers above 8 GeV and negligible hadron calorimeter energy deposits were considered as electrons.
Tau candidates were obtained by clustering the events into a three-jet configuration using the LUCLUS [14] algorithm.The jet with the lowest charged multiplicity was chosen as the tau candidate.As the tau lepton predominantly decays into a final state with one or three charged particles, only jets containing between one and four charged tracks were selected.

Selection
The event selection was based upon a multi-layer perceptron neural network [15].The network was optimised separately for eν e qq ′ , µν µ qq ′ , τ ν τ qq ′ candidates containing only one charged particle and τ ν τ qq ′ candidates with several charged particles.
Having removed the lepton candidate in eν e qq ′ and µν µ qq ′ events, the LUCLUS jet clusterization algorithm (with a d join of 7.5 GeV/c) was used to cluster the remaining particles.Events containing more than three jets were re-clustered, forcing them into a three-jet configuration.τ ν τ qq ′ events were clustered as the tau candidate and a two-jet system.The events were reconstructed using a constrained fit imposing conservation of four-momentum and equality of the two W masses in the event.As the energy of the tau lepton is unknown, due to the emission of at least one neutrino in its decay, the mass in the τ ν τ qq ′ channel is entirely determined by the jet system.
The neural network relied upon the characteristic event properties in each decay channel.The input variables included lepton momentum, polar angle of the missing momentum, and the isolation of the lepton candidate from the hadronic system of the event.The electron and muon identification was obtained as a strong or loose tag and the network was optimised separately in each decay channel for these two categories of identification.
The network was tuned on samples of signal and background simulation events, and its performance estimated on independent samples of events.After applying a cut on the network output the selection performance is as shown in Tab. 1. Events that passed the cut in the muon channel were selected, the remaining events were considered as electron channel candidates and, if they were again rejected, were then analysed under the tau channel hypothesis.The neural network output in the τ ν τ qq ′ channel is shown in Fig. 1.
The tau selection sample contains a significant proportion of other semi-leptonic decays: the composition was estimated from simulation to be 67% τ ν τ qq ′ , 19% eν e qq ′ , 5% µν µ qq ′ with the remaining fraction dominated by qq(γ) background events.This corresponds to a 39% selection efficiency for τ ν τ qq ′ events.Further information on a selection of W + W − events in Delphi with a similar performance is available in [16].
The fraction of semi-leptonic WW events in the sample was extracted from simulation as a function of the neural network output: this is referred to below as the event purity P e .This feature is particularly useful for the tau selection, where the proportion of background events is highest.

Likelihood Function
The following likelihood function was evaluated for each selected events with a reconstructed mass in the range 68 − 92 GeV/c 2 : where P e is the event purity, discussed above, S ′′ is the signal function that describes the reconstructed mass distribution of the semi-leptonic W decays, and B is used to describe background processes.The reconstructed event mass m f it and its estimated error σ f it are both obtained from the constrained fit.The distribution of background events is extracted from simulation as a function of m f it .
The signal function S ′′ is defined in terms of S,S ′ as discussed below.The function S relies on the convolution of three components, using x and m as the dummy integration variables: (2) BW P S is a phase-space corrected relativistic Breit-Wigner distribution (representing the W mass distribution) which is convoluted with the Gaussian function G describing the detector resolution.The width of the Gaussian depends upon the reconstructed mass error obtained in the constrained fit for that event.Details of the BW P S and G terms are given in [2].A recent addition to the analysis is the description of the ISR spectrum, which is parametrised as where x γ is the ratio of the photon energy to the centre-of-mass energy and β is calculated from the electromagnetic constant (α), the centre-of-mass energy squared (s) and the electron mass (m e ): Including this ISR term decreases the bias on the fitted W mass by approximately 400 MeV and improves the expected error by 2 ± 1%.
The event selection contains a significant fraction of τ ν τ qq ′ events in the electron and muon channel samples, and of eν e qq ′ events in the tau sample (see Tab. 1).In the tau channel the mass of the event is determined from the jet system.The behaviour of true τ ν τ qq ′ and eν e qq ′ events in this fit are found to be similar.However, in the electron and muon channel samples the behaviour of the τ ν τ qq ′ events is somewhat different to that of the eν e qq ′ , µν µ qq ′ events.The τ ν τ qq ′ events have a worse mass resolution and a small negative bias on the mass.The fraction of tau events, which have been wrongly classified and are contained in the electron and muon channel samples, has been parametrised in bins of the lepton energy and the measured missing mass.This event impurity P τ e was then taken into account in the likelihood function for the electron and muon samples, by defining the signal function S ′′ as where S ′ is analogous to S, but with the width of the Gaussian resolution function increased according to simulation studies.All remaining biases in the analysis due to using this approximate likelihood description are corrected for in the calibration procedure as described in Sec.3.3.

Fully-Hadronic Decay Channel
The analysis of the fully-hadronic channel was based on that applied in [3].However, the implementation now relies on kinematic fits with four rather than six constraints and includes a new ISR treatment.

Event Selection
A sample of hadronic events was selected by requiring more than 13 charged particles and a total visible energy exceeding 1.15 E BEAM .
q q(γ) events were suppressed by demanding an effective centre-of-mass energy [17], after ISR emission, of greater than 161 GeV.The algorithm for assessing the e + e − collision energy considers both the emission of unobserved ISR photons in the beam-pipe and photon candidates detected in the electromagnetic calorimeters.
The DURHAM jet clustering algorithm [18] with y cut of 0.002 was applied to the event.If any of the resulting jets contained less than three particles or had an invariant mass smaller than 1 GeV/c 2 , clustering was continued to a higher value of y cut .Events with less than four jets were then rejected, while events containing six or more jets were reclustered into five objects representing four quarks plus one hard gluon jet.
Events containing b-quarks were rejected using the Delphi b-tag package [19], this removes 17% of ZZ events and 6% of q q(γ) background while reducing the signal efficiency by only 0.2% .A four constrained kinematic fit [3] was applied to the remaining events, enforcing conservation of energy and momentum.
A variable to discriminate between q q events with hard gluon radiation and signal events was constructed.This compound variable relied upon the fitted jet energies and the inter-jet angles.The expected fraction of qq ′ qq ′ events in the selected sample, the event purity, was parametrised as a function of this variable.Events with an estimated purity below 25% were rejected.
The performance of the event selection is shown in Tab. 1.Further information on a selection of W + W − events in Delphi with a similar performance is available in [16].

Likelihood Function
For each of the selected events an event likelihood was constructed : As in [3] the signal function S(m x , m y , M W , Γ W ) consists of Breit-Wigner terms for the WW and the ZZ contribution and a phase space correction factor.A flat distribution B accounts for background processes and wrong jet pairings in the signal events.Both S and B were normalised to 1 over the integration area.The fraction of the signal and background likelihoods used for each event depend upon the event purity P ef f e .This purity was parametrised as a function of a discriminating variable as described above.
The sum i w i,e • p i,e (m x , m y ) is a weighted sum of the probability densities p i,e that the event e corresponds to two heavy objects with mass m x and m y .
The probability density p i,e (m x , m y ) ∝ exp − 1 2 χ 2 i,e (m x , m y ) was determined for all jet pairings (three possibilities for a four-jet event and ten for a five-jet event) and with three different clustering algorithms (DURHAM [18], CAMJET [20] and DICLUS [21]).The relative probabilities w i,e that the corresponding jet pairing was the correct one were estimated using jet charge information and the transverse momentum of the gluon candidate (see [3]).The three jet clustering algorithms were given the same weight.
The probability was calculated using a Gaussian approximation for the χ 2 : The masses m fit x , m fit y , their errors σ mx and σ my and the correlation between them, ρ xy , are taken from a four constrained kinematic fit.When the χ 2 4C is larger than the number of degrees of freedom (NDF = 4), the χ 2 i,e (m x , m y ) is rescaled with a factor NDF/χ 2 4C in order to compensate for non-Gaussian resolution effects.This procedure decreases the computing time taken by an order of magnitude compared with the full six constrained fit [3], while resulting in only a minimal reduction in the W mass precision obtained (2 ± 1%).
A new feature of this analysis is a treatment of events under the collinear ISR hypothesis.A kinematic fit was performed with modified constraints to simulate the emission of an ISR photon of momentum p z inside the beam pipe: The probability that the missing momentum in the z direction was indeed due to an unseen ISR photon was extracted from the simulation as a function of |p fit z |/σ pz , where σ pz is the estimated error on the fitted z momentum component; only events with this ratio greater than 1.5 were treated with the mechanism described below.
Additional probability density p i,e terms were then included in the likelihood sum for these events, with a relative weight factor derived from the ISR hypothesis probability.An example of the effect of including the ISR hypothesis is shown in Fig. 2.This treatment was applied to 16% of the events and resulted in an improvement of the expected W mass error for these events of 15%.

Mass and Width Extraction
The distribution of the reconstructed invariant masses of the selected events after applying a kinematic fit, imposing four-momentum conservation and the equality of the two di-jet masses, are shown in Fig. 3.This plot is provided for illustrative purposes only, the mass and width fitting procedure is described below.
The combined likelihood of the data can be obtained from the product of the event likelihoods described above.The W mass and width were extracted from maximum likelihood fits.The W mass fit is performed assuming the standard model value for the W width.The W width was obtained assuming a mass of 80.35 GeV/c 2 .The correlation between M W and Γ W was found to have a negligible impact on the extracted width value.
The mass and width analyses have been calibrated separately in each of the decay channels (qq ′ qq ′ ,eν e qq ′ ,µν µ qq ′ ,τ ν τ qq ′ ).The biases of the analyses were estimated by re-weighting generated simulation samples to obtain the fitted mass and width values.The re-weighting was performed using the extracted matrix element of the EXCALIBUR generator.The linearity of the mass analysis was estimated using independent simulation samples generated at three W mass values, while the re-weighting procedure was used for the width analysis.The analyses were corrected with the calibration results, and the statistical error on the bias is included in the systematic error.
The analyses were checked by performing fits to a large number of samples of simulation events.Each sample was comprised of a mixture of signal and background simulation events to represent the expected distribution in the data.The pull distribution was demonstrated to be compatible with a Gaussian of width one to an accuracy of better than 1%.The mean expected statistical error in the W mass was 262 MeV/c 2 for eν e qq ′ , 203 MeV/c 2 for µν µ qq ′ , 311 MeV/c 2 for τ ν τ qq ′ and 104 MeV/c 2 for the qq ′ qq ′ channel.

Systematic Uncertainties
The sources of systematic error that have been considered for the W mass and width determinations are described in the subsections below.The results of these studies are summarised in Tabs.8 and 9.

Calibration
The accuracy with which the bias of the analysis can be determined is limited by the size of the generated simulation samples.Sufficient events were generated to limit this error to less than 10% of the statistical error in any given channel.The calibration procedure is described in Sec.3.3.

Detector Effects
The data taken at the Z peak were used to study, and limit, possible errors in the detector simulation model.
Muon studies were performed on a selected sample of Z → µ + µ − events.From the dimuon sample corrections to the inverse momentum scale, 1/p, were calculated separately for positive and negative muons as a function of the lepton polar angle.The systematic error on this correction was estimated by varying it by half of its value.The momentum resolution (typically 0.001 in 1/p) was found to be slightly better in simulation than in the data (a maximum difference of 10%).This was corrected by smearing the simulation with a Gaussian.An extra smearing of 0.0005 in 1/p was used to estimate the systematic error coming from this correction.The combined systematic error from these corrections is quoted for the µν µ qq ′ channel as the lepton correction systematic error in Tabs.8 and  9.
The correction of the energy scale of electrons was determined from Bhabha events at the Z peak in different polar angle regions.The residual systematic error on this absolute energy scale was estimated to be 0.5%.In each of these polar angle regions, the energy resolution of simulation events was degraded by applying a Gaussian smearing, and the residual error on this smearing was estimated to be 1%.The dependence of the energy calibration as a function of the electron energy was checked using low energy electrons from Compton events at the Z peak, and high energy electrons from radiative Bhabha scattering at high centre-of-mass energy.In these cases the true energy of the lepton was deduced from 3-body kinematics using only the angular information and assuming that the unseen particle was along the beam axis.The absolute energy calibration was found to be compatible with requiring no additional corrections in all energy ranges.The systematic error coming from this source was estimated assuming a 1% change of slope in the energy calibration in the range of interest (between 25 and 70 GeV).The lepton correction uncertainty for the eν e qq ′ channel, quoted in Tabs.8 and 9, is the quadratic sum of the errors from these three error sources.
The lepton correction error is not quoted for the τ ν τ qq ′ channel as the reconstructed tau lepton carries no information for the reconstructed W mass.
Jet energies were studied in Z → q q events as a function of the polar angle and reconstructed energy of the jet.The comparison of data and simulation showed agreement within a band of ±2% over most of the Delphi detector's angular coverage and an overall uncertainty of 1% was estimated.The simulation was smeared by a Gaussian function to improve the description of the observed energy spread in the data and the residual error on this smearing was estimated to be ±4%.The dependence of the energy calibration as a function of the jet energy was checked using low energy jets from q q + gluon events at the Z peak and high energy jets from radiative Z's at high centre-of-mass energy.The true jet energy was estimated from 3-body kinematics; in the radiative events the unseen photon was assumed to be along the beam axis.No additional energy calibration slope was necessary over the relevant energy range (25 to 75 GeV) and a 1% change in slope was used to calculate the systematic error.A study of the acollinearity of jets in Z → q q events was performed and an appropriate smearing to the simulation of the jet angular direction was estimated.A systematic error was estimated by applying an extra 5 mrad angular smearing.The jet correction uncertainty, quoted in Tabs.8 and 9, is the quadratic sum of these four errors.
A possible source of angular distortion in Delphi is the uncertainty on the length to width ratio of the detector.The detector is aligned relative to the vertex detector, the largest uncertainty being the radius of this detector which is known to a precision of ±0.1%.This error is listed in the M W systematics table under aspect ratio.The corresponding uncertainty for the width measurement is less than 10 MeV.

Background Description
The background level was changed by ±10% in the simulation; this easily covers the expected uncertainty in the accepted cross-section, as discussed in [16].The dominant background source, Z → q q(γ), was generated using both JETSET and HERWIG fragmentation models and mass fits performed using both these samples.The 4-jet rate in data and simulation was also studied using events collected at the Z peak.It is concluded that the background description is a relatively small component of the systematic error for the mass measurement, as shown in Tab. 8 and somewhat more important for the width (see Tab. 9).

Fragmentation
A study of the possible effects on M W and Γ W of the simulation of the event fragmentation was performed, the results are provided in Tabs. 2 and 3 .
Events were produced using a fast simulation package where the value of Λ QCD and σ q were changed with respect to the standard Delphi JETSET tuned values [12] by twice their estimated errors2 .The estimated errors are ±0.018GeV for Λ QCD and ±0.007 GeV for σ q .
The event samples generated by Aleph in the context of the LEP WW workshop are used to compare results from the HERWIG and JETSET fragmentation models.The HERWIG events were produced with a recent tuning [22] which provides a better description of the data than previous HERWIG versions.
The results on the Delphi tuning of JETSET are also compared with those of the Aleph collaboration.This comparison cross-checks several effects as the Aleph events were produced with a different generator (KORALW) from the Delphi events which includes a different ISR and final-state radiation (FSR) treatment.
As the results in Tabs. 2 and 3 are all compatible with zero, we quote a systematic error from fragmentation reflecting twice the statistical precision of the JETSET tuning parameter studies in the combined semi-leptonic and fully-hadronic channels.

Mixed Lorentz Boosted Z's
The agreement between data and simulation can be cross-checked using the method of Mixed Lorentz Boosted Z's.Z events were selected from the Z peak data sample collected during 1998 and the corresponding simulation sample.Through a suitable choice of Lorentz boost and superimposing two Z events a WW event may be emulated.The angular distribution of the Z events used was chosen to match that expected in WW events.
This technique has been applied to the qq ′ qq ′ mass measurement and the hadronic system in semi-leptonic events.The differences between data and simulation obtained are 2 ± 2 MeV/c 2 and 3 ± 10 MeV/c 2 in the two topologies respectively, where the errors are statistical.These results reinforce the view that the quoted systematic errors for detector and fragmentation effects are conservative.A study of possible intrinsic uncertainties of the MLBZ method [23] in the qq ′ qq ′ channel estimates the accuracy of this technique to be 5 MeV/c 2 , and demonstrates excellent agreement between data and simulation as a function of relevant event variables.

ISR
A comparison of two independent models of ISR was performed.WW events produced with the EXCALIBUR generator were re-weighted as a function of the total ISR energy in the event.The weights were obtained for each WW decay channel from a generator level study of the KORALW [13] ISR treatment (based on the YFS exponentiation approach) and that of the standard algorithm used in Delphi QEDPS [9] (based on a parton shower approach).The results of the W mass and width fits are shown in Tabs. 4 and 5.We conservatively choose to quote the largest deviation observed in any of the channels.

LEP Beam Energy
The average Lep beam energy at Delphi is evaluated by the energy working group [24] at 15 minute intervals of running or after significant changes in the beam energy.The measured centre-of-mass energy is imposed as a constraint in the kinematic fit, and hence the relative error on the beam energy translates to the same fractional error on the W mass determination.The spread of energies of the electrons and positrons in the LEP beams was found to have a negligible impact on the mass and width measurements.

Bose-Einstein Correlations
Bose-Einstein statistics dictate that the production amplitude for final state particles should be symmetrical under the exchange of identical bosons.The omission of these correlations between particles from different W bosons in our standard simulation could lead to a systematic error on the W mass and width measurement in the qq ′ qq ′ channel.
A clear picture has yet to emerge from the experimental study of this phenomenon [25].
To evaluate the possible size of the effect on the mass and width measurements we have considered several phenomenological models; the results are given in Tab. 6.The relevant value for the systematic uncertainty is the difference between the shifts obtained from Bose-Einstein correlations inside individual W's and that between W's.The models used are: • The LUBOEI algorithm [26] in JETSET changes the momentum vectors of identical final state bosons to model the two particle correlation function, and then offers a range of options to restore energy and momentum conservation.• Global re-weighting aims to reproduce the enhancement of identical bosons close in phase-space by giving weights to events.This procedure does ensure energy and momentum conservation, but may adversely affect other event distributions.Delphi results were obtained using the Kartvelishvili/Kvatadze/Møller (KKM) re-weighting scheme [27] and reported on in [3].
• In [3] we also reported results for a study based on a modification of the JETSET fragmentation model introducing quantum mechanical 2-particle and 3-particle interference (ST) for identical bosons using a local re-weighting technique [28].
The version of LUBOEI model studied here is BE 32 , in which a local energy and momentum conservation procedure is applied [26].Six sets of fully simulated Z events were generated at a range of values of correlation strength (λ) and radius (r), and the four momentum difference, Q, between all selected same charge particle pairs was calculated.An interpolation was performed on the basis of the Q distribution to obtain λ = 1.35 and r = 0.6 fm (P ARJ(93) = 0.34 GeV) which provide the optimal description of the Delphi Z peak data.The simulated event samples were produced with the standard Delphi JETSET tuning.The W mass and width shifts were evaluated using samples produced with fast simulation.

Colour Reconnection
The hadronization of two W bosons in the qq ′ qq ′ channel may not occur independently.The colour flow between the W's may lead to a shift in the measured W mass.Although experimental work is progressing [29], a suitable sensitivity has not yet been reached by the measurements to limit the effect on M W .
In the previous publication [3] we reported results with the ARIADNE colour reconnection model, the results are repeated in Tab. 7.However, we note that the current version of the ARIADNE 2 model is disfavoured by Lep1 data [30], and that we do not consider the ARIADNE 3 model for the systematic error assessment as it allows perturbative phase reconnection where calculations have shown the effect to be small [31].
The Sjöstrand/Khoze models of colour reconnection are available in the JETSET framework.We have used the SK1 and SK2 models.Results for the SK1 model are quoted for 30% of reconnected events, this is the same reconnected fraction as in the SK2 model.Events were produced using DELSIM and a fast detector simulation.By processing the same event sample through the full detector simulation and the fast simulation, the reliability of the fast detector simulation for this study was clearly demonstrated.The simulated event samples were produced with the standard Delphi JETSET tuning.
However, we report that the standard implementation of SK1 shows numerical instabilities that reduce the reliability of the model used.In this model the reconnection probability is a function of the string overlap.This overlap is calculated by numerical integration through sampling.The accuracy of this calculation is improved by using 1000 sampling points rather than the default value of 100.In addition, we also report results from a more efficient sampling of the string overlap, in which the sampling is performed along strings taking into account their life-time and the total overlap is calculated as the sum of the overlap of pairs of string pieces.
The observed W mass and width shifts are given in Tab.7 and in Fig. 4 the observed W mass and width shifts are shown as functions of the reconnected fraction of events in the SK1 model.
We conservatively choose to quote as a systematic error the largest effect observed in our studies; this is approximately 50 MeV for both the mass and width analyses.

Correlations
The components of the systematic error arising from the jet energy scale corrections, aspect ratio, ISR, fragmentation and LEP beam energy are taken as correlated between the analyses in the different decay channels for the 189 GeV data.The background description and lepton modelling uncertainties are treated as uncorrelated between WW decay channels.
In Sec.5.2 a combination with the previously published Delphi results is performed.The LEP energy correlation matrix is used [24] in this process.The calibration statistics error is uncorrelated between years, while all other systematic errors are conservatively assumed to be fully correlated between years.

189 GeV Results
The W mass and width results of the analyses described in this paper are presented in Tabs. 10 and 11.The error is divided into its statistical component, indicated (stat), the main systematic component (syst) and the systematic from the beam energy uncertainty (LEP ).In the qq ′ qq ′ channel an error from final state interaction effects (F SI) is also included.

Combined DELPHI Results
These results are combined with the previously published Delphi results on the W mass determination by direct reconstruction [2,3] and the result from the threshold crosssection measurement [1] to obtain: M W = 80.359 ± 0.074(stat) ± 0.032(syst) ± 0.017(LEP ) ± 0.033(F SI) GeV/c 2 .
The mass difference between the measurements for the fully-hadronic and semi-leptonic channels, ∆M W (qq ′ qq ′ − ℓν ℓ qq ′ ), has been determined: A significant non-zero value for ∆M W could indicate that FSI effects are biasing the value of M W determined from the qq ′ qq ′ events.As ∆M W is primarily of interest as a check of the possible effects of final state interactions, the errors from FSI effects are set to zero in this determination: all other errors and correlations were as described above.
The combination for the W width with the 1997 data [3]

Figure 1 :
Figure 1: The output distribution of the event selection neural network for τ ν τ qq ′ candidates.The shaded areas indicate the contribution from the various simulated states, the data are shown as points with statistical error bars.The value of the selection cut applied is indicated by the dashed line.Note that the order of the simulation contributions in the figure follows that in the key, the τ ν τ qq ′ signal is the first distribution shown and background from two-photon diagrams the last.

Figure 3 :
Figure 2: An example of the reconstructed probability density function i w i,e • p i,e (m + , m − ) (see text) for the invariant masses in a simulated 4-jet qq ′ qq ′ event without (a) and with (b) the hypothesis of collinear ISR.The first 3 sigma contours are shown.The normalization of the different solutions prevents the high mass contours from reaching the 1 sigma probability level, while the small difference in the low mass solutions originates from the jet charge information.The generated masses of the two W bosons in the event are marked with a cross.

Table 2 :
Results of a study of fragmentation effects on the M W measurement (see text).All results are given for the fitted mass analyses on the samples studied with respect to the standard Delphi JETSET sample, unless otherwise stated.The statistical error on the observed difference is also given.

Table 3 :
Results of a study of fragmentation effects on the Γ W ISR study (MeV/c 2 ) Model eν e qq ′ µν µ qq ′ τ ν τ qq ′ qq ′ qq ′ KORALW re-weighted -QEDPS +2 ± 7 0 ± 5 −12 ± 8 −16 ± 3 W measurement (see text).All results are given for the fitted mass analyses on the samples studied with respect to the standard Delphi JETSET sample, unless otherwise stated.The statistical error on the observed difference is also given.M

Table 5 :
A comparison of two ISR treatments for the Γ W measurement, the KORALW form and the QEDPS treatment.The statistical error on the result is also provided.

Table 6 :
Results of the Delphi studies on Bose-Einstein Correlations, see text for details.The BE 32 samples were produced with λ = 1.35, r = 0.6 fm

Table 7 :
Results of the Delphi studies on Colour reconnection effects, see text for details.For the SK1 model results are given for 30% of reconnected events.M W Systematic Errors (MeV/c 2 ) Sources of systematic error eν e qq ′ µν µ qq ′ τ ν τ qq ′ ℓν ℓ qq ′ qq ′ qq ′

Table 11 :
189GeV Γ W results.The error is divided into its statistical component, indicated stat, the main systematic component syst and the systematic from the beam energy uncertainty LEP .In the qq ′ qq ′ channel an error from final state interaction effects F SI is also included.The ℓν ℓ qq ′ results represents the combination of the results obtained in the three semi-leptonic channels.