Measurement of Triple Gauge-Boson Couplings at LEP up to 189GeV

The triple gauge-boson couplings involving the W are determined using data samples collected with the ALEPH detector at mean centre-of-mass energies of 183 GeV and 189 GeV, corresponding to integrated luminosities of 57 pb^-1 and 174 pb^-1, respectively. The couplings, g^Z_1, Kappa_gamma and lambda_gamma, are measured using W-pair events, single-W production and single-gamma production. Each coupling is measured individually with the other two coupling fixed at their Standard Model value. Including ALEPH results from lower energies, the 95% confidence level intervals for the deviation to the Standard Model are -0.087<Dg^Z_1<0.141 -0.200<DKappa_gamma<0.258 -0.062<Lambda_gamma<0.147. Fits are also presented where two or all three couplings are allowed to vary. In addition, W-pair events are used to set limits on the C- or P-violating couplings g^V_4, g^V_5, Kappa_V, and Lambda_V, where V denotes either gamma or Z. No deviations from the Standard Model expectations are observed.


Introduction
The existence of the triple gauge-boson couplings (TGC) is a direct consequence of the SU(2) L × U(1) Y gauge theory. The measurement of the TGCs represents a fundamental test of the non-Abelian nature of the Standard Model. The triple WWγ and WWZ couplings have been studied at LEP in e + e − collisions at energies above the W-pair production threshold, using direct W-pair production (e + e − → W + W − ) [1,2], single-W production (e + e − → Weν) and single-γ production (e + e − → ννγ(γ)) [3][4][5]. Measurements of the TGCs have also been made at the Tevatron from studies of di-boson production [6]. This paper presents new results for the TGCs from analyses of W-pair, single-W, and single-γ final states using data recorded in 1997 and 1998 with the ALEPH detector. In 1997 and 1998 ALEPH recorded total integrated luminosities of 56.81 pb −1 and 174.20 pb −1 , at mean centre-of-mass energies of 182.66 GeV and 188.63 GeV, denoted as 183 and 189 GeV.
The most general Lorentz invariant parametrisation of the WWγ and WWZ vertices can be described by 14 independent complex couplings [7][8][9], 7 for each vertex: g V 1 , g V 4 , g V 5 , κ V , λ V ,κ V andλ V , where V denotes either γ or Z. Assuming electromagnetic gauge invariance, C-and P-conservation, the set of 14 couplings can be reduced to 5 parameters: g Z 1 , κ γ , κ Z , λ γ and λ Z , with Standard Model values g Z 1 = κ Z = κ γ = 1 and λ Z = λ γ = 0. Precision measurements at the Z resonance at LEP and SLC also provide bounds on the couplings [10,11]. However, local SU(2) L × U(1) Y gauge invariance reduces the relevance of these bounds [10] and introduces the constraints: where ∆ denotes the deviation of the respective quantity from its non-zero Standard Model value, and θ W is the weak mixing angle. Hence, only three parameters remain: ∆g Z 1 , ∆κ γ , and λ γ [9].
Using data from e + e − → W + W − final states all three couplings ∆g Z 1 , ∆κ γ and λ γ can be tested, whereas the single-W and single-γ final states allow measurements of only the WWγ-couplings, ∆κ γ and λ γ . Although the contribution from W-pair production dominates the combined limits, the single-W and single-γ events provide complementary information, which enhances the sensitivity especially for ∆κ γ .
In this analysis the three couplings ∆g Z 1 , ∆κ γ and λ γ are measured individually with the two other couplings fixed at zero, their Standard Model value. Fits are also presented, where two or all three couplings are allowed to vary.
The C-or P-violating sector of the TGCs is weakly bound. Indirect limits onκ γ , λ γ ,κ Z and g Z 4 exist, while there are no direct or indirect limits on the parametersλ Z , g γ 4 , g γ 5 and g Z 5 [11]. Only the parameterλ γ is tightly constrained by precision low-energy measurements [12]. This paper includes, for the first time, single-parameter fits to the unconstrained real and imaginary parts of the 8 couplings g V 4 , g V 5 ,κ V , andλ V , all zero in the Standard Model, based on an analysis of semileptonic (eνqq and µνqq) W-pair events.
The paper is organised as follows. In Section 2, a brief description of the ALEPH detector is given. The Monte Carlo event generators used in the analyses are presented in Section 3. The analysis of the single-γ final states is described in Section 4 and Section 5 is devoted to the single-W analysis. The description of the two analyses is rather concise, as they are presented in earlier publications [3,4]. In Section 6 the measurement of TGCs from W-pair 1 events is discussed in detail. Finally, all measurements are combined with ALEPH results from W + W − production at 172 GeV [1], single-W production at 183 GeV [3] and single-γ production at 183 GeV [4]. The resulting single-and multi-parameter fits are discussed in Section 7, followed by a summary and conclusions in Section 8.

The ALEPH Detector
A detailed description of the ALEPH detector can be found in Ref. [13,14]. The central part of the ALEPH detector is dedicated to the reconstruction of the trajectories of charged particles. Following a charged particle from the interaction point outwards, the trajectory is measured by a two-layer silicon strip vertex detector (VDET), a cylindrical drift chamber and a large time projection chamber (TPC). The three tracking detectors are immersed in a 1.5 T axial field provided by a superconducting solenoidal coil. Altogether they measure charged particle momenta with a resolution of δp T /p T = 6 × 10 −4 p T ⊕ 0.005 (p T in GeV/c). In the following, charged particle tracks reconstructed with at least one hit in the VDET, at least four hits in the TPC, and originating from within a cylinder of length 20 cm and 2 cm radius centred on the nominal interaction point and parallel with the beam, are referred to as good tracks.
Photons and electrons are identified in the electromagnetic calorimeter (ECAL), situated between the TPC and the coil. It is a lead-proportional-wire sampling calorimeter segmented in 0.9 • ×0.9 • towers read out in three sections in depth. It has a total thickness of 22 radiation lengths and yields a relative energy resolution of 0.18/ √ E+0.009, with E in GeV, for isolated photons. At low polar angles, the ECAL is supplemented by two calorimeters, LCAL and SiCAL, principally used to measure the integrated luminosity collected by the experiment. Electrons are identified by their transverse and longitudinal shower profiles in ECAL and their specific ionisation in the TPC. A detailed description of the photon identification can be found in [14].
The iron return yoke is equipped with 23 layers of streamer tubes and forms the hadron calorimeter (HCAL). The latter provides a relative energy resolution of charged and neutral hadrons of 0.85/ √ E, with E in GeV. Muons are distinguished from hadrons by their distinct pattern in HCAL and by the muon chambers composed of two double-layers of streamer tubes outside HCAL.
The information from the tracking detectors and the calorimeters are combined in an energy flow algorithm [14]. For each event, the algorithm provides a set of charged and neutral reconstructed particles, called energy flow objects, which are used in the analysis. Studies of Z → qq events show that the angular resolution of jets reconstructed from energy flow objects is typically 30 mrad in space and the energy resolution approximately σ E = (0.6 √ E + 0.6)(1 + cos 2 θ) GeV, where E is the jet energy in GeV and θ is the polar angle with respect to the z-axis along the e − beam direction.

luminosity.
The efficiency for the single-γ cross section measurement is estimated using a modified version of the KORALZ [15] Monte Carlo program. The KORALZ generator simulates initial state photons using YFS exponentiation [16]. The generator is modified to include the effects from photons produced as bremsstrahlung off the exchanged virtual W. This treatment includes the expected Standard Model contribution and possible anomalous couplings together with the interference of the two. The effect on the overall cross section is found to be small (∼0.2%) for Standard Model couplings. However, it can be as large as a few percent in certain kinematical regions. The predictions of the modified KORALZ Monte Carlo are confirmed by an independent generator NUNUGPV [17], which is based on exact lowest order amplitudes for the production of up to three photons in the final state, modified for higher order QED effects using transverse momentum dependent structure functions.
For the single-W study the GRC4F program [18] is used to simulate the four-fermion signal process final state eνff . The effective QED coupling constant is fixed to be α QED = 1/130.2 as suggested in [19]. For initial state radiation, the photon structure function approach is utilised. Final state radiation and tau decays are simulated with PHOTOS [20] and TAUOLA [21].
For the analysis of W-pair final states, the KORALW [22] generator, which includes all four-fermion diagrams contributing to W + W − -like final states, is used to produce the primary reference sample with a W mass of 80.35 GeV/c 2 . The KORALW generator is interfaced with JETSET [23], PHOTOS [20], and TAUOLA [21] for fragmentation, final state radiation and τ decays, respectively. In addition, several samples are generated using the double resonant CC03 [9] diagrams with non-standard values for one coupling at a time, to check the reconstruction and TGC determination. Finally, a sample generated with the double resonant CC03 diagrams is used to optimise selection efficiencies and parametrise the corrections used in the kinematic fitting.
In order to include the effects from various background processes, Monte Carlo samples are generated with a corresponding integrated luminosity of each background sample of at least 20 times that of the data. PYTHIA [23] is used to generate e + e − → qq(γ), ZZ, Zee, and e + e − → Weν event samples. In the ZZ sample, events with W + W − -like final states are discarded to avoid double counting. Two-photon processes are simulated with the PHOT02 [24] generator. The KORALZ [15] and UNIBAB [25] generators are used for the di-lepton final states.

Single-γ production analysis
Events with one or more photons and missing energy can be used to probe the anomalous WWγ coupling parameters ∆κ γ and λ γ . Although the single-γ channel is less sensitive to the couplings compared to the W-pair and single-W channels, it provides complementary information. A detailed description of the Standard Model processes involved in the reaction e + e − → ννγ(γ) and the modelling of the measured triple gauge-boson couplings can be found in [4].
The sensitivity to the WWγ couplings in the single-γ channel comes from the W-W fusion diagram. The W's exchanged in this t-channel diagram are predominantly at low momentum transfer. The single-γ channel is therefore mainly sensitive to ∆κ γ because contributions from λ γ contain higher powers of the W momenta. Furthermore, the effect of anomalous TGCs depends on the energy of the photon. For low energy photons, below the radiative return to the Z peak the sensitivity arises from the interference between the Standard Model and the anomalous contribution; this interference has a linear dependence on the TGCs. In the region around the radiative return to the Z peak, the sensitivity is minimal. For high energy photons, above the radiative return to the Z peak, the dependence on the TGCs is quadratic.

Event selection and determination of the TGCs
The events are selected from the 189 GeV data sample using the procedure described in [26]. In summary, single-γ events are selected by requiring at least one photon candidate with θ γ > 20 • and p γ T /E beam > 0.1 and no additional activity in form of reconstructed charged tracks or energy deposits in the forward regions (below 14 • ). Events where a photon has converted into a e + e − pair are not considered.
Anomalous contributions to the WWγ vertex increase the total cross section and lead to characteristic energy and angular distributions of the final state photons. For the singleγ channel the TGCs are extracted from the data by performing a maximum likelihood fit based on the overall number of observed photons, their polar angles θ γ and scaled energies where P i is the probability density function of observing event i with a given value of x E and θ γ and N exp is the expected number of events including background. The probability density function and the expected number of events for different values of the couplings are constructed by reweighting fully simulated single-γ events. Distributions of the polar angle and the scaled energy for single-γ events are shown in Figure 1. Two separate kinematic regions are used in the fit, excluding a region around the Z peak return, where the sensitivity is small.
The total numbers of photons used in the fit are 120 (128 expected) and 260 (258 expected) below and above the excluded region.

Results
At present energies, the cross section and the shape contribute equally in the likelihood function for ∆κ γ , whereas the result for λ γ is dominated by the sensitivity to the shape above the excluded region. The estimation of the systematic uncertainties follows the procedure described in Ref. [4] and the different contributions are summarised in Table 1. The fitted results for the 189 GeV data for single parameter fits, where each coupling is determined setting the other coupling to its Standard Model value, are where the first error is the statistical error and the second is the systematic uncertainty. The 95% confidence level limits including systematic errors are: The validity of these 95% C.L. limits and the error from the likelihood fit have been checked using many Monte Carlo samples corresponding to the data luminosity as described in [4].
Combining with the previous measurement for centre-of-mass energies between 161 and 183 GeV [4], the 95% C.L. limits on ∆κ γ and λ γ from single-γ production are The negative log-likelihood functions curves are shown in Figure 2 for the 189 GeV data, 161 -183 GeV data [4] and the combined results. In the combination, the systematic errors from acceptance and theoretical prediction are assumed to be fully correlated, while all other sources are taken as uncorrelated.

Single-W production analysis
Single-W production, e + e − → Weν, is sensitive to the WWγ vertex. This sensitivity comes from the γ-W fusion diagram, where the momentum transfer is low. As for the single-γ channel, the single-W channel is therefore mostly sensitive to ∆κ γ [27].

Selection
The analysis of single-W production is performed on the 189 GeV data sample. All W decay modes are used and the selection of each W decay final state, described in the following, has been optimised for the single-W signal definition used in the previous analysis at lower centre-of-mass energies [3]: θ e < 34 mrad, E ℓ > 20 GeV and | cos θ ℓ | < 0.95 for leptonic decays, M qq ′ > 60 GeV/c 2 for hadronic decays, where θ e is the polar angle of the scattered electron, E ℓ and θ ℓ are the energy and polar angle of leptons from the W decay. M qq ′ is the invariant mass of the quark pair. The cut angle at 34 mrad corresponds to the lower edge of the acceptance of the ALEPH detector. As single-W production is dominated by t-channel processes, the outgoing electron is predominantly emitted at small polar angles. Another specific feature is the large missing momentum carried away by the electron-neutrino, and therefore a common selection criterium for all single-W final states is the requirement of the missing momentum direction to be within the detector acceptance, | cos θ miss | < 0.9.

Leptonic Selection
The leptonic W decay is characterised by a high energy isolated lepton. Allowing for a multiprong decay of the tau, events with one or three good tracks (| cos θ| < 0.95) are accepted. The selection cuts are the same as in the analysis in Ref. [3] and are summarised here.
In addition to the cut on the missing momentum direction, tagged two-photon events are rejected by requiring that no energy be detected within a cone of 12 • around the beam axis (E 12 = 0).
The remaining backgrounds, mainly untagged two-photon events and two-fermion events, are eliminated by requiring that the transverse missing momentum be greater than 0.06 √ s. This threshold is increased to 0.1 √ s if the missing momentum direction points to within 10 • in azimuth to the boundaries between the two LCAL halves or between the six inner sectors for the TPC. It is required that no energy is found within a wedge of 10 • opposite to the direction of the lepton transverse momentum. To reduce the background from e + e − → Zee with Z decaying to neutrinos, events are rejected if an electron candidate track is identified and its energy, including the neutral energy in a 10 • cone around the track, is less than 20 GeV.
The selection efficiencies for the three final states are 75% (eνeν), 77% (eνµν) and 43% (eντ ν). The main background source is Zee where the Z decays to ν µνµ or ν τντ (the ν eνe case is a four-fermion final state which is e + e − → Weν like and is part of the signal).
In the data, 23 events are observed in agreement with the expectation from the Standard Model of 26.5 events (17.7 signal events). The composition is 15 events with an electron (8.4 signal and 6.2 background expected), 4 events with a muon (6.6 signal and 0.4 background expected) and 4 events with a tau (2.7 signal and 2.2 background expected). The distributions of the lepton transverse momentum, p l T , for single-W events passing the final selection cuts for the leptonic electron and muon W decay, are shown in Figure 3.

Hadronic selection
For the hadronic W decay, the event topology is characterised by two acoplanar jets with an invariant mass around that of the W boson. The selection is the same as in the analysis in Ref. [3]. In addition to the cut on missing momentum direction, at least seven good tracks are required. Similarly to the leptonic selection, tagged two-photon events and two-fermion events with initial state radiation are rejected by demanding that the energy E 12 be less than 0.025 √ s. The visible mass is required to exceed 60 GeV/c 2 and to be less than 90 GeV/c 2 to reject untagged two-photon events at the low end of the mass spectrum and ZZ events at the high end.
Events for which the energy in a wedge of 30 • centred on the transverse missing momentum direction is greater than 0.1 √ s are rejected. The acollinearity angle between the two hemisphere (defined by the event thrust axis) momentum directions is required to be less than 165 • . The semileptonic final state (ℓνqq) of W-pair production is efficiently rejected by requiring that no identified electron or muon with an energy of more than 0.05 √ s be reconstructed. The tau jet reconstruction algorithm of Ref. [3] is used in order to further reject semileptonic decays of W-pairs which contain a tau lepton. After all cuts, the semileptonic W-pair production, primarily events with one tau lepton, remains as the dominant background.
The efficiency for the hadronic W channel is about 43%. In the data, 53 events are observed, in agreement with the Standard Model expectation of 63.1 events (23.5 signal events). The visible mass distribution of the selected events is displayed in Figure 3.

Results
Limits on ∆κ γ and λ γ are derived from the total rate of single-W events, which is sensitive to the WWγ couplings. The upper limit on the single-W signal cross section has been calculated while varying only one coupling at a time, and the 95% C.L. limits on ∆κ γ and λ γ for the 189 GeV data are including the systematic uncertainties. The different contributions to the systematic errors are summarised in Table 2. The total systematic error amounts to 7% for the leptonic and 12% for the hadronic channel on the predicted numbers of signal events. The overall systematic errors are small compared to the statistical precision, which amounts to 33% for the leptonic channel and 55% for the hadronic channel.

W-Pair production analysis
The large number of W-pair events produced yield the dominant sensitivity to the TGCs. The process is sensitive to both the WWγ and the WWZ couplings via the s-channel Wpair production diagrams and the sensitivity to the coupling λ γ is higher than that of the single-W and the single-γ processes.

Event selection and kinematic reconstruction
In this section the event selections for the three distinct W + W − event topologies, ℓνqq, qqqq, and ℓνℓν, are described. Selected events are exclusively classified in the following order of priority: µνqq, eνqq, qqqq, τ νqq, and ℓνℓν. The expected numbers of events after all cuts used in the TGC results for signal and background processes at both centre-of-mass energies are summarised in Table 3 for each channel, along with the corresponding selection efficiencies and purities.
The event selection procedure for semileptonic W + W − events is similar to that used for the W mass measurement at the corresponding energy [28,29]. At 183 GeV, events are reconstructed such that they contain a high energy lepton candidate and two jets [28]. The charged particle with the highest momentum component anti-parallel to the missing momentum is chosen as lepton candidate. At 189 GeV the selection criteria for the lepton track are slightly changed, using the lepton track isolation [29]. The DURHAM-PE [28] clustering algorithm is applied to all energy flow objects not used to construct the lepton four-momentum, and these are forced into two jets. After this preselection, the probability for the event being signal is determined using the momentum of the lepton, the total missing transverse momentum and the lepton isolation from the closest jet. At this stage events passing a cut on the probability are subjected to a kinematic fit in order to improve the resolution on the reconstructed four-momenta of the W decay products. The kinematic fit and additional reconstruction cuts are described in the context of the specific TGC analysis.
A W-pair event can be characterised by the five measured angles, θ W , the W − production angle between the W − and initial e − in the W + W − rest frame, the polar and azimuthal angles of the lepton, θ * l and φ * l , in the rest frame of its parent W and the polar and azimuthal angles of a quark jet, θ * jet and φ * jet , in the rest frame of its parent W. The distributions of the five angles cos θ W , cos θ * l , φ * l , cos θ * jet and φ * jet , for eνqq and µνqq events at 189 GeV after selection and reconstruction is represented in Figure 5.
The event selection is based on two complementary algorithms, developed for the cross section measurement at 161 [30] and 172 GeV [31], but modified to account for the change in event kinematics with centre-of-mass energy. In summary, events passing a set of preselection cuts are selected as semileptonic τ candidate events if they fulfil either a global or topological 8 selection. The DURHAM-PE clustering algorithm is then applied to all energy flow objects that are not used to construct the tau four-momentum, and these are forced into two jets. More detailed descriptions of the selection and tau reconstruction can be found in the publications on the W mass determination at the corresponding energies [28,29].
To improve the resolution of the angular observables a 3-constraint kinematic fit is applied, requiring four-momentum conservation and reference mass constraints. In the kinematic fit the direction of the τ is approximated by its visible decay products and the extra energy loss is compensated by correction coefficients obtained from Monte Carlo simulated W + W − →τ νqq events. For single prong τ decays the charge of the τ is directly accessible, but in the case of three-prong τ decays ambiguities arise due to mis-assigned particles from the jets to the τ . For three-prong τ decays the charge of the τ is therefore determined from the sign of the pseudorapidity-weighted average jet charge of the τ decay products (Section 6.1.3), where the pseudo-rapidity is defined with respect to jet-axis. The charge mis-assignment in τ νqq events is 5% for one-prong and 41% for three-prong τ decays.
The distribution of the cosine of the W − production angle from τ νqq events can be seen in Figure 6.
To extract the hadronic W + W − signal with high purity and efficiency, the selection is based on a neural network [32]. Events passing a preselection designed to remove the qq(γ) background, are assigned a neural network output, based on global event properties, heavy quark flavour tagging, jet properties and WW kinematics. A detailed description of the selection algorithm at 183 and 189 GeV is given in [32].
For the hadronic W + W − events the reconstruction of the relevant information is more complicated since the W − direction is not known and the information on the particle flavours in either W system is not discriminant. In this case the four jets can be paired in three different ways. To select the best pairing, a 6-constraint kinematic fit is applied to all three possible pairings. The kinematic fit requires four-momentum conservation and reference mass constraints. The four-momenta obtained in the kinematic fit for the pairing with the lowest χ 2 value are then used in the final determination of the TGCs, while the other combinations are discarded. The efficiency of this algorithm to find the correct combination was found to be 78% at 183 and 75% at 189 GeV.
To assign a jet pair to the W + or W − a jet charge algorithm is used. The jet charge, Q jet , is obtained from the pseudorapidity-weighted average charge of jet particles. The jet pair charge is defined by the sum of the two jets assigned to a W, Q W = Q jet1 + Q jet2 . A jet pair is then assigned to the W + based on the charge difference between the two jet pairs, ∆Q, with a probability P + . The probability P + is given by where N + is the probability density function for the charge difference between the two W systems for true W + jet pairs obtained from Monte Carlo event samples [1]. Figure 7 shows the distribution of the di-jet charge of the two W systems for true W + and W − decays, obtained from Monte Carlo generator information. The distribution of the di-jet charge from semileptonic events is also shown for both data and Monte Carlo. The data are well reproduced by the Monte Carlo simulation. The charge assignment efficiency for correctly paired hadronic W + W − events amounts to approximately 76% for P + > 0.5.
The distribution of the cosine of the W − production angle from qqqq events can be found in Figure 6.

W + W − →ℓνℓν events
The selection of W + W − →ℓνℓν events (where ℓ denotes an electron or muon) is mainly based on variables used in Ref. [30], namely missing transverse momentum, missing mass and kinematic properties of the lepton candidate. The information of these and other variables is combined in a neural network. A detailed description of all the variables used in the neural network can be found in the appendix A.
For purely leptonic W + W − events the momenta of the two neutrinos are unknown. However, in the absence of ISR and neglecting the W width, the constraint that the two lν systems should have the W mass (M 1 = M 2 = 80.35 GeV/c 2 ) in combination with the usual four-momentum conservation allows a reconstruction of the neutrino momenta. The quadratic nature of the mass constraint results in a two-fold ambiguity, corresponding to flipping both neutrinos with respect to the plane defined by the charged leptons. As detector resolution, ISR and the finite W width are not included in this reconstruction hypothesis, 28.0% of the events have no physical solution and a zero-constrained kinematic fit is employed. The fit determines a set of values for the reference masses, M 1 and M 2 , in the mass-constraints for the two lν systems, for which a physical solution exists. By this method 92% of the events without a solution are recovered, resulting in a 97.7% reconstruction efficiency for signal events. The majority of the events which fail to have a solution are purely leptonic W + W − events with at least one leptonically decaying τ , which is the dominant background. In addition to the selection by a cut on the neural network output, purely leptonic events are only accepted if a physical solution is found inside the mass window of 55 − 105 GeV/c 2 .
The distribution of the cosine of the W − production angle for ℓνℓν events at 189 GeV after selection and reconstruction is represented in Figure 6.

Determination of the TGCs
Three different methods, described in the following, are used to extract the couplings in the different W-pair final states.

The optimal observable methods
The general idea of optimal observables (OO) [33] is to project the sensitive kinematic information for a given coupling g i onto the one-dimensional distribution of a suitably defined variable O (1) i . The coupling g i can then be extracted from a fit to this distribution or, equivalently, from the measurement of the mean value O (1) i of the optimal observable. Since the amplitudes are linear in the TGCs the differential cross section can be expanded in these couplings g i containing no terms beyond the second order where g i denotes any type of couplings and Ω denotes phase space variables taking into account reconstruction ambiguities for the individual W + W − channels. The zero-order term, S 0 (Ω), represents the Standard Model contribution. Using the first order term, a given set of couplings, g, can be determined by minimising is re-optimised.
• Adding the information contained in the second order term of the expansion in Eq. 3 (called OO 2 in the further text). This is achieved by including the second optimal observable and adding additional terms of the same structure as the ones in Eq. 4 to the χ 2 , including new terms describing the correlation between O ij . The second order observable increases the sensitivity when the information contained in the first order observable decreases [34]. The covariance matrix for the mean values are obtained by reweighting fully simulated W + W − events.
In both cases the information from the measured cross section is included by adding a Poisson term to the likelihood function.
By construction, these methods are bias-free and take into account any experimental effect, provided that the Monte Carlo simulation describes the data correctly. For a given channel, contributions from any other channel are considered as coupling dependent backgrounds. In addition the efficiencies and purities of each selection are parameterised as function of the couplings.
For semileptonic events, both OO analyses, OO 1 and OO 2 , apply a two-constraint kinematic fit using four-momentum conservation, equal W-mass hypothesis and including a massless neutrino. The corresponding χ 2 -probability of the fit is required to be larger than 0.005 for an event to be selected for the TGC extraction. This cut improves the purity of the sample and discards poorly reconstructed events. In addition, an OO-window cut is applied, optimised to improve the sensitivity to the TGCs by reducing the contamination of background events with OO values incompatible with W-pair production.
For hadronically decaying W's, there remains a twofold ambiguity since the quark flavours are undetermined. Hence, for semileptonic eνqq, µνqq and τ νqq events the contributions are averaged over the quark and anti-quark directions. For hadronic events, considering the W charge as undetermined, there is an eightfold ambiguity. This is included in the final extraction of the TGCs, where each contribution is weighted with the corresponding di-jet charge probability, Eq. 2. For leptonic events, ℓνℓν, the contributions are averaged over the two solutions for the neutrino momenta.

Maximum likelihood-fit
A maximum likelihood analysis (LL), Eq. 1, of the channels eνqq and µνqq is used to measure the C-or P-violating couplings and as a cross-check for the CP-conserving couplings. The measured variables are the five angles described in Sect. 6.1.1. As in the optimal observable analyses, no quark flavour tagging is performed and the quark and anti-quark directions are averaged. The probability density function, P , is given by where g denotes a set of couplings and the anglesΩ = (θ W , θ * l , φ * l , θ * jet , φ * jet ) are calculated using the charged lepton, neutrino, and quark jet four-vectors. The quantity b(Ω) is the background distribution as predicted by Monte Carlo. The W + W − signal distribution, s(Ω,ḡ), is defined by whereŝ is the squared invariant mass of the W + W − system, r(Ω,Ω true ) is the detector resolution function, ǫ(ŝ,Ω true ) is the detection efficiency, F (ŝ) is an initial state radiation function [35], and dσ dΩ ′ (ŝ,Ω true ,ḡ) is the lowest-order narrow-width differential cross-section for W + W − production and decay [7].
The normalisation factors B and S(ḡ) are the integrals of b(Ω) and s(Ω,ḡ), where S(ḡ) is evaluated by reweighting W + W − Monte Carlo events in order to include detector resolution and efficiency. The proper evaluation of the normalisation S(ḡ) is crucial for the success of the likelihood method.
Some approximations to s(Ω,ḡ) are made when evaluating the numerator of Eq. 5. In particular, most of the information in the detection efficiency function ǫ(ŝ,Ω) arises from the charged lepton momentum p l and polar angle θ l in the laboratory reference frame. This dependence has been parameterised with a two-dimensional efficiency function.
The event selection for this method is the same as for the optimal observables. A kinematic fit using the equal mass hypothesis and four-momentum conservation is applied, where the corresponding χ 2 -probability is required to be larger than 0.02. The efficiencies, purities and numbers of events for data and Monte Carlo simulation are shown in Table 4.

Comparison of methods
The three methods (OO 1 , OO 2 and the maximum likelihood) discussed above are used to determine the TGCs from WW events. The OO 1 method is used to measure the couplings ∆g Z 1 , ∆κ γ and λ γ for all five WW final states considered in this analysis, namely eνqq, µνqq, τ νqq, qqqq and ℓνℓν. The OO 2 method is employed to measure the couplings ∆g Z 1 , ∆κ γ and λ γ for the eνqq, µνqq and qqqq channel. The maximum likelihood method is used in the two semileptonic channels, eνqq and µνqq, to measure the standard set of couplings ∆g Z 1 , ∆κ γ and λ γ and the real and imaginary parts of the C-or P-violating couplings g V 4 , g V 5 ,κ V , andλ V , where V denotes either γ or Z. A comparison of the three couplings ∆g Z 1 , ∆κ γ and λ γ , using the three methods OO 1 , OO 2 and maximum likelihood, is given in Table 5. For OO 1 and OO 2 , the combined results at 183 and 189 GeV are extracted by adding up the corresponding χ 2 terms, while for the maximum likelihood method results are extracted by summing up the corresponding log L functions. The error intervals for each coupling are defined as the 68% confidence intervals obtained by integration of the likelihood functions, to accommodate cases with non-parabolic behaviour of the log-likelihood function.
The linearity of the three fitting procedures is checked by repeating the fits using Monte The consistency of the Monte Carlo simulations with the data is verified by comparisons of the distributions of the input quantities to the selections for data and Monte Carlo. In addition, the stability of the analysis with respect to the event selection is tested by varying the main selection criteria within reasonable limits. In no case significant discrepancies are found.
The reliability of the errors from each fitting procedure is investigated by performing fits to a large number of independent Monte Carlo samples, each corresponding to the integrated luminosity of the data. These samples, typically 300, are then reweighted to non-zero values for the TGCs and passed through the analysis chain. The expected 68% confidence levels, obtained from the distributions of the fit values, show good correspondence with the 68% confidence intervals obtained for data.
The consistency of the data results from the three different methods used for the eνqq and µνqq channels, has been checked by performing fits to a large number of independent Monte Carlo samples. The results obtained with the three different methods are compatible and the expected spread between the methods show a good agreement with the observed differences in the data.
The expected 68% confidence level intervals, obtained from the distributions of the fit values for the three methods at 189 GeV for the three couplings ∆g Z 1 , ∆κ γ and λ γ , are listed in Table 6. The expected errors of the OO 2 and the maximum likelihood method are very similar, whereas the OO 1 method is slightly worse in the case of ∆κ γ . For the final results, the OO 2 method is therefore employed in the analysis of the eνqq, µνqq and qqqq final states and the OO 1 analysis is used for the remaining τ νqq and ℓνℓν final states.

Systematic uncertainties
In the following the different sources of systematic errors for each decay channel and their determination for the various methods are briefly described. The different contributions of each source to the total systematic error for the three couplings ∆g Z 1 , ∆κ γ and λ γ , as 14 obtained with the OO methods, are given for each channel in Table 7. The systematic uncertainties for the combined eνqq and µνqq channels are listed in Table 8 for the C-or P-violating couplings, determined with the maximum likelihood method.
The following sources, listed in the approximate relative importance, have been considered to be fully correlated between the channels: -Fragmentation: The effect of fragmentation in hadronic W decays is estimated by fitting the couplings in samples of events generated with KORALW where the default JETSET fragmentation is replaced by HERWIG [36] fragmentation. The HERWIG fragmentation parameters are tuned at the Z using hadronic events with flavour tagging [29].
-W + W − cross section: The uncertainty due to the theoretical error on the expected W + W − cross section predicted by KORALW is estimated by changing the W + W − cross section by ±2% [9].
-Luminosity: The effect of the error on the integrated luminosity is estimated by varying the measured integrated luminosity by ±0.7% [37].
-LEP energy: The uncertainty on the LEP energy affects the determination of the couplings via the kinematic fitting procedure and the cross section measurement. The values of the LEP energy are varied in the range ±0.050 GeV [38], which has a negligible effect on the results.
-W mass: The analysis is repeated using Monte Carlo samples generated at different values of the W mass to investigate the effects due to the uncertainty of ±62 MeV in the W mass measured at hadron colliders [39].
-Calorimeter absolute scale: The absolute energy scale of the electromagnetic and hadronic calorimeters is determined using hadronic Z events. The uncertainties in the absolute scale are found to be ±0.9% and ±2% for the electromagnetic and hadronic calorimeter, respectively. The effect of a possible miscalibration of the calorimeters is evaluated on Monte Carlo samples by scaling the electromagnetic and hadronic part of the measured energy independently by these amounts. The largest of the observed shifts for each calorimeter is combined in quadrature.
-Particle tracking: The definition of a good charged track has been tightened in this analysis in order to minimise possible effects from residual tracking distortions primarily in the forward regions of the detector. Corrections for the distortions are determined by studying Z→ µ + µ − events, and possible remaining distortions have been estimated using Bhabha events. The systematic uncertainty related to tracking is assessed by applying the corrections and adding the remaining distortions independently to Monte Carlo event samples and repeating the analysis [28].
-Jet energy corrections: Detailed comparisons of reconstructed jets in Monte Carlo and data are used to parametrise small corrections to Monte Carlo jet energies as function of the jet polar angle to the beam axis [28]. In order to evaluate the effect from the uncertainty in the Monte Carlo jet energy corrections, two alternative correction functions, corresponding to ±1σ errors of the discrepancies, are used. The largest shift with respect to the nominal correction is taken as the systematic error [28].
-Higher order terms: The effect from missing higher order terms, O(α 3 ), in the simulation of initial state radiation in the KORALW generator, is assessed following the procedure described in Ref. [28]. In summary, the error on the couplings is determined by comparing fits of Monte Carlo samples with events weighted to O(α)/O(α 2 ) with fits to the corresponding unweighted samples, which have been generated in the second order leading-log approximation. Recently, new improved calculations with the Double Pole Approximation [40] have appeared. The improved CC03 cross section and angular distributions predicted by two independent Monte Carlo programs, RacoonWW [41] and YFSWW [42], could introduce small changes on the couplings. For the time being these effects have not been included.
Errors assumed to be uncorrelated between channels include: -Bose-Einstein correlations: The effect of Bose-Einstein correlations in the qqqq channel is investigated by repeating the analysis on Monte Carlo events generated with KORALW and fragmented using JETSET with Bose-Einstein correlations for all particles, following the implementation in LUBOEI [43]. The scheme for restoring four-momentum conservation denoted BE 3 , which has been tuned to the LEP1 Z data, is considered [29].
-Colour reconnection: The uncertainty arising from possible colour reconnection effects is assessed by studying Monte Carlo implementations of different colour reconnection scenarios in the parton evolution scheme in JETSET [44]. The analysis is repeated with qqqq events generated with the EXCALIBUR generator and hadronised with and without colour reconnection in the model referred to as SK1, as described in [28]. The systematic error is taken as the difference in fitted couplings from samples without colour reconnection and with colour reconnection in about 30% of the events.
-Background estimation: The error on the couplings from the uncertainties in the background estimation is evaluated by varying the normalisation of the main background processes. The background from QCD is changed by ±5% based on comparisons between data and Monte Carlo simulation. The background from γγ, Zee and ZZ processes is varied by ±30% ±20% and ±2%, respectively, to account for the theoretical uncertainty in the prediction for those processes [9].
-Monte Carlo statistics: The effect of the Monte Carlo statistics is included in the systematic uncertainty.
-Jet charge assignment: To investigate the effects from the uncertainties on the jet charge, the reconstructed W charge is shifted by 0.01. This number is based on comparisons between data and Monte Carlo simulation (Z peak data) [45]. The effect on the couplings is found to be negligible.
The systematic uncertainties listed above as fully correlated between channels are also assumed to be fully correlated between years. In addition, the systematic errors from Bose-Einstein correlations and colour reconnection are taken to be to fully correlated between years.
For both optimal observable methods, OO 1 and OO 2 , the systematic errors have been calculated based on the changes in the mean values of the respective observables. This is incorporated in the TGC extraction by including the corresponding covariance matrix for   the systematic uncertainties. By this procedure the systematic uncertainties are folded with the proper statistical correlations between the optimal observables and the results of the fits include both the statistical and systematic errors. The systematic uncertainties listed in Table 7 are derived from the changes in the optimal observables mean values. They are not used as such in the analysis but serve only as a representation of the systematic contributions from the different sources. The systematic uncertainties for the maximum likelihood method are convoluted into the log L functions by assuming parabolic behaviour of the systematic errors around the fitted TGC value.

Final results from W-pair production
The combined results from all W + W − decay channels at 183 and 189 GeV for the three couplings ∆g Z 1 , ∆κ γ and λ γ , are obtained by combining the OO 2 analysis of the eνqq, µνqq and qqqq final states with the OO 1 analysis of the τ νqq and ℓνℓν final states. The correlation of the systematic errors between the different channels and energies are included as described in Section. 6.3.2. The results for ∆g Z 1 , ∆κ γ and λ γ , including systematic uncertainties, are listed in Table 9. The final 68% and 95% combined W-pair result for the three couplings ∆g Z 1 , ∆κ γ and λ γ is summarised in Table 10. The corresponding log L curves, including systematic uncertainties, are shown in Figure 8.
The maximum likelihood method is used in the two semileptonic channels, eνqq and µνqq, to measure the real and imaginary parts of the C-or P-violating couplings g V 4 , g V 5 ,κ V , andλ V , where V denotes either γ or Z. The combined 183 and 189 GeV results, including the systematic uncertainties, for the C-or P-violating couplings are summarised in Table 10.
In all cases described above, each coupling is determined fixing the other couplings to their Standard Model values. The error intervals for each coupling are defined as the 68% confidence level intervals obtained by integration of the likelihood functions, to accommodate cases with non-parabolic behaviour of the log-likelihood function.

Combined TGC results
The measurements from single-γ, single W and WW production, are combined with previous ALEPH results from W + W − production at 172 GeV [1], single-W production at 183 GeV [3] and single-γ production at 183 [4]. The combined results are listed in Table 11. In Figure 9 the corresponding one-parameter log L curves are shown.
To study the full correlation between the parameters, two-and three-parameter fits, where two or all three couplings are allowed to vary, are also presented. The fits use the combined information from W-pair production, single-W production and single-γ production at 183-189 GeV.
For the three parameter fit the results and the errors computed from a variation from the minimum of the log L functions of 0.5, are summarised in Table 12 including the systematic uncertainties. The correlation matrix of the three-parameter fit is also given in Table 12. This correlation matrix is evaluated at the local minimum, and the correlations vary substantially depending on the exact value of the minimum. The projections onto the two dimensional plane of the three dimensional envelope of the 95% confidence level volume, representing the integration of the confidence over the corresponding third coupling, are shown in Figure 10. The 95% confidence limits of the respective 2-parameter fits of the three pairs of couplings (∆g Z 1 , ∆κ γ ), (∆g Z 1 , λ γ ) and (∆κ γ , λ γ ) are shown as full lines. The systematic uncertainties are included in the limits shown. No deviations from the Standard Model expectations are observed.

Summary and conclusions
The triple gauge-boson couplings have been measured using W-pair events at 183 and 189 GeV, single-W production at 189 GeV and single-γ production at 189 GeV. Combining with previous ALEPH results from W + W − production at 172 GeV, single-W production and single-γ production at 183 GeV, the three couplings ∆g Z 1 , ∆κ γ and λ γ have been measured individually, assuming the two other couplings to be fixed at their

A Leptonic Neural Network Input Variables
The Neural Network (NN) calculates an approximation of the multidimensional probability density function in the following 13 input variables, for signal and backgrounds. The NN is applied, after preselection, to events with at least two opposite charged tracks with momentum -after bremsstrahlung correction -in excess of 15 GeV. The NN uses variables related to the lepton candidates, to the missing momentum, global event variables, and WW kinematics. They are listed here together with their relative discriminating power, namely the statistical correlation with the neural network output: • missing mass squared (13.7%); • missing transverse momentum (11.5%); • angle between the two most energetic tracks (9.6%); • energy of the second most energetic track, (8.7%); • total energy found in a 12 o cone around the beam axis (8.2%); • number of identified leptons with an energy greater than 15 GeV(8%); • missing transverse momentum with respect to the plane defined by the beam axis and the 3D-thrust axis (7.8%); • energy of the most energetic track (7.4%); • invariant mass of the two most energetic tracks (6.6%); • missing longitudinal momentum (6.3%); • scalar sum of the transverse components of the two most energetic tracks with respect to a 2D-thrust axis, built from the projection of the track momenta on the transverse plane (5.4%); • number of isolated neutral clusters with energy more than 4 GeV outside a cone of 10 o around each of the two most energetic tracks and forming an invariant mass with each of them of more than 2 GeV (4.1%); • cosine of the angle between the most energetic track and the axis perpendicular to the plane defined by the second most energetic track and the z-axis (2.6%).  Figure 10: Multi-parameter fits using the combined data from single-γ, single-W and W-pair production at energies up to 189 GeV. The two-dimensional 95% confidence level contours for the three pairs of couplings, a) (∆g Z 1 , ∆κ γ ), b) (∆g Z 1 , λ γ ) and c) (∆κ γ , λ γ ). The solid lines show the 95% confidence level contours of the two-parameter fit. The shaded area is a projection onto the two-dimensional plane of the three-dimensional envelope of the 95% confidence level volume. The Standard Model point is represented by a star.