Measurement of W boson polarisations and CP-violating triple gauge couplings from W+W- production at LEP

Measurements are presented of the polarisation of W+W- boson pairs produced in e+e- collisions, and of CP-violating WWZ and WWGamma trilinear gauge couplings. The data were recorded by the OPAL experiment at LEP during 1998, where a total integrated luminosity of 183 pb^-1 was obtained at a centre-of-mass energy of 189 GeV. The measurements are performed through a spin density matrix analysis of the W boson decay products. The fraction of W bosons produced with longitudinal polarisation was found to be sigma_L/sigma_total = (21.0 +- 3.3 +- 1.6)% where the first error is statistical and the second systematic. The joint W boson pair production fractions were found to be sigma_TT/sigma_total = (78.1 +- 9.0 +- 3.2) %, sigma_LL/sigma_total = (20.1 +- 7.2 +- 1.8) % and sigma_TL/sigma_total = (1.8 +- 14.7 +- 3.8) %. In the CP-violating trilinear gauge coupling sector we find kappa_z = -0.20 +0.10 -0.07, g^z_4 = -0.02 +0.32 -0.33 and lambda_z = -0.18 +0.24 -0.16, where errors include both statistical and systematic uncertainties. In each case the coupling is determined with all other couplings set to their Standard Model values except those related to the measured coupling via SU(2)_LxU(1)_Y symmetry. These results are consistent with Standard Model expectations.


Introduction
We report measurements of the properties of W pair production in e + e − collisions using data recorded by the OPAL detector at LEP at a centre-of-mass energy of 189 GeV with a total integrated luminosity of 183 pb −1 . We perform a spin density matrix (SDM) analysis [1,2] of the production and decay properties of the W bosons using the semi-leptonic WW → qqeν, WW → qqµν and WW → qqτ ν final states. Using suitable summations of SDM elements we present measurements of the inclusive production cross-sections for each of the transverse and longitudinal polarisation states of the W. The SDM elements are also sensitive to triple gauge couplings. We present measurements of CP-violating couplings involving the W ± , Z 0 and photon.
The doubly resonant e + e − → W + W − production process proceeds via s-channel Z 0 or photon exchange, or via t-channel neutrino exchange, collectively known as the CC03 diagrams. The s-channel processes contain the WWZ and WWγ triple gauge boson vertices. The most general Lorentz invariant Lagrangian describing these vertices [3,4] contains 14 independent couplings. In order to facilitate measurements the number of parameters is often reduced to three by assuming SU (2) L × U (1) Y gauge invariance and charge conjugation (C) and parity (P) invariance. The resulting independent couplings are conventionally taken as: ∆κ γ , ∆g z 1 and λ [5]. Measurements of these couplings using data recorded at LEP2 [6][7][8][9][10] and the Tevatron [11] have been reported elsewhere.
If C, P and CP-invariance are not assumed then several additional couplings may be present. The CP-violating ones can be taken asκ V andλ V which violate P and conserve C, and g V 4 which violates C but conserves P ( V = Z 0 or γ) [3]. A further parameter, g V 5 , violates both C and P but conserves CP so is not considered here. SU (2) L × U (1) Y symmetry imposes the following relations [12][13][14] which are assumed in this analysis:κ z = − tan 2 θ wκγ λ z =λ γ (1) All of the triple gauge boson couplings (TGCs) can in principle be measured through observations of the W production angular distribution, and distributions of the W decay products. One way of realising such an analysis is through the spin density matrix (SDM) approach. In this approach the individual contributions to the W production angular distribution arising from each of the different possible helicity states of the W bosons can in principle be determined exclusively. These exclusive SDM distributions exhibit different behaviour with respect to each of the TGCs. The TGCs can therefore be determined from the SDM elements in a second step. In W boson pair production, the SDM analysis is particularly suited to the extraction of the CP-violating couplings. Indeed, the imaginary parts of the off-diagonal elements of the SDM are completely insensitive to the CPconserving couplings, and will only deviate from their Standard Model predictions in the presence of CP-violation at the triple gauge boson vertex.
Before proceeding to the TGC measurements we investigate the exclusive production cross-sections for each of the transverse and longitudinal polarisations of W bosons. These can be made in the context of the SDM analysis by suitable summations of SDM elements which are described in detail in the following section. The study of the longitudinal cross-section is particularly interesting as this degree of freedom of the W only arises in the Standard Model through the electro-weak symmetry breaking mechanism. Previous measurements of the proportion of longitudinally polarised W bosons produced at LEP2 have been been reported by the OPAL [9] and L3 [15] experiments.
In addition we present limits on the CP-violating parameters. Previous limits on the CP-violating TGCs have been obtained at the DELPHI experiment using data at 161 GeV and 172 GeV [8]. Limits on the CP-violating WWγ couplings,κ γ andλ γ , have been obtained at the D0 experiment [16] from the process pp → ℓνγ + X. Measurements of the neutron electric dipole moment show that the electromagnetic interaction is CP-conserving to very high accuracy [17].

The Spin Density Matrix
The two-particle joint spin density matrix (SDM) [4] describes the contribution of each of the helicity states of the W bosons to the W + W − production cross-section.
The amplitude for the production of a W − with helicity τ − (= −1,0 or 1) and a W + with helicity τ + is denoted as F (λ) [1,4], where λ (= ± 1 2 ) denotes the helicity of the e − , s is the square of the centre-of-mass energy and cos θ W is the W production angle in the centre-of-mass frame. The SDM (ρ) elements are normalised products of these amplitudes given by [4]: The normalisation ensures that τ − τ + ρ τ − τ − τ + τ + = 1. The ρ matrix is Hermitian, giving 80 independent real coefficients which may be experimentally measured. The diagonal elements are defined as the subset of elements where τ − = τ ′ − and τ + = τ ′ + . These diagonal elements are purely real and are equivalent to the probability of producing a final W + W − state with helicities τ + and τ − respectively. The off-diagonal elements represent interference terms between the different helicity states of each W.
In the context of a limited sample of events it may not be possible to measure all of the components independently. It is then useful to consider the single particle density matrix for either the W − or W + , formed by summation over the helicity states of the other W. For example the W − matrix elements ρ W − τ − τ ′ − are obtained by summation over τ + = τ ′ + : The Hermitian matrix ρ W − has eight independent real coefficients, and satisfies the normalisation con- The three real diagonal elements represent the relative production probabilities for a final state W − with a particular helicity.
The constraints of CPT and CP-invariance impose additional symmetries on the density matrix at tree level [2]. CPT-invariance imposes the conditions: CP invariance imposes the condition: Thus CPT and CP-conservation together dictate that all coefficients from ρ W τ 1 τ 2 are real. Deviations from the validity of equation (6) would represent an unambiguous signal for CP-violation at tree level.
Deviations from the equality of equation (5) represent a signal of higher-order (loop) effects beyond tree level [2]. The exclusive differential cross-sections for the production of a W with transverse (T) or longitudinal (L) helicity are obtained by weighting the total differential cross-section by the relevant elements of the single particle ρ matrix: and the corresponding total cross-sections are given by: The single W SDM elements can be obtained from measurements of the properties of the W decay products by the application of suitable projection operators Λ τ τ ′ [1,2] which assume the V−A coupling of the W to fermions. The W decays are characterised by the polar and azimuthal angles of the decay fermion in the W rest frame, θ * and φ * respectively. In this analysis we consider only the semi-leptonic event class, WW → qqℓν, where one W decays to a lepton (e, µ or τ ) and a neutrino and the other to two hadronic jets. In this case the values of θ * and φ * of the lepton may be determined unambiguously and the corresponding single W SDM elements are given by: In section 4 we describe how this expression is realised as a sum over the events observed in the data sample.
In the case of the hadronically decaying W differentiation between the particle and anti-particle decay products is very difficult. However, certain projection operators are symmetric under the transformation cos θ * → − cos θ * , φ * → φ * +π, i.e. the interchange of one of the decay particles from the W boson with the other decay particle, for example the u-type quark with the d-type in a hadronically decaying W boson. This means that a number of combinations of the SDM elements may still be extracted from W → qq decays: these include ρ ++ + ρ −− and ρ 00 . Both W bosons in the event can therefore be used to measure the polarised cross-section, since only these terms appear in equations (7).
Returning to the full two particle SDM, the analogous expressions to (7) but describing the differential cross-section for both W bosons in the pair being transversely polarised (TT), both being longitudinally polarised (LL), or one of each polarisation (TL) are: dσ TL dcos θ W = (ρ ++00 + ρ −−00 + ρ 00++ + ρ 00−− ) dσ dcos θ W and the analogous expression to (9) is where the integral is now over the decay angles of both W bosons.

The Data Sample
The W-pair data used in this analysis were collected by the OPAL [18] detector at LEP. The accepted integrated luminosity in 1998, evaluated using small angle Bhabha scattering events observed in the silicon tungsten forward calorimeter [19], is 183.14 ± 0.55 pb −1 [20]. The luminosity-weighted mean centre-of-mass energy for the data sample is √ s = 188.64 ± 0.04 GeV.
In this analysis we use only the W-pair events decaying to the qqeν, qqµν and qqτ ν channels. These qqℓν events were first selected using the likelihood selection described in [9,20,21]. The selection is designed to optimise the rejection of Z 0 /γ → qq and four-fermion backgrounds. A total of 1252 qqℓν candidates was selected at this stage. Monte Carlo studies show this selection is about 88% efficient.
Further cuts are applied in order to obtain a sample of well reconstructed events. A full description can be found in [10]. A brief overview of the procedure is given here. For the qqeν and qqµν events a well reconstructed lepton track is required. For the qqτ ν events, either one track or a narrow jet consisting of three tracks is assigned as the τ decay product.
Kinematic fits are now applied to the events. For the qqeν and qqµν events a one-constraint fit is applied that requires energy-momentum conservation and neglects any initial state radiation. Events are accepted if the fit converges with a probability larger than 0.001. An improvement is made to the resolution of the angular observables by performing a second kinematic fit which constrains each reconstructed W mass to the average value measured at the Tevatron [22]. Events for which the fit converges with a probability of at least 0.001 are accepted; for other events we revert to the previous fit.
For the qqτ ν events a different kinematic fit is applied [9]; here the reconstructed mass of the two W bosons is required to be equal. In order to obtain a W mass from the leptonic part of these events, it is necessary to assume that the direction of the visible part of the τ decay approximates the direction of the τ lepton. The fit is required to converge with a probability of at least 0.025. This cut rejects 41% of the background. It also rejects about 14% of the signal events, but preferentially rejects events where the τ decay products are not correctly identified.
After application of all selection cuts 1065 qqℓν candidates remain, with 359 in the qqeν channel, 386 in the qqµν channel and 320 in the qqτ ν channel. The W production and decay angles for these events are shown in figure 1. The background contributions to the data sample are estimated using Monte Carlo samples. The resulting contaminations are expected to be four-fermion (non-CC03) 3.1%, Z 0 /γ → qq 2.1%, two-photon (multiperipheral) 0.2% and CC03 four-fermion (WW → qqqq and ℓνℓν) 0.5%.

Monte Carlo
A number of Monte Carlo models are used to provide estimates of efficiencies and backgrounds as well as the expected W-pair production and decay angular distributions for different TGC values. The Monte Carlo samples are generated at √ s = 189 GeV. All Monte Carlo samples mentioned below are passed through the full OPAL simulation program [23] and then subjected to the same selection and reconstruction procedure as the data. The ERATO [13] Monte Carlo program is used in this analysis because it can generate samples of four-fermion Monte Carlo events with non-Standard Model values of CP-violating anomalous couplings. The EXCALIBUR [24] Monte Carlo is also used extensively in this analysis for correction of detector effects and as the main tool to compare to the 189 GeV data. A comparison between ERATO and EXCALIBUR angular distributions and spin density matrix elements was undertaken and no statistically significant differences were seen between them. The other four-fermion generator used to calculate background contributions and systematic uncertainties is grc4f [25]. The WW generators used to calculate theoretical predictions for the W-pair polarisations and to calculate systematic uncertainties are EXCALIBUR, HERWIG [26], PYTHIA [27] and KORALW [28]. The background Z 0 /γ → qq samples are generated using PYTHIA. For samples of two-photon (multiperipheral) background PYTHIA, PHOJET [29] and HERWIG were used.

Experimental Method
Measurements of the production cross-sections for each of the transverse and longitudinal states of the W bosons are obtained from summations of the SDM elements obtained from the qqℓν data sample, after correcting for detector acceptance, resolution effects and background contamination.
To calculate the SDM elements the events are divided into eight equal bins, k, of cos θ W . In each bin the SDM elements are obtained by a summation over events, i, weighted by the relevant SDM projection operators as shown in equation (12). This is effectively the realisation of equation (9) as a summation over observed events. N k is the number of events in bin k. In the extraction of all the individual W SDM elements CPT invariance is assumed, so all leptonic (hadronic) W + and W − decays can be combined.
Estimates of the production cross-sections are derived from the diagonal SDM elements; therefore the projection operators only involve the polar decay angle cos θ * . It is not required that the SDM elements be positive definite. Both the leptonic and hadronic W decays can be used since the SDM combinations occurring in equation (7) are symmetric with respect to the use of the fermion or antifermion polar decay angle. Each event is also weighted by a correction factor, f , for detector acceptance and resolution effects: (14) where N cor k is the corrected number of events in bin k: The factors f k (θ * i ) are obtained from fully simulated Monte Carlo events, and are calculated as a function of cos θ W and cos θ * . They are defined to be the ratio of the number of reconstructed events to the number of generated events in each bin. The correction factors have an effect of between 5% and 20%. The typical resolution of the measured cos θ W is found to be 0.05, which is less than 20% of the bin width used in correcting this distribution. The typical resolution of the measured cos θ * is found to be 0.07, which is 70% of the bin width used in the calculation of the correction factors.
The distribution of uncorrected events, shown in figure 1, also has to be corrected for detector acceptance and resolution effects. This is done in a similar manner to the correction of the SDM elements by calculating a correction factor from Monte Carlo events in a bin-wise fashion and applying it to the data. The estimated background contribution is subtracted from the SDM elements and the dσ/dcos θ W distribution.
The polarised differential cross-sections, dσ L /dcos θ W and dσ T /dcos θ W , are then obtained by multiplying the binned unpolarised differential cross-section by the relevant SDM combinations, following equations (7). Figure 2 shows the individual W boson differential cross-sections obtained from the data by this method. The transverse and longitudinal components obtained from the leptonically decaying W and the hadronically decaying W are shown separately, together with the values obtained by combining the two. The polarisations of the W bosons in the W-pair event are correlated. The correlation is estimated to be 0.07, and is taken into account in the errors on the combined cross-sections.

Experimental Results
The fraction of each polarisation state, obtained by integrating over cos θ W and then dividing by the total cross-section, is given in table 1. Once again the correlated polarisations of the two W bosons are taken into account when deriving the errors. The estimation of the systematic errors shown in tables 1 and 2 is described in section 6.
The joint polarised cross-sections are obtained by extracting the joint SDM elements in a similar way to equations (12) where in this case the ρ and Λ are those corresponding to the particular SDM combinations appearing in equation (10). Note that there is now a projection operator for each W, and the correction factor is binned in terms of both polar decay angles. Figure 3 shows the joint differential cross-sections obtained by this method, these cross-sections are not constrained to be positive definite. between each helicity fraction, calculated from the data, are found to be: σ TT :σ LL = 0.67±0.02, σ TT :σ TL = −0.93±0.01 and σ LL :σ TL = −0.89±0.01.
The individual W polarised differential cross-sections show good agreement with the Standard Model predictions, as do the overall fraction of each individual W polarisation, as seen in table 1. The W pair polarised cross-sections show less good agreement, but σ LL /σ total is within two standard deviations of the Standard Model prediction and the other two are just over two standard deviations away. The χ 2 value for the three measurements compared to the Standard Model expectations, including systematic uncertainties, is 4.7. This corresponds to a χ 2 probability of 10%.

Measured
Expected σ TT /σ total 0.781 ± 0.090 ± 0.033 0.572 ± 0.010 σ LL /σ total 0.201 ± 0.072 ± 0.018 0.086 ± 0.008 σ TL /σ total 0.018 ± 0.147 ± 0.038 0.342 ± 0.016 Table 2: The fraction of each helicity combination of WW pairs. The expected values are calculated from Monte Carlo studies. The first errors on the measured fractions are statistical and the second are systematic. The errors on the expected fractions are statistical only.

Measurement of Anomalous Couplings and Test of CP-Invariance
Measurements of anomalous coupling parameters are obtained through a comparison of distributions of SDM elements obtained from the leptonically decaying W bosons in the qqℓν data sample with distributions obtained from fully simulated Monte Carlo samples. In contrast to the polarised crosssection measurements, neither data nor Monte Carlo events are corrected for detector or acceptance effects. Instead the experimentally observed SDM distributions are compared directly with those for fully detector simulated Monte Carlo events. Monte Carlo samples for arbitrary TGC values are obtained by a re-weighting technique applied to a large Standard Model sample. A χ 2 minimisation procedure is then used to find the simulated data set which best fits the data and hence obtain the best fit TGC parameter value.
The events are again divided into eight equal bins of cos θ W . Equation (12)  This method is used to calculate both the CP-conserving and the CP-violating couplings. However, when performing the fit for CP-conserving couplings only the six real SDM coefficients are used, and in the case of CP violating couplings all nine SDM coefficients are included. The effect of correlations between SDM elements is included. Although there is no correlation between cos θ W bins, within each bin the measurements of the SDM elements are highly correlated. The correlations are obtained directly from the data.
The Monte Carlo sample used in the fitting procedure is generated for Standard Model couplings using EXCALIBUR. All TGC-dependent four-fermion diagrams are included in this sample. SDM distributions for arbitrary couplings are obtained by re-weighting the fully simulated EXCALIBUR Monte Carlo events. Technically this was achieved by re-weighting events with the matrix elements from [1]. The five characteristic W production and decay angles are constructed from the original four-momenta of the primary fermions in the simulated EXCALIBUR events. This treatment neglects the effects of four-fermion background, but these are checked explicitly for CP-conserving couplings, and found to be negligible. The expected contribution from four-fermion background events produced in channels unaffected by triple gauge couplings, and those from other backgrounds, are added to the EXCALIBUR sample.
As the SDM elements are normalised to the number of events in each cos θ W bin, information from the production angle of the W is effectively removed from the SDM fit. In order to include this an independent χ 2 is derived from the comparison of the shape of the cos θ W distributions and this is incorporated into the fit. No information from the overall cross-section is included.
The method of re-weighting and fitting is tested by performing fits to large samples of fully simu- These latter results are consistent with those measured in [10]. All couplings are set to their Standard Model values except the coupling being measured and those related to it via the SU (2) L × U (1) Y symmetry. A full breakdown of the results for the CP-violating couplings is shown in table 3.
The χ 2 plots for all the fits including systematics can be seen in figure 5. The double minima in the χ 2 curves for the CP-violating couplings derived from the cos θ W distribution alone reflect that the cos θ W distribution is sensitive only to the absolute magnitude of the coupling. The same is true for the real SDM elements resulting in the double minima in the SDM χ 2 curve forκ Z . It is only the imaginary parts that lift the degeneracy. It is evident that the CP-violating couplings are better constrained by the SDM elements than by the W boson production angular distribution, whereas for the CP-conserving couplings the converse is true. A further test of tree level CP violation is given by comparing the imaginary SDM elements from the W + to those from the W − , as in equation (6). These comparisons, which give a model independent test of CP-violation in the triple gauge coupling, are shown in figure 6. Deviations from zero would only be due to CP-violation at tree level. No deviations are seen, complementing the results of the measurements of the CP-violating TGCs.
Also shown in figure 6 are the combinations of imaginary SDM elements that test for effects beyond tree level, as described by equation (5). Any deviation from zero in these plots could only be caused by either loop effects or CPT-violation; all CP-violating tree level effects cancel out. The data are consistent with no effect.

Systematic Uncertainties
Systematic uncertainties in the measurements are calculated as described below. The individual contributions to the errors on the polarised cross-section fractions are given in table 4. For the measurements of the TGCs, all error sources listed below are included except that from the detector correction. For all TGC measurements the systematic errors are included as extra uncertainties on the contents of each bin of the SDM element and cos θ W distributions. Each systematic uncertainty is taken to be uncorrelated between bins in both cos θ W and the different SDM element distributions.
• Jet reconstruction: Uncertainties in the modelling of jet reconstruction are estimated by varying the reconstruction in a Monte Carlo sample. The resolutions of the three jet parameters (energy, cos θ j , φ j ) are varied by 10%, and the energy is shifted by 5% to account for systematic uncertainties. This is done for both the quark jets and the τ jets. The size of the variations for the quark jets is determined from extensive studies of back-to-back jets at Z 0 energies. A possible systematic shift in the reconstructed direction of the boson has been estimated using radiative Z 0 /γ → qq events. The possible shift in | cos θ W | was found to be less than 0.01 [10]. The same uncertainties are taken on the τ jets as for the quark jets.
• Hadronisation: Uncertainties due to the hadronisation model are estimated by comparing Monte Carlo fragmented with HERWIG5.9 to Monte Carlo using the JETSET7.4 [27] hadronisation scheme. Both Monte Carlos have been tuned to OPAL data. Variations in the calculated polarised crosssections, SDM elements and cos θ W distribution between the different samples are taken as the systematic uncertainty.
• Monte Carlo generators: The modelling of the four-fermion production processes is assessed by comparing SDM elements and cos θ W distributions from the ERATO and grc4f Monte Carlo programs to those from EXCALIBUR. All these generators have a different calculation of the matrix elements and a different treatment of ISR. The calculated values of the individual W and W-pair polarised cross-section fractions are also compared. Variations in the calculated polarised cross-sections, SDM elements and cos θ W distribution between the different samples are taken as the systematic uncertainty.
• Background simulation: Possible systematic effects due to the simulation of the dominant Z 0 /γ background are accounted for by replacing the PYTHIA Monte Carlo with a HERWIG sample. The two-photon background is removed and doubled as well. Variations in the extracted polarised crosssections, SDM elements and cos θ W distribution are taken as the systematic uncertainty.
• Lepton ID: The efficiency of the lepton identification in the data compared to that in the Monte Carlo is investigated as a function of the polar angle and energy of the lepton. The Monte Carlo is weighted to adjust these distributions to the data. Variations in the calculated polarised cross-sections, SDM elements and cos θ W distribution between the samples before and after weighting are taken as the systematic uncertainty.
• Detector effect correction: The use of Standard Model Monte Carlo to correct for detector effects could introduce some bias. This is tested for by comparing the helicity fractions calculated from generator level non-Standard Model Monte Carlo with those calculated from the same sample after full detector simulation and correction using correction factors determined from Standard Model Monte Carlo. This test is done with six samples of Monte Carlo each with one of the couplings ∆κ γ , ∆g z 1 , λ, set to ±1.
• Lepton energy scale and charge misassignment: Uncertainties on the lepton energy scale are tested by shifting the lepton energy by 0.3%. These are found to be negligible. The momentum resolution of electrons and muons and their charge misassignment are investigated by varying the resolution in Q/p t by 10%. Here Q is the lepton charge and p t is the transverse momentum with respect to the beam direction.

Conclusion
W-pair events with one leptonically and one hadronically decaying W are analysed to extract the polarisation properties of the W boson and test CP invariance. Both the individual W, and W-pair, polarised cross-sections have been measured. The individual W polarised cross-sections are welldescribed by the Standard Model predictions, both inclusively and separately for the leptonically and hadronically decaying W bosons (table 1). The results are consistent with those measured by L3 at the same centre-of-mass energy [15], with the results presented here being more precise.
The W-pair polarised cross-sections are measured for the first time at LEP. All the results are found to be consistent with the Standard Model expectations (table 2).
Triple gauge couplings are extracted using information from the leptonically decaying W, together with the W production polar angle. They are found to be consistent with the Standard Model predictions. The CP conserving couplings have been measured elsewhere with the same data sample [10], and serve as a consistency check for this analysis. The CP violating couplingsκ z ,λ z and g z 4 are measured for the first time with OPAL data, table 3, and constraints are thus placed on possible new CP-violating contributions to the W production and decay processes.     Figure 1: Distributions of the kinematic variables cos θ W , cos θ * l , cos θ * jet , φ * l and φ * jet , as obtained from the qqℓν events. The points represent the data. The histograms show the expectation of the Standard Model and the cases ∆g z 1 = ±0.5. The shaded histogram shows the non-qqℓν background. In the case of the W + decaying to the lepton, the value of φ * l is shifted by π in order to overlay the W + and W − distributions.