Measurement of the hadronic form factor in D-0 -> K(-)e(+)nu(e) decays 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, 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. 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, 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, 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, M. J. Kobel, 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, K. A. George, F. Di Lodovico, 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, G. Calderini, J. Chauveau, P. David, L. Del Buono, Ch. de la Vaissière, O. Hamon, Ph. Leruste, J. Malclès, J. Ocariz, A. Perez, L. Gladney, M. Biasini, R. Covarelli, E. Manoni, C. Angelini, G. Batignani, S. Bettarini, M. Carpinelli, R. Cenci, A. Cervelli, F. Forti, M. A. Giorgi, A. Lusiani, PHYSICAL REVIEW D 76, 052005 (2007)

by the BABAR detector at the PEP-II electron-positron collider. The corresponding decay branching fraction, relative to the decay D 0 ! K ÿ , has also been measured to be R D BD 0 ! K ÿ e e =BD 0 ! K ÿ 0:927 0:007 0:012. From these results, and using the present world average value for BD 0 ! K ÿ , the normalization of the form factor at q 2 0 is determined to be f 0 0:727 0:007 0:005 0:007 where the uncertainties are statistical, systematic, and from external inputs, respectively. DOI

I. INTRODUCTION
Measurements of exclusive semileptonic D decays provide an accurate determination of the hadronic form factors entering in these decays. Assuming that the CKM matrix is unitary, the elements jV cs j and jV cd j can be determined: In D 0 ! K ÿ e e decays [2], with a pseudoscalar hadron emitted in the final state, and neglecting the electron mass, the differential decay rate depends only on one form factor f q 2 , dÿ dq 2 G 2 F 24 3 jV cs j 2 jp K q 2 j 3 jf q 2 j 2 ; (2) where G F is the Fermi constant, q 2 is the invariant mass squared of the two leptons, e and e , andp K q 2 is the kaon three-momentum in the D 0 rest frame [3]. In this paper we present measurements of the q 2 variation and absolute value of the hadronic form factor at q 2 0 for the decay D 0 ! K ÿ e e . The data consist of D mesons produced in e e ÿ ! c c continuum events at a center of mass energy near the 4S mass, and were recorded by the BABAR detector at the Stanford Linear Accelerator Center's PEP-II collider. A semi-inclusive reconstruction technique is used to select charm semileptonic decays with high efficiency. As a result of this approach, events with a photon radiated during the D 0 decay are included in the signal. The systematic uncertainties are kept as low as possible by using control samples extracted from data where possible.
Measurements of D ! K ' ' , based on smaller signal events samples, have been published by the CLEO [4], FOCUS [5], and Belle [6] Collaborations. FOCUS (E831) is a Fermilab heavy flavour photoproduction experiment whereas CLEO and Belle have analyzed events from e e ÿ annihilation at an energy close to the 4S mass. This paper is organized as follows. A general description of the hadronic form factor, f q 2 , is given in Sec. II, where the different parameterizations considered in this analysis are explained. In Section III a short description of the detector components that are important to this measurement is given. The selection of signal events and the rejection of background are considered in Section IV. In Section V, the measured q 2 variation of the hadronic form factor is discussed and compared with previous measurements. In Sec. VI the measured decay rate is given and in Sec. VII these measurements are combined to obtain the value of f 0.

II. THE F Q 2 HADRONIC FORM FACTOR
The amplitude for the decay D 0 ! K ÿ ' ' depends on two hadronic form factors: where V s c. The constraint f 0 f 0 0 ensures that there is no singularity at q 2 0. When the charged lepton is an electron, the contribution from f 0 is proportional to m 2 e and can be neglected in decay rate measurements.
The parameterizations of f q 2 which have been compared with present measurements and a few examples of theoretical approaches, proposed to determine the values of corresponding parameters, are considered in the following.

A. Form-factor parameterizations
The most general expressions of the form factor f q 2 are analytic functions satisfying the dispersion relation: Imf t t ÿ q 2 ÿ i : (4) The only singularities in the complex t q 2 plane originate from the interaction of the charm and the strange quarks in vector states. They are a pole, situated at the D s mass squared and a cut, along the positive real axis, starting at threshold (t m D m K 2 ) for D 0 K ÿ production.

Taylor expansion
This cut t-plane can be mapped onto the open unit disk with center at t t 0 using the variable: t ÿ t 0 p : (5) In this variable, the physical region for the semileptonic decay (0 < t < t ÿ q 2 max m D ÿ m K 2 ) corresponds to a real segment extending between z max 0:051. This value of z max is obtained for The z expansion of f is thus expected to converge quickly. The most general parameterization [7], consistent with constraints from QCD, is which is based on earlier considerations [8]. The function Pt zt; m 2 D s has a zero at the D s pole mass and jPj 1 along the unit circle; is given by where V can be obtained from dispersion relations using perturbative QCD and depends on u m s =m c [9]. At leading order, with u 0 [10], The choice of P and is such that Having measured the first coefficients of this expansion, Eq. (9) can constrain the others. This constraint, which depends on V , may have to be abandoned in the case of charm decays as the charm-quark mass may not be large enough to prevent the previous evaluation of V from receiving large 1=m c and QCD corrections. However the parameterization given in Eq. (6) remains valid and it has been compared [7] with available measurements. The first two terms in the expansion were sufficient to describe the data.

Model-dependent parameterizations
A less general approach assumes that the q 2 variation of f q 2 is governed mainly by the D s pole and that the other contributions can be accounted for by adding another effective pole at a higher mass [11]: with pole 1= pole ÿ pole =1 ÿ pole and pole 1= pole . If in addition, the form factors f and f 0 must obey a relation, valid at large recoil and in the heavy quark limit, then pole 1= pole [11] ( pole pole and pole 0 in this case). Equation (10) becomes: known as the modified pole ansatz. Initially an even simpler expression, the simple pole ansatz, was proposed which considered only the contribution from the D s pole. In the following, the pole mass entering in is fitted. Note that such an effective pole mass value has no clear physical interpretation and that the proposed q 2 variation does not comply with constraints from QCD. The obtained value may nonetheless be useful for comparison with results from different experiments.

B. Quantitative expectations
Values of the parameters that determine f q 2 were obtained initially from constituent quark models and from QCD sum rules. These two approaches have an intrinsically limited accuracy. In this respect, results from lattice QCD computations are more promising because their accuracy is mainly limited by available computing resources.

Quark models
Quark model calculations estimate meson wave functions and use them to compute the matrix elements that appear in the hadronic current. There are a large variety of theoretical calculations [12]. Among these models we have selected the ISGW model [13], simply because it is widely used to simulate heavy hadron semileptonic decays. This model was expected to be valid in the vicinity of q 2 max , a region of maximum overlap between the initial and final meson wave functions. In ISGW2 [14] the exponential q 2 dependence of the form factor has been replaced by another parameterization, with a dipole behavior, expected to be valid over a larger q 2 range: The predicted values of the parameters are f q 2 max 1:23 and r 1:12 GeV ÿ1 for the charge radius [14].

QCD sum rules
QCD sum rules [15] and their extension on the light cone [16], are expected to be valid at low q 2 . Using a value of 150 MeV for the strange quark mass, one obtains [16]: f 0 0:78 0:11 and pole ÿ0:07 0:15 ÿ0:07 ; (14) using the modified pole ansatz. The uncertainty of f 0) is estimated to be of order 15%, and the q 2 dependence is expected to be dominated by a single pole at the D s mass because the value of pole is compatible with zero.

Lattice QCD
Lattice QCD computation is the only approach able to compute f q 2 from first principles. Current results must be extrapolated to physical values of light-quark masses and corrected for finite lattice size and discritization effects. There have been several evaluations of f q 2 for different values of the momentum transfer in the quenched approximation [17,18]. These results have been combined [17], giving f 0 0:73 0:07. The first unquenched calculation has been published recently [19]: f 0 0:73 0:03 0:07 and pole 0:50 0:04, using the modified pole ansatz to parameterized the q 2 dependence of the form factor.

C. Analyzed parameterizations
The different parameterizations of f q 2 considered in this analysis are summarized in Table I, along with their corresponding parameters and expected values, where available.

III. THE BABAR DETECTOR AND DATASET
A detailed description of the BABAR detector [20] and of the algorithms used for charged and neutral particle reconstruction and identification [21] is provided elsewhere. Charged particles are reconstructed by matching hits in the 5-layer double-sided silicon vertex tracker (SVT) with track elements in the 40-layer drift chamber (DCH), which is filled with a gas mixture of helium and isobutane. Low momentum particles which do not leave enough hits in the DCH due to the bending in the 1.5 T magnetic field, are reconstructed in the SVT. Charged hadron identification is performed combining the measurements of the energy deposition in the SVT and in the DCH with the information from the Cherenkov detector (DIRC). Photons are detected and measured in the CsI(Tl) electro-magnetic calorimeter (EMC). Electrons are identified by the ratio of the track momentum to the associated energy deposited in the EMC, the transverse profile of the shower, the energy loss in the DCH, and the Cherenkov angle in the DIRC. Muons are identified in the instrumented flux return, composed of resistive plate chambers interleaved with layers of steel and brass.
The results presented here are obtained using a total integrated luminosity of 75 fb ÿ1 registered by the BABAR detector during the years 2000 -2002. Monte Carlo (MC) simulation samples of 4S decays, charm and other light-quark pairs from continuum equivalent, respectively, to 2.8, 1.2 and 0.7 times the data statistics, respectively, have been generated using GEANT4 [22]. These are used mainly to evaluate background components. Quark fragmentation, in continuum events, is described using the JETSET package [23]. The MC distributions have been rescaled to the data sample luminosity, using the expected cross sections of the different components (1.3 nb for c c, 0.525 nb for B B ÿ and B 0 B 0 , 2.09 nb for light u u, d d and s s quark events). Dedicated MC samples of pure signal events, equivalent to 7 times the data statistics, are used to correct measurements for efficiency and finite resolution effects. They have been generated using the modified pole parameterization ansatz for f q 2 with pole 0:50. Radiative decays (D 0 ! K ÿ e e ) are modeled by PHOTOS [24]. To account for one of the most important sources of background, a special event sample with, in each event, at least one cascade decay D ! D 0 , D 0 ! K ÿ 0 e e (or its charge conjugate) has been generated with a parameterization of the form factors in agreement with measurements from the FOCUS Collaboration [25]. Events with a D and a D 0 decaying into K ÿ or K ÿ 0 have been reconstructed in data and simulation. These control samples have been used to adjust the c-quark fragmentation distribution and the kinematic characteristics of particles accompanying the D meson in order to better match the data. They have been used also to measure the reconstruction accuracy on the missing neutrino momentum. In addition, K ÿ events are used as a reference channel in the D 0 ! K ÿ e e branching fraction measurement.

IV. SIGNAL RECONSTRUCTION
We reconstruct D 0 ! K ÿ e e decays in e e ÿ ! c c events where D 0 originates from the D ! D 0 . The

Modeling
Parameters Expected values z expansion [8] a 0 , r k , k 1 No prediction General two poles [11] f 0, pole , pole No prediction Modified pole [11] f 0, pole pole 0 Simple pole f 0, m pole m pole m D s ISGW2 [14] f t ÿ , I f t ÿ 1:23 I 0:104 GeV ÿ2 main sources of background arise from events with a kaon and electron candidate. Such events come from 4S decays and the continuum production of charmed hadrons. Their contribution is reduced using variables sensitive to the particle production characteristics that are different for signal and background events.

A. Signal selection
Charged and neutral particles are boosted to the center of mass system (c.m.) and the event thrust axis is determined. The direction of this axis is required to be in the interval j cos thrust j < 0:6 to minimize the loss of particles in regions close to the beam axis. A plane perpendicular to the thrust axis is used to define two hemispheres, equivalent to the two jets produced by quark fragmentation. In each hemisphere, we search for pairs of oppositely charged leptons and kaons. For the charged lepton candidates we consider only electrons or positrons with c.m. momentum greater than 0:5 GeV=c.
Since the e momentum is unmeasured, a kinematic fit is performed, constraining the invariant mass of the candidate K ÿ e e system to the D 0 mass. In this fit, the D 0 momentum and the neutrino energy are estimated from the other, charged and neutral, particles measured in the event.
Masses of charged particles are assigned according to the information provided by particle identification algorithms. The D 0 direction is taken as the direction opposite to the sum of the momenta of all reconstructed particles in the event, except for the kaon and the positron associated with the signal candidate. The energy of the jet is determined from the total c.m. energy and from the measured masses of the two jets. The neutrino energy is estimated as the difference between the total energy of the jet containing the kaon and the electron and the sum of the particle energies in that hemisphere. A correction, which depends on the value of the missing energy measured in the opposite jet, is applied to account for the presence of missing energy due to particles escaping detection, even in the absence of a neutrino from the D 0 decay.
The D 0 candidate is retained if the 2 probability of the kinematic fit exceeds 10 ÿ3 . Detector performance for the reconstruction of the D 0 direction and for the missing energy are measured using events in which the D 0 decays into K ÿ . Corrections are applied to account for observed differences between data and simulation. Each D 0 candidate is combined with a charged pion, with the same charge as the lepton, and situated in the same hemisphere. The mass difference m mD 0 ÿ mD 0 is evaluated and is shown in Fig. 1. This distribution contains events which in addition pass the requirements on the Fisher discriminant F B B suppressing B B background and also give a satisfactory kinematic fit constraining the invariant D mass which is described in Section IV C. This last requirement is the reason of the slow decrease of the m distribution. At large m values, there is a small excess of background in data and the simulation is rescaled accordingly, in Sec. VA 5, to determine the signal component. Only events with m < 0:16 GeV=c 2 are used in the analysis.

B. Background rejection
Background events arise from 4S decays and hadronic events from the continuum. Three variables are used to reduce the contribution from B B events: R 2 (the ratio between the second and zeroth order Fox-Wolfram moments [26]), the total charged and neutral multiplicity and the momentum of the soft pion ( s ) from the D .
These variables exploit the topological differences between events with B decays and events with c c fragmentation. The particle distribution in 4S decay events tends to be isotropic as the B mesons are produced near threshold, while the distribution in c c events is jetlike as the c.m. energy is well above the charm threshold. This also results in a softer D momentum spectrum in 4S decays compared to c c events. Corresponding distributions of these variables for signal and background events are given in Fig. 2. These variables have been combined linearly in a Fisher discriminant. The requirement F B B > 0:5 retains 65% of signal and 6% of B B-background events. Background events from the continuum arise mainly from charm particles as requiring an electron and a kaon reduces the contribution from light-quark flavors to a low level. Because charm hadrons take a large fraction of the charm-quark energy charm decay products have higher average energies and different angular distributions (relative to the thrust axis or to the D direction) compared with other particles in the hemisphere emitted from the hadronization of the c and c quarks. These other particles are referred to as ''spectator'' in the following; the ''leading'' particle is the one with the largest momentum. To reduce background from c c events, the following variables are used: (i) the D 0 momentum; (ii) the spectator system mass, m sp , which has lower values for signal events; (iii) the direction of the spectator system momentum relative to the thrust axis cos sp-thrust ; (iv) the momentum of the leading spectator track; (v) the direction of the leading spectator track relative to the D 0 direction; (vi) the direction of the leading spectator track relative to the thrust axis; (vii) the direction of the lepton relative to the kaon direction, in the e e rest frame, cos e ; (viii) the charged lepton momentum, jp e j, in the c.m.
frame. The first six variables depend on the properties of c-quark hadronization whereas the last two are related to decay characteristics of the signal. Distributions for four of the most discriminating variables are given in Fig. 3. D 0 ! K ÿ events have been used to tune the simulation parameters so that distributions of the variables used to reject background agree with those measured with data events. These eight variables have been combined linearly into a Fisher discriminant variable (F cc ) and events have been kept for values above 0. This selection retains 77% of signal events that were kept by the previous selection requirement and rejecting 66% of the background (Fig. 4).
The remaining background from c c-events can be divided into peaking (60%) and nonpeaking (40%) candidates. Peaking events are those background events whose distribution is peaked around the signal region. These are mainly events with a real D in which the slow is included in the candidate track combination. Backgrounds from e e ÿ annihilations into light u u, d d, s s quarks and B B events are nonpeaking. These components, from simulation, are displayed in Fig. 1.

C. q 2 measurement
To improve the accuracy of the reconstructed D 0 momentum, the nominal D mass is added as a constraint in the previous fit and only events with a 2 probability higher than 1% are kept ( Fig. 1 is obtained requiring only that the fit has converged). It is verified that the sample composition of the nonpeaking background does not depend on this last condition. The measured q 2 r distribution, where q 2 r p D ÿ p K 2 , is given in Fig. 5. There are 85 260 selected D 0 candidates containing an estimated number of 11280 background events. The nonpeaking component comprises 54% of the background.
To obtain the true q 2 distribution, the measured one has to be corrected for selection efficiency and detector resolution effects. This is done using an unfolding algorithm based on MC simulation of these effects.
The variation of the selection efficiency as a function of q 2 is given in Fig. 6. The resolution of the q 2 measurement for signal events is obtained from MC simulation. The resolution function can be fitted by the sum of two Gaussian functions, with standard deviations 1 0:066 GeV 2 and 2 0:219 GeV 2 , respectively. The narrow component corresponds to 40% of the events.
To obtain the unfolded q 2 distribution for signal events, corrected for resolution and acceptance effects, the Singular Value Decomposition (SVD) [27] of the resolution matrix has been used. This method uses a twodimensional matrix which relates the generated q 2 distribution to the detected distribution, q 2 r , as input. After subtracting the estimated background contribution, the measured binned q 2 r distribution is linearly transformed values have zero mean and standard deviation equal to unity [27]. Using toy simulations, we find that seven SV have to be kept with events distributed over ten bins in q 2 . Because the measurement of the form-factor parameters relies on the measured q 2 r distribution, it does not require unfolding, and is independent of this particular choice.

V. RESULTS ON THE q 2 DEPENDENCE OF THE HADRONIC FORM FACTOR
The unfolded q 2 distribution, normalized to unity, is presented in Fig. 7 and in Table II. Also given in this table are the statistical and total uncertainties and the correlations of the data in the ten bins. Figure 7 shows the result of fits to the data for two parameterizations of the form factor with a single free parameter, the simple pole and the modified pole ansatz. Both fitted distributions agree well with the data.
A summary of these and other form-factor parameterizations is given in Table III. These results will be discussed in detail in Sec. V B.
The fit to a model is done by comparing the number of events measured in a given bin of q 2 with the expectation from the exact analytic integration of the expression TABLE II. Statistical and total uncertainty matrices for the normalized decay distribution (corrected for acceptance and finite resolution effects) in ten bins of q 2 from 0 to 2 GeV 2 , and for the ratio R D (see Section VI). The total decay distribution has been normalized to unity for q 2 varying over ten 0:2 GeV 2 intervals. The uncertainty matrices are provided for both the statistical (upper half ) and total (lower half ) uncertainties. The uncertainty on each measured value ( i ) is given along the diagonal. Off-diagonal terms correspond to the correlation coefficients ( ij ). From this information, the error covariance matrix elements correspond to: ii 2 i and ij ij i j and the matrix is symmetric relative to the diagonal. jp K q 2 j 3 jf q 2 j 2 over the bin range, with the overall normalization left free. The result of the fit corresponding to the parameterization of the form factor using two parameters [see Eq. (10)] is given in Fig. 8.

A. Systematic uncertainties
Systematic uncertainties of the form-factor parameters are likely to originate from imperfect simulation of c-quark fragmentation and the detector response, from uncertainties in the background composition and the individual contributions for the selected signal sample, the uncertainty in the modeling of the signal decay and the measurement of the q 2 distribution. We study the origin and size of various systematic effects, correct the MC simulation, if possible, assess the impact of the uncertainty of the size of correction on the fit results, and adopt the observed change as a contribution to the systematic uncertainty on the fitted parameters for the different parameterizations under study. Some of these studies make use of standard BABAR evaluations of detection efficiencies, others rely on special data control samples, for instance hadronic decays D 0 ! K ÿ or K ÿ 0 .

c-quark hadronization tuning
The signal selection is based on variables related to c-quark fragmentation and decay properties of signal events. Simulated events have been weighted to agree with the distributions observed in data. Weights have been obtained using events with a reconstructed D 0 decaying into K ÿ . After applying these corrections, the distribution of the Fisher discriminant that contains these variables is compared for data and simulation. The remaining relative difference, which is below 5% over the range of this quantity, is used to evaluate the corresponding systematic uncertainty. It corresponds to the variations on fitted quantities obtained by correcting or not for this difference.

Reconstruction algorithm
It is important to verify that the q 2 variation of the selection efficiency is well described by the simulation. This is done by analyzing D 0 ! K ÿ 0 as if they were K ÿ e e events. The two photons from the 0 are removed and events are reconstructed using the algorithm applied to the semileptonic D 0 decay. The ''missing'' 0 and the charged pion play, respectively, the roles of the neutrino and the electron. To preserve the correct kinematic limits, it is necessary to take into account that the ''fake'' neutrino has the 0 mass and that the fake electron has the mass.
Data and simulated events, which satisfy the same analysis selection criteria as for Ke e , have been compared. For this test, the cos e and p e variables are removed from the Fisher discriminant, because distributions for these two variables in K ÿ 0 events are different from the signal events.
The ratio of efficiencies measured in data and simulation is fit with a linear expression in q 2 . The corresponding slope, 0:71 0:68%, indicates that there is no significant bias when the event selection criteria are applied. The measured slope is used to define a correction and to estimate the corresponding systematic uncertainty.

Resolution on q 2
To measure possible differences between data and simulation on the q 2 reconstruction accuracy, D 0 ! K ÿ 0 events are used again. Distributions of the difference q 2 r ÿ q 2 , obtained by selecting events in a given bin of q 2 are compared. These distributions are systematically slightly narrower for simulated events and the fraction of events in the distant tails are higher for data (see Fig. 9).
With the D 0 ! K ÿ sample we study, in data and simulation, the accuracy of the D 0 direction and missing energy reconstruction for the D 0 ! K ÿ e e analysis. This information is used in the mass-constrained fits and thus influences the q 2 reconstruction. Once the simulation is tuned to reproduce the results obtained on data for these parameters, the q 2 resolution distributions agree very well, as shown in Fig. 9. One half of the measured variation on the fitted parameters from these corrections has been taken as a systematic uncertainty.

Particle identification
Effects from a momentum-dependent difference between data and simulated events on the charged lepton and on the kaon identification have been evaluated. Such differences, which are typically below 2%, have been measured for selected, high purity samples of electrons and kaons. These corrections have been applied and the observed variation has been taken as the estimate of the systematic uncertainty.

Background estimate
The background under the D signal has two components that have, respectively, nonpeaking and peaking behavior.
The nonpeaking background originates from non-c c events and from continuum charm events in which the s candidate does not come from a cascading D . By comparing data and simulated event rates for m > 0:18 GeV=c 2 (see Fig. 1), a correction of 1.05 is determined from simulation for the nonpeaking background. This correction is applied and an uncertainty of 0:05 is used as the corresponding systematic uncertainty. The effect of a different variation of the nonpeaking background level with q 2 r has been evaluated by comparing corresponding distributions in data and simulation for events selected with m between 0.18 and 0:20 GeV=c 2 . The measured difference, of 5% over the q 2 r range, is used to determine the corresponding systematic uncertainty.
Events which include a slow pion originating from D decay contribute in several ways to the peaking background. The production rate of D mesons in the simulation is in agreement with expectations based on measurements from CLEO [28]. The uncertainty of 0:06 on this comparison is dominated by the systematic uncertainty from the CLEO result.
To study the remaining effects, the peaking background components have been divided according to the process from which they originate and have been ordered by decreasing level of importance: (i) the K ÿ and the electron originate from a D 0 decay (54%). The main source comes from D 0 ! K ÿ 0 e e . We have corrected the decay branching fraction used for this channel in the MC (2.02%) using recent measurements (2:17 0:16% [1]). The uncertainty on this value has been used to evaluate the corresponding systematic uncertainty. (ii) the electron comes from a converted photon or a Dalitz decay (24%). It has been assumed that the simulation correctly accounts for this component; (iii) the K ÿ does not originate from a charm hadron (14%). This happens usually when there is another negative charged kaon accompanying the D . We have studied the production of charged kaon accompanying a D using D 0 ! K ÿ events and measure a correction factor of 0:87 0:02 and 0:53 0:02, respectively, for same sign and opposite sign K-D pairs. The simulation is modified accordingly and the remaining systematic uncertainty from this source becomes negligible; (iv) fake kaon candidate (mainly pions) (6%) or fake electrons (1%). Differences between data and MC on the evaluation of fake rates have been studied in BABAR. As this affects small components of the total peaking background rate, the effect of these differences has been neglected.

Fitting procedure and radiative events
To fit form-factor parameters we compare the number of expected events in each bin with the measured one after all corrections. In this approach it is always assumed that the q 2 variation of f q 2 is given exactly by the form-factor parameterization. This hypothesis is not correct, a priori, for radiative decays as q 2 p D ÿ p K 2 p e p p 2 is not (perhaps) equal to the variable that enters in f for such decays. PHOTOS is used to generate decays with additional photons and the modified pole ansatz is taken to parameterize the hadronic form factor in signal events. To quantify possible distortion of the fit we compare the fitted value of a form-factor parameter with the one obtained from a fit to the generated q 2 distribution (see Table IV).
Corresponding corrections, given in the second column of Table IV, have been applied and quoted uncertainties enter in the systematic uncertainty evaluation.
To evaluate the importance of corrections induced by radiative effects, we have compared also the fitted value of a parameter on q 2 distributions generated with and without using PHOTOS. Measured differences are given in the last column of Table IV. They have been applied to the values quoted in Table III for the different parameters. We measure also that radiative effects affect mainly the fraction of the decay spectrum in the first bin in Table II which has to be increased by 0.0012 to correct for this effect.

Control of the statistical accuracy in the SVD approach
Once the number of SV is fixed, one must verify that the statistical precision obtained for each binned unfolded value is correct and if biases generated by removing information are under control. These studies are done with toy simulations. One observes that the uncertainty obtained from a fit of the unfolded distribution is underestimated by a factor which depends on the statistics of simulated events and is 1:06 in the present analysis. Pull distributions indicate also that the unfolded values, in each bin, have biases which are below 10% of the statistical uncertainty. Similar studies are done for the determination of form-factor parameters.

Summary of systematic errors
The systematic uncertainties for determining formfactor parameters are summarized in Table V. The systematic error matrix for the ten unfolded values is computed by considering, in turn, each source of uncertainty and by measuring the variation, i , of the corresponding unfolded value in each bin (i). The elements of this matrix are the sum, over all sources of systematic uncertainty, of the quantities i j . The total error matrix is evaluated as the sum of the matrices corresponding, respectively, to statistical and systematic uncertainties. Quoted uncertainties correspond to MC statistics. The last column gives the impact of the radiative effects on the form-factor measurements as predicted by PHOTOS.

Parameter
Measured difference (true-fitted) Bias from radiation

B. Comparison with expectations and with other measurements
The summary of the fits to the normalized q 2 distributions are presented in Table III. As long as we allow the form-factor parameters to be free in the fit, the fitted distributions agree well with the data and it is not possible to reject any of the parameterizations.
However, if the form-factor parameters are constrained to specific predicted values, the agreement is not good. For the ISGW2 model, the predicted dependence of the form factor on q 2 disagrees with the data (see Table III) and the fitted value of the parameter I differs from the predicted value, I 0:104 GeV ÿ2 by more than a factor two.
As observed by previous experiments, the simple pole model ansatz, with m pole m D s 2:112 GeV=c 2 does not reproduce the measurements. This means that the contribution from the continuum DK interaction cannot be neglected. If one introduces a second parameter pole to account for contributions from an effective pole at higher mass [see Eq. (10)] the two parameters are fully correlated and there is no unique solution, as illustrated in Fig. 8. The modified pole ansatz corresponds to pole 0.
In Table VI the fitted parameters for the simple pole ansatz and the modified pole [11] ansatz are compared for different experiments. The fitted pole masses are all well below the mass of the D s meson. The results presented here are consistent within the stated uncertainties with earlier measurements. Except for the BELLE measurement, all other measurements appear to favor a value of pole that is lower than the value predicted by lattice QCD, namely pole 0:50 0:04.
In Fig. 10, the dependence of the form factor on q 2 is presented. The data are compared to earlier measurements by the FOCUS experiment, as well as with predictions from lattice QCD calculations [19]. As stated above, the data favor a somewhat lower value for pole .
The data have also been mapped into the variable z. Figure 11 shows the product P f as a function of z. By convention, this quantity is constrained to unity at z z max , which corresponds to q 2 0. We perform a fit to a polynomial, P f 1 r 1 z r 2 z 2 . The data are compatible with a linear dependence, which is fully consistent with the modified pole ansatz for f q 2 , as illustrated in Fig. 11.  [5]. The band corresponds to lattice QCD [19] with the estimated uncertainty.

VI. BRANCHING FRACTION MEASUREMENT
The D 0 ! K ÿ e e branching fraction is measured relative to the reference decay channel, D 0 ! K ÿ . Events with a radiated photon are included in the signal sample. Specifically, we compare the ratio of rates for the decay chains D ! D 0 , D 0 ! K ÿ e e , and D 0 ! K ÿ in data and simulated events; this way, many systematic uncertainties cancel, (to simplify these expressions, apart in the first fraction, the notation for radiated s has been omitted.) The first line in this expression is the ratio of the branching fractions for the two channels used in the simulation: The second line is the ratio of the number of c c simulated events and the integrated luminosities for the two channels:  (17) The third line corresponds to the ratios of measured numbers of signal events in data and in simulation, and the last line gives the ratios of the efficiencies to data and simulation.

A. Selection of candidate signal events
The selection of D 0 ! K ÿ e e candidates is explained in Secs. IVA and IV B. For the rate measurement, the constraint on the D mass is not applied and also the momentum of the soft pion candidate is not included in the Fisher discriminant variable designed to suppress B B background. Since generic simulated signal events used in this measurement have been generated with the ISGW2 model, they have been weighted so that their q 2 distribution agrees with the measurement presented in this paper. Furthermore, we require for the Fisher discriminant F B B > 0 and restrict m < 0:16 GeV=c 2 . After background subtraction, there remain 76283 323 and 95302 309 events in data and simulation, respectively. This gives ND 0 ! K ÿ e e data ND 0 ! K ÿ e e MC 0:8004 0:0043: To select D 0 ! K ÿ candidates, the same data samples are used (in practice the overlap between the two data sets is higher than 95%) and particles, in each event, are selected in the same way. The same selection criteria on the Fisher discriminant to suppress B B-events, on the thrust axis direction and on other common variables are applied. Events are also analyzed in the same way, with two hemispheres defined from the thrust axis. In each hemisphere a D 0 candidate is reconstructed by combining a charged K with a pion of opposite sign. These tracks have to form a vertex and the K mass must be within the range 1:77; 1:95 GeV=c 2 . Another charged pion of appropriate sign is added to form a D candidate.
In addition, the following selection criteria are used: (i) the fraction of the beam momentum, in the c.m. frame, taken by the D candidate must exceed 0.48 to remove contributions from B B events; (ii) the measured D 0 mass must be in the range between 1.83 and 1:89 GeV=c 2 . This requirement eliminates possible contributions from remaining D 0 ! K ÿ K or ÿ decays [see Fig. 12(a)]; (iii) the vertex fits for the D 0 and D have to converge.
The m distribution for candidate events is shown in Fig. 12(c). The following components contribute to the D signal [see Fig. 12  The m distribution corresponding to other event categories does not show a peaking component in the D signal region. The total background level, is normalized to data using events in the m interval between 0.165 and 0:200 GeV=c 2 [see Fig. 12(c)]. This global scaling factor is equal to 1:069 0:011. After background subtraction, the m distributions obtained in data and simulation can be compared in Fig. 12(d). Since the D signal is narrower in the simulation, we use a mass window such that this difference has a small effect on the measured number of D events in data and in the simulation. There are 166 960 409 and 134 537 374 candidates selected in the interval m 2 0:142; 0:149 GeV=c 2 for simulated and data events, respectively. This means:

B. Efficiency corrections
The impact of the selection requirement on the reconstructed K mass has been studied. The K mass distri-bution signal for simulated events (not including radiative photons) is compared with the corresponding distribution obtained with data after background subtraction. The background contributions are taken from the simulation. The fraction of D 0 candidates in the selected mass range (between 1.83 and 1:89 GeV=c 2 ) is 97:64 0:25% in MC and 97:13 0:29% in data events. The ratio of efficiencies is equal to Since D 0 ! K ÿ e e events have been selected using a selection requirement on m, we need to confirm that the distribution of this variable is similar in data and simulation. This is checked by comparing the distributions obtained with D 0 ! K ÿ 0 events analyzed as if they were semileptonic decays. The m distributions are compared in Fig. 13. Below 0:16 GeV=c 2 , there are 0:93552 0:000 66 of the D candidates in the simulation and 0:93219 0:000 78 for data. The corresponding ratio of efficiencies (MC/data) is equal to 1:0036 0:0010. Using D 0 ! K ÿ 0 events, we also measure the difference between the fraction of events retained after the mass-constrained fits. Namely, it is 0:980 38 0:000 37 in the simulation compared to 0:97438 0:000 49 in data. The relative efficiency (MC/data) for this selection is 1:0062 0:0006. Based on these two measured corrections the ratio of efficiencies are: D 0 ! K ÿ e e MC D 0 ! K ÿ e e data 1:0098 0:0011: The quoted uncertainties, in this section, are of statistical origin and will be included in the statistical uncertainty on R D . Other differences between the two analyzed channels are considered in the following section and contribute to systematic uncertainties.

C. Systematic uncertainties on R D
A summary of the systematic uncertainties on R D are given in Table VII. They originate from selection criteria which are different for the two channels. Some of these uncertainties are the same as those already considered for the determination of the q 2 variation of f .

Correlated systematic uncertainties
Systematic uncertainties on the decay rate coming from effects that contribute in the measurement of the q 2 dependence of f q 2 are evaluated in Sec. VA and the full covariance matrix for the measurements of the number of D 0 ! K ÿ e e signal events and the fraction of the decay spectrum fitted in each of the ten bins is determined. Among the sources of systematic uncertainties, listed in Table V, those corresponding to: (i) the reconstruction algorithm, (ii) the tuning of the resolution on q 2 , (iii) the corrections applied on electron identification, and (iv) the background normalization are taken as common sources. Corresponding relative uncertainties on R D are given in Table VII.
Other systematic uncertainties contributing to the formfactor measurement also affect the reference channel and so their effects on R D cancel. They are related to the c-hadronization tuning and to the corrections applied on the kaon identification.

Selection requirement on the Fisher discriminant
The stability of the fraction of D 0 ! K ÿ e e events selected in data and in simulation as a function of the Fisher discriminant, F cc , designed to suppress c c background has been examined. This is done by comparing the distributions of this variable measured in data and in simulation as given in Fig. 4 for two selected intervals in m.
The value corresponding to F cc > 0 and for events selected in the range m < 0:16 GeV=c 2 is used as the central result and half the difference between the measurements corresponding to F cc greater than ÿ0:25 and 0:25 is taken as systematic uncertainty. This range corresponds to a relative change of 40% of the efficiency for signal events, and gives an uncertainty of 0:0061 on the ratio of data and simulated signal candidates given in Eq. (18).
3. D counting in D 0 ! K ÿ D candidates are selected in the range m 2 0:142; 0:149 GeV=c 2 . From the simulation it is expected that the fraction of signal events outside this interval is equal to 1.4%. Even though the D signal is slightly narrower in the simulation, there is not a large discrepancy in the tails. The fraction of signal events measured in the  sidebands m 2 0:140;0:142 0:149;0:150 GeV=c 2 is 0.4% and 0.5%, respectively, for simulation and data. An uncertainty of 0:004, corresponding to 30% uncertainty on the total fraction of events outside the selected m interval is assumed.

VII. SUMMARY
The decay rate distribution for the channel D 0 ! K ÿ e e has been measured in ten bins of q 2 in Table II. Several theoretical expectations for the variation of this form factor with q 2 have been considered and values for the corresponding parameters have been obtained (see Table III). The q 2 variation of the form factor can be parameterized with a single parameter using different expressions. The ISGW2 model with expected values for the parameters is excluded, as is the pole mass parameterization with m pole m D s . The value of the decay branching fraction has been also measured independent of a model in Section VI D.
Combining these measurements and integrating Eq. (2), the value of the hadronic form factor is obtained: where B is the measured D 0 ! K ÿ e e branching fraction, D 0 410:1 1:5 10 ÿ15 s [1] is the D 0 lifetime and I R q 2 max 0 jp K q 2 j 3 jf q 2 =f 0j 2 dq 2 . To account for the variation of the form factor within one bin, and, in particular, to extrapolate the result at q 2 0, the pole mass and the modified pole ansatze have been used; the corresponding values obtained for f 0 differ by 0.002.
Taking the average between these two values and including their difference in the systematic uncertainty, this gives f 0 0:727 0:007 0:005 0:007; (24) where the last quoted uncertainty corresponds to the accuracy on BD 0 ! K ÿ , D 0 and jV cs j. It agrees with expectations and, in particular, with LQCD computations [19]. Using the z expansion of Eq. (6) and assuming that m c 1:2 GeV=c 2 to evaluate V from Eq. (8), we find a 0 2:98 0:01 0:03 0:03 10 ÿ2 . The high accuracy of the present measurement will be a reference test for improved lattice determinations of the q 2 variation of f .