Combined results on b-hadron production rates, lifetimes, oscillations and semileptonic decays

Combined results on b-hadron lifetimes, b-hadron production rates, B^0_d - \bar{B^0_d} and B^0_s -\bar{B^0_s} oscillations, the decay width difference between the mass eigenstates of the B^0_s -\bar{B^0_s} system, and the values of the CKM matrix elements | V_cb | and \Vert V_ub | are obtained from published and preliminary measurements available in Summer 99 from the ALEPH, CDF, DELPHI, L3, OPAL and SLD Collaborations.

Accurate determinations of b-hadron decay properties provide constraints on the values of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [1]. The |V cb | and |V ub | elements can be obtained from semileptonic decay rates into charmed and non-charmed hadrons, and measurements of the oscillation frequency in B 0 d − B 0 d and B 0 s − B 0 s systems give access to |V td | and |V ts |.
Elements of the CKM matrix govern weak transitions between quarks. Experimental results are obtained from processes involving b-hadrons. Effects from strong interactions have thus to be controlled; and b-hadrons are also a good laboratory in this respect. Lifetime differences between the different weakly decaying hadrons can be related to interactions between the heavy quark and the light quark system inside the hadron. The polarization of Λ 0 b baryons produced by b-quarks of known polarization, emitted from Z decays at LEP or SLC, indicates how polarization is transmitted from the heavy quark to the baryon(s) in the hadronization process. Rates and decay properties of excited D * * states 1 produced in b-hadron semileptonic decays are needed to obtain accurate determinations of |V cb |. Decays of D * * states are also an important source of background in other channels, and their properties have to be monitored. Finally, it is mandatory to measure the production rates of the different weakly decaying b-hadrons produced during the hadronization of b-quarks created in high energy collisions, because all these states have different properties so the study of any one of them requires control of the background from the others.
Results obtained on b-hadron lifetimes, b-hadron production rates, B 0 d −B 0 d and B 0 s −B 0 s oscillations, and the values of the CKM matrix elements |V cb | and |V ub | made available during Summer 1999 are presented here. A limit on the decay width difference of the mass eigenstates of the B 0 s − B 0 s system is also given. These quantities have been obtained by averaging published and preliminary measurements released publicly by the ALEPH, CDF, DELPHI, L3, OPAL and SLD experiments. Whenever possible, the input parameters used in the various analyses have been adjusted to common values, and all known correlations have been taken into account. Close contacts have been established between representatives from the experiments and members of the different working groups in charge of the averages, to ensure that the data are prepared in a form suitable for combinations. Working group activities are coordinated by a steering group 2 .
Section 2 presents the values of the common input parameters that contribute to the systematic uncertainties given in this note. Published results obtained by other experiments, mainly operating at the Υ(4S) resonance, have been included to obtain most accurate values. In particular, it contains studies on some characteristics of D * * mesons in semileptonic b-decays, and on the Λ 0 b polarization. As some measured quantities are needed for the evaluation of others, an iterative procedure has been adopted to obtain stable results.
Section 3 describes the averaging of the b-hadron lifetime measurements. The combined values are compared with expectations from theory. The inclusive b-hadron and B 0 d meson 3 lifetimes are needed in the determination of |V ub | and |V cb | in order to convert measured branching fractions into partial widths that can then be compared with theory.
In Section 4, oscillations of neutral B mesons are studied. Also the production rates of the different b-hadrons in jets induced by a b-quark are determined using direct measurements together with the constraints provided by B mixing. The B 0 d production rate is an important input for the |V cb | measurement using B 0 d → D * + ℓ − ν ℓ decays, and the sensitivity to B 0 s -B 0 s oscillations depends on the B 0 s production rate. In Section 5, a limit on the decay width difference between mass eigenstates of the B 0 s -B 0 s system is given. As this average is obtained for the first time, more details on the adopted procedure have been given than for lifetimes or oscillations measurements.
The determination of |V cb | and |V ub | presented in Sections 6 and 7 includes only LEP results. These averages are also released for the first time. The determination of |V ub | is based on a novel technique using the lepton momentum and the mass of the hadronic system. The accuracy of these results (especially on |V cb |) depends mostly on theoretical uncertainties. Appendix E gives details on the theoretical inputs used.
A summary of all results obtained by the different working groups is given in Section 8. In addition, Appendices A-D contain, respectively, the individual measurements of the production rates of narrow D * * states, Λ 0 b polarization, b-hadron lifetimes, and direct measurements of b-meson and b-baryon production rates, that have been used in the present averages.

Common input parameters
The b-hadron properties used as common input parameters in these averages are given in Table 1. Most of the quantities have been taken from results obtained by the LEPEWWG [2,3] or quoted by the PDG [4]. The others, which concern the production rates and decay properties of D * * mesons in b-hadron semileptonic decays and the value of the Λ 0 b polarization, are explained later in this section.
Note the following: These values apply only to b-hadrons produced in Z decays. R b and R c are the respective branching fractions of the Z boson into bb and cc pairs in hadronic events, < x E > is the mean fraction of the beam energy taken by a weakly decaying bhadron.
• shape of the b-quark fragmentation function.
The value of < x E > obtained assuming the Peterson function [7], is given in Table  1. To evaluate the corresponding systematic uncertainty, parameter(s) governing b-quark fragmentation functions, taken from other models, have been varied in accordance with the uncertainty quoted for < x E >. Fragmentation functions taken from two models [8,9] have been chosen to estimate the systematic uncertainties coming from the shape of the function. These models typically yield results on either side of those obtained using the Peterson function which is commonly used by the experiments. For analyses which are rather insensitive to the fragmentation  Table 1: Common set of input parameters used for the derivation of the various measurements presented in this paper. The first set of results has been taken from those obtained by the LEPEWWG [2,3] or quoted by the PDG [4], the second set corresponds to averages obtained in the present report and, for the sake of completeness, in the last set, branching fractions in charm decays needed for the production rates of these particles are listed as well. The abbreviated notations sl., r.s. and. w.s. correspond, respectively, to semileptonic, right sign and wrong sign.
uncertainty, it is considered adequate to use only the Peterson model and to inflate the uncertainty on < x E > to ±0.02.
• the inclusive semileptonic branching fraction of b-hadrons.
The largest contribution to the systematic error comes from the uncertainty on the semileptonic decay model. This error is reduced from ±0.084×10 −2 to ±0.065×10 −2 when the fit is performed including the forward-backward asymmetry measurements in addition to all heavy flavour measurements performed at the Z. This happens because only asymmetry measurements obtained using leptons depend on the semilep-tonic decay model, measurements using a lifetime tag combined with jet-charge or Dmeson reconstruction do not. To achieve consistency between these measurements, the fit effectively constrains the size of the error attributed to the semileptonic decay model, thus reducing the corresponding systematic error on the semileptonic branching fraction. Including asymmetry measurements gives an inclusive b-hadron semileptonic branching fraction of (10.62 ± 0.17)%. To be conservative, we have not included the asymmetry measurements in the fit to extract the semileptonic branching fraction used in the present averages. The average value was obtained with the average B 0 − B 0 mixing parameter and R b value given in Table 1.
In the absence of direct measurements of the semileptonic branching fractions for the different B meson states, it has been assumed, when needed in the following analyses, that all b-hadron semileptonic widths are equal. This hypothesis is strictly valid for B − and B 0 d → cℓ − ν ℓ decays because of isospin invariance originating from the b → c transition, which is ∆I = 0 (similar considerations apply also to D 0,+ → sℓ + ν ℓ decays). It is not valid in B − and B 0 d → uℓ − ν ℓ decays, but the induced difference between B − and B 0 d total semileptonic decay rates can be neglected. It has been assumed valid also for B 0 s mesons and b-baryons, but for b-baryons an uncertainty of 15% has been added, estimated by comparing the lifetime ratios and semileptonic branching fraction ratios for b-mesons and b-baryons [10]. Exclusive semileptonic branching fraction averages, given in the following for the B 0 d meson, have been obtained using this hypothesis. Results for other b-hadron flavours can be obtained using, in addition, the corresponding lifetime ratios. The latter are obtained from b-hadron lifetimes given in Section 3.
• gluon splitting to heavy quarks.
The quantities P(g → cc) and P(g → bb) are defined as the ratios BR(Z→qqg, g→QQ)

BR(Z→hadrons)
in which, respectively, Q is a c or a b quark.
• b-hadron decay multiplicity. Table 1 is an average of DELPHI [11] and OPAL [12] measurements which does not include charged decay products from the long lived particles K 0 s and Λ 0 .

Inclusive D * * production rate in b-hadron semileptonic decays
The inclusive b-hadron semileptonic branching fraction into D * * mesons has been measured in three different ways: -by subtracting the contributions of B 0 d → (D + + D * + )ℓ − ν ℓ from the total semileptonic branching fraction of B 0 d mesons, yielding: using the inclusive semileptonic branching fraction given in Table 1 and published values for the exclusive rates [4].
-from a measurement at the Υ(4S) of the rate of final states with a D * + and using models to account for the other channels [13]: -from an inclusive measurement of final states in which a D or a D * + meson is accompanied by a charged hadron and assuming that non-strange D * * decay channels involve only Dπ and D * π final states: As the χ 2 of these four measurements is equal to three per degree of freedom, the uncertainty on the weighted average has been multiplied by √ 3, giving: 2.2 D * * decays to D * mesons in semileptonic b-decays Semileptonic decays to excited charm states which subsequently decay to a D * + are a source of correlated (physics) background in studies of B 0 d meson properties. It is appropriate to express the different measurements in terms of the parameter b * * , defined as the branching fraction of B 0 d semileptonic decays involving D * * final states in which a D * , charged or neutral, is produced: This is to differentiate such D * * decays from those into a D meson directly. Throughout this section it is assumed that decays of D * * mesons involve at most one pion (or one kaon for strange states).
The following measurements have been interpreted in terms of the quantity b * * and the production fractions (f B i ) and lifetimes (τ (B i )) of the different types of weakly decaying b-hadrons: -semi-inclusive measurements of semileptonic decays in which a D * + and a charged pion have been isolated: -the ARGUS measurement [13]: where the value of 0.77 corresponds to the modelling used for D * * decays in that analysis; -the inclusive production rate of charged D * mesons in semileptonic b decays: +f which has been obtained by averaging results quoted by the PDG [4] for corresponding decays of B 0 d and B − mesons, and the scaling factor α = 0.75 ± 0.25 has been introduced in (8) to account for a possible SU(3) flavour violation when comparing D * * + s → D * + K 0 and D * * 0 → D * + π − decays.
In the above, the production fractions (f B i ) and lifetimes (τ (B i )) are taken from Sections 4 and 3 respectively.
These measurements form a coherent set of results yielding: From this result, and using the same hypotheses, it is possible to derive other quantities of interest in several analyses presented in this paper such as: It is also possible to determine the fraction of B 0 d semileptonic decays which contain a D * + : (14) and the fraction of b-quark jets in which the D * + comes from a D * * decay: As these results (10) to (15) are highly correlated they are represented in Table 1 by a single entry, chosen to be result (11).

Other semileptonic decays to D * + mesons
Other mechanisms giving a D * + accompanied by a lepton are important because they are further sources of background for the exclusive channel B 0 d → D * + ℓ − ν ℓ .
2.3.1 Charged D * production in semileptonic decays of b-hadrons involving τ leptons Table 1 is obtained from the average measured value of BR(b → τ ν τ X) taken from the PDG [4] and assuming that a D * + is produced in (50 ± 10)% of the cases, as measured in semileptonic decays involving light leptons. The uncertainty on this last number has been increased, as compared to the value obtained in Equation (14), to account for a possible different behaviour of decays when a τ is produced instead of a light lepton because of differences in the masses involved in the final state.

Charged
The double charm production rate in b-decays was measured by the ALEPH [17] and CLEO [18] Collaborations, isolating the contributions of different D meson species. From these measurements the inclusive semileptonic branching fraction, involving a wrong sign charmed meson, has been obtained [3]: The product BR(b → D * + X c Y) × BR(X c → ℓ − X) quoted in Table 1 is then determined assuming again that a D * + is produced in (50 ± 10)% of the cases.

Production of narrow D * * states in b-hadron semileptonic decays
A modelling of b-hadron semileptonic decays requires detailed rate measurements of the different produced states which can be resonant or non-resonant D ( * ) nπ systems 4 , each having characteristic decay properties. As an example, there are four orbitally excited states with L=1. They can be grouped in two pairs according to the value of the spin of the light system, j = L ± 1/2 (L=1).
States with j = 3/2 can have J P = 1 + and 2 + . The 1 + state decays only through D * π, and the 2 + through Dπ or D * π. Parity and angular momentum conservation imply that in the 2 + the D * and π are in a D wave but allow both S and D waves in the 1 + . However, if the heavy quark spin is assumed to decouple, conservation of j( = 3/2) forbids S waves even in the 1 + . An important D wave component, and the fact that the masses of these states are not far from threshold, imply that the j = 3/2 states are narrow. These states have been observed with a typical width of 20 MeV/c 2 , in accordance with the expectation.
On the contrary, j = 1/2 states can have J P = 0 + and 1 + , so they are expected to decay mainly through an S wave and to be broad resonances with typical widths of several hundred MeV/c 2 . The first experimental evidence for broad 1 + states has been obtained by CLEO [19].
At present, information on the composition of b-hadron semileptonic decays involving D * * is rather scarce. ALEPH [14] and CLEO [24] have reported evidence for the production of narrow resonant states (D 1 and D * 2 ), which can be summarised as follows: = 0.37 ± 0.16 or < 0.6 at the 95% CL (19) where R * * is the ratio between the production rates of D * 2 and D 1 in b-meson semileptonic decays. The original measurements are listed in Appendix A. Absolute rates for the D 1 and D * 2 mesons have been obtained assuming that the former decays always into D * π and the latter decays into Dπ in (71 ± 7)% of the cases, as measured by ARGUS [25] and CLEO [26].
It is of interest to note that these results are very different from naive expectations from HQET in which R * * ≃ 1.6 [20]. This implies large 1/m c corrections in these models, as was noted for instance in [21]. These corrections not only reduce the expected value for R * * , they also enhance the expected production rate of j = 1/2 states which can become of similar importance as the j = 3/2 rate.
Comparing Equations (17)(18) and (4), it can be deduced that narrow D * * states account for less than one third of the total D * * rate. If, in addition, model expectations for the production rates of j = 1/2 states are used, it can be concluded that non resonant production of D ( * ) π states is as important as, or even dominates over, the production of resonant states. As a consequence, models like ISGW [22,23], which contain only resonant charmed states, do not provide a complete description of b-hadron semileptonic decays.

Λ 0 b polarization
Even for unpolarized e ± beams, the b-quarks produced in Z decays are highly polarized: P b = −0.94. Hard gluon emission and quark mass effects are expected to change P b by only 3% [27]. However, during the hadronization process of the heavy quark, part or all of the initial b-quark polarization may be lost by the final weakly decaying b-hadron state. The b-mesons always decay finally to spin zero pseudoscalar states, which do not retain any polarization information. In contrast, Λ 0 b baryons are expected to carry the initial b-quark polarization since the light quarks are arranged, according to the constituent quark model, in a spin-0 and isospin-0 singlet. However, b-quark fragmentation into intermediate Σ The Λℓ final state is used to select event samples enriched with Λ 0 b semileptonic decays, and the polarization is measured using the average values of the lepton and neutrino energies. Samples of events enriched in b-mesons, which carry no polarization information, are used to calibrate the measurements.
Measurements from ALEPH [29], DELPHI [30] and OPAL [31], collected in Appendix B, have been averaged, assuming systematic uncertainties to be correlated, yielding: where the uncertainty is still dominated by the statistics. The following sources of systematic uncertainties, common to all measurements, have been considered: -Λ production in semileptonic and inclusive decays of the Λ + c , -Λ + c polarization, -theoretical uncertainties in the modelling of polarized Λ 0 b semileptonic decays.

Averages of b-hadron lifetimes
Best estimates for the various b-hadron lifetimes, for the ratio of the B + and B 0 d lifetimes, and for the average lifetime of a sample of b-hadrons produced in b-jets have been obtained by the B lifetime working group 5 . Details on the procedure used to combine the different measurements can be found in [32].
The measurements used in the averages presented here are listed in Appendix C. Possible biases can originate from the averaging procedure. They depend on the statistics and time resolution of each measurement and in the way systematic uncertainties have been included. These effects have been studied in detail using simulations and have been measured to be at a level which can be neglected as compared to the statistical uncertainty.

Dominant sources of systematic uncertainties
Dominant sources of systematic uncertainties leading to correlations between measurements are briefly reviewed (more details can be found in [32]). These are the background estimation, the evaluation of the b-hadron momentum, and the decay length reconstruction.
The background can be due either to physics processes leading to a final state similar to that used to tag the signal, or to accidental combinations of tracks which simulate the decay of interest. When "physics" background is present, the experimental uncertainties on the branching fractions of the background processes and on the lifetimes of the background particles lead to a systematic error which is correlated between different experiments. When the background is combinatorial, the amount and/or the lifetime of the background "particles" is normally extracted from the data using the sidebands of mass distributions or wrong sign combinations. In this case the related systematic uncertainty is usually not correlated between experiments. But, in the measurement of the b-baryon  lifetime using Λℓ correlations, the amount of accidental background is obtained from the wrong sign combinations and a correction, common to all analyses, has to be applied to this number to take into account the production asymmetry of accidental Λℓ pairs.
Most of the exclusive B lifetime measurements are based on the reconstruction of B decay length and momentum. In most analyses the B particles are only partially reconstructed, and their energies are estimated from the energies of the detected decay products. In all cases, systematic uncertainties have been evaluated for the following effects: • determination of the b-quark fragmentation function, • branching fractions of B and charmed hadrons, • b-hadron masses, • b-baryon polarization and • modelling of neutral hadronic energy.
Finally there are uncertainties correlated within an experiment. They are due to primary and secondary vertex reconstruction procedures, detector resolution, tracking errors, B flight direction reconstruction and detector alignment.

Measurements of B 0 d and B + lifetimes
Apart from the measurements by CDF using large samples of exclusive B 0 d → J/ψK ( * )0 and B + → J/ψK + decays, the most precise measurements of B 0 d and B + lifetimes originate from two classes of partially reconstructed decays. In the first class the decay B → D ( * ) ℓ + ν ℓ X is used in which the charge of the charmed meson distinguishes between neutral and charged B mesons. In the second class the charge attached to the B decay vertex is used to achieve this separation.
The following sources of correlated systematic uncertainties have been considered: background composition (includes D * * branching fraction uncertainties as obtained in Section 2), momentum estimation, lifetimes of B 0 s and b-baryons (as obtained in Sections 3.3 and 3.4), and fractions of B 0 s and b-baryons produced in Z decays (as measured in Section 4.1).
The world average lifetimes of B 0 d and B + mesons are: The several measurements and the average values are shown in Figure 1.

B 0 s lifetime measurements
The most precise measurements of the B 0 s lifetime originate from partially reconstructed decays in which a D − s meson has been completely reconstructed (see Figure 2-left). The following sources of correlated systematic uncertainties have been considered: average b-hadron lifetime used for backgrounds, B 0 s decay multiplicity, and branching fractions of B and charmed hadrons.
The world average lifetime of B 0 s mesons is equal to: This result has been obtained neglecting a possible difference, ∆Γ B 0 s , between the decay widths of the two mass eigenstates of the B 0 s − B 0 s system. The measurement of ∆Γ B 0 s is explained in Section 5. The several measurements and the average value are shown in Figure 2-left.

b-baryon lifetime measurements
The most precise measurements of the b-baryon lifetime originate from two classes of partially reconstructed decays (see Figure 2-right). In the first class, decays with an exclusively reconstructed Λ + c baryon and a lepton of opposite charge are used. In the second class, more inclusive final states with a baryon (p, p, Λ or Λ) and a lepton have been used.
The following sources of correlated systematic uncertainties have been considered: experimental time resolution within a given experiment, b-quark fragmentation distribution into weakly decaying b-baryons, Λ 0 b polarization, decay model, and evaluation of the b-baryon purity in the selected event samples. As the measured b-hadron lifetime is proportional to the assumed b-hadron mass, the central values of the masses are scaled to m(Λ b ) = (5624 ± 9) MeV/c 2 and m(b-baryon) = (5670 ± 100) MeV/c 2 , before computing the averages. Uncertainties related to the decay model include mostly assumptions on the fraction of n-body decays. To be conservative it is assumed that they are correlated whenever given as an error. Furthermore, in computing the average, results have been corrected for the effect of the measured value of the Λ 0 b polarization and it has been assumed that b-baryons have the same fragmentation distribution as all b-hadrons (Section 2).
The world average lifetime of b-baryons is then: This value and the single measurements are given in Figure 2-right.
Keeping only Λ ± c ℓ ∓ and Λℓ − ℓ + final states, as representative of the Λ 0 b baryon, the following lifetime is obtained: Averaging the measurements based on the Ξ ∓ ℓ ∓ final states gives a lifetime value for a sample of events containing Ξ 0 b and Ξ − b baryons:

Average b-hadron lifetime
Available results have been divided into three different sets and a separate average has been computed for each set (see Figure 3-left): (1) measurements at LEP and SLD which accept any b-hadron decay: (2) measurements at LEP based on the identification of a lepton from a b-hadron decay: (3) measurement at CDF based on the identification of a J/ψ from a b-hadron decay: The reason for this division is that, since the lifetimes of the individual b-hadron species are different, a meaningful average lifetime can only be computed for samples in which the composition is the same.
The following sources of correlated systematic uncertainties have been considered when evaluating the averages: b-and c-quark fragmentation and decay models, BR(b → ℓ), BR(b → c → ℓ), τ c , BR(c → ℓ) and the charged track multiplicity in b-hadron decays (see Table 1).
The measurement in Set 1 is the average b-hadron lifetime for a sample of weakly decaying b-hadrons produced in Z decays.
The measurements of Set 2 are based on a sample enriched with B + with respect to other production fractions, and thus are expected to correspond to a higher average lifetime. In practice this difference is expected to be small. Assuming that all b-hadrons have the same semileptonic partial width, the average lifetimes in the inclusive and semileptonic samples can be related using the measured values of the production rates (Section 4) and lifetimes of the different types of b-hadrons, giving: Finally the CDF J/ψ measurement probably corresponds to a different composition in b-hadron species compared to any of the previous sets of measurements.
As a consequence, the average b-lifetime, τ b , is obtained by combining the values of τ incl. b and τ sl b only, and neglecting the small correction expected in Equation (31): The difference between the values of τ incl. b and τ sl b amounts to less than two standard deviations.

b-hadron lifetime ratios
Ratios of other b-hadron lifetimes to the B 0 d lifetime have also been determined, using the previous averages of the individual lifetimes and the direct measurements of the ratio between B + and B 0 d lifetimes shown in Figure 3-right. The results are given in Table 2.
Lifetime ratio Measured value Table 2: Ratios of b-hadron lifetimes relative to the B 0 d lifetime.
It can be seen that the predictions [71] for the meson lifetimes based on factorization, given in Equation (33), which had encountered theoretical criticism [72], are on the mark. A recent lattice study [73] finds a result quite consistent with [71].
However a large discrepancy is observed for baryons. It remains to be clarified if this problem is related only to the validity of the quark models used to determine the parameters in the case of baryons, as was suggested in reference [71], or to more basic defects in the whole theoretical picture. Lattice evaluations for b-baryons are on the way [74].

B 0 d and B 0 s oscillations and b-hadron production fractions
The four LEP Collaborations and CDF and SLD have published or otherwise released measurements of ∆m d [75,76,77,78,79,80] and lower limits for ∆m s [81,82,83,84,85]. Combined results, as well as estimates of b-hadron fractions, have been prepared by the B oscillation working group 6 .
The estimates of the b-fractions, described in Section 4.1, are important inputs needed for the combination of the B 0 d and B 0 s oscillation results. The procedure used to combine ∆m d results is explained in Section 4.2, followed by a description of the common systematic uncertainties in the analyses, and the result for the mean value of ∆m d . Then combined results on the B 0 s oscillation amplitude, as well as an overall lower limit on ∆m s are presented in Section 4.3. More details on these procedures can be found in reference [86].

Measurements of b-hadron production rates in b-jets
In this analysis, the relative production rates of the different types of weakly-decaying b-hadrons are assumed to be similar in b-jets originating from Z decays and in high transverse momentum b-jets in pp collisions at 1.8 TeV. This hypothesis can be justified by considering the last steps of jet hadronization to be a non-perturbative QCD process occurring at a scale of order Λ QCD .
Direct information on production rates is available from measurements of branching fraction products using channels with characteristic signatures. All the measurements used in the present analysis are listed in Appendix D.
At CDF and LEP, the B 0 s production rate has been evaluated using events with a D + s accompanied by a lepton of opposite sign in the final state (Table 21). Since the rate of these events is given by . This has been done by assuming, from SU(3) symmetry, that the partial semileptonic decay widths into D, D * and D * * final states are the same for all B mesons. A lower value for the B 0 s semileptonic branching fraction with a D − s in the final state is obtained by assuming that all D * * s states decay into a non-strange D meson. The branching fraction corresponding to B 0 . This corresponds to a lower limit on D − s production because the possibility that D * * − s states decay into D − s KK or D − s η ( ′ ) has been neglected. In addition, assuming that only the measured D ( * ) π final state in B 0 d or B + semileptonic decays corresponds to D ( * ) K transitions for the B 0 s , an upper value can be obtained for the B 0 s decay of interest using the relation: based on isospin symmetry. In practice, the lower and upper values are compatible within uncertainties. This allows f Bs = (11.2 +3.1 −2.7 )% to be extracted from the LEP measurements. Using the same assumptions, the CDF Collaboration has measured the ratio f Bs /(f B + + f B d ) = 0.213 ± 0.068 from final states with an electron and a charm meson In a similar way, the fraction of b-baryons is estimated from the measured production rates of Λ + c ℓ − [93,94] and Ξ − ℓ − [95,96] final states (Table 22) yielding, respectively, has been obtained considering that there could be one Λ + c produced in every decay or, at the other extreme, that all excited charmed baryons, assumed to be produced with a rate similar to D * * final states in B mesons, decay into a D meson and a non-charmed baryon. Similar considerations have been applied to Ξ b semileptonic decays. The semileptonic decay width Γ(Ξ b → Ξ c Xℓ − ν ℓ ) has been taken to be equal to Γ(Λ 0 b → Λ + c Xℓ − ν ℓ ) and it has been assumed that, within present uncertainties, BR(Ξ 0 c → Ξ − X) = BR(Λ + c → Λ 0 X). The total b-baryon production rate is then: f b-baryon = (12.0 +3.9 −2.7 )%, assuming the same production rates for Ξ 0 b and Ξ − b baryons. This is then averaged with a direct measurement of f b-baryon = (10.2 ± 2.8)% from the number of protons in b-events [97]. Finally, the CDF Collaboration has measured the production rate of b-baryons relative to non-strange B mesons, [90]. The fraction of B + mesons has also been measured [98] using the charge of a large sample of inclusively reconstructed secondary vertices: f B + = (41.4 ± 1.6)% (Table 23).
All these measurements have been combined to obtain average production rates for B 0 s , b-baryons, B + and B 0 d 8 , imposing that: and The results obtained are given in Table 3. Correlated systematics between the different measurements, coming mainly from the poorly measured branching fractions of D + s and Λ + c charmed hadrons, have been taken into account. It may be noted that constraint (36) is expected to be true for b-mesons even though it is known to be wrong for c-mesons. In bjets, and also in c-jets, isospin invariance of strong interactions predicts similar production rates of mesons in which the heavy quark is associated to a u or d antiquark. Strong and electromagnetic decays of these states result in different rates for weakly decaying mesons with a u or a d flavour because the D * → Dπ decays occur very close to threshold, and the threshold prevents the transition D * + → D + π − . However for B 0 d and B + mesons, no asymmetry is expected: B * mesons decay electromagnetically, leaving the flavour of the light spectator quark in the b-meson unchanged; non-strange B * * mesons decaying through strong interactions and having masses away from threshold will not induce any asymmetry either; B * * 0 s mesons decay also by strong interactions into B 0( * ) d K 0 or B +( * ) K − with equal probabilities, although a possible tiny difference between these two rates can be expected for decays of narrow states occurring very close to threshold, because of the mass difference between K 0 and K + mesons 9 .
b-hadron fractions Table 3: Average values of b-hadron production rates and their correlation coefficients obtained from direct measurements.
b-hadron fractions correlation coefficients Table 4: Average values of b-hadron production rates and their correlation coefficients obtained from direct measurements and including constraints from results on B mixing.
Additional information on the production rates can be obtained from measurements of the time-integrated mixing probability of b-hadrons. For an unbiased sample of semileptonic b-hadron decays in b-jets, with fractions g B d and g Bs of B 0 d and B 0 s mesons, this mixing probability is equal to: where χ d is the time-integrated mixing probability for B 0 d mesons (see Section 4.2) and χ s =1/2 is the corresponding quantity for B 0 s mesons 10 . As already mentioned in Section 2, the semileptonic width is assumed to be the same for all b-hadron species, implying are the lifetime ratios. This leads to the relation: . This is used to extract another determination of f Bs from the b-baryon fraction of Table 3, the lifetime ratio averages of Table 2, the world average value of χ d from Equation (44) of Section 4.2, and the χ average from the LEPEWWG (see Table 1). This new estimate of f Bs , f mixing Bs = (10.3 ± 1.5)%, is then combined with the b-hadron rates from direct measurements, taking into account correlations 11 and imposing the conditions (35) and (36). The final b-hadron fractions are displayed in Table 4.
The production rate for bs states is expected to be rather similar to the probability (γ s ∼ 12%) for creating an ss pair during the hadronization process. But, because of their masses, non strange B * * mesons are not expected to decay into B 0( * ) s K, whereas strange B * * mesons are expected to decay mainly into non-strange weakly decaying b-mesons. As 10 The assumption χ s = 1 2 can be justified by the existence of limits on ∆m s obtained from D s -lepton analyses, which have negligible dependence on f Bs . 11 There is a small statistical correlation between χ d and χ, arising from the fact that a few ∆m d analyses at LEP are based on the same samples of dilepton events as the ones used to extract χ. This correlation is ignored, with a negligible effect on the final result.   a result, the fraction of B 0 s in b-jets is expected to be given by f Bs ≃ γ s (1 −P(b → B * * )) ≃ 10%, as observed.

Combination method for ∆m d
The measurements of ∆m d are now quite precise and it is important that the correlated systematic uncertainties are correctly handled. Furthermore, many results depend on physics parameters for which different values were used in the original analyses. Before being combined, the measurements of ∆m d are therefore adjusted on the basis of a common consistent set of input values.
For each input parameter the ∆m d measurement and the corresponding systematic uncertainty are linearly rescaled in accordance with the difference between the originally used parameter value and error, and the new common values. This is done for the systematic uncertainties originating from the b-hadron lifetimes, the b-hadron fractions and the mixing parameters χ and χ d . The statistical and systematic uncertainties of each individual measurement are symmetrized. An alternative approach using asymmetric uncertainties (when quoted as such) would produce negligible differences in the combined result.
The combination procedure makes a common fit of ∆m d and the common input parameters. It is assumed that ∆m d may be expressed as a function X(Y 1 . . . Y Nsys ) with a weak dependence on the systematic sources Y α . Expanding X then gives: where the quantities σ syst α are the correlated errors on ∆m d from systematic sources α, and Y α = (z α − z f it α )/δ α . In this last expression, z α and z f it α are the input and fitted values of the systematic parameter α, and δ α is the variation used to calculate σ syst α . The following χ 2 is then constructed: where X i are the measurements of ∆m d and σ uncor i the quadratic sum of the statistical and uncorrelated systematic uncertainties on X i . This χ 2 is minimized with respect to the parameters X 0 and Y α ; the result for ∆m d is taken as the value of X 0 at the minimum (and the values of Y α at the minimum are ignored). This method gives the same results as a χ 2 minimization with inversion of a global correlation matrix. As several measurements also have a statistical correlation, the first sum in Equation (40) is then generalized to handle an N × N error matrix describing the statistical and uncorrelated systematic uncertainties.
The following sources of systematic uncertainties, common to analyses from different experiments, have been considered: • b-lifetime measurements (Section 3). Different measurements use the b-hadron lifetimes in different ways as input: some use the actual lifetimes, and others use ratios of lifetimes. As this leads to complicated correlations between the analyses, the following procedure is adopted. For each measurement a single combined value of all b-lifetime-related systematic uncertainties is computed. This systematic is then treated as fully correlated with all other such numbers from the other measurements. Tests show that this procedure gives a negligible bias and results in a conservative evaluation of the uncertainty.
• b-quark fragmentation (Section 2). In some analyses the b-hadron momentum dependence is treated through the variation of the ǫ parameter in the Peterson fragmentation function, and in others through the mean scaled b-hadron energy < x E >.
In such a case, the corresponding systematic uncertainties are treated as fully correlated, without attempting to adjust the individual results to a common value of ǫ or < x E >.
• fraction of B + mesons remaining in B 0 d enriched samples for analyses using D * mesons. This depends on decay branching fractions and on the mistag probability.
Common systematics due to purely experiment-dependent factors (i.e. common to different results in a particular experiment) have not been included in the above list, but nonetheless are treated as correlated in the fit.
Following this procedure, a combined value: is obtained. Using the relation 12 : and the B 0 d lifetime of Section 3, the above ∆m d result can be converted to: Averaging with the time-integrated mixing results obtained by ARGUS and CLEO at the Υ(4S) (see Table 1), yields finally: The individual measurements of ∆m d and their combined value are shown in Figure 4. The determination of ∆m d as described above and of the b-hadron fractions described in Section 4.1 cannot be performed sequentially: the values of the fractions are needed to perform the ∆m d fit, and the best estimates of the fractions can be obtained only once the final ∆m d average is known. This circular dependence has been handled by including the calculation of the fractions in the ∆m d fitting procedure, in such a way that the final results quoted for ∆m d and the b-hadron fractions form a consistent set. 3 ps −1 is derived from this spectrum. Points with error bars are the fitted amplitude values and corresponding uncertainties, including systematics. The dark shaded area is obtained by multiplying these uncertainties by 1.645 such that the integral of the probability distribution, assumed to be Gaussian, is equal to 5% above this range. The clear shaded area is obtained in a similar way but including only statistical uncertainties in the fit; as a consequence, central values of the two regions do not necessarily coincide.
No experiment has yet directly observed B 0 s oscillations, so the task of a combination procedure is to calculate an overall limit or to quantify the evidence for a signal from the information provided by each measurement. This is done using the amplitude method of reference [99]. At each value of ∆m s , in the range of interest, an amplitude is measured in each analysis, where the expected value of the amplitude is unity at the true frequency. An overall limit on ∆m s is then inferred from the combined amplitude spectrum by excluding regions of ∆m s where the amplitude is incompatible with unity.
Studies with Monte Carlo simulation show that the measured amplitude has a Gaussian distribution around its expected value, which is zero for frequencies much lower than the true one, and unity at the true frequency. In addition, the error on the amplitude depends on the statistical power of the analysis but not on the true value of the oscillation frequency. Therefore, if at a given frequency (∆m s ) a measured amplitude (µ) with error σ is found, this value of the frequency can be excluded at the 95% C.L. if the probability of measuring an amplitude equal or smaller than µ, for a true amplitude of unity, is smaller than 5%, that is µ −∞ G(A − 1, σ)dA < 0.05, or equivalently µ + 1.645σ < 1 where G is the Gaussian function.
Before combining the measured amplitudes, the individual central values and systematic uncertainties are adjusted, using the same procedure as for ∆m d measurements, to common values of the b-hadron fractions (Table 4), b-hadron lifetimes (Section 3), and ∆m d (Equation (45)). The adjustment to a common value of f Bs is performed first, and needs a special treatment for certain analyses, as explained below. Although the final value of f Bs is calculated under the assumption that B 0 s mixing is maximal, studies indicate that the effect on the combined ∆m s result is negligible.
The statistical uncertainty on the measured amplitude is expected to be inversely proportional to the B 0 s purity of the analysed sample [99]. This purity is of the order of f Bs for inclusive analyses. For analyses where a full or partial B 0 s reconstruction is performed, the purity is much less dependent on the assumed value of f Bs . Therefore, the statistical uncertainties on the amplitudes measured in inclusive analyses have been multiplied by f used Bs /f new Bs , where f used Bs is the B 0 s fraction assumed in the corresponding analysis, and f new Bs is the B 0 s fraction from Table 4. This correction is also applied on the central values of the amplitude measured in inclusive analyses, after having checked that the relative uncertainty on the amplitude is essentially independent of f Bs for large enough values of ∆m s .
The amplitudes (measured at each value of ∆m s ) are averaged using the same procedure as for the combination of ∆m d results. In addition to systematics which are correlated within the individual experiments, the following sources of correlated systematic uncertainties have been taken into account: • b-hadron lifetime measurements, • b-quark fragmentation, • direct and cascade semileptonic branching fractions of b-hadrons, • b-hadron fractions, • ∆m d measurements, • ∆Γ B 0 s measurements.  Figure 6: Measurements of the B 0 s oscillation amplitude. The amplitudes are given at ∆m s = 10 ps −1 , along with the relevant statistical and systematic errors. The sensitivities are given within parentheses in the column on the right. Note that the individual measurements are as quoted in the original publications, but the average, corresponding to the hatched area, includes the effects of adjustments to a common set of input parameters. The dashed lines correspond to amplitude values of 0 and 1.
Data using the amplitude method have been provided by ALEPH, DELPHI, OPAL (for their lepton-jet analysis only 13 ), CDF and SLD. The combined amplitude spectrum is shown in Figure 5. All values of ∆m s below 14.3 ps −1 are excluded at the 95% CL. The combined expected sensitivity to ∆m s , defined as the largest value of ∆m s that would have been excluded if the combined amplitudes were zero at all values of ∆m s , is 14.5 ps −1 . In Figure 6 the values of the amplitude at ∆m s = 10 ps −1 from the various analyses, are shown. In addition, the sensitivities are given.
Positive measured amplitudes are found for frequency values near and above the sensitivity limit. Some values deviate from zero by more than twice the total estimated error. The likelihood profile obtained from the amplitude spectrum, using the prescription of [99], shows a minimum at 16.3 ps −1 . The likelihood value at this minimum is 1.9 units below the asymptotic likelihood value for ∆m s → ∞. Because the measurements at dif-ferent frequencies are correlated, it is not possible to calculate analytically the probability that, in the explored frequency range, a fluctuation as or more unlikely occurs in a sample where the true frequency is far beyond the sensitivity. Therefore, an estimation procedure based on fast Monte Carlo experiments has been developed [100]: the above probability is found to be (10 ± 2)%, where the quoted uncertainty is due to the assumptions and approximations adopted in the fast Monte Carlo generation. to first approximation [103]), thus offering the possibility of measuring ∆m s via the lifetime difference rather than the oscillation frequency. This could be particularly important if the oscillation frequency is too fast to be measured with the present experimental proper time resolution. In addition, if ∆Γ B 0 s does turn out to be sizable, the observation of CP violation and the measurement of CKM phases from untagged B 0 s samples can be imagined [104].
The existing experimental constraints on the width difference and the combination of these constraints are reported here 14 .   has also been proposed recently [106]. The available results are summarised in Table 5. The values of the limit on ∆Γ B 0 s /Γ B 0 s quoted in the last column of this table have been obtained by the working group.

Experimental constraints on ∆Γ
Methods based on double exponential lifetime fits to samples containing a mixture of CP eigenstates have a quadratic sensitivity to ∆Γ B 0 s (inclusive, semileptonic, D + s -hadron), whereas methods based on isolating a single CP eigenstate have a linear dependence on ∆Γ B 0 s (φφ, J/ψφ). The latter are therefore, in principle, more sensitive to ∆Γ B 0 s ; but they tend to suffer from reduced statistics.
In order to obtain an improved limit on ∆Γ B 0 s , the results based on fits to the proper time distributions are used to apply a constraint on the allowed range of 1/Γ B 0 s . The world average B 0 s lifetime is not used, as its meaning is not clear if ∆Γ B 0 s is non-zero. Instead, it is chosen to constrain 1/Γ B 0 s to the world average τ B 0 d lifetime (1/Γ B 0 s ≡ 1/Γ B 0 d = τ B 0 d = (1.562 ± 0.029) ps). This is well motivated theoretically, as the total widths of the B 0 s and B 0 d mesons are expected to be equal within less than one percent [71,108] and ∆Γ B 0 d is expected to be negligible.
Further information on the various individual measurements is now given.
• L3 inclusive b-sample: in an unbiased inclusive B sample, all B 0 s decay modes are measured, including decays into CP eigenstates. An equal number of B short s and B long s mesons are therefore selected, and the proper time dependence of the B 0 s signal is given by: If the proper time dependence of this sample is fitted assuming only a single exponential lifetime then, using the definitions of Γ B 0 s and ∆Γ B 0 s and assuming that ∆Γ B 0 s /Γ B 0 s is small, the measured lifetime is given by: L3 effectively incorporates P incl (t) into the proper time fit of an inclusive b-sample and applies the constraint 1/Γ B 0 s = (1.49 ± 0.06) ps.
• DELPHI B 0 s → D + s ℓ − ν ℓ X: in a semileptonic selection, the ratio of short and long B 0 s 's selected is proportional to the ratio of the decay widths and the proper time dependence of the signal is: If this sample is fitted assuming only a single exponential lifetime for the B 0 s , then the measured lifetime is, always in the limit that ∆Γ B 0 s /Γ B 0 s is small, given by: The single lifetime fit is thus more sensitive to the effects of ∆Γ B 0 s in the semileptonic than in the fully inclusive case. Information on ∆Γ B 0 s is obtained by scanning the likelihood as a function the two parameters 1/Γ B 0 s and ∆Γ B 0 s /Γ B 0 s and applying the τ B 0 d constraint. • OTHERS: other analyses of the B 0 s semileptonic lifetime [105] have not explicitly considered the possibility of a non-zero ∆Γ B 0 s value. Nevertheless, the fact that the single exponential lifetime for this case (Equation (50)) is sensitive to ∆Γ B 0 s allows information on ∆Γ B 0 s to be extracted. The average B 0 s semileptonic lifetime has been recalculated excluding the DELPHI result just discussed, and information on ∆Γ B 0 s has been obtained by using Equation (50) and applying the τ B 0 d constraint. The validity of Equation (50) in the presence of background contributions has been verified using a toy Monte Carlo [109].
• ALEPH B 0 s → φφX (counting method): only those B 0 s decays which are CP eigenstates can contribute to a width difference between the CP even and CP odd states. An analysis [111] This lifetime is related to ∆Γ B 0 s and Γ B 0 s via the expression: ALEPH [106] has measured the lifetime of φφX events and extracted information on ∆Γ s /Γ s with the help of the world average B 0 s lifetime obtained from semileptonic B 0 s decays. The result listed in Table 5 has been obtained by the working group with the assumption 1/  (54) in which P semi (t) and P single (t) are defined in Equations (49) and (52) where P single has been defined in Equation (52). CDF measures the lifetime of J/ψφ events and information on ∆Γ B 0 s /Γ B 0 s is obtained after applying the 1/Γ B 0 s ≡ τ B 0 d constraint and including the experimental uncertainty on f short .

Combined limit on ∆Γ B 0 s
In order to combine the analyses summarised in Table 5, the result of each analysis has been converted to a two-dimensional log-likelihood in the (1/Γ B 0 s , ∆Γ B 0 s /Γ B 0 s ) plane. This log-likelihood has either been provided by each experiment or reconstructed from the measured lifetimes using the expected dependence of this quantity on 1/Γ B 0 s and ∆Γ B 0 s /Γ B 0 s . The latter procedure was necessary for the OTHERS and CDF entries of Table 5. The L3 analysis is not included in the average as the two-dimensional likelihood was not provided and could not be reconstructed from the available information. The τ B 0 d constraint is not applied on 1/Γ B 0 s at this stage. Systematic uncertainties are included in the individual log-likelihood distributions.
The log-likelihood distributions have been summed and the variation of the global negative log-likelihood function has been measured with respect to its minimum (∆L). The 68%, 95% and 99% C.L. contours of the combined negative log-likelihood are shown in Figure 7a. The corresponding limit on ∆Γ B 0 s /Γ B 0 s is: An improved limit on ∆Γ B 0 s /Γ B 0 s can be obtained by applying the τ B 0 d = (1.562 ± 0.029) ps constraint. When expressed as a probability density, this constraint is: The value of the uncertainty σ τ B 0 d includes, eventually, the theoretical uncertainty on the . The two-dimensional log-likelihood obtained, after including the constraint, supposed to be exact, is shown in Figure 7b. The resulting probability density distribution for ∆Γ B 0 s /Γ B 0 s is shown in Figure 7c. The corresponding limit on ∆Γ B 0 s /Γ B 0 s is: If an additional 2% uncertainty, assumed to be Gaussian distributed, is incorporated to account for the theory assumption τ B 0 d = 1/Γ B 0 s , the effect on the result is small: There are two methods to measure |V cb |: the inclusive method, which uses the semileptonic decay width of b-decays and the OPE; and the exclusive method, where |V cb | is extracted by studying the exclusive B 0 d → D * + ℓ − ν ℓ decay process using HQET. The B → Dℓ − ν ℓ channel has not been averaged to date.
In this note, both methods are used to determine values for |V cb |, which are then combined to produce a single average 16 . The semileptonic b-decay width, determined by the LEP heavy flavour electroweak fit to ALEPH, DELPHI, L3 and OPAL data, is used to determine BR(b → ℓ − ν ℓ X). Results from ALEPH [112] 17 , DELPHI [113] 18 and OPAL [114] are used to perform a LEP |V cb | average in the B 0 d → D * + ℓ − ν ℓ decay channel. These measurements are combined using a method similar to that used by the B oscillations working group.
Theoretical input parameters needed to extract |V cb | from actual measurements are detailed in Appendix E. In both analyses, the validity of the quark/hadron duality hypothesis has been assumed.

Inclusive |V cb | determination
In the inclusive method, the partial width for semileptonic B meson decays to charmed mesons is related to |V cb | using the following expression (Appendix E): where X c represents all final states containing a charmed quark. The rest of the expression, within parentheses, represents the correction to the muon decay formalism, depending on the b-and c-quark masses and on the strong coupling constant; µ 2 π is the average of the square of the b-quark momentum inside the b-hadron. The three last contributions to the uncertainty correspond respectively, to the uncertainties coming from the QCD perturbative expansion (scale dependence, truncation at finite order, ..), from the b-quark mass determination and from neglected terms of order m −n b,c , with n ≥ 3, appearing in the O.P.E. formalism.
Experimentally, the semileptonic width of b-hadrons is determined from its semileptonic branching fraction and lifetime. In Z decays, a mixture of B 0 d , B − , B 0 s and b-baryons is produced, such that the inclusive semileptonic branching fraction measured at LEP is an average over the different hadrons produced: where τ b is the average b-hadron lifetime. Therefore the semileptonic width of b-hadrons can be obtained using the inclusive semileptonic branching fraction and the average bhadron lifetime. The two last equalities in Equation (62) assume that all b-hadrons have the same semileptonic width. This hypothesis may be incorrect for b-baryons. Taking into account the present precision of LEP measurements of b-baryon semileptonic branching fractions and lifetimes, an estimate of the correction to Equation (62) is about 1.5% (see Section 2). The average LEP value for BR(b → ℓ − ν ℓ X) = (10.58 ± 0.07(stat.) ± 0.17(syst.))% is taken from the global fit which combines the heavy flavour measurements performed at the Z (see Section 2). The BR(b → ℓ − ν ℓ X u ) contribution is subtracted from BR(b → ℓ − ν ℓ X), using the LEP average value given in Section 7. For the average b-hadron lifetime, the world average value of τ b is used, as obtained in Equation (32).

Sources of systematic errors
The systematic errors assumed at present in the determination of the inclusive semileptonic decay width can be grouped into the following categories: • errors related to the efficiency and purity of the b-tagging algorithm.
As a lifetime b-tag is involved, effects due to the uncertainties in the sample composition in terms of different heavy hadrons, and uncertainties in the heavy hadron lifetimes, are considered.
• input parameters influencing the signal and background normalisation.
The values of the production fractions: , and the rates of gluon splitting P(g → cc) and P(g → bb) have been fixed to the values of [2]. R b , R c , BR(b → c → ℓ) and BR(c → ℓ) are parameters of the LEPEWWG fit. Values for all these quantities are listed in Table 1.
• the average fraction of the beam energy carried by the weakly decaying b-hadron.
Different models have been considered for the shape of the fragmentation function and the free parameters of the models have been determined from the data.
These effects on the lepton spectra have been included.
The average LEP value for BR(b → ℓ − ν ℓ X) is taken from the global LEPEWWG fit which combines the heavy flavour measurements performed at the Z [2], but removing the forward-backward asymmetry measurements (see Section 2).
• detector specific items.
The errors listed in the last item are uncorrelated among the different experiments. The others have been split into their uncorrelated and correlated parts. The error on the average b-hadron lifetime is assumed to be uncorrelated with the error on the semileptonic branching fraction. The propagation of these errors to the error on |V cb | is done assuming that they are Gaussian in the branching fraction and the lifetime, respectively.

Inclusive |V cb | average
Using the expression given in Equation (61), the following value is obtained: where the first error is experimental and the second is from theory. The experimental contributions due to the semileptonic branching fraction and the lifetime are ±0.37 ×10 −3 and ±0.18×10 −3 , respectively. The dominant systematic uncertainty, of theoretical origin, comes from the determination of the kinetic energy of the b-quark inside the b-hadron as explained in Appendix E.

Exclusive |V cb | determination
In the exclusive method, the value of |V cb | is extracted by studying the decay rate for the process B 0 d → D * + ℓ − ν ℓ as a function of the recoil kinematics of the D * + meson. The decay rate is parameterized as a function of the variable w, defined as the product of the four-velocities of the D * + and the B 0 d mesons. This variable is related to the square of the four-momentum transfer from the B 0 d to the ℓ − ν ℓ system, q 2 , by: and its values range from 1.0, when the D * + is produced at rest in the B 0 d rest frame, to about 1.50. Using HQET, the differential partial width for this decay is given by: where K(w) is a known phase space term and F D * (w) is the hadronic form factor for the decay. Although the shape of this form factor is not known, its magnitude at zero recoil, w = 1, can be estimated using HQET. In the heavy quark limit (m b → ∞), F (w) coincides with the Isgur-Wise function [117,118] which is normalised to unity at the point of zero recoil. Corrections to F D * (1) have been calculated to take into account the effects of finite quark masses and QCD corrections [119]. Calculations of this correction yield F D * (1) = 0.88 ± 0.05 (Appendix E). Since the phase space factor K(w) tends to zero as w → 1, the decay rate vanishes at w = 1 and the accuracy of the extrapolation relies on achieving a reasonably constant reconstruction efficiency in the region close to w = 1. The unknown function F D * (w) is approximated with an expansion around w = 1 due to Caprini, Lellouch and Neubert (CLN) [115]: where ρ 2 is the slope parameter at zero recoil and z = The ratio between the axial and vector form factors is included in K(w). Theoretical predictions restrict values of ρ 2 to be in the range: −0.14 < ρ 2 < 1.54. An alternative parametrization, obtained earlier, can be found in [116].

Sources of systematic uncertainties
The systematic uncertainties in the determination of |V cb | using the semileptonic decay B 0 d → D * + ℓ − ν ℓ can be grouped into the following categories: • normalisation: B 0 d meson production rate, D ( * ) branching fraction to the tagged final states (including topological BR), B 0 d lifetime (this is needed to obtain the B 0 d → D * + ℓ − ν ℓ decay partial width), and the b → B fragmentation function (which influences the reconstruction efficiency).
• background from physical processes: comprising B 0 d → D * + τ − ν τ , B 0 d → D * + X c (followed by the semileptonic decay τ /X c → ℓ − ν ℓ X) and, particularly, the intermediate production of excited charm mesons D * * which subsequently decay to a D * + ; • detector specific items: selection efficiency (lepton identification, tracking, vertexing), non physics background (combinatorial, hadron mis-identification), resolution, fitting, etc. This last set is treated as uncorrelated among experiments, and therefore will be ignored in the following discussion.

Normalisation
The B 0 d production rate at LEP is given by the product: where R b is taken from [2] (Table 1) and f B d from Section 4.1.
The values for the charm meson decay branching fractions are taken from the PDG [4]. Correlations among some of them (e.g. D 0 → K − π + with D 0 → Knπ, etc.) are included. Analyses based on the inclusive reconstruction of pions from D * + cascade decays may be affected by the knowledge of the topological D 0 branching fractions 19 ; they are taken from MARKIII measurements [120].
The B 0 d lifetime determined in Section 3.2 is used. Knowledge of the B 0 d fragmentation function is necessary in order to compute the fraction of B 0 d which were not reconstructed because they did not have enough energy to be detected. B 0 d hadrons produced in e + e − annihilations carry on average a large fraction, < x E >, of the beam energy (Table 1); consequently only a small fraction of them are outside the selection acceptance.

Physics background
Charged D * mesons can be accompanied by a lepton of opposite charge in two other b-hadron decay processes: 19 A topological branching fraction corresponds to the branching fraction into a fixed number of charged particles emitted in the final state. Their respective branching fractions have been evaluated in Section 2.3. The production rates of the different D * * states (see Sections 2.2 and 2.4) and the variation of their corresponding form factors as a function of w, have also to be considered. Published results on |V cb | are based on the old Isgur-Wise model [22], which predicts a sizeable D * * background rate near the end-point spectrum. As a consequence, in this model the error on the overall amount of D * * has a large effect on |V cb | while having a negligible contribution on the slope parameter ρ 2 . However, HQET predicts that, in the infinite charm mass limit, the rate near w = 1 is suppressed by a further factor (w 2 − 1) when compared with the signal [20,121]. In this case, the D * * rate uncertainty would have a large effect on the slope with only a small influence on |V cb | [113]. However, models in this extreme case fail to predict the ratio R * * (Equation (19)) between the production rates of the two narrow states.
A treatment which accounts for O(1/m c ) corrections is proposed in [21]. Several possible approximations of the form factors are provided, depending on five different expansion schemes and on three different input parameters. To be conservative, at the present stage of the analysis, each proposed scheme was tested in turn and the input parameters were varied over their full range. The quoted central value corresponds to the arithmetic average of the values obtained with the two extreme models. The systematic error due to the modelling of the D * * background was computed as half the difference between the two extreme results.

Corrections applied to the measurements
Since the three LEP measurements have been performed using different methods and inputs, they must be put on the same footing before being averaged. ALEPH and DELPHI measurements have been updated to use the CLN [115] extrapolation method (66) and the Ligeti [21] B 0 d → D * * + ℓ − ν ℓ model. Since OPAL uses the Caprini-Neubert extrapolation method and the JETSET B 0 d → D * * + ℓ − ν ℓ model, corrections to the OPAL measurement have been estimated using the ALEPH analysis, which is similar.
Corrections for changing to the standard input parameters, listed in Table 7, have been calculated as for the ∆m d measurement (see Section 4.2). The central value of each analysis is adjusted according to the difference between the used and desired parameter values and the associated systematic error. The systematic error itself is then scaled to reflect the desired uncertainty on the input parameter. Table 6 lists the corrected results for the three experiments.

Exclusive |V cb | average
The combination method for F D * (1)|V cb | and ρ 2 is the same as the method used for ∆m d (see Section 4.2) generalized to the combination of two or more correlated parameters. 34.2 ± 1.9 ± 2.6 1.10 ± 0.24 ± 0.21 Table 6: Experimental results on F D * (1)|V cb | and ρ 2 corrected for common inputs.
The LEP average (see Figure 9-right) gives: The parameters used are listed in Table 1. The uncertainty on the variation with w of the B 0 d → D * * + ℓ − ν ℓ form factors is taken to be fully correlated between experiments. The dominant systematic uncertainties on F D * (1)|V cb | are listed in Table 7. The largest comes from the B 0 d → D * * + ℓ − ν ℓ contribution. The confidence level of the fit is 15%. The error ellipses of the corrected measurements and of the LEP average are shown on Figure 8.

Overall |V cb | average
The combined |V cb | average (see Figure 9-right) can be extracted taking into account correlations between the inclusive and exclusive methods. The most important source of correlations comes from theoretical uncertainties in the evaluation of µ π , the average momentum of the b-quark inside the b-hadron (Appendix E). In the determination of experimental systematic uncertainties, theoretical uncertainties in the modelling of b → ℓ decays and the exact amount of b → D * * decays are taken as fully correlated. If these correlations were neglected, the central value for the |V cb | average would change by only 0.3%. Uncertainties from lepton identification and background also contribute, but to a much lesser extent. All other sources provide negligible contributions to the correlated error. The various contributions to the uncertainty are categorised in Table 8. The combined value is: where, within the total error of 1.9, 1.5 comes from uncorrelated sources and 1.2 from correlated sources.  7 Average of LEP |V ub | measurements The first LEP combined determinations of BR(b → ℓ − ν ℓ X u )(ℓ − = e − or µ − ) and the derivation of |V ub | have been obtained by combining the results reported by the ALEPH [122], DELPHI [123] and L3 [124] Collaborations 20 . The three analyses rely on different techniques to measure the inclusive yield of b → u transitions in semileptonic b-hadron decays. The experimental details can be found in the original publications.
All three experiments reported evidence for the b → ℓ − ν ℓ X u transition and measured its rate BR(b → ℓ − ν ℓ X u ). DELPHI fitted |V ub | / |V cb | directly to the fraction of candidate b → ℓ − ν ℓ X u decays in the selected data sample. For this averaging, the corresponding value of BR(b → ℓ − ν ℓ X u ) has been derived by using the value of |V cb | obtained in Section 6.3 and the relationship between |V ub | and BR(b → ℓ − ν ℓ X u ) given below. In order to average these results, the sources of systematic uncertainties have been divided into two categories. The first contains (a) uncorrelated systematics due to experimental systematics, such as lepton identification, b-tagging, vertexing efficiency and energy resolution, and (b) uncorrelated systematics from signal modelling and background description. The second contains correlated systematic uncertainties deriving from the simulation of signal b → u and background b → c transitions. The contributions from the statistical, experimental, and uncorrelated and correlated modelling uncertainties are summarised in Table 9.  The correlated systematics are summarised in Table 10. Differences in the analysis techniques adopted by the three experiments are reflected by the different sizes of the systematics uncertainties estimated from each common source. Important common systematics arise from the D topological branching fractions and the rate of D → K 0 decays. D decays represent a potential source of background for b → u decays because both are characterized by a small hadronic mass and a low charged multiplicity. The sensitivity to the topological branching fractions is reduced in the DELPHI analysis by applying a rescaling of the mass M X of the hadronic system, based on the reconstructed ℓ − ν ℓ X mass, and by the use of identified kaons for separating signal from background events. This explains the different sensitivity of the ALEPH and DELPHI analyses to these two important sources of systematic uncertainty. ALEPH and L3 are sensitive to the uncertainties in the b fragmentation function due to the kinematical variables used for discriminating b → ℓ − ν ℓ X u from b → ℓ − ν ℓ X c decays. The DELPHI result is sensitive to the assumed production rates of b hadron species due to the use of kaon anti-tagging to reject b → c, thus rejecting also B 0 s and Λ 0 b decays; and it is sensitive to the contribution of D ( * ) π and D * * states in semileptonic decays because the resulting difference in the vertex topology is also used for discriminating b → u from b → c decays.
The systematic uncertainties on the b → u signal have been grouped into inclusive model and exclusive model classes, which are assumed to be fully correlated. The first corresponds to the uncertainty in the modelling of the kinematics of the b-quark in the heavy hadron. It has been estimated from the spreads of the results obtained with the ACCMM [125], the Dikeman-Shifman-Uraltsev [126], and the parton [127] models in the ALEPH and DELPHI analyses. In the case of the L3 analysis, the uncertainties in the single pion and in the lepton energy spectra were evaluated from the discrepancies between the model of Ref [128] and the ISGW [22] model respectively to the JETSET 7.4 [129] prediction. The exclusive model uncertainty corresponds to the modelling of the hadronic final state in the b → ℓ − ν ℓ X u decay. These uncertainties have been estimated by replacing the parton shower fragmentation model in JETSET [129] with the fully exclusive ISGW2 [23] model by ALEPH and DELPHI, and by propagating a 100% uncertainty on the B → πℓ − ν ℓ rate by L3. ALEPH and L3 have also taken into account the uncertainty from the modelling of the charmless semi-leptonic decay of b-baryons. This has not been considered by DELPHI as they remove decays containing identified protons and kaons, thus suppressing the contribution of b-baryons as mentioned above. In addition to these sources, ALEPH has estimated a b → u uncertainty from the energy cut-off value Λ for the hybrid model adopted [130]; DELPHI allowed for the b-quark pole mass m b and the expectation value of the kinetic energy operator < p 2 b > uncertainties which have been assumed to be uncorrelated.
The three measurements of BR(b → ℓ − ν ℓ X u ) have been averaged using the Best Linear Unbiased Estimate (B.L.U.E.) technique [131]. Using the inputs from Table 9 and Table 10, the LEP average value for BR(b → ℓ − νX u ) was found to be: with a confidence level for the combination of 0.70 (see Figure 10). The value of the |V ub | element has been extracted using the following relationship derived in the context of (OPE) [132,133]: by assuming m b = (4.58 ± 0.06) GeV/c 2 (Appendix E). From the LEP average of BR(b → ℓ − ν ℓ X u ) obtained above, the b-lifetime value, τ b , obtained in Section 3, and the quoted |V ub | uncertainty coming from OPE (±4.7 % relative), a probability density function for |V ub | has been calculated. The resulting distribution is shown in Figure 10-right, where all errors are convoluted together assuming that they are Gaussian in BR(b → ℓ − ν ℓ X u ), with the exception of the OPE error assumed to be Gaussian in |V ub |. The small part of this function in the negative, unphysical region, corresponding to only 0.14%, has been discarded and the probability density function renormalised accordingly. The median of this function has been chosen as the best estimate of |V ub |, and ±34.135% and ±47.725% of the integral of the probability density function around this value have been used to define the 1σ and 2σ confidence regions, denoted hereafter as the 68% and 95% confidence levels, obtaining:  Figure 10: Left: the determinations of BR(b → ℓ − ν ℓ X u ) by ALEPH, DELPHI and L3 and the resulting LEP average. Right: the probability density function for |V ub | corresponding to the LEP average value of BR(b → ℓ − ν ℓ X u ) with the value of the median and two confidence intervals indicated. Unphysical negative entries have been discarded and the probability density function renormalised accordingly.
All the uncertainties have been included in these estimates. The application of the above procedure for each error source separately yields the following detailed result for the 68% confidence level:

Summary of all results
This paper provides precise combined results (see Table 11), from measurements submitted to the 1999 Summer Conferences, on parameters which govern production and decay properties of b-hadrons emitted in high energy b-quark jets. At the Z pole, the polarization of b-baryons is found to be significantly different from zero but it is reduced as compared to the initial b-quark polarization of −0.94. Production of Σ ( * ) b baryons has been invoked [28] to explain this difference. The accuracy on B 0 d and B − lifetimes is better than 2%. The B − lifetime is significantly longer than the B 0 d lifetime, in agreement with the original expectations based on the OPE and parton-hadron duality [71]. But the clear difference between measured and expected lifetimes of b-baryons remains to be explained. It may point to a failure of parton-hadron duality in inclusive B decays, or to a failure of quark models used to evaluate the mean values of operators contributing to b-baryon decays.   Table 11: Summary of the results obtained on b-hadron production rates and decay properties using data available by mid 1999. More details on systematic uncertainties can be found in the corresponding sections.
The accuracy on B 0 d and B − production rates in b-quark jets has improved by a factor two as compared to published values [4]. In conjunction with the better determination of the B 0 d lifetime, this measurement has given a more precise determination of |V cb | from B 0 d → D * + ℓ − ν ℓ decays. The B 0 s fraction in jets is close to 10% with an uncertainty which has been reduced by 30%. Its determination is an important input to searches for B 0 s -B 0 s oscillations. The accuracy on the b-baryon rate has improved by about a factor two because of new measurements using spectator baryons.
The mass difference between mass eigenstates of the B 0 d -B 0 d system is measured with a relative precision better than 4%. The corresponding quantity for the B 0 s -B 0 s system is still unmeasured, in spite of impressive progress in the sensitivity reached by the experiments.
A first combined result on the decay width difference between mass eigenstates of the B 0 s -B 0 s system has been obtained. Because of the very accurate measurements obtained on the inclusive b-lifetime and semileptonic branching fraction, the accuracy of |V cb | determined from inclusive semileptonic decays is entirely limited by theoretical uncertainties. An experimental control of these uncertainties is thus needed. The measurement of |V cb | from exclusive B 0 d → D * + Xℓ − ν ℓ decays is limited by uncertainties related to D * * production and decay mechanisms. Present experimental results on production characteristics of these states in b-hadron semileptonic decays have been detailed in Section 2 of the present report.
The fact that the theoretical uncertainties on the determinations of |V cb | from inclusive and exclusive measurements are largely uncorrelated has been used in evaluating a global average.
The new measurement of |V ub | from LEP experiments has reached an accuracy similar to previous determinations at the Υ(4S) but with a better sensitivity over a larger fraction of the phase space than in measurements using the lepton energy end-point region.
New measurements from LEP and SLD are still expected in the year 2000 and will improve present determinations of b-hadron production and decay properties, which already provide stringent constraints on the shape of the CKM unitarity triangle just before the era of B-factories. semileptonic b-decays 3.73 ± 0.85 ± 0.52 ± 0.24 [24] Table 12: Measured values of the D 1 production rate in b-semileptonic decays. The third systematic uncertainty, quoted in the CLEO analysis, comes from the variation of the detection efficiency when changing the parameters of the signal form factors in the ISGW2 model [23].
5.9 ± 6.6 ± 1.0 ± 0.4 [24] Table 13: Measured values of the D * 2 production rate in b-semileptonic decays. In the original publication from ALEPH, as these values differ from zero only by one or two standard deviations, only upper limits were quoted. The third systematic uncertainty, quoted in the CLEO analysis, comes from the variation of the detection efficiency when changing the parameters of the signal form factors in the ISGW2 model [23].
µ 2 π is the average of the square of the heavy quark momentum inside the B hadron. ∆M 2 includes all terms of order m −n Q , with n ≥ 2 and its absolute value is expected to be smaller than 15 MeV/c 2 . The main difference between the theoretical analyses given in [142] and [144] comes from the evaluation of the uncertainties on non-perturbative corrections. A different approach can be found in [145].
The detailed balance of uncertainties given in [142], which has been revised for consistency with the values adopted at present for the different parameters, is the following: Uncertainties originating from non-perturbative QCD dominate over the perturbative contribution (pert.). The hypotheses or the results on which their contributions were based are reviewed in the following. The label excit. denotes the contribution of excited charmed final states.
The value µ 2 π = (0.5±0.1) GeV 2 has been used and it gives ±0.024 uncertainty on F D * (1). This evaluation has been justified considering that, using QCD sum rules, the following value has been obtained [146]: and that a model-independent lower bound has been established [147]: µ 2 G is the expectation value of the chromomagnetic operator.

E.2.2 High mass excitations
In QCD sum rules, used in [142,144], the effect expected from high mass hadronic states has been introduced in a rather arbitrary way. To be conservative, in [142], it has been assumed that their effect, on the correction terms which behave as m −2 Q , can vary between 0 and 100% of ∆ A 1/m 2 . Such a variation corresponds to ±0.035 on F D * (1).

E.2.3 Higher order non-perturbative corrections
The value of ±0.025 uncertainty related to the contribution from terms of order at least m −3 Q is obtained in [142].

E.2.5 Expected improvements in the control of theoretical errors
The uncertainty attached to "excited states" is expected to be decreased if better measurements of D * * production rates and hadronic mass distributions are obtained 26 .
The uncertainty related to µ 2 π is expected to decrease by measuring moments of the lepton momentum in the B rest frame.
E.2.6 Form factors for D * * production near to w = 1 It has been realized that the w-dependence of the form factors used to describe D * * production in b-hadron semileptonic decays is important for evaluating the systematic uncertainty on |V cb | related to the fraction of D * mesons coming from D * * decays. It has been proposed to use the model described in [21], which predicts a rate for the ratio of 2 + over 1 + narrow states more in agreement with the experimental results given in Section 2. Additional information on D * * production which may bring constraints on such models would be welcome.
E.3 Measurement of |V cb | using the inclusive semileptonic decay b → ℓ − X rate The expression relating |V cb | to the inclusive semileptonic branching fraction can be found in [142]: The central value of reference [142] has been lowered by 2% to account for a different choice of the b-quark mass. A very similar result for the central value and the uncertainties can be found in [133].

E.3.1 Uncertainties related to quark masses
It has been assumed that m b (µ) is known with an uncertainty of ±60 MeV/c 2 , as discussed in Section E.1, and that the difference between the b-and c-quark masses is known with an uncertainty of ±40 MeV/c 2 , related to the error on µ 2 π . These variations induce ±0.010 and ±0.012 uncertainties on |V cb |, respectively. 26 D * * denotes all decay modes which are not B → Dℓ − ν ℓ and B → D * ℓ − ν ℓ .

E.3.2 Adopted value
Adding in quadrature the quoted uncertainties in Equation (87) leads to a ±2.5% relative uncertainty on |V cb |. Since these errors have not been cross-checked by other theoretical approaches or experimental measurements and because most of these uncertainties depend on the evaluation of m c (or µ 2 π ) which is less under control than m b , it has been decided to inflate the total error by an arbitrary factor of two, leading to a theoretical uncertainty of ±5% (see reference [134] for a more detailed discussion). This implies, in practice, that uncertainties of ±120 MeV/c 2 , ±80 MeV/c 2 , and ±0. E.5 Sources of theoretical errors entering into the measurement of the ratio |V ub | |V cb | In the evaluation of the ratio |V ub | |V cb | , several common experimental and model systematics have a reduced effect.
For the ratio |V ub | incl |V cb | incl , uncertainties originating from the determination of m b have to be considered as fully correlated. Leading effects of order O(1/m 3 ) are expected to be uncorrelated between |V ub | and |V cb | determinations [134]. This is also expected for perturbative uncertainties.

E.6 Conclusions and summary
The conclusions of the workshop which are given in this Appendix have left aside several subtleties concerning the exact meaning of the parameters entering into the different expressions for |V ub | (incl.), |V cb | (incl.) and |V cb | (excl.). The interested reader has therefore to consult the documents quoted in the references for more information.
Two groups (at least) have obtained consistent results on central values and uncertainties for |V ub | (incl.) and |V cb | (incl.).
Uncertainties on |V cb | (incl.) have been enlarged, in a rather arbitrary way, by a factor two to have some margin because of possible additional contributions (reliability of 1/m c expansion, need for other techniques to evaluate m b ,...).
Dedicated experimental studies on D * * production and on the distributions of moments of the lepton momentum in the B rest frame are needed to evaluate with better accuracy and greater confidence the most important sources of systematic errors contributing in the theoretical expressions.
The following central values and theoretical uncertainties have been used in the measurements of |V ub | and |V cb | obtained from combined LEP analyses: All uncertainties have been added in quadrature. The correlated uncertainty between |V cb | (incl.) and |V cb | (excl.) measurements stem from the effect of µ 2 π . The correlated uncertainty between |V ub | (incl.) and |V cb | (incl.) measurements is limited to the effect of m b . Hence no correlated uncertainty between |V ub | (incl.) and |V cb | (excl.) measurements is assumed.