Measurements of the absolute branching fractions of B-+/--> K-+/- Xc(c)over-bar-$ RID C-2728-2008 RID C-5223-2009 RID C-5719-2008 RID D-1055-2009

B. Aubert, R. Barate, D. Boutigny, F. Couderc, Y. Karyotakis, J. P. Lees, V. Poireau, V. Tisserand, A. Zghiche, E. Grauges, A. Palano, M. Pappagallo, A. Pompili, J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu, G. Eigen, I. Ofte, B. Stugu, G. S. Abrams, M. Battaglia, A. B. Breon, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, J. Kadyk, L. T. Kerth, Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, P. J. Oddone, T. J. Orimoto, M. Pripstein, N. A. Roe, M. T. Ronan, W. A. Wenzel, M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, S. E. Morgan, A. T. Watson, M. Fritsch, K. Goetzen, T. Held, H. Koch, B. Lewandowski, M. Pelizaeus, K. Peters, T. Schroeder, M. Steinke, J. T. Boyd, J. P. Burke, N. Chevalier, W. N. Cottingham, T. Cuhadar-Donszelmann, B. G. Fulsom, C. Hearty, N. S. Knecht, T. S. Mattison, J. A. McKenna, A. Khan, P. Kyberd, M. Saleem, L. Teodorescu, A. E. Blinov, V. E. Blinov, A. D. Bukin, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko, A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, A. N. Yushkov, D. Best, 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, C. Buchanan, B. L. Hartfiel, A. J. R. Weinstein, S. D. Foulkes, J. W. Gary, O. Long, B. C. Shen, K. Wang, L. Zhang, D. del Re, H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, S. Rahatlou, V. Sharma, J. W. Berryhill, C. Campagnari, A. Cunha, B. Dahmes, T. M. Hong, M. A. Mazur, J. D. Richman, W. Verkerke, 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, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, J. S. Minamora, I. Narsky, T. Piatenko, F. C. Porter, A. Ryd, A. Samuel, R. Andreassen, G. Mancinelli, B. T. Meadows, M. D. Sokoloff, F. Blanc, P. Bloom, S. Chen, W. T. Ford, J. F. Hirschauer, A. Kreisel, 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, Q. Zeng, D. Altenburg, E. Feltresi, A. Hauke, B. Spaan, T. Brandt, J. Brose, M. Dickopp, V. Klose, H. M. Lacker, R. Nogowski, S. Otto, A. Petzold, J. Schubert, K. R. Schubert, R. Schwierz, J. E. Sundermann, D. Bernard, G. R. Bonneaud, P. Grenier, S. Schrenk, Ch. Thiebaux, G. Vasileiadis, M. Verderi, D. J. Bard, P. J. Clark, W. Gradl, F. Muheim, S. Playfer, Y. Xie, M. Andreotti, V. Azzolini, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinetto, E. Luppi, M. Negrini, L. Piemontese, F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, P. Patteri, I. M. Peruzzi,* M. Piccolo, A. Zallo, A. Buzzo, R. Capra, R. Contri, M. Lo Vetere, M. Macri, M. R. Monge, S. Passaggio, C. Patrignani, E. Robutti, A. Santroni, S. Tosi, G. Brandenburg, K. S. Chaisanguanthum, M. Morii, E. Won, J. Wu, R. S. Dubitzky, U. Langenegger, J. Marks, S. Schenk, U. Uwer, G. Schott, W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, R. L. Flack, J. R. Gaillard, J. A. Nash, M. B. Nikolich, W. Panduro Vazquez, X. Chai, M. J. Charles, W. F. Mader, U. Mallik, A. K. Mohapatra, V. Ziegler, J. Cochran, H. B. Crawley, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, J. Yi, N. Arnaud, M. Davier, X. Giroux, G. Grosdidier, A. Höcker, F. Le Diberder, V. Lepeltier, A. M. Lutz, A. Oyanguren, T. C. Petersen, S. Plaszczynski, S. Rodier, P. Roudeau, M. H. Schune, A. Stocchi, G. Wormser, C. H. Cheng, D. J. Lange, M. C. Simani, D. M. Wright, A. J. Bevan, C. A. Chavez, I. J. Forster, J. R. Fry, E. Gabathuler, R. Gamet, K. A. George, D. E. Hutchcroft, R. J. Parry, D. J. Payne, K. C. Schofield, C. Touramanis, C. M. Cormack, F. Di Lodovico, W. Menges, R. Sacco, C. L. Brown, G. Cowan, H. U. Flaecher, M. G. Green, D. A. Hopkins, P. S. Jackson, T. R. McMahon, S. Ricciardi, F. Salvatore, D. Brown, C. L. Davis, J. Allison, N. R. Barlow, R. J. Barlow, C. L. Edgar, M. C. Hodgkinson, M. P. Kelly, G. D. Lafferty, M. T. Naisbit, J. C. Williams, C. Chen, W. D. Hulsbergen, A. Jawahery, D. Kovalskyi, C. K. Lae, D. A. Roberts, G. Simi, G. Blaylock, C. Dallapiccola, S. S. Hertzbach, R. Kofler, V. B. Koptchev, X. Li, T. B. Moore, S. Saremi, H. Staengle, S. Willocq, R. Cowan, K. Koeneke, G. Sciolla, S. J. Sekula, M. Spitznagel, F. Taylor, R. K. Yamamoto, H. Kim, P. M. Patel, S. H. Robertson, A. Lazzaro, V. Lombardo, F. Palombo, J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Kroeger, J. Reidy, D. A. Sanders, D. J. Summers, H. W. Zhao, S. Brunet, D. Côté, P. Taras, B. Viaud, H. Nicholson, N. Cavallo, G. De Nardo, F. Fabozzi, C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, C. Sciacca, M. Baak, H. Bulten, G. Raven, H. L. Snoek, L. Wilden, C. P. Jessop, J. M. LoSecco, T. Allmendinger, G. Benelli,

Several exclusive decays of B mesons of the form B ! K X c c (where X c c is one of the charmonium states c , J= , c0 , c1 , 0 c , 0 , 00 ) have been observed by reconstructing the charmonium state from its decay to some known final state, f [1,2].In principle, such B decays provide a direct probe of charmonium properties since the phase space is large for all known states and all should be produced roughly equally, in the absence of a strong selection rule [3].However, with this technique only the product of the two branching fractions BB !K X c c BX c c ! f is measured, thereby reducing the precision of BB !K X c c when the daughter branching fraction is poorly known.
We describe here a complementary approach, based on the measurement of the kaon momentum spectrum in the B center-of-mass frame, where two-body decays can be identified by their characteristic monochromatic line, allowing an absolute determination of BB !K X c c .Knowledge of the B center-of-mass system is obtained by exclusive reconstruction of the other B meson from a 4S decay.In addition to obtaining new information on known charmonium states, this method is used to search for the X3872 state, recently observed in B ! K X3872 decays by Belle [4] and BABAR [5], in the subsequent decay X3872 !J= ÿ .The same method allows a search for charged partners of the X3872 in B 0 decays, independent of the X3872 decay mode.The nature of X3872 resonance is still unclear; different interpretations [6] have been proposed but more experimental data will be needed to discriminate between them.
For this analysis we use a data sample of 210:5 fb ÿ1 integrated luminosity, corresponding to 231:8 10 6 B B pairs.The data have been collected with the BABAR detector at the SLAC PEP-II asymmetric-energy collider, where 9 GeV electrons and 3.1 GeV positrons collide at a center-of-mass energy 10.58 GeV, corresponding to the mass of the 4S resonance.A detailed description of the BABAR detector can be found in [7].Charged tracks are reconstructed with a 5 layer silicon vertex tracker (SVT) and a 40 layer drift chamber (DCH), located in a 1.5 T magnetic field generated by a superconducting solenoid.The energy of photons and electrons is measured with an electromagnetic calorimeter made up of CsI(Tl) crystals.Charged hadron identification is done with ionization measurements in the SVT and DCH and with an internally reflecting ring imaging Cherenkov detector.The instru-mented flux return of the solenoid is used to identify muons.
The analysis is performed on a sample of events where a B meson is fully reconstructed (B recon ).For these events, the momentum of the other B (B signal ) can be calculated from the momentum of B recon and the beam parameters.We select events with a K not used for the reconstruction of B recon and calculate its momentum (p K ) in the B signal center-of-mass system.
B recon mesons are reconstructed in their decays to exclusive D H final states, where H is one of several combinations of , K , 0 , and K 0 S hadrons; a detailed description of the method can be found in [8].
The number of B events in the data is determined with a fit to the distribution of the beam energy substituted mass , where E CM is the total center-ofmass energy, determined from the beam parameters, and p B is the measured momentum of B recon in the center-ofmass frame.The fit function is the sum of a Crystal Ball function [9] describing the signal and an ARGUS function [10] for each background component (e e ÿ !q q where q is u, d, s, or c or misreconstructed Bs), the relative weights of which are obtained from a Monte Carlo simulation (MC), while the total normalization factor is determined from the data.A total of 378 580 1110 events with a fully reconstructed B is obtained.
Fifteen variables related to the B recon decay characteristics, its production kinematics, the topology of the full event, and the angular correlation between B recon and the rest of the event are used in a neural network (NN1) to reduce the large background, mainly due to non-B events.The network has 80% signal efficiency while rejecting 90% of the background.The m ES distribution after this selection is shown in Fig. 1.Only events with 5:275 < m ES < 5:285 GeV=c 2 are used in the analysis.
We now consider only tracks not associated with B recon .Most K produced in B decays originate from D mesons and their spectrum, although broad, peaks at low p K .In the B rest frame, these K are embedded in a ''minijet'' of D decay products, while signal K recoil against a massive (3-4 GeV=c 2 ) state and therefore tend to be more isolated.A second neural network (NN2) rejects background from secondary K by using 15 input variables describing the energy and track multiplicities measured in the K hemisphere, the sphericity of the recoil system, and the angular PRL 96, 052002 ( 2006) week ending 10 FEBRUARY 2006 052002-4 correlations between the K and the recoil system.These variables have been chosen to be independent of the particular decay topology of the recoil system.Since the topology of the event changes with the recoil mass, we have considered separately two recoil mass regions in the training of this neural network: the ''high-mass'' region, corresponding to 1:0 < p K < 1:5 GeV=c and the ''lowmass'' region, for 1:5 < p K < 2:0 GeV=c.The signal training sample is B ! K X c c MC simulation while the background sample consists of simulated K from D meson decays in the same momentum range.The chosen cuts on the NN2 outputs correspond to 85% signal efficiency; the background rejection factor varies between 2.5 in the X3872 and 0 region and 1.5 in the J= region.The selection criteria are optimized for MC signal significance with the high-mass region blinded.
The kaon momentum distribution shows a series of peaks due to the two-body decays B ! K X c c corresponding to the different X c c masses, superimposed on a smooth spectrum due to K coming from multibody B decays, or non-B background.The mass of the X c c state (m X ) can be calculated directly from , where m B and m K are the B and K masses and E K is the K energy.The resonance width ÿ X can be obtained from the Breit-Wigner width of the peak in the p K spectrum ÿ K , obtained after deconvolution with the momentum resolution function, using ÿ X ÿ K K m B =m X , where K p K =E K .
We determine the number of B ! K X c c events (N X ) from a fit to the p K distribution.The branching fraction for the decay channel is calculated as: where X is the efficiency determined from the MC simu-lation and N B the number of B mesons in the sample.An alternative method, which we use to improve the branching fraction measurement in the case of c , is to normalize to the channel B ! K J= , which is well measured in the literature [11], according to: In this relative measurement, the systematic errors that are common to both resonances cancel in the ratio.The two methods are combined to extract BB !K c , taking into account the correlations between them.We fit the p K spectrum using an unbinned maximum likelihood method.The background is well modeled by a third degree polynomial and each signal is a Breit-Wigner function folded with a resolution function.The masses and widths of the c and 0 c mesons are left free; all others are fixed to values from Ref. [11].The resolution function has two parts: a Gaussian with varying from 6 MeV=c at p K ' 1:1 GeV=c to 12 MeV=c at p K ' 1:7 GeV=c describes the 72.5% of the signal where B recon is correctly reconstructed; if B recon is incorrect, but has m ES within our range, the p K resolution is a bifurcated Gaussian with 78 and 52 MeV=c on the left-and right-hand side of the peak, respectively.
The spectrum in the low-mass region is expected to exhibit two peaks, at p K 1:683 GeV=c corresponding to the J= , and at p K 1:754 GeV=c for the c meson.These two peaks are clearly seen in Fig. 2(a); both have a significance of 7. The number of events under each peak obtained from the fit is NJ= 259 41 and N c 273 43.
The spectrum in the high-mass region is fitted with a background and seven signal functions, corresponding to the following states: 0 , c0 , c1 , c2 , 00 , 0 c , and X3872.The resulting fit is shown in Fig. 2(b), with the yields given in Table I.The h c charmonium state lies near the c1 , and it is difficult to distinguish the peaks from these two decays.A fit including the h c yields a number of h c events consistent with zero, and a fit performed with free c1 mass and width gives values consistent with a narrow c1 ; therefore, we have no evidence for h c production.
Several sources of systematic error affecting these measurements have been evaluated.The relative errors on absolute measurements are the same for all states; many of these cancel partially in relative measurements, and all are summarized in Table II.''B counting'' refers to uncertainties in the fit parametrization used to determine the number of fully reconstructed B recon .It is one of the largest errors in absolute measurements, and cancels in ratios.The mass scale is verified to a precision of 1:5 MeV=c in p K by floating the masses of the well-measured J= , c1 , and 0 peaks; we assign a systematic error corresponding to this shift.We also consider variations in the background and signal model parametrizations, which partially cancel in the case of ratios.Errors in the K track reconstruction and identification efficiency are evaluated by comparing data and MC control samples.The systematic error in the NN1 and NN2 selections is evaluated by comparing efficiencies and distributions in data and MC, and studying efficiency variation with p K .We verified that the NN2 selection is not dependent on visible energy or multiplicity of the recoil part of the B meson decay.Adding in quadrature, the total relative error on an absolute measurement is 9.0%.The total is reduced to 3.3% for the relative measurement of J= and c , and to 5.9% for states in the high-mass region relative to J= .For the extraction of relative branching fractions, an additional 4% error, labeled (ext) in the following, comes from the present uncertainty of BB !K J= 10:0 0:4 10 ÿ4 [11].
In the high-mass region, clear signals are found for c1 and 0 (with significance 6.0 and 3:2, respectively), an excess of events is present for 0 c and 00 [12], while no signal is found for c0 , c2 , and X3872.The branching fractions and upper limits are summarized in Table I.
In the low-mass region, our J= measurement is consistent with the world average.From the c and J= yields and the reference branching fraction we can derive the result with the relative measurement method BB !K c rel 10:6 2:3stat 0:4sys 0:4ext 10 ÿ4 .We combine this result with the absolute measurement of Table I, taking the correlated errors into account, to obtain BB !K c 8:7 1:5 10 ÿ4 .
We obtain from our fits the c and 0 c masses and widths and find m c 2982 5 MeV=c 2 , ÿ c < 43 MeV, and m 0 c 3639 7 MeV=c 2 , ÿ 0 c < 23 MeV, where the width limits are both at 90% C.L.
Taking BB !K X3872 < 3:2 10 ÿ4 , and using an average of the Belle [4] and BABAR [5] measurements of BB !K X3872 BX3872 !J= ÿ we set a lower limit BX3872 !J= ÿ > 4:2% at 90% C.L.This branching fraction, for which there are not yet any predictions, is sensitive to the distribution of charm quarks inside the X3872.A search for charged partners of the X3872 is performed by examining K recoiling from a sample of 245.6 k reconstructed B 0 decays.No signal is seen and we find BB 0 !K X3872 < 5 10 ÿ4 at 90% C.L.   We combine our BB !K c with a previous BABAR measurement of BB !K c B c ! K K [13] to obtain B c ! K K 8:5 1:8%, significantly improving the precision of the world average [11].Since this branching fraction is used as a reference for all c yield measurements, our result will lead to more precise c partial widths and more stringent comparisons with theoretical models.For example, from an average of BJ= !c B c ! K K measured by Mark-III [14], DM2 [15], and BES [16], we obtain BJ= !c 0:79 0:20%, and using the value ÿ c ! B c ! K K 0:48 0:06 keV [11] we calculate ÿ c ! 5:6 1:4 keV.Both results are more precise than the world average [11].Similarly, we obtain B 0 c ! K K 8 5% and ÿ 0 c ! 0:9 0:5 keV.
In conclusion, a novel technique is used to measure directly the absolute branching fractions of the various charmonium states X c c in two-body decays B ! K X c c (Table I).The results for X c c c ; J= ; 0 are in agreement with previous measurements, and the c result significantly improves the present world average.Upper limits are set for c0 and c2 , confirming factorization suppression [17].Measurements of B ! K 0 c and B ! K 00 branching fractions are reported, although with poor significance.Upper limits are given for X3872 and for production of a possible charged partner in B 0 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.This work is supported by DOE and NSF (USA), NSERC (Canada), IHEP (China), CEA and CNRS-IN2P3 (France), BMBF and DFG (Germany), INFN (Italy), FOM (The Netherlands), NFR (Norway), MIST (Russia), and PPARC (United Kingdom).Individuals have received support from CONACyT (Mexico), A. P. Sloan Foundation, Research Corporation, and Alexander von Humboldt Foundation.

FIG. 1 .
FIG.1.The m ES distribution of all B recon after the NN1 selection.The solid line represents a fit described in the text; cumulative background contributions from e e ÿ !q q where q is u, d, s, or c (dashed line), B 0 (dotted line), B (dash-dotted line) events are shown.The arrows indicate the cuts used in the analysis (see text).

FIG. 2 (
FIG. 2 (color online).Kaon momentum spectrum for the (a) low-mass and (b) high-mass regions.The lines represent the fit described in the text.Arrows show the expected positions of known charmonium states.

TABLE II .
Summary of systematic errors in percent for absolute and the J= : c relative measurement.

TABLE I .
Event yields and absolute branching fractions BB !K X c c from the fits to the p K spectrum.The first error is statistical, the second systematic, and B upper limits are given at 90% C.L., taking into account the 9.0% systematic error.