Measurements of branching fraction, polarization, and charge asymmetry of B-+/-->rho(+/-)f(0) and a search for B-+/-->rho(+/-)f(0)(980) RID C-2728-2008 RID C-5223-2009 RID C-5719-2008 RID D-1055-2009 RID A-2675-2009

B. Aubert, M. Bona, D. Boutigny, F. Couderc, Y. Karyotakis, J. P. Lees, V. Poireau, V. Tisserand, A. Zghiche, E. Grauges, A. Palano, J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu, G. Eigen, I. Ofte, B. Stugu, G. S. Abrams, M. Battaglia, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, M. S. Gill, Y. Groysman, R. G. Jacobsen, J. A. Kadyk, L. T. Kerth, Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, T. J. Orimoto, M. Pripstein, N. A. Roe, M. T. Ronan, W. A. Wenzel, P. del Amo Sanchez, M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, A. T. Watson, T. Held, H. Koch, B. Lewandowski, M. Pelizaeus, K. Peters, T. Schroeder, M. Steinke, J. T. Boyd, J. P. Burke, W. N. Cottingham, D. Walker, D. J. Asgeirsson, T. Cuhadar-Donszelmann, B. G. Fulsom, C. Hearty, N. S. Knecht, T. S. Mattison, J. A. McKenna, A. Khan, P. Kyberd, M. Saleem, D. J. Sherwood, 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, M. Bruinsma, M. Chao, S. Curry, I. Eschrich, D. Kirkby, A. J. Lankford, P. Lund, M. Mandelkern, R. K. Mommsen, W. Roethel, D. P. Stoker, S. Abachi, C. Buchanan, S. D. Foulkes, J. W. Gary, O. Long, B. C. Shen, K. Wang, L. Zhang, H. K. Hadavand, E. J. Hill, 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, G. Nesom, T. Schalk, B. A. Schumm, A. Seiden, P. Spradlin, D. C. Williams, M. G. Wilson, J. Albert, E. Chen, A. Dvoretskii, F. Fang, D. G. Hitlin, I. Narsky, T. Piatenko, F. C. Porter, A. Ryd, 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, W. O. Ruddick, J. G. Smith, K. A. Ulmer, S. R. Wagner, J. Zhang, A. Chen, E. A. Eckhart, 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, T. Brandt, V. Klose, H. M. Lacker, W. F. Mader, R. Nogowski, J. Schubert, K. R. Schubert, R. Schwierz, J. E. Sundermann, A. Volk, D. Bernard, G. R. Bonneaud, E. Latour, 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, G. Cibinetto, E. Luppi, M. Negrini, A. Petrella, L. Piemontese, E. Prencipe, 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, G. Brandenburg, K. S. Chaisanguanthum, M. Morii, J. Wu, R. S. Dubitzky, J. Marks, S. Schenk, U. Uwer, D. J. Bard, W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, 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, A. G. Denig, M. Fritsch, G. Schott, N. Arnaud, M. Davier, G. Grosdidier, A. Höcker, F. Le Diberder, V. Lepeltier, A. M. Lutz, A. Oyanguren, S. Pruvot, S. Rodier, P. Roudeau, M. H. Schune, A. Stocchi, W. F. Wang, G. Wormser, C. H. Cheng, D. J. Lange, D. M. Wright, C. A. Chavez, I. J. Forster, J. R. Fry, E. Gabathuler, R. Gamet, K. A. George, D. E. Hutchcroft, D. J. Payne, K. C. Schofield, C. Touramanis, A. J. Bevan, F. Di Lodovico, W. Menges, R. Sacco, G. Cowan, H. U. Flaecher, D. A. Hopkins, P. S. Jackson, T. R. McMahon, S. Ricciardi, 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, M. T. Naisbit, J. C. Williams, J. I. Yi, C. Chen, W. D. Hulsbergen, A. Jawahery, C. K. Lae, D. A. Roberts, G. Simi, G. Blaylock, C. Dallapiccola, S. S. Hertzbach, X. Li, T. B. Moore, S. Saremi, H. Staengle, R. Cowan, G. Sciolla, S. J. Sekula, M. Spitznagel, F. Taylor, R. K. Yamamoto, H. Kim, S. E. Mclachlin, P. M. Patel, S. H. Robertson, A. Lazzaro, V. Lombardo, 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, N. Cavallo, G. De Nardo, F. Fabozzi, C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, C. Sciacca, M. A. Baak, G. Raven, H. L. Snoek, C. P. Jessop, J. M. LoSecco, T. Allmendinger, G. Benelli, L. A. Corwin, K. K. Gan, K. Honscheid, D. Hufnagel, P. D. Jackson, H. Kagan, R. Kass, 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, PRL 97, 261801 (2006) P H Y S I C A L R E V I E W L E T T E R S week ending 31 DECEMBER 2006

We measure the branching fraction (B), polarization (f L ), and CP asymmetry (A CP ) of B ! 0 decays and search for the decay B ! f 0 980 based on a data sample of 231:8 10 6 4S ! B B decays collected with the BABAR detector at the SLAC PEP-II asymmetric-energy B factory. In B ! 0 decays we measure B 16:8 2:2 2:310 ÿ6 , f L 0:905 0:042 0:023 ÿ0:027 , and A CP ÿ0:12 0:13 0:10, and find an upper limit on the branching fraction of B ! f 0 980! ÿ decays of 1:9 10 ÿ6 at 90% confidence level. DOI The measurement of the CP-violating phase of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1] is an important part of the present program in particle physics. Violation of CP symmetry is manifested as a nonzero area of the CKM unitarity triangle [2]. In this Letter we report the measurement of the branching fraction, polarization, and CP asymmetry of the B ! 0 decay mode, which is needed for the isospin analysis used to extract argÿV td V tb =V ud V ub [3]. We also set an upper limit on the unknown branching fraction of B ! f 0 980! ÿ , which is measured to control this background to the B ! 0 analysis.
In B 0 B 0 ! ÿ decays [4] the interference between the B B oscillations which depend on V td and the dominating tree-level amplitude b ! u ud causes a time-dependent CP asymmetry that depends on sin2. The presence of loop (penguin) amplitudes leads to a shift j ÿ eff j, between the physical weak phase and the effective one eff , experimentally measured in B 0 ! ÿ decays [5,6]. However, the penguin amplitudes in these decays are known to contribute at a very low level because of the small upper limit of 1:1 10 ÿ6 at 90% confidence level (C.L.) [7], obtained from the branching fraction of the penguin dominated mode B 0 ! 0 0 . The size of can be extracted from the full isospin analysis combining all B ! modes [3].
In B ! decays, a spin zero particle decays into two spin one particles. The final state is therefore a superposition of two transversely polarized modes (helicity 1) and one longitudinal mode (helicity 0), which can be measured through an angular analysis. The longitudinal polarization fraction f L is defined as the fraction of decays to the helicity zero state f L ÿ L =ÿ, where ÿ is the total decay rate and ÿ L is the decay rate to the longitudinally polarized final state. The transverse polarization is a mixed CP state, while the longitudinal state is pure CP even. The previous measurements of f L [8,9] showed the decay is consistent with being fully longitudinally polarized.
Our analysis is performed in the helicity frame [10] as a function of the two helicity angles and 0 where the helicity angle of a 0 meson is defined as the angle between its daughter and the direction opposite to the B meson in the 0 rest frame. The polarization f L can be extracted from the differential decay rate: Here we integrate over the angle between the -meson decay planes. The measurements presented in this Letter are based on data collected with the BABAR detector [11] at the SLAC PEP-II asymmetric-energy e e ÿ collider. The analyzed data sample of 231:8 2:6 10 6 B B pairs produced at the 4S resonance corresponds to an integrated luminosity of 210:5 fb ÿ1 .
To reconstruct B ! 0 and B ! f 0 decays, we select events with at least three charged tracks and one neutral pion candidate. Charged tracks are required to originate from the interaction point and have particle identification information inconsistent with kaon, electron, and proton hypotheses. We form 0 ! candidates from pairs of calorimeter showers, each with a photonlike lateral spread and a minimum energy of 50 MeV. The invariant mass of 0 candidates is required to fall in the range 0:10 < m < 0:16 GeV=c 2 .
The mass of charged candidates must satisfy 0:396 < m 0 < 1:146 GeV=c 2 where the lowside re-quirement on the 0 mass is chosen to exclude K 0 S ! ÿ decays. Neutral final state meson candidates ( 0 , f 0 ) must satisfy 0:520 < m ÿ < 1:146 GeV=c 2 . In order to suppress backgrounds with low momentum pions, the helicity angles are required to fall in the ranges ÿ0:8 < cos < 0:95 and j cos 0 j < 0:95. Backgrounds from D 0 ! K ÿ 0 and D 0 ! ÿ 0 decays are reduced by requiring the candidate D 0 invariant mass to be at least 40 MeV=c 2 away from the D 0 mass. About 20% of the selected events have multiple B candidates and the one that has the reconstructed 0 mass closest to the 0 mass is kept. In the case that more than one candidate has the same reconstructed 0 mass, we select one at random.
Continuum decays represent the largest source of background and are reduced by requiring j cos T j < 0:8, where T is the cosine of the angle between the B thrust axis and that from the rest of the event (ROE). To further discriminate signal from continuum, we also use a neural network built out of five event-shape variables: a Fisher discrimi-nant combining the 0th and 2nd order monomials [12] for charged particles and neutral clusters of the ROE; the cosine of the angle between the direction of the B and the collision axis (z) in the center-of-mass (c.m.) frame; the cosine of the angle between the B thrust axis and the z axis; the variable j cos T j defined above; and the sum of transverse momenta in the ROE relative to the z axis. The output is transformed into a variable x NN which has roughly Gaussian signal and background distributions. We select candidates in a range of x NN that removes 54% of continuum background events while retaining 90% of the signal. After these selections, about 85% of the remaining events are from continuum decays.
Signal event candidates are further identified based on two kinematic variables: the beam-energy-substituted After the selection criteria are applied, the efficiency is 8.4% for longitudinal and 18.6% for transverse polarized B ! 0 decays. The selection efficiency is 16.6% for B ! f 0 decays. Any possible interference effects between the B ! 0 and B ! f 0 are neglected.
An unbinned extended maximum likelihood fit is applied to the selected sample of N tot 74293 events in order to measure the B ! 0 event yield, polarization, and charge asymmetry as well as the B ! f 0 event yield. The likelihood function is where M is the number of hypotheses (signal, misreconstructed signal, continuum, and B-background classes), and n k (n j ) represents the number of measured events for each hypothesis determined by maximizing the likelihood function. P j x i is the product of the probability density functions (PDFs) of hypothesis j evaluated at the ith event's measured variables,x i fm ES ; E; m 0 ; m ÿ ; cos ; cos 0 ; x NN g. In addition, the charge asymmetry, obtained from the measured B ÿ and B signal candidate decay yields, A CP Each discriminating variable in the likelihood function is modeled with a PDF extracted either from the data, or from high statistics Monte Carlo (MC) simulated data samples. The correlations between the variables are assumed to be small and the PDFs independent. This is checked with systematic error studies, and corrections are applied where necessary.
The continuum background E, m ES , and x NN distributions are modeled with one-dimensional parametrized distributions taken from fits to the data. Correlations are observed between the m and cos distributions for both -meson candidates, which are taken into account with two-dimensional PDFs. The signal component is modeled with one-dimensional parametrized distributions for each of six variables; m ES is modeled with a nonparametric PDF [13]. The signal PDF shapes are obtained from fits to signal MC sample after the selection is applied. Events with a true B ! 0 decay but with wrong tracks or calorimeter clusters assigned to the final state are referred to as self cross feed (SCF) events. They make up 35% and 14% of the selected longitudinally and transversely polarized signal samples, respectively. The longitudinal and transverse SCF components and B-background PDFs are determined in a similar manner using high statistics MC samples and modeled with nonparametric PDFs [13] for each variable.
To , and BK 0 ! K ÿ 2=3. The decays B ! a 1 and B ! a 1 have few experimental constraints [17,18]. We adopt the following B branching ratios, in units of 10 ÿ6 , and assume a 100% systematic uncertainty: a 0 1 12, a 1 0 6, a 0 1 a 1 0 48. Table I  We considered systematic effects from biases in the fit model, which are due to the imperfect modeling of some correlations between fitted variables in the likelihood function. Tests of the fit are made by using a large number of MC samples containing the amounts of signal, continuum, and B-background events measured or fixed in the data fit, and where correlations between variables are modeled for signal and B backgrounds. The fits to these samples should reproduce the number of MC events generated. A shift of the fitted values with respect to the generated ones indicates a bias of ÿ49 events on the 0 yield, 0:009 on f L , and ÿ4:6 events on the f 0 yield in the fitting procedure. We use the same technique to study the effects of correlations between the neural net and helicity variables in the q q continuum and observe a fit bias of 24 events on the 0 yield, ÿ0:001 on f L , and ÿ19 events on the f 0 yield. No other significant correlations were observed between the other discriminating variables. These two biases are corrected from the fit measurements, and half of each separate fit bias is taken as the systematic error (cf. Table II).
Many of the B-background rates are poorly known. The effect of uncertainties in these values is evaluated by varying the number of events in each background category within the range allowed by the error on the branching fraction. Fourteen nonresonant backgrounds that are not in the default fit are tested by adding them singly to the fit with a yield that is allowed to vary. The only shift seen was associated with the mode B ! 0 0 , and is taken as a symmetric systematic uncertainty.
The systematic error associated with misreconstructed signal is evaluated by taking the difference between the default fit and the one for which these events are not modeled, and therefore mostly absorbed into the ''other charmless'' background category. We consider the error due to the uncertainty on the signal, B background, and continuum PDF shapes and estimate a systematic error by varying the parameters obtained from MC calculations that govern these shapes within their statistical uncertainty. The impact of the uncertainty on the measurement of the f 0 mass and width [19] has also been evaluated. The values of the systematic errors described above are given in Table II. Systematic uncertainties in the reconstruction and calibration procedure introduce a systematic error of 3% after a correction of ÿ2:5% on the 0 reconstruction efficiency, 3.9% after a correction of ÿ1:5% on the track reconstruction efficiency, and a systematic error of 1.1% from the particle identification. The uncertainty on the efficiency ratio between longitudinal and transverse events is found to be negligible. The error on A CP includes a 0.45% uncertainty in the charged track reconstruction asymmetry, a 4% uncertainty from the detector's intrinsic charged particle identification asymmetry, and a 9% uncertainty which is the largest single shift obtained when assuming a uniform probability for the charge asymmetry of every B-background individually.  In summary, we measure the branching fraction, longitudinal polarization, and CP asymmetry of the decay B ! 0 , using a dataset of about 231: 8  The measurement of the branching fraction has improved by a factor of about 2 with respect to the previous BABAR measurement [8], and supersedes it. The isospin relations between branching ratios are consistent between this measurement and those of ÿ and 0 0 [14], validating the approach [3] used to constraint . Moreover, our measurements confirm that this mode is largely longitudinally polarized. They also confirm that the charge asymmetry is consistent with zero as expected for decays proceeding through one decay channel only; this suggests the contributions of electroweak penguins are small in the B ! system.
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. *Also with Università di Perugia, Dipartimento di Fisica, Perugia, Italy †