Measurement of |Vcb| using B0 ->D* l nu_l Decays

A preliminary measurement of |Vcb| and the branching fraction BR(B0 ->D* l nu_l) has been performed based on a sample of about 55,700 B0 ->D* l nu_l decays recorded with the BaBar detector. The decays are identified in the D* ->D0 pi final state, with the D0 reconstructed in three different decay modes. The differential decay rate is measured as a function of the relativistic boost of the D* in the B0 rest frame. The value of the differential decay rate at `zero recoil', namely the point at which the D* is at rest in the B0 frame, is predicted in Heavy Quark Effective Theory as a kinematic factor times F(1)|Vcb|, where F(1) is the unique form factor governing the decay. We extrapolate the measured differential decay rate to the zero recoil point and obtain F(1)|Vcb|=(34.03+-0.24+-1.31) x 10^-3. Using a theoretical calculation for F(1) we extract |Vcb| = (37.27+-0.26(stat.)+-1.43(syst.)+1.5-1.2(theo.) x 10^-3. From the integrated decay rate we obtain BR(B0 ->D* l nu_l) = (4.68+-0.03+-0.29)%.


Introduction
In the framework of the Standard Model, the weak coupling between quarks of different flavors is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix.Its elements are not predicted by theory, but are constrained only to give a unitary matrix.
A precise measurement of |V cb |, the element corresponding to the b → c transitions, is needed to determine whether the CKM matrix provides a quantitatively accurate description of the CP violation observed in B 0 mesons.Progress in the phenomenological description of heavy flavor semileptonic decay allows the determination of |V cb | with small theoretical uncertainty, either from the inclusive process b → cℓν, or from an analysis of the form factors in the decay B0 → D * + ℓ − νℓ .The present measurement is based on the second approach [1].
The decay rate for the process is proportional to |V cb | 2 and to the square of the hadronic matrix elements describing the transition from a B 0 to a D * + meson.In the context of the Heavy Quark Effective Theory (HQET) [2], the matrix element is proportional to a single form factor F(w), where w is the product of the B 0 and D * + four-vector velocities.For w = 1, the D * + is produced at rest in the B 0 rest frame.Heavy-flavor symmetry predicts the normalization F(1) = 1 in the limit of an infinitely massive b-quark.Corrections to this prediction due to perturbative QCD have been computed up to second order in α S [3] and the effect of finite b and c quark masses has been calculated in the framework of HQET, yielding F(1) = 0.913 +0.030 −0.035 [4].A measurement of the differential decay rate near w = 1 determines |V cb | with small theoretical uncertainty.The rate at w = 1 is suppressed by phase space.Consequently, the differential rate dΓ/dw is measured, where F(w) is parametrized using several different functional forms [5,6] (see also discussion below), and is extrapolated to w = 1.Results based on this approach have been reported by the ARGUS [7], Belle [8], and CLEO [9] collaborations operating at the Υ (4S) and by ALEPH [10,11], DELPHI [12,13] and OPAL [14] at LEP.This paper presents preliminary results for the measurements of |V cb | and of the branching fraction B( B0 → D * + ℓ − νℓ ) performed on a sample of B0 → D * + ℓ − νℓ decays selected from events with a high momentum charged lepton ℓ − and a D * +1 .The D * + is reconstructed from its decay to a charged pion and a D 0 .This pion is produced with small momentum, and is commonly referred to as the slow pion (π slow ).
This paper is organized as follows.The next section is dedicated to a brief description of the BABAR detector and of the data sets employed.Section 3 describes the event selection and the composition of the sample used for the measurement.The fit method and its results are described in Section 4; the evaluation of the systematic error is discussed in Section 4.4.Finally, Section 5 presents the conclusions and the comparison with other measurements of |V cb |.
The BABAR detector is described in detail elsewhere [16].The momenta of charged particles are measured using a combination of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) in a 1.5 T solenoidal magnetic field.A detector of internally-reflected Cherenkov radiation (DIRC) is used for charged-particle identification.Kaons are identified with a neural network based on the likelihood ratios calculated from dE/dx measurements in the SVT and DCH, and from the information from the DIRC.A finely segmented CsI(Tl) electromagnetic calorimeter (EMC) is used to detect photons and neutral hadrons, and to identify electrons.Electron candidates are required to have a ratio of EMC energy to track momentum, EMC cluster shape, DCH dE/dx, and DIRC Cherenkov angle all consistent with the electron hypothesis.The instrumented flux return (IFR) contains resistive plate chambers for muon and long-lived neutral hadron identification.Muon candidates are required to have IFR hits located along the extrapolated DCH track, and energy deposition in the EMC and penetration length in the IFR consistent with a minimum ionizing particle.

Event Selection
Events are selected that contain a fully-reconstructed D * + and an oppositely charged lepton.Those D * + ℓ − combinations not coming from the signal decay are suppressed by means of kinematic cuts.
A candidate lepton is searched for among all charged tracks with momentum greater than 1.2 GeV/c in the Υ (4S) rest frame.Electron candidates are selected with an efficiency of about 90% and a hadron misidentification probability of less than 0.2%.Muon candidates are selected with an efficiency of ≃ 60% and hadron misidentification probability is about 2%.
The D * + candidate is selected in the decay mode D * + → D 0 π + and the D 0 meson is reconstructed in the three modes K − π + , K − π + π + π − , and K − π + π 0 .The π 0 is reconstructed from two photons, each with energy larger than 30 MeV, and must have a total energy larger than 200 MeV and an invariant mass between 119.2 and 150.0 MeV/c 2 .The invariant mass of the two photons is constrained to the π 0 mass and the pair is kept as a π 0 candidate if the χ 2 probability of the fit exceeds 1%.Charged kaon candidates satisfy loose identification criteria for the K − π + mode and tighter criteria for the K − π + π + π − and K − π + π 0 modes.For the K − π + π 0 mode, the D 0 candidate is retained if the square of the decay amplitude in the Dalitz plot for the three-body candidate, based on measured amplitudes and phases [17], is larger than 10% of its maximum value.D 0 candidates are accepted if they have an invariant mass within 17 MeV/c 2 of the D 0 mass for the K − π + and K − π + π + π − modes, and within 34 MeV/c 2 for the K − π + π 0 mode.The invariant mass of the decay products is then constrained to the D 0 mass [18] and the tracks are constrained to a common vertex by a simultaneous fit.The D 0 candidate is retained if the χ 2 probability of the fit is more than 0.1%.The low-momentum pion candidates for the D * + → D 0 π + slow decay are selected from among tracks with total momentum in the laboratory frame less than 450 MeV/c, and momentum transverse to the beam line in the laboratory frame greater than 50 MeV/c.The momentum of the D * + candidate in the Υ (4S) rest frame is required to be between 0.5 and 2.5 GeV/c.These requirements retain essentially all signal events and reject higher momentum D * + from continuum events.Candidates from D * + are preselected with the cut on the mass difference ∆M = M D 0 π slow − M D 0 < 165 MeV/c 2 .The ∆M distribution has a kinematic threshold at the mass of the π + , and a peak at 145.5 MeV/c 2 with a resolution of about 1 MeV/c 2 .The combinatoric background below the signal is evaluated using events in the ∆M sideband.A harder ∆M cut is applied later.
The χ 2 probability of the fit to a common vertex of the lepton, the π slow , and the D 0 candidate, constrained to the beam spot, must be greater than 1%.
Candidates for D * + ℓ − are selected if | cos ∆θ * thrust | < 0.85, where ∆θ * thrust is the angle between the thrust axis of the D * + ℓ − candidate and the thrust axis of the remaining charged and neutral particles in the event.The distribution of | cos ∆θ * thrust | is peaked at 1 for jet-like continuum events, and is flat for BB events, where the two B mesons decay independently and practically at rest with respect to each other.This cut retains about 85% of signal and rejects about 47% of continuum events.
Finally, the angle between the lepton and the D * + in the Υ (4S) rest frame must satisfy cos θ D * ℓ < 0. This cut reduces significantly the events due to random combinations of a lepton and a D * + produced by two different B mesons, while keeping most of the signal events.

Sample Composition
Several processes contribute to the selected sample: 1. Signal decays.
2. Combinatoric background from any random combination of tracks mimicking a true D 0 or from a correctly identified D 0 , combined with a low-momentum charged track not originating from D * + decay.This category includes events from B B and from continuum production.Because of its combinatorial origin, the ∆M distribution does not exhibit a peak.All the other backgrounds listed below peak in ∆M .
3. Fake-lepton background, in which a charged hadron is misidentified as a lepton and is combined with a correctly reconstructed D * + .Due to the very good performance of the electron identification, the rate of fake electrons is negligible.
4. Continuum background, where a true D * + and a lepton (true or fake) are produced from e + e − → cc events.
5. "Uncorrelated" background, in which the lepton and the D * + originate from two different B mesons.Most of these combinations are removed by the requirement that the lepton and D * + have opposite charge.Opposite charge combinations arise mainly from events with B 0 B 0 oscillations where both B mesons decay into the same flavor eigenstate, or events in which the D * + comes from the hadronization of the virtual W from B decay, and thus has a charge opposite to the one expected from its parent B.
6. B → D * + ℓ − νℓ X background, where one (or more) hadrons are produced in addition to the D * + , for either neutral or charged B. A considerable fraction of this background is due to the intermediate production of orbitally excited D states, which then decay to a D * + π final state; this is referred to as D * * background in the following.
7. "Correlated" background, in which the B decays to a D * + and a heavy particle (either a charm meson or a τ ), which then decays to a lepton.These events have the same charge correlation as signal events, but are suppressed by their low branching fraction and by the kinematic cuts.Such events correspond to a few percent of the peaking sample.Their contribution is fixed to the value computed in the simulation.(right) candidates in on-resonance data.Points correspond to data, and the curve is the result of unbinned maximum likelihood fits to this signal sample, the fake-lepton control sample and the rescaled off-resonance events.The shaded distributions correspond to the continuum, fake-lepton, and combinatorial backgrounds described in the text.The correlated-lepton, the B → D * + ℓ − νℓ X and the uncorrelated backgrounds are included in the signal distribution.
The fraction of combinatoric events below the D * + peak is determined from a fit to the ∆M distribution in the region 0.138 < ∆M < 0.165 GeV/c 2 (see Figure 1).In the fit, D * + events are described by the sum of two Gaussian functions while the shape of combinatoric events is described by the function 1 − exp − ∆m−mπ c 1 ∆m mπ c 2 , where the overall normalization, c 1 , and c 2 are determined by the fit.The width of the ∆M distribution for signal events is dominated by the experimental resolution.In one third of events the π slow track is reconstructed in the SVT and in the DCH; in the remainder, the track is reconstructed only in the SVT.As the resolution in the first sample is better, the two samples are fitted separately.
The ∆M distribution of a fake-lepton control sample is fitted simultaneously to determine the fake-lepton background.The control sample is obtained by selecting events with the kinematic cuts described above, except for the requirement that the candidate lepton fails a loose lepton selection criterion.Lepton identification efficiencies and hadron misidentification probabilities, estimated from data control samples of pure leptons and hadrons, are used to scale the number of peaking events in this hadron-substituted sample to the expected amount of fake-lepton background in the signal sample.
The continuum contribution is fitted from the off-resonance event sample and scaled according to the ratio between the on-resonance and off-resonance luminosities.
For each D 0 final state, the signal, fake-lepton and continuum samples are fitted simultaneously.The mean values, widths and the relative normalization of the two Gaussian functions describing the signal are common to the three data sets, while the parameters describing the shape of the combinatoric background are fitted independently to each data sample.Cuts on ∆M are then applied to reduce the amount of combinatorics in the subsequent stages of the analysis: 0.143 < ∆M < 0.148 GeV/c 2 if the π slow is reconstructed by the SVT alone, 0.144 < ∆M < 0.147 GeV/c 2 otherwise.
The fraction of uncorrelated and D * * events is determined by exploiting the decay kinematics.The angle between the B 0 and the pseudo-particle obtained by adding the D * + and the ℓ − fourmomenta is computed according to the following expression: where M miss is the invariant mass of all the other particles produced in association with the lepton and the D * + in the decay of the B 0 .We compute cos θ B 0 ,D * ℓ by fixing M miss to zero, under the hypothesis that this is a true B0 → D * + ℓ − νℓ decay, for which only a neutrino is missing.Neglecting the broadening due to experimental resolution, the obvious constraint −1 < cos θ B 0 ,D * ℓ < 1 applies to signal events.Due to the production of additional particles, events from D * * background have positive M 2 miss values and therefore produce large tails below -1 in cos θ B 0 ,D * ℓ distribution, while uncorrelated background events preferentially populate the region with cos θ B 0 ,D * ℓ > 1.This allows for the determination of the fraction of D * * , uncorrelated, and signal events in the sample, by means of a fit to the cos θ B 0 ,D * ℓ distribution in data (see for instance Figure 2).The shapes of all components in the fit are determined from Monte Carlo simulation, while the number of events for each of the other background categories is determined from the fit to the ∆M distributions.The same procedure is applied to 18 independent samples, identified according to lepton flavors (×2), period of data acquisition (×3), D 0 decay mode (×3).The data sets are further split into bins of w (defined in Section 4.1) in order to reduce the systematic error corresponding to variations in background composition and normalization with w.The resulting large statistical uncertainties are reduced by requiring that the fractions of uncorrelated and D * * background with respect to the signal be the same for electrons and muons, and for all D 0 decay modes.The validity of these hypotheses is verified by Monte Carlo simulation.The distributions in cos θ B 0 ,D * ℓ for the B0 → D * + e − νe , D 0 → K − π + and B0 → D * + µ − νµ , D 0 → K − π + π 0 samples (integrated over w ) are shown in Figure 2.
Only events with | cos θ B 0 ,D * ℓ | < 1.2 are further analyzed.This selects a final sample of 55700 signal events for the measurement.

Parametrization of the Decay Width
Assuming the connection between the various form factors provided by HQET, the expected number of signal events can be expressed as a function of w by the relation The factor of 4 accounts for the fact that two B mesons are produced in each event, and that both electrons and muons are used; N Υ is the number of Υ (4S) produced, f 00 is the fraction Γ(Υ (4S) → is the branching fraction for the decay D * + → D 0 π + , B(D 0 → K − nπ) is the branching fraction for the decay of the D 0 to the final state considered, τ B 0 is the B 0 lifetime and r is the ratio of meson masses, r = M D * + /M B 0 .The values employed for these parameters, as determined by independent measurements, are reported in Table 1.
Table 1: Values of inputs to the decay width calculation, as obtained from Ref. [18].
The w-dependent reconstruction efficiency ǫ(w) is determined by means of the detailed detector simulation.The w dependence of the form factor F(w) is unknown.Since only a small range of w is allowed by phase space, a Taylor series expansion limited to second order has typically been used: Apart from F(1), the theory does not predict the values of the higher order coefficients, which must be determined experimentally.The first measurements of |V cb | were performed assuming a linear expansion, i.e., setting c = 0 ( [7,9,10,13]).However, the requirement of analyticity and positivity of the QCD functions describing the local currents leads to the prediction that c should be positive and it should be related to the heavy meson radius ρ 2 F (see [5]) by the relation: Results employing this analyticity bound were obtained by the ALEPH, DELPHI and OPAL collaborations (see [11,12,14]).An improved parametrization based on the above consideration is proposed in [6], accounting for higher order terms and so reducing to ±2% the relative uncertainty on |V cb | due to the form factor parametrization.A new function A 1 (w) is introduced in that calculation, connected to F(w) by the following relation: where R 1 (w) and R 2 (w) are ratios of axial and vector form factors also given in Ref. [6].The following parametrization, depending only on a single unknown parameter ρ 2 A 1 , is proposed for A 1 (w): It should be noted that, in the limit w → 1, A 1 (w) → F(w), so that A 1 (1) = F(1).Monte Carlo events employed for this analysis were produced with a linear parametrization for F(w), while experimental data are fitted using expressions of Ref. [6].
Experimentally the variable w can be expressed as: where q 2 ≡ (p B 0 − p D * + ) 2 .The magnitude of the B 0 momentum is known; its direction is obtained from Equation 1 with an azimuthal ambiguity about the direction of the D * + ℓ − pair.The two extreme solutions corresponding to the minimal and maximal angles between the B 0 and the D * + are used to define the quantity The simulation shows that w is a good estimator of w, with a resolution σ( w − w) ∼ 0.02, corresponding to about 4% of the full physical range.

Fit Method
Data and Monte Carlo simulated events are collected in ten equal-size w bins.A few events, which, due to resolution, have w > w max = 1.503 are discarded.A least-squares fit is then performed comparing the number of events observed in each bin to the sum of signal and background events.The number of background events in each bin is determined as explained above and is fixed in the fit.The signal contribution is obtained at each step in the minimization procedure by properly weighting each generated event surviving the selection.The χ 2 is: where the index i runs over the ten w bins; N i data(MC) is the number of data (signal Monte Carlo) events found in the i th bin; N i bck , σ i bck are the numbers of background events and their errors.The weight for each Monte Carlo event is computed as the product of three terms where • W L is an overall fixed scale factor, that accounts for the relative luminosity of data and signal Monte Carlo events, for the difference in the actual branching fractions relevant to this analysis3 , and any possible overall efficiency scale factor.
• W ǫ i accounts for the efficiency correction, computed event per event as the ratio of the efficiency in the real data and in the Monte Carlo.It accounts for differences in tracking, particle identification, π 0 reconstruction, and depends on the three-momenta of the lepton, the soft pion and of all the particles from the D 0 decay.It remains unchanged at all steps of the minimization.(If the Monte Carlo simulation were perfectly tuned, this factor would be identically one).It should be noted that only the factor due to the slow pion tracking efficiency introduces a net dependence of this term on w.
is the term accounting for the theoretical function describing the decay.It depends on the parameters to be determined (A 1 (1)|V cb | and ρ 2 A 1 ) and varies for each step of the minimization.
Here the function f theo corresponds to the expressions in Equations 3 and 5, while the function used to generate Monte Carlo events, f M C , contains a simpler, linear parametrization of the form factor dependence on w.
The value of the branching ratio is then computed by integrating the differential expression in Equation 3.
The fit has been performed as a blind analysis, i.e., the values of A 1 (1)|V cb | and ρ 2 A 1 were hidden until the study of the systematic errors was completed and all consistency checks were performed.

Fit Results
The fit is performed separately for each data set, providing 18 statistically independent determinations of A 1 (1)|V cb | and of ρ 2 A 1 .The results are reported in Table 2.In the same Table, the branching fraction, as obtained by integrating the differential decay width, is also shown.The agreement between the data and the fit function is good.The distribution of the χ 2 probabilities is uniform.The measured w distributions are shown in Figure 3, where the three distributions for the D 0 decay modes are summed.
It should be noted that, because the statistical errors are small, the measurements are dominated by the systematic errors, which introduce considerable correlations among the samples.For this reason, first the systematic error is discussed, and then the average is presented.

Systematic Errors and Consistency Checks
The individual sources of systematic errors are reported in Table 4 and are described in detail below.They can be grouped as • Factors affecting the overall normalization (see Equation 3): Each of these factors is varied by its uncertainty, and the resulting variation in the parameters is added to the systematic error.The uncertainties on N Υ , f 00 , τ B 0 and B(D * + ) are completely correlated for all the samples.The uncertainties on the D 0 branching fractions depend on the decay mode and are only partially correlated [18].
• Experimental uncertainties in the efficiency: These include track reconstruction, lepton and kaon identification, π 0 selection, vertex reconstruction and the cut on the Dalitz plot for the K − π + π 0 mode.For high-momentum tracks, the tracking efficiency is determined comparing the independent information from the SVT and DCH.This approach results in ±0.8% systematic uncertainty in the reconstruction efficiency for each track; the errors from each track are added linearly (i.e., the error on the branching fraction is N times the single track error, where N is the total number of tracks employed, including the lepton).Most of the π slow do not reach the DCH, and are reconstructed in the SVT alone.The efficiency for these low-momentum tracks is computed from the angular distribution of the π slow in the D * + rest frame.To evaluate this, a large set of D * + → D 0 π + slow , D 0 → K − π + decays is selected from generic hadronic events.For fixed values of the D * + momentum, the observed angular distribution is compared to the one predicted for the decay of a vector meson to two pseudoscalar mesons, and the deviation from the expected shape is attributed to the inefficiency of track reconstruction.This inefficiency is parametrized as a function of the particle momentum in the laboratory frame.The values of the parameters describing the efficiency are computed in different bins of the polar angle of the track, both in the data and in the simulation.The systematic error is computed by varying each parameter of the efficiency function by its uncertainty, including correlations.
The efficiencies and systematic errors for lepton and kaon identification and π 0 reconstruction lated background by a factor of about 2 in the full cos θ B 0 ,D * ℓ range, and by 30% in the signal region.No appreciable change in the fit result is observed.As an alternative, the fraction of uncorrelated events is determined by counting the events in the rejected region cos θ D * ℓ > 0 and propagating this number to the complementary region using the Monte Carlo.The fraction obtained is then fixed in the cos θ B 0 ,D * ℓ fit, which is used to determine the fraction of D * * and signal events only.Once again, no appreciable change in the result is observed.Finally, the hypothesis that the ratio of the amount of D * * and uncorrelated background over the signal be the same for all D * + decay modes is checked with the simulation.While the test is satisfactory for the uncorrelated background, a slight discrepancy is observed for D * * events.For K − π + π + π − decay the fitted fraction is higher by 30% than for the other modes.The effect is the same for electron and muon events and does not depend on w.Therefore the w fit is repeated by increasing the D * * amount in the K − π + π + π − sample and decreasing it correspondingly in the K − π + and K − π + π 0 .The difference in the result is propagated as systematic error.
• Theoretical uncertainties: A considerable fraction of the allowed phase-space is removed by the lower limit on the lepton momentum.This induces a model-dependence in the computation of the branching fraction.The lepton spectrum is determined partly by the fitted shape of the Isgur-Wise function, and partly by the angular profile of the decay, which is governed by the parameters R 1 and R 2 defined above.Even with perfect acceptance, the uncertainty on R 1 and R 2 directly affects ρ 2 A 1 and A 1 (1)|V cb |, because these parameters enter in the fit function.The values of these parameters determined by the CLEO collaboration [21] are used in this analysis.The systematic error is taken as the observed variation of A 1 (1)|V cb |, ρ 2 A 1 and of B( B0 → D * + ℓ − νℓ ), when R 1 and R 2 are floated by their uncertainties, accounting for their correlation.
• Fit method: The fit procedure is validated using an independent sample of Monte Carlo simulated events, of the same size as the data, and following the procedure described above, using however a linear extrapolation for the Isgur-Wise function.The fitted values of ρ 2 A 1 and A 1 (1)|V cb | are consistent within their statistical errors with the input values.The statistical error on this test is added to the systematic error.In addition, to test that the fit method does not bias the result, a set of 500 toy experiments are performed, in which random events are generated with realistic efficiency and resolution.Each toy experiment has the same sample size for data and Monte Carlo as in the actual measurement.Within errors, no difference is observed between the average of the fitted parameters and their generated values.Also, it is observed that the pull width is consistent with one.

Results
The results from all subsamples are combined by using the COMBOS package [22], taking into account the correlations between the samples.COMBOS was originally developed by the LEP-B-Oscillation working group to combine the results of ∆M d from different experiments, and was then adapted by the LEP-V cb working group to compute the world average for |V cb |, ρ 2 A 1 and B( B0 → D * + ℓ − νℓ ).Besides computing averages, errors, and confidence levels, COMBOS also provides a breakdown of the single error sources, which is reported in Table 4.
All the parameters affecting the normalization are treated as fully correlated.The errors on the branching fractions of the D 0 to each final state are fully correlated for each state, but, at present, the correlations among different decay modes are neglected.The systematic errors on particle identification are considered as fully correlated among data taking periods.The error on tracking is determined by systematic effects common to all years and therefore is considered as fully correlated among all the samples.In contrast the error on π slow tracking efficiency is still dominated by the statistical uncertainty from the control samples; therefore it is completely correlated for samples collected in the same year of data acquisition, but uncorrelated for different years.Model errors are common to all data sets.The mean values are obtained by first averaging results for different running periods, and then computing the average result for the three years.The confidence level of the A 1 (1)|V cb | and ρ 2 A 1 average values is 30%.The results are reported in Table 5.For the sum of all 18 data subsamples the following result is obtained: ρ 2 A 1 = 1.23 ± 0.02 ± 0.28, (10) B( B0 → D * + ℓ − νℓ ) = (4.68 ± 0.03 ± 0.29) %, (11) where the first error is statistical and the second is systematic.The statistical correlation between A 1 (1)|V cb |and ρ 2 A 1 is 92%.

Conclusions
A sample of 55,700 signal events from the decay process B0 → D * + ℓ − νℓ is selected from a set of 86 × 10 6 Υ (4S) → BB decays collected by the BABAR detector.The product of |V cb | and the form factor at zero recoil, F(1) = A 1 (1), is measured, as well as the derivative of the form factor, ρ 2 A 1 , again at zero recoil.The integrated branching fraction is also computed.These results are in agreement with those obtained by Belle [8] and with the LEP average [18], but differ by more than two standard deviations from the CLEO result [9].Using the value A(1) = 0.913 +0.030 −0.035 reported in [4], based on a lattice QCD calculation, the value: Both the experimental and theoretical errors are independent.By interpreting the theoretical errors as ±68% confidence interval, the two values differ by 2.0 standard deviations.

Figure 1 :
Figure1: ∆M distribution for events passing all selection criteria: D * + e − (left) and D * + µ − (right) candidates in on-resonance data.Points correspond to data, and the curve is the result of unbinned maximum likelihood fits to this signal sample, the fake-lepton control sample and the rescaled off-resonance events.The shaded distributions correspond to the continuum, fake-lepton, and combinatorial backgrounds described in the text.The correlated-lepton, the B → D * + ℓ − νℓ X and the uncorrelated backgrounds are included in the signal distribution.

Figure 2 :
Figure 2: Distributions in cos θ B 0 ,D * ℓ for the B0 → D * + e − νe , D 0 → K − π + (left) and B0 → D * + µ − νµ , D 0 → K − π + π 0 samples (right).Points with error bars show the data.Shaded areas represent the Monte Carlo distributions.The fit determines the normalization of the Monte Carlo distributions.The tail to the left in the signal electron sample is due to bremsstrahlung.Below the histograms are shown the fractional deviations of the data from the fit results.

Figure 3 :
Figure 3: The experimental distribution of w compared to fit results separately for electrons (left) and muons (right).The plots from all the different D 0 decay modes and years of data taking are summed together.The points are the data and the histograms are the results of the fit to the sum of the signal and the different sources of background.The fit residuals are also shown below.

Table 5 :
Fit results for the three running periods.