Study of B0→ρ+ρ- decays and constraints on the CKM angle α

We present results from an analysis of B0→ρ+ρ- decays using (383.6±4.2)×106 BB pairs collected by the BABAR detector at the PEP-II asymmetric-energy B factory at SLAC. The measurements of the B0→ρ+ρ- branching fraction, longitudinal polarization fraction fL, and the CP-violating parameters Slong and Clong are as follows: B(B0→ρ+ρ-)=(25.5±2.1(stat)-3.9+3.6(syst))×10-6, fL=0.992±0.024(stat)-0.013+0.026(syst), Slong=-0.17±0.20(stat)-0. 06+0.05(syst), Clong=0.01±0.15(stat)±0.06(syst). We determine the unitarity triangle angle α, using an isospin analysis of B→ρρ decays. One of the two solutions, α=[73.1,117.0]° at 68% confidence level, is compatible with standard model-based fits of existing data. Constraints on the unitarity triangle are also evaluated using an SU(3) symmetry-based approach. © 2007 The American Physical Society.


I. INTRODUCTION
Charge conjugation-parity (CP) violation was first seen in the decay of neutral kaons [1].It was shown some 40 years ago that CP violation is a necessary but insufficient condition required to generate a net baryon-antibaryon asymmetry in the Universe [2].The standard model (SM) of electroweak interactions describes CP violation as a consequence of a complex phase in the three-generation Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [3,4]: (1) V CKM describes the couplings of the u, c, and t quarks to d, s, and b quarks, which are mediated by the exchange of a W boson.In B-meson decays the CP violating parameters of the SM are most directly related to the angles and sides of the so-called unitarity triangle (UT), shown in Fig. 1.The angles , , and are defined as argÿV cd V cb =V td V tb ; (3) Any nontrivial phase in V ij is CP violating.CP violating phases originating from the CKM matrix are related to weak interactions, and therefore referred to as weak phases.In the Wolfenstein convention [5], the angle is the phase of V ub , is the phase of V td , and is the phase difference between V ub and V td constrained to satisfy ÿ ÿ through the unitarity of V CKM .In Fig. 1 the UT is shown in the complex ; plane, where the apex is given by the phase-convention indepen-dent definition [6]: The quest to understand CP violation remains, despite its observation in the B meson system by both the BABAR [7] and Belle experiments [8], since the SM does not, through the CKM phase, incorporate enough CP violation to explain the current matter-antimatter asymmetry [9,10].The CKM angle is measured to a 1 precision with b !c cs decays [11,12] and is consistent with current predictions [13][14][15].A significant deviation from results of SMbased fits of existing data for any of the UT angles would be a clear indication of new physics [16,17].
In principle, direct experimental measurements of can be made from decays that proceed mainly through a b !u ud tree diagram such as B 0 !ÿ , , ÿ , and a 1 [18,19].Interference between the direct decay and decay after B 0 B 0 mixing in each of these decays results in a time-dependent decay-rate asymmetry that is sensitive to the angle .Figure 2 shows the leading order tree and gluonic penguin loop contributions to the decay B 0 !ÿ .The presence of penguin contributions complicates the extraction of from these decays.Using isospin relations [20], measurements of the B ! 0 [21,22] and B 0 ! 0 0[23] branching fractions show that the penguin contribution in B ! is smaller than the leading tree diagram.The use of SU3 flavor symmetry to increase the precision on the weak phase constrained using B 0 !ÿ and B ! K 0 decays has also been proposed [24].Both of these approaches are discussed in Sec.X. Section X and Ref. [25] describe a number of possible sources of theoretical uncertainty.In B 0 !ÿ decays, a spin 0 particle (the B 0 meson) decays into two spin 1 particles ( mesons), as shown in Fig. 3. Subsequently, each meson decays into a 0 pair.As a result, the CP analysis of B 0 !ÿ is complicated by the presence of one amplitude with longitudinal polarization and two amplitudes with transverse polarization.The longitudinal amplitude is CP even, while the transverse amplitudes contain CP-even and CP-odd states.The decay is observed to be dominated by the longitudinal polarization [26,27], with a fraction f L defined as the fraction of the helicity zero state in the decay.Integrating over the angle between the decay planes , the angular decay rate is 1 ÿ f L sin 2 1 sin 2 2 ; (6) where the helicity angles i1;2 are the angles between the 0 momentum and the direction opposite to that of the B 0 in the rest frame.
We identify (tag) the initial flavor of the reconstructed B candidate, B rec , using information from the other B meson, B tag , in the event.The decay rate f f ÿ for a neutral B meson decaying into a CP eigenstate, where the B tag is a B 0 ( B 0 ), is given by Here t is the proper time difference between the decay of the B rec and B tag mesons in an event, B 0 is the mean B 0 lifetime, m d is the B 0 -B 0 mixing frequency [28], and f is the CP eigenvalue of the decay.This assumes that there is no difference between the B 0 and B 0 lifetimes, ÿ 0. The parameters S and C describe B-mixing induced and direct CP violation, respectively.For the longitudinal polarization ( f 1) S S long and C C long are defined as where CP q p A A [29], q and p are complex constants that relate the B meson flavor eigenstates to the mass eigenstates, and A=A is the ratio of amplitudes of the decay of a B 0 or B 0 to the final state under study.CP violation is probed by studying the time-dependent decay-rate asymmetry where ÿ ( ÿ) is the decay rate for B 0 ( B 0 ) tagged events.This asymmetry has the form The transverse polarization is an admixture of CP-even ( f 1) and CP-odd ( f ÿ1) parts, where each CP eigenstate has a decay-rate distribution of the form of Eq. ( 7).As the longitudinal polarization dominates this decay, we describe the CP admixture of the transverse polarization with common parameters S tran and C tran .We set S tran C tran 0, and vary these parameters when evaluating systematic uncertainties.
In the absence of penguin contributions in B 0 !ÿ , S long sin2 and C long 0. The presence of penguin contributions with different weak phases to the tree-level amplitude shifts the experimentally measurable parameter eff away from the value of .In the presence of penguin contributions eff , where and C long can be nonzero.
For B ! 0 one measures a time-integrated charge asymmetry, which is analogous to Eq. (10) and is defined as TABLE I. Measurements of the B ! decays.Branching fractions (B) are in units of 10 ÿ6 .The column marked C long (A CP ) denotes the time-dependent charge asymmetry C long in the case of the B 0 !ÿ decay, and the time-integrated charge asymmetry A CP N ÿ N N N ; (13) where N ( N) is the number of B ( B) decays to this final state.Table I summarizes the most recent measurements of the complete set of experimental inputs needed to extract from the B ! decays.
In this article, we present an update of the timedependent analysis of B 0 !ÿ and measurement of the CKM angle reported in Ref. [26] and the branching fraction reported in Ref. [30].Improvements to the analysis reported here include an increase in the data sample analyzed, a tighter constraint on the proper time difference between the two B meson decays in selected events, an improved algorithm to determine the flavor of B mesons, a modified multivariate analyzer for background suppression, and an improved description of the signal and background model.

II. THE DATA SET AND BABAR DETECTOR
The results presented in this paper are based on data collected with the BABAR detector at the PEP-II asymmetric e e ÿ storage ring [32] operating at the Stanford Linear Accelerator Center.At PEP-II, 9.0 GeV electrons and 3.1 GeV positrons are collided at a center-of-mass energy of 10.58 GeV which corresponds to the mass of the 4S resonance.The asymmetric energies result in a boost from the laboratory to the center-of-mass (CM) frame of 0:56.The data set analyzed has an integrated luminosity of 349 fb ÿ1 corresponding to 383:6 4:2 10 6 B B pairs recorded at the 4S resonance (on peak).An additional 27:2 fb ÿ1 of data were recorded about 40 MeV below this energy (off peak) for the study of continuum background, where light quarks are produced in the process e e ÿ !q q (q u, d, s, c).
The BABAR detector is described in detail elsewhere [33].Surrounding the interaction point is a five doublesided layer silicon vertex tracker (SVT) which measures the impact parameters of charged particle tracks in both the plane transverse to and along the beam direction.A 40layer drift chamber (DCH) surrounds the SVT and provides measurements of the momenta for charged particles.Both the SVT and DCH are surrounded by a solenoid magnet that provides a 1.5 T magnetic field.Charged hadron identification is achieved through measurements of particle energy loss in the tracking system and the Cherenkov angle obtained from a detector of internally reflected Cherenkov light.A CsI(Tl) electromagnetic calorimeter (EMC) provides photon detection, electron identification, and 0 reconstruction.Finally, the instrumented flux return of the magnet allows discrimination of muons from pions.For the most recent 134:0 fb ÿ1 of data, a portion of the resistive plate chambers constituting the muon system has been replaced by limited streamer tubes [34 -36].
We use a right-handed coordinate system with the z axis along the electron beam direction and the y axis upward, with the origin at the nominal beam interaction point.Unless otherwise stated, kinematic quantities are calculated in the laboratory rest frame.The other reference frame which we commonly use is the CM frame of the colliding electrons and positrons.
We use Monte Carlo (MC) simulated events generated using the GEANT4 [37] based BABAR simulation.

III. RECONSTRUCTION OF B CANDIDATES
A. Photon and 0 reconstruction Photons are reconstructed from localized energy deposits in the EMC that are not associated with a charged track.We require photon candidates to have an energy greater than 50 MeV, and a lateral shower profile [38] to be consistent with the photon hypothesis.We reconstruct neutral pions from combinations of two distinct photon candidates where the invariant mass m satisfies 0:10 < m < 0:16 GeV=c 2 .

B. reconstruction
We combine reconstructed 0 mesons with charged tracks that are consistent with the hypothesis to form candidates.The invariant mass m 0 of the reconstructed is required to lie between 0.5 and 1:0 GeV=c 2 , to select events in the vicinity of the resonance.We require that the helicity angle of each meson satisfies ÿ0:90 < cos i < 0:98.This selection criteria suppress continuum and B backgrounds.

C. B 0 reconstruction and event selection
We combine two oppositely charged candidates to form the B candidate B rec .All particles in the rest of the event (ROE) are combined to form the other B meson in the event B tag .In addition to the two charged tracks in the B rec , we require that there is at least one track in the B tag .
In order to suppress potential backgrounds from e e ÿ ! e e ÿ , ÿ events, we require the second-tozeroth Fox-Wolfram moment R 2 [39] of the event to be less than 0.98.Continuum events are the dominant background which is reduced by requiring the absolute value of the cosine of the angle between the B rec thrust axis (TB) and that of the ROE (TR) to satisfy j cosTB; TRj < 0:8.We retain 17.1% and 20.1% of longitudinal and transverse signal, respectively, by requiring the aforementioned selection criteria.
We calculate t z=c from the measured separation z between the B rec and B tag vertices [40].We determine the B rec vertex from the two charged-pion tracks in its decay.The B tag decay vertex is obtained by fitting the other tracks in the event with constraints from the B rec momentum and the beam-spot location.The RMS resolution on t is 1.1 ps.We only use events that satisfy jtj < 15 ps and for which the error on t ( t ) is less than 2.5 ps.
Signal events are identified using two kinematic variables, the difference E between the CM energy of the B rec , E B , and and the beam-energy-substituted mass, where s p is the total CM energy.The B rec momentum p B and four-momentum of the initial state E i ; p i are defined in the laboratory frame.We accept candidates that satisfy 5:25 <m ES < 5:29 GeV=c 2 and ÿ0:12 < E < 0:15 GeV.An asymmetric E selection is used in order to reduce backgrounds from higher-multiplicity B decays.The resolution on m ES is dominated by the knowledge of the energy of the e and e ÿ beams, and the resolution on E is dominated by the reconstruction performance of the EMC.
After the selection criteria mentioned above have been applied, the average number of candidates per event is approximately 1.5.In events with more than one candidate, we select the candidate that minimizes the sum of m ÿ m 0 2 where m 0 is the true 0 mass [28].In 0.3% of events, the same 0 mesons are used by multiple B candidates.In such an event we randomly select the candidate to keep.

IV. CONTINUUM BACKGROUND SUPPRESSION
In addition to the cuts on cos i , R 2 , and j cosTB; TRj that directly remove background events, we use an artificial neural network in order to discriminate between continuum background and signal events.For this purpose we combine the following eight variables into a single output, N .
(i) The coefficients, L 0 , L 2 , split into sums over the ROE for neutral and charged particles: L 0;n , L 2;n and L 0;c , L 2;c .The coefficients are defined as L k P p j j cos j j k , where k 0, 2, p j is the particle momentum, and j is the angle of the particle direction relative to the thrust axis of the B candidate.Both p j and j are defined in the CM frame.(ii) j cosB; Zj, the absolute value of the cosine of the angle between the direction of the B and z axis in the CM frame.This variable follows a sine squared distribution for B B events, whereas it is almost uniform for q q. (iii) j cosTB; TRj.This variable, previously defined in Sec.III C, is strongly peaked at unity for q q events.B B events are more isotropic as the B mesons are produced close to the kinematic threshold.(iv) j cosTB; Zj, the absolute value of the cosine of the angle between the B thrust and the z axis.(v) The scalar sum of the transverse momenta p t in the ROE.This sum includes neutral and charged particles.The distributions of these input variables are shown in Fig. 4. Figure 5 shows N for signal MC simulated events and continuum background samples (off-peak data) and the efficiency for signal and q q background as a function of cut on the minimum value of N .We require N to be greater than ÿ0:4.Note that later we use this variable in the maximum likelihood fit described in Sec.VII.
The samples used to train the neural net were correctly reconstructed MC simulated events and off-peak data.To avoid overtraining, we used an independent sample of these data (i.e.distinct from the sample used for the training) to check the performance of the network.The training is stopped when the separation between the signal and continuum background is stable.

V. MISRECONSTRUCTED SIGNAL AND SELECTION EFFICIENCY
Misreconstructed signal candidates, referred to as selfcross-feed (SCF) signal, may pass the selection requirements even if one or more of the pions assigned to the ÿ state belongs to the other B in the event.These SCF candidates constitute 50.7% (27.9%) of the accepted longitudinally (transversely) polarized signal.The majority of SCF events have both charged pions from the ÿ final state and unbiased CP information.These correct (right) track SCF events are denoted by RT SCF.There is a SCF component (13.8% of the signal) where at least one track in B rec is from the ROE.These wrong track (WT) events are used to determine the signal yield N signal and f L but have biased CP information, and are treated separately for the CP result.A systematic error is assigned to the CP results from this type of signal event.The total selection efficiency for longitudinally (transversely) polarized signal is 7.7% (10.5%).

VI. B MESON FLAVOR TAGGING
A key ingredient in the measurement of time-dependent CP asymmetries is to determine whether, at the time of decay of the B tag , corresponding to t 0, the B rec was a B 0 or a B 0 .This ''flavor tagging'' is achieved with the analysis of the decay products of the recoiling B meson B tag .The overwhelming majority of B mesons decay to a final state that is flavor specific, i.e.only accessible from either a B 0 or a B 0 , but not from both.The purpose of the flavor tagging algorithm is to determine the flavor of B tag with the highest possible efficiency tag and lowest possible probability w of assigning a wrong flavor to B tag .It is not necessary to fully reconstruct B tag in order to determine its flavor.In the presence of a finite mistag probability w, the CP asymmetry is reduced by a dilution factor 1 ÿ 2w, so Eq. ( 7) becomes The figure of merit for the performance of the tagging algorithm is the effective tagging efficiency which is related to the statistical uncertainty in the coefficients S and C through We use a neural network based technique [11,40] that isolates primary leptons, kaons, and pions from B decays to final states containing D mesons, and high momentum charged particles from B decays, to determine the flavor of the B tag .The output of this algorithm is divided into seven mutually exclusive categories.These are (in order of decreasing signal purity) Lepton, Kaon I, Kaon II, Kaon-Pion, Pion, Other, and Untagged.The performance of this algorithm is determined using fully reconstructed neutral B decays to D ÿ ; ; a 1 and is summarized in Table II.The categories assigned correspond to events with leptons, kaons, and pions in the decay products of B tag .The untagged events contain no flavor information and therefore carry no weight in the time-dependent analysis.In addition, these events are not considered useful for the branching fraction measurement since they are dominated by continuum background.The total Q of this algorithm is 32:7 0:7%.

VII. LIKELIHOOD FIT METHOD
On applying the selection criteria described above, we obtain a sample of 37 424 events that enter the fit.These events are dominated by backgrounds from q q (81.4%) and B B (16.6%) events.The remaining 2% of events are considered as signal.We distinguish between the following components in the fit: (i) correctly reconstructed signal, (ii) SCF signal, split into RT and WT parts, (iii The dominant B backgrounds come from components (iii) and (iv).The signal, continuum, and B background models are described in Secs.VII A, VII B, and VII C, respectively.We use an unbinned, extended maximum likelihood (ML) fit to extract N signal , f L , S long , and C long .The likelihood function incorporates the following previously defined eight discriminating variables to distinguish signal from the backgrounds: m ES , E, t, N , and the m k and cos k values of the two mesons, where k 1, 2. For each of the aforementioned components j, we construct a probability density function (PDF) that is the product of one-dimensional PDFs for each of the variables.The PDFs do not account for all possible correlations among the discriminating variables, and we account for possible biases as a result of neglecting these correlations as discussed in Sec.VIII.For each event i, the PDFs can be written as We determine the PDFs for signal and each of the B background components from MC distributions for each discriminating variable.The PDFs for the continuum background are determined from on-peak and off-peak data.The likelihood function is where n j are the event yields of hypotheses j (determined from the ML fit) and N is the observed number of events in the sample.The normalized exponent takes into account Poisson fluctuations in the total number of events.We simultaneously fit events in the mutually exclusive flavor tag categories listed in Table II, excluding events in the untagged category.There are six continuum background yields, one for each flavor tag category, and we use a single yield for the signal and each of the B backgrounds, while accounting for the relative fractions of events expected in each of the flavor tag categories.The fit results are obtained by minimizing the value of ÿ lnL.

A. Signal model
The signal has contributions from true and SCF longitudinally (long) and transversely (tran) polarized events.In addition to this, the longitudinally polarized SCF signal PDF is further subdivided into the two categories of misreconstructed signal: RT and WT SCF signal, and all signal PDFs are subdivided according to the B tag flavor tag category.The total signal PDF is given by where f long RTWT SCF is the fraction of longitudinally polarized RT (WT) SCF signal and f tran SCF is the fraction of transversely polarized SCF signal.The PDFs P long RTWT SCF and P tran SCF are defined accordingly.In order to extract the observable f L from the fit, we account for the different reconstruction efficiencies for longitudinally and transversely polarized signal using TABLE II.Tagging efficiency , average mistag fraction w, mistag fraction difference between B 0 and B 0 tagged events w, and Q for B 0 !ÿ events.where f obs L is the fraction of signal from longitudinally polarized events in the data, and R is the ratio of longitudinal to transverse polarized signal efficiencies.
All of the SCF PDFs are further subdivided into parts that contain correctly reconstructed and misreconstructed mesons.Four subdivisions are used, corresponding to both mesons being correctly reconstructed (TT), one being correctly reconstructed and the other being misreconstructed (TF or FT), and both mesons being misreconstructed (FF).Here the first (second) meson is positively (negatively) charged.The SCF PDFs have the form where T and F are correctly and misreconstructed mesons, respectively, f TT , f TF , f FT , f FF are the fractions of each type of misreconstructed event, and the PDFs for m ES , E, and the helicity and mass for each reconstructed meson are the products of one-dimensional PDFs.The signal PDFs for N and t have different shapes according to the B tag flavor tag category assigned to an event.
The longitudinally polarized true and RT SCF signal m ES distributions are described by a Gaussian with a power law tail [41][42][43] which takes the following form: where m 0 is the mean, is the width of the Gaussian part, and the parameters and n describe the exponential tail.The remaining PDFs for the signal m ES distributions are the sum of a Gaussian and a Gaussian with an exponential tail.The signal E distributions for the longitudinally polarized SCF events are described by first and second order polynomials, with the exception of TT events, which are modeled with the sum of a second order polynomial and a Gaussian.The other E PDFs are described by the sum of a Gaussian and a Gaussian with an exponential tail.
The N distributions are modeled using smoothed histograms of MC simulated events.The signal t distribution is described by Eq. ( 16) convolved with a triple Gaussian resolution function given by R sig t; t f core Gt; core t ; core t f tail Gt; tail t ; tail t f outlier Gt; outlier ; outlier ; (25) where G is a Gaussian with mean l and width l for l core, tail, and outlier.The parameters tail , outlier , and outlier are set to 3.0 ps, 8.0 ps, and 0.0 ps, respectively.The remaining parameters of the resolution function are determined from signal MC simulated events scaled by the differences between data and MC simulated events for large samples of exclusively reconstructed B decays to D ÿ ; ; a 1 final states described in Sec.VI.The values of the mean and width of the core Gaussian contribution to the resolution function depend on the flavor tagging category of an event.There are three signal t PDFs, one for the true and RT SCF longitudinally polarized signal, one for the WT SCF longitudinally polarized signal, and one for the transversely polarized signal.Each of these PDFs has distinct values of S and C as described in Secs.I and V.The WT SCF longitudinally polarized signal has a different resolution function with respect to the rest of the signal.The cos i distribution for true signal events is described by the expected cosine squared or sine squared distribution multiplied by a polynomial acceptance function.The cos i PDFs for T mesons in longitudinally (transversely) polarized SCF signal are described by sixth order polynomials (smoothed histograms of MC simulated events).The cos i PDFs for misreconstructed mesons are described by smoothed histograms of MC simulated events.The longitudinally polarized signal mass PDFs are described using relativistic Breit-Wigner function and third order polynomials for T and F contributions, respectively.The transversely polarized signal SCF mass PDF distributions are described using smoothed histograms of MC simulated events.
In addition to N signal , f L , S long , and C long , we determine the mean of the core Gaussian part of the m ES and E distributions from the fit.We obtain parameters consistent with expectations from MC simulated events.

B. Continuum background model
The continuum background PDFs are also subdivided according to the misreconstruction of mesons in the final state and have the form P P m ES P EP N P t f TT P TT cos 1 ; m 1 ; cos 2 ; m 2 f TF P TF cos 1 ; m 1 ; cos 2 ; m 2 f FT P FT cos 1 ; m 1 ; cos 2 ; m 2 f FF P FF cos 1 ; m 1 ; cos 2 ; m 2 ; where helicity and mass PDFs for each reconstructed meson are the products of one-dimensional PDFs.
The continuum distribution for m ES is described by a phase-space-motivated distribution [44] with the following form: where .The E and N shapes are modeled with third and fourth order polynomials, respectively.The parameters of the m ES , E, and N shapes are allowed to vary in the fit to the on-peak data.The continuum mass distribution is described using a relativistic Breit-Wigner function and a third order polynomial PDF for T and F contributions, respectively.The cos i distribution is described by a third order polynomial.The continuum t distribution has a prompt lifetime component convolved with a triple Gaussian resolution function.The parameters of the mass and helicity distributions are obtained from a fit to the off-peak data, and the remaining parameters are determined in the fit.

Charm B backgrounds
Combinatorial events from b !c transitions are the dominant B backgrounds.These components have shapes similar to continuum and do not peak in the signal region for the discriminating variables.The functional form used for the PDFs of these background components is given by Eq. ( 26).We parametrize the m ES and E distributions of these backgrounds using the phase-space-motivated distribution of Eq. ( 27), and a third order polynomial, respectively.The mass distribution is described using a relativistic Breit-Wigner function and a third order polynomial PDF for T and F contributions, respectively.The remaining PDFs are described using smoothed histograms of MC simulated events.Each of the PDF parameters for these backgrounds are determined from samples of MC simulated events, and the yields of these components are determined in the fit.When studying systematic uncertainties, we modify the PDF used for t so that it has a form similar to the signal which uses the resolution function of Eq. ( 25).An effective lifetime that is smaller than B 0 is used to account for misreconstruction of these events and the finite charm meson lifetime.The value of this parameter is obtained by fitting MC simulated events for this category of events.

Charmless B backgrounds
Some of the charmless B backgrounds have PDFs similar to the signal for one or more of the discriminating variables, so it is important to correctly account for such events in the fit.We consider the 20 exclusive and two inclusive components of this type of B background listed in Table III.If a charmless B background decay contributes an event yield more than 1% of the expected signal yield, we model that mode exclusively.The remaining 140 charmless background components considered were combined with the appropriate weightings to form the neutral and charged inclusive charmless B background components.Where possible, we use branching fractions from existing measurements.Where measurements are not available, we have either tried to use SU2 and SU3 flavor symmetries to relate an unmeasured decay to a measured one or, where this is not possible, we have assigned a branching fraction of 10 10 ÿ6 to a decay mode.An uncertainty of 100% is assumed on all extrapolated branching fractions.We assume that the a 1 meson decays into a three pion final state via .The decay B 0 ! a 0 1 0 is penguin dominated and is expected to have a small branching fraction, as is the case for the penguin dominated decays B 0 ! 0 0and B 0 ! 0 0.As a result, we have neglected possible contributions from this potential source of background.Although we do not see evidence for B 0 ! 0, we do consider this channel as a potential source of background.Contribution from B 0 ! 4 is expected to be even smaller and is neglected.When considering the systematic uncertainty arising from interference between signal and other ÿ 0 0 final states, we assume that the nonresonant ÿ 0 0 final state has the same branching fraction as that of B 0 ! 0, where BB 0 ! 0is calculated from the yield determined in the B 0 !ÿ nominal fit.This is discussed further in the Appendix.
The functional form used for the PDFs of these charmless B background components is given by Eq. ( 26).The E distributions are described by third order polynomials, except for nonresonant B 0 ! 0which uses smoothed histograms of MC simulated events.The m distributions for true mesons are parametrized using a relativistic Breit-Wigner function, and the fake m (combinatorial 0 ) distribution is described using third order polynomials.The remaining background shapes are described using smoothed histograms of MC simulated events.The yield for B 0 ! 0decays is allowed to vary in the fit.All other charmless background yields are fixed to expectations.This constraint is relaxed when studying possible sources of systematic uncertainties.When studying systematic uncertainties from possible CP violation in the B background, we modify the PDF used for t so that it has a form similar to the signal which uses the signal resolution function in Eq. ( 25).

VIII. LIKELIHOOD FIT VALIDATION
Before applying the fitting procedure to the data, we subject it to various tests.The aim of these tests is to verify that one can extract the signal observables N signal , f L , S long , and C long in a controlled way.Consistency of the likelihood fit is checked with ensembles of experiments simulated from the PDFs.The event yields are generated according to a Poisson distribution with mean n j for each category in these ensembles.In each of these tests we verify that the values generated for the signal observables are reproduced with the expected resolution.The distribution of ÿ lnL for the ensemble of experiments in comparison to that obtained when fitting the data provides an additional, but limited, cross-check of the consistency.
The PDFs used in the likelihood do not account for all possible correlations among the discriminating variables.We account for possible biases as a result of neglecting these correlations by fitting ensembles of experiments obtained from samples of signal and the charmless B background MC simulated events combined with charm backgrounds and q q background events simulated from the PDFs.The MC simulated events used in these ensembles have these correlations modeled.We find a positive bias of 38:5 5:3 events on N signal , and a negative bias of 0:016 0:002 on f L , and we do not observe a significant bias on S long and C long .
As continuum events are the dominant background, we apply the fitting procedure to the off-peak data (after correcting for the difference in s p ) to confirm that we do not find a fake signal in this control sample of events.We fit 790 off-peak events and extract signal and continuum yields of 8 7 and 782 28 events, respectively.
A blind analysis technique has been used for the extraction of S long and C long , where the actual values of these observables have been hidden by an offset.With the values of S long and C long hidden, we perform the following fit cross-checks.We first verify that the uncertainties on the signal observables, and the value of ÿ lnL obtained from the fit to data, are compatible with the ensembles of simulated experiments described above.We then validate the stability of our results by introducing a variety of modifications to the fitting procedure.In addition to the nominal set of variables determined from data, we allow the physics parameters B 0 and m d to vary in turn.The signal observables do not change significantly when doing this, and the results obtained for B 0 and m d are 1:72 0:16 ps and 0:36 0:22@=ps, respectively, consistent with the reported world average [28].
The Lepton and Kaon I tagged events have the highest signal purity, and lowest w probability, and dominate our CP results.We check that the results obtained from fitting only these categories of events are compatible with our nominal blind results.When doing this, we observe shifts of ÿ0:08 ( ÿ 0:05) on S long (C long ) relative to the nominal result.
Given that there are a number of B backgrounds that contribute to the data that have yields fixed to expectations, we validate this assumption by allowing each fixed B background yield to vary in the fit.We obtain background yields consistent with our expectations, and observe the shifts on signal parameters listed in Table IV.A systematic uncertainty corresponding to the largest observed deviation is assigned on our results.Once these checks have been completed, the fitting procedure is frozen and we extract the values of S long and C long by removing the offset.

A. Fit results
From the ML fit described above, we obtain the following results: Nsignal 729 60stat; f L 0:992 0:024stat; S long ÿ0:17 0:20stat; C long 0:01 0:15stat; after correction for a 38:5 5:3 event fit bias (see Sec. VIII), a ÿ76 event bias from SCF on the signal yield (see Sec. IX B), and a correction for a ÿ0:016 0:002 fit bias on f L .The correlation between S long and C long is ÿ0:035.We discuss the origin of these fit biases in Sec.IX B. We calculate the branching fraction of this decay using B N signal =N b b , where is the efficiency for signal corresponding to the observed f L , and N B pairs is the number of B B pairs analyzed.We obtain B B 0 !ÿ 25:5 2:1stat 10 ÿ6 : The B 0 ! 0background yield obtained from the fit is 9:2 53:6 events.Figure 6 shows distributions of m ES , E, cos i , and m 0 for the Lepton and Kaon I tagged events with a loose requirement on N .Relative to the total number of events in the fit, the plot of m ES contains 15.6% of the signal and 1.1% of the total background.For the other plots there is an added constraint that m ES > 5:27 GeV=c 2 ; these requirements retain 13.9% of the signal and 0.4% of the total background.Figure 7 shows the t distribution for B 0 and B 0 tagged events, as well as the time-dependent decay-rate asymmetry of Eq. ( 10).Here we apply the same selection criteria as in Figs.6(b)-6(d).

B. Systematic uncertainty studies
Table V lists the possible sources of systematic uncertainties on the values of N signal , f L , S long , and C long that have been studied, and are described in the following.
(i) The uncertainty from PDF parametrization is obtained by varying PDF shape parameters by 1, in turn.The deviations obtained are added in quadrature to give the quoted uncertainty.(ii) The systematic uncertainty from the fraction of SCF events is obtained from the difference between the nominal fit result and that obtained when fitting the data and removing the SCF from the fit model.In the case of the signal yield we correct for half of the difference observed, and use 100% of the correction as an uncertainty.The uncertainty on the other signal observables comes from the difference observed between the two sets of fit results.(iii) The kinematic endpoint position in m ES is extracted from the fit.Changes in beam energy in the data can affect the endpoint position.To account for possible effects of this, we vary the kinematic endpoint position in m ES by 0:45 MeV=c 2 which is determined from control samples of fully reconstructed B 0 decays.(xii) We estimate the systematic error on our results coming from neglecting the interference between B 0 !ÿ and other 4 final states: B ! a 1 , 0 , and 0 0 .Details of this calculation are given in the Appendix.(xiii) As the PDFs used in the ML fit do not account for all of the correlations between discriminating variables used in the fit, the results have a small bias.We calculate the fit bias on the signal observables as described in Sec.VIII.The observed bias on N signal and f L is corrected, and 100% of the correction is assigned as a systematic uncertainty.We do not observe a significant bias on S long and C long , and we conservatively assign a systematic uncertainty from fit bias based on the statistical precision of this test.(xiv) Small imperfections in the knowledge of the geometry of the SVT over time can affect the measurement of S long and C long .We vary the alignment according to the results obtained from the study of e e ÿ ! e e ÿ , ÿ events in order to estimate the magnitude of this systematic error on our CP results.The branching fraction has multiplicative systematic uncertainties from the reconstruction of 0 mesons in the detector (6%), uncertainties in the reconstruction of charged particles (0.8%), and the discrimination of from other types of charged particles (1%).In addition to these uncertainties, there is a 1.1% uncertainty on the number of B B pairs in the data sample.The statistical uncertainties arising from the MC samples used in this analysis are negligible.

A. The SU2 isospin analysis
In SU2 isospin symmetry, the amplitudes of B decays to final states [20] are related by 1 2 p A ÿ A 0 ÿ A 00 ; (28) for the longitudinal polarization and each of the CP eigenstates of the transverse polarization, where A ij ( A ij ) are the amplitudes of B ( B) decays to the final state with charge ij, where i ; ÿ; 0 and j ÿ; 0. These two relations correspond to triangles in a complex plane as shown in Fig. 8.In the usual phase convention [20], the amplitudes A are rotated to Ã in order to align the base of the triangles.After aligning A 0 and A ÿ0 , the phase difference between A ÿ and A ÿ is 2.There are two sources of ambiguities on the measurement of .There is a fourfold ambiguity coming from the orientation of the isospin triangles, and a twofold ambiguity from the trigonometric conversion in Eq. (12).In order to measure , one must measure the branching fractions and charge asymmetries of B decays to ÿ , 0 , 0 0 .There are several assumptions implicitly used in the isospin-based direct measurements of : (i) This approach only considers tree and gluonic penguin contributions.Possible contributions from electroweak penguins (EWP) are neglected as they do not obey SU2 isospin symmetry.These EWPs have the same topology as the gluonic penguin diagram in Fig. 2, with the gluon replaced by or Z 0 bosons.In the absence of EWP contributions jA 0 j j A ÿ0 j, and A CP is zero for B ! 0 .Several groups have estimated the correction due to the SU2 breaking effect of EWP contributions to be 1.5-2.0degrees [54,55].These estimates consider contributions from the two EWP operators assumed to be dominant in the effective Hamiltonian.(ii) The possible effect on the isospin analysis from 0 -! mixing [55] is neglected.(iii) Other SU2 symmetry breaking effects are neglected.Estimates of the magnitude of these effects are much less than the current experimental precision [56,57].Possible isospin 5=2 amplitudes also break the SU2 triangle construction [58].(iv) The isospin analysis outlined above neglects possible I 1 amplitudes [59].The presence of I 1 amplitudes in B 0 !ÿ can be tested by measuring S long and C long for different ranges of the invariant 0 mass.We constrain the CKM angle and penguin contribution from an isospin analysis of B ! decays.The inputs to the isospin analysis are the amplitudes of the CP-even longitudinal polarization of the final state, as well as the measured values of S long and C long for B 0 !ÿ .We use the following numerical inputs in the isospin analysis: (i) The measurements of BB 0 !ÿ , f L , S long , and C long presented here.(ii) The measurements of B and f L for B ! 0 from Ref. [21].(iii) The measurement of BB 0 ! 0 0from Ref. [23].
To interpret our results in terms of a constraint on from the isospin relations, we construct a 2 that includes the 8.The isospin triangle for B ! decays.

XI. CONCLUSIONS
We report the measurement of the branching fraction, f L , and CP violation parameters, S long and C long , for the decay B 0 !ÿ using a data sample of 383:6 4:2 10 6 B B pairs.We obtain the following results: BB 0 !ÿ 25:5 2:1stat 3:6 ÿ3:9 syst 10 ÿ6 ; f L 0:992 0:024stat 0:026 ÿ0:013 syst; S long ÿ0:17 0:20stat 0:05 ÿ0:06 syst; C long 0:01 0:15stat 0:06syst: Using these results, and experimental knowledge of the other B ! final states, we perform an isospin analysis to obtain a measurement of the CKM angle .The confidence level distribution for is shown in Fig. 9.The solution obtained that is compatible with the results of SMbased fits of existing data is 73:1; 117:0 at 68% CL.The uncertainty on the measurement of from the isospin analysis is dominated by penguin pollution.If one uses the flavor SU3 approach described in the text to constrain , one obtains the constraint 83:3; 105:8 at 68% CL.

ACKNOWLEDGMENTS
We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the excellent luminosity and machine conditions that have made this work possible.The success of this project also relies critically on the expertise and dedication of the computing organizations that support BABAR.The collaborating institutions wish to thank SLAC for its support and the kind hospitality extended to them.An extensive study of the B backgrounds, and the associated systematic error, has been presented earlier in this paper.A number of B backgrounds decay into the same final state as the signal.These are (i) B 0 ! a 1 , a 1 ! 0, (ii) B 0 ! 0, nonresonant, (iii) B 0 !ÿ 0 0 , nonresonant.The systematic error associated with these modes was evaluated by propagating the uncertainty on the branching ratio and CP asymmetries to the final measurements, with the appropriate description of acceptance, SCF, resolutions, and other reconstruction biases (see Sec. IX B).The likelihood described in Sec.VII does not account for possible interference between B 0 !ÿ and other B 0 ! 4 final states.
Acceptance, SCF, and other reconstruction biases are not taken into account in this study, and instead, the ÿ measurements are averaged in a perfect region of the phase space referred to as the '' ÿ band,'' on which the analysis technique is assumed to have a uniform sensitivity.This ÿ band corresponds to the kinematic selection of (see Sec. III B): 0:5 < m 0 < 1:0 GeV=c 2 and ÿ0:9 < cos i < 0:98.
The remainder of this appendix describes the final state wave function, decay amplitudes, effective CP asymmetries, and the estimate of the systematic uncertainty from neglecting interference between the signal and other B 0 ! 4 final states.

FIG. 3 .FIG. 4 (
FIG.3.A schematic view of the decay of a B meson via two mesons to a four-pion final state.The meson final states are shown in their rest frames, and is the angle between the decay planes of the mesons.

1 FIG. 5 (
FIG.5(color online).The left plot shows the distribution of N for off-peak data (dashed line) and longitudinally polarized signal MC simulated events (solid line).The right plot shows the signal efficiency (solid line) and continuum background efficiency (dashed line) distribution as a function of cut on the minimum value of N .

FIG. 6 (FIG. 7 (
FIG.6(color online).The distributions for the highest purity tagged events for the variables (a) m ES , (b) E, (c) cosine of the helicity angle, and (d) m 0 .The dashed lines are the sum of backgrounds, and the solid lines are the full PDF. FIG.
This work is supported by the US Department of Energy and National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), the Commissariat a `l'Energie Atomique and Institut National de Physique Nucle ´aire et de Physique des Particules (France), the Bundesministerium fu ¨r Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale di Fisica Nucleare (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), the Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, Ministerio de Educacio ´n y Ciencia (Spain), and the Science and Technology Facilities Council (United Kingdom).Individuals have received support from the Marie-Curie IEF program (European Union) and the A. P. Sloan Foundation.APPENDIX: INTERFERENCE OF B 0 !ÿ WITH OTHER B 0 ! 4 MODES

1 FIG. 10 (
FIG. 10 (color online).Confidence level on obtained from the flavor SU3 approach.The horizontal dashed lines correspond to the 68% (top) and 90% (bottom) CL intervals.The constraint shown in this figure uses some of the same inputs as the SU2 isospin constraint of Fig.9, and the two cannot be averaged.

TABLE III .
The components of charmless B backgrounds considered, along with the branching fraction (B) and number of events (N bg ) expected in the selected data sample, where y indicates a longitudinally polarized final state, and NR denotes a nonresonant final state.Where appropriate, branching fractions of these decays are quoted, including the branching fraction of subdecay modes and measured or expected values of f L .

TABLE IV .
The shifts on signal yield[Nsignal], fraction of longitudinally polarized events (f L ), S long (S long ), and C long (C long ) obtained when floating the yield of each fixed B background in turn.The symbol y indicates a longitudinally polarized final state, and NR denotes a nonresonant final state.The yields of decay modes not listed in this table are allowed to vary in the nominal fit.

TABLE V .
Summary of additive systematic uncertainty contributions on the signal yield N signal , f L , S long , and C long .

TABLE VI .
Systematic error coming from interference effects due to three contaminating amplitudes.