A Measurement of the Tau Topological Branching Ratios

Using data collected in the DELPHI detector at LEP-1, measurements of the inclusive (cid:28) branching ratios for decay modes containing one, three, or (cid:12)ve charged particles have been performed, giving the following results:


Introduction
Historically, there has been an inconsistency in the measured τ exclusive branching ratios: the measurements of exclusive branching ratios to decay modes containing one charged particle did not sum up to the inclusive branching ratio for a charged multiplicity of one.A number of analyses have attempted and succeeded in resolving this question [1][2][3][4].However the current world average values from direct measurements of the topological branching ratios in the Particle Data Group listings [5] have uncertainties which are significantly larger than the values obtained through combined fits to all the τ decay data.In the case of the 1-prong branching ratio (Γ 2 in [5]), the difference between the average and the fit value is more than two standard deviations.
This paper presents a dedicated simultaneous measurement of the decay rates of the τ lepton to different final states as a function of the charged particle multiplicity.
The relevant components of the DELPHI detector and the data-set are described in Section 2. The definition of the measured branching ratios and the method used to derive them are introduced in Section 3. In Section 4 the preselection of the sample of e + e − → τ + τ − events is described.Section 5 contains a description of the reconstruction of charged particles, as well as the reconstruction algorithms for photons and hadrons which interact with the detector material before the tracking subdetectors.The τ decays are classified according to their charged particle multiplicity.This is described in Section 5.4.The results and systematic studies are presented in Sections 6.1 and 6.2 respectively.

Experimental apparatus and data sample
The DELPHI detector and its performance are described in detail elsewhere [6,7].The subdetector units particularly relevant for this analysis are summarised here.These detector components covered the full solid angle considered in the analysis except where specified, and sat in a 1.2 Tesla solenoidal magnetic field parallel to the z-axis 1 .
The reconstruction of charged particles in the barrel region of DELPHI used a combination of the measurements in four different cylindrical subdetectors: a silicon Vertex Detector (VD); the Inner Detector (ID), consisting of a jet chamber tracker and a wire chamber used for trigger purposes; the Time Projection Chamber (TPC); the Outer Detector (OD).
The VD had three layers of silicon micro-strip modules, at radii of 6.3, 9.0 and 11.0 cm from the beam axis.The space point precision was about 8 µm in r-φ.For data from 1993 onwards, a measurement of r-z was produced in the outermost and innermost layers.This had a precision of about 15 µm.The-two track resolution was 100 µm in r-φ and 200 µm in r-z.
The ID had an inner radius of 12 cm and an outer radius of 28 cm.The inner jet chamber part lay at radii below 22 cm.It had a two-track resolution in r-φ of 1 mm and a precision in r-φ of 50 µm.The outer triggering part was a five-layer wire chamber.This was replaced for the 1995 data with a straw detector containing less material.
The TPC, extending from 30 cm to 122 cm in radius, was the main detector for charged particle reconstruction.The main track reconstruction information was provided by 16 concentric circles of pads which supplied up to 16 three-dimensional space points on a track.In addition, ionisation information was extracted from up to 192 wires.This was used for particle identification purposes.Every 60 • in φ there was a boundary region between read-out sectors about 1 • wide which had no instrumentation.At cosθ = 0 there was a cathode plane which caused a reduced tracking efficiency in the polar angle range | cosθ| < 0.035.The TPC had a two-track resolution of about 3 mm in r-φ and 1.5 cm in z.
The OD, with 5 layers of drift cells at a radius of 2 m from the beam axis, was important for the momentum determination of energetic particles.
The principal device for electron and photon identification was the High density Projection Chamber (HPC), located outside the OD, which allowed full reconstruction of the longitudinal and transverse shower components.It covered the polar angle region | cosθ| < 0.75.In the forward region, 0.800 < | cosθ| < 0.985, the calorimetry was performed by a lead-glass array.For very small polar angles, the calorimetry was performed by the luminosity monitors, the SAT in 1992 and 1993, and the STIC in 1994 and 1995.
Outside the magnet solenoid lay the hadron calorimeter and muon chambers, which were used for hadron and muon identification.
The data were collected by the DELPHI detector from the LEP electron-positron collider, in the years 1992 through 1995, at centre-of-mass energies √ s of the e + e − system between 89 and 93 GeV, on or near to the Z resonance.It was required that the VD, ID, TPC and HPC were in good operational condition for the data sample analysed.The integrated luminosity of the data sample was 135 pb −1 of which about 100 pb −1 was taken at √ s of 91.3 GeV, close to the maximum of the Z boson resonance production cross-section.
The selection procedures were studied with samples of simulated events which had been passed through a detailed simulation of the detector response [7] and reconstructed with the same program as the real data.Separate samples were produced corresponding to the detector configurations in different years.The Monte Carlo event generators used included: KORALZ 4.0 [8] for e + e − → τ + τ − events; DYMU3 [9] for e + e − → µ + µ − events; BABAMC [10] and BHWIDE [11] for e + e − → e + e − events; PYTHIA 5.7 [12] for e + e − → qq events.Four-fermion final states were produced using two different generators.BDK [13] was used for reactions with four leptons in the final state.This included two-photon events where one or two e + or e − were not observed in the detector.TWOGAM [14] was used to generate e + e − → (e + e − )qq events.The KORALZ 4.0 generator incorporated the TAUOLA 2.5 [15] package for modelling the τ decays.A value of 1.2% was used for the branching ratio of the Dalitz decay π 0 → γe + e − .

Method
From a high purity sample of Z → τ + τ − events, τ decays were identified with charged multiplicity one, three or five, corresponding to the topological branching ratios defined as: where h is either a charged π or K meson.In this definition, the K 0 S meson was treated as a neutral particle, even when it decayed to π + π − before the tracking detectors.In τ decays only about 5% of K 0 S mesons decay before the VD.The h − K 0 (neutrals)ν τ component was therefore treated as signal for the 1-prong τ decays and background for the 3-prong τ decays.Although a π 0 meson which decayed to the γe + e − final state produced two charged particles observed in the detector, it was also treated as a neutral particle, following the τ branching ratio definition given by the Particle Data Group [5].Thus the branching ratios B 1 , B 3 and B 5 correspond to the parameters Γ 2 , Γ 55 and Γ 101 respectively in [5].
No attempt was made to measure τ decays to more than five charged particles and it was assumed that such decays had a negligible branching ratio.A 90% CL upper limit on the inclusive branching ratio for the τ decay to seven charged particles has been set at 2.4 × 10 −6 by the CLEO experiment [16].
At LEP it is possible to separate cleanly τ + τ − events from other final states, permitting the efficient selection of a high purity τ sample.One can measure the branching ratio B(τ → X) for the decay of the τ to a final state X using the expression where N X is the number of identified decays of type X found in the sample of N τ τ -decay candidates, preselected with efficiency τ with a background fraction of b τ .X is the total (preselection × identification) efficiency for selecting the decay mode τ → X, with a background fraction of b X , including background from other τ decay modes.In a τ + τ − event selection without specific requirements on either of the two τ candidates in the event, τ will be identical to the τ + τ − event selection efficiency.
In the case where several branching ratios are measured simultaneously, candidate τ decays can be classified as a function of the detector components used in the charged particle track reconstruction.These classes (detailed in Section 5.4) typically have different signal to background ratios for different types of decays and thus contain different levels of information.To obtain the branching ratios, a fit can then be performed to the predicted number N i,pred of decays in a class i: where i,m is the probability of a τ decay of true multiplicity m being attributed to class i and N i,bkg is the estimated background in class i due to non-τ + τ − events.This can be extended to a simultaneous fit for both τ decays in a candidate τ + τ − event.The predicted number N i,j,pred of τ + τ − events with one τ decay in class i and the other in class j is where N i,j,bkg is the estimated number of non-τ + τ − background events in class i, j and mn ij is the probability of an event with decays of true multiplicity m and n being attributed to class i, j.N τ τ is the number of τ + τ − events and can be left free and obtained from a fit to the different classes.This method takes into account correlations between the two observed τ decays in an event and has the advantage that non-τ + τ − backgrounds tend to populate event classes which have correlations between the two candidate τ decays (e.g.dimuons in the 1-versus-1 topology but not in 1-versus-3) and thus exhibits less sensitivity to backgrounds.It was used as the principal method in the following analysis.

e + e − → τ + τ − preselection
In e + e − → Z → τ + τ − events at √ s = M Z , the τ + and τ − are produced back-to-back, ignoring radiative effects.Each τ decays, producing one, three, or more charged particles in addition to one or two neutrinos and, possibly, neutral mesons.All particles apart from the neutrinos can be detected in DELPHI.A τ + τ − event is thus characterised by two low multiplicity jets which appear approximately back-to-back in the laboratory.The loss of the neutrinos implies that not all the energy in the event is seen in the detector.
Except where explicitly stated, all quantities calculated using charged particles used only charged particle tracks which had hits in at least the TPC or OD, to ensure good momentum reconstruction from the long track length, and had an impact parameter with respect to the centre of the interaction region of less than 1.5 cm in r-φ and 4.5 cm in z.
The basic e + e − → τ + τ − preselection described below has some differences compared with that used previously in studies of the τ leptonic branching ratios [17].There are slightly modified cuts to reduce non-τ backgrounds and the addition of subsamples of events containing τ decays with a nuclear reinteraction in the detector material and τ decays where one charged particle was reconstructed only in the VD and ID tracking subdetectors.

Preselection criteria
To ensure that the τ decay products were in the region of DELPHI corresponding to the HPC acceptance, the thrust axis of the event, calculated using only charged particle tracks, was required to lie in the polar angle region defined by | cos θ| < 0.732.The event was split into two hemispheres, each associated to a candidate τ decay, by the plane perpendicular to the thrust axis and containing the point at the centre of the interaction region.Each hemisphere had to contain at least one charged particle.It was required that at least one charged particle lie in the polar angle region defined by | cos θ| > 0.035 in order to reduce the effects of the TPC tracking inefficiency near cos θ = 0.
Most hadronic Z decays were rejected by requiring that the event contain at most eight charged particle tracks with hits in either the TPC or OD and satisfying the impact parameter cuts described above.
Four-fermion events were suppressed by requiring that the event isolation angle be greater than 160 • .This was defined as the minimum angle between any pair of charged particles in opposite τ decay hemispheres.
Backgrounds from µ + µ − and e + e − final states and cosmic rays were reduced by requiring that the isolation angle be less than 179.5 • for events with only two charged particles.
The µ + µ − and e + e − contamination was reduced further by requiring that both 1/2 and E rad = ( ) 1/2 be less than unity.The variables p 1 and p 2 are the momenta of the highest momentum charged particles in hemispheres 1 and 2 respectively.The quantity p 1 was obtained from the formula and p 2 by analogy with the indices 1 and 2 interchanged.The angles θ 1 and θ 2 are the polar angles of the highest momentum charged particle in hemispheres 1 and 2 respectively.
The variables E 1 and E 2 are the total electromagnetic energies deposited in cones of halfangle 30 • about the momentum vectors p 1 and p 2 respectively, while E j = cp j , for j = 1, 2. These definitions of p rad and E rad are identical to those used for the e + e − → τ + τ − measurements at LEP-2 [18].They take into account more correctly the radiative events where the radiative photon is close to the beam axis than the definition used previously in LEP-1 analyses [17,19].In the limit of no radiation the new and old definitions are equivalent.The distributions of p rad and E rad are shown in Fig. 1a and Fig. 1b respectively, after all cuts have been imposed except the one on the displayed variable.There are some small data/simulation discrepancies in the momentum resolution, visible in the dimuon peak at p rad ≈ 1.4.This had a negligible effect on the analysis.Much of the remaining background from the dileptonic channels came from events containing hard radiation lying far from the beam.These events should lie in a plane.Where two charged particles and a photon were visible in the detector, such events were removed by requiring that the sum of the angles between the three particles was greater than 359.8 • .
A further reduction in four-fermion background was achieved by requiring a minimum visible energy in the event.The visible energy E vis is the sum of the energies of all charged particles and the electromagnetic energies of all neutral particles, neglecting energy deposits recorded by the very forward calorimeters (SAT and STIC) at angles of less than 12 • from the beam axis.Events were accepted if E vis was greater than 0.09× √ s.For events with only two charged particles, an additional condition was applied that the vectorial sum of the charged particle momenta in the r-φ plane be greater than 0.4 GeV/c.Four-fermion events typically have very low values of of this quantity compared with τ + τ − events.
Most cosmic rays were rejected by the upper cut on isolation angle.Further rejection was carried out by requiring at least one charged particle in the event have an impact parameter with respect to the interaction region of less than 0.3 cm in the r-φ plane and that in the τ decay hemisphere opposite to this particle at least one charged particle have an impact parameter in the r-φ plane of less than 1.5 cm.It was also required that both event hemispheres have a charged particle track whose perigee point lay within 4.5 cm of the interaction region in z.Additional τ + τ − events were selected where one of the τ 's produced a 1-prong decay and the other τ decay produced one charged hadron which interacted with the detector material inside the TPC and where the interaction was reconstructed by the nuclear reinteraction reconstruction algorithm described in Section 5.3.The cuts in p rad and E rad were not applied.This increased the selection efficiency by about 1% with a relative background of 0.8%.
A further class of events was selected where one of the τ decays left one charged particle track with hits in only the VD and the ID due to inefficiencies in the TPC and OD.Many of these tracks contained no polar angle information and had poor momentum resolution, so they were excluded from the calculation of event quantities such as the thrust or p rad .The p rad requirement was replaced with a cut on the momentum of the highest momentum charged particle in the opposite hemisphere, requiring it to be less than 90% of the beam momentum.This cut was tightened to 70% if there was an identified muon in the hemisphere.This subsample of events contributed an extra 4% of events with a relative background of 5%.
The total efficiency, including all classes of events, was estimated from simulation to be (51.74± 0.04)%.This corresponds to an efficiency of approximately 85% within the angular acceptance cuts.

Backgrounds
Knowledge of the background from qq events is of particular importance in this analysis as it constitutes a serious background in the higher multiplicity topologies.Uncertainties in the low multiplicity backgrounds have a reduced effect on the measured branching ratios due to the large B 1 value and the resultant cancellation in Eqn. 1.After the cuts described above, the qq background estimated from simulation is (0.80 ± 0.05)%, where the error is purely statistical.To reduce systematic effects due to qq background, a further level of rejection was performed.A number of observables display separation power between τ + τ − and qq events while being only weakly dependent on the topology classification.These were: the isolation angle; the event thrust; the number of neutral electromagnetic showers in an event; the invariant masses reconstructed using charged particles in each candidate τ decay hemisphere of an event, and the invariant masses reconstructed using both charged and neutral particles in each candidate τ decay hemisphere.These variables were used as inputs to a neural network with one hidden layer of eight neurons and a single output neuron.The output neuron distribution is shown in Fig. 2. A cut at 0.05 in the output neuron gave a rejection factor of about three, leaving a background of 0.29% from e + e − → qq, while rejecting only 1.3 × 10 −3 of the remaining τ + τ − events.Discrepancies between data and simulation in this variable are discussed in Section 6.2.1.
Background levels were estimated as accurately as possible using the data.Variables sensitive to a particular type of background were chosen.The background levels were extracted by fitting for the background contribution in these variables assuming the background shape from the simulation.Typically a relative precision of order 10% or better was achieved, ensuring systematic uncertainties on the topological branching ratios below the expected statistical precision.
The background from e + e − → qq was estimated using the data by performing a fit to the distribution of the neural network output neuron.The shape of the background distribution was taken from simulation and the fit was performed for the normalisation factor.A factor 0.92 ± 0.03 was obtained compared to the simulation.
The background due to the fully leptonic four-fermion final states was estimated from the data by extrapolating from the observed background in the low isolation angle region for events with a low visible energy, as well as for events tagged by the presence of a muon (for e + e − µ + µ − ) or an electron (for e + e − e + e − ).The relative uncertainty obtained on the normalisation for these final states was ±10%.Fig. 3 shows the isolation angle distributions for different classes of events.The background from dimuon events in the τ + τ − sample was estimated by studying the momentum distribution for high momentum identified muons, and by studying the p rad distributions for events with two identified muons.This yielded a correction factor relative to the simulation of 0.96 ± 0.03.The p rad distribution, after all other τ + τ − preselection cuts, of events containing two identified muons is shown in Fig. 4a.Similar studies were carried out for Bhabha background using the associated electromagnetic energy and anti-electron tagging in the hadron calorimeter and muon chambers.This, combined with studies [17] comparing the BHWIDE and BABAMC Bhabha generators indicated that a rescaling upwards by 1.15 ± 0.15 of the Bhabha background estimated from simulation was necessary.In addition, further cross-checks were performed by studying the p rad and E rad distributions in the region outside the one occupied by τ + τ − events.The E rad distribution, after all other τ + τ − preselection cuts, for events with two identified electrons is shown in Fig. 4b.
The cosmic ray contamination was estimated from the data to be (50 ± 8) × 10 −5 by extrapolating from the number of events seen outside the cuts in impact parameter.

Charged particle tracking
The charged particle track pattern recognition was based on a robust algorithm designed to minimise the number of bad associations between hits produced by different particles in the different elements of the tracking system.In general, most tracks contain hits in all the tracking subdetectors, the VD, ID, TPC and OD.Most importantly, for the rejection of photon conversion products, the attachment of VD hits to a track must be as efficient and as unambiguous as possible.The VD efficiency and misassociation probability in simulation have been tuned to the data with e + e − → µ + µ − events and with large samples of hadronic events for the purposes of heavy quark tagging [20].
The TPC was the most important subdetector for charged particle reconstruction.The section of a track reconstructed in the TPC had to contain at least three space points out of a maximum possible 16.In dense topologies, where tracks lay closer together than the inherent two-track resolution of the detector, the clustering algorithm used to reconstruct space points used the information contained in the φ profile of a single cluster to resolve two tracks down to distances of typically 2 to 3 mm in r-φ.Where tracks were closer than this, a space point was associated simultaneously to each of the tracks with which it was compatible and assigned a large error so as not to introduce systematic biases in the reconstructed parameters for these tracks.For regions away from the TPC sector boundaries, the probability to include either the TPC or the OD on an isolated track was close to 100%.
After the data had been passed through the standard algorithm, a second algorithm was run which produced VD-ID tracks from combinations of hits in the VD and the ID jet chamber.This algorithm recuperated charged particle tracks around the boundary regions of the TPC or cases where the two-track resolution of the TPC prevented it from resolving both charged particles.The two-track resolution of the VD is more than one order of magnitude better in r-φ than the TPC.However the VD-ID tracks had poorly reconstructed momenta and particle charge sign determination due to their short length.This had no consequences for the charged particle multiplicity reconstruction.

Photon conversions
A good understanding of the photons converting in material before the tracking subdetectors is important in order to measure correctly the true charged particle multiplicity in a τ decay.About 7% of photons interact with the material before the TPC gas volume, giving an e + e − pair detected in some of the tracking chambers.In particular, unreconstructed conversions form an important potential background from lower multiplicity topologies in higher multiplicity samples.Furthermore, to measure in data both the conversion reconstruction efficiency and the amount of material, it was necessary to study the rate of both reconstructed and unreconstructed converted photons.
Converted photons were reconstructed using the algorithms described below.These are described in more detail in [21].Photons converting in the beam-pipe, VD, ID and TPC inner wall were reconstructed using charged particle tracks observed in the TPC.For each TPC track the position of the tangent to the helix which passed through the interaction region was calculated.This point gave the estimated position of the conversion.All tracks for which this point was compatible with the interaction region within one standard deviation were neglected as conversion candidates.If two conversion candidates were found with compatible conversion points, they were accepted as a converted photon provided that: the two tracks had opposite curvature; the mean conversion radius of the two tracks was greater than 5 cm; at least one of the tracks had no associated point inside the mean conversion radius; the angles of the two tracks at the conversion radii agreed within 30 mrad in φ and 15 mrad in θ.This algorithm had an efficiency, estimated from simulation, of about 60% for photon energies below 6 GeV falling to about 20% at 10 GeV.It was complemented by an algorithm using a fit to a joint origin for two charged particle tracks assuming zero opening angle, with similar requirements on the position of the reconstructed conversion radius to the previous algorithm.This had an efficiency of about 40% for photon energies above 8 GeV.The average efficiency for the combination of both algorithms was (68.1 ± 0.2)% for the photon energy spectrum in simulated 1-prong τ decays and (59.8 ± 0.4)% for 3-prong decays.
In 3-prong decays, the probability of reconstructing a π + π − pair incorrectly as a photon conversion was estimated from simulation to be 1 × 10 −4 .For systematic studies the uncertainty on this was taken as 100% of the size of the effect.The probability of reconstructing a K 0 S → π + π − decay incorrectly as a photon conversion was estimated from simulation to be 1.5%, however this did not alter the observed charged particle multiplicity.
The rates of reconstructed and unreconstructed converted photons were studied as a function of the radial coordinate of the estimated conversion point.The conversions were separated into those taking place in the beam-pipe, VD and ID jet chamber inner wall, the ID jet chamber sensitive volume, the ID jet chamber outer wall and ID Trigger Layers, and the TPC inner wall.The bulk of the material was in the TPC inner wall.In order to be as independent as possible of the τ + τ − presample and modelling of the photon production rate in τ decays, the study was performed using e + e − → µ + µ − γ and e + e − → e + e − γ radiative events.The events were selected using electron or muon identification and used kinematic constraints to predict the direction and energy of the photon.The rates were recorded for events with an HPC shower, reconstructed conversion, or no reconstructed conversion but some charged particle tracks consistent with an electron-positron pair arising from a converted photon.In the latter case the radius of the conversion point was taken as being that of the material layer lying just inside the innermost first measured point on either of the two tracks believed to be due to the conversion.The rates were normalised to the total number of observed photons, including conversions.The identical procedure was applied to the simulation.Correction factors for the reconstructed conversion rate were obtained from the ratio A The subscripts d and s stand for data and simulation respectively.R 0 d,s is the total rate of photons, X EM d,s is the probability of a conversion taking place (proportional to the amount of material in terms of radiation lengths) and conv d,s is the reconstruction efficiency.The ratio A conv unrec is given by  The ratio R 0 d /R 0 s of the rates of photons produced in the e + e − → µ + µ − γ and e + e − → e + e − γ radiative events was taken to be unity with a systematic uncertainty of ±2% [22].
The factors A conv rec and A conv unrec obtained as a function of interaction radius are given in Table 1.These factors were used to reweight the simulated events used in the estimation of the selection efficiencies to account for data/simulation discrepancies in the material and reconstruction efficiency.The distributions of the conversion radius and energy for reconstructed converted photons are shown in Fig. 5, where the simulation has been corrected with the factors obtained in the e + e − → µ + µ − γ and e + e − → e + e − γ events.There are still some localised discrepancies between data and simulation due to simplifications in the detector material description in the simulation program.For the τ presample, the rate of reconstructed conversions observed in data was 1.004±0.008times that predicted by the simulation after the application of the correction factors, in excellent agreement.
The data taken in different years were studied separately, in particular before and after the beginning of 1995, the point after which the material in the ID trigger layers was substantially reduced.Consistency was observed between the different years and the change in conversion rate observed in the ID trigger layers between 1992-1994 and 1995 was consistent with the expectation.
The conclusion, integrating up to a radius of 50 cm, is that the material is underestimated in the detector simulation by about 10% in terms of radiation lengths, while the conversion reconstruction algorithm efficiency is similar in data and simulation.
The kinematical conversion reconstruction was complemented by conversion rejection based on an electron identification algorithm [7] which combined information from the TPC ionisation and the HPC electromagnetic shower reconstruction.This algorithm also rejected many π 0 Dalitz decays, which was particularly important for rejecting backgrounds in the 5-prong τ decay sample.A conversion was flagged if there was at least one tightly tagged electron in a hemisphere with greater than one track in the TPC or if there were a pair of oppositely charged particles both of which were loosely tagged as electrons.These criteria were chosen to minimise the probability of misidentifying a π + π − pair as e + e − rather than to maximise the conversion identification efficiency.In simulation, the efficiency to reject the remaining conversions was (59.6±0.6)% in 1-prong decays and (54.5±0.8)% in 3-prong decays.The electron-positron pair rejection also identified (72.7±1.2)% of Dalitz pairs from π 0 decays in simulated 1-prong decays, giving an improved classification for these decays.The efficiency was (68.3±3.1)% for 3-prong decays.The efficiency to identify Dalitz pairs was higher than for conversions due to biases arising from the kinematic conversion reconstruction algorithm which had already identified and reconstructed most converted photons.
The conversion rejection efficiency of the electron identification algorithm was crosschecked in τ decays containing three charged particle tracks, from the rate of oppositely charged pairs of particles which satisfied the electron tagging requirements outlined above and where neither particle had a hit in the VD.The ratio of the efficiency in data to that in simulation was estimated to be 0.96±0.05, in good agreement with unity.The relative systematic uncertainty on the efficiency of this algorithm was conservatively taken to be 10%, to account for any systematic effects in the estimation of the ratio and in its application to all classes of charged particle tracks.
Fig. 6 shows the invariant mass distribution for pairs of oppositely charged particles in τ decay hemispheres containing three charged particles, where one or both of the pair were tagged as an electron candidate and where both tracks had hits in the VD.The electron mass was assumed for both particles.The VD association requirement rejected most photon conversions, which occurred after the VD, while not affecting the signal of π 0 Dalitz decays, or of hadrons misidentified as electrons.The contribution from Dalitz decays peaks at low masses, below the π 0 mass.The tail from Dalitz decays in the region of masses around 200 MeV to 600 MeV arises from decays where only one of the electron or positron was identified and there were two possible combinations of oppositely charged particles.In such decays only half the possible combinations were due to the e + e − pair.In these cases each combination has been given a weight of one half.The region of mass greater than 200 MeV was dominated by π + π − pairs where one of the pions was misidentified as an electron.This permitted a determination from the data of the probability to misidentify a 3-prong τ decay as a 1-prong decay containing a conversion candidate.The misidentification probability was estimated to be (1.29 ± 0.09)%.The uncertainty was dominated by data statistics and contained a contribution from the VD association efficiency.From the ratio of data to simulation in the low mass spike it was estimated that the efficiency to reject a π 0 Dalitz decay in data was 0.90±0.07times the efficiency obtained in simulation.The uncertainty included contributions from the error on the Particle Data Group best estimation of the π 0 → γe + e − branching ratio [5], the uncertainty in the VD efficiency and the limited simulation statistics in addition to the statistical uncertainty.No correction was applied to the simulation and a relative uncertainty of 10% on the Dalitz rejection efficiency was taken.Most δ-rays were rejected by the requirement of VD hits on a track.Events in simulation containing a δ-ray were reweighted with the radius dependent factors describing the ratio of material in terms of radiation lengths 1.There is some ambiguity in the choice of this factor.The associated systematic uncertainties are discussed in Section 6.2.3.

Nuclear reinteractions
About 3% of hadrons reinteract inelastically with the material before the TPC gas volume, and produce typically up to 10 charged secondary particles.These reinteractions were reconstructed by an algorithm which was designed to find secondary reinteraction vertices using the tracks from outgoing charged particles produced in nuclear reinteractions.
The algorithm was robust with respect to the number of charged particles produced and to the position of the reinteraction with the DELPHI detector material.Only tracks which had impact parameters inconsistent with coming from the interaction region were considered.The algorithm iterated over all pairs of tracks, and attempted to fit a vertex in three dimensions if neither of the tracks had already been included in an earlier successfully fitted vertex.However, if one of the tracks had already been included in an earlier fitted vertex, an attempt was made to include the other track in that vertex.This was the main part of the algorithm permitting the reconstruction of vertices with an arbitrary number of outgoing charged particle tracks.In the rare case where both tracks had already been associated to different fitted vertices, an attempt was made to merge those two vertices, if their positions were compatible.In all cases a candidate reinteraction vertex had to pass a χ 2 cut and all outgoing tracks had to have a first observed hit consistent with the particle originating from the reconstructed vertex.Candidate vertices with two outgoing tracks and kinematically compatible with the decay of a K 0 S or Λ produced in the interaction region were not considered.
A search was then performed for the incoming particle which caused the interaction.Tracks which had been reconstructed in only the VD or the VD and ID were extrapolated to the candidate nuclear reinteraction vertex, and associated to the vertex if consistency was found.This prevented double counting of VD-ID tracks and reconstructed nuclear reinteractions.If there was no such track the momentum and charge of the incoming hadron were estimated from the sum over the outgoing particles associated to the vertex.For τ decays it is a good approximation to assume all such nuclear reinteraction vertices are caused by charged particles, as K 0 mesons are the only source of neutral hadrons producing nuclear reinteractions and are infrequently produced.
In simulation, the algorithm was (82.4±0.2)%efficient in 1-prong τ decays for incoming particle momenta above 1 GeV/c, and (83.2±0.4)%efficient in 3-prong τ decays.The efficiency had little dependence on the τ decay multiplicity or on the charged secondary particle multiplicity in the nuclear reinteraction.The position of the reinteraction vertex was reconstructed with a precision of better than 1 mm in three dimensions.
In a similar manner as for the converted photons, the reconstructed and unreconstructed nuclear reinteractions were studied as a function of the reinteraction radius.The rates of reconstructed nuclear reinteractions were obtained from the τ + τ − sample in data and simulation.The radial distribution of reconstructed hadronic reinteractions is shown in Fig. 7. Some localised discrepancies, more marked than in the case of the photon conversions, are visible in the detailed description of the ID trigger layers and TPC in the simulation.The flat tail for radii greater than 35 cm is due to reinteractions in the TPC cathode plane at θ ≈ 90 • .
The unreconstructed nuclear reinteractions were measured by studying the multiplicity distribution for charged particle tracks with impact parameters in the r-φ plane greater than 1.5 cm in the hemisphere opposite a well reconstructed leptonic τ decay candidate.Simulation studies showed that this was dominated by τ decays with an unreconstructed nuclear reinteraction for multiplicities of four and greater.The τ decay hemispheres used to estimate this effect had to have no reconstructed nuclear reinteraction.This removed the risk of a false measurement arising from potential discrepancies in the charged particle multiplicity due to data/simulation differences in the track inclusion efficiency of the nuclear reinteraction reconstruction algorithm.The reinteraction radius was taken as the .The region of the ID jet chamber gas volume, (13 cm to 22 cm) has been included together with the region at lower radius as it has a very small number of nuclear interaction lengths, of order 10 −4 .These results indicate that for nuclear interactions the material at lower radii than the TPC gas volume is over-estimated by about 9% in the simulation program while the reconstruction efficiency is consistent in data and simulation.There is a significant difference in the estimated data/simulation material ratios for nuclear and electromagnetic interactions.This can be attributed to the complexity of the composite material structure together with the means by which each of these quantities is introduced in the simulation program, as well as to possible weaknesses in the modelling of the hadronic reinteractions in the simulation program.The same effects are observed in e + e − → qq events.The uncertainty in the initial production rate of hadrons in τ decays is relatively small and does not give an important systematic uncertainty in the measurement of A nucl rec and A nucl unrec .As for the photon conversions, the nuclear reinteractions were studied for each year separately, particularly before and after the upgrade of the ID trigger layer in 1995.
Simulation studies indicate that the probability of reconstructing a fake nuclear reinteraction from τ decay secondary particles where none was present was 5×10 −4 in 3-prong τ decays.However most such events were still classified as 3-prongs.The associated uncertainty was taken to be 100% of the size of the effect.The probability of reconstructing a fake nuclear reinteraction from a K 0 S decay was negligible.The probability that a hadron produced in a τ decay would undergo an elastic nuclear reinteraction before the TPC sensitive volume was estimated from simulation to be (3.29±0.05)%.The typical signature of such an interaction was a charged particle track with a kink at the scattering point in the detector material.For large scatters the reconstruction algorithm would not reconstuct one but two tracks, one before and one after the scatter.This led typically to a signature of one VD-ID track and one TPC-OD track, and did not alter the observed track multiplicity in the TPC and OD.

Topology selection criteria
By judicious application of the VD hits requirement, and taking into account that the τ must decay into an odd number of charged particles, algorithms were constructed to minimise the sensitivity to the VD association efficiency and the knowledge of the material reinteractions; these are described below.In the following a "good" track is defined as a charged particle track which has associated TPC or OD hits and is not identified as originating from a conversion or nuclear reinteraction in the detector material.For simplicity, the VD-ID track classification, in addition to including tracks reconstructed using only the VD and ID detectors, includes tracks reconstructed as the ingoing charged hadron to a nuclear reinteraction.In τ decay hemispheres where the presence of a converted photon was signalled by the electron identification algorithm, the observed number of good tracks was reduced by two (or by one in the rare case of an even number of good tracks) for the purposes of the topology classification.
Individual τ decay hemispheres were allotted to five different classes dependent on the number of reconstructed tracks in the TPC, VD, etc.The first three of these classes correspond closely to events with clean reconstruction of a 1-prong, 3-prong or 5-prong topology.The other two classes contain the small number of τ decay hemispheres where there was more ambiguity in the reconstruction of the topology, but where some discrimination was still possible.
A 1-prong τ decay was defined as a τ decay hemisphere satisfying any of the following criteria: • only one good track with associated VD hits, and no other tracks with associated VD hits; • only one good track, without VD or ID hits, and one VD-ID track; • no good tracks, and only one VD-ID track.
3-prong τ decays were isolated by demanding τ decay hemispheres satisfying at least one of the following sets of criteria: • three, four or five good tracks, of which either two or three had VD hits; • two good tracks with associated VD hits, plus one VD-ID track; • one good track with associated VD hits, plus one or two VD-ID tracks pointing within 3 • in azimuth of a TPC sector boundary.
Candidate 5-prong τ decays were selected if they satisfied at least one of the following topological criteria: • five good tracks of which at least four had two or more associated VD hits; • four good tracks with VD hits, and one other VD-ID track.
Additional criteria were applied in the selection of 5-prong τ decays due to the large potential background from hadronic Z decays and mis-reconstructed 3-prong τ decays.The background originating from ν τ 3h ± ≥ 1π 0 final states with a Dalitz decay was expected to occur at a similar level to the signal.The electron rejection criteria described in Section 5.2 reduced this background by about 70%, and it was further suppressed by requiring that all good tracks had a momentum greater than 1 GeV/c.To reject Z → qq events it was required that the total momentum of the 5-prong system be greater than 20 GeV/c.Only good tracks were included in the calculation of this quantity.These three classes accounted for 97.6% of candidate τ decays in the τ + τ − sample.The remaining 2.4% of candidate τ decays were mostly 1-prong and 3-prong τ decays with some pattern recognition failure or detector inefficiency.These were classified into two categories, 1 and 3 , corresponding to 1-prong and 3-prong τ decays respectively.Candidate τ decay hemispheres of type 1 had to satisfy at least one of the following criteria: • one good track, with associated VD hits, plus two VD-ID tracks not pointing within 3 • in azimuth of a TPC sector boundary; • one good track, plus a conversion candidate selected with the electron identification algorithm.
The 3 category contained all other τ decay hemispheres failing to pass any of the previous classification requirements.In particular it contained those events with either four or five good tracks of which either four or more had associated VD hits, but not satisfying the 5-prong class requirements.This criteria had a higher efficiency for 5-prong decays than for 3-prong decays, but due to the much greater 3-prong branching ratio the 3 category was dominated by 3-prong decays.Table 3 contains the efficiencies of these selection requirements for the different exclusive τ decay modes and the inclusive single-hemisphere topological selections, as obtained from simulation and after corrections for observed discrepancies between data and simulation in the rate and reconstruction efficiency of material reinteractions.Events were classified according to both the hemisphere classifications into 14 event topology classes: 1-1, 1-1 , 1-3, 1-3 , 1-5, 1 -1 , 1 -3, 1 -3 , 1 -5, 3-3, 3-3 , 3-5, 3 -3 and 3 -5.The class 5-5 was not included as these events were removed by the multiplicity cut in the τ + τ − preselection and in any case would give a negligible signal contribution.The selection of event classes up to multiplicity five assumes that the inclusive branching ratio of the τ to seven charged particles is negligible.While inclusion of higher multiplicity classes is possible, the DELPHI sample size is insufficient to reach the level of precision obtained by CLEO [16].
The various sources of non-τ + τ − background are detailed in 6 The fit and systematics

Fitting procedure
A reweighting technique was used to take into account the observed data/simulation discrepancies in the rates and reconstruction efficiencies of nuclear reinteractions, photon conversions and δ-rays.Each secondary or tertiary particle i in a simulated τ + τ − event k was given a weight W part i .For particles interacting with the detector material this weight was obtained from the studies described in Section 5 and summarised in Tables 1 and 2. These weights were a function of the radius and type of the reinteraction, as follows:
• reconstructed converted photons: It was also necessary to reweight particles which did not reinteract so as to maintain the internal normalisation of the simulation sample with respect to parameters such as angular distributions, momentum distributions, or the τ exclusive branching ratios.Failure to achieve this could lead to biases within the simulation sample as a function of photon multiplicities, charged hadron multiplicities, angular regions with different amounts of material, or could bias the preselection of the τ + τ − sample with respect to certain types of decay mode.Non-interacting photons were given a weight W part i given by where N γ is the total number of photons produced by τ decays or by interactions with material of τ decay products in the simulation sample.N γ int is the number of photons with a material reinteraction.The sum is over all photons with a material reinteraction.To maintain the angular and momentum distributions the above weight calculation was performed separately for individual bins in a three-dimensional space of the photon polar angle, azimuthal angle and momentum.A similar procedure was used to reweight noninteracting hadrons, particles not emitting a δ-ray (a very small correction) and electrons without any bremsstrahlung emission.
With every secondary and tertiary particle i in a simulated event k given a weight W part i , the weight W event k of the event was given by the product of all these weights: Class Observed Fit Output Table 5: The second and third columns contain the number of observed events in each class and number predicted from the maximum likelihood fit of the τ topological branching ratios.The fourth column contains the χ 2 contribution for each class to the overall χ 2 .
The numbers of selected events in each event topology class are shown in Table 5.A maximum likelihood fit assuming Poissonian probabilities was performed to the reweighted numbers of events estimated using Eqn. 3. The number of τ + τ − events (N τ τ in Eqn. 3) and the branching ratios B 1 and B 5 were allowed to vary in the fit while B 3 was constrained by the relation B 1 + B 3 + B 5 = 1.The output of the fit is also shown in Table 5.The results of this fit were: B 1 = (85.316± 0.093)%; B 3 = (14.569± 0.093)%; B 5 = (0.115±0.013)%,where only the statistical error from the fit is quoted.An estimate of the consistency of the fit was made by calculating a χ 2 .This took into account only the statistical uncertainties.It gave a χ 2 /n.d.f. of 16.4/11, indicating good consistency.The contribution to the χ 2 from each class is shown in Table 5.

Systematics
In general, the systematic uncertainties on the topological branching ratios due to any particular effect were estimated simultaneously for B 1 , B 3 and B 5 by repeating the analysis, including the τ + τ − preselection, after modifying the relevant variable in the simulation.This accounted for correlations in the systematic uncertainties between the different branching ratios and correlations in efficiencies and backgrounds between the different event classes and the τ + τ − preselection.

Preselection
The main systematic effects of the τ + τ − selection criteria on the result can arise through the mis-calibration of the quantities used in the selection.These quantities can be classified into cuts related to energy or momentum measurements, such as E rad , p rad and E vis , or those cuts related to multiplicity, such as the isolation angle, which can be incorrectly estimated if extra secondary particles are produced in reinteractions with detector material.This second effect is taken account of by the systematics in the reweighting of the secondary interactions, discussed in Section 6.2.3.A separate systematic uncertainty is included for the energy and momentum scales below.Other effects such as tracking efficiency, trigger efficiency and uncertainties on the τ exclusive branching ratios can have an effect on the τ + τ − preselection efficiency.These sources of systematic uncertainty are discussed below, and in general the systematic uncertainties on the topological branching ratios take into account effects in the preselection.
Any remnant discrepancies due to τ + τ − preselection criteria were checked by studying the agreement between data and simulation for the distributions of the τ + τ − selection variable for each of the different event topology classes.With the corrections to the background levels discussed in Sections 4.2 and 5, good agreement was found, in particular in the regions of the cuts.
The uncertainties due to the backgrounds from non-τ + τ − sources were estimated by varying the background normalisations by their uncertainties obtained as described in Section 4, and are listed in Table 6.The distribution of the neural network variable used for qq rejection (see Fig. 2) displays a 15% excess in the data compared with simulation in the region [0.05;0.8] of the output neuron distribution.This region was dominated by τ + τ − events and was studied to see if there were any discrepancies in the input variables, and by checking the multiplicity distributions of the events.No discrepancies were found.The systematic uncertainty was estimated by rescaling the numbers of rejected signal events in all classes by 15%.

Tracking
Within the angular acceptance of this analysis, there are four tracking detectors which can contribute to the track reconstruction of all charged particles: the VD, ID, TPC and OD.The reconstruction of "good" tracks is strongly dependent on the TPC with its full three-dimensional readout.The redundancy in the track reconstruction obtained by the inclusion of τ decay hemispheres with only a VD-ID track in the τ + τ − sample and 1-prong subsample reduces greatly the sensitivity of the measurement to inefficiencies within a given subdetector, and allows direct cross-checks to be made.
In simulation, for non-interacting particles, the efficiency for the TPC to reconstruct an isolated charged particle which passes through a sensitive gas volume far from dead regions is 99.95%.However, even with the redundancy of the different tracking subdetector components in DELPHI, to measure this in data directly is difficult because of material reinteractions between different subdetectors which can cause the particle to be lost before entering the TPC sensitive volume.These effects tend to reduce the measured efficiency.For the data, the efficiency of the TPC in its sensitive regions was estimated using the redundancy of the tracking system to be (99.35± 0.05)%.The identical procedure was applied to the simulation, yielding an efficiency of (99.31 ±0.01)%, in excellent agreement with the data, and implying that the modelling of the TPC efficiency in the simulation was accurate.The true inefficiency in simulation is a factor of 12 lower than the inefficiency estimated by this method.The uncertainty on the TPC reconstruction efficiency in the TPC sensitive region was taken conservatively to be 0.05% to account for any systematic effects in its estimation.
The TPC efficiency for isolated tracks in the sector boundary region was studied using inclusive low multiplicity events but excluding the τ + τ − events.This sample consists of radiative dilepton events and the high energy part of the two-photon leptonic event spectrum.The number of reconstructed tracks in the TPC boundary region was compared to that expected by normalisation of the sensitive region of the TPC sectors, in bins of momentum.The loss of tracks was compatible with zero below about 5 GeV/c, rising to 3.5% at 45 GeV/c.Data and simulation agree well in their behaviour.The inefficiency estimated from data for the τ + τ − momentum distribution was (2.73±0.04)%,compatible with the result in simulated τ + τ − events.The TPC efficiency was varied by throwing away tracks containing TPC hits in simulation and repeating the analysis, including the τ + τ − preselection stage.The error was scaled by a factor two to account for any systematic effects in the procedure.
The attachment of VD hits to a track is the main criterion used in this analysis to determine the multiplicity of a τ decay.The association efficiency and mis-association probability of a VD hit to be attached to a charged particle track has been studied [20] on large samples of e + e − → qq events, and on isolated topologies such as µ + µ − final states.These studies indicate that the efficiency in simulation and data agree to within ±2%.The associated systematic uncertainties were obtained from the shifts in the results observed when randomly removing 2% of associated VD hits and repeating the analysis, including the τ + τ − preselection stage.
Within the τ sample itself the rate of association of different subdetectors to a reconstructed charged particle track was studied for different τ decay topologies.In candidate 1-prong decays, the hit association probabilities for the VD, ID and OD were studied for tracks which had an associated TPC track segment.The probabilities obtained in data and simulation were compared and the relative differences calculated.These differences were 0.0% for the VD, +2.3% for the ID and −1.2% for the OD.The equivalent numbers for 3-prong decays were −0.7%, +2.5% and +0.6% respectively.The results for the ID were unchanged if in addition the VD was required to be associated to a track, and vice-versa.The small fraction of tracks without a TPC track segment showed a level of agreement for the proportions of tracks with different subdetectors which was better than 3% in all cases.
In simulation, the preselection efficiency was found to be independent of the combination of the subdetectors attached to a track at a level below 0.1%.Given the observed discrepancies between data and simulation, any effect on the preselection efficiency due to this source was O(10 −5 ).The uncertainties arising from the probability of including an ID hit on a track were estimated to be 2.2 × 10 −5 for both B 1 and B 3 with a full anticorrelation and has been included in the tracking systematic uncertainty.The analogous uncertainty for the OD was negligible.
It was observed that the level of Bhabha background was correlated with the existence of an OD hit on a track.This was due to bremsstrahlung, and resolution effects in the p rad variable.The level of agreement observed between the data and the simulation implied that the Bhabha background was consistent with the estimation made in Section 4.2.
The two-track resolution has been studied by data-simulation comparison of the minimum opening angle in three dimensions between tracks in τ decay hemispheres with more than one track.Fig. 8a shows the efficiency in simulation to reconstruct a 3-prong τ decay as a function of the minimum opening angle.It is flat within 1% except for a fall-off of about 5% in reconstruction efficiency for a minimum opening angle below 3 mrad.As is visible in the distribution of the minimum opening angle in candidate 3-prong τ decays shown in Fig. 8b, less than 1% of 3-prong τ decays lie in this region and data and simulation are compatible within the statistical precision.A systematic uncertainty was attributed to the two-track resolution by varying the efficiency in the region with minimum opening angle below 3 mrad by the uncertainty allowed from the statistical uncertainty of the data.There was no observable fall-off in efficiency for low values of the minimum differences between the azimuthal or polar angles of tracks in a 3-prong hemisphere.The minimum opening angle distribution in candidate 5-prong τ decays is shown in Fig. 9.
The association efficiency of the TPC in 3-prong τ decays and its dependence on the minimum opening angle were studied by comparing the rate of candidate τ decays containing three tracks with associated TPC hits with the rate of decays containing either two tracks with TPC hits plus one track without TPC hits or one track with TPC hits plus two tracks without TPC hits.Fig. 10 shows the data over simulation ratio of the fraction of candidate 3-prong τ decays containing three tracks with associated TPC hits, as a function of the minimum opening angle.It shows compatibility with unity in all regions except for a small discrepancy of about +3% for minimum opening angles in the region of 4 to 7 mrad.The average of the ratio for all minimum opening angles is 1.0016 ± 0.0008, consistent with unity within two standard deviations.The difference from unity can be directly related to the efficiency for a charged particle track to have associated TPC hits in 3-prong τ decays.The uncertainties on the topological branching ratios were estimated by varying the probability in simulated 3-prong τ decays to have TPC hits associated to a track by 1.6 × 10 −3 and repeating the analysis.The same procedure was carried out for 5-prong decays.
Fig. 8b shows some differences between data and simulation in the region between 3 and 15 mrad, with the simulation tending to lie slightly above the data.This can be due to the decay modelling, for which the attendant systematic uncertainties are discussed in Section 6.2.4.However track reconstruction effects cannot be excluded.While systematics in the TPC association efficiency are accounted for above, the deviation from unity of the ratio shown in B 1 and B 5 induced via the correlations in the fit.This accounted for potential differences between data and simulation in the VD-ID track reconstruction in the small minimum opening angle region of 3-prong τ decays.The reconstructed charge in a τ decay hemisphere was studied for different event and decay topologies.The charge was not used to classify τ decays or select τ + τ − events, so this provided an indirect cross-check of the track reconstruction.Use was made of the constraint that the two τ -decay hemispheres arose from τ 's of opposite charge.For events in the data belonging to the 1-1 event class, an estimated (99.88±0.01)% of τ decays had the charge correctly reconstructed.The rate in the simulation was (99.91±0.01)%.For 3-prong decays in data, (98.57±0.09)%had the correctly signed unit charge, compared with (98.60±0.03)% in simulation, while (0.54±0.06)% had a charge of three with the correct sign as compared with (0.52±0.02)% in simulation.The remainder had the wrongly signed charge.In 5-prong decays an estimated (90.5±0.7)% had the correctly reconstructed unit charge in data compared with (89.7±0.3)% in simulation.

Reinteractions and K 0
S reconstruction Uncertainties from the photon conversion reconstruction were estimated by varying by their uncertainties the correction factors for the reconstructed and unreconstructed conversions, given in Table 1, which were obtained from data test samples as described in Section 5.2.A contribution for the uncertainty in the incident photon rate in the dileptonic test samples was included.The resultant uncertainties are dominated by the contribution from the unreconstructed conversions.A similar approach was taken for the nuclear reinteractions.Here again the uncertainty arising from the unreconstructed nuclear reinteractions was the most significant contribution.
Both the photon conversion and hadronic reinteraction errors contain contributions to account for the uncertainty on the rate of fake reconstructions which cause the topology to be misidentified.The 1-prong τ decay class has a number of different subclasses which cover the different reconstruction possibilities for an elastic nuclear reinteraction in a 1-prong hadronic τ decay.In simulation, of the 3.3% of hadrons in τ decays undergoing an elastic scatter before the TPC sensitive volume 99.4% were correctly attributed to the 1-prong class, leading to a classification inefficiency of 2.0 × 10 −4 for 1-prong hadronic τ decays.A cross-check of the elastic scatters was performed using a kink algorithm which attempted to link up pairs of tracks consistent with both tracks having been produced by a single charged particle experiencing a large scatter in detector material.Comparison of the rates in data and simulation showed agreement within 20%, which was taken as a relative uncertainty on the effect.The effects of fake kink reconstructions by the algorithm had a negligible effect.
The uncertainty attributed to the electron identification had contributions from the uncertainties on the efficiencies for conversion rejection and for rejection of π 0 Dalitz decays, and from the probability of misidentifying a hadron as an electron.The study of the electron identification algorithm and its systematic uncertainties is described in Section 5.2.
The δ-ray rate is proportional to the mean atomic number Z of the material, while the factors which are applied correct for the number of radiation lengths, which is approximately proportional to the mean Z 2 .This gives an ambiguity in the correction for δ-rays depending on whether it is due to the wrong mean Z, or to the wrong quantity of material, but with the correct mean Z.The correction factor used in the analysis, X EM d /X EM s , assumes the latter case.The former case would imply a factor of (X EM d /X EM s ) 1/2 .This was applied and the observed variations on the branching ratios taken as systematic uncertainties.
Only the fraction (≈5%) of K 0 S → π + π − decays occuring inside the VD, and which the K 0 S reconstruction algorithm failed to identify, gave fake charged primary τ decay products and hence the incorrect topology assignment.For the region inside the outer layer of the VD, the rate of successfully reconstructed K 0 S 's was 0.91 ± 0.09 times the simulation prediction, showing good consistency.The quoted error is purely statistical.
By comparison of data and simulation for candidate K 0 S mesons reconstructed as decaying beyond the VD, the reconstruction efficiency and relative rate of photon conversions misidentified as K 0 S were estimated to agree within 15%.The K 0 S reconstruction algorithm can reconstruct a fake K 0 S from two primary τ decay particles thus reducing the τ decay charged multiplicity by two.Simulation studies indicate that the scale of this effect was 5 × 10 −5 of the 3-prong rate.The related systematic uncertainty was taken to be equal to the size of the effect.

Exclusive branching ratios and τ decay modelling
The exclusive branching ratios were varied within the uncertainties quoted in the Particle Data Listings [5].To take account of any hidden correlations the quoted uncertainties have been scaled up by a factor of 1.5.The largest single contribution is from the decay mode 3π ± 3π 0 ν τ with a high π 0 multiplicity and a large relative uncertainty on the branching ratio.The uncertainty also included a small contribution due to the decay modes τ − → K − π − π + π 0 ν τ and τ − → K − K + π − π 0 ν τ which were not included in the simulation.
The uncertainties associated with the modelling of the 3-and 5-prong decays were estimated by correcting the efficiencies taking into account differences between data and simulated invariant mass distributions.In addition, the hadronic structure of the 3π final state was varied between the default TAUOLA [15] model and that obtained in the DELPHI analysis of the 3π structure in τ decays [23].For the 3ππ 0 structure the parameterisation of Model 1 of [24] was used and, as a cross-check, the parameterisation of 3ππ 0 used in [23] was used to reweight the distributions of the minimum opening angle.

Trigger, energy scale, τ polarisation and simulation statistics
The trigger efficiency for τ + τ − final states was (99.98 ± 0.01)% for events within the polar angle acceptance.Studies indicated that the inefficiency was due to events where both τ 's decayed via the τ → µνν mode [17].It was assumed that the full inefficiency of (2 ± 1) × 10 −4 was contained in the 1-1 event class and the efficiency of this class and of the τ + τ − preselection were modified accordingly.The associated systematic uncertainty was obtained by varying the inefficiency by its error.
The energy and momentum scales and resolution can affect many of the quantities used in the analysis, such as the τ + τ − selection variables, invariant masses or the thrust (used in the qq rejection).The momentum scale was varied by 0.2%, the electron energy by 0.5% and the neutral electromagnetic energy by 0.2%.These variations were obtained in a study carried out for the τ polarisation measurement [25].The associated systematic uncertainty on the topological branching ratio results was small.
The 1-prong selection efficiency has a slight sensitivity to the average τ polarisation because of acceptance effects in τ → πν τ , Kν τ decays due to the p rad cut used in the τ + τ − selection.The analysis used the result and uncertainty from the DELPHI analysis on τ polarisation [25].
The systematic uncertainty due to the limited simulation statistics was also included.

Other cross-checks and summary
A cross-check of the fitting procedure was performed by using the simulation itself as input to the fit; the results of the fit were compared with the input branching ratios and agreement was observed.This was repeated with different input branching ratios, again

Figure 2 :
Figure2: The output neuron distribution for the neural network used for the final qq rejection.The dots are data, the line is simulation and the shaded region shows the simulated qq contribution.The selected events lie in the region with output neuron values above 0.05.

Figure 3 :
Figure3: Distributions of isolation angle (in degrees) in candidate τ + τ − events: a) events containing two identified muons; b) events containing at least one tagged electron; c) all events in the τ + τ − sample.In all cases data are the dots and simulation is the line.In (a) and (b) the shaded area is the four-fermion background, while in (c) it contains all background contributions.

Figure 4 :
Figure4: Distributions used in the study of the dileptonic background to the τ + τ − preselection: a) p rad distribution for events containing two identified muons; b) E rad distribution for events containing two identified electrons.In both cases the dots are the data and the line is simulation.In (a) the shaded area is the background contribution due to the reaction e + e − → µ + µ − (γ), while in (b) it is the background contribution from e + e − → e + e − (γ) events.
conversions in data and simulation respectively.Analogous correction factors for the unreconstructed conversion rate were obtained from the ratio A conv unrec of the rates R unrec d and R unrec s of unreconstructed conversions.The measured correction factors were used to estimate the reconstruction efficiency and amount of material.The ratio A conv rec is given by The conversion reconstruction efficiency in data conv d can be obtained from the data alone via the relation conv d of any assumptions on the material or, if R rec d and R unrec d are measured on the same data sample, on the incident photon flux.If the ratio R 0 d /R 0 s is known, measurements of A conv rec , A conv unrec and conv d , combined with conv s from the simulation, can provide an estimation of the material ratio X EM d /X EM s for electromagnetic interactions.

Figure 5 :
Figure 5: a) Distribution in e + e − → τ + τ − events of the estimated conversion radius for photons converted in front of the TPC.Data are dots, the line is simulation after applying correction factors from e + e − → µ + µ − γ and e + e − → e + e − γ events.b) Reconstructed conversion photon energy in e + e − → τ + τ − events.Data are dots, the line is simulation after applying correction factors.

Figure 6 :
Figure 6: Invariant mass m e + e − distributions for pairs of oppositely charged particles with one or both particles identified as electron candidates and where both tracks have hits in the VD.Data are dots, the simulation is the line.The shaded area shows the contribution from π 0 → γe + e − decays.The region m e + e − > 200 MeV is dominated by misidentified π + π − pairs.

Figure 7 :
Figure 7: Distribution in e + e − → τ + τ − events of reconstructed radius of the inelastic nuclear reinteractions.Data are dots, the line is simulation after applying correction factors.

1 bFigure 8 :
Figure 8: a) The selection efficiency in simulation for candidate 3-prong decays as a function of the minimum angle in three dimensions between any two charged particles.b) The distribution of the minimum angle in three dimensions between any two charged particles in candidate 3-prong τ decays.Dots are data, line is simulation.

1 Figure 9 :
Figure 9: The minimum angle in three dimensions between any two charged particles in candidate 5-prong τ decays.Dots are data, line is simulation.

1 Figure 10 :
Figure 10: For candidate 3-prong τ -decays, the ratio of data to simulation of the fraction of candidate decays where three tracks had associated TPC hits, as a function of the minimum opening angle.

Table 1 :
Factors used to correct simulated photon conversions and data/simulation discrepancies as a function of the radius of the interaction point.Column 5 shows the data/simulation ratio of reconstruction efficiency and column 6 the equivalent ratio for the material.

Table 2 :
Factors used to correct simulated nuclear reinteractions for data/simulation discrepancies as a function of the radius of the reinteraction point.Column 5 shows the data/simulation ratio of reconstruction efficiency and column 6 the equivalent ratio for the material.

Table 3 :
Estimates of the selection and topology classification efficiencies, in percent, for different exclusive decay modes, as obtained from simulation.The efficiencies are corrected for observed discrepancies between data and simulation in the rate and reconstruction efficiency of material reinteractions.The quoted uncertainties are from the simulation statistics only.