A Study of One-Prong Tau Decays with a Charged Kaon

From an analysis of the ionisation energy loss of charged particles selected from 110326 e+e- ->tau+tau- candidates recorded by the OPAL detector at e+e- centre-of-mass energies near the Z0 resonance, we determine the one-prong tau decay branching ratios: Br(tau- ->nu_tau K->=0h0) = 1.528 +- 0.039 +- 0.040 % Br(tau- ->nu_tau K-) = 0.658 +- 0.024 +- 0.029 % where the h0 notation refers to a pi0, an eta, a K^0_S, or a K^0_L, and where the first uncertainty is statistical and the second is systematic.


Introduction
Precise studies of the decays of the tau lepton have been made possible by the availability of large, low background tau-pair samples, such as those produced in e + e − collisions at CESR and at LEP I. This paper reports on an analysis of one-prong tau decay modes containing a charged kaon, using the complete set of data collected by the OPAL experiment between 1990 and 1995 at e + e − centre-of-mass energies near the Z 0 resonance. The excellent charged particle identification capability of the OPAL detector is exploited to obtain precise measurements of the τ − → ν τ K − and τ − → ν τ K − ≥ 0h 0 branching ratios, where the charge conjugate decays are implied in these interactions and throughout this paper. The h 0 notation refers to a π 0 , an η, or a K 0 , where the K 0 can be either a K 0 S or a K 0 L .
Measurements of the τ − → ν τ K − branching fraction may be used to determine the charged kaon decay constant, f K , and the CKM matrix element, |V us |. In addition, precise measurements of this branching ratio may also be used to set limits on physics beyond the Standard Model, such as charged Higgs effects in tau decays [1].

The OPAL Detector
The OPAL detector is described in detail in reference [2]. The component of the detector most important to this analysis is a large volume jet chamber, which measures the momentum and specific ionisation energy loss, dE/dx, of charged particles. The jet chamber provides the only means of distinguishing between charged pions and kaons in the momentum range of interest to this study.
The jet chamber is a cylinder 4 m long and 3.7 m in diameter, and is divided by cathode wire planes into 24 azimuthal sectors. Each sector contains one radial plane of anode wires, which are staggered to resolve left-right ambiguities. The chamber is contained in a solenoidal magnetic field of 0.435 T, and is filled with an argon-methaneisobutane gas mixture at a pressure of 4 atmospheres. This arrangement provides an intrinsic transverse spatial hit resolution of σ xy = 130 µm, and a two-hit resolution of 2.5 mm. 1 In the barrel region of the jet chamber (| cos θ| < 0.68), the ionisation energy loss of a charged particle is sampled up to 159 times. A truncated mean is formed by discarding the highest 30% of the measurements, resulting in a dE/dx resolution of 3.1% for isolated tracks in the chamber [3].
Other tracking detectors include a high precision silicon microvertex detector surrounding the beam pipe, a precision gas vertex detector, and a layer of wire chambers with drift direction parallel to the z-axis lying immediately outside the jet chamber. These z-chambers extend to | cos θ| = 0.80, which encompasses the acceptance of this analysis, and accurately determine the polar angle of charged particles traversing the central detector. A lead-glass electromagnetic calorimeter and presampler chambers are located outside the magnetic coil and jet chamber pressure vessel. The electromagnetic calorimeter is primarily used in this analysis to identify tau decays that include a π 0 or a K 0 in the final state. The return yoke of the OPAL magnet is instrumented for hadron calorimetery, and is surrounded by external muon chambers. The hadron calorimeter and muon chambers are mainly used in this analysis to select control samples containing muons and to veto non-tau backgrounds in the tau-pair candidate sample.

Parameterisation of Ionisation Energy Loss
The mean ionisation energy lost per unit distance by a particle of charge q and speed β traversing a gas is approximately described by the Bethe-Bloch equation [4,5]: where A 1 and A 2 are constants which depend on the gas composition, and δ(β) is a function which describes polarisation effects in the gas.
The ionisation energy loss parameterisation used by OPAL, D pred , is a variation of the Sternheimer-Peierls modification to the Bethe-Bloch formula [3,6], and uses six input parameters. Figure 1(a) shows the dependence of dE/dx on the momentum of various particles in the OPAL jet chamber, as predicted by D pred . Figure 1(b) shows the particle resolving power, R ij , versus momentum for various pairs of particle species i and j, where: pred | / σ, and where σ is the resolution of the measured dE/dx. Note that the OPAL jet chamber yields a pion/kaon separation of at least 2σ for particles between about 2 to 30 GeV/c. The parameters used in the calculation of D pred and σ are tuned to the measured energy loss of charged particles in hadronic Z 0 → qq decays, yielding a dE/dx parameterisation that is accurate enough for nearly all analyses of OPAL data.
However, hadronic Z 0 → qq decays tend to have a much higher track multiplicity than one-prong tau decays, leading to small systematic differences between the measured dE/dx in the two environments. Although these differences are less than one σ in magnitude, they must be corrected for accurate particle identification in one-prong tau decays. Based on a study of the measured energy loss of muons in e + e − → µ + µ − , τ − → µ − ν µ ν τ , and e + e − → e + e − µ + µ − control samples, the following corrections are applied to the the OPAL dE/dx parameterisation, D pred , and to the measured energy loss and its uncertainty, D meas and σ: Multiplicative correction to D pred : To obtain a parameterisation which correctly predicts the measured energy loss of particles from tau decays, D pred must be corrected with a multiplicative factor, s(β). Parameterising s(β) as a second order polynomial in − log(1 − β 2 ) yields satisfactory results. The magnitude of the correction is of order 1% for tracks in one-prong tau decays.
Additive correction to D meas : The measured energy loss of charged particles in lowmultiplicity events is found to depend strongly on the azimuthal separation, φ, between the track and anode plane in the jet chamber cell. To improve the dE/dx resolution, this behaviour is corrected with a function, f (φ), with seven parameters tuned to the measured energy loss of e + e − → µ + µ − candidates. The magnitude of this correction is also of order 1%.
Multiplicative correction to σ: To correctly predict the spread in the dE/dx of charged particles from one-prong tau decays, the OPAL parameterisation of σ must be corrected with a multiplicative factor, s res ≈ 0.9.

Monte Carlo Generated Event Samples
For this analysis, tau lepton decays are simulated with the KORALZ 4.0 Monte Carlo generator and the TAUOLA 2.4 decay package [7,8]. A Monte Carlo sample of 600000 tau-pair events is generated with input branching ratios determined from world averages or theoretical expectations [9]. Background contributions from non-τ sources were evaluated using Monte Carlo samples obtained from the following generators: hadronic Z 0 → qq events were simulated using JETSET 7.4 [10], e + e − → µ + µ − events using KORALZ, Bhabha events using BABAMC [11] and BHWIDE [12], and fourfermion events using VERMASEREN 1.01 [13] and FERMISV [14].
All Monte Carlo samples are passed through a detailed simulation of the OPAL detector [15] and are subjected to the same analysis chain as the data.

Data Selection
This analysis uses the full data set collected by the OPAL detector in the years 1990 through 1995 at e + e − centre-of-mass energies close to the Z 0 resonance. Only data for which the tracking chambers and the electromagnetic calorimeter were fully operational are retained. The topology of e + e − → τ + τ − events is characterised by a pair of backto-back, narrow jets with low particle multiplicity. Jets are defined in this analysis by grouping tracks and electromagnetic clusters into cones with a half-angle of 35 • , where each cone is assumed to contain the decay products of one of the tau leptons. Taupair candidates are selected by requiring two low-multiplicity jets (N track ≤ 6) with an average polar angle of | cos θ ave | < 0.68. Background events from other lepton pairs, hadronic Z 0 → qq decays, and four-fermion events are reduced with cuts on the event topology and total visible energy. These selections produce a sample of 110326 tau-pair candidates with background f non−τ = 1.56 ± 0.10%, estimated by data control samples and Monte Carlo background samples. The tau-pair selection procedure is described in detail in reference [16].
One-prong tau decays are selected from this sample by choosing tau decay cones containing only one well-reconstructed charged track, consistent with coming from the origin. Tracks in this sample are required to have at least 40 jet chamber hits used in the calculation of the measured dE/dx, to have at least 3 hits in the z−chambers, and to have reconstructed momentum between 2 and 50 GeV/c. A sample of 143528 candidates is selected by this procedure.
The number of charged kaons in the sample is determined by maximising a likelihood function based on dE/dx. The τ − → ν τ K − ≥ 0h 0 branching ratio is then obtained by correcting the number of kaons for backgrounds and selection inefficiencies. As a cross-check of the general analysis procedure, the τ − → ν τ π − ≥ 0h 0 branching ratio is also determined in an analogous fashion, and compared to the world average.
The tracks in the one-prong tau decay candidate sample consist of electrons, muons, pions, and kaons from both tau decays and non-tau sources. To determine the number of charged kaons in the sample, a maximum likelihood fit to the ionisation energy loss distribution of the tracks in the sample is performed. The likelihood function used in the fit is: where W kj is the dE/dx weight of charged particle j under particle hypothesis k, and where f k is the fraction of particle type k in the sample, where k is either pion, kaon, muon, or electron. The values of f k are constrained to be non-negative, and their sum is constrained to be 1.
f ′ µ and ∆f ′ µ are the fraction, and uncertainty on the fraction, of muons in the sample, as estimated by Monte Carlo generated events. The estimated muon fraction is corrected using information from data and Monte Carlo muon control samples that are selected using information from the OPAL muon chambers. This constraint is added to the fit to help distinguish between pions and muons in the sample, which have quite similar dE/dx. N is the total number of particles in the sample. D j meas is the measured dE/dx of the j th charged particle. D kj pred is the predicted dE/dx for charged particle j, calculated with the OPAL parameterisation under particle hypothesis k, as obtained from the measured dE/dx of charged particles in hadronic Z 0 → qq events.
σ kj is the dE/dx uncertainty, calculated using the OPAL parameterisation, as obtained from the measured dE/dx of hadronic Z 0 → qq events.
s res is the multiplicative correction to σ kj , as determined from the one-prong control samples as described in Section 2.1.
s(β) is the β dependent multiplicative correction to D pred , as determined from the one-prong control samples.
f (φ) is the φ dependent correction to the measured dE/dx, as obtained from the one-prong control samples.
Efforts are made to obtain a dE/dx parameterisation for the tau decay environment that is optimal for many particle species over a wide range of momenta. However, it is possible that the dE/dx corrections described in Section 2.1 may be somewhat more (or less) optimal for pions than they are for kaons in the momentum range of interest. Thus, to correct for any possible species-dependent quality differences in the parameterisation of dE/dx, an extra factor, C K , is allowed to multiply the kaon predicted energy loss, D K pred , and is allowed to vary freely in the fit using the likelihood function found in Equation 2. The value of C K returned by the fit is C K = 0.9943 ± 0.0009. Independent likelihood fits are performed in 13 momentum bins of variable size in the range of 2 to 50 GeV/c. A test of the fit with Monte Carlo generated samples verifies that the resulting estimates for the kaon fraction are efficient and unbiased.
The number of charged kaons and pions found in the sample, summed over all momentum bins, is given in Table 1, and the stretch dE/dx distribution of tracks in all momentum bins of the sample is shown in Figure 2. Stretch energy loss, S i , is defined as: where i is the particle hypothesis used to calculate the predicted dE/dx. The normalisation of the predicted distributions of the kaons, pions, muons, and electrons in Figure 2 is obtained from the results of the likelihood fit. The momentum distributions of the charged kaons and pions in the sample, as estimated by the likelihood fit, are shown in Figure 3.

The Branching Ratio Calculation
To determine the τ − → ν τ K − ≥ 0h 0 branching ratio, the number of charged kaons found by the dE/dx likelihood fit in the one-prong candidate sample is corrected for backgrounds and selection inefficiencies: where N τ is the number of tau decay candidates. There were 220652 tau decay candidates recorded in the barrel region of the OPAL detector between the years 1990 and 1995.
f non−τ is the estimated background from non-τ sources in the tau decay candidates (f non−τ = 1.56 ± 0.10%).
N K TOTAL is the total number of kaons in the one-prong tau decay sample, as estimated by the likelihood fit. The number of kaons summed over all momentum bins in the sample is shown in Table 1.
N K bkgnd is the number of background kaons in the one-prong tau decay sample, as estimated by Monte Carlo generated tau decay and hadronic Z 0 → qq samples. The number of background kaons summed over all momentum bins in the sample is shown in Table 1.
0h 0 decays in the tau-pair candidate sample to contribute kaons to the one-prong tau decay sample, and is shown in Table 1. The efficiency has been corrected for biases introduced by the tau-pair selection procedure, and the uncertainty on the efficiency includes the systematic uncertainty arising from this correction. The bias corrections are all consistent with 1 to within one standard deviation of their Monte Carlo statistical uncertainties.
The result of the branching ratio calculation is: where the first uncertainty is statistical and the second is systematic. A summary of the systematic uncertainties, all of which will be described in subsequent sections, is shown in Table 3.
As a cross-check of the general analysis procedure, we also determine the τ − → ν τ π − ≥ 0h 0 branching ratio. Thus, in a similar manner to the τ − → ν τ K − ≥ 0h 0 branching ratio calculation, the total number of charged pions found in the oneprong tau decay sample by the likelihood fit is corrected for backgrounds and the τ − → ν τ π − ≥ 0h 0 efficiency to yield B π − ≥0h 0 = 48.26 ± 0.25%, where the uncertainty is statistical only. The linear correlation coefficient between the OPAL τ − → ν τ K − ≥ 0h 0 and τ − → ν τ π − ≥ 0h 0 branching ratios is approximately −20%. The OPAL τ − → ν τ π − ≥ 0h 0 branching ratio is consistent with the world average, derived from the branching ratios of the tau to the various exclusive final states that include a charged pion listed in reference [17], B π − ≥0h 0 = 48.36 ± 0.23%, and thus we conclude that biases introduced by the analysis procedure are negligible.

dE/dx Systematic Studies
A significant source of systematic uncertainty in the measurement of the inclusive branching ratio is the uncertainty in the parameterisation of the predicted energy loss. To assess this systematic, we determine the sensitivity of the likelihood estimates of the number of kaons within the data sample to the uncertainties in the dE/dx correction factors obtained from the one-prong control samples. The likelihood function shown in Equation 2 is therefore modified such that: n is a vector containing the dE/dx correction factors used in the likelihood fit.
V is the covariance matrix for the dE/dx correction factors, as obtained from the fits to the one-prong control samples.
A 'prior-belief' method is used to assess the systematic uncertainty; in the first iteration of the fit, the correction factors are allowed to vary, and the returned values are found to be consistent with the values obtained from the control samples. In the second iteration, the likelihood fit is repeated, this time keeping the dE/dx correction factors fixed to the values from the first iteration. The systematic uncertainty in N K TOTAL (and N π TOTAL ) due to the parameterisation of dE/dx is then obtained from the square root of the quadrature difference of the fit uncertainties in N K TOTAL (N π TOTAL ) from the two iterations, and is shown in Table 1.

Efficiency Estimation
The efficiencies for kaons and pions from signal events in the tau-pair candidate sample to enter the data sample are estimated using signal events generated with the KORALZ 4.0 Monte Carlo generator and the TAUOLA 2.4 decay package, as described in Section 3. The efficiency for τ − → ν τ K − ≥ 0h 0 decays in the tau-pair candidate sample to contribute to the one-prong sample, ǫ K − ≥0h 0 , is calculated via: where f is the ratio of the exclusive τ − → ν τ K − branching ratio to the inclusive τ − → ν τ K − ≥ 0h 0 branching ratio, as obtained from OPAL data (see Section 6.2). The ǫ K − ≥1h 0 efficiency is calculated from the average of the relevant efficiencies appearing in Table 2, weighted by the associated world average branching ratios. The inclusive and exclusive branching ratios determined by this analysis are all less than 15% correlated with the world average branching ratios used in this efficiency calculation. All efficiencies have been corrected for biases introduced by the tau-pair selection procedure, and, in addition, the Monte Carlo efficiency estimates for the signal events to pass the jet chamber and z-chamber hit requirements are corrected using a control sample of data one-prong tau decays. The uncertainty on each efficiency includes the systematic uncertainties arising from these corrections.
The Monte Carlo modelling of the probability for photons to convert to an e + e − pair in the OPAL detector material has been found by previous studies of data control samples to be correct to within an absolute uncertainty of about ±0.5% [18]. To assess the systematic uncertainty in the branching fraction associated with this modelling, we determine the change in the branching fraction when the efficiency for each exclusive final state is changed by +0.005 for each associated photon in the final state; thus the efficiency for τ − → ν τ K − π 0 decays is changed by 2 × 0.005 = 0.01, and the efficiency for τ − → ν τ K − π 0 π 0 decays is changed by 4 × 0.005 = 0.02. We also determine the change in the branching fraction when the efficiency for each exclusive final state is changed by −0.005 for each associated photon in the final state. The resulting systematic uncertainties in the branching fractions are shown in Table 3.
The variation of the branching ratio due to alternative intermediate resonant structure scenarios for the τ − → ν τ K − π 0 π 0 final state is assessed using efficiency estimates from signal events generated by the modified version of TAUOLA 2.4 through non-resonant production only. The efficiency calculated from this sample for τ − → ν τ K − π 0 π 0 to contribute to the data candidates is within 0.4σ of the central value, where σ represents the combined Monte Carlo statistical uncertainties. Since the two different scenarios produce consistent results for the efficiency, and since there is no a priori reason to expect a strong efficiency dependence on the intermediate resonant structure, no systematic uncertainty for this effect is assigned.
In a similar fashion, the variation of the branching ratio due to alternative intermediate resonant structure scenarios for the τ − → ν τ K − K 0 π 0 final state is explored using events generated through non-resonant production only. The τ − → ν τ K − K 0 π 0 efficiency calculated from these events is within 0.6σ of the central value, thus no additional systematic uncertainty is assigned.
Since the observed momentum distributions for both the pions and the kaons in the sample are well described by the Monte Carlo generated events, no systematic uncertainty is assigned for possible error in either the estimate of the total efficiency, or the momentum dependence of the efficiency estimates, since it is apparent that such errors are small relative to the other systematic uncertainties, as listed in Table 3.

Background Estimation
The number of background kaons in N K TOTAL is estimated by applying the same selection criteria to Monte Carlo samples of hadronic Z 0 → qq decays and three-prong tau decays. From the number of selected events that contain kaons, the estimated number of background kaons is derived, as summarised in Table 1. There are very few background kaons in N K TOTAL , most of which result from the small amount of low-multiplicity hadronic Z 0 → qq events contaminating the tau-pair sample. Due to this small contamination, the systematic uncertainty associated with the background correction is neglible.
The number of background pions in N π TOTAL is estimated in a comparable manner.

Exclusive Branching Ratios
To distinguish between one-prong tau final states with and without accompanying neutral particles, information from the electromagnetic calorimeter is used. Tau decays that include neutral particles tend to have a larger relative electromagnetic energy deposition, E/p, associated with the tau decay cone than tau decays without accompanying neutrals. Such decays also tend to have a larger number of associated clusters, N clus , in the electromagnetic calorimeter than tau decays without neutrals. To ensure accurate modelling of the E/p and N clus distributions of one-prong tau decays to charged kaons, corrections to the Monte Carlo simulation are obtained from a data τ − → ν τ π − ≥ 0h 0 control sample. The following sections describe these corrections, and how they are applied in the determination of the τ − → ν τ K − branching ratio.

Corrections to the Modelling of E/p and N clus
The two-dimensional (p, E/p) and (p, N clus ) distributions of data τ − → ν τ π − ≥ 0h 0 decays are used to determine corrections to the Monte Carlo modelling of the E/p and N clus distributions of τ − → ν τ K − ≥ 0h 0 decays. To this end, the data one-prong tau decay sample is divided into 13 bins of momentum between 2 to 50 GeV/c and 6 bins of E/p, and the number of charged pions and kaons within each bin is determined by maximising the dE/dx likelihood function shown in Equation 2. After background correction, this procedure yields the (p, E/p) distributions of τ − → ν τ π − ≥ 0h 0 and τ − → ν τ K − ≥ 0h 0 decays. The data one-prong sample is also divided into 13 bins of momentum and 4 bins of N clus , and the (p, N clus ) distributions of τ − → ν τ π − ≥ 0h 0 and τ − → ν τ K − ≥ 0h 0 decays are determined in a comparable manner.
The momentum distribution of the charged pions in the one-prong tau decay sample, as shown in Figure 3, appears to be adequately modelled by Monte Carlo generated τ − → ν τ π − ≥ 0h 0 events. However, the E/p and N clus distributions, obtained by summing the (p, E/p) and (p, N clus ) distributions over momentum bins, are not as well modelled by the Monte Carlo, as is shown in Figure 4.
∆N π ij is the statistical uncertainty on N π ij , as obtained from the dE/dx likelihood fit.
B l is the world average branching ratio of tau decay l, taken from reference [17].
ǫ l is the efficiency for events from tau decay l in the tau-pair candidate sample to contribute to the one-prong tau decay sample.
S l ij is the probability, calculated using Monte Carlo generated events, for tau decays of type l in the one-prong tau decay sample to contribute to the (i, j) th bin of the (p, E/p) distribution.
S l j is the correction to the j th bin of the Monte Carlo E/p distribution of tau decays of type l. The S l j are the only parameters allowed to float freely in the fit.
s l i is a normalisation factor used to ensure that the correction to E/p does not affect the shape of the momentum distribution of decays of type l, and is determined from the condition j s l i S l ij S l j = P l i , where P l i is the probability, as calculated from Monte Carlo generated events, for tau decays of type l to contribute to the i th bin of the momentum distribution.
The absolute linear correlation coefficients between N clus and p, and N clus and E/p are less than 20% in all classes of signal. Thus, the N clus and the (p, E/p) distributions are used as discriminators in this analysis, and N clus is treated as a variable uncorrelated to both p and E/p.

The τ − → ν τ K − Branching Ratio
The τ − → ν τ K − branching ratio, B K − , is determined by a χ 2 fit to the (p, E/P ) and N clus distributions of one-prong tau decays with a charged kaon, where the modelling of the distributions is corrected using the results of the previous section. The following χ 2 is minimised: where l is an index referring to either and the indices i, j, and k run over the p, E/p, and N clus bins, respectively. The number of background kaons within each bin is N bkgnd , while ∆N K is the statistical uncertainty on the number of kaons within each bin, N K . Also: and f is the ratio of the OPAL τ − → ν τ K − branching ratio to the OPAL τ − → ν τ K − ≥ 0h 0 inclusive branching ratio, and is allowed to float in the fit.
N K TOTAL is the total number of charged kaons in the one-prong tau decay sample, as determined in Section 5, and shown in Table 1.
ǫ l is the efficiency for events from tau decay l in the tau-pair candidate sample to contribute to the one-prong tau decay sample. The ǫ K − h 0 , ǫ K − ≥2h 0 , and ǫ K − ≥1h 0 efficiencies are calculated from the average of the relevant efficiencies appearing in Table 2, weighted by the associated world average branching ratios. All efficiencies have been corrected for biases introduced by the tau-pair selection procedure.
S l ij is the probability, calculated using Monte Carlo generated events, for tau decays of type l in the one-prong tau sample to contribute to the (i, j) th bin of the (p, E/p) distribution.
S π − j , S π − π 0 j , and S π − π 0 π 0 j are the corrections to the j th bin of the Monte Carlo E/p distribution of the various one-prong tau decays to pions, as determined from the χ 2 fit using Equation 5.
s l i is a normalisation factor used to ensure that the correction to E/p does not affect the shape of the momentum distribution of decays of type l.
T l k is the probability, calculated using Monte Carlo generated events, for tau decays of type l in the one-prong tau sample to contribute to the k th bin of the N clus distribution.
T π − j , T π − π 0 j , and T π − π 0 π 0 j are the corrections to the k th bin of the Monte Carlo N clus distribution of the various one-prong tau decays to pions, as determined in the previous section.
The value of f derived from the result of the fit using Equation 6 is f = 0.436 ± 0.013, where the uncertainty is statistical only. The τ − → ν τ K − branching ratio derived from this result is: where the first uncertainty is statistical and the second is systematic, and the exclusive branching ratio is 60% correlated with the OPAL inclusive τ − → ν τ K − ≥ 0h 0 branching ratio. A summary of the systematic uncertainties on B K − is shown in Table 3.
The statistical uncertainty in B K − due to the statistical uncertainty in the number of kaons within each bin of the (p, E/p) and N clus distributions is assessed by repeating the fit many times, each time varying the central value of the number of kaons within a particular bin, a bin at a time, by plus, then minus, one statistical standard deviation. The square root of half of the quadrature sum of the net variations produced in the branching ratio is taken as the statistical uncertainty.
The total χ 2 per degree of freedom of the fit using Equation 6 is 63.0/50, where 56.6/46 results from the χ 2 per degree of freedom between the data and predicted (p, E/p) distributions. Figure 5 shows the E/p and N clus distributions of one-prong tau decays with a charged kaon, where the normalisation of the contributions to the predicted distribution comes from the OPAL τ − → ν τ K − and τ − → ν τ K − ≥ 0h 0 branching ratios. If the E/p and N clus corrections are not applied, the total χ 2 per degree of freedom from the fit is 79.9/50.
To check for undue variation that may be produced by the choice of binning used in the fit, the above procedure is repeated using various binning schemes. The RMS variation of the results returned by the different fits is less than the statistical uncertainty from the original fit. In addition, a test of the entire fitting procedure (including the calorimeter variable correction procedure described in Sections 6.1 and 6.2) is performed with Monte Carlo generated samples, and verifies that the resulting estimates of f are efficient and unbiased.

dE/dx Systematic Error
The uncertainty in the N K within each N clus or (p, E/p) bin due to uncertainty in the dE/dx correction factors is determined as in Section 5.2. The dE/dx uncertainties are highly correlated between the bins. Thus the uncertainty in B K − due to the dE/dx systematic uncertainty in the number of kaons within each bin of the (p, E/p) and N clus distributions is assessed by repeating the fit many times, each time varying the the central value of the number of kaons within a particular bin, a bin at a time, by plus, then minus, one standard deviation of the dE/dx systematic uncertainty. Half of the linear sum of the absolute variations produced in the branching ratio is taken as the dE/dx systematic uncertainty in B K − , and is shown in Table 3.

Systematic Uncertainty Associated with the Efficiency Estimation
The uncertainty in the τ − → ν τ K − branching ratio due to the limited Monte Carlo statistics used to determine the ǫ K − , ǫ K − h 0 , and ǫ K − ≥2h 0 efficiencies is assessed via the prior-belief method by modifying Equation 6 such that: where ǫ ′ k and ∆ǫ ′ k are the efficiency and associated uncertainty, due to limited Monte Carlo statistics, for kaons from signal channel k in the tau-pair candidate sample to enter the data sample.
ǫ k is the efficiency used in the fit for kaons from signal channel k to contribute to the data sample.
In the first iteration of the fit the ǫ k are allowed to float. In the second, they are fixed to the values from the first iteration. The square root of the quadrature difference in the errors in the branching ratio is taken as the systematic uncertainty associated with the limited Monte Carlo statistics used in the assessment of the efficiencies, and is shown in Table 3.
Various world average exclusive branching ratios are used in the calculation of ǫ K − h 0 and ǫ K − ≥2h 0 (for instance, ǫ K − h 0 is the average of ǫ K − π 0 and ǫ K − K 0 , weighted by the world average branching ratios B K − π 0 and B K − K 0 ). In addition, these same branching ratios are used to assess the inclusive τ − → ν τ K − h 0 and τ − → ν τ K − ≥ 2h 0 branching ratios appearing in Equation 6. The uncertainty in the τ − → ν τ K − branching ratio due to the uncertainties in these world average branching ratios is assessed by varying the branching ratios, one at a time, by plus, then minus, one standard deviation. The square root of half of the quadrature sum of the net variations produced in the branching ratio is taken as the systematic uncertainty, and is shown in Table 3.
The systematic uncertainty in the efficiency and branching fraction due to the uncertainty associated with the Monte Carlo modelling of the probability for photons to convert to an e + e − pair in the OPAL detector material is assessed as described in Section 5.3.

Systematic Uncertainty Associated with the Modelling of E/p and N clus
The uncertainty in the τ − → ν τ K − branching ratio due to the statistical uncertainties in the E/p correction factors is assessed using the prior-belief method by modifying Equation 6 such that: where n ′ is a vector containing the central values of the E/p corrections, as obtained from the fit to the τ − → ν τ π − ≤ 2π 0 sample.
n is a vector containing the values of the E/p corrections used in the χ 2 fit.
V is the covariance matrix for the S correction factors, as obtained from the fit to the τ − → ν τ π − ≤ 2π 0 (p, E/p) and N clus distributions.
In the first iteration, all bins of the three sets of S corrections are allowed to vary in the fit, and the returned values are found to be consistent with the input values. In the second iteration, the fit is repeated, this time keeping the E/p corrections fixed to the values from the first iteration. The systematic uncertainty in B K − due to the modelling of E/p is then obtained from the square root of the quadrature difference of the fit uncertainties from the two iterations.
The uncertainty in the τ − → ν τ K − branching ratio due to the statistical uncertainties in the N clus correction factors is assessed in an analogous fashion. The combined systematic uncertainty in B K − due to the corrections to the modelling of E/p and N clus is shown in Table 3.

Summary and Discussion
From an analysis of the ionisation energy loss of charged particles selected from a sample of 220652 tau decay candidates recorded in the barrel region of the OPAL detector at e + e − centre-of-mass energies close to the Z 0 resonance, we determine the branching ratios: where the first uncertainties are statistical and the second are systematic. The τ − → ν τ K − and τ − → ν τ K − ≥ 1h 0 branching ratios are 60% and 80% correlated to the τ − → ν τ K − ≥ 0h 0 branching ratio, respectively.
Both the inclusive and exclusive branching ratios are in agreement with previous measurements, shown in Table 5, and both contribute significantly to the reduction of the uncertainty on the world averages of the respective quantities.
The theoretical prediction for the τ − → ν τ K − branching ratio is calculated via [20]: where the values of f K = 0.1598 ± 0.0015 GeV and |V us | = 0.2196 ± 0.0023 are both taken from reference [17]. The OPAL measurement is in agreement with this prediction to within 1.3 standard deviations of the combined uncertainties.
The OPAL measurement of the τ − → ν τ K − branching fraction may be used to determine either f K or |V us | through Equation 7. Assuming f K = 0.1598 ± 0.0015 GeV, a value of |V us | = 0.2121 ± 0.0063 is obtained. Assuming a value of |V us | = 0.2196 ± 0.0023, a value of f K = 0.1544 ± 0.0046 GeV is obtained. In both cases, the uncertainty on the determined quantity is dominated by the uncertainty on the τ − → ν τ K − branching fraction.
The theoretical prediction for the τ − → ν τ K − ≥ 0h 0 branching ratio is obtained by summing the theoretical predictions for the five states appearing in Table 4: The OPAL measurement is in agreement with this prediction to within 0.7 standard deviation of the combined uncertainties.
N τ 220652 f non−τ 1.56 ± 0.10% # in one-prong sample 143528 N π TOTAL 75566 ± 389 ± 246 N π bkgnd 89 ± 16 ǫ π − ≥0h 0 0.720 ± 0.001 ± 0.001 0.759 ± 0.006 ± 0.007 Table 1: Pion and kaon composition of the one-prong tau decay candidate sample, as estimated by the likelihood fit to the measured dE/dx of tracks in the sample. The first uncertainty on N π TOTAL and N K TOTAL is the statistical uncertainty from the fit, and the second is due to the systematic uncertainty in the dE/dx correction factors. Also shown are the estimated backgrounds in N π TOTAL and N K TOTAL , along with the efficiencies for signal events to contribute to the one-prong tau decay sample. The first uncertainty on each efficiency estimate is due to the limited statistics of the Monte Carlo generated samples, while the second is due to the uncertainty in the various branching ratios used to calculate the efficiency, as listed in Table 2.