Branching fraction measurement of B̄0→D(*)+π- and B-→D(*)0π- and isospin analysis of B̄→D(*)π decays

Using 65 (cid:6) 10 6 (cid:1) (cid:1) 4 S (cid:3) ! B (cid:1) B events collected with the BABAR detector at the PEP-II e (cid:4) e (cid:5) storage ring at the Stanford Linear Accelerator Center, we measure the color-favored branching fractions B (cid:1) (cid:1) B 0 ! D (cid:4) (cid:2) (cid:5) (cid:3) (cid:7) (cid:1) 2 : 55 (cid:8) 0 : 05 (cid:8) 0 : 16 (cid:3) (cid:6) 10 (cid:5) 3 , B (cid:1) (cid:1) B 0 ! D (cid:2)(cid:4) (cid:2) (cid:5) (cid:3) (cid:7) (cid:1) 2 : 79 (cid:8) 0 : 08 (cid:8) 0 : 17 (cid:3) (cid:6) 10 (cid:5) 3 , B (cid:1) B (cid:5) ! D 0 (cid:2) (cid:5) (cid:3) (cid:7) (cid:1) 4 : 90 (cid:8) 0 : 07 (cid:8) 0 : 22 (cid:3) (cid:6) 10 (cid:5) 3 and B (cid:1) B (cid:5) ! D (cid:2) 0 (cid:2) (cid:5) (cid:3) (cid:7) (cid:1) 5 : 52 (cid:8) 0 : 17 (cid:8) 0 : 42 (cid:3) (cid:6) 10 (cid:5) 3 , where the ﬁrst error is statistical and the second is systematic. With these results and the current world average for the branching fraction for the color-suppressed decay (cid:1) B 0 ! D (cid:1)(cid:2)(cid:3) 0 (cid:2) 0 , the cosines of the strong phase difference (cid:3) between the I (cid:7) 1 = 2 and I (cid:7) 3 = 2 isospin amplitudes are determined to be cos (cid:3) (cid:7) 0 : 872 (cid:4) 0 : 008 (cid:4) 0 : 031 (cid:5) 0 : 007 (cid:5) 0 : 029 for the (cid:1) B ! D(cid:2) process and cos (cid:3) (cid:7) 0 : 924 (cid:4) 0 : 019 (cid:4) 0 : 063 (cid:5) 0 : 017 (cid:5) 0 : 054 for the (cid:1) B ! D (cid:2) (cid:2) process. Under the isospin symmetry, the results for cos (cid:3) suggest the presence of ﬁnal-state interactions in the D(cid:2) system.

B events collected with the BABAR detector at the PEP-II e e ÿ storage ring at the Stanford Linear Accelerator Center, we measure the color-favored branching fractions B B 0 ! D ÿ 2:55 0:05 0: 16 10 ÿ3 , B B 0 ! D ÿ 2:79 0:08 0:17 10 ÿ3 , BB ÿ ! D 0 ÿ 4:90 0:07 0:22 10 ÿ3 and BB ÿ ! D 0 ÿ 5:52 0:17 0:42 10 ÿ3 , where the first error is statistical and the second is systematic. With these results and the current world average for the branching fraction for the color-suppressed decay B 0 ! D 0 0 , the cosines of the strong phase difference between the I 1=2 and I 3=2 isospin amplitudes are determined to be cos 0:872 0:0080:031 ÿ0:007ÿ0:029 for the B ! D process and cos 0:924 0:0190:063 ÿ0:017ÿ0:054 for the B ! D process. Under the isospin symmetry, the results for cos suggest the presence of final-state interactions in the D system. DOI The B ! D and B ! D processes provide very good opportunities to test the theories of hadronic B-meson decays due to their clean and dominant hadronic decay channels. With the development of heavy quark effective theory (HQET) [1,2] and soft collinear effective theory (SCET) [3,4], the theoretical description for these hadronic decays has improved considerably, and the factorization hypothesis in heavy quark hadronic decay has been put on a more solid basis. The three decay amplitudes A for B ! D can be expressed in terms of two isospin amplitudes, A 1=2 and A 3=2 , under the isospin symmetry of the strong interaction: where isospin amplitudes A 1=2 and A 3=2 correspond to the transitions into D final states with pure I 1=2 and I 3=2 isospin eigenstates [5,6]. An identical decomposition holds for B ! D decays. The isospin amplitudes are not necessarily the same in the B ! D and B ! D systems. In the context of QCD factorization [6], A 1=2 and where m b is the b-quark mass and QCD is the QCD scale. The deviation of the ratio A 1=2 = 2 p A 3=2 from unity is a measure of the departure from the heavy-quark limit. The QCD factorization implies that the relative phase of A 1=2 and A 3=2 is O QCD =m b . Final-state interactions (FSI) in the I 3=2 and I 1=2 channels can lead to a nonzero . A large value of will substantially suppress the destructive interference for the color-suppressed decay B 0 ! D 0 0 , thereby increasing the associated branching fraction.
Recent experimental results on the color-suppressed decay B 0 ! D 0 0 [7-9] provide evidence for a sizable relative strong interaction phase between color-favored and color-suppressed B 0 ! D decay amplitudes. It has been suggested [5] that improved measurements of the color-favored hadronic two-body decay of the B meson will lead to a better understanding of these QCD effects. Further experimental results on the color-favored decay B ! D suggest the presence of final-state interactions in the B ! D process [10]. This paper presents new measurements of the branching fractions of B ÿ ! D 0 ÿ and B 0 ! D ÿ (charge conjugation is implied throughout this paper) and of the relative phase .
This analysis uses 65:2 0:7 10 6 B B pairs collected at the 4S resonance with the BABAR detector [11] at the PEP-II asymmetric-energy storage ring during the 2001-2002 data taking period. Charged tracks are detected by a 5-layer silicon vertex tracker and a 40-layer drift chamber. Hadrons are identified by measuring the ionization energy loss dE=dx in the tracking system and the opening angle of the Cherenkov radiation in a ring-imaging detector. Photons are identified by an electromagnetic calorimeter. These systems are mounted inside a 1.5-T solenoidal superconducting magnet.
Kaon and pion candidates are selected from chargedparticle tracks using dE=dx and the Cherenkov light signature. Each charged track, except the track used as the soft pion to reconstruct D ! D 0 , is required to have at least 12 hits in the drift chamber and a transverse momentum greater than 100 MeV=c. D 0 and D candidates are reconstructed in the K ÿ and K ÿ channels, respectively. In each case, D meson candidates are required to have a mass within 3 of the mean reconstructed mass value, where the mass resolution is approximately 7 MeV=c 2 for D 0 and 6 MeV=c 2 for D . A vertex fit is performed on D 0 (D ) candidates with the mass constrained to the nominal value [12]. A D 0 candidate is combined with a low momentum or 0 to form a D or D 0 candidate, where the 0 candidate is formed from two photon candidates and must have an invariant mass between 120 and 145 MeV=c 2 . Combinations with an invariant mass difference M m D 0 ÿ m D 0 between 143 and 148 MeV=c 2 for D and between 138 and 146 MeV=c 2 for D 0 , corresponding to 3 about the m peak, are retained. Each B meson candidate is reconstructed using the selected D or D candidate and an additional charged track that is not consistent with the kaon hypothesis.
To reject jetlike continuum background events, the normalized second Fox-Wolfram moment R 2 [13], computed with charged tracks and neutral clusters, is required to be less than 0.5. We also require j cos T j to be less than 0.85, where T is the angle between the thrust axis of the B candidate and the thrust axis of the rest of the event in the e e ÿ center-of-mass (CM) frame.
B candidates are identified using the beam-energy- where E and p are the energy and momentum of the reconstructed B candidate and s p is the total energy in the e e ÿ CM frame. B signal candidates have m ES m B , the B meson mass, and E ' 0, within their respective resolutions. The resolution in E, E , for various B modes ranges from 15.7 to 18.1 MeV. We require that jE ÿ hEij < 3 E . For events with more than one B candidate, a 2 is defined with the D mass m D , m and their resolutions as and the candidate with the smallest 2 is chosen. The event yield n for each mode of B ! D ÿ is extracted by fitting the m ES distribution of the selected B candidates with an unbinned extended maximum likelihood fit. The m ES distribution is fit to the sum of a signal component, modeled as a Gaussian, and a background shape. The background shape is parameterized as the sum of a Gaussian, representing the peaking background events that peak in m ES , and a phase space parameterization function [14] representing nonpeaking combinatorial background and continuum events. The parameters describing the background shape, including the relative normalization of the peaking component, are determined by fitting Monte Carlo (MC) simulated samples, with the signal events removed. The total signal and background event yields, as well as the shape parameters describing signal events, are free parameters in the fit. The fitted m ES distributions for each of the B meson decay modes are presented in Fig. 1. The peaking background yield n pb is about (2 -4)% of the observed B signal yield, as shown in Table I.
For each studied B decay mode of B ! D , the branching fraction is calculated as Here N B B is the total number of B B pairs; " is the efficiency determined from signal Monte Carlo events; f represents f ÿ or f 00 , the charged or neutral B meson production ratios at the 4S, which we assume to be f ÿ f 00 0:5; and BD is the branching fraction of D or D decaying to its reconstructed final state [12]. The branching fractions we obtain are reported in Table I.
The final states D selected by this analysis are, in general, accompanied by some small amount of final-state radiation (FSR). We model final-state radiation in our experiment with PHOTOS [15], which predicts that 6%-7% of our selected events, varying slightly with decay mode, are accompanied by an average FSR energy of about 17 MeV. Approximately two-thirds of this energy is produced in the initial B decay, while the remainder is generated in the D decay.
We summarize systematic uncertainties on the measurements from various sources in Table II. N B B is the uncertainty on the total number of B B pairs in data. The error on the efficiency, ", is due to signal Monte Carlo sample statistics. The uncertainty from combinatoric back- The uncertainties in tracking, particle identification, and 0 reconstruction efficiencies are due to potential residual inaccuracies in the Monte Carlo simulation, after correcting for known differences. The dominant uncertainty is from the D branching fractions BD and the tracking efficiency.
With the branching fractions of the four color-favored decay modes B 0 ! D ÿ and B ÿ ! D 0 ÿ , as well as the two color-suppressed modes B 0 ! D 0 0 , one can calculate cos. Following Ref. [16] (equations have been modified to use the notation from Ref. [5]), cos for B ! D (similarly for B ! D ) can be expressed as jA 3=2 j 2 1 3 ÿD 0 ÿ ; jA 1=2 j 2 ÿD ÿ ÿD 0 0 ÿ 1 3 ÿD 0 ÿ : (9) Using the measured branching fractions in this analysis, the ratio of the B lifetimes B ÿ = B 0 1:071 0:009 [12], and the branching fractions B To estimate the systematic error on cos for B ! D (and, similarly, B ! D ), we use a Monte Carlo technique [10]. We simulate 10 6 experiments, varying the measured branching fractions, the used color-suppressed decay branching fraction, and B ÿ = B 0 about their central values according to Gaussian distributions where their errors are taken as the sigmas of the Gaussian distributions, to calculate the cos. The correlation of the systematic errors between the two color-favored decay modes in the cos calculation is taken into account. We assume the errors are uncorrelated between the colorfavored and color-suppressed modes. The statistical error on cos is estimated in a similar fashion, with only the statistical errors on the branching fractions of colorfavored modes are used in the procedure. The resulting normalized distribution of cos, i.e., the estimated likelihood function of cos, is obtained. Figure 2 shows the likelihood function of cos from the described experiments in which both the statistical and systematic errors are taken into account.
We define 1 confidence interval of cos as the integral of its likelihood function over the region around the nominal value of cos, which is calculated from the central values of the branching fractions, to 68.27% (half below and half above the nominal value) of the total area.

Systematic error
for the B ! D and B ! D system, respectively, where the first error is statistical and the second is systematic. The likelihood function from the simulated experiments, with both statistical and systematic errors are taken into account, is shown in Fig. 3.
In summary, we have measured the branching fractions for the color-favored B 0 ! D ÿ and B ÿ ! D 0 ÿ decays. Using these measurements together with the current world averages for B B 0 ! D 0 0 and B B 0 ! D 0 0 , we extract the cosines of the relative strong phase in the D and D systems, and the ratios of the I 3=2 and I 1=2 isospin amplitudes. Our results for the B ! D branching fractions, except for B ÿ ! D 0 ÿ , are consistent with the current world average values [12] but have a better precision. The branching fraction of B ÿ ! D 0 ÿ from this measurement is greater than the world average by about 2. Our results for cos differ from unity by about 4:3 for B ! D decays and 1:1 for B ! D decays. The result of cos for B ! D decays is consistent with the result in Refs. [9,10], and under the isospin symmetry it suggests the presence of final-state interactions in B ! D decays.
We are grateful for the excellent luminosity and machine conditions provided by our PEP-II colleagues, and for the substantial dedicated effort from the computing organizations that support BABAR. The collaborating institutions wish to thank SLAC for its support and kind hospitality.