Measurement of charged current triple gauge boson couplings using W pairs at LEP

Triple gauge boson couplings are measured from W-pair and single photon events recorded by the OPAL detector at LEP at center-of-mass energies between 183 - 209 GeV with a total integrated luminosity of 680 inverse picobarns. Only CP-conserving couplings are considered and SU(2)xU(1) relations between the WWZ and the WWgamma couolings are used, resulting in four independent couplings. Each coupling is determined in a separate fit, assuming the other couplings to take their Standard Model values. Fits are also done allowing some of the couplings to vary simultaneously. The results are compared with the Standard Model predictions.

k now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, U.S.A. l now at University of Toronto, Dept of Physics, Toronto, Canada m current address Bergische Universität, Wuppertal, Germany n now at University of Mining and Metallurgy, Cracow, Poland o now at University of California, San Diego, U.S.A. p now at Physics Dept Southern Methodist University, Dallas, TX 75275, U.S.A. q now at IPHE Université de Lausanne, CH-1015 Lausanne, Switzerland r now at IEKP Universität Karlsruhe, Germany s now at Universitaire Instelling Antwerpen, Physics Department, B-2610 Antwerpen, Belgium t now at RWTH Aachen, Germany u and High Energy Accelerator Research Organisation (KEK), Tsukuba, Ibaraki, Japan v now at University of Pennsylvania, Philadelphia, Pennsylvania, USA w now at TRIUMF, Vancouver, Canada x now at University of Florence, Florence, Italy.y now at IASF-CNR Sezione di Bologna, 40129 Bologna, Italy.* Deceased

Introduction
The triple gauge boson vertices, WWγ and WWZ, are expected in the Standard Model due to its non-Abelian nature.The Standard Model predicts exact values of the associated parameters, the triple gauge boson couplings (TGCs).Any significant deviation of the measured values from this prediction would be a signature of new physics beyond the Standard Model.Several measurements of these couplings have already been performed at LEP using W-pair [1][2][3][4][5], single W [6] and single photon [7] production.Limits on TGCs also exist from studies of di-boson production at the Tevatron [8].This paper reports results on TGCs using W-pair events from the full OPAL LEP2 data sample.
Assuming a general Lorentz-invariant Lagrangian which satisfies electromagnetic gauge invariance, charge conjugation (C) and parity (P) invariance, five independent TGC parameters are expected to describe the WWγ and WWZ vertices.These can be taken as g z 1 , κ z , κ γ , λ z and λ γ [9][10][11][12].In the Standard Model g z 1 =κ z =κ γ =1 and λ z =λ γ =0.When the restriction on C and P invariance is dropped, but their product, CP, is required to be invariant, a sixth coupling, g z 5 [9][10][11][12], which violates both C and P, is added.This coupling vanishes in the Standard Model.Couplings which violate CP [13] are not considered in this analysis.
Precision measurements on the Z 0 resonance and lower energy data motivate the following SU(2)×U (1) relations between the five C and P conserving couplings [9,11], Here ∆ indicates a deviation of the respective quantity from its Standard Model value and θ w is the weak mixing angle at tree level, defined by cos θ w = M W /M Z .These two relations leave three independent C and P conserving couplings, ∆κ γ , ∆g z 1 and λ(=λ γ =λ z ), plus the C and P violating coupling g z 5 .These couplings are not significantly restricted [14,15] by measurements from data collected at the Z 0 pole at LEP and SLC.
While the single W and single photon processes can also be used to measure TGCs, W-pair production is by far the most sensitive process for this measurement.Its sensitivity comes from the triple gauge boson vertex connecting the intermediate s-channel Z 0 /γ to the outgoing W-bosons.The TGC contribution depends on the helicity states of the outgoing W bosons, determining the angular distributions of the W bosons and of their decay products.The total W-pair cross-section is also affected by the TGCs, but for centre-of-mass energies well above the threshold for W-pair production and for small values of anomalous couplings its sensitivity, compared with the angular distributions, is much lower.Very often, however, the angular distribution analysis yields two ambiguous solutions for the TGCs, one being quite far away from the Standard Model values [16], and the information from the total cross-section can help in resolving this ambiguity.
The production and decay of W-pair events can be characterized by five angular variables.These are the W − production polar angle 1 , θ W , the polar and azimuthal angles, θ * and φ * , of the decay fermion from the W − in the W − rest frame 2 , and the corresponding polar and azimuthal angles of the decay anti-fermion from the W + .The experimental accessibility of these angles depends on the final state produced when the W bosons decay.
All W-pair final states are used in this study, namely the leptonic, ℓν ℓ ℓ ′ ν ℓ ′ , the semileptonic, qqℓν ℓ , and the hadronic, qqqq final states, with branching fractions of 10.5%, 43.9% and 45.6% respectively in the Standard Model.
The next section discusses the OPAL data and Monte Carlo samples.The following three sections then present the first part of this study using the W-pair angular distributions.For this purpose, the kinematic observables of each event are reconstructed as described in Section 3 and then used in Section 4 to construct optimal observables and analyse them in terms of the TGCs.The study of systematic errors for this part of the analysis is summarised in Section 5. Section 6 describes the second part of this analysis, where the total W-pair event rate is used to extract additional information on the TGCs.The TGC results are presented in Section 7 and summarised in the last section.

Data and Monte Carlo samples
The OPAL detector is described in detail elsewhere [17].The data were collected during 1997-2000 around eight different centre-of-mass energies.The integrated luminosity at each energy is evaluated using small angle Bhabha scattering events observed in the silicon tungsten forward calorimeter.The luminosity-weighted average values of the centre-of-mass energies and the corresponding luminosities are listed in Table 1.
The main Monte Carlo generator used to simulate the W-pair signal is KandY, which is a combination of KORALW1.51 [18] and YFSWW3 [19] running concurrently.These programs generate all four-fermion final states using the full set of electroweak diagrams including the W + W − production diagrams (class 3 CC03) and other four-fermion graphs, such as e + e − → Weν e , e + e − → Z 0 e + e − and e + e − → Z 0 Z 0 .All diagrams are corrected for initial and final state radiation, while a more complete correction of O(α) radiative effects is applied to the CC03 diagrams.For each centre-of-mass energy, 500,000 events have been generated using these Monte Carlo programs, and these events are used as our reference samples for the TGC analysis.The TGC dependence is obtained by reweighting these samples using the fourfermion matrix element calculation taken from the EXCALIBUR [20] Monte Carlo program.This program has also been used to generate samples with anomalous TGCs which are used to cross-check our TGC extraction method.Although EXCALIBUR includes less complete O(α) radiative effects, it can still be used for small anomalous TGC values, where the effect of the missing O(α) corrections on the anomalous TGC contribution can be neglected.
Background sources to the W-pair signal are four-fermion final states, such as qqℓℓ, which are 2 The axes of the right-handed coordinate system in the W rest-frame are defined such that the z-axis is along the parent W flight direction in the overall centre-of-mass system, and the y-axis is in the direction − → e − × − → W where − → e − is the electron beam direction and − → W is the parent W flight direction. 3In this paper, the three lowest order W-pair production diagrams, i.e. t-channel ν e exchange and s-channel Z 0 /γ exchange, are referred to as "CC03", following the notation of [9].
produced only by diagrams that do not involve the triple gauge boson vertex, two-fermion final states and two-photon processes.The four-fermion background is generated using KORALW1.42 [21] and grc4f [22].As an alternative for systematic studies we replace KORALW1.42 by EXCALIBUR.Two-fermion final states are generated using KK2F [23].For systematic studies, KORALZ [24], for e + e − →µ + µ − , e + e − →τ + τ − and e + e − →νν, and BHWIDE [25], for e + e − →e + e − , are also used.Background to the semileptonic final state from single-tag two-photon processes, where the outgoing electron or positron is detected, is evaluated using HERWIG [26], and F2GEN [27] is used as an alternative for systematic studies.Background from leptonic final states from untagged two-photon processes, where both outgoing electron and positron escape undetected into the beam pipe, is calculated using the Vermaseren [28] and BDK [29] generators.
To estimate the fragmentation and hadronisation systematics, Monte Carlo samples of W + W − events were produced by the KORALW1.42 generator, with the fragmentation and hadronisation stages generated separately by either PYTHIA5.7 [30], HERWIG6.2 [26] or ARIADNE4.08[31].These hadronic simulation generators have been tuned to Z 0 hadronic decays [32].Similarly, for the Z 0 /γ → qq background, the same hadronic simulation programs have been used to fragment and hadronise events produced by the KK2F generator.
All Monte Carlo samples mentioned above were processed by the full OPAL simulation program [33] and then reconstructed in the same way as the data.

W-pair event selection and reconstruction
W-pair events are selected with the same procedure used in the W-pair cross-section measurement as described in references [3,34].There are three independent selections corresponding to the three final states.Each candidate is then reconstructed in order to obtain the maximum possible information on the W production and decay angles, which are needed to extract the couplings.Events which cannot be well reconstructed are rejected from the sample, as will be described below.The numbers of candidates left at the different centre-of-mass energies for the three final states are listed, along with the expected values, in Table 1.
Three kinematic fits with different sets of requirements are used in the event reconstruction: by 3 due to the invisible neutrino.For qqτ ν τ events there is at least one additional unobserved neutrino from the τ decay, resulting in a loss of one constraint.The momentum sum of the track(s) assigned to the τ can still be used as an approximation of the τ flight direction, relying on its high boost, but the τ energy is left unknown.Finally for ℓν ℓ ℓ ′ ν ℓ ′ events, where none of the leptons is a τ , there are two invisible neutrinos.Hence, six constraints are lost and only requirement C is applied.
In the following we discuss the reconstruction of each final state separately.

Reconstruction of qqℓν ℓ final states
Starting from the sample used for the W-pair cross-section measurement, we select candidate qqeν e and qqµν µ events with a reconstructed lepton track.The track charge is needed to reconstruct the W − polar angle.For the qqτ ν τ events, either one track or a narrow jet consisting of three tracks is assigned as the τ decay product.
The OPAL tracking detectors are used to reconstruct the muon momentum and the electron direction, whereas the electromagnetic calorimeters are used to give a more accurate measurement of the electron energy.As explained above, the direction of τ candidates is given directly by the τ decay products, whilst the energy of the τ is estimated using a kinematic fit.
The remaining tracks and calorimeter clusters in the event are grouped into two jets using the Durham k ⊥ algorithm [36].The total energy and momentum of each of the jets are calculated using the charged track momenta and calorimeter energies and correcting for double counting employing the method described in [37].
The different kinematic fits described above are used to improve the resolution in the five kinematic variables used for the TGC analysis.This resolution is correlated with the kinematic fit probabilities which indicate the event reconstruction quality.Therefore we use the results of the kinematic fits to classify the events into nine quality classes which are used later in the TGC analysis.The exact definitions of these classes and the fraction of accepted events in each class To optimize the resolution in the kinematic variables used in the TGC analysis they are taken, according to the quality class, either from the directly measured values or from one of the three kinematic fits, as detailed in Table 2.The qqτ ν τ events which fail kinematic fit B or pass it with fit probability below 0.025 are rejected.This cut suppresses those qqτ ν τ events which are correctly identified as belonging to this decay channel but where the τ decay products are not identified correctly, leading to an incorrect estimate of the τ flight direction or its charge.The fraction of such events in the qqτ ν τ sample is thus reduced from 18% to 12%.
Finally, the kinematic variables are used to calculate the optimal observables (see Section 4) which are used to extract the TGCs.Those events (typically 1% of the sample) with optimal observable values in the far tails of their distributions are interpreted as being badly measured and discarded.
Performing the complete analysis on Monte Carlo samples generated with different anomalous couplings shows that our event selection does not introduce any bias on the TGC results.
The purities of the final samples are about 96%, 98% and 91% for the qqeν e , qqµν µ and qqτ ν τ final states respectively.The background consists mainly of TGC-independent four-fermion final states and Z 0 /γ → qq events.Cross migration between the qqτ ν τ and each of the qqeν e and qqµν µ decay channels is at the level of 4-5% and is considered as signal.
In the reconstruction of the qqℓν ℓ events we obtain cos θ W by summing the kinematically fitted four-momenta of the two jets.The decay angles of the leptonically decaying W are obtained from the charged lepton four-momentum, after boosting back to the parent W rest frame.For the hadronically decaying W it is not possible to distinguish between the jets of the quark and anti-quark.This ambiguity is taken into account in the analysis described below.
In Figure 1 we show the distributions of the five angles corresponding to the qqℓν ℓ data sample and the expected distributions for λ = ±0.5 and 0. These expected distributions are obtained from Standard Model Monte Carlo samples.To produce the λ = ±0.5 distributions the events have been appropriately reweighted using the EXCALIBUR matrix element calculation.All Monte Carlo distributions are normalised to the number of events observed in the data.The ℓ is shifted by 180 • in order to overlay W + and W − distributions in the same plot, which is valid so long as CP is conserved.
shaded histograms show the non-qqℓν ℓ background and their normalisation relative to the signal is according to the Standard Model.For the cos θ * jet and φ * jet distributions the jet with 0 ≤ φ * jet ≤ 180 • is arbitrarily chosen as the quark jet from the decay of the W − , or the anti-quark jet from the decay of the W + .Sensitivity to λ is apparent mainly for cos θ W .The contribution of cos θ * ℓ , φ * ℓ , cos θ * jet and φ * jet to the overall sensitivity enters mainly through their correlations with cos θ W .

Reconstruction of qqqq final states
Using the Durham k ⊥ algorithm [36] each selected qqqq candidate is reconstructed as four jets whose energies are corrected for double counting of charged track momenta and calorimeter energies [37].These four jets can be paired into W-bosons in three possible ways.To improve the resolution on the jet four-momenta we perform the five-constraint kinematic fit B for each possible jet pairing.In the following we consider only good jet pairings, namely those pairings with a successful fit B yielding a W-mass between 70 and 90 GeV.We accept only events with at least one good jet pairing.
After this cut the overall efficiency is around 80%, varying slightly with the centre-of-mass energy, and the contamination from other W + W − decays is below 0.3%.The major background contribution is from Z 0 /γ → qq where additional quarks or gluons are radiated off the primary quarks.This contribution is between 12% to 17%, decreasing with centre-of-mass energy.The contribution of TGC-independent four-fermion final states is 4.5%, except at 183 GeV where it drops to 3%.
We use a jet charge method as an estimator for the W charge.For the two jets coming from one W candidate we calculate the momentum weighted charge of all tracks contributing to this W candidate as, , where N W and N tot are the numbers of tracks in the W candidate and in the total event respectively, q i is the charge of the i th track and p ||,i is the momentum component parallel to the jet axis.We take the W candidate with the more positive charge as the W + .Around 38% of the events have more than one good jet pairing, but the fraction of events where all three possible jet pairings are good is below 5%.To select the correct jet pairing, the good pairings in the event are processed through a likelihood algorithm [38] which takes into account correlations between the input variables.We use as input to our likelihood the W candidate mass as obtained from the kinematic fit B, the absolute value of the difference between the W candidate masses as obtained from the kinematic fit A and the absolute value of the jet charge difference between the two W candidates.The algorithm has been tuned on Monte Carlo events at each centre-of-mass energy in order to give the optimal performance.The di-jet pairing with the highest likelihood output is taken as the correct pairing.The distribution of the highest jet pairing likelihood output is shown in Figure 2(a).The distribution of the W charge separation for all highest likelihood pairings is shown in Figure 2(b).
According to the W + W − KORALW Monte Carlo, the probability of selecting the correct di-jet combination is about 80%, and the probability of correct assignment of the W charge, once the correct pairing has been chosen, is about 77%.Both probabilities vary slightly with the centre-of-mass energy but depend strongly on the jet pairing likelihood output.
In order to reflect this dependence on the jet pairing likelihood output we use the output as a quality variable to classify the events into four quality classes which are later used in the TGC analysis.The definitions of these classes and the fraction of accepted events in each class are listed in Table 3.For classes with higher jet pairing likelihood output the probability for selecting the correct di-jet combination and the correct W charge assignment increases.The numbers are given in The distributions of the variables cos θ W , cos θ * jet and φ * jet are shown in Figure 3, along with the expected distributions for the Standard Model (λ=0) and for λ=±0.5.In the cos θ * jet and φ * jet distributions there are two entries per event, and according to our convention (see Section 3.1), the jet with 0 ≤ φ * jet ≤ 180 • is chosen as the quark jet from the decay of the W − , or the anti-quark jet from the decay of the W + .The Monte Carlo histograms are normalised to the number of data events.The shaded histograms represent the separate contributions from wrong pairings, correct pairing but wrong charge and background.Their relative normalisation is according to the Standard Model.The data are again described well by the Standard Model prediction.

Reconstruction of ℓν ℓ ℓ
′ ν ℓ ′ final states The reconstruction of ℓν ℓ ℓ ′ ν ℓ ′ final states is possible, using kinematic fit C, only if there is no τ -lepton in the final state.To reject events with τ -leptons we use the lepton identification algorithm designed for the W branching ratio measurement [34].The event rate analysis selection sometimes selects events with only one high momentum charged lepton reconstructed in the detector: such events were removed from the angular distribution analysis as the momenta of both charged leptons must be measured.As a result of these two cuts, the efficiency for ℓν ℓ ℓ ′ ν ℓ ′ events with ℓ, ℓ ′ = e, µ drops from ∼88% to ∼73%, varying slightly with centre-of-mass energy.
The contamination left in the sample from ℓν ℓ ℓ ′ ν ℓ ′ final states with a τ -lepton is 14% and the background from other final states is 2%.
Kinematic fit C has zero constraints and reduces to solving a quadratic equation.In the ideal case of no measurement errors and satisfying the conditions where fit C is valid, namely no ISR and narrow W width, one expects to obtain two real solutions corresponding to a two-fold ambiguity in the angle set cos θ W , φ * 1 and φ * 2 .There is no ambiguity in the angles θ * 1 and θ * 2 which have a one-to-one correspondence with the lepton momenta.In realistic conditions, however, including ISR, finite W width and measurement errors, one may obtain no real solution, but a pair of complex conjugate solutions.The fraction of events with complex solutions falls with increasing centre-of-mass energy from 30% to 24%, in agreement with Monte Carlo prediction.In these events, the imaginary parts of the complex solutions are set to zero, yielding a single real angle set.The total weight of all events is the same, independent of whether one or two solutions are employed.In the latter case, the two solutions each are assigned a weight of one half.
Figure 4 shows the distributions of the five reconstructed angles and the expected distributions for λ=0, ±1.The Monte Carlo histograms are normalized to the number of data events.The shaded histograms represent the contribution of background and their normalisation relative to the signal is according to the Standard Model.The data distributions agree with the expectation of the Standard Model.

Optimal observable analysis
The angular distributions from all final states are consistent with the Standard Model expectation and no evidence is seen for any significant contributions from anomalous couplings.For a quantitative study we use the method of optimal observables in the same way as in [4].This method relies on the linear dependence of the triple gauge vertex Lagrangian on the TGCs, corresponding to a second-order polynomial dependence of the differential cross-section, where S (0) , S (2) ij are functions of the five phase-space variables, Ω=(cos . For this kind of dependence it has been shown [39] that all information contained in the phase-space variables Ω, which is relevant to the four couplings, is retained in the whole set of 14 observables O 3,4).In this analysis we use the mean values of O In the one-parameter fit of α i we only use the corresponding first order, O i , and second order, O i , observables.In this way, the other couplings are not used in the fit.Similarly, in two-parameter fits where two couplings are allowed to vary simultaneously, we use only the five observables corresponding to these couplings and their mutual correlations.When three couplings are allowed to vary, nine observables are used.
The calculation of the mean optimal observables and their expected values is done separately for each centre-of-mass energy and for each of the qqℓν ℓ , qqqq and ℓν ℓ ℓ ′ ν ℓ ′ final states.The optimal observables are constructed for each event k with the set of phase-space variables Ω k using the analytic expression for the CC03 Born differential cross-section to calculate the values of ii (Ω k ) and S (2) ij (Ω k ).This calculation takes into account the reconstruction ambiguities for each particular final state.For the qqqq channel, the ambiguity in the W-charge determination is taken into account by weighting the S (0) , S analysis, where cos θ * 1 , φ * 1 (cos θ * 2 , φ * 2 ) are the decay angles of the charged lepton (anti-lepton) from the W − (W + ) in the parent W rest frame.All events enter with a total weight of one, but in the case of events with two ambiguous solutions, each solution enters with a weight of 0.5.The histograms, normalized to the number of data events, show the expectation of the Standard Model and the cases of λ = ±1.The shaded histograms represent the contribution of the background.
for each charge hypothesis with the corresponding probability.This probability is determined from the Monte Carlo as a function of the estimated charge difference |Q W 1 − Q W 2 | between the two W candidates in the event.
For qqℓν ℓ and qqqq events the mean optimal observables are calculated using weights assigned to the event according to its quality class (see Tables 2 and 3).For qqℓν ℓ events, the weight for a particular class is inversely proportional to the variance of the optimal observable distribution of the events in that class.For qqqq events, the weights are calculated in a similar way, using the variances of the optimal observable resolution distributions.In this way, events where the quality of the reconstruction is poor obtain a lower weight.Using this weighted mean has been found, as expected, to enhance our sensitivity.
The expected mean values of each observable as a function of the TGCs are calculated from four-fermion reference Monte Carlo samples generated according to the Standard Model.The EXCALIBUR matrix element calculation is used to reweight the Monte Carlo events to any TGC value required.The contribution of background events is also taken into account in the calculation of the expected mean values using corresponding Monte Carlo samples.
A χ 2 fit of the measured mean values O to the corresponding expectations E[O](α) is performed to extract the couplings from the data.The covariance matrix for the mean values is calculated from the Monte Carlo events and scaled to the number of data events.
In some of our fits the χ 2 function has a double minimum and therefore it is necessary to validate the errors obtained from the χ 2 fits.This is performed with a large number of Monte Carlo "experiments".In each of these "experiments", we use for each centre-of-mass energy and each final state a subsample of Standard Model Monte Carlo events.The size of each subsample corresponds to the luminosity for that final state at the particular centre-of-mass energy.Background events are also included in these subsamples in the appropriate proportions.For each "experiment", the various subsamples are analysed in the same way as the real data events.The distributions of the fit results are centred around the expected values, but there are some non-Gaussian tails.We therefore test the reliability of the error interval, as given by the region where χ 2 is no more than 1 above its minimum value, by counting the fraction of "experiments" where the correct value falls within this region.The error estimate is considered to be reliable if the calculated fraction is consistent with 68%; otherwise the corresponding elements of the covariance matrix for the O values are scaled up by the necessary factor to obtain 68% of the subsamples within the error interval.The resulting scaled covariance matrix is also used to analyse the real data.This is done separately for each fit.The resulting scale factors vary between 1 and 1.076 for the one-parameter fits.However, in the three-parameter fit of ∆κ γ , ∆g z 1 and λ, where double minima occur more frequently, a larger scale factor of 1.28 is obtained.These scale factors have been found to be appropriate also for the 95% C.L. intervals.a) There is a theoretical uncertainty in the O(α) correction due to missing higher orders in the YFSWW3 Monte Carlo generator program.As a conservative estimate of this uncertainty, this O(α) correction is removed, degrading the events to KORALW1.42 level.The efffect of this change is taken as a systematic uncertainty.
b) KORALW includes ISR up to O(α 3 ).As a conservative estimate of missing higher orders, we replace it with ISR up to O(α) only and take the difference as a systematic uncertainty.c) Another uncertainty in the signal Monte Carlo generator is related to the PYTHIA fragmentation and hadronisation model which is used.To study the effect of this uncertainty we compare results of the analysis when applied to W + W − Monte Carlo samples generated with KORALW1.42 using either PYTHIA, HERWIG or ARIADNE for the fragmentation and hadronisation phase.The effect of changing to HERWIG is the larger and is used to estimate the uncertainty.The main effect is on the qqqq final state, where HERWIG predicts a higher probability of correct jet pairing and correct W charge assignment than PYTHIA.This effect is partially due to the average charged multiplicity for light quarks being lower in our tuned version of HERWIG than in PYTHIA.Since the data are in agreement with PYTHIA [41], we weight the HERWIG Monte Carlo events so as to reproduce the same charged multiplicity distribution as PYTHIA.The effect of this weighting is a reduction of the systematic error by about 15-30%.
d) To assess the uncertainty in the TGC matrix element calculation which is used to weight our Monte Carlo events, we replace the calculation from the EXCALIBUR Monte Carlo generator program with the one from grc4f.
e) The Monte Carlo simulation of the OPAL detector has been studied using back-to-back jets and leptons in Z 0 events collected during calibration runs which were taken regularly every year between the high energy runs.Some scaling and smearing had to be applied to the reconstructed jets and leptons in the Monte Carlo in order to achieve the best agreement with data.The uncertainties in these corrections correspond to variations in our results which are taken as systematic errors.These corrections are usually applied after the event selection.To estimate their effect on the selection of our sample, they were applied as a systematic check before the event selection, and the difference from the standard result was added in quadrature with the effects of the uncertainties on the corrections.
f) Uncertainties in the background estimation are determined by varying both its shape and normalisation.To study the effect of the shape we use alternative samples for each background source, keeping the normalisation fixed.For the TGC-independent fourfermion background we use EXCALIBUR and grc4f instead of KORALW.For the Z 0 /γ → qq background we use KK2F with HERWIG or ARIADNE rather than PYTHIA fragmentation.Finally, the HERWIG two-photon samples are replaced with F2GEN.For the background normalisation uncertainty the contribution from each source is scaled, one at a time, by the factors 1.2, 1.3 and 2.0 respectively.These scale factors are determined from studies with control data samples and from the differences between the alternative Monte Carlo generators, as explained in [34].
g) The consequences of the uncertainties in the centre-of-mass energy and the W mass are estimated as follows.For the centre-of-mass energy uncertainty of 40 MeV, Monte Carlo samples at different centre-of-mass energies are used as pseudo-data.For the W mass we assume M W =80.456±0.059GeV/c 2 [35] which is 0.126 GeV/c 2 higher than the one used in our Monte Carlo generators.We do not correct our TGC fit result for this shift in M W , but it is accounted for in our systematic error, along with the M W measurement error.The combined effect is assessed using Monte Carlo samples generated with different W masses.
h) Bose-Einstein correlations (BEC) in qqqq events might affect the reconstruction of the W bosons and measured W charge distribution.We investigate this effect with a Monte Carlo program simulating BEC via re-adjustment of final state momenta using the LUBOEI model [42].When both inter-W and intra-W correlations are taken into account, the bias is found to be larger than the case where only intra-W boson correlations are present.
i) The jet reconstruction and the measured W charge distribution in qqqq events might also be affected by colour reconnection.This effect is investigated with several MC samples generated according to the Sjöstrand-Khoze [43] models.We consider model I with colour reconnection probabilities of 34% and 97%.Models II and II' are also considered.Model I with colour reconnection probability of 97% produces the largest bias.
The Monte Carlo samples are larger than the data by a factor of 150 -900, depending on the centre-of-mass energy.Therefore, the uncertainty due to Monte Carlo statistics is neglected.
The effect of each of these sources on the mean values of the optimal observables is estimated separately for each centre-of-mass energy and for each final state.The deviations obtained in the mean observables are used to construct a systematic error covariance matrix, where the systematic deviations from each source are assumed to be fully correlated between all centreof-mass energies.For the combination of all final states, the covariance matrix is extended to incorporate the three final states.All uncertainties, except for those associated with the background, are assumed to be fully correlated between the relevant final states.This matrix is added to the statistical covariance matrix, and the fit is redone to obtain the TGC parameters with their total (statistical and systematic) errors.For demonstration, we also calculate the separate contribution from each source of systematic error to the error on the TGC parameters, after combining the different centre-of-mass energies and final states.These are listed in Table 4.

Event rate TGC analysis
The event rate study, unlike the angular distribution analysis, requires a detailed investigation of the systematic uncertainties associated with the overall selection efficiency.This investigation is part of the W-pair cross-section analysis and has to be done separately for each centre-ofmass energy and for each final state.So far, it has been completed for centre-of-mass energies of 183 GeV [3] and 189 GeV [34].The TGC event rate analysis is therefore also restricted to these two centre-of-mass energies.In the case where the angular distribution analysis has two TGC solutions, one being far away from the Standard Model value, the results from this restricted event rate analysis would still be sufficient to reject that solution.The data from higher centre-of-mass energies would make a negligible additional contribution to our results.Compared with the angular distributions, the sensitivity of the event rates to TGCs around the Standard Model values is much smaller and its dependence on the overall luminosity is weaker.
The event rate analysis has been described already in our previous publications [3,4].The numbers of expected events are re-evaluated for this analysis using the KandY Monte Carlo samples to obtain the Standard Model expected number of events.
The same selection, efficiency and background estimates as in the total cross-section analysis are used here.The numbers of selected events in each final state and for each of the two centre-ofmass energies are quoted in Table 5 along with the Monte Carlo expectations.The systematic errors on the expected values are due to Monte Carlo statistics, 0.5% uncertainty in the total cross-section, 0.2% uncertainty in the luminosity, uncertainties due to data/MC differences, tracking losses, detector occupancy and fragmentation.A detailed description of all these sources can be found in [34].The expected numbers of signal events include contributions from all four-fermion final states which can be produced by diagrams involving the WWγ and WWZ triple gauge boson vertices.The other four-fermion final states are considered as background along with the two-fermion final states and two-photon processes.The total numbers of expected events are consistent within the statistical and systematic errors with the observed numbers.Therefore, there is no evidence for any significant contribution from anomalous couplings.A quantitative study of TGCs from these results is performed by comparing the numbers of observed events in each of the three event selection channels at each of the two centre-of-mass energies with the expected numbers of signal and background events, where the number of signal events is parametrised as a second-order polynomial in the TGCs.This parametrisation is based on the second-order polynomial dependence of the cross-section, as was used also in the optimal observable analysis (see Section 4).The polynomial coefficients in terms of event rate for each final state are calculated by weighting our Monte Carlo events to correspond to various TGC values using a four-fermion matrix element calculation procedure taken from the EXCALIBUR Monte Carlo generator program.

√ s
The probability to observe the measured number of candidates, given the expected value, is calculated using a Gaussian distribution where the errors also include the systematic uncertainties.We assume full correlation between the systematic uncertainties for the three final states except for efficiency and background which are taken to be uncorrelated.All systematic uncertainties are assumed to be fully correlated between the two different centre-of-mass energies.

Combined TGC results
The results of the various fits are listed in Table 6.We present first the results of the angular distribution analysis without systematic errors obtained by performing one-parameter fits for each final state simultaneously for the eight centre-of-mass energies.For comparison, the average fit errors of the Monte Carlo "experiments" are also quoted in Table 6, demonstrating the different sensitivities for the various final states and TGC parameters.In most cases, the fit errors are rather close to the expected ones, except for ∆κ γ , where the fit errors for some final states are also very asymmetric due to the proximity of a second minimum in the χ 2 function.
Next we present the results of the angular distribution analysis with the systematic errors included as explained in Section 5.This modifies the relative contributions to the fit of the different optimal observables and different centre-of-mass energies, as well as the correlations between them.Therefore, not only the fit errors are affected, but also the fit results themselves.The corresponding log L curves5 for the different channels are shown in Figure 5.The results from the three final states are combined, taking into account the common systematic errors as explained in Section 5, and the combined results are listed in Table 6.
The results of the event rate analysis are presented by log L curves shown by the thin solid lines in Figure 5.They are combined with the results of the angular distribution analysis and the combined log L plots are shown by the thick solid lines.The correlation between the results of the two analyses, which is due to uncertainties in the background normalisation, centre-of-mass energy and W mass, is below 7%.Therefore this correlation is neglected in the combination.The numerical results of the combination are listed in Table 6 along with the corresponding χ 2 /d.o.f.As expected, the event rate information has relatively little impact on the final results.

Figure 1 :
Figure 1: Distributions of the kinematic variables cos θ W , cos θ * ℓ , φ * ℓ , cos θ * jet and φ * jet , obtained from the qqℓν ℓ candidates.The solid points represent the data.The histograms show the Monte Carlo expectation of the Standard Model (solid line) and the cases of λ = +0.5 and λ = −0.5 (dashed and dotted lines respectively).The Monte Carlo histograms are normalised to the number of data events.The shaded histograms show the non-qqℓν ℓ background.In the case of W + → lν decays the value of φ *ℓ is shifted by 180 • in order to overlay W + and W − distributions in the same plot, which is valid so long as CP is conserved.

Figure 2 :
Figure 2: Distributions for qqqq data events (points) for the KandY four-fermion Monte Carlo (histogram) of (a) jet-pairing likelihood corresponding to the most likely combination; (b) charge separation between the two W candidates.

Figure 3 :
Figure 3: Distribution of the kinematic variables cos θ W , cos θ * jet , and φ * jet for qqqq events.In the cos θ * jet , and φ * jet distributions there are two entries per event.The solid points represent the data.The histograms show the Monte Carlo expectation of the Standard Model (solid line) and the cases of λ=+0.5 and λ=-0.5 (dashed and dotted lines respectively).The Monte Carlo histograms are normalised to the number of data events.The open histograms show the contribution from correct pairing and correct charge assignment.The single hatched histograms represent the contribution from correct pairing, but wrong charge assignment.The crosshatched histograms correspond to the contribution from wrong pairing, and the dark histograms show the contribution of background.

Figure 5 :Figure 6 :
Figure 5: Negative log-likelihood curves obtained from the angular distribution information of W-pair events taken at the centre-of-mass energies of 183 -209 GeV, separately for the three different final states: qqℓν ℓ (dotted lines), qqqq (dashed lines) and ℓν ℓ ℓ ′ ν ℓ ′ (dash-dotted lines).The thin solid lines describe the contribution of the event rate information from all final states corresponding to centre-of-mass energies of 183 and 189 GeV.The systematic errors are included.The thick solid lines are obtained by combining all sources of information.The curves for each TGC parameter are obtained by setting the other three parameters to their Standard Model values.

Table 1 :
Observed and expected numbers of data candidates selected for the angular distribution analysis after all cuts in the different final states and for different centre-of-mass energies.

Table 2 :
Definition of the nine quality classes of qqℓν ℓ events.P A , P B and P C are the probabilities of kinematic fits A, B and C.
are listed in Table2.The event population in each class is well modelled by the Monte Carlo.

Table 3 .
The population of each class is well modelled by Monte Carlo.

Table 3 :
Definition of the four quality classes of qqqq events.The correct charge assignment probability is calculated only for those events with correct jet pairing.

Table 4 :
Contributions to the systematic uncertainties in the determination of the TGC parameters.The total systematic errors are obtained by combining all contributions at the optimal observable level and then calculating the effect on the TGC parameters.

Table 5 :
Observed and expected numbers of events in each W + W − final state for the centreof-mass energies of 182.7 GeV and 188.6 GeV.The separation between signal and background is explained in the text.All the quoted errors are systematic.