Measurement of the Absolute Branching Fractions for D − s → l − ν̄ l and Extraction of the Decay Constant f

The absolute branching fractions for the decays D (cid:1) s ! D (cid:2)(cid:1) s ! D (cid:1) s (cid:4) and X represents additional pions from fragmentation. The D (cid:1) s ! ‘ (cid:1) (cid:1) ‘ events are detected by full or partial reconstruction of the recoiling system DKX(cid:4)‘ . The branching fraction measurements are combined to determine the D (cid:1) s decay constant f D s ¼ ð 258 : 6 (cid:3) 6 : 4 (cid:3) 7 : 5 Þ MeV , where the ﬁrst uncertainty is statistical and the second is systematic.

s ! D À s and X represents additional pions from fragmentation. The D À s ! ' À ' events are detected by full or partial reconstruction of the recoiling system DKX'. The branching fraction measurements are combined to determine the D À s decay constant f D s ¼ ð258:6 AE 6:4 AE 7:5Þ MeV, where the first uncertainty is statistical and the second is systematic. DOI The D À s meson can decay purely leptonically via annihilation of the " c and s quarks into a W À boson [1]. In the standard model (SM), the leptonic partial width ÀðD À s ! ' À " ' Þ is given by where M D s and m ' are the D À s and lepton masses, respectively, G F is the Fermi coupling constant, and V cs is an element of the Cabibbo-Kobayashi-Maskawa quark mixing matrix. These decays provide a clean probe of the pseudoscalar meson decay constant f D s .
Within the SM, f D s has been predicted using several methods [2]; the most precise value by Follana [4] and leptoquarks [5], may explain this difference. In addition, f D s measurements provide a cross-check of QCD calculations which predict the impact of NP on B and B s meson decay rates and mixing. High precision determinations of f D s , both from experiment and theory, are necessary in order to discover or constrain effects of NP.
We present absolute measurements of the branching fractions of leptonic D À s decays with a method similar to the one used by the Belle Collaboration [6,7]. An inclusive sample of D À s 's is obtained by reconstructing the rest of the event in reactions of the kind e þ e À ! c " Here, D represents a charmed hadron (D 0 , D þ , D Ã , or Ã þ c ), K represents the K 0 S or K þ required to balance strangeness in the event, and X represents additional pions produced in the c " c fragmentation process. When the charmed hadron is a Ã þ c an additional antiproton is required to ensure baryon number conservation. No requirements are placed on the decay products of the D À s so that the selected events correspond to an inclusive sample. The 4-momentum of each D À s candidate, p r , is measured as the difference between the momenta of the colliding beam particles and the fully reconstructed DKX system: p r ¼ p e þ þ p e À À p D À p K À p X À p . The inclusive D À s yield is obtained from a binned fit to the distribution in the recoil mass m r ðDKXÞ ffiffiffiffiffi ffi p 2 r p . Within this inclusive sample, we determine the fraction of events corresponding to D À s ! À " , D À s ! e À " e , and D À s ! À " decays. In the SM, ratios of the branching fractions for these decays are e À " e : À " : À " ¼ 2 Â 10 À5 :1:10, due to helicity and phase-space suppression.
The analysis is based on a data sample of 521 fb À1 , which corresponds to about 677 Â 10 6 e þ e À ! c " c events, recorded near ffiffi ffi s p ¼ 10:58 GeV by the BABAR detector at the SLAC PEP-II asymmetric-energy collider. The detector is described in detail in Refs. [8,9]. Charged-particle momenta are measured with a 5 layer, double-sided silicon vertex tracker and a 40 layer drift chamber inside a 1.5 T superconducting solenoidal magnet. A calorimeter consisting of 6580 CsI(Tl) crystals (EMC) is used to measure electromagnetic energy. Measurements from a ringimaging Cherenkov radiation detector, and of specific ionization (dE=dx) in the silicon vertex tracker and drift chamber, provide particle identification (PID) of charged hadrons. Muons are mainly identified by the instrumented magnetic flux return, and electrons are identified using EMC and dE=dx information. The analysis uses Monte Carlo (MC) events generated with EVTGEN and JETSET [10,11] and passed through a detailed Geant4 [12] simulation of the detector response. Final state radiation from charged particles is modeled by PHOTOS [13]. Samples of MC events for e þ e À annihilation to q " q (q ¼ u, d, s, c, b) (generic MC) are used to develop methods to separate signal events from backgrounds. In addition, we use dedicated samples for D À s production and leptonic decays (signal MC) to determine reconstruction efficiencies and the distributions needed for the extraction of the signal decays.
We reconstruct D candidates using the following 15 modes: All 0 's and K 0 S 's used in this analysis are reconstructed from two photons or two oppositely charged pions, respectively, and are kinematically constrained to their nominal mass values [14]. The K 0 S in a D candidate must have a flight distance from the e þ e À interaction point (IP) greater than 10 times its uncertainty. For each D candidate we fit the tracks to a common vertex, and for each mode, we determine the mean and of the reconstructed signal mass distribution from a fit to data. We then simultaneously optimize a set of selection criteria to max- , where S refers to the number of D candidates after subtraction of the background B within a mass window defined about the signal peak, and where B is estimated from the sideband regions of the mass distribution. In addition to the size of the mass window, several other properties of the D candidate are used in the optimization: the center-of-mass (CM) momentum of the D, PID requirements on the tracks, the probability of the D vertex fit, and the minimum lab energy of 0 photons. The CM momentum must be at least 2:35 GeV=c in order to remove B meson backgrounds. After the optimization the relative contributions to the total signal sample are 74.0% D 0 , 22.6% D þ , and 3.4% Ã þ c . Multiple candidates per event are accepted.
To identify D mesons originating from D Ã decays we reconstruct the following decays: The photon energy in the laboratory frame is required to exceed 30 MeV for 0 ! and 250 MeV for D Ã0 ! D 0 decays. The invariant mass must be within 3 sigma of the 0 peak. For all D Ã decays, the mass difference mðD Ã Þ À mðDÞ is required to be within 2.5 sigma of the peak value.
A K candidate is selected from tracks not overlapping with the D candidate. PID requirements are applied to each K þ candidate, and a K 0 S candidate must have a flight distance greater than 5 times its uncertainty.
An X candidate is reconstructed from the remaining AE 's and 0 's not overlapping with the DK candidate. In the laboratory frame, a AE must have a momentum greater than 100 MeV=c and each photon from a 0 decay must have energy greater than 100 MeV. We reconstruct X modes without 0 's with up to three charged pions, and modes with one 0 with up to two charged pions. The total charge of the X candidate is not checked at this stage.
Finally, we select a candidate for the signal D ÃÀ s decay by requiring a minimum energy of 120 MeV in the laboratory frame, and an angle with respect to the direction of the D candidate momentum in the CM frame greater than 90 degrees. This photon cannot form a 0 or candidate when combined with any other photon in the event. In addition, the cluster must pass tight requirements on the shower shape in the EMC and a separation of at least 15 cm from the impact of any charged particle or the position of any other energy cluster in the EMC.
Only DKX candidates with a total charge of þ1 are selected to form a right-sign (RS) sample, from which we extract the D À s signal yield. The charm and strange quark content of the DKX must be consistent with recoiling from a D À s . The RS sample includes candidates for which consistency cannot be determined due to the presence of a K 0 S . We define a wrong-sign (WS) sample with the same charge requirement above, but by requiring that the charm and strange quark content of the DKX be consistent with a recoil from a D þ s . The WS sample contains a small fraction of signal events due mainly to DKX candidates for which the total charge is misreconstructed. The generic MC shows that the WS sample, after subtraction of the signal contribution, correctly models the backgrounds in the RS sample.
A kinematic fit to each DKX candidate is performed in which the particles are required to originate from a common point inside the interaction point region, and the D mass is constrained to the nominal value [14]. The 4-momentum of the signal D ÃÀ s is extracted as the missing 4-momentum in the event. We require that the D ÃÀ s candidate mass be within 2:5 of the signal peak. For MC signal events, the mean is found to be consistent with the nominal value and varies between 37 and 64 MeV=c 2 depending on the number of pions in X.
We perform a similar kinematic fit with the signal included and with the mass recoiling against the DKX constrained to the nominal D ÃÀ s mass [14] in order to determine the D À s 4-momentum. We require that the D À s CM momentum exceed 3:0 GeV=c, and that its mass be greater than 1:82 GeV=c 2 . After the final selections, there remain on average 1.7 D À s candidates per event, due mainly to multiple photons that can be associated with the D ÃÀ s decay. In order to properly count events in the fits described below, we assign weight 1=n to each D À s candidate, where n is the number of D À s candidates in the event. We define n R X and n T X to be the number of reconstructed and true pions in the X system, respectively. The efficiency for reconstructing signal events depends on n T X . However, the n T X distribution is expected to differ from the MC simulation due to inaccurate fragmentation functions used by JETSET. To correct for these inaccuracies, we extract the D À s signal yields from a fit to the two-dimensional histogram of m r ðDKXÞ versus n R X . The probability distribution function (PDF) for the signal distribution is written as a weighted sum of the MC distributions for j ¼ n T X , Sðm; n R X Þ ¼ X 6 j¼0 w j S j ðm; n R X Þ: The weights w j have to be extracted from this fit. To constrain the shape of the weights distribution, we introduce the parametrization w j / ðj À Þ e Àj together with the condition P j w j ¼ 1. This parametrization is motivated by the distribution of weights in the MC. The value ¼ À1:32 is taken from a fit to MC, whereas and are determined from the fit to data.
The RS and WS samples are fitted simultaneously to determine the background. The fit to the WS sample uses a signal component similar to that used in the RS fit, except that due to the small signal component, the weights are fixed to the MC values and the signal yield is determined from signal MC to be 11.8% of the RS signal yield. The shapes remaining after the signal component is removed from the WS sample, B i ðmÞ (i ¼ n R X ), are used to model the RS backgrounds. A shape correction is applied to B 0 to account for a difference observed in the MC. We add these components with free coefficients (b i ) to construct the total RS background shape: Bðm; n R X Þ ¼ P 3 i¼0 b i B i ðmÞði À n R X Þ. Thus in addition to , , and the total signal yield, there are 3 additional free parameters b i ði ¼ 0; 1; 2Þ in the RS fit. Having constructed the inclusive D À s sample, we proceed to the selection of D À s ! À " events within that sample. We use the m r ðDKXÞ range between 1.934 and 2:012 GeV=c 2 , which contains an inclusive D À s yield (N D s ) of ð67:2 AE 1:5Þ Â 10 3 . We require that there be exactly one more charged particle in the remainder of the event, and that it be identified as a À . In addition, we require that the extra neutral energy in the event, E extra , be less than 1.0 GeV; E extra is defined as the total energy of EMC clusters with individual energy greater than 30 MeV and not overlapping with the DKX candidate. Since the only missing particle in the event should be the neutrino we expect the distribution of E extra to peak at zero for signal events. We determine the 4-momentum of the " candidate through a kinematic fit similar to that described earlier in the determination of the D À s 4-momentum, but with the À included in the recoil system. In this fit we constrain the mass recoiling against the DKX system to the nominal value for the D À s [14]. To extract the signal yield, we perform a binned maximum likelihood fit to the m 2 r ðDKXÞ distribution using a signal PDF determined from reconstructed signal MC events that contain the signal decay chain D ÃÀ s ! D À s with D À s ! À " . The background PDF is determined from the reconstructed generic MC events with signal events removed. The fit is shown in Fig. 3(a), and the number of signal events extracted, N , is listed in Table I.
The D À s ! À " branching fraction is obtained from where the D À s ! À " reconstruction efficiency, " j , is determined using the signal MC sample with j ¼ n T X , and " j D s is the corresponding inclusive D À s reconstruction efficiency. The efficiency ratios " j =" j D s decrease from 87% to 33% as j increases from 0 to 6. The weighted average, " " , and the value determined for BðD À s ! À " Þ are listed in Table I. The statistical uncertainty includes contributions from N D s , " " , and N (with correlations taken into accounted). The systematic uncertainty is determined by varying the parameter values in the inclusive D À s fit which were fixed to MC values, by varying the resolution   on the D À s signal PDF (for both mass and n R X ), and by estimating how well the MC models the nonpeaking component of the signal PDF observed in Figs. 1 and 2. The nonpeaking signal component in the m r ðDKXÞ distribution arises from DKX candidates in events that contain the signal decay D ÃÀ s ! D À s , but for which the photon candidate is misidentified and is due to other sources such as 0 or decays, or tracks or K 0 L interacting in the calorimeter. Uncertainties are assigned for possible mismodeling of the signal or background m 2 r ðDKXÞ distributions due to possible differences in the position or resolution of the mass distribution, or mismodelings of different D À s decays. Uncertainties in the efficiencies due to tracking and À identification are included. This measurement supersedes our previous result [15].
Using a procedure similar to that for D À s ! À " we search for D À s ! e À " e events. The fit to the m 2 r ðDKXeÞ distribution, shown in Fig. 3(b), gives a signal yield N e consistent with 0. We obtain an upper limit on BðD À s ! e À " e Þ by integrating a likelihood function from 0 to the value of BðD À s ! e À " e Þ corresponding to 90% of the integral from 0 to infinity. The likelihood function consists of a Gaussian function written in terms of the variable BN D s " " e with mean and sigma set to N e and its total uncertainty, respectively. To account for the uncertainties on N D s " " e , the main Gaussian is convolved with another Gaussian function centered at the measured value of N D s " " e with sigma set to the N D s " " e total uncertainty. The value obtained for the upper limit is listed in Table I. We find D À s ! À " decays within the sample of inclusively reconstructed D À s events by requiring exactly one more track identified as an e À or À , from the decay À ! e À " e or À ! À " . We remove events associated with D À s ! À " decays by requiring m 2 r ðDKXÞ > 0:5 GeV2=c 4 . Since D À s ! À " events contain more than one neutrino we use E extra to extract the yield of signal events; these are expected to peak towards zero, while the backgrounds extend over a wide range. The signal and background PDFs are determined from reconstructed MC event samples. The fits are shown in Figs. 3(c) and 3(d); the signal yields are listed in Table I. We determine BðD À s ! À " Þ from the e À and À samples using Eq. (3) and accounting for the decay fractions of the À [14]. The values obtained are listed in Table I and are consistent with the previous BABAR result [16]. The error-weighted average [17] of the branching fractions is BðD À s ! À " Þ ¼ ð5:00 AE 0:35ðstatÞ AE 0:49ðsystÞÞ Â 10 À2 . The weights used in the average are computed from the total error matrix and account for correlations. As a test of lepton flavor universality we determine the ratio BðD À s ! À " Þ=BðD À s ! À " Þ ¼ ð8:27 AE 0:77ðstatÞ AE 0:85ðsystÞÞ, which is consistent with the SM value of 9.76.
As a cross-check of this analysis method, we measure the branching fraction for the hadronic decay D À s ! K À K þ À . Within the inclusive D À s sample, we require exactly three additional charged-particle tracks that do not overlap with the DKX candidate. PID requirements are applied to the kaon candidates. The mass of the K À K þ À system must be between 1.93 and 2:00 GeV=c 2 , and the CM momentum above 3:0 GeV=c. We combine the K À K þ À system with the signal and extract the signal yield from the mðKKÞ distribution. For this mode we choose the loose selection m r ðDKXÞ > 1:82 GeV=c 2 , because this variable is correlated with mðKKÞ; this corresponds to an inclusive D À s yield of N D s ¼ ð108:9 AE 2:4Þ Â 10 3 . We model the signal distribution using reconstructed MC events that contain the decay chain D ÃÀ s ! D À s and D À s ! K À K þ À . In the generic MC and a high statistics control data sample (for which the inclusive reconstruction was not applied), the background was found to be linear in mðKKÞ. From a fit to the mðKKÞ distribution, shown in Fig. 3(e), we determine a signal yield of N KK ¼ 1866 AE 40 events.
We compute the D À s ! K À K þ À branching fraction using Eq. (3). The efficiency for reconstructing signal events is determined from the signal MC in three regions of the K À K þ À Dalitz plot, corresponding to À , K À K Ã0 , and the rest. A variation of $8% is observed across the Dalitz plot, leading to a correction factor of 1.016 on " j KK . The weighted efficiency ratio is found to be " " KK ¼ 29:5%, and we obtain BðD À s ! K À K þ À Þ ¼ ð5:78 AE 0:20ðstatÞ AE 0:30ðsystÞÞ%. The first uncertainty accounts for the statistical uncertainties associated with the inclusive D À s sample and N KK . The second accounts for systematic uncertainties in the signal and background models, and the inclusive D À s sample, as well as the reconstruction and PID selection of the K À K þ À candidates. This result is consistent with the value ð5:50 AE 0:23 AE 0:16Þ% measured by CLEO-c [18].
Using the leptonic branching fractions measured above, we determine the D À s decay constant using Eq. (1) and the known values for m ' , m D s , jV ud j (we assume TABLE I. Average efficiency ratios, signal yields, branching fractions, and decay constants for the leptonic D À s decays. The first uncertainty is statistical and the second is systematic.