Measurement of |Vcb| using B0 to D*l nu decays

The magnitude of the Cabibbo-Kobayashi-Maskawa matrix element Vcb has been measured using B0 to D*l nu decays recorded on the Z0 peak using the OPAL detector at LEP. The D* to D0 pi+ decays were reconstructed both in the particular decay modes D0 to K- pi+ and D0 to K- pi+ pi- and via an inclusive technique. The product of |Vcb| and the decay form factor of the B0 to D* l nu transition at zero recoil F(1) was measured to be F(1)|Vcb| = (37.1+-1.0+-2.0)x10^-3, where the uncertainties are statistical and systematic respectively. By using Heavy Quark Effective Theory calculations for F(1), a value of |Vcb| = (40.7+-1.1+-2.2+-1.6)x 10^-3 was obtained, where the third error is due to theoretical uncertainties in the value of F(1). The branching ratio Br(B0 to D* l nu) was also measured to be (5.26+-0.20+-0.46)%.


Introduction
Within the framework of the Standard Model of electroweak interactions, the elements of the Cabibbo-Kobayashi-Maskawa mixing matrix are free parameters, constrained only by the requirement that the matrix be unitary.The values of the matrix elements can only be determined by experiment.Heavy Quark Effective Theory (HQET) provides a means to extract the magnitude of the element V cb from particular semileptonic b decays, with relatively small theoretical uncertainties [1].
In this paper, the value of |V cb | is extracted by studying the decay rate for the process B0 → D * + ℓ − ν as a function of the recoil kinematics of the D * + meson1 [2][3][4].The decay rate is parameterised as a function of the variable ω, defined as the scalar product of the four-velocities of the D * + and B0 mesons.This is related to the square of the four-momentum transfer from the B0 to the ℓ − ν system, q 2 , by and ranges from 1, when the D * + is produced at rest in the B0 rest frame, to about 1.50.Using HQET, the differential partial width for this decay is given by where K(ω) is a known phase space term and F(ω) is the hadronic form factor for this decay [3].
Although the shape of this form factor is not known, its magnitude at zero recoil, ω = 1, can be estimated using HQET.In the heavy quark limit (m b → ∞), F(ω) coincides with the Isgur-Wise function [4] which is normalised to unity at the point of zero recoil.Corrections to F(1) have been calculated to take into account the effects of finite quark masses and QCD corrections, yielding the value and theoretical uncertainty F(1) = 0.913 ± 0.042 [5].Since the phase space factor K(w) tends to zero as ω → 1, the decay rate vanishes at ω = 1 and the accuracy of the extrapolation relies on achieving a reasonably constant reconstruction efficiency in the region near ω = 1.Previous measurements of |V cb | have been made using B mesons produced on the Υ(4S) resonance [6] and in Z 0 decays [7][8][9].These analyses used a linear or constrained quadratic expansion of F(ω) around ω = 1.An improved theoretical analysis, based on dispersive bounds and including higher order corrections, has since become available [10].This results in a parameterisation for F(ω) in terms of F(1) and a single unknown parameter ρ 2 constrained to lie in the range −0.14 < ρ 2 < 1.54, ρ 2 corresponding to the slope of F(ω) at zero recoil.
The previous OPAL measurement [9] used the decay chain D * + → D 0 π + , with the D 0 meson being reconstructed in the exclusive decay channels D 0 → K − π + and D 0 → K − π + π 0 .In this paper, a new analysis is described in which only the π + from the D * + decay is identified, and no attempt is made to reconstruct the D 0 decay exclusively.This technique, first employed by DELPHI [8,11], gives a much larger sample of B0 → D * + ℓ − ν decays than the previous measurement, but also larger background, requiring a rather more complex analysis.The measurement of [9] is also updated to use the new parameterisation of F(ω) [10], and improved background models and physics inputs.In both cases, the initial number of B 0 mesons is determined from other measurements of B 0 production in Z 0 decays.
The new reconstruction technique is described in Section 2, the determination of ω for each event in Section 3, the fit to extract F(1)|V cb | and ρ 2 in Section 4 and the systematic errors in Section 5.The updated exclusive measurement is discussed in Section 6 and the measurements are combined and conclusions drawn in Section 7.
2 Inclusive reconstruction of B0 → D * + ℓ − ν events The OPAL detector is well described elsewhere [12].The data sample used in this analysis consists of about 4 million hadronic Z 0 decays collected during the period 1991-1995, at centre-of-mass energies in the vicinity of the Z 0 resonance.Corresponding simulated event samples were generated using JETSET 7.4 [13] as described in [14].
Hadronic Z 0 decays were selected using standard criteria [14].To ensure the event was well contained within the acceptance of the detector, the thrust axis direction2 was required to satisfy | cos θ T | < 0.9.Charged tracks and electromagnetic calorimeter clusters with no associated tracks were then combined into jets using a cone algorithm [15], with a cone half angle of 0.65 rad and a minimum jet energy of 5 GeV.The transverse momentum p t of each track was defined relative to the axis of the jet containing it, where the jet axis was calculated including the momentum of the track.A total of 3 117 544 events passed the event selection.
The reconstruction of B0 → D * + ℓ − ν events was then performed by combining high p and p t lepton (electron or muon) candidates with oppositely charged pions from the D * + → D 0 π + decay.Electrons were identified and photon conversions rejected using neural network algorithms [14], and muons were identified as in [16].Both electrons and muons were required to have momenta p > 2 GeV, transverse momenta with respect to the jet axis p t > 0.7 GeV, and to lie in the polar angle region | cos θ| < 0.9.
The event sample was further enhanced in semileptonic b decays by requiring a separated secondary vertex with decay length significance L/σ L > 2 in any jet of the event.The vertex reconstruction algorithm and decay length significance calculation are described fully in [14].Together with the lepton selection, these requirements result in a sample which is about 90 % pure in bb events.
An attempt was made to estimate the D 0 direction in each jet containing a lepton candidate.Each track (apart from the lepton) and calorimeter cluster in the jet was assigned a weight corresponding to the estimated probability that it came from the D 0 decay.The track weight was calculated from an artificial neural network, trained to separate tracks from b decays and fragmentation tracks in b jets [14].The network inputs are the track momentum, transverse momentum with respect to the jet axis, and impact parameter significances with respect to the reconstructed primary and secondary vertices (if existing).The cluster weights were calculated using their energies and angles with respect to the jet axis alone, the energies first being corrected by subtracting the energy of any charged tracks associated to the cluster [17].
Beginning with the track or cluster with the largest weight, tracks and clusters were then grouped together until the invariant mass of the group (assigning tracks the pion mass and clusters zero mass) exceeded the charm hadron mass, taken to be 1.8 GeV.If the final invariant mass exceeded 2.3 GeV, the jet was rejected, since Monte Carlo studies showed such high mass D 0 candidates were primarily background.For surviving jets, the momentum p D 0 of the group was used as an estimate of the D 0 direction, giving RMS angular resolutions of about 45 mrad in φ and θ.The D 0 energy was calculated as The selection of pions from D * + decays relies on the small mass difference of only 145 MeV [18] between the D * + and D 0 , which means the pions have very little transverse momentum with respect to the D 0 direction.Each track in the jet (other than the lepton) was considered as a slow pion candidate, provided it satisfied 0.5 < p < 2.5 GeV and had a transverse momentum with respect to the D 0 direction of less than 0.3 GeV.If the pion under consideration was included in the reconstructed D 0 , it was removed and the D 0 momentum and energy recalculated.The final selection was made using the reconstructed mass difference ∆m between D * + and D 0 mesons, calculated as where the D * + energy is given by E D * = E D 0 + E π and momentum by p D * = p D 0 + p π .
A new secondary vertex was then iteratively reconstructed around an initial seed vertex formed by the intersection of the lepton and slow pion tracks.Every other track in the jet was added in turn to the seed vertex, and the vertex refitted.The track resulting in the lowest vertex fit χ 2 was retained in the seed vertex for the next iteration.The procedure was repeated until no more tracks could be added without reducing the vertex fit χ 2 probability to less than 1 %.The decay length L ′ between the primary vertex and this secondary vertex, and the associated error σ L ′ , were calculated as in [14].The pion candidate was accepted if the decay length satisfied −0.1 cm < L ′ < 2 cm and the decay length significance satisfied The resulting distributions of ∆m for opposite and same sign lepton-pion combinations are shown in Figure 1(a) and (b).The predictions of the Monte Carlo simulation are also shown, broken down into contributions from signal B0 → D * + ℓ − ν events, 'resonant' background containing real leptons and slow pions from D * + decay, and combinatorial background, made up of events with fake slow pions, fake leptons or both.
In Monte Carlo simulation, about 45 % of opposite sign events with ∆m < 0.17 GeV are signal B0 → D * + ℓ − ν events, 14 % are resonant background and 41 % are combinatorial background.The resonant background is made up mainly of B − → D * + π − ℓ − ν, B0 → D * + π 0 ℓ − ν and Bs → D * + K 0 ℓ − ν decays.These are expected to be dominated by b semileptonic decays involving orbitally excited charm mesons (generically referred to as D * * ), e.g.B − → D * * 0 ℓ − ν followed by D * * 0 → D * + π − .These decays will be denoted collectively by B → D * + h ℓ − ν.Small contributions are also expected from b → D * + τ νX decays with the τ decaying leptonically, and b → D * + D − s X with the D − s decaying semileptonically (each about 1 % of opposite sign events).For same sign events with ∆m < 0.17 GeV, there is a small resonant contribution of about 6 % from events with a real D * + → D 0 π + where the D 0 decays semileptonically, and the rest is combinatorial background.
The most important background, from B → D * + h ℓ − ν decays, comes from both charged B + and neutral B 0 and B s decays, whereas the signal comes only from B 0 decays.Therefore the charge Q vtx of the reconstructed secondary vertex containing the lepton and slow pion and its estimated error σ Qvtx were calculated, using where w i is the weight for track i of charge q i to come from the secondary vertex, and the sums are taken over all tracks in the jet [19,20].The weights were calculated in a similar way to those used for the D 0 direction reconstruction, using a neural network with the track momentum, transverse momentum and impact parameter significances with respect to the reconstructed primary and secondary vertices as inputs.The weights for the lepton and slow pion candidate tracks were set to one.The vertex charge distributions for opposite sign events with ∆m < 0.17 GeV are shown in Figure 1(c) and (d).
The reconstructed vertex charge and error were used to divide the data into different classes enhanced or depleted in B + decays, thus reducing the effect of this background and increasing the statistical sensitivity.Three classes c were used-class 1 where the charge is measured poorly (σ Qvtx > 0.9), class 2 where the charge is measured well and is compatible with a neutral vertex (σ Qvtx < 0.9, |Q vtx | < 0.5) and class 3 where the charge is measured well and is compatible with a charged vertex (σ Qvtx < 0.9, |Q vtx | > 0.5).

Reconstruction of ω
The recoil variable ω was estimated in each event using the reconstructed four-momentum transfer to the ℓν system: together with equation 1.The B 0 and D * + energies E B 0 and E D * were estimated directly, whilst the momentum vectors p B 0 and p D * were estimated using the energies together with the reconstructed polar and azimuthal angles, as described in more detail below.Since the slow pion has very little momentum in the rest frame of the decaying D * + , the momentum (and hence energy) of the D * + was estimated by scaling the reconstructed slow pion momentum by m D * + /m π , as for the D 0 → K − π + π 0 channel in [9].A small correction (never exceeding 12 %) was applied, as a function of cos θ * , the angle of the slow pion in the rest frame of the D 0 .This procedure gave a fractional D * + energy resolution of 15 %.The polar and azimuthal angles of the D * + were reconstructed using weighted averages of the slow pion and D 0 directions, giving resolutions of about 22 mrad on both φ and θ.
The energy of the B 0 was estimated using a technique similar to that described in [21], exploiting the overall energy and momentum conservation in the event to calculate the energy of the unreconstructed neutrino.First, the energy E bjet of the jet containing the B 0 was inferred from the measured particles in the rest of the event, by treating the event as a two-body decay of a Z 0 into a b jet (whose mass was approximated by the B 0 mass) and another object making up the rest of the event.Then, the total fragmentation energy E frag in the b jet was estimated from the measured visible energy in the b jet and the identified B 0 decay products: Finally, the B 0 energy was calculated as E B 0 = E bjet − E frag , giving an RMS resolution of 4.4 GeV.
The b direction was estimated using a combination of two techniques.In the first, the momentum vector of the B 0 was reconstructed from its decay products: the neutrino energy being estimated from the reverse of the event visible momentum vector: p ν = −p vis .The visible momentum was calculated using all the reconstructed tracks and clusters in the event, with a correction for charged particles measured both in the tracking detectors and calorimeters [17].This direction estimate is strongly degraded if a second neutrino is present (e.g. from another semileptonic b decay in the opposite hemisphere of the event), and its error was parameterised as a function of the visible energy in the opposite hemisphere.The resulting resolution is typically between 40 and 100 mrad for both φ and θ.
The second method of estimating the B 0 direction used the vertex flight direction-i.e. the direction of the vector between the primary vertex and the secondary vertex reconstructed around the lepton and slow pion as described in Section 2. The accuracy of this estimate is strongly dependent on the B 0 decay length, and was used only if the decay length significance L ′ /σ L ′ exceeded 3.After this cut, the angular resolution varies between about 15 and 100 mrad, and is worse for θ in the 1991 and 1992 data where accurate z information from the silicon microvertex detector was not available.The B0 → D * + ℓ − ν candidate was rejected if the two reconstruction methods gave θ or φ angles disagreeing by more than three standard deviations, which happened in 7 % of Monte Carlo signal events.Finally, the two B 0 direction estimates were combined according to their estimated uncertainties, giving average resolutions of 35 mrad on φ and 43 mrad on θ, including events where only the first method was used.
The estimate of q 2 derived from the B 0 and D * + energies and angles was improved by applying the constraint that the mass µ of the neutrino produced in the B0 → D * + ℓ − ν decay should be zero.The neutrino mass is given from the reconstructed quantities by: The constraint was implemented by calculating q 2 using a kinematic fit, incorporating the measured values and estimated uncertainties of the B 0 and D * + energies and angles (the B 0 angular uncertainties varying event by event).This procedure improved the average q 2 resolution from 2.78 GeV 2 to 2.57 GeV 2 .In Monte Carlo simulation, 11 % of signal events were reconstructed with q 2 < 0, corresponding to an unphysical value of ω larger than ω max ≈ 1.5, and were rejected.
The resulting distributions of reconstructed ω for various ranges of true ω (denoted ω ′ ) in simulated B0 → D * + ℓ − ν decays are shown in Figure 2(a-e).The average RMS ω resolution is about 0.12, but there are significant non-Gaussian tails.The resolution was parameterised (separately for the 1991-2 and 1993-5 data) as a continuous function R(ω, ω ′ ) giving the expected distribution of reconstructed ω for each true value ω ′ .The resolution function was implemented as the sum of two asymmetric Gaussians (i.e. with different widths either side of the peak) whose parameters were allowed to vary as a function of ω ′ .The convolution of this resolution function with the Monte Carlo ω ′ distribution is also shown in Figure 2(a-e), demonstrating that the resolution function models the ω distributions well.
The efficiency to reconstruct B0 → D * + ℓ − ν decays ǫ(ω ′ ) is shown in Figure 2(f), together with a second order polynomial parameterisation.The efficiency varies with ω ′ , but is reasonably flat in the critical region near ω ′ = 1 where the extrapolation to measure F(1)|V cb | is carried out.

Fit and results
The values of F(1)|V cb | and ρ 2 were extracted using an extended maximum likelihood fit to the reconstructed mass difference ∆m, recoil ω and vertex charge class c of each event.Both opposite and same sign events with ∆m < 0.3 GeV and ω < ω max were used in the fit, the high ∆m and same sign events serving to constrain the combinatorial background normalisation and shapes in the opposite sign low ∆m region populated by the B0 → D * + ℓ − ν decays.Using the ∆m value from each event in the fit, rather than just dividing the data into low ∆m 'signal' and high ∆m 'sideband' mass regions, increases the statistical sensitivity as the signal purity varies considerably within the low ∆m region.
The logarithm of the overall likelihood was given by where the sums of individual event log-likelihoods ln L a i and ln L b j are taken over all the observed M a opposite sign and M b same sign events in the data sample, and N a and N b are the corresponding expected numbers of events.
The likelihood for each opposite sign event was given in terms of different types or sources of event by where N a s is the number of expected events for source s, f a s,c is the fraction of events in source s appearing in vertex charge class c, M s,c (∆m) is the mass difference distribution for source s in class c and P s (ω) the recoil distribution for source s.For each source, the vertex charge fractions f a s,c sum to one and the mass difference M s,c (∆m i ) and recoil P s (ω) distributions are normalised to one.The total number of expected events is given by the sum of the individual contributions: N a = 4 s=1 N a s .There are four opposite sign sources: (1) signal B0 → D * + ℓ − ν events, (2) B → D * + h ℓ − ν events where the D * + is produced via an intermediate resonance (D * * ), (3) other opposite sign background involving a genuine lepton and a slow pion from D * + decay and (4) combinatorial background.The sum of sources 2 and 3 are shown as 'resonant background' in Figure 1.A similar expression to equation 4 was used for L b j , the event likelihood for same sign events.In this case, only sources 3 and 4 contribute.
The mass difference distributions M c,s (∆m) for sources 1-3 were represented by analytic functions, whose parameters were determined using large numbers of simulated events, as were the recoil distributions P s (ω) for sources 2 and 3.The fractions in each vertex charge class for sources 1-3 were also taken from simulation.For the signal (source 1), the product of the expected number of events N a  (0.76 ± 0.16) % [22], see text Br(b → D * + τ − νX) (0.65 ± 0.13) % [18], see text Br( B0 (13.9 ± 0.9) % [18] Table 1: Input quantities used in the fits for F(1)|V cb | and ρ 2 .The values marked 'see text' are derived using methods explained in Section 4.
and recoil distribution P 1 (ω) was given by convolving the differential partial decay width (equation 2) with the signal resolution function and reconstruction efficiency: where N Z is the number of hadronic Z 0 decays passing the event selection, R b ≡ Γ b b/Γ had is the fraction of hadronic Z 0 decays to bb, f B 0 the fraction of b quarks hadronising to a B0 and τ B 0 the B 0 lifetime.The factor of four accounts for the two b hadrons produced per Z 0 → bb event and the two identified lepton species (electrons and muons).The form factor F(ω ′ ) is given in [10] in terms of the normalisation F(1) and slope parameter ρ 2 .The efficiency function ǫ(ω ′ ) and resolution function R(ω, ω ′ ) were described in Section 3, and the known phase space factor K(ω ′ ) is given in [9].The assumed values of the numerical quantities are given in Table 1.
Since the data are divided into different vertex charge classes enhanced and depleted in B + decays, the fit gives some information on the amount of B → D * + h ℓ − ν background.The predicted level of this background in the fit was therefore allowed to vary under a Gaussian constraint corresponding to the branching ratio of Br(b → D * + h ℓν) = (0.76 ± 0.16) %.The latter has been calculated from the measured branching ratio Br(b → D * + π − ℓνX) = (0.473 ± 0.095) % [22], assuming isospin and SU(3) flavour symmetry to obtain the corresponding b → D * + π 0 ℓν and b → D * + K 0 ℓν branching ratios.A scaling factor of 0.75 ± 0.25 was included for the last branching ratio to account for possible SU(3) violation effects reducing the branching ratio D * * + s → D * + K 0 compared to the expectation of 3 2 Br(D * * + → D * + π 0 ).The P 2 (ω) distribution for these events was taken from simulation, using the calculation of Leibovich et al. [23] to predict their recoil spectrum.
The numbers N a,b 3 and P 3 (ω) distributions for the small background contributions covered by source 3 (both opposite and same sign) were taken from Monte Carlo simulation, with branching ratios adjusted to the values given in Table 1.The branching ratio for b → D * + τ − ν was derived using the inclusive branching ratio Br(b → τ X) = (2.6 ± 0.4) % [18], assuming a D * + is produced in a fraction 0.25 ± 0.03 of the time, as seen for the corresponding decays b → D * + ℓX and b → ℓX (ℓ = e, µ) [18].The rate of b → D * + D − s X was assumed to be dominated by the two body decay B0 → D * + D ( * )− s .The rate of real D * + decays in the same sign background depends on the production fractions of D * + in bb and cc events, which were taken from [24].
The parameters of the analytic functions describing the combinatorial background ∆m) and P 4 (ω)) were fitted entirely from the data, with only the choice of functional forms motivated by simulation.The shapes of the mass and recoil functions (including a small correlation between ∆m and ω) are constrained by the same sign sample (which is almost entirely combinatorial background), and are the same for each vertex charge class c.The opposite sign high ∆m region serves to normalise the number of combinatorial background events in the low ∆m region for each vertex charge class.
The values of F(1)|V cb | and ρ 2 were extracted by maximising the total likelihood given by equation 3. The values of F(1)|V cb | and ρ 2 were allowed to vary, together with the level of B → D * + h ℓ − ν background and 13 auxiliary parameters describing the combinatorial background distributions.A result of F(1)|V cb | = (37.5 ± 1.2) × 10 −3 , ρ 2 = 1.12 ± 0.14 was obtained, where the errors are only statistical.The correlation between F(1)|V cb | and ρ 2 is 0.77.The distributions of reconstructed ω for opposite and same sign events with ∆m < 0.17 GeV, together with the fit results, are shown in Figure 3.The fit describes the data well, both in this region and the high ∆m region dominated by combinatorial background, for all three of the vertex charge classes.
By integrating the differential partial decay width (equation 2) over all values of ω, the branching ratio Br( B0 → D * + ℓ − ν) was also determined to be where again the error is statistical only.This result is consistent with the world average once systematic errors are included.
Many previous results have been obtained using a constrained quadratic expansion for the form factor: 2 , where a is a slope parameter to be determined by the fit, and b is constrained to b = 0.66a 2 − 0.11 [25].To allow comparison with such measurements, the fit was also performed with this parameterisation of F(ω), giving the results F(1)|V cb | = (36.9±1.2)×10−3 and a 2 = 0.88 ± 0.14, the correlation between F(1)|V cb | and a 2 being 0.79.The difference in the two curvature parameters ρ 2 and a 2 is in good agreement with the expectation of ρ 2 − a 2 ≈ 0.21 [10].

Systematic Errors
Systematic errors arise from the uncertainties in the fit input parameters given in Table 1, the Monte Carlo modelling of the signal ω resolution, the recoil spectrum of b → D * * ℓν decays and selection efficiencies, and possible biases in the fitting method.The resulting systematic errors on the values of F(1)|V cb |, ρ 2 and Br( B0 → D * + ℓ − ν) are summarized in Table 2 and described in more detail below.
Input quantities: The various numerical fit inputs were each varied according to the errors given in Table 1 and the fit repeated to assess the resulting uncertainties.
b → D * * ℓν decays: The calculation of Leibovich et al. [23] was used to simulate the recoil spectrum of B → D * + h ℓ − ν decays, assumed to be produced via the semileptonic decay B 0 → D * * ℓν.
Here D * * represents a P-wave orbitally excited charm meson.The calculation predicts the recoil spectra and relative rates of semileptonic decays involving both the narrow D 1 and D * 2 states and wide D * 0 and D * 1 states.All these decays are suppressed close to ω = 1 by an extra factor of (ω 2 − 1) when compared with the signal B0 → D * + ℓ − ν decays.This reduces the uncertainty due to the rate of B → D * + h ℓ − ν decays in the extrapolation of the signal recoil spectrum to ω = 1.Non-resonant B → D * + h ℓ − ν decays are not included in the model, but are not expected to contribute close to ω = 1.The differential decay rates in [23] are given in terms of five possible expansion schemes and several unknown parameters: a kinetic energy term η ke and the slopes of the Isgur-Wise functions for the narrow and wide D * * states τ1 and ζ1 .These parameters were varied within the allowed  Error Source ∆(F(1) for the inclusive analysis.The fractional errors on F(1)|V cb | and Br( B0 → D * + ℓ − ν) are given, whereas the errors on ρ 2 are absolute.
ranges −0.75 < η ke < 0.75 GeV, −2 < τ1 < −1 and −2 < ζ1 < 0 subject to the constraint that the ratio R = Γ( B → D * 2 ℓν)/(Γ( B → D 1 ℓν) lie within the measured range R = 0.37 ± 0.16 [22,26].This excludes the expansion schemes A ∞ and B ∞ of [23] and constrains the allowable values of η ke in the others.The fraction of B → D * * ℓ − ν decays involving the narrow D 1 and D * 2 states, which is not precisely predicted, was varied in the range 0.22 ± 0.06, obtained by comparing the measured rates for B 0 semileptonic decays involving D + , D * + , D 1 and D * 2 with the inclusive semileptonic decay rate [18,22,26].The systematic errors were determined as half the difference between the two parameter sets giving the most extreme variations in F(1)|V cb | and ρ 2 , and the central values were adjusted to half way between these two extremes.The values of both F(1)|V cb | and ρ 2 are most sensitive to variations in η ke , which is constrained by the measured value of R.
b fragmentation: The effect of uncertainties in the average b hadron energy x E = E b /E beam was assessed in Monte Carlo simulation by reweighting the events so as to vary x E in the range 0.702 ± 0.008 [28], and repeating the fit.The largest of the variations observed using the fragmentation functions of Peterson, Collins and Spiller, Kartvelishvili and the Lund group [27] were taken as systematic errors.
D 0 decay multiplicities: The signal reconstruction efficiency and vertex charge distributions are sensitive to the B 0 decay multiplicity, which depends only on the D 0 decay for the B 0 decay channels of interest.The systematic error was assessed by varying separately the D 0 charged and π 0 decay multiplicities in Monte Carlo simulation according to the measurements of Mark III [29].The branching ratio D 0 → K 0 , K0 was also varied according to its uncertainty [18].The resulting uncertainties on F(1)|V cb | and ρ 2 from each variation were added in quadrature to determine the total systematic errors.
ω reconstruction: The modelling of the ω resolution depends on the description of the D * + and B 0 energy and angular distributions in the simulation.The reconstructed D * + and B 0 energy distributions in data and simulation were compared, and the means were found to differ by 0.04 and 0.13 GeV respectively.The opposite hemisphere missing energy was found to agree within 5 %.The corresponding systematics were assessed by shifting or scaling the data distributions and repeating the ω reconstruction and fit.The modelling of the angular resolution was checked by studying the agreement of the two angular estimators-the slow pion and D 0 directions for the the D * + , and the missing energy vector and vertex flight directions for the B 0 .The angular resolutions were found to be up to 5 % worse in data, and the systematic error was assessed by degrading the simulated resolution appropriately.Finally, the fraction of events with ω reconstructed in the physical region (ω < ω max ) was found to be 3.5 % smaller in the data, in both opposite and same sign charge samples.The reconstruction efficiency was corrected for this effect, and an additional systematic error of half the correction (1.7 %) assumed.The final systematic errors due to ω resolution modelling are dominated by the B 0 θ resolution.
Lepton identification efficiency: The electron identification efficiency has been studied using control samples of pure electrons from e + e − → e + e − events and photon conversions, and found to be modelled to a a precision of 4.1 % [14].The muon identification efficiency has been studied using muon pairs producted in two-photon collisions and Z 0 → µ + µ − events, giving an uncertainty of 2.1 % [30].
Vertex tag efficiency: The fraction of hemispheres with identified leptons which also had a selected secondary vertex was found to be about 4 % less in data than in simulation.The overall fraction of vertex tagged hemispheres was also found to be about 4 % lower in data.These discrepancies were translated into systematic errors on the efficiency to tag a semileptonic b decay with a secondary vertex in either the same or the opposite hemisphere, in each case attributing the whole discrepancy to a mismodelling of b hadron decays.The resulting errors on the same and opposite hemisphere tagging efficiencies were taken to be fully correlated.
Track reconstruction efficiency: The overall track reconstruction efficiency is known to be modelled to a precision of 1 % [14], and a similar uncertainty was found to be appropriate for the particular class of slow pion tracks from D * + decays.The systematic error was assessed by randomly removing 1 % of tracks in the simulation and repeating the fit.
Tracking resolution: Uncertainties in the tracking detector resolution affect the efficiency, ω reconstruction and vertex charge distributions.The associated error was assessed in the simulation by applying a global 10 % degradation to all tracks, independently in the r-φ and r-z planes, as in [14].To verify the correctness of the statistical errors returned by the fit, it was performed on many separate subsamples, and the distribution of fitted results studied.Further checks on the data included performing the analysis separately for B 0 decays involving electrons and muons, dividing the sample according to the year of data taking, and varying the lepton transverse momentum cut.In all cases, consistent results were obtained.
Including all systematic uncertaintities, the final result of the inclusive analysis is where the first error is statistical and the second systematic in each case.
6 Measurement using exclusively reconstructed D * + decays In this analysis, the D 0 from the B0 → D * + ℓ − ν, D * + → D 0 π + decay is reconstructed exclusively in the decay modes D 0 → K − π + ('3-prong') and D 0 → K − π + π 0 ('satellite'-where the π 0 is not reconstructed).The event selection, reconstruction and determination of ω are exactly the same as described in [9].The determination of the signal and background fractions, and the fit to extract F(1)|V cb | and ρ 2 are performed as in [9], but have been updated using the improved form factor parameterisation [10], the updated input parameters given in Table 1 and the b → D * * ℓν decay model [23] discussed in Section 5.
The fit is performed on 814 3-prong and 1396 satellite candidates, of which 505 ± 44 and 754 ± 72 are attributed to B0 → D * + ℓ − ν signal decays.The result of the fit is where again the first errors are statistical and the second systematic.The statistical correlation between F(1)|V cb | and ρ 2 is 0.95.The distribution of reconstructed ω for selected candidates (both 3-prong and satellite) is shown in Figure 4.The branching ratio Br( B0 → D * + ℓ − ν) has also been determined to be Br( B0 → D * + ℓ − ν) = (5.11± 0.19 ± 0.49) %.
The systematic errors arise from uncertainties in the background levels in the selected samples, as well as uncertainties in the Monte Carlo simulations.They have been evaluated following the procedures described in [9], and are summarised in Table 3.The selection efficiency error includes contributions from lepton identification efficiency, b fragmentation and detector resolution uncertanties, described in detail in [31].The largest change with respect to the previous result is due to the improved b → D * * ℓν modelling, with the suppression of this background at values of ω close to one.This reduces the statistical error, the systematic uncertainty due to the rate of such decays, and shifts the central value of F(1)|V cb | upwards as compared to [9].

Conclusions
The CKM matrix element |V cb | has been measured by studying the rate of the semileptonic decay B0 → D * + ℓ − ν as a function of the recoil kinematics of both inclusively and exclusively reconstructed D * + mesons.The two results are combined, taking into account the statistical correlation of 18 % and correlated systematic errors from physics inputs and detector resolution.The results are: where the first result is statistical and the second systematic in each case.The statistical and systematic correlations between F(1)|V cb | and ρ 2 are 0.90 and 0.54 respectively.These results supersede our previous publication [9].They are consistent with other determinations of F(1)|V cb | at LEP [7,8] and the Υ(4S) resonance [6].The branching ratio is consistent with the world average result of (4.60 ± 0.27) % [18].The result for F(1)|V cb | is the most precise to date from any single experiment.where the uncertainties are statistical, systematic and theoretical respectively.

Figure 1 :
Figure 1: Reconstructed ∆m distributions for selected (a) opposite sign and (b) same sign lepton-pion combinations; reconstructed vertex charge Q vtx distributions for opposite sign events with ∆m < 0.17 GeV and (c) σ Qvtx > 0.9 and (d) σ Qvtx < 0.9.In each case the data are shown by the points with error bars, and the Monte Carlo simulation contributions from signal B0 → D * + ℓ − ν decays, other resonant D * + decays and combinatorial background are shown by the open, single and cross hatched histograms respectively.The vertex charge classes 1 (poorly reconstructed), 2 (well reconstructed neutral vertex) and 3 (well reconstructed charged vertex) are indicated.

Figure 2 :
Figure 2: Reconstructed ω distributions for various ranges of true ω (denoted ω ′ ) in Monte Carlo B0 → D * + ℓ − ν events, together with the prediction from the resolution function.Events reconstructed in the unphysical region with ω > ω max are rejected.The reconstruction efficiency as a function of true ω is also shown.

Figure 3 :
Figure 3: Distributions of reconstructed ω for (a) opposite sign and (b) same sign events with ∆m < 0.17 GeV.The data are shown by the points with error bars and the expectation from the fit result by the histograms.The contributions from signal B0 → D * + ℓ − ν, resonant and combinatorial backgrounds are indicated.

Fit method:
The entire fitting procedure was tested on a fully simulated Monte Carlo sample seven times bigger than the data, with true values of F(1)|V cb | = 32.5 × 10 −3 and ρ 2 = 1.3.The fit gave the results F(1)|V cb | = (31.8± 0.5) × 10 −3 and ρ 2 = 1.25 ± 0.07, consistent with the true values.For each variable, the larger of the deviation of the result from the true value and the statistical error were taken as a systematic errors due to possible biases in the fit.Additionally, the large Monte Carlo sample was reweighted to change the values of F(1)|V cb | and ρ 2 , and the fit correctly recovered the modified values.

Figure 4 :Table 3 :
Figure 4: Distribution of reconstructed ω for B0 → D * + ℓ − ν candidates in the exclusive analysis.The data are shown by the points with statistical error bars, and the fit result by the histogram.The contributions from signal B0 → D * + ℓ − ν, resonant and combinatorial backgrounds are indicated.