Determination of the form factors for the decay B(0)-> D*(-)l(+)nu(l) and of the CKM matrix element vertical bar Vcb vertical bar RID C-5223-2009 RID C-5719-2008 RID D-1055-2009

B. Aubert, M. Bona, D. Boutigny, Y. Karyotakis, J. P. Lees, V. Poireau, X. Prudent, V. Tisserand, A. Zghiche, J. Garra Tico, E. Grauges, L. Lopez, A. Palano, G. Eigen, I. Ofte, B. Stugu, L. Sun, G. S. Abrams, M. Battaglia, D. N. Brown, J. Button-Shafer, R. N. Cahn, Y. Groysman, R. G. Jacobsen, J. A. Kadyk, L. T. Kerth, Yu. G. Kolomensky, G. Kukartsev, D. Lopes Pegna, G. Lynch, L. M. Mir, T. J. Orimoto, M. Pripstein, N. A. Roe, M. T. Ronan,* K. Tackmann, W. A. Wenzel, P. del Amo Sanchez, C. M. Hawkes, A. T. Watson, T. Held, H. Koch, B. Lewandowski, M. Pelizaeus, T. Schroeder, M. Steinke, W. N. Cottingham, D. Walker, D. J. Asgeirsson, T. Cuhadar-Donszelmann, B. G. Fulsom, C. Hearty, N. S. Knecht, T. S. Mattison, J. A. McKenna, A. Khan, M. Saleem, L. Teodorescu, V. E. Blinov, A. D. Bukin, V. P. Druzhinin, V. B. Golubev, A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, K. Yu Todyshev, M. Bondioli, S. Curry, I. Eschrich, D. Kirkby, A. J. Lankford, P. Lund, M. Mandelkern, E. C. Martin, D. P. Stoker, S. Abachi, C. Buchanan, S. D. Foulkes, J. W. Gary, F. Liu, O. Long, B. C. Shen, L. Zhang, H. P. Paar, S. Rahatlou, V. Sharma, J. W. Berryhill, C. Campagnari, A. Cunha, B. Dahmes, T. M. Hong, D. Kovalskyi, J. D. Richman, T. W. Beck, A. M. Eisner, C. J. Flacco, C. A. Heusch, J. Kroseberg, W. S. Lockman, T. Schalk, B. A. Schumm, A. Seiden, D. C. Williams, M. G. Wilson, L. O. Winstrom, E. Chen, C. H. Cheng, A. Dvoretskii, F. Fang, D. G. Hitlin, I. Narsky, T. Piatenko, F. C. Porter, G. Mancinelli, B. T. Meadows, K. Mishra, M. D. Sokoloff, F. Blanc, P. C. Bloom, S. Chen, W. T. Ford, J. F. Hirschauer, A. Kreisel, M. Nagel, U. Nauenberg, A. Olivas, J. G. Smith, K. A. Ulmer, S. R. Wagner, J. Zhang, A. M. Gabareen, A. Soffer, W. H. Toki, R. J. Wilson, F. Winklmeier, Q. Zeng, D. D. Altenburg, E. Feltresi, A. Hauke, H. Jasper, J. Merkel, A. Petzold, B. Spaan, K. Wacker, T. Brandt, V. Klose, H. M. Lacker, W. F. Mader, R. Nogowski, J. Schubert, K. R. Schubert, R. Schwierz, J. E. Sundermann, A. Volk, D. Bernard, G. R. Bonneaud, E. Latour, V. Lombardo, Ch. Thiebaux, M. Verderi, P. J. Clark, W. Gradl, F. Muheim, S. Playfer, A. I. Robertson, Y. Xie, M. Andreotti, D. Bettoni, C. Bozzi, R. Calabrese, A. Cecchi, G. Cibinetto, P. Franchini, E. Luppi, M. Negrini, A. Petrella, L. Piemontese, E. Prencipe, V. Santoro, F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, S. Pacetti, P. Patteri, I. M. Peruzzi, M. Piccolo, M. Rama, A. Zallo, A. Buzzo, R. Contri, M. Lo Vetere, M. M. Macri, M. R. Monge, S. Passaggio, C. Patrignani, E. Robutti, A. Santroni, S. Tosi, K. S. Chaisanguanthum, M. Morii, J. Wu, R. S. Dubitzky, J. Marks, S. Schenk, U. Uwer, D. J. Bard, P. D. Dauncey, R. L. Flack, J. A. Nash, M. B. Nikolich, W. Panduro Vazquez, P. K. Behera, X. Chai, M. J. Charles, U. Mallik, N. T. Meyer, V. Ziegler, J. Cochran, H. B. Crawley, L. Dong, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, A. V. Gritsan, Z. J. Guo, C. K. Lae, A. G. Denig, M. Fritsch, G. Schott, N. Arnaud, J. Béquilleux, M. Davier, G. Grosdidier, A. Höcker, V. Lepeltier, F. Le Diberder, A. M. Lutz, S. Pruvot, S. Rodier, P. Roudeau, M. H. Schune, J. Serrano, V. Sordini, A. Stocchi, W. F. Wang, G. Wormser, D. J. Lange, D. M. Wright, C. A. Chavez, I. J. Forster, J. R. Fry, E. Gabathuler, R. Gamet, D. E. Hutchcroft, D. J. Payne, K. C. Schofield, C. Touramanis, A. J. Bevan, F. Di Lodovico, K. A. George, W. Menges, R. Sacco, G. Cowan, H. U. Flaecher, D. A. Hopkins, P. S. Jackson, T. R. McMahon, F. Salvatore, A. C. Wren, D. N. Brown, C. L. Davis, J. Allison, N. R. Barlow, R. J. Barlow, Y. M. Chia, C. L. Edgar, G. D. Lafferty, T. J. West, J. I. Yi, J. Anderson, C. Chen, A. Jawahery, D. A. Roberts, G. Simi, J. M. Tuggle, G. Blaylock, C. Dallapiccola, S. S. Hertzbach, X. Li, T. B. Moore, E. Salvati, S. Saremi, R. Cowan, P. H. Fisher, G. Sciolla, S. J. Sekula, M. Spitznagel, F. Taylor, R. K. Yamamoto, S. E. Mclachlin, P. M. Patel, S. H. Robertson, A. Lazzaro, F. Palombo, J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Kroeger, D. A. Sanders, D. J. Summers, H. W. Zhao, S. Brunet, D. Côté, M. Simard, P. Taras, F. B. Viaud, H. Nicholson, G. De Nardo, F. Fabozzi, L. Lista, D. Monorchio, C. Sciacca, M. A. Baak, G. Raven, H. L. Snoek, C. P. Jessop, J. M. LoSecco, G. Benelli, L. A. Corwin, K. K. Gan, K. Honscheid, D. Hufnagel, H. Kagan, R. Kass, J. P. Morris, A. M. Rahimi, J. J. Regensburger, R. Ter-Antonyan, Q. K. Wong, N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, M. Lu, R. Rahmat, N. B. Sinev, D. Strom, J. Strube, E. Torrence, N. Gagliardi, A. Gaz, M. Margoni, M. Morandin, A. Pompili, M. Posocco, M. Rotondo, F. Simonetto, R. Stroili, C. Voci, E. Ben-Haim, H. Briand, J. Chauveau, P. David, L. Del Buono, Ch. de la Vaissière, O. Hamon, B. L. Hartfiel, Ph. Leruste, J. Malclès, J. Ocariz, PHYSICAL REVIEW D 77, 032002 (2008)


I. INTRODUCTION
The study of the semileptonic decay B 0 !D ÿ ' ' [1] is interesting in many respects.In the standard model, the rate of this weak decay is proportional to the square of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V cb , which is a measure of the weak coupling of the b to the c quark.This decay is also influenced by strong interactions.Their effect can be parametrized by two axial form factors, A 1 and A 2 , and one vector form factor, V, each of which depends on the momentum transfer squared q 2 of the B meson to the D meson.The form of this dependence is not known a priori.In the framework of heavy-quark effective field theory (HQET) [2,3], these three form factors are related to each other through heavy-quark symmetry (HQS), but HQET leaves three free parameters, which must be determined by experiment.
The extraction of jV cb j relies on the measurement of differential decay rates.HQS predicts the normalization of the decay rate at the maximum q 2 , and jV cb j is determined from an extrapolation of the form factors to this value.The precise determination requires corrections to the HQS prediction for the normalization, as well as a measurement of the variation of the form factors near the maximum q 2 , where the decay rate goes to zero as the phase space vanishes.
Several experiments have measured jV cb j based on studies of the differential decay width for B 0 !D ÿ ' ' decays [4 -10].These analyses of the one-dimensional differential decay rate resulted in the measurement of only one of the form-factor parameters, and they relied on a measurement by the CLEO Collaboration [11] for the other two.The uncertainty in these two parameters introduces the largest systematic uncertainty in all previous measurements of jV cb j using this method.
Furthermore, for measurements of jV ub j, based on both inclusive and exclusive B ! X u ' decays, improved knowledge of all form factors is important to correctly describe the dominant B ! X c ' background.
In this paper, we present a simultaneous measurement of jV cb j, of the branching fraction for B 0 !D ÿ ' ' , and of the three form-factor parameters, based on measurements of three one-dimensional decay distributions.Thus we extend the earlier BABAR measurement [10] where F1jV cb j and one of the form-factor parameters are measured, fully accounting for correlations between these onedimensional distributions.
We combine the results of this analysis with another BABAR measurement of the form factors [12], which employs a fit to the full four-dimensional decay distribution on a partial sample of the data.
The D ÿ candidates are reconstructed from the D ÿ !D 0 ÿ decays, and the D 0 mesons are reconstructed in three different decay modes, K ÿ , K ÿ ÿ , and K ÿ 0 .Electrons or muons are paired with the D ÿ to form signal candidates.The large data sample permits a precise determination of the background contributions, largely based on data, and thus results in smaller experimental uncertainties.
This leads to a further reduction of the form-factor errors.

II. FORMALISM A. Kinematic variables
The lowest-order quark-level diagram for the decay B 0 !D ÿ ' ' is shown in Fig. 1.This decay is completely characterized by four variables, namely, three angles and the Lorentz-invariant variable w, which is linearly related to q 2 and defined as where m B and m D are the masses of the B and the D mesons (2.010 and 5.2794 GeV, respectively [13]), and P B and P D are their four-momenta.In the B rest frame the expression for w reduces to the Lorentz boost D E D =m D .The ranges of w and q 2 are restricted by the kinematics of the decay, with q 2 0 corresponding to and w min 1 corresponding to q 2 max m B ÿ m D 2 10:69 GeV 2 : (3) The three angular variables, shown in Fig. 2, are as follows: (i) ' , the angle between the direction of the lepton in the virtual W rest frame and the direction of the W in the B rest frame; (ii) V , the angle between the direction of the D in the D rest frame and the direction of the D in the B rest frame; (iii) , the angle between the plane formed by the D and the plane formed by the W decay.

B. Four-dimensional decay distribution
The Lorentz structure of the B 0 !D ÿ ' ' decay amplitude can be expressed in terms of three helicity amplitudes (H , H ÿ , and H 0 ), which correspond to the three polarization states of the D , two transverse and one longitudinal.For low-mass leptons, these amplitudes are expressed in terms of the three functions h A 1 w, R 1 w, and R 2 w [2,3]: where The functions R 1 w and R 2 w are defined in terms of the axial and vector form factors as By convention, the function h A 1 w is defined as For w ! 1, the axial form factor A 1 w dominates, and in the limit of infinite bÿ and cÿ quark masses, a single form factor describes the decay, the so-called Isgur-Wise function [14].
The fully differential decay rate in terms of the three helicity amplitudes is By integrating this decay rate over all but one of the four variables, w, cos ' , cos V , or , we obtain the four onedimensional decay distributions from which we extract the form factors.The differential decay rate as a function of w is where FIG. 2. Definition of the angles ' , V , and for the B 0 !D ÿ ' ' decay, mediated by a vector boson W; s refers to the low-momentum pion from the decay D ÿ !D 0 ÿ s .
and Gw is a known phase space factor, It is important to note that h A 1 1 F 1 corresponds to the Isgur-Wise function [14] at w 1.In the infinite quark-mass limit, the HQS normalization gives F 1 1.Corrections to this HQS prediction have been calculated in the framework of lattice QCD.A recent calculation, performed in a quenched approximation, predicts (including a QED correction of 0.7%) F 1 0:919 0:030 ÿ0:035 [15].This value is compatible with estimates based on nonlattice methods [16,17].

C. Form-factor parametrization
Since HQET does not predict the functional form of the form factors, a parametrization is needed for their extraction from the data.Perfect heavy-quark symmetry implies that R 1 w R 2 w 1, i.e., the form factors A 2 and V are identical for all values of w and differ from A 1 only by a simple kinematic factor.Corrections to this approximation have been calculated in powers of QCD =m b and the strong coupling constant s .Various parametrizations in powers of w ÿ 1 have been proposed.Among the different predictions relating the coefficients of the higher order terms to the linear term, we adopt the following expressions derived by Caprini, Lellouch, and Neubert [18], The three parameters, 2 , R 1 1, and R 2 1, cannot be calculated; they must be extracted from data.

III. DATA SAMPLE, RECONSTRUCTION, AND SIMULATION
The data used in this analysis were recorded by the BABAR detector, operating at the PEP-II asymmetricenergy e e ÿ collider.The data sample corresponds to a luminosity of 79 fb ÿ1 recorded on the 4S resonance (on-resonance sample), and 9:6 fb ÿ1 recorded at a centerof-mass energy 40 MeV lower (off-resonance sample).
The BABAR detector and event reconstruction have been described in detail elsewhere [19,20].The momenta of charged particles are measured by a tracking system consisting of a five-layer silicon vertex tracker (SVT) and a 40layer drift chamber (DCH).Charged particles of different masses are distinguished by their energy loss in the tracking devices and by a ring-imaging Cerenkov detector (DIRC).Electromagnetic showers from electrons and photons are measured in a CsI(Tl) calorimeter (EMC).These detector components are embedded in a 1.5-T magnetic field of the solenoid.Electron candidates are selected on the basis of the ratio of the energy detected in the calorimeter to the track momentum, the calorimeter shower shape, the energy loss in the drift chamber, and the angle of the photons reconstructed in the DIRC.Muons are identified in a set of resistive plate chambers inserted in the steel flux-return of the magnet (IFR).Information from the IFR is combined with the track momentum measurement and energy deposition in the EMC to improve the separation of muons from charged hadrons.
The electron and muon identification efficiencies and the probabilities to misidentify a pion, kaon, or proton as an electron or muon have been measured as a function of the laboratory momentum and the angles with clean samples of tracks selected from data [21].Within the acceptance of the calorimeter, defined by the polar angle in the laboratory frame, ÿ0:72 < cos lab < 0:92, and above 1.0 GeV, the average electron efficiency is 91%, largely independent of the electron momentum.The average hadron misidentification rate is less than 0.2%.The muon detection extends to polar angles of ÿ0:91 < cos lab < 0:95.For a hadron misidentification rate of typically 2.0%, the average efficiency is close to 65%.
The criteria for distinguishing charged kaons from charged pions are chosen to maximize the efficiency while controlling the background, and thus differ for the decay modes under study.Consequently, the efficiency varies from 87% to 97%.The uncertainties are typically 1%.
We determine the tracking efficiency for highmomentum tracks by comparing the independent information from SVT and DCH.We compute the efficiency for low-momentum tracks reconstructed in the SVT alone from the angular distribution of the ''slow'' pion, ÿ s , in the D ÿ rest frame.We use a large sample of D ÿ !D 0 ÿ s , D 0 !K ÿ decays selected from hadronic B decays.For fixed values of the D ÿ momentum, we compare the observed angular distribution to the one expected for the decay of a vector meson to two pseudoscalar mesons.This study is performed in several bins of the polar angle.We define the relative efficiency as the ratio of the observed to the expected distribution and parametrize its dependence on the laboratory momentum of the ÿ s .Below 100 MeV, this efficiency drops very steeply, and reaches zero at about 60 MeV.
Neutral pions are reconstructed from pairs of photon candidates of more than 30 MeV detected in the EMC and assumed to originate from the interaction point.For photon pairs with an invariant mass within 15.75 MeV of the nominal 0 mass, we perform a kinematic fit, constraining the mass.We require that the probability of the fit exceeds 1%.The efficiencies, including the EMC acceptance, vary between 55% and 65% for 0 energies ranging from 0.3 to 2.5 GeV, for a mass resolution between 5.5 and 7.5 MeV.
We use Monte Carlo (MC) simulation of the production and decay of B mesons at the 4S resonance and of the detector response [22] to estimate signal and background efficiencies, and to extract the observed signal and background distributions to fit the data.We assume that the 4S decays exclusively to B B pairs.The simulated sample of generic B B events corresponds to roughly 3 times the B B data sample.Information from studies of selected control data samples on efficiencies and resolutions is used to improve the accuracy of the simulation.Comparisons of data with the MC simulations have revealed small differences in the tracking and particle detection efficiencies, which have been corrected for.The MC simulations include radiative effects such as bremsstrahlung in the detector material.QED final state radiation is modeled by PHOTOS [23], and decays with radiative photons are included in the signal sample.
In the MC simulations the branching fractions for hadronic B and D decays are based on average values reported in the Review of Particle Physics [13].B 0 !D ÿ ' ' signal events are generated with the HQET-based form factors, using the specific parametrization by Caprini, Lellouch, and Neubert [18].Values of the form-factor parameters are taken from measurements by the CLEO Collaboration [11].B ! D ' decays, involving orbitally excited charm mesons, are generated according to the ISGW2 model [24], and decays to nonresonant charm states are generated following the prescription of Goity and Roberts [25].

IV. EVENT SELECTION AND BACKGROUND SUBTRACTION A. Event selection
The reconstruction of the events and the selection of candidate B 0 !D ÿ ' ' decays are largely common to the earlier BABAR analysis [10].
We select events that contain a D ÿ candidate and an oppositely charged electron or muon with momentum in the range 1:2 < p ' < 2:4 GeV.Unless explicitly stated otherwise, momenta are measured in the 4S rest frame, which is boosted relative to the laboratory frame, 0:56.We reconstruct D ÿ in the decay channel D ÿ !D 0 ÿ s , with the D 0 decaying to K ÿ , K ÿ ÿ , or K ÿ 0 .The tracks of the charged hadrons from the D 0 candidate are fitted to a common vertex, and the candidate is rejected if the fit probability is less than 0.1%.We require the invariant mass of the hadrons to be compatible with the D 0 mass within 2:5 times the experimental resolution, corresponding to 34 MeV for D 0 !K ÿ 0 decays and 17 MeV for the other decays.For the decay D 0 !K ÿ 0 , we accept only candidates from portions of the Dalitz plot where the square of the decay amplitude [26] exceeds 10% of the maximum.We select D ÿ candidates with a momentum in the range 0:5 < p D < 2:5 GeV.For the ÿ s from the D ÿ decay, the momentum in the laboratory frame must be less than 450 MeV, and the momentum transverse to the beam must be greater than 50 MeV.Finally, the lepton, the ÿ s , and the D 0 are fitted to a common vertex using a constraint from the beam-beam interaction point.The probability for this fit is required to exceed 1%.
In semileptonic decays, the presence of an undetected neutrino complicates the separation of the signal from background.For a signal B 0 decay, the D ÿ and the charged lepton originate from the B 0 and the only missing particle is a massless neutrino.The absolute value of the B momentum, p B , is known from the total energy in the event, and its direction is constrained to lie on a cone centered on the D ÿ ' momentum vector.The opening angle of this cone, B;D ' , is computed for each event, Here m, E, and p refer to the mass, the energy, and the absolute value of the momentum.The condition j cos B;D ' j 1:0 should be fulfilled in a perfectly reconstructed decay.The value of w depends on the azimuthal angle of the B 0 direction, which cannot be determined.We therefore approximate w by the average of the four values of w corresponding to the azimuthal angles 0, =2, , and 3=2, as was done in [12].This approximation results in an average resolution for w of 0.04.

B. Background subtraction
The background subtraction is performed separately for each of the four kinematic variables (w, cos ' , cos V , and ), to be used for the extraction of the form-factor parameters and jV cb j.For each variable, we divide the events into ten subsets, each corresponding to one of ten bins of the distribution, and the fits are performed separately for the ten subsets.
The selected events are divided into six signal samples by separating decays into electrons and muons and the three D 0 decay modes.In addition to signal events, each subsample contains background events from six different sources: (i) combinatorial background (events from B B and continuum q q production in which at least one of the hadrons assigned to the D ÿ does not originate from the D ÿ decay); (ii) continuum background (D ÿ ' combinations from e e ÿ !c c); (iii) fake-lepton background (a true D ÿ combined with a hadron misidentified as a lepton); (iv) uncorrelated background (' and D ÿ originating from the decay of two different B mesons); (v) B background involving higher-mass charm states, Except for the combinatorial background, all background sources contain a true D ÿ !D 0 ÿ s decay and thus are expected to exhibit a peak in the m m D ÿ ÿ m D 0 distribution, where m D ÿ and m D 0 are the reconstructed masses of the D ÿ and D 0 candidates.We determine the background distributions from data, except for the correlated background, which amounts to less than 1.8% of the total sample of selected candidates.
We determine the signal and background composition in two steps.First, we estimate the combinatorial, the continuum, and the fake-lepton background from fits to the m distributions (Fig. 3).Second, we fix the background levels for these three sources and determine the uncorrelated background and the B ! D X' ' background from fits to the cos B;D ' distributions (Fig. 4).The shape and normalization of the small correlated background is fixed to the MC predictions and based on measured branching fractions [13].

Fits to m distributions
To estimate the shape and normalization of the combinatorial, the continuum, and fake-lepton backgrounds from the measured m distributions, we use, in addition to the on-resonance data, off-resonance data, as well as a set of on-resonance events in which no lepton is identified and a charged hadron is selected to take its place.We refer to this data sample as the fake-lepton sample.For each of these data sets, we perform an unbinned maximum-likelihood fit to the m distributions.The data are fitted to a sum of a peak, due to correctly reconstructed D ÿ !D 0 ÿ s decays, and a combinatorial background.
The peak is described as a sum of three Gaussian resolution functions with three mean values, three different widths, and two parameters that specify the contributions of continuum and fake-lepton backgrounds relative to the total number of events above the combinatorial background.The combinatorial background is described, as in [10], by an empirical function, where N is the normalization, m 0 refers to the kinematic threshold equal to the pion mass, and c 1 and c 2 are free parameters.
Since the m resolutions and background yields depend on the D 0 decay mode, and on whether the low-momentum pion track is reconstructed in the SVT alone or in both SVT and DCH, the parameters describing the peak contributions are determined separately for the three subsamples corresponding to the D 0 decay modes, each divided into two classes of events, depending on the detection of the slow pion.
The off-resonance distributions are scaled to the onresonance luminosity, and the simulated lepton signal and fake-lepton samples are adjusted to reproduce the detection efficiencies and misidentification probabilities determined from independently selected data control samples.For the on-and off-resonance data and for the fake-lepton sample, the parameters describing the shape and normalization of the combinatorial background differ, and they are therefore determined separately, while the FIG. 3 (color online).m distributions for the decays D ÿ !D 0 ÿ s with D 0 !K ÿ for events in which the ÿ s track is reconstructed in the SVT only, (a) for on-resonance, (b) for offresonance data in which an electron or muon was identified, and (c) for on-resonance data in which no charged lepton was found in the event.The data (points with statistical errors), integrated over the full w range, are compared to the result of the simultaneous fit (solid line) to all three distributions.The shaded area shows the fitted combinatorial background; the remainder (white area) represents the sum of the signal and peaking backgrounds.
mean and the widths of the Gaussian functions are common.
The fits extend in m from 138 to 165 MeV.Given the very large number of parameters and the many data subsamples and various subsets, the fits to the m distributions are performed in several steps.Initially, the parameters describing the combinatorial background contributions are fixed to values determined from binned 2 fits to simulated distributions.With these starting values for the combinatorial background shapes, an unbinned maximum-likelihood fit is performed to all subsets of the data, both on-and off-resonance, to determine the eight peak-shape variables and three additional relative peak yields for each data subset.To improve the agreement with the data, the parameters describing the combinatorial background distributions are then refitted, together with the three peak yield parameters, with the remaining parameters describing the components of the peaking signal and background fixed to the results of the previous fit.
As an illustration for these background fits, Fig. 3 shows the m distributions for the decays D ÿ !D 0 ÿ s with D 0 !K ÿ for events in which the ÿ s track is reconstructed in the SVT only.The data are compared to the results of the combined fit to three distributions, for se-lected samples of on-resonance data, off-resonance data, and on-resonance data without an identified lepton.
The fitted fractions of combinatorial, fake-lepton, and continuum background events are determined for the peak regions in m, which are defined as 144 to 147 MeV for events with the slow charged pion detected in the SVT and DCH, and 143 to 148 MeV for decays with the ÿ s detected in the SVT alone.

Fits to cos B;D ' distributions
In a second step, a binned 2 fit is performed to the cos B;D ' distributions in the range ÿ10 < cos B;D ' < 5, to determine the signal contribution and the normalization of the uncorrelated and B ! D X' ' background events.Neglecting resolution and radiative effects, signal events meet the constraint j cos B;D ' j 1, while the distribution of B ! D X' ' events extends below ÿ1, and the uncorrelated background events are spread over the entire range considered.Because of final state radiation, the events in the electron samples also contribute to lower values of cos B;D ' .This fit is performed separately for the six signal samples, corresponding to the D 0 mesons reconstructed in three different decay modes, divided into events with electrons or muons.Each signal sample is further divided into ten subsets corresponding to ten bins in one of the four kinematic variables.We perform the fits separately for each bin, with the individual shapes for the signal and for each of the six background sources taken from MC simulation.The fraction of the combinatorial, fake-lepton, and continuum events are taken from the m fits and fixed.
To reduce the sensitivity to statistical fluctuations, we require that the ratio of B ! D X' ' background and of uncorrelated background to the signal be the same for all three D 0 decay modes, with either electrons or muons.The cos B;D ' distributions for the six signal samples and the results of the fits are shown in Fig. 4. In total, there are 68 840 decay candidates in the range j cos B;D ' j < 1:2.The number of selected events and the fractions of background events are given in Table I.As expected, the fake rate is about a factor 10 higher for decays to muons than to electrons.Except for the combinatorial backgrounds, which vary significantly for the three D 0 decays modes, most of the other background fractions are similar for the six data samples.

V. EXTRACTION OF jV cb j AND FORM-FACTOR PARAMETERS
We determine F 1jV cb j and the three form-factor parameters by extending the one-dimensional least-squares fit to the w distribution used previously [10] to a combined fit of three one-dimensional binned distributions, with binby-bin background subtraction.We have chosen this approach to avoid the statistical limitations of a fit to a binned four-dimensional decay distribution.
In principle, any kinematic observable that is sensitive to the form-factor parameters can be used.We have examined the sensitivity of the four kinematic variables, w, cos ' , cos V , and (see Sec. II), to the form-factor parameters and found that the distribution is practically insensitive, and thus we select the remaining three variables.As done for the background estimate, we again divide the distributions into ten bins with equal bin size for w and cos V , and varying bin size for cos ' in order to have a more similar population in all the bins for this variable.
We perform a least-squares fit to these three projected one-dimensional distributions to extract the form factors and F 1jV cb j, using the independently determined background estimates.We account for the correlations by noting that the statistical covariance between the content of two bins in two different one-dimensional distributions is determined by the common number of events in these bins, while it is zero for bins in the same distribution.Since each of the three distributions that are included in the fit contain the same events and thus have the same normalization, we reduce the total number of bins used in the fit by 2, from 30 to 28.The choice of the bins that are left out is arbitrary, and it has been verified that the fit result does not depend on this choice.

A. The least-squares fit
The concept for this fit is an extension of the one introduced in the previous BABAR analysis [10].For a given bin with index i, N data i is the total number of observed events, N bkg i is the estimated number of background events, and N MC i refers to the number of MC-simulated signal events.The fit function can be written as where the indices i and j run over the 28 bins, and the index k runs over all MC-simulated events, N MC i , in bin i. W k i is a weight assigned to the kth simulated signal event in bin i to evaluate the expected signal yield as a function of the free parameters of the fit, and C ij is the covariance matrix element for a pair of bins i and j.
Each weight W k i is the product of four weights, is the product of the correction factors for efficiencies, applied on a particle-byparticle basis, which accounts for the residual differences in reconstruction and particle-identification efficiencies between the data and the Monte Carlo simulation, as a function of particle momentum and polar angle.(4) The factor W S;k i accounts for potential small differences in efficiencies among the six data subsamples and allows the adjustment of the D 0 branching fractions, properly accounting for their correlated systematic uncertainties.W S;k i is the product of several scale factors that are free parameters in the fit, each constrained to its expected value within the estimated experimental uncertainty.Specifically, to account for the uncertainty in the multiplicitydependent tracking efficiency, we introduce a factor W S trk 1 N trk 1 ÿ trk , where N trk is the number of charged tracks of the D ' candidates and trk is constrained to 1.0 within the estimated uncertainty of trk 0:8% in the single-track efficiency.
Similarly, multiplicative correction factors ' , K , and 0 are introduced to adjust the efficiencies for leptons ( e for e or for ), kaons, and 0 mesons, each within their estimated uncertainties, ' ( e for e or for ), K , and .Likewise, BK , BK , BK 0 are introduced to adjust the individual D 0 branching fractions within their current measurement uncertainties [13].Correlations between the branching fraction measurements are taken into account by the covariance matrix C B .As a result, W S;k i can be expressed as The corrections to the kaon and 0 efficiency, the decay multiplicity, N trk , and the D 0 branching fraction depend on the particular event k in bin i.The complete ansatz for the X 2 function used in the fit is The indices n, m refer to the three D 0 decay modes.The addition of these extra terms allows us to fit all subsamples simultaneously, while taking into account the correlated systematic uncertainties and effectively propagating these uncertainties, via the weights, to the uncertainties on the free parameters of the fit.
The fit procedure has been tested on a large variety of Monte Carlo generated signal samples.The fits are performed for multiple samples, comparable in size to the data.The resulting pulls are consistent with a Gaussian distribution, with no evidence for systematic biases.The width of the pull distribution has been found to be consistent with 1, giving confidence in the uncertainty extracted from the fit.

B. The covariance matrix
A direct consequence of this method is that the covariance matrix for the measurements is not diagonal, since the events in different bins are not all statistically independent.The total covariance matrix is the sum of three separate matrices: one for the measured data yields, and one each for the estimated signal and background yields.The diagonal elements of the matrices are the uncertainties of the bin contents.The covariance of bins belonging to the same distribution is zero, and the covariance of bins from different distributions is the variance of the number of events that is common for the two bins.
For the data, the covariance matrix is determined under the assumption that the data obey Poisson statistics, and therefore the variance of the bin or of the intersection of two bins is simply the number of events.The estimated signal matrix is built in an analogous way, where the variance of a number of weighted events n is approximated by the sum of the squares of the weights, P i1;n w 2 i .The calculation of the background covariance matrix is less straightforward.The diagonal elements are simply the estimated variances of the measured background, according to the procedure described in Sec.IV.However, the background extraction procedure does not directly determine the number of common events in two bins, because this procedure is based on a rather complex sequence of fits to shapes and event yields.The solution adopted here is to use the number of common background events in two bins as predicted by the simulation, corrected for tracking and particle-identification (PID) efficiencies and adjustments to account for the background estimates.This is done in such a way that for each bin of one of the kinematic observables the background is equal to the data-based background estimate.The choice of the observable (w, cos ' , cos V , or ) is arbitrary, and the variation of the results with this choice is used to evaluate the systematic uncertainty introduced by this procedure.
Since the total number of background events estimated for the four different distributions is not exactly the same, we average the results obtained for the four distributions.The spread of the background normalization values is found to be almost twice as large as the estimated uncertainty from the error propagation.This can be explained by the fact that the uncertainties of the assumed signal and background shapes in the m fits are not accounted for in this propagation.We derive an estimate of this additional uncertainty using the number of background events determined from the fits to the background distributions.They are very similar for w and cos ' , and for cos V and .The estimate is given by the average of the minimum and maximum difference between these two groups.

VI. THE FIT RESULTS AND SYSTEMATIC UNCERTAINTIES A. Results of the fit
The results of the simultaneous fit to the three onedimensional distributions are presented in Table II.The stated uncertainties on 2 , R 1 1, R 2 1, and F 1jV cb j are taken from the MINUIT minimization program [27].Among the form-factor parameters, the correlations are quite large, but their correlation with F 1jV cb j is less than 0.23.All the parameters are found to be consistent with their nominal value of 1 within their uncertainties.
Figure 5 compares the one-dimensional projections for the kinematic observables with the results of the fit, detailing the signal and background contributions.To provide an additional check on the background estimation, we also show the projection.The goodness of the fit, ignoring that the measurements used are not all independent, can be stated as 2 =d:o:f 23:8=24, corresponding to a probability of 47.3%.
As a cross check, we perform the fit separately for the six subsamples.The quality of the fits is generally good, and the results, shown in Table III, agree within the uncertainties obtained, with the possible exception of the Ke sample.Detailed investigations of the background estimates and fits for this sample did not reveal any anomalies.

B. Systematic uncertainties
A summary of statistical and systematic uncertainties on the measured parameters is presented in Table IV, including the breakdown of those for the measurement of the B 0 !D ÿ ' ' branching fraction.

Uncertainties included in the fit
The uncertainty of the parameters resulting from the fit is not purely statistical, since the systematic uncertainty sources that are not common to all events are accounted for in the fit through the parameters in the weights W S;k i .As described above, these weights account for residual uncertainties in the lepton and hadron identification, the charged particle tracking and 0 efficiencies, and the individual D 0 branching fractions.We can determine the statistical uncertainties of the fit by repeating the fit with all parameters fixed at their fitted value.The systematic uncertainties are then obtained by subtracting the statistical covariance matrix from the total covariance matrix.While the uncertainties related to the detector performance are relatively TABLE II.The results of the fit: The parameters, their uncertainties, and the off-diagonal elements of the correlation matrix.The stated uncertainties include systematic ones that are not common to all events and are therefore included in the fit (see Sec. VI B 1).

Soft pion efficiency
A major source of uncertainty on F 1jV cb j is the reconstruction efficiency for the low-momentum pion from the D ÿ decay, since it is highly correlated with the D ÿ momentum and thereby with w.The functions parametrizing the efficiency for data and MC simulation are consistent within the statistical uncertainties.To assess the systematic uncertainty on jV cb j, we vary the parameters of the efficiency function by their uncertainty, including correlations.We add in quadrature the uncertainty in the absolute scale as determined from higher-momentum tracks reconstructed in both the SVT and the DCH.The resulting systematic uncertainty on jV cb j is 1.3%.

D l vertex reconstruction efficiency
The uncertainties from the D l vertex reconstruction have been evaluated by observing the impact of changes  III.Results of fits performed separately for the six subsamples corresponding to each combination of three D 0 decay modes and the charged lepton.The uncertainties represent the total uncertainty of the fit, except for F 1jV cb j, where it is split into the statistical and the systematic contribution included in the fit. in the standard vertex fit procedure.First, we remove the constraint on the average position of the beam-beam interaction point, and second, we remove the lepton track from the vertex fit.We take the larger of the observed variations of the parameters as the systematic uncertainty.

B momentum
The B momentum is impacted by small changes in the energies of the two colliding beams.In the Monte Carlo simulation the beam energies are assumed to be constant.We have examined the impact of these variations on the measured distributions and have concluded that we can account for this difference by rescaling the cos B;D ' values in the simulation by a factor of 0.97.Half of the observed relative change of the fitted parameters coming from the adjustment is assumed as the systematic uncertainty due to this effect.This is the largest systematic uncertainty on R 1 1 and R 2 1.

Radiative corrections
Radiative corrections to the B 0 !D ÿ ' ' decays are simulated by PHOTOS 2.0 [23], which describes the final state photon radiation (FSR) up to O 2 .In the event reconstruction, no attempt is made to recover photons emitted in the decay.The simulated prediction of the reconstructed kinematic variables is sensitive to the details of the radiative corrections.This is particularly important for electrons, for which FSR results in the long tail below ÿ1 in the cos B;D ' distribution.
At present, no detailed calculation of the full O radiative corrections to the B 0 !D ÿ ' ' decay is available.Recently, a new O calculation of radiative corrections in K 0 !ÿ e e decays has become available [28], which allows detailed comparisons of the radiated photons with data and with PHOTOS calculations.These new calcu-lations agree well with kaon data.From the comparison with PHOTOS, it is evident that the radiated photon energy spectrum is reproduced quite well by PHOTOS, but the photon emission angle with respect to the electron differs significantly.
To assess the systematic uncertainty due to the imperfect treatment of the radiative corrections in B decays, we have used the comparison presented in Ref. [28] and reweighted the simulated decays to reproduce the photon angular distribution for photons above 10 MeV in the B rest frame.The effect of the reweighting has been used as an estimate of the systematic uncertainty.The possible impact of the radiative effects on the hadronic interactions, i.e. the form factors, is unknown and therefore not considered.

B ! D X' ' background description
The semileptonic branching fraction and form factors for the four higher-mass charm states, D , and for nonresonant production of D and D are not well known.The shapes of these different components of the B ! D X' ' background are taken from simulation.Their relative yield is obtained from the fit to the cos B;D ' distribution.To account for the uncertainty of the composition of this background, the fits have been repeated using only one of the contributions at a time.The study was done only for the D 0 !K ÿ subsample; it is assumed to be valid also for the other subsamples.The estimated uncertainty on the fit parameters is taken as half of the biggest change observed.

Background estimates
As explained in Sec.V B the background covariance matrix is built using the measured background shape in one observable, and the simulation information for the others.The choice of the observable is arbitrary.The systematic uncertainty is evaluated by comparing results for the four kinematic observables.The maximum observed variation with respect to the standard fit result is taken as the estimate for this uncertainty.

Global normalization factors
There are several quantities that affect only the overall normalization of the data, and thus not the form-factor parameters.Their contributions to the systematic uncertainty on F 1jV cb j and the branching fraction are listed in the bottom half of Table IV.They are as follows: the B 0 lifetime B 0 1:530 0:009 [13], the D ÿ !D 0 ÿ branching fraction BD ÿ !D 0 ÿ 67:7 0:5%, and N B B , the number of B B events in the total data sample.The systematic uncertainty on N B B is 1.1%.The effect of the uncertainties in the D 0 branching fractions has already been discussed; it is subsample specific, but it affects all parameters because it changes the fraction of signal events from different D 0 decays.The uncertainty on the ratio f ÿ =f 00 B4S !B B ÿ =B4S !B 0 B 0 1:037 0:029 [13] affects both the absolute number of measured B 0 !D ÿ ' ' decays and the ratio of background events from B 0 and B decays in the simulation.This second aspect influences the cos B;D ' distributions and therefore the background determination.For this reason, this uncertainty also affects the form-factor parameters.The systematic uncertainty is equated with the observed change in the fit parameters for a 1 standard deviation change in the value of f ÿ =f 00 .

A. Summary of results of this analysis
A sample of about 52 800 fully reconstructed B 0 !D ÿ ' ' decays recorded by the BABAR detector has been analyzed to extract both F 1jV cb j and the formfactor parameters, 2 , R 1 1, and R 2 1, in the Caprini-Lellouch-Neubert parametrization [18].The D ÿ candidates are reconstructed from the D ÿ !D 0 ÿ decays, and the D 0 mesons are reconstructed in three different decay modes, K ÿ , K ÿ ÿ , and K ÿ 0 .Electrons or muons are paired with the D ÿ to form signal candidates.The large data sample has permitted a more precise determination of the background contributions, largely based on data, and thus has resulted in smaller experimental uncertainties.

B. Combination of results with the previous BABAR measurement of the form-factor parameters
The BABAR Collaboration recently published a measurement [12] of the same form-factor parameters for B 0 !D ÿ ' ' decays based on an unbinned maximumlikelihood fit to the four-dimensional decay distribution [Eq.(10)].This fit is sensitive to the interference of the three helicity amplitudes and thus results in significant smaller uncertainties on the form-factor parameters.The fit does not attempt an absolute normalization of the decays, and thus is not sensitive to F 1jV cb j.It resulted in 2 1:145 0:066 0:035, R 1 1 1:396 0:070 0:027, and R 2 1 0:885 0:046 0:013.
We combine the two BABAR measurements of the formfactor parameters, taking into account the correlation between them, and obtain F 1jV cb j 34:4 0:3 1:1 10 ÿ3 ;  tainties of the form-factor parameters.The event sample and the sample of Monte Carlo simulated events used in Ref. [12] are a subset of the one used in the present analysis, namely, about 15 000 selected B 0 candidates with D 0 !K ÿ decays combined with electrons.Except for the selection of the D 0 decay, the event selection and the determination of the background shapes and the signal extraction are almost identical for the two analyses.Therefore, all the detector-related systematic uncertainties should be the same, as well as the uncertainties from the background models and those due to input parameters as the branching fractions.Thus, we retain the systematic measurement uncertainties established in this paper.The combined statistical errors are still larger than the total systematic uncertainties, but not by a large factor.An upper limit for the correlation between the two measurements has been estimated on the basis of the ratio of the uncertainties, and is found to be less than 0.45.
The correlation coefficients for the combined measurements are 2 ; R 1 1 71%; 2 ; R 2 1 ÿ83%; 2 ; F 1jV cb j 27%; R 1 1; R 2 1 ÿ84%; R 1 1; F 1jV cb j ÿ39%; R 2 1; F 1jV cb j 22%: Figure 7 shows the correlations between the fitted variables and their uncertainties, both for the present analysis and for the combined result with Ref. [12].The contours correspond to 2 1, i.e. 39% CL.The correlations between the form-factor parameters are quite large, but their correlation with F 1jV cb j is less than 0.4, and the sign of the coefficients differs, resulting in a much reduced overall dependence of F 1jV cb j on these form factors.
The corresponding branching fraction of the decay B 0 !D ÿ ' ' is B B 0 !D ÿ ' ' 4:69 0:04 0:34%: The combined results of the two BABAR analyses supersede all previous BABAR measurements of the form-factor parameters, of the exclusive branching fraction for the B 0 !D ÿ ' ' decay, and of jV cb j extracted from this decay.
The value of the branching fraction presented here is smaller than the average of previous measurements [13].This measurement combined with BB 0 !D ÿ ' ' 2:08 0:18% [13] represents only 65 7% of the total branching fraction for the B 0 !X c ' ' decays.The remaining fraction of 35% is expected to involve highermass charm states.The branching fractions for decays to these individual higher-mass states are not well known, in particular, those involving broad resonances or nonresonant D ; states [29,30].
The combination of the two BABAR measurements results in a further reduction of the form-factor uncertainties compared to the previous BABAR analysis [12], for which the uncertainties on R 1 1 and R 2 1 had already been reduced by a factor of 4 or more, compared to the CLEO measurement [11].The uncertainty on 2 has also been reduced, by a factor of 5 with respect to the BABAR measurement in Ref. [10].The correlation between F 1jV cb j and 2 , which was sizable for all previous measurements, has been reduced significantly, and this also leads to a smaller uncertainty for jV cb j.
The resulting value of jV cb j is fully compatible with the earlier BABAR measurement [10], and most earlier measurements [13], but it is significantly smaller than the CLEO measurement [31].
jV cb j can also be extracted from inclusive B ! X c ' ' decays.Recent measurements of jV cb j rely on moments of the electron energy and the hadron mass spectrum, combined with moments of the energy photon spectrum in inclusive B ! X s decays.Here X c and X s refer to charm and strange hadronic states, resonant or nonresonant.Global fits to such moments, as a function of the minimum lepton and photon energy, have been performed in terms of two different QCD calculations, one in the so-called 1S scheme [32] and the other in the kinetic scheme [33].The results are in very good agreement.The weighted average is jV cb j 41:7 0:7 10 ÿ3 [13]; the experimental and theoretical uncertainties are comparable in size.The experimental techniques and the theoretical calculations for exclusive and inclusive decays are completely different and uncorrelated.Given the sizable experimental error and the large uncertainty of the form-factor normalization of the exclusive measurement, the results are consistent.This should give us confidence in the experimental methods employed and the theoretical calculations, but also an incentive to improve on both.

FIG. 4 (
FIG. 4 (color).The cos B;D ' distributions (points with statistical errors) for selected events in the six subsamples (integrated over the full range in w): (a,b) D !K, (c,d) D !K, and (e,f) D !K 0 , compared to the result of the fits to signal and background contributions.

FIG. 5 (
FIG. 5 (color).Comparison of the measured distributions (data points), (a) w, (b) cos ' , (c) cos V , and (d) , with the result of the fit, shown as the sum of the fitted signal yield and the estimated background distributions.The statistical uncertainties of the data are too small to be visible.

3 FIG. 6 .
FIG.6.The measured w dependence of F wjV cb j (data points) compared to the theoretical function with the fitted parameters (solid line).The experimental uncertainties are too small to be visible.

TABLE I .
Number of selected candidates and background fractions, separately for the subsamples identified by the D 0 decay mode and the charged lepton.

TABLE IV .
Breakdown of statistical and systematic uncertainties.