Measurements of BR (b ->tau- nubar_tau X) and BR (b ->tau- nubar_tau D*+- X) and Upper Limits on BR (B- ->tau- nubar_tau) and BR (b ->s nu nubar)

Inclusive branching ratios involving b to tau transitions are measured in approximately four million hadronic Z decays collected by the ALEPH detector at LEP. The fully-inclusive branching ratio b ->tau nu X and the semi-inclusive branching ratio b ->tau nu D*+/- X are measured to be (2.43 +/- 0.20 +/- 0.25)% and (0.88 +/- 0.31 +/- 0.28)%, in agreement with the standard model predictions. Upper limits on the branching fractions b ->tau nu and b ->s nu nubar are set to 8.3 10**-4 and 6.4 10**-4 at the 90% C.L. These results allow a 90% C.L. lower limit of 0.40 (GeV/c**2)**-1 to be set on the tan(beta)/mH+/- ratio, in the framework of type-II two-Higgs-doublet mode


Introduction
Third-generation fermions are involved in both the initial and final states of b → τ transitions.A measurement of branching ratios pertaining to these transitions can be compared to the standard model predictions, yielding direct tests of heavy fermion interactions.These interactions are especially sensitive to the mechanism underlying the electroweak symmetry breaking, i.e., to the origin of mass, and allow extensions of the standard model to be constrained in this respect.In the study presented here, the branching fractions of the following processes are considered.
• The inclusive branching fraction BR(b → τ − ντ X) can be compared to the standard model prediction of 2.30±0.25%[1], as computed in the framework of the Heavy Quark Effective Theory (HQET).The standard model transition, illustrated in Fig. 1a, could be modified by the exchange of a new charged boson, as shown in Fig. 1b with a charged Higgs boson.In any type-II two-Higgs-doublet model, such as the minimal supersymmetric extension of the standard model, the corresponding contribution to the transition amplitude is proportional to (tan β/m H ± ) 2 [2,3,4], where tan β is the ratio of the vacuum expectation values of the two Higgs doublets, and m H ± is the mass of the charged Higgs boson.
• A similar test can be performed with the semi-inclusive decay b → τ − ντ D * + X.In the standard model, this branching fraction is expected to be approximately 1% [5].
The proportion of D * relative to D final states would be reduced in the presence of a charged Higgs boson [6].
The measurements of the branching fractions of these four processes, b → τ − ντ X, b → τ − ντ D * ± X, B − → τ − ντ and b → sν ν, are performed with all data collected by the ALEPH detector between 1991 and 1995, at centre-of-mass energies at and around the Z resonance, corresponding to approximately four million hadronic Z decays.As described in previous ALEPH studies applied to smaller data samples [10,11], a single technique is used to extract the four branching ratios.The initial-state b quark is identified by means of standard b-tagging algorithms, and the final-state τ lepton (or the two neutrinos in the case of b → sν ν) is identified by the missing energy carried away by the neutrinos.
The energy-flow and b-tagging algorithms are briefly described in Section 2 along with properties of the ALEPH detector relevant for the present study.The selection algorithms aimed at rejecting background events with a large missing energy are presented 1 in Section 3. The method to estimate the residual background yields is detailed in Section 4. The branching ratios BR(b → τ − ντ X) and BR(b → τ − ντ D * ± X) are determined in Section 5, and the limits on BR(B − → τ − ντ ) and BR(b → sν ν) are extracted in Section 6.Finally, an alternative measurement of BR(b → τ − ντ X) with opposite-sign di-leptons in the final state is presented in Section 7. The results are interpreted in the framework of type-II two-Higgs-doublet models in Section 8 and summarized in Section 9.

The ALEPH Detector
A detailed description of the ALEPH detector can be found in Ref. [12], and of its performance in Ref. [13].Charged particles are detected in the central part, consisting of a precision silicon vertex detector, a cylindrical drift chamber and a large time projection chamber, measuring altogether up to 31 space points along the charged particle trajectories.A 1.5 T axial magnetic field is provided by a superconducting solenoidal coil.Charged particle transverse momenta are reconstructed with a 1/p T resolution of (6 × 10 −4 5 × 10 −3 /p T ) (GeV/c) −1 .In the following, good tracks are defined as charged particle tracks reconstructed with at least four hits in the time projection chamber, originating from within a cylinder of length 20 cm and radius 2 cm coaxial with the beam and centred at the nominal collision point, and with a polar angle with respect to the beam such that | cos θ| < 0.95.
Jets originating from b quarks are identified with a lifetime b-tagging algorithm [14], which takes advantage of the three-dimensional impact parameter resolution of charged particle tracks.For tracks with two space points in the silicon vertex detector (i.e., | cos θ| < 0.7), this resolution can be parametrized as (25 + 95/p) µm, with p in GeV/c.
In addition to its rôle as a tracking device, the time projection chamber also measures the specific energy loss by ionization dE/dx.It allows low momentum electrons to be separated from other charged particle species by more than three standard deviations.
Electrons (and photons) are also identified by the characteristic longitudinal and transverse developments of the associated showers in the electromagnetic calorimeter, a 22 radiation length thick sandwich of lead planes and proportional wire chambers with fine read-out segmentation.The relative energy resolution achieved is 0.18/ √ E (E in GeV) for isolated electrons and photons.
Photon conversions to e + e − in the detector material are identified as a pair of oppositely charged identified electrons satisfying the following conditions: (i) the distance of closest approach to the beam of the two reconstructed tracks is larger than 2 mm; (ii) the distance between the two tracks at their point of closest approach is smaller than 2 cm in space; (iii) the position of this point is consistent with a material boundary; and (iv) the invariant mass is smaller than 20 MeV/c 2 , when calculated as for an e + e − pair coming from this point of closest approach.
Muons are identified by their characteristic penetration pattern in the hadron calorimeter, a 1.5 m thick yoke interleaved with 23 layers of streamer tubes, together with two surrounding double-layers of muon chambers.In association with the electromagnetic calorimeter, the hadron calorimeter also provides a measurement of the hadronic energy with a relative resolution of 0.85/ √ E (E in GeV).
Taus are identified by the missing energy carried away by their decay neutrinos.The total visible energy is measured with an energy-flow reconstruction algorithm which combines all the above measurements, supplemented by the energy detected at low polar angle (down to 24 mrad from the beam axis) by two additional electromagnetic calorimeters, used for the luminosity determination.The relative resolution on the total visible energy varies between 0.60/ √ E for high multiplicity final states and 0.25/ √ E for final states of low multiplicity without neutral hadrons.In addition to the total visible-energy measurement, the energy-flow reconstruction algorithm also provides a list of reconstructed objects, classified as charged particles, photons and neutral hadrons, and called energy-flow particles in the following.

Event selections
A b → τ transition followed by the decay of the τ always produces two energetic ν τ 's.Such transitions can therefore be identified in e + e − → b b events on the basis of a large measured missing energy.The argument holds as well for b → sν ν decays, which can be selected along the same lines as b → τ transitions.In this section, the algorithms aimed at rejecting background events leading to large measured missing energy are described.An alternative selection, based on the presence of two identified leptons with opposite electric charge, is described in Section 7.

Preselection
Approximately four million hadronic Z decays are selected in the data collected between 1991 and 1995 at energies at and around the Z resonance when the standard criteria (at least five good tracks carrying at least 10% of the centre-of-mass energy [15]) are applied.To keep only two-jet events well contained in the detector acceptance, the polar angle is required to satisfy | cos θ thrust | < 0.7 to match the acceptance of the vertex detector, and the thrust value must exceed 0.85.
Each event is then divided in two hemispheres with respect to the plane perpendicular to the thrust axis.In each hemisphere, the missing energy E 1,2  miss is defined as the difference between the expected true energy E 1,2  true and the measured visible energy E 1,2 vis .The latter is determined from the total energy of all energy-flow particles contained in that hemisphere, while the former is estimated from the centre-of-mass energy √ s and with energy-momentum conservation: E 1,2 true = (s + m 2 1,2 − m 2 2,1 )/2 √ s, where m 1 and m 2 are the measured invariant masses in the two hemispheres [16].
The main background to the final states searched for in this analysis consists of Z → b b or cc events followed by a semi-leptonic b or c decay into an electron or a muon (hereafter called lepton and denoted ℓ), with a large missing energy carried away by the neutrino ν ℓ .This background can be considerably reduced by rejecting hemispheres in which a lepton is identified.The standard lepton identification [17] is not used here, but is replaced by much looser criteria.A good track is tagged as an electron either if the dE/dx is compatible with that of an electron and incompatible with that of a pion, or if the transverse and longitudinal profiles of the associated electromagnetic shower are compatible with those of an electron.Electrons and positrons originating from identified photon conversions are not considered in this process.Similarly, a good track is tagged as a muon if it is associated to a few hits either in the last layers of the hadron calorimeter or in the two layers of muon chambers.A lepton identification efficiency over 95% is achieved for electrons (muons) with a momentum in excess of 1 (2) GeV/c.In addition, to reduce the contamination from Z decays into lighter quark pairs (uū, d d, ss and cc) in which a large missing energy is faked due to finite detector resolution effects, the final sample is enriched in Z → b b events with b tagging.The corresponding criteria are dependent on the signal final state and are described in the following subsections.

The b → τ − ντ X final state
For the b → τ − ντ X final state selection, a hemisphere is kept only if the opposite hemisphere is tagged as arising from a b quark.Each good track in the opposite hemisphere is assigned a probability of originating from the primary interaction point, on the basis of its impact parameter significance.The confidence level α hemi that all good tracks in that hemisphere come from the primary interaction point, determined under the assumption that the individual probabilities are uncorrelated [14], is required to be smaller than 1%.
Residual backgrounds like Z → τ + τ − decays, two-photon processes or beam-gas interactions are likely to yield a large missing energy in the final state and might therefore bias the analysis.These background sources are reduced down to negligible levels by requiring at least seven good tracks and a total missing energy smaller than 50 GeV, with almost no additional loss of the signal.

The B − → τ − ντ and b → sν ν final states
The same b-tagging technique as for the b → τ − ντ X selection is applied but the residual backgrounds are rejected in a slightly different manner.The opposite hemisphere is required to contain at least six good tracks, and its missing energy must be smaller than 25 GeV.This latter cut reduces the effect on the missing-energy distribution in the signal hemisphere, which is expected to extend towards larger values than in the previous case.

The b
The b → τ − ντ D * ± X final state selection is primarily based on an exclusive reconstruction of a D * ± , very similar to that performed in Ref. [18] although with slightly tighter cuts as a result of the larger statistics accumulated until 1995.
Candidate D * ± 's are searched for in the decay channel D * ± → D 0 π soft , followed by Neutral pions are identified as photon pairs with an invariant mass compatible with the π 0 mass and a χ 2 probability in excess of 1%.Candidate K 0 S 's are identified as pairs of oppositely charged particle tracks (assumed to be charged pions) forming a secondary vertex at least 2 cm away from the primary interaction point, and with a reconstructed invariant mass within 5 MeV/c 2 of the K 0 S mass.Candidate K ± 's are identified as good tracks of momentum in excess of 3 GeV/c and with a dE/dx compatible with that of a charged kaon within two standard deviations.In addition the angle θ * between the kaon momentum direction, evaluated in the D 0 rest frame, and the D 0 line-of-flight must satisfy | cos θ * | < 0.9.
The reconstructed invariant mass of the candidate D 0 's is required to be within 15 MeV/c 2 (50 MeV/c 2 ) of the D 0 mass for the first (last) two decay channels.The ratio of the D 0 energy to the beam energy is required to be between 0.25 and 0.50.In addition, the reconstructed D 0 vertex must be incompatible with the primary interaction point by at least two standard deviations.
Each reconstructed D 0 is associated with a low momentum pion (p π soft < 4.2 GeV/c) to form a D * ± .The reconstructed invariant mass difference m D * ± − m D 0 is required to be within 20 MeV/c 2 of the nominal mass difference.
Finally, the D * ± hemispheres are tagged as arising from a b quark by requiring the presence of an additional good track in a cone of half-angle 30 • around the D 0 momentum direction.This track must be incompatible with originating from the primary interaction point by more than three standard deviations, and its electric charge must be opposite to that of the π soft .No b tagging needs therefore be applied to the opposite hemisphere.

Event yields
The expected event yields were determined with the following samples of events, simulated with JETSET [19,20], and subsequently processed through the full GEANT [21] simulation of the detector: • approximately four million hadronic Z decays for the background estimate; • over 40 000 b b events, each containing a b → τ − ντ X, with the τ polarization determined in the limit of the free-quark model [1] for m b = 4.8 GeV/c 2 and m c = 1.4 GeV/c 2 ; • about 20 000 b b events with at least one B − → τ − ντ and 100% polarized τ 's ; • almost 20 000 b b events with at least one b → sν ν, generated with the B → X s ν ν decay probability given in Ref. [8], in which the new-physics parameters were conservatively chosen to minimize the expected selection efficiency.The model was improved with a realistic X s invariant mass distribution, described as a Gaussian of mean and variance 1.35 and 0.60 GeV/c 2 [9], supplemented by two peaks at the K and K * masses [22] with branching fractions of 0.06 ± 0.02 and 0.29 ± 0.01.In all these samples, the polarization of the τ 's was properly taken into account in the simulation of their decay kinematics [23], and events containing a Λ b were reweighted according to the measured Λ b polarization [24], so that the correct distribution of the ν τ momenta, and therefore the missing energy, be obtained.
The numbers of hemispheres expected to be selected from the various background and signal processes and the numbers of hemispheres selected in the data are listed in Table 1.
The measurement of the various branching fractions relies on the observation of an excess of events at large E miss in the missing-energy distribution of the data with respect to that of the simulated background.Given the limited number of events expected from the various signal processes (Table 1), it is essential to have a detailed understanding of all uncertainties affecting the background distribution both in shape and normalization (Section 4) and of the signal selection efficiencies (Section 5).

Background estimate
The missing-energy distribution of the background to b → τ and b → sν ν transitions depends on detector performance in three areas, (i) the visible-energy reconstruction for hemispheres with a priori no missing energy such as Z hadronic decays into light-quark pairs; (ii) the effectiveness of the lepton veto to reject background from b, c → ℓν ℓ X; and (iii) the performance of the b-tagging.To minimize the influence of possible inaccuracies in the detector simulation, these three quantities are derived from the data and subsequently used to re-weight the simulated events.

The visible-energy calibration
As mentioned above, the accuracy of the visible-energy reconstruction can be studied with a sample of events in which no significant missing energy is expected.Such a sample is selected in a way identical to that followed in Sections 3.1 and 3.2 except that the btagging criterion is replaced by a b-rejection one requiring α hemi > 0.5.It yields fractions of 88.7%, 10.5% and 0.8% for (uū,d d,ss), cc and b b events, respectively.
The visible energy E vis is the sum of the energy E cha of the charged particles reconstructed in the central tracker, the energy E pho of the photons detected in the electromagnetic calorimeter, and the energy E neu of the neutral hadrons, the reconstruction of which makes use of all identification and energy measurement capabilities of the detector (Section 2).These three contributions to the visible energy may therefore be affected by different systematic uncertainties and must be calibrated separately.The hemispheres were thus binned in a grid according to the fractional contributions to the visible energy of the charged, photon and neutral hadron components (E cha /E vis , E pho /E vis , E neu /E vis ).The accuracy of the detector simulation was examined by comparing the observed and simulated missing energy in each of the bins.
A noticeable inaccuracy of the simulation was found in hemispheres with a large proportion of neutral hadronic energy.This disagreement with the data could lead to a large systematic uncertainty in the final result.Moreover, the neutral hadronic energy is expected to be smaller in b b events than in other hadronic Z decays.It was therefore decided to reject all hemispheres for which E neu is in excess of 7 GeV, both in this calibration procedure and in the selections described in the previous section, preserving 69.6% of the events in the data and 68.7% in the simulation.
The residual differences between the data and the simulation were corrected by scaling E cha , E pho and E neu by f cha , f pho and f neu , respectively, in the simulated events.The values of these calibration factors were obtained by minimizing the following χ 2 , where the mean energy values are computed in each bin i, E data vis,i is the total visible energy measured in that bin, the E MC xxx,i values are the charged, photonic and neutral hadronic energies expected in that bin, and σ i is the uncertainty due to the limited statistics of the event samples.The result of the fit is shown in Table 2.While f cha and f pho are compatible with unity, a sizeable calibration is found to be still necessary for the neutral hadronic energy, despite the cut at 7 GeV.This effect, already reported in Ref. [10], results from an inadequacy of the simulation of nuclear interactions.f cha 1.002 ± 0.001 f pho 1.000 ± 0.002 As another consequence, the E neu resolution is better in the simulation than in the data.In addition to the rescaling, the simulated neutral hadronic energy was therefore smeared on an event-by-event basis, leading to a (1.1 ± 1.6)% relative deterioration of the E neu resolution.This smearing was performed in order to equalize the observed and expected numbers of events with a reconstructed E miss in excess of 14 GeV.Indeed, in events where no real missing energy is expected, such a large reconstructed E miss value is often due to a loss of neutral hadronic energy.
The distributions of the reconstructed missing energy in the control sample are displayed in Fig. 2, both for data and simulated events, before and after the calibration procedure.The agreement is found to be adequate over the whole spectrum, and in particular at large missing energy where the background to the signal of interest in this paper has to be evaluated.
Finally, due to the specific hadronization and fragmentation of b-quark jets, the sharing between the charged, photonic and neutral hadronic energy may lead to an overall visibleenergy shift different from that observed for lighter quark species.The visible energy of hemispheres tagged as arising from a b quark was therefore rescaled by an overall factor f b , determined from the b-tagged sample of Section 3.2 so as to match the peak positions (which are not affected by events with truly missing energy) of the expected and observed visible-energy distributions.This b-specific calibration factor is found to be compatible with unity (Table 2).

The lepton rejection effectiveness
In order to estimate the effectiveness of the lepton veto, aimed at rejecting Z → b b or cc events followed by a semi-leptonic b or c decay into an electron or a muon, the veto was applied to unbiased, pure, lepton data and simulated samples.
A 97.6% pure electron sample was obtained by selecting hadronic events with a photon conversion, as described in Section 2, in which only one of the two particles has to be identified as an electron.Photon conversions are well suited to study the electron veto since the hadronic environment is very similar to that of semi-leptonic heavy-flavour decays.The electron-veto effectiveness is thus given by the probability, after background subtraction, that the second particle be identified as an electron by the lepton-veto criteria.However, this probability depends on (i) the electron momentum, the spectrum of which is different in photon conversions and in semi-leptonic b/c decays; and (ii) the number of wires used for the dE/dx measurement, which tends to be smaller in photon conversions than in semi-leptonic b/c decays due to shared hits with the first electron of the converted pair.Therefore, the electron identification probability was mapped as a function of the electron momentum and the number of hits associated to the corresponding track, both for the data and the simulated events with photon conversions.The ratio of these two maps was then used to re-weight the simulated b/c events with semi-leptonic decays.
For the muons, a sample originating from the γγ → µ + µ − process was selected.The absence of hadronic environment is less important because the identification is achieved with the external layers of the hadron calorimeter and the muon chambers.A 97.1% pure sample of low-momentum muons was obtained by requiring exactly two good tracks (accompanied by no neutral energy-flow particle) with opposite electric charge, one of which being identified as a muon, a total momentum less than 30 GeV/c and an invariant mass smaller than 2 GeV/c 2 .Above 3 GeV/c, muons are expected to traverse entirely the hadron calorimeter and the muon chambers, making the identification efficiency independent of momentum.The γγ → µ + µ − sample was therefore supplemented by Z decays into muon pairs to determine the latter.As above, and after background subtraction, the muon-veto effectiveness was given by the probability that the second particle be identified as a muon by the lepton-veto criteria as a function of momentum, both in the data and in the simulation.The ratio was then used to re-weight the simulated b/c events with semi-leptonic decays.
The identification probability for electrons and muons, determined from these samples, is displayed as a function of the lepton momentum in Figs.3a and 3b, for data and simulation, and is compared to the identification efficiency for simulated b/c events with semi-leptonic decays.The data-to-simulation ratio, applied as a correction factor to the latter, is consistently found to be close to unity.

The b-tagging efficiency
The efficiency of the b-tagging criterion can be determined directly from the data using the double tag method [14], so as to minimize the systematic uncertainties related to the limited knowledge of b-hadron production and decay.The small correlation between the b-tagging probabilities in the two hemispheres is negligible for the present purpose.
From the knowledge of the fractions R b and R c of hadronic Z decays into b b and cc, the efficiencies ε b and ε c of the b-tagging criterion on b and c hemispheres can be measured in the data by comparing the fraction f 1 of hemispheres which pass a given α hemi cut with the fraction f 2 of events in which both hemispheres pass the same cut.These fractions are given by where ε x is the hemisphere tagging efficiency for light-quark hemispheres, the value of which is determined directly with the data as explained in Ref. [14].The difference between the observed and expected ε x values for a given α hemi cut arises from from different impact parameter resolutions in the data and in the simulation.This difference was estimated from the fraction of hemispheres satisfying α hemi neg < 0.01, where α hemi neg is the hemisphere probability recomputed with only those tracks having negative impact parameter (and reversing the impact parameter sign).Such tracks nearly always originate from the primary vertex, and account for resolution differences.The relatively small contribution of long-lived particles (such as K 0 S or Λ 0 ) is taken from the simulation.With the cut α hemi < 0.01, ε x was found to be 2.46%, i.e., slightly larger in the data than predicted by the simulation by a factor of 1.17 ± 0.11.With the ε x value obtained this way, Eq. ( 2) was solved for ε b and ε c .The tagging efficiencies for b (56.7%) and c (14.8%) hemispheres were also found to be larger in the data by factors of 1.03 ± 0.01 and 1.03 ± 0.03, respectively, and were used to re-weight the simulated events accordingly.
The b → τ − ντ X and b → τ − ντ D * ± X branching fractions can be determined by evaluating the excess of events over the background expected at large missing energy from the two selections (Sections 3.2 and 3.4).The E miss distributions in the b → τ − ντ X and b → τ − ντ D * ± X final states are displayed in Figs.4a and 4b, for the data and for the re-weighted simulation (Section 4), normalized to the number of observed events.

Branching ratio extraction
An excess of events at large E miss is indeed observed over the background in both distributions.The corresponding branching ratios were extracted with a binned-likelihood fit of the expected to the observed missing-energy distributions, keeping the normalization to the number of events observed.The fit was performed in E miss intervals chosen so as to (i) minimize their total uncertainty, defined as the quadratic sum of the statistical and the systematic contributions (Section 5.2); and (ii) make the measurements statistically independent of that of the B − → τ − ντ branching fraction.The optimal intervals, indicated in Fig. 4, are found to be 18-35 GeV and 12-35 GeV, respectively.
The numbers of events in these intervals, observed in the data and expected from both the background and the signal processes, are displayed in Table 3.In particular, the cascade decay b → D − s X with D − s → τ − ντ , which yields an E miss spectrum similar to that of the decay b → τ − ντ X, is included in the backgrounds to this signature.These numbers yield measured branching fractions of BR(b → τ − ντ X) = [2.43 ± 0.20(stat.)± 0.25(syst.)]%, and BR(b → τ − ντ D * ± X) = [0.88 ± 0.31(stat.)± 0.28(syst.)]%, in agreement with the standard model predictions of (2.30 ± 0.25)% and approximately 1%, respectively.The systematic uncertainties, also indicated above, are discussed in the next section.Table 3: Numbers of events observed in the data at large E miss in the b → τ − ντ X and b → τ − ντ D * ± X final states.Also indicated are the events expected from the various background and signal processes.For the latter, the fitted branching ratio value was assumed.The corresponding selection efficiencies are indicated in brackets.

Systematic uncertainties
The measured branching ratios BR(b → τ − ντ X) and BR(b → τ − ντ D * ± X) can be affected by several sources of systematic effects on the signal selection efficiencies and by residual inaccuracies of the background simulation.The relevant sources common to b → τ − ντ X and b → τ − ντ D * ± X are addressed in turn below, and their effects on the branching fractions are summarized in Table 4.
1.The missing-energy distributions of b → τ − ντ X and b, c → ℓν ℓ X depend on the experimentally determined quantities x b , x c , BR(b → ℓν ℓ X) and BR(b → c → ℓν ℓ X).However, the small, residual differences between the E miss distributions in the data and the simulation, obtained by reversing the lepton-veto criteria (Fig. 5), were found to be well accounted for by the measured uncertainties on their values.The whole extraction procedure of the branching fractions was therefore repeated by varying these quantities within their uncertainties, so as to evaluate the corresponding systematic effects.Similarly, the branching ratio dependence on the D * * content in b decays was accounted for according to Ref. [25].
2. The b-fragmentation modelling affects the b-energy spectrum, and therefore also the missing-energy spectrum in b → τ − ντ X and b, c → ℓν ℓ X.The effect of its knowledge was estimated in a similar manner by using the four different fragmentation functions described in Ref. [26], adjusted to reproduce the measured value of x b .The corresponding systematic uncertainty was defined as the largest change in the measured branching fractions.3. The ντ energy spectrum in the b rest frame also depends on whether it is computed with HQET or with the spectator model [1].To reflect this difference, the neutrino spectrum was distorted with a weight depending linearly on the neutrino energy such that the averaged energy changes by ±10%, and the resulting change in the measured branching fractions noted.
7. The visible-energy calibration procedure (Section 4.1) is tuned with light-quark events, and may not be entirely accurate for b b events in the signal region.In this region, the missing energy is found in the simulation to mainly originate from misreconstructed neutral hadrons.It is observed that b b events yield more hemispheres with a large neutral hadronic energy in the data than in the simulation.This excess is expected to increase the residual background in the signal region by ∼ 20%, treated as an additional systematic uncertainty.
Two additional systematic effects, specific to the b → τ − ντ D * ± X final state were also identified.First, the combinatorial background was estimated in that case from the sidebands around the D 0 peak, and used to derive the corresponding uncertainty on BR(b → τ − ντ D * ± X).Second, the b-tagging criterion specific to this channel (Section 3.4) yields a different efficiency in data and simulation, which affects BR(b → τ − ντ D * ± X).
The missing-energy distribution of the events selected as described in Section 3.3 is displayed in Fig. 6, and is compared to that of the background.Also indicated in Fig. 6 is the expected enhancement at large missing energy, should either the B − → τ − ντ or the b → sν ν branching ratio be 1%.No background subtraction was performed, making the validity of the results presented here unaffected by possible systematic uncertainties related to the knowledge of the residual background, largely suppressed by a lower cut on the missing energy at 35 GeV.The latter cut was optimized so as to maximize the expected 90% C.L. upper limit, evaluated with simple event counting and in the absence of new physics [30], on the B − → τ − ντ and b → sν ν branching ratios.The numbers of events, observed in the data and expected from background and signal, are displayed in Table 5 in three missing-energy intervals.Two events with a missing energy in excess of 35 GeV were observed, with 2.5 ± 1.6 events expected from all background processes.In absence of any systematic uncertainty, the 90% C.L. upper limits on the B − → τ − ντ and b → sν ν branching fractions are found to be 8.1 × 10 −4 and 6.2 × 10 −4 , respectively.
However, these limits are affected by the uncertainty on the expected fraction of B − → τ − ντ and b → sν ν events with such a large missing energy or, almost equivalently, with such a large value of x b .(Most of these events are characterized by a value of x b in excess of 0.9.)This fraction was determined [31] to be 0.146 +0.025 −0.021 , which translates Table 5: Numbers of events observed in the data in three different E miss intervals.Also indicated are the events expected from the various background and signal processes.For the latter, a branching ratio of 1% was assumed.
E miss interval [30,35] GeV [35,40]  to an uncertainty of 15% on the number of events expected.In addition, the fraction of B − in Z → b b events is known [32] to be (38.9± 1.3)%.The latter uncertainty of 3% on the number of events expected affects only the limit on the B − → τ − ντ branching ratio.These two uncertainties were taken into account following the method of Ref. [33], yielding the limits at the 90% confidence level.
7 Measurement of BR(b → τ − ντ X) with di-leptons An alternative method of measuring BR(b → τ − ντ X) was developed with events where both the τ and the accompanying D decay to e or µ.Hence, the signature used to tag the signal events is a pair of leptons (e, µ) of opposite sign in a jet.The background, originating from b → cℓ − νℓ followed by c → qℓ + ν ℓ , is about 20 times larger than the signal.Signal and background are therefore separated on the basis of their different kinematic properties.Although this method is statistically less powerful than that based on missing energy, it represents an interesting cross-check since it is based on a complementary sample of events and sources of systematic uncertainties are largely different.Indeed, the main contribution to the systematic error comes from the uncertainty on the product BR(b → ℓν ℓ c, c → ℓν ℓ q) and on the double charm decay rates, B → D s D(X) and B → D 0 D(X).
The presence of three neutrinos in the decay chain of signal events is the main difference between signal and background.As a consequence, b → τ − ντ X decays give larger missing energy, softer lepton spectrum and smaller charged multiplicity to the jet containing the lepton candidates.These different kinematic properties of the various categories of events are used to separate the signal from the background.A multivariate analysis technique with a multilayered neural network (NN) is used to obtain the best discriminating power.

Event selection
Electron and muon identification follows the standard criteria [17], with two refinements for muons.Firstly, the muon momentum cut is lowered from 3 GeV/c to 2.2 GeV/c to increase the acceptance for the signal (2.2 GeV/c is the minimum momentum for a muon to reach the muon chambers) and secondly, any track "shadowed" by another track is rejected in the di-lepton selection.Two tracks are said to be shadowing if they share in common either hits in the last ten layers of the hadron calorimeter or at least one three-dimensional muon chamber hit.This latter cut allows a good description of the background to be obtained in the sample of events with two muons in the same jet, as shown in Ref. [34].To obtain the di-lepton sample used for the training of the neural network, two additional cuts on the invariant mass M ℓℓ of the lepton pair are applied.
• Electron pair invariant masses are required to be greater than 0.3 GeV/c 2 in order to reject γ conversions or Dalitz decays of η, π 0 .
• Electron pairs and muon pairs are required to have an invariant mass smaller than 2.5 GeV/c 2 to avoid pairs coming from J/ψ decays.
After these selection cuts, the simulated di-lepton sample consists of 451 signal events, 7834 b → ℓ background events and 5249 di-leptons coming from hadrons mis-identified as leptons and light-hadron decays.

Analysis method and background estimation
Jets are reconstructed by the JADE clustering algorithm [35] with a y cut of 0.0044 [17].
The following variables are selected as input to the NN: the missing energy in the hemisphere of the jet containing the lepton pair, the momenta of the leptons boosted to the reconstructed b-hadron rest frame [34], the invariant mass of the lepton pair, the transverse momenta of the leptons with respect to the jet direction, the total energy and the charged energy of the jet containing the lepton pair, the fraction of the jet charged energy carried by the leptons and the number of charged particles in the jet.
The main background consists of di-leptons coming from b → cℓ, c → ℓ ′ ν ℓ ′ decays.To determine the product BR(b → ℓν ℓ c)×BR(c → ℓ ′ ν ℓ ′ ), the semi-leptonic branching fraction BR(b → ℓν ℓ X) is fixed to the LEP average value [29], while the b → c → ℓ ′ fraction is determined from the same sample of di-lepton events as follows: a neural network with the input variables described above is trained to separate b → cℓ → ℓℓ ′ X decays from all other processes; and the b → c → ℓ fraction is extracted by fitting the simulated From the dependence between BR(b → τ − ντ X) and r of Eq. 3, and from the value of BR SM (b → τ − ντ X), a 90% C.L. upper limit on r was extracted [3,4,38]: tan β/m H ± < 0.49 (GeV/c 2 ) −1 at 90% CL. ( The τ polarization actually plays here an important rôle since, for r = 0.49 (GeV/c 2 ) −1 , it amounts to −0.28 while the standard model prediction is −0.735 in b → τ − ντ X decays.This effect had been neglected in previous ALEPH analyses [11].
The dependence between BR(B − → τ − ντ ) and r of Eq. 4 and the standard model value of the branching ratio, BR SM (B − → τ − ντ ) = 7.4 × 10 −5 (f B /160 MeV) 2 (|V ub |/0.004) 2 , were interpreted in terms of a 90% C.L. upper limit on r: tan β/m H ± < 0.40 (GeV/c 2 ) −1 at 90% CL. ( In the limit extraction, the systematic uncertainties on f B [39] and on V ub [25] were taken into account following the method of Ref. [33].The combination of the two results in Eqs. 5 and 6 cannot improve on the latter, since the corresponding r value does not enhance the b → τ − ντ X branching ratio.

Figure 1 :
Figure 1: Feynman graphs for the processes b → τ − ντ X, B − → τ − ντ and b → sν ν, in the standard model (a,c,e) and in two-Higgs-doublet models with the exchange of a charged Higgs boson (b,d,f).

Figure 2 :
Figure 2: Reconstructed missing-energy distributions using the light-quark tag and e/µ veto (a) before calibration; (b) after missing-energy recalibration and applying the cut E neu < 7 GeV.

Figure 3 :
Figure 3: Lepton identification efficiencies (a) for electrons and (b) for muons as a function of momentum in the control samples (see text) for the data (open circles) and the simulation (full circles), and their ratio (triangles).The histograms show the identification efficiency for simulated b/c events with semi-leptonic decays.

Figure 4 :
Figure 4: Missing-energy distributions in the b → τ − ντ X (a) and the b → τ − ντ D * ± X (b) final states,for the data (dots) and the simulation (histograms).The latter is subdivided in (i) the fitted signal contribution (empty histogram); (ii) the contribution from b and c semi-leptonic decays (light hatching); and (iii) the residual background (dark hatching).Also indicated are the E miss intervals considered for the branching ratio measurements (Section 5.1).

Figure 5 :
Figure 5: E miss distributions obtained for the data (dots) and the simulation (histogram) from the selection described in Section 3.2 with the lepton-veto criteria reversed.

Figure 6 :
Figure 6: Missing-energy distributions in the B − → τ − ντ and b → sν ν final state selection, for the data (dots) and for the simulated background (full histogram).Also indicated are the contributions of the B − → τ − ντ (dashed histogram) and b → sν ν (dotted histogram) processes with a branching ratio of 1%.The arrow shows the region in which the two limits are calculated.

Figure 7 :Figure 8 :
Figure 7: Data and normalized Monte Carlo distributions for (a) the missing energy of the leptonic hemisphere and (b) the momentum of one lepton boosted to the reconstructed b-hadron rest frame.The normalization of the simulated distribution for signal events corresponds to the BR(b → τ − ντ X) value obtained from the fit.

Table 1 :
Numbers of hemispheres selected in the data in the four different final states.Also indicated are the numbers of hemispheres expected from the various background and signal processes.For the latter, the standard model branching ratio values were assumed.The standard model branching ratios and the selection efficiencies are indicated in brackets.

Table 2 :
Calibration parameters (see text for a detailed description) for E vis .