Measurement of the Tau Polarisation at LEP

The polarisation of $\tau$'s produced in Z decay is measured using 160 pb$^{-1}$ of data accumulated at LEP by the ALEPH detector between 1990 and 1995. The variation of the polarisation with polar angle yields the two parameters ${\cal A}_e = 0.1504 \pm 0.0068 $ and ${\cal A}_{\tau} = 0.1451 \pm 0.0059$ which are consistent with the hypothesis of $e$-$\tau$ universality. Assuming universality, the value ${\cal A}_{e{-}\tau} = 0.1474 \pm 0.0045$ is obtained from which the effective weak mixing angle $\sin^2 {\theta_{\mathrm{W}}^{\mathrm{eff}}} =0.23147 \pm 0.00057 $ is derived.


Introduction
Owing to parity violation in Z production and decay, the τ leptons produced in the e + e − → Z → τ + τ − reaction are polarised. Except for very small O(m 2 τ /m 2 Z ) corrections, the τ + and τ − have opposite helicities, hence opposite polarisations. In the present paper, P τ is defined as the τ − polarisation.
The P τ dependence on the angle θ formed by the τ − direction and the e − beam is given, in the improved Born approximation at the Z peak, by P τ (cos θ) = − A τ (1 + cos 2 θ) + A e (2 cos θ) (1 + cos 2 θ) where A f b is the forward-backward charge asymmetry of τ production and In the standard model, the vector (g l V ) and axial vector (g l A ) couplings of the Z to lepton l are independent of the lepton flavour and related to the effective weak mixing angle by Therefore, the measurement of P τ (cos θ) allows a test of the e-τ universality prediction A τ = A e and, assuming universality, gives a determination of the weak mixing angle.
The analyses presented here use the complete set of data (160 pb −1 ) accumulated by the ALEPH detector at LEP I from 1990 to 1995. They supersede the previously published results [1][2][3] The τ polarisation is extracted from five decay channels which amount to a total branching ratio of about 90%, namely τ → eνν, τ → µνν, τ → πν, τ → ρν, τ → a 1 ν. The channels τ → Kν and τ → Kπ 0 ν are included in τ → πν and τ → ρν respectively, and the a 1 decay is observed as both a 1 → π + π − π ± (3h) and a 1 → π 0 π 0 π ± (h2π 0 ). The methods employed for the P τ measurement in each channel are described in the following section. At each energy, a fit of Eq. 1 allows to determine uncorrected parameters A τ and A e . Final corrections are needed to get the effective couplings at the Z, whose combinations 2g l V g l A /[(g l V ) 2 + (g l A ) 2 ] are written A l to distinguish them from the uncorrected A l . They are computed with the help of the ZFITTER program [4].
2 Methods for τ polarisation measurement 2.1 Information from a polarised τ decay The general method to measure the τ polarisation in an optimal way [5] takes advantage of the linear dependence on the polarisation of the distribution of decay products.
In each decay channel, the decay products are described by a set of n observables ξ and their distribution reads 1 with the normalization and positivity conditions f (ξ)d n ξ = 1, g(ξ)d n ξ = 0, f ≥ 0 and |g| ≤ f .
With the optimal variable for the channel defined by the distribution in Eq. 4 can be cast, with no loss of information, in the general reduced formŴ where W + and W − are the distributions for positive and negative helicity respectively.
For events in which both τ 's decay into hadrons, it is possible to get information on the line of flight of the two τ 's [6]. It is shown in Ref. [5] that, in an ideal measurement, using the knowledge of the τ direction to construct the ω variable gives for ρ and a 1 decays the same sensitivity as that for π decay ( Table 1).
The polarisation is obtained by means of a fit ofŴ + (ω) andŴ − (ω) functions to the experimentalŴ (ω) distribution. The W ± functions are constructed from Monte Carlo events and normalized to the acceptances. Non-tau background is taken into account in the fit.
The validity of the standard model for the description of τ decays has been checked elsewhere [7] by the measurement of the τ decay parameters.
In the standard model, which is assumed here, the decay distributions for the τ → lνν, τ → πν, and τ → ρν channels are completely determined by Lorentz invariance and there is no hadronic decay model dependence in the definition of ω. Furthermore, because the τ polarisation always appears through the product χP , where χ is the handedness of the neutrino [2], the definition of the ω variable is the same for τ + and τ − , as are its distributions when expressed in terms of the τ − polarisation P τ = P τ − = −P τ + .
The situation is more intricate in the case of the a 1 decay. There, in addition to the decay angles, the two-pion masses are required to describe the hadronic system. Their dependence is embodied in four "structure functions", W A , W C , W D , and W E [8], whose computation requires a model of the a 1 decay. The model [8] that is used in the present study assumes a sequential a 1 → ρπ decay. The sign of the W E function Table 1: Ideal sensitivities for the polarisation measurement in the τ decay channels without and with the τ direction τ [5]. The sensitivity is defined as S = 1/∆ √ N , where ∆ is the statistical error expected for a sample of N events.

Channel
Sensitivity without τ with τ πν 0.58 0.58 ρν 0.49 0.58 a 1 ν 0.45 0.58 lνν 0.22 − 2 depends on the charge of the τ so that the expression of the ω variable is not exactly the same for τ + and τ − .
For the τ → lνν and τ → πν decay modes, the optimal variable is the ratio x = E/E beam of the energy of the charged particle (l or π) to the beam energy.
For the other channels, the τ decay is described by one angle and the mass of the hadronic system. The subsequent decay of the hadron is described by two angles in the case of the ρ, and by two effective masses and three Euler angles in the case of the a 1 . If the τ direction is not available, one of the angles used to describe the hadron decay is not measurable and the analytical form of the decay distributions is made more complex by the Wigner rotation between the laboratory and τ rest frame helicity axes. The expressions used for the construction of the ω variables can be found in Refs. [2] and [8].
Although all the information relevant to the polarisation measurement is contained in the sole ω variable, the distributions of the other parameters are useful when checking the quality of the simulations and the understanding of energy calibrations. For example, the two angles used to describe the τ → ρν decay when the τ direction is not available are [2] cos ψ τ ∝ (2(E π ± + E π 0 )/E beam − 1) and cos ψ ρ ∝ (E π ± − E π 0 )/(E π ± + E π 0 ). The comparison of their distributions in data and Monte Carlo is the best test of the quality of the measurement.

The analyses
Two complementary analyses have been performed, which emphasize differently the sensitivity of the estimators and the reduction of the systematic uncertainties.
In the first one, the information from the opposite hemisphere of the event is used for background rejection when studying a τ decay. This approach, which was followed in previous analyses [1][2][3], is called hereafter the "single τ method". In this method, the information exploited to construct the ω variable comes from only one of the two hemispheres defined by the plane perpendicular to the thrust axis.
As mentioned in the previous section, for events in which both τ 's decay into hadrons it is possible to get information on the line of flight of the two τ 's and construct an ω variable giving an improved sensitivity for the τ → ρν and τ → a 1 ν decays. This more global approach is followed in the second analysis and is called hereafter the "τ direction method".
These two methods use the same set of data and many common standard ALEPH tools. The two analyses have nevertheless been performed independently and thus are described in detail in Sections 4 and 5. The combination of their results is presented in Section 6.
3 Detector and data sets 3

The ALEPH detector
The ALEPH detector and its performance are described elsewhere [9,10]. The study of τ polarisation requires a very special care in the alignment and calibration procedures to master the differences between real and simulated data at the required level.
The charged track measurement rests on three elements, the Vertex Detector (VDET), the Inner Tracking Chamber (ITC), and the Time Projection Chamber (TPC).
The VDET, a double-sided silicon strip detector, provides precise measurements close to the interaction point and, therefore, plays an important role in determining the τ direction.
The ITC is a drift chamber with wires parallel to the beam and a short drift time. Besides providing a trigger, it is an efficient tool against cosmic rays. Because its complete drift time is about 200 ns, to be compared to the 11 µs interval between bunch collisions, its sensitive time window is small. Therefore, there are few if any ITC hits in most triggers on cosmic rays. If present, they are not reconstructed with the correct timing and the reconstructed track is distorted leading to erroneous distances of approach to the interaction point.
The TPC is a very large volume drift chamber. Together with the VDET and ITC, it brings a resolution on the track momentum which reaches 6 × 10 −4 p (p in GeV), but to get this accuracy a precise calibration is needed. For these analyses, the distortions from the TPC are derived from the study of nonradiative µ + µ − events in φ and cos θ bins, where φ and θ are the azimuthal and polar angles of the µ − with respect to the e − beam. The departures from the corresponding observations on Monte Carlo are corrected for, considering that the distortions on the sum of positive and negative track momenta are due to field effects and that the distortions on the difference are due to sagitta measurement effects, i.e. alignment effects.
Apart from momentum and angle measurements, the tracking system is used to measure the impact parameter relative to the beam axis, denoted d 0 .
The tracking detectors are surrounded by the electromagnetic calorimeter (ECAL). All these elements are installed inside a superconducting coil which provides a 1.5 T axial field. The return yoke of the magnetic field is instrumented to form a hadron calorimeter (HCAL).
The ECAL is a lead/proportional wire chamber sampling device of 22 radiation length thickness. The insensitive regions between its modules (cracks) represent 2% of the barrel and 6% of the endcap areas. The anode wire signal is read plane by plane. The cathodes are divided in pads, making 74,000 towers pointing to the interaction point. These towers are read out in three sections in depth ("storeys"or "stacks"). The use of the energies measured with the pad towers is mandatory whenever a precise direction is needed but, in the case of τ events, the energy collected on the wires may be advantageous. It does not suffer from disconnected channels ensuring that the total energy left in each module is properly known. This comes from the fact that the lack of a plane is, to first order, corrected for by the energy calibration. It has a very low level of noise (about 10 MeV for a complete module). The energy resolution therefore does not suffer from threshold effects and is slightly better than the resolution of the pad signal. Another interesting feature is that the 45 plane signals provide a complete depth profile of the energy deposit.
The calibration of the ECAL was performed independently in the two analyses presented below but the salient features are common and the technique is the same for pad and wire signals. First a correction map is applied which takes care of local gain variations related to mechanical distortions, then gain dependencies as a function of φ and θ are applied to each module, finally the modules are inter-calibrated and normalized to the response observed in Monte Carlo. A problem in this procedure comes from the nonlinearity with energy of the calorimeter response which reaches 4% at 50 GeV. Bhabha events provide a calibration at the beam energy but not a measurement of the nonlinearity which can be parametrized in first order by a constant α in the expression giving the true energy as a function of the measured one, E t = E m (1 + αE m ). To extract α, clean electrons are selected and the E/p ratio is compared, as a function of energy, between real data and Monte Carlo where the saturation effects have not been simulated. This is done independently for barrel and endcaps.
The HCAL is made of 23 layers of streamer tubes inserted between 5 cm thick iron plates. The individual tube hits are recorded digitally and the energy is read out in 4,788 projective towers. This calorimeter is used primarily for discriminating between pions and muons. This does not require a specific energy calibration. For the current purpose, use is made essentially of the digital pattern provided by the firing of the streamer tubes. In particular, the pattern recognition described in Ref. [2] is used to get variables describing the shape of the shower.
The HCAL is surrounded by two double layers of streamer tubes, the muon chambers, providing additional information for muon identification.
The luminosity calorimeters are not utilized in the analyses presented here.

The data and Monte Carlo samples
The data used in these analyses were accumulated with the ALEPH detector during the years 1990 to 1995 at centre-of-mass energies close to the Z mass. This set, often referred to as the LEP I data set, corresponds to about 160 pb −1 of data. A detailed description of the characteristics of each period can be found in Ref. [11]. For the present study no requirement on the luminosity measurement quality is made.
Monte Carlo simulation samples were generated for the τ + τ − events as well as for the different backgrounds through the complete detector simulation chain. A sample of 1.9 × 10 6 τ + τ − events was generated using KORALZ [12]. For each year and energy of data taking a sample of about six times the real data sample was simulated with the same energy and the effective detector geometry of the year.

Single τ analysis 4.1 Introduction
The aim of this analysis is to reduce the systematic uncertainties, even at the price of some inefficiency. Data have been processed separately for each year of LEP operation in order to take as much care as possible of time-variations in the behaviour of the detector.
Particle and channel identification, as well as non-τ background rejection, are achieved by means of likelihood estimators, in order to optimize the efficiencies. These estimators are constructed from Monte Carlo distributions of the relevant variables, but their efficiencies are directly measured on selected samples of data, which are used also to study the properties of the remaining backgrounds.
Specific procedures have been designed for the γγ and cosmic ray background rejection. They are described in detail in the following.
To cross-check the photon identification and reconstruction, an inclusive one-prong hadronic analysis, which infers the polarisation from the charged pion momentum distribution, has been devised, in addition to the standard study of the previously listed channels.

Tools
The general tools applied in the analysis are presented here. Additional information on their performances is given in Section 4.4 dedicated to the study of the systematic uncertainties.

Charged particle identification
The charged particle identification procedure (PID) described in Ref. [11] is used for the non-τ background rejection and for the decay channel selection.
Charged particles are classified as electron, muon, or hadron (pions and kaons are treated as a single class thereafter referred to as "pions"). Since there is negligible probability of mistaking electron and muon (∼ 2.5×10 −4 ), the PID is based on two likelihood estimators which give the relative probability for a track to be an electron versus a pion (e/π), and to be a muon versus a pion (µ/π), respectively. The µ/π PID is restricted to charged tracks with momentum above 1.3 GeV/c. An identification efficiency of about 99% with a rate of misidentification of about 1% is achieved for the µ's in the previously defined region and for the electrons inside a fiducial region which excludes about twice the size of the physical ECAL cracks. A study of selected data samples [11] has shown that the efficiencies are reproduced by the Monte Carlo at the two per mil level. The small resultant corrections are applied in the polarisation measurements.

Photon reconstruction
The reconstruction of photons from clusters of cells in the electromagnetic calorimeter is performed with the ALEPH standard algorithm [10] which searches for local energy maxima in the ECAL pad clusters. Fluctuations of a shower can generate "fake photons" [17] which are artefacts of the clustering algorithm or true photons produced by secondary interactions in the ECAL.
In order to reduce the number of such fake photons, an estimator of fakeness (P f ak ) is built from the relevant parameters [17]: the oriented distance to the closest charged track, the fraction of the shower energy deposited in the first ECAL stack, and the total ECAL shower energy. The sign of the distance is computed depending of the position of the shower with respect to the track bending in the r − φ projection. The reference distributions are taken from Monte Carlo, with the ECAL barrel and end-cap being handled independently. A cut on the fakeness likelihood estimator provides a partition of the photon candidate sample in two subsamples: the "fake photon" sample (P f ak >0.75), enriched in fake photons and the "genuine photon" sample (P f ak <0.75), enriched in genuine photons. The distributions of the variables used to construct the P f ak estimator are shown in Fig. 1 for each of the two subsamples of photon candidates. A good agreement between data and τ + τ − Monte Carlo can be observed.

Non-τ background rejection
Non-τ background in the selected samples has four main sources: Bhabha events, Z → µ + µ − decays, Z → qq decays, and γγ processes. As stated before, likelihood estimators exploiting the information of a single hemisphere are used for the suppression of the backgrounds that arise from the first three categories.
The Bhabha estimator (E e + e − ) provides the likelihood for one hemisphere to be part of a Bhabha event rather than of a τ + τ − . It is based on three variables: the e/π probability of the track, its momentum, and the ECAL wire energy. Specific reference distributions are used for tracks pointing to an ECAL crack.
The Z → µ + µ − estimator (E µ + µ − ) is similar to the previous one with the µ/π track probability replacing the e/π probability.
A third estimator (E qq ) tags a Z hadronic decay in one hemisphere. It makes use of the opening angle between the charged tracks, the number of charged tracks, the number of objects reconstructed from the pattern of the energy deposited in the calorimeters, and the mass and energy of the jet. It is very similar to the estimator used for the leptonic cross section measurement and is described in Refs. [11] and [18]. Different reference distributions are taken for each ECAL region (barrel, endcaps, and their overlap).
The γγ background is rejected by cuts on the three pertinent variables: visible energy, missing transverse momentum, and acollinearity of the two jets. The cut values are optimized for each final state.
Two less important backgrounds can also contribute: ℓ + ℓ − V and cosmic ray events. The ℓ + ℓ − V events are eliminated by the same cuts that reject Bhabha, Z → µ + µ − , and Z → qq events. Cosmic rays can contaminate the τ → hν and τ → µνν channels. They appear as events with two collinear tracks of the same momentum coming from the vertex region and a low energy deposit in the electromagnetic calorimeter. A first set of cuts is applied to select them: ECAL energy below 2.5 GeV, |p 1 | − |p 2 | < 5.0 GeV/c, min(|d 1 0 |, |d 2 0 |) > 0.1 cm, and |d 1 0 + d 2 0 | < 0.5 cm, where p i and d i 0 are respectively the momentum and signed impact parameter of the track i. The short gating time of the ITC is then used [19]. On the sample of cosmic ray candidates selected by the preceding cuts, the distribution of the number of ITC hits is measured. It is found that 68% of the events have no hits and 88% less than eight. Since the number of ITC hits is related only to the arrival time of the particles, the comparison of its distributions in the two samples, cosmic ray candidates and τ decays, is used first to estimate quantitatively the residual background in the τ samples. Finally, the events in the previously defined sample are rejected, along with events satisfying the same criteria on ECAL energy and momenta, and N IT C < 8.

Channel selection
The ALEPH standard leptonic pre-selection described in Ref. [11], which essentially requires a multiplicity of charged tracks between two and eight, is used first. The different channels are then selected in the angular acceptance (| cos θ| < 0.9) with the help of the charged particle and photon identification. The "π ± "(γ(s)) mass distribution in the muon sample before the hadronic veto. The shaded area shows the contribution of τ → µνν taken from the τ + τ − Monte Carlo. The black dots are the data, the solid (dashed) histogram the τ + τ − Monte Carlo prediction taking (not taking) into account the correction to the hadron misidentification rate described in Section 4.2.1.
The leptonic τ decay sample is selected from hemispheres containing only one charged track coming from the interaction region and identified as electron or µ. In the case of electron, it is required that the charged track does not point to an ECAL crack. To reduce the hadronic τ decay background possibly arising from particle misidentification, it is required that there is no genuine γ pair with an invariant mass compatible with the π 0 mass, and that the jet (hemisphere) mass is smaller than 0.5 GeV/c 2 . Since bremsstrahlung is only significant for electrons, muon hemispheres with one or more photons are also rejected when the hemisphere mass lies between 0.5 GeV/c 2 and 1.6 GeV/c 2 . This "hadronic veto" is illustrated by Fig. 2 which shows the γγ mass in electron hemispheres, and the jet mass in muon hemispheres containing at least one photon.
The non-τ background is rejected by cuts on the above described E e + e − and E µ + µ − estimators. To further suppress the Bhabha event contamination in the electron channel, the total ECAL wire energy is required to be smaller than 0.9 √ s. The ECAL wire energy distribution, before this last cut, is shown in Fig. 3. The agreement observed between the data and the sum of signal and background Monte Carlo shows the good understanding of the non-τ background contamination in the electron channel. The two-photon background is rejected by the cuts defined in the previous section.

τ → hν
For the selection of hadronic decays without a neutral hadron, the presence of only one charged track coming from the vertex region, with a momentum larger than 1.3 GeV/c (limit of the validity of µ-PID), is required. It is also required that this track is not identified as an electron or a muon by the PID. Two further cuts are imposed on the same hemisphere: no genuine reconstructed photon and a calorimetric energy deposit consistent with the momentum of the charged hadron. This last cut reduces the K 0 L π background.
The rejection of the hadronic Z decays and other non-τ backgrounds is performed by cuts on the values of the previously defined likelihood estimators computed on the opposite hemisphere. The two-photon background is rejected by the cuts described in Section 4.2.3.

τ → ρν
This selection is similar to the hν channel selection as far as the characteristics of the charged track and the non-τ background rejection are concerned. The requirements on the photons in the charged track hemisphere are opposite. The presence of one or two genuine photons is required. When a single photon is reconstructed, its energy must be greater than 0.5 GeV and a cut on the πγ mass (0.4 GeV/c 2 < m πγ < 1.6 GeV/c 2 ) is imposed. When two photons are reconstructed, no cut on the π2γ mass is applied but the γγ mass has to be compatible with a π 0 mass (80 MeV/c 2 < m γγ < 250 MeV/c 2 ).
The quality of the simulation is illustrated in Fig. 4, which shows the π 0 (one and two γ(s)) energy distribution in the selected ρ sample along with the single photon fraction, and in Fig. 5, where the jet mass in the one γ sample (before the above described cut on this mass) is displayed.
Here also, the selection criteria on the charged track and the non-τ background rejection are the same as in the hν case. The additional requirements on the photons in the selected hemisphere are as follows. The number of genuine reconstructed photons, n γ , ranges from two to four. For n γ = 2, the γγ mass has to be larger than the π 0 mass and inconsistent with it. For n γ = 3, the reconstruction of at least one π 0 is required and, for n γ = 4, the reconstruction of two π 0 's is required. The choice of the best association of the photons and the π 0 reconstruction are performed by a kinematic fit [10] which also improves the energy determination.
The purity of the selection and the quality of the reconstruction for the four photon sample are shown in Fig. 6 where the mass of a γγ system is displayed when the mass of the two remaining photons is consistent with the π 0 hypothesis.
The presence of three charged tracks coming from the interaction region is required. At least two of the three tracks must be identified as charged hadrons by the PID. The following cuts against the presence of π 0 's in the hemisphere are applied: no genuine photon with energy larger than 1 GeV, no genuine γ pair with E γγ > 3 GeV, and m γγ compatible with m π 0 . Hemispheres containing a V 0 , as defined by the ALEPH V 0 package [10], whose mass is compatible with a K 0 S are rejected. In addition, the effective mass of at least one of the π − π + combinations is required to be larger than 0.55 GeV/c 2 . The non-τ background rejection is the same as in the other hadronic decay channels.
The charged track selection and the non-τ background rejection are those used in the selection of the hν channel. Only the cuts against the two-photon event background are tighter, because the pion momentum distribution used for the polarisation measurement is very sensitive to this background. Since the purpose of the study of the inclusive channel is to cross-check the photon and π 0 treatment, no cut using information from the electromagnetic calorimeter in the selected hemisphere is applied.
A summary of the errors in the single-τ method is given in Table 2 and their estimation is discussed in the following.
The dominant contributions are photon/π 0 reconstruction, ECAL calibration, nonτ background, and modelling of the decay for the a 1 . This last effect is evaluated [2] by changing the model and its parameters within the limits allowed by the data.
The systematic effects are evaluated for each of the previously described tools using the difference between the data and Monte Carlo plus one standard deviation as input. The errors are propagated through all the steps of the analysis, taking into account the fact that a given tool can be used in several places. For that, A τ and A e are measured by means of the standardŴ ± (ω) functions on Monte Carlo events weighted for the change of the tool response. The shifts of the values of A τ and A e give the associated errors. Cross-checks of the estimations are presented for some of the important contributions.

Charged particle identification
The comparison of the PID performances on data and Monte Carlo is reported in Ref. [11]. The uncertainty on the data/Monte Carlo efficiency-ratio computed for each year of data taking is used to estimate systematic error.

Photon reconstruction
Clustering algorithm The first step of photon reconstruction is the "standard " ALEPH clustering algorithm [10]. Its parameters are the distance to the closest charged track and the threshold energy. It has been shown in Ref. [17] that the distance parameter can be changed by 1 mm, while preserving within one standard deviation the agreement between data and Monte Carlo on the distance distribution. This parameter shift introduced in the simulation induces negligible effects on the polarisation, even in channels containing photon(s). The possible effect of the threshold value is investigated by modifying the reconstructed shower energy in the Monte Carlo by the clustering correction [10]. The resulting variation of the polarisation measurement is small but not completely negligible. It is part of the value given in the line "photon" of Table 2.

Fake photon estimator
The second step is the genuine/fake classification for which the P f ak estimator is used. In the first place, the efficiency of the identification of genuine photons is measured on a selected sample made of photons forming a π 0 with one "high quality"photon (converted γ or isolated shower with 4 GeV< E < 20 GeV). For energies greater than 4 GeV, the efficiency is not significantly different from one. The efficiencies for E < 4 GeV are shown in Fig. 7. The ratio between data and Monte Carlo is (99.50±0.19)×10 −2 and an efficiency variation of 0.69×10 −2 is used as input in the error computation. Since the Monte Carlo is able to reproduce the characteristics of the fake photon production in a hadronic shower but not its precise normalization [17], a weighting procedure is used to correct it. The weight for an event is written w(p) n , where n is the fake photon multiplicity and p the hadron momentum. The function w(p) is obtained by a fit of the likelihood estimator distributions. A check of the procedure is given by the comparison of data and Monte Carlo for the γγ mass spectra shown in Fig. 8. The related systematic error is estimated in the standard way described at the beginning of Section 4.4, comparing the measured weight to one.

Bhabha, µ + µ − , and qq event backgrounds
The E e + e − , E µ + µ − , and E qq estimators use information from only one of the hemispheres of an event, so correlations with the other hemisphere can only be due to geometrical correlations in the detector construction (ECAL cracks, overlap, etc.). The measurement of the estimator efficiencies in the non-correlated regions is described first, and then an evaluation of the residual background, which is not affected by the possible correlations, is presented. This last point concerns only the hadronic channels because the selection criteria for leptonic τ decays eliminate the correlated, insensitive parts of the detector.

Efficiency of the estimators
In order to measure the efficiencies, almost pure samples of data and Monte Carlo events are selected for the τ + τ − , e + e − , µ + µ − , and qq final states. A clean τ decay (three charged π's or a ρ), a high energy electron, a high energy muon or a hadronic jet is required in one hemisphere and the opposite one is used to study the estimators. Examples of their responses in the case of E e + e − and E µ + µ − are shown in Fig. 9 and 10.
A cut at a value of 0.8 on these two estimators allows a good separation between τ 's and leptonic background. The measured efficiencies, taking into account the small residual contamination of the samples, are given in Table 3. These numbers are used to determine the systematic uncertainties which contribute to the line "selection" of   Figure 10: (a) The E µ + µ − distribution after tagging a τ in the opposite hemisphere.
(b) Efficiency on τ + τ − events as function of the cut value. Black dots are the data, empty boxes and histogram, the τ + τ − Monte Carlo. Table 2 and to the line "non-τ background" for the two leptonic channels.
A similar analysis, already described in Ref. [11], is performed for the E qq estimator.

Residual backgrounds
The estimation of the residual background, chiefly Bhabha events, in the hadronic channels is an essential task on account of its impact on the A e measurement. For that reason, procedures were developed to measure directly the background contributions from the data themselves.
The method is exemplified by the case of the τ → hν channel. Applying all the cuts used for the channel selection but the e/π identification cut on the τ side, and the E e + e − cut on the opposite one, a set of data is constituted which is divided in two subsets by the response of the e/π PID: (a) for e PID and (b) for π PID. The distributions of E e + e − in the two subsets and the τ + τ − Monte Carlo expectations shown in Fig. 11 clearly display the Bhabha event background. The events with E e + e − > 0.8 in sample (b) are used to determine the kinematic characteristics of this background. The normalization of the residual background after the E e + e − cut is given by the number of events with E e + e − > 0.8 in sample (b) multiplied by the ratio of the numbers of events with E e + e − < 0.8 and E e + e − > 0.8 in sample (a). This background is subtracted for the polarisation measurement and the statistical error on its estimation enters the systematic error. The small difference of the efficiencies of the cut on E e + e − in samples (a) and (b), mainly due to geometrical effects and initial state radiation, can be computed by Monte Carlo and is added to the systematic uncertainty.
As a cross-check, an independent study of the Bhabha event background has been performed, using the tools developed for the measurement of the τ leptonic branching ratios [19]. Its conclusions are consistent with the above presented results.
Similar procedures are used to study the µ + µ − contamination. They show that such a background is reduced to a negligible level by the selection process.

Two-photon background
In a similar way, the residual γγ background can be studied in the data from the excesses of events in the distributions of the variables (acollinearity, visible energy, etc.) used for its rejection.
Such a background appears mostly in the leptonic channels. The contaminations evaluated by a study of the opposite hemispheres are (0.2 ± 3.2) × 10 −4 for the electron channel and (5.5 ± 3.6) × 10 −4 for the muon. To compute the tiny associated uncertainties the FERMISV Monte Carlo [16] is used.

Cosmic ray background
The contaminations of the two relevant channels, τ → hν and τ → µνν are estimated using the independent information from the track parameters and the number of ITC hits. Their values are (6.9±1.9)×10 −4 and (1.5±0.5)×10 −4 respectively, leading to very small contributions to the uncertainties. Out-of-time cosmic rays are normally identified as hadrons, owing to the HCAL and muon chamber inefficiency in such events [19].

Momentum calibration
The uncertainty on the momentum calibration described in Section 3.1 is used, for each cos θ bin, to estimate the contribution to the polarisation systematic errors.

ECAL calibration
The ECAL energy calibration (Section 3.1) is obtained from electrons taken in Bhabha events, γγ processes, or τ decays to electrons. Uncertainties coming from the statistics of these samples and the handling of the ECAL saturation are introduced as systematic errors, leading to an uncertainty of 2.5× 10 −3 on the calibration, which is used as input for the computation of polarisation errors.
An independent check of the ECAL calibration has been performed using the µ + µ − γ events for which the photon energy is deduced, with a small error, from the angles, by the formula γ between the kinematically calculated and ECAL energies, and σ Eγ the expected resolution, are shown in Fig. 12 for data and Monte Carlo and a very good agreement can be observed.
Since some calorimetric cuts are used in the selection of the two hadronic final states without photons, a tuning of the ECAL response to charged hadrons is also performed. The related uncertainties are independently evaluated.  Table 4 gives the numbers of selected data events, the selection efficiencies in the angular acceptance, and the contaminations from τ and non-τ backgrounds.

Results
The polarisation is measured independently for each channel and for each year of data taking. Some of the systematic errors are common to all the years and some are common to different channels.
In a first step, the measurements for each channel in the different LEP operation periods are combined, taking care of the common systematic errors and adding the small corrections dependent on the centre-of-mass energy described in Section 6.1. The values obtained are given in the corresponding lines of Table 5. Figure 13 shows the x distributions for the hadron, muon, and electron channels; and the ω distributions for the ρν and a 1 ν channels (3h and h2π 0 ). The distributions for data and the contributions (Monte Carlo) from the two τ helicities and non-τ background are superimposed.
The measurements from all the channels are then combined taking into account their correlations evaluated by Monte Carlo. The inclusive hadronic analysis has been devised as a cross-check of the exclusive ones since their systematic uncertainties are essentially uncorrelated. The values given in Table 5 show that this test is satisfactory. Statistically, on the contrary, the exclusive and inclusive analyses are Distributions of the ω variable for τ → ρν (d) and τ → a 1 ν: a 1 → 3h (e) and a 1 → h2π 0 (f). The dotted and dash-dotted lines correspond to the contributions of left-and right-handed τ 's respectively. The shaded area is the non-τ background contribution.  The final values of A e and A τ are given in the last line of Table 5. Their errors are corrected for the correlation between opposite hemispheres.
As an illustration of the understanding of the ρ channel, Fig. 14 shows the τ and ρ decay angle distributions with good agreement between data and Monte Carlo, even for low energy ρ's (cos ψ τ < −0.9).
5 Analysis with the τ direction method 5

.1 Philosophy of the method
The τ polarisation analysis described in this section is based on two ideas: use of proven analysis tools and enhancement of the measurement sensitivity. On one hand, event selection, particle identification, photon and π 0 reconstruction, and decay classification are taken from previous ALEPH analyses of τ leptonic and hadronic branching ratios [17,19], and of the hadronic spectral functions [20,21]. On the other hand, the maximum sensitivity to the τ polarisation can be achieved in the hadronic decay channels through a complete set of observables including the determination of the initial τ direction. This new information, as compared to standard analyses, is obtained through a kinematic reconstruction complemented by the precision measurement of the charged particles using the microvertex detector in order to remove part of the intrinsic ambiguity of the method.

τ -pair selection
This analysis is based on a global event selection that retains τ pair candidates from Z decays. Each event is divided into two hemispheres along the thrust axis. The τ event selection is described in [17] and references therein. The method is based on the removal of non-τ background identified as Bhabha, mu-pair, hadronic, photon-photon induced, and cosmic ray events. A well-tested procedure of reconstructing the spatial energy flow [10] is used in order to include all the energy measurements and their topology in the detector.
However, to enhance the selection efficiency and to reduce the biases in the polarisation measurement, some minor modifications of the usual cuts [17] are introduced: • cuts against γγ events and the cut on the energy of the leading tracks in both hemispheres against Bhabha and di-muon events are not applied when both hemispheres are identified to be hadronic.
• in order to reduce the contamination from γγ → τ τ and γγ → qq background, events with one hemisphere with a single charged hadron and no π 0 are removed if | cos θ| > 0.792 and x < 0.088. • hadronic cuts are loosened taking advantage of the more precise hadronic mass determination after π 0 reconstruction, rather than relying on the energy-flow variables used in the selection.

Particle identification
A likelihood method for charged particle identification incorporates the relevant information from the detector. In this way, each charged particle is assigned a set of probabilities from which a particle type is chosen. A detailed description of the method can be found in Refs. [19] and [22].
Briefly, the particle identification is based on the following information: (i) the dE/dx in the TPC, (ii) the longitudinal and transverse shower profile in the ECAL near the extrapolated track, and (iii) the energy and the average shower width in the HCAL, together with the number of planes, out of the last ten of the HCAL, which had fired and the number of hits in the muon chambers. Because of the calorimetric requirements, the complete identification method is only applied to particles with momenta larger than 2 GeV/c.
The performance of the identification method has been studied in detail using samples of leptons from Bhabha, µ-pair, and two-photon events, and hadrons from π 0 -tagged hadronic τ decays. All data sets cover the full angular and momentum range [17,19]. Ratios of particle identification efficiencies in the data and the corresponding Monte Carlo samples are obtained as a function of momentum. They are then used to correct the Monte Carlo efficiencies in the τ + τ − sample, thus taking into account the different momentum distributions and particle environment. Misidentification probabilities are obtained in the same way. The results of these analyses were already given for the 1990-1993 data [19]. Figures 15 and 16 show the corresponding efficiencies and misidentification probabilities measured in the 1994-1995 data. Table 6 gives the efficiency matrix for one-prong τ decays used for this analysis, i.e., for particles with a momentum larger than 2 GeV/c and not in a crack region between ECAL modules. The 2 GeV/c cut is not applied to electron candidates because of the good dE/dx separation at low momentum, while the crack veto is not required for muons which are essentially identified in the HCAL. The numbers refer to the 1991-1995 data sample.  . The values are obtained directly from data using a π 0 -tagged hadron sample.

Non-τ background
A new method is developed to directly measure in the final data samples the contributions of the major non-τ backgrounds: Bhabhas, µ + µ − pairs, and γγ → e + e − , µ + µ − events. The procedure does not rely on an absolute normalization for the simulated channels, which is difficult to obtain reliably after the large necessary rejection is applied. It only makes use of a qualitative description of the distribution of the discriminating variables. The basic idea is to apply cuts on the data in order to reduce as much as possible the τ τ population while keeping a high efficiency for the background source under study. The background fraction is directly fitted in the cut data sample and then extrapolated to the standard sample.
For Bhabha background, different event topologies are considered in turn. For e-e and e-crack topologies, the acoplanarity angle (defined by the two tracks azimuthal difference) is required to be larger than 179 • or the acollinearity angle (i. e. 180 • minus the angle between the directions of the two tracks) should be less than 175 • with the difference of transverse energies being less than 3 GeV. The corresponding Bhabha efficiency (compared to the normal sample) is 88% as determined from the simulation. The angular distribution of the restricted sample is then fitted to τ τ and Bhabha components (also including a small contribution from the other non-τ backgrounds) from the simulation. Therefore, the Bhabha Monte Carlo input is only used to determine the (large) selection efficiency and the cos θ shape inside the final sample, not relying on any determination of the absolute Monte Carlo normalization. The derived Bhabha contribution has a statistical uncertainty which is assigned as a systematic error. Several combinations of variables have been tried, showing a good stability of the result within its error. For the more numerous e-h and h-h samples, a total energy cut of 55 GeV is applied and the main reduction of τ events is achieved by suppressing true hadrons as compared to electrons misidentified as hadrons: this is done by restricting the opposite hadron in an e-h event to have an electron identification probability larger than 0.01 (most of the true hadrons are below this value). A similar approach is used for the h-h samples with even tighter identification probability cuts for hemispheres with only one charged hadron (in this case 86% of the hadrons have an electron probability smaller than 0.0002). The Bhabha signal then appears on a distribution of the acollinearity angle between the two hemispheres (δθ) at values close to 180 • .
A similar technique is used to estimate the µ-pair background. An efficiency of 90% is obtained by requiring an acoplanarity angle larger than 178 • for µµ events in the final sample. Here, the fitted distribution is that of the calculated photon energy along the beam for a postulated µµγγ beam kinematics to take the most general case compatible with only the information on the two muons. For µ-h and h-h final states, the procedure follows the one used for Bhabha background.
The measurement of the remaining background induced by γγ interactions is based on cuts on the total event energy (typically E < 35 GeV, depending on the channels) and the acollinearity angle between the two hemispheres (δθ < 175 • ). The γγ background shows up clearly at small overall transverse momentum well above the small τ surviving contribution.
All the above non-τ background contributions are determined on the polarisation 27 data sample with the help of special cuts. The corresponding Monte Carlo generators are only used to compute the efficiency of the cuts. This method is illustrated in Fig. 17 showing the evidence for Bhabha and γγ → τ τ or qq backgrounds in the ρν channel, the contribution from µµ being negligible in this case. As an important bonus, this procedure leads to a determination from data alone of the background distributions which are relevant for the polarisation measurement. This method is preferable to using the Monte Carlo predictions which are questionable in view of the large rejection factors achieved. Distributions of the ω variable obtained in the ρν channel are given in Fig. 18.
As in the analysis of Section 4, cosmic ray background is determined to be very small after the additional cuts [19]. The contamination from hadronic Z decays affects only the hadronic channels. Its small contribution is estimated with the JETSET 7.4 generator [24]. Tests were previously made to ascertain the reliability of the prediction for the rate of surviving low multiplicity events by comparing to data [17]. A 30% uncertainty was derived from these tests. The non-τ background contributions are given in Table 7, separately for the considered channels.

Photon and π 0 reconstruction
As in the previous analysis, the photon and π 0 reconstruction is performed according to a likelihood method in which several estimators and probabilities are computed in order to distinguish between genuine and fake photons produced by hadronic interactions in ECAL or by electromagnetic shower fluctuations [17].
An algorithm is used to form photon pairs in each hemisphere [17] using a π 0 estimator D π 0 i,j = P γ i P γ j P π 0 i,j for two photons i and j, where P π 0 i,j is the probability from a kinematic π 0 -mass constrained fit. A good agreement between data and Monte Carlo is observed for the distributions of P π 0 i,j and D π 0 i,j . Figure 19 shows the corresponding π 0 mass distributions with the expected contribution from photon pairs including at least one fake photon.
High-energy π 0 's with overlapping photon showers are reconstructed through an analysis of the spatial energy deposition in the ECAL towers ("unresolved π 0 's"). All the remaining single ("residual") photons are then considered and those with 29 P γ > 0.5 are selected as π 0 candidates. The P γ distributions are given in Fig. 20 for all photons and for the residual ones after renormalizing the amount of fake photons in the simulation by a factor depending on the τ decay hadronic final state and ranging from 1.17 to 1.28. The shapes in data and Monte Carlo are in good agreement. The energy dependence of the data distributions is also well reproduced by the simulation: hence the photon efficiencies resulting from the D π 0 ij and P γ cuts are adequately described by the Monte Carlo, as already observed in Section 4.4.2 and Fig. 7.
Finally, photon conversions are identified following the procedure described in Ref. [17]. They are added to the list of good photons and included in the π 0 reconstruction. The fractions of resolved, unresolved and residual π 0 's are presented in Fig. 21 as a function of the π 0 energy. It is seen that the data distributions are well reproduced by the simulation. The fraction of resolved "calorimetric" π 0 's remains significant out to large energies in this analysis as compared to the previous one (see Section 4.3.3 and Fig. 4); this effect is the result of a looser π 0 definition here, where highly asymmetric γγ pairs are more likely to be counted as resolved π 0 's, whereas they would be identified as unresolved in the other method because of the tight invariant mass cut. Figure 19: Photon-pair mass distribution for reconstructed π 0 candidates. Data (points) are compared to the simulation (histogram). The shaded histogram gives the simulated contribution where at least one of the photons is fake.

Decay classification
τ decays are classified according to the number of charged tracks and their identification, and the number of reconstructed π 0 's. Beyond the classification achieved in the previous analyses [17,[19][20][21], some additional cuts have been developed in order to reduce the ρ and K * background level in the h channel, as explained below.
In order to reject ρ decays where the π 0 is close to the charged pion and undetected as such, the sum of the energies in the first two stacks of ECAL is required to be less than 75% of the charged pion momentum. The energy is measured including all the pads in a cone of 30 • around the track. To further reduce the ρ contamination, the hemispheres with one charged hadron and a single residual photon with a probability P γ greater than 0.2 are rejected.
The K * (decaying into K 0 L π) component is reduced by a cut on the HCAL energy deposit (E hcal ) and the azimuthal offset (δφ) between the track impact and the energy weighted barycentre of the energy deposit in HCAL [25]. The cut is applied on the following variables: where P h is the hadron momentum and σ h the corresponding HCAL energy resolution. The value of the parameter ξ is +1 if φ barycentre is in the direction of the charged hadron bending, and −1 otherwise. The following region in the (δ E , δφ) plane is cut : In order to verify that those supplementary cuts do not introduce biases into the measurement of the τ polarisation, the ratio of the efficiencies for each of the three cuts on data and on Monte Carlo, ǫ data /ǫ MC , is studied as a function of the pion momentum (Fig. 22). The computation of these efficiencies and of their statistical errors takes into account the non-pion background estimated by Monte Carlo. The small bias introduced by these cuts is evaluated by a linear fit and corrected in the final polarisation measurement. The corresponding statistical error is taken as a systematic uncertainty.
For the ρ channel, a cut on the invariant mass is required, m ρ > 0.4 GeV/c 2 , to reduce the feedthrough from the h mode. For the a 1 channel, the invariant mass is required to be less than the τ mass. In the case of a decay into three charged hadrons (3h), three good tracks are required in the hemisphere.
The channel identification efficiencies, the τ background, and the non-τ background contributions as measured in the data (except for qq) are shown in Table 7. For events where both τ 's decay into hadrons (called hadronic events), it is possible to reconstruct the τ + τ − direction, thereby enhancing the sensitivity of the polarisation measurement in the ρ and a 1 channels as seen in Section 2. This method is used here for the first time.
In hadronic events, although the momentum of each neutrino is unknown, the τ direction can be determined up to a twofold ambiguity. The two solutions lie at the intersection of the cones having for axis the hadronic direction in each hemisphere and an opening angle φ 0 computed from the measured hadron momentum P h and energy E h , assuming the pion mass for each individual hadron, where P τ and E τ are the energy and the momentum of the parent τ assumed to be produced at the beam energy, and m h the hadronic mass.
A procedure is used for events where the cones do not intersect, mainly because of detector resolution effects. The hadron momenta are allowed to fluctuate within the expected resolutions of charged and neutral pions. The τ direction is taken as the average of the found solutions among the random trials, provided at least 10 solutions are found for 500 trials. This procedure attributes a τ direction to half of the events where the two cones are not initially intersecting.
As a result, 80% of the hadronic events are available for the polarisation analysis using the τ direction. The remaining 20% without a reconstructed direction are analysed in the standard way. On the whole 52% of the τ → ρν decays benefit from this improvement.

Handling the two τ directions
Since the typical τ decay length is 2 mm at LEP I, the precise determination of the secondary charged tracks brings some information allowing in principle to lift the twofold ambiguity and to choose the actual τ direction. This procedure stems from the 3-dimensional method developed for the measurement of the τ lifetime [6]. In practice one can only separate the two solutions on a statistical basis with the help of an estimator, D h , and assign to each of the two directions a probability P 1,2 to be the true one. This procedure is explained in Appendix A, while more details can be found in Ref. [26].
Reference distributions for the estimator are set up with Monte Carlo for various event configurations (π-π, ρ-ρ, ρ-π, etc.) and probabilities P 1,2 for each of the two found τ directions to be the true one are obtained.
If the highest probability is retained, it is found by simulation that the closest of the two directions to the true one is chosen in 65% of the cases for all channels combined. A more relevant way to characterize the procedure is to consider the precision achieved on the τ direction. While the mean value of the angle between the two directions is 24 mrad, the angle between the direction with the higher probability and the true one is 12 mrad on average. The corresponding value for the direction closest to the true one is 7 mrad (dominated by resolution effects) and a value of 15 mrad is obtained if a random choice between the two solutions is made. In practice, the two solutions are retained, with their respective probabilities.

Analysing the polarisation with the τ direction
The polarisation is analysed in each hemisphere using the proper optimal variables, ω 1,2 , calculated for the observed decay for each choice of direction using the description of the hadronic system given in Ref. [8]. Both directions are entered into the expected decay distribution depending on P τ : with the new optimal observable The ideal sensitivity given in Table 1 is naturally degraded by detector resolution effects and mostly by the imperfect determination of the τ direction. The expected sensitivities are given in Table 8 from a Monte Carlo study. The gain in sensitivity using the calculated probabilities instead of affecting each direction with an a priori 50% weight is quite small, but it provides the ultimate gain in precision achievable within the ALEPH detector capabilities.
For events where one of the τ 's decays leptonically and for which the information on the τ direction cannot be retrieved, the acollinearity angle between the decay products in the two hemispheres contains additional information as compared to the separate analysis of single decays. This method [2,23] is also used in this analysis.

Summary of systematic effects
Possible sources of systematic effects come from the τ selection, particle identification, and π 0 reconstruction. The dominant systematic uncertainties are due to photon  identification and π 0 reconstruction, the non-τ background, and, in the case of the a 1 modes, the dynamics of the decay. Table 9 shows the summary of the main components of the systematic uncertainties for A τ and A e .
The particle identification efficiency matrix is measured as a function of momentum using data samples of kinematically identified particles, as described in Section 5.2.2. The procedure yields both the relevant efficiencies and their associated uncertainties which are dominated by the statistics of the control samples. The effect of this systematic uncertainty on the polarisation measurement is significant only for the e or µ channels.
To study the effect of the τ pair event selection, the difference in the selection efficiency between data and the Monte Carlo is analysed for each cut, as discussed in detail in Refs. [17,19]. The systematic error on the selection has a significant effect only on the h channel where the nearly flat momentum distribution is distorted at low values by the γγ cuts and for large momenta by the Bhabha and µ-pair cuts.
The branching ratios used in the Monte Carlo are corrected using the measured values of ALEPH [17,19] and the uncertainty in the amount of feedthrough in a given τ decay channel from the other modes is obtained by varying the branching ratios within the errors given in Refs. [17,19], taking into account the correlations between them. The impact on the τ polarisation measurement is negligible.

36
The tracking resolution affects only the a 1 decay since three reconstructed charged hadrons are required. In this decay the tracks are generally close to each other and therefore partially overlap. A reconstruction problem could occur if the two same-sign tracks have nearly the same transverse momentum. Data and Monte Carlo distributions of the angle between the two tracks at the vertex are in fact in agreement within statistics and a systematic uncertainty for the polarisation is derived. The systematic error on the momentum calibration is also estimated with the largest effect observed on the µ channel.
A complete investigation of the observables entering photon and π 0 reconstruction is undertaken. For each variable a comparison between data and Monte Carlo is performed in order to search for possible systematic deviations. These differences are then injected in the complete reconstruction procedure, affecting both the number of events in a given channel and the corresponding ω distributions. Each data-Monte Carlo comparison yields a systematic uncertainty limited by the statistics of the data. As discussed in full detail in Ref. [17], these studies include the effect of the clusterization algorithm, the photon probabilities, the threshold energy cut for low energy photons or the minimal distance between a charged track and a photon in the electromagnetic calorimeter. Similarly, the effect of the cuts on the π 0 estimators is also determined. All the contributions to the total systematic uncertainties for this source are added in quadrature.
Another source of uncertainty in the photon and π 0 treatment is due to the problem of fake photons generated by hadron interactions in ECAL or electromagnetic shower fluctuations. This source of photon candidates is underestimated in the Monte Carlo simulation with respect to the data. This deficit produces a bias on the polarisation measurement and needs to be corrected for with the systematic uncertainty carefully estimated. In general, the effect has been studied by removing fake photons in the simulation and repeating the analysis: decay classification and polarisation determination in all channels. This method establishes the sensitivity of the measurement to the problem. Then the fake photon deficit in the simulation is measured in all the needed hadronic channels by fitting probability distributions in data and Monte Carlo as discussed in Section 5.2.4 and the corresponding systematic uncertainties are deduced.
A more direct method is used in the h channel for the events with a residual photon where a tighter probability cut was applied to reject more ρ background (see Section 5.2.5). Since the ρ contribution can be well estimated by Monte Carlo, it is possible to use the results of this cut both on data and Monte Carlo to estimate the relative effect of fake photons as a function of momentum. Any momentum dependence of the Monte Carlo deficit would produce a bias for the polarisation. Such a bias is measured (Fig. 22) and the statistical uncertainty in the linear momentum fit is taken as systematic error. This procedure in the h channel does not take into account the loss of events with multiple fake photons entering the definition of π 0 candidates. A special study was made to estimate this effect as the simulation tends to underestimate the fake photon multiplicity, with the result that the effect on the polarisation is small compared to that of residual photons.
Since some additional calorimetric cuts are introduced in order to further reduce the level of the ρ and K * background in the hadron channel, it is important to study their effect on the hadron momentum distribution. The results have been presented in Section 5.2.5 and the corresponding systematic uncertainties derived.
Finally, the effect on the A τ and A e measurements from the different sources of nonτ background has been studied in detail directly on the data with the method presented in Section 5.2.3. This procedure allows the energy distribution of the remaining Bhabha events in the selected sample to be investigated and a direct measurement of their cos θ distribution to be made. The latter information is particularly relevant for the A e determination because of the asymmetric Bhabha contamination. Using as input the distributions derived from data with their statistical and systematic uncertainties, their impact on A e and A τ is assessed. It is important to note that the dominant statistical uncertainty affects both the normalization of the background and the shape of the ω and cos θ distributions. The most seriously affected channel is clearly that of the electron, both for A τ and A e . Table 10 gives the results for A e and A τ in the seven channels considered in this analysis. The statistical and systematic uncertainties are given separately. The last line presents the combination of these results taking into account their correlations, including the polarisation correlation between the two hemispheres in an event. Figures 23,24, and 25 show the ω distribution for the ρ and a 1 decay modes. The ω distributions for the decays where it was possible to reconstruct the τ direction are also shown separately.

Standard model corrections
The A τ and A e parameters are obtained by fitting the integral, in each cos θ bin, of the function (1) to the measured polarisation. To obtain the combinations of effective couplings (A l ) related to sin 2 θ eff W , corrections have to be applied which take care of the fact that the polarisation was not measured at the Z pole and that there are radiative corrections and Z-γ interference effects. These corrections are computed by the ZFITTER program [4] in the following way.
For a set of the standard model parameters (Z mass, top mass, Higgs mass, α s ), ZFITTER provides the value of the effective couplings in the framework of the standard model. Cuts are defined on acollinearity and momentum which reproduce the ones used in the analysis. Then the A f b asymmetry, and the polarisation in each cos θ bin are calculated at the given energy. The same fitting procedure as for the data is applied and the difference between the effective couplings and the fitted parameters is taken as the correction.
The study of the stability of the corrections against the variation of the input parameters within their errors shows that the related error is negligible compared to the systematic errors involved in the analyses. The same is true for the values chosen for the cuts. With this procedure, the inaccuracy of (1) in describing the polarisation does not introduce any bias but may simply reduce very slightly the sensitivity. The level of these corrections is 0.04% with an uncertainty of 0.01%. It has been checked that using as input the measured values of A f b gives the same results.

Consistency and combination of the results
To investigate the consistency of the two analyses, both statistical and systematic effects have to be taken into account.
As a result of the different selection procedures, there is a significant number of events which belong to only one of the two samples. The fractions of unshared events in the different channels, averaged over the full angular distribution, are 24% (π), 27% (ρ), 24% (a 1 → 3π), 47% (a 1 → π2π 0 ), 24% (e), and 13% (µ). Typically, these fractions are by a factor 1.6 larger in the small-angle region (| cos θ| > 0.7) than in the central part. The distribution of the events has to be taken into account in the evaluation of the consistency and in the combination of the measurements since this small-angle region has a significant contribution to the A e determination. Another statistical effect comes from the channels where the τ direction is taken into account since different information is used in the two analyses. Finally, the Monte Carlo samples do not completely overlap.
The systematic effects of the two methods have in common the detector behaviour and the basic event reconstruction; however the analysis tools used are largely independent. Furthermore, the procedures developed to measure efficiencies, and various corrections are different in the two analyses and thus a part of the estimated systematic uncertainties is uncorrelated.
The A τ values given by the two analyses -for each channel and globally -are in agreement within the estimated uncertainties. The combination of all channels, gives ∆A τ = (0.14 ± 0.52)%, where the statistical and systematic part of the error are comparable. The agreement is not as good in the case of the A e measurement where the deviation between the two analyses, mainly due to the hadron channel, is ∆A e = (0.92 ± 0.48)%. However, the systematic uncertainties on A e , dominated by the effect of the non-τ background (chiefly Bhabha events), are smaller by a factor of 5-10 than the statistical errors, and the many cross-checks performed show that the background levels and distributions are known within the quoted uncertainties. This lead to the conclusion that the observed difference is the result of a statistical fluctuation.

Final results
The consistency of the two sets of results having been checked, they are combined to obtain the final results. Since the two analyses have comparable statistical and systematic uncertainties, the combination procedure is not a critical one. It was chosen to average the results using equal weights. For the computation of the final errors the statistical correlation between the two samples is taken into account together with the systematic uncertainty correlations.  The results are given in Table 11 for all the analysed modes and in Table 12 for the global values. Figure 26 shows the average tau polarisation as a function of cos θ.

Discussion
The A τ and A e parameters are measured by means of similar methods by the three other LEP collaborations [27][28][29]. The A e asymmetry is also measured by the SLD collaboration [30] using the beam polarisation.     Comparing these numbers with the values of the couplings presented in Ref. [11] shows the important weight of the τ polarisation measurement in the Z leptonic coupling determination.

Conclusion
From the A τ and A e values obtained through the τ polarisation measurement: Using both the τ polarisation and the leptonic forward-backward asymmetries measured by ALEPH, the value sin 2 θ eff W = 0.23130 ± 0.00048 (14) is obtained. consequence of the τ mass constraint, where E, p, and m are energy, momentum, and mass of τ 's and hadrons (π's). Accordingly, the two kinematically found τ − directions, defined by the unit vectorsτ 1 andτ 2 , are the intersections of two cones having −ĥ + andĥ − for axes and φ ± for opening angles. The vectorτ 1 +τ 2 lies then in the hadronic plane and is, therefore, orthogonal to D h (Fig. 30). This implies The vector A h , projection ofτ on the hadronic plane, is the internal bisector ofτ 1 and τ 2 . These properties, with Eq. 16, lead to since l is, by definition, a positive length. The positivity of l and Equations 16 and 17 imply that the physical τ direction is determined by the relation For hemispheres with a 3-prong decay, the vectorĥ ± is defined as equal to the sum of the three momenta and originating from the 3-track common vertex.

A.2 The general case
When π 0 's are produced, the vector A h , which is kinematically determined, is still measurable, but the exact reconstruction of D h given in the previous section is not possible and an approximation must be made. 49 The vectorsĥ ± are now unit vectors along the sums of the hadron momenta in each hemisphere. Two vectors, which make with A h an orthogonal system, can be unambiguously measured: which is normal to the hadronic plane, and N =ĥ ⊥ ×τ =ĥ + cos φ − +ĥ − cos φ + , which is normal to the ambiguity plane spanned byτ 1 andτ 2 . In terms of h ⊥ and N, the analytical expressions of A h and D h are

The vector
−→ ∆H is equal to −−→ H + H − , where H ± are arbitrary points on the lines drawn along the vectorsĥ ± from the τ ± decay points, lines which coincide with the reconstructed tracks in the ideal case of the previous section.
The charged tracks reconstructed in each hemisphere are used to define an approximation of D h . Letĉ ± be unit vectors along the momenta of the charged particles, C ± points on the charged tracks in the neighbourhood of the decay points, c ⊥ =ĉ − ×ĉ + , and −→ ∆C = −−→ C + C − . If the vectorsĉ ⊥ andĥ ⊥ are collinear, the space planes parallel to the hadronic plane and containing the τ ± decay points contain also the charged tracks. The vector D h is then equal to the minimum approach vector between the charged tracks. In the most general case, an equation similar to Eq. 16 can be written [6]: