Observation of the decay B^{0}_{s}\rightarrow K^{0}\overline{K}{}^{\,0}


Observation of the decay

B. Pal University of Cincinnati, Cincinnati, Ohio 45221    A. J. Schwartz University of Cincinnati, Cincinnati, Ohio 45221    A. Abdesselam Department of Physics, Faculty of Science, University of Tabuk, Tabuk 71451    I. Adachi High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    H. Aihara Department of Physics, University of Tokyo, Tokyo 113-0033    D. M. Asner Pacific Northwest National Laboratory, Richland, Washington 99352    T. Aushev Moscow Institute of Physics and Technology, Moscow Region 141700 Institute for Theoretical and Experimental Physics, Moscow 117218    R. Ayad Department of Physics, Faculty of Science, University of Tabuk, Tabuk 71451    T. Aziz Tata Institute of Fundamental Research, Mumbai 400005    V. Babu Tata Institute of Fundamental Research, Mumbai 400005    I. Badhrees Department of Physics, Faculty of Science, University of Tabuk, Tabuk 71451 King Abdulaziz City for Science and Technology, Riyadh 11442    S. Bahinipati Indian Institute of Technology Bhubaneswar, Satya Nagar 751007    A. M. Bakich School of Physics, University of Sydney, NSW 2006    E. Barberio School of Physics, University of Melbourne, Victoria 3010    P. Behera Indian Institute of Technology Madras, Chennai 600036    V. Bhardwaj University of South Carolina, Columbia, South Carolina 29208    B. Bhuyan Indian Institute of Technology Guwahati, Assam 781039    J. Biswal J. Stefan Institute, 1000 Ljubljana    A. Bobrov Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    A. Bozek H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342    M. Bračko University of Maribor, 2000 Maribor J. Stefan Institute, 1000 Ljubljana    T. E. Browder University of Hawaii, Honolulu, Hawaii 96822    D. Červenkov Faculty of Mathematics and Physics, Charles University, 121 16 Prague    V. Chekelian Max-Planck-Institut für Physik, 80805 München    A. Chen National Central University, Chung-li 32054    B. G. Cheon Hanyang University, Seoul 133-791    R. Chistov Institute for Theoretical and Experimental Physics, Moscow 117218    K. Cho Korea Institute of Science and Technology Information, Daejeon 305-806    V. Chobanova Max-Planck-Institut für Physik, 80805 München    Y. Choi Sungkyunkwan University, Suwon 440-746    D. Cinabro Wayne State University, Detroit, Michigan 48202    J. Dalseno Max-Planck-Institut für Physik, 80805 München Excellence Cluster Universe, Technische Universität München, 85748 Garching    N. Dash Indian Institute of Technology Bhubaneswar, Satya Nagar 751007    Z. Doležal Faculty of Mathematics and Physics, Charles University, 121 16 Prague    Z. Drásal Faculty of Mathematics and Physics, Charles University, 121 16 Prague    A. Drutskoy Institute for Theoretical and Experimental Physics, Moscow 117218 Moscow Physical Engineering Institute, Moscow 115409    D. Dutta Tata Institute of Fundamental Research, Mumbai 400005    S. Eidelman Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    H. Farhat Wayne State University, Detroit, Michigan 48202    J. E. Fast Pacific Northwest National Laboratory, Richland, Washington 99352    B. G. Fulsom Pacific Northwest National Laboratory, Richland, Washington 99352    V. Gaur Tata Institute of Fundamental Research, Mumbai 400005    A. Garmash Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    R. Gillard Wayne State University, Detroit, Michigan 48202    Y. M. Goh Hanyang University, Seoul 133-791    P. Goldenzweig Institut für Experimentelle Kernphysik, Karlsruher Institut für Technologie, 76131 Karlsruhe    D. Greenwald Department of Physics, Technische Universität München, 85748 Garching    O. Grzymkowska H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342    J. Haba High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    T. Hara High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    K. Hayasaka Kobayashi-Maskawa Institute, Nagoya University, Nagoya 464-8602    H. Hayashii Nara Women’s University, Nara 630-8506    X. H. He Peking University, Beijing 100871    W.-S. Hou Department of Physics, National Taiwan University, Taipei 10617    K. Inami Graduate School of Science, Nagoya University, Nagoya 464-8602    A. Ishikawa Department of Physics, Tohoku University, Sendai 980-8578    Y. Iwasaki High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801    W. W. Jacobs Indiana University, Bloomington, Indiana 47408    I. Jaegle University of Hawaii, Honolulu, Hawaii 96822    H. B. Jeon Kyungpook National University, Daegu 702-701    D. Joffe Kennesaw State University, Kennesaw GA 30144    K. K. Joo Chonnam National University, Kwangju 660-701    T. Julius School of Physics, University of Melbourne, Victoria 3010    K. H. Kang Kyungpook National University, Daegu 702-701    E. Kato Department of Physics, Tohoku University, Sendai 980-8578    T. Kawasaki Niigata University, Niigata 950-2181    C. Kiesling Max-Planck-Institut für Physik, 80805 München    D. Y. Kim Soongsil University, Seoul 156-743    H. J. Kim Kyungpook National University, Daegu 702-701    K. T. Kim Korea University, Seoul 136-713    M. J. Kim Kyungpook National University, Daegu 702-701    S. H. Kim Hanyang University, Seoul 133-791    K. Kinoshita University of Cincinnati, Cincinnati, Ohio 45221    P. Kodyš Faculty of Mathematics and Physics, Charles University, 121 16 Prague    S. Korpar University of Maribor, 2000 Maribor J. Stefan Institute, 1000 Ljubljana    P. Križan Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana J. Stefan Institute, 1000 Ljubljana    P. Krokovny Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    T. Kuhr Ludwig Maximilians University, 80539 Munich    R. Kumar Punjab Agricultural University, Ludhiana 141004    T. Kumita Tokyo Metropolitan University, Tokyo 192-0397    A. Kuzmin Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    Y.-J. Kwon Yonsei University, Seoul 120-749    I. S. Lee Hanyang University, Seoul 133-791    C. H. Li School of Physics, University of Melbourne, Victoria 3010    H. Li Indiana University, Bloomington, Indiana 47408    L. Li University of Science and Technology of China, Hefei 230026    L. Li Gioi Max-Planck-Institut für Physik, 80805 München    J. Libby Indian Institute of Technology Madras, Chennai 600036    D. Liventsev CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801    P. Lukin Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    T. Luo University of Pittsburgh, Pittsburgh, Pennsylvania 15260    M. Masuda Earthquake Research Institute, University of Tokyo, Tokyo 113-0032    D. Matvienko Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    K. Miyabayashi Nara Women’s University, Nara 630-8506    H. Miyata Niigata University, Niigata 950-2181    R. Mizuk Institute for Theoretical and Experimental Physics, Moscow 117218 Moscow Physical Engineering Institute, Moscow 115409    G. B. Mohanty Tata Institute of Fundamental Research, Mumbai 400005    S. Mohanty Tata Institute of Fundamental Research, Mumbai 400005 Utkal University, Bhubaneswar 751004    A. Moll Max-Planck-Institut für Physik, 80805 München Excellence Cluster Universe, Technische Universität München, 85748 Garching    H. K. Moon Korea University, Seoul 136-713    T. Mori Graduate School of Science, Nagoya University, Nagoya 464-8602    R. Mussa INFN - Sezione di Torino, 10125 Torino    E. Nakano Osaka City University, Osaka 558-8585    M. Nakao High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    T. Nanut J. Stefan Institute, 1000 Ljubljana    Z. Natkaniec H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342    M. Nayak Indian Institute of Technology Madras, Chennai 600036    N. K. Nisar Tata Institute of Fundamental Research, Mumbai 400005 Aligarh Muslim University, Aligarh 202002    S. Nishida High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    S. Ogawa Toho University, Funabashi 274-8510    S. Okuno Kanagawa University, Yokohama 221-8686    P. Pakhlov Institute for Theoretical and Experimental Physics, Moscow 117218 Moscow Physical Engineering Institute, Moscow 115409    G. Pakhlova Moscow Institute of Physics and Technology, Moscow Region 141700 Institute for Theoretical and Experimental Physics, Moscow 117218    C. W. Park Sungkyunkwan University, Suwon 440-746    H. Park Kyungpook National University, Daegu 702-701    S. Paul Department of Physics, Technische Universität München, 85748 Garching    T. K. Pedlar Luther College, Decorah, Iowa 52101    L. Pesántez University of Bonn, 53115 Bonn    R. Pestotnik J. Stefan Institute, 1000 Ljubljana    M. Petrič J. Stefan Institute, 1000 Ljubljana    L. E. Piilonen CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061    C. Pulvermacher Institut für Experimentelle Kernphysik, Karlsruher Institut für Technologie, 76131 Karlsruhe    J. Rauch Department of Physics, Technische Universität München, 85748 Garching    E. Ribežl J. Stefan Institute, 1000 Ljubljana    M. Ritter Ludwig Maximilians University, 80539 Munich    A. Rostomyan Deutsches Elektronen–Synchrotron, 22607 Hamburg    S. Ryu Seoul National University, Seoul 151-742    H. Sahoo University of Hawaii, Honolulu, Hawaii 96822    Y. Sakai High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    S. Sandilya Tata Institute of Fundamental Research, Mumbai 400005    T. Sanuki Department of Physics, Tohoku University, Sendai 980-8578    Y. Sato Graduate School of Science, Nagoya University, Nagoya 464-8602    V. Savinov University of Pittsburgh, Pittsburgh, Pennsylvania 15260    T. Schlüter Ludwig Maximilians University, 80539 Munich    O. Schneider École Polytechnique Fédérale de Lausanne (EPFL), Lausanne 1015    G. Schnell University of the Basque Country UPV/EHU, 48080 Bilbao IKERBASQUE, Basque Foundation for Science, 48013 Bilbao    C. Schwanda Institute of High Energy Physics, Vienna 1050    Y. Seino Niigata University, Niigata 950-2181    K. Senyo Yamagata University, Yamagata 990-8560    O. Seon Graduate School of Science, Nagoya University, Nagoya 464-8602    I. S. Seong University of Hawaii, Honolulu, Hawaii 96822    V. Shebalin Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    T.-A. Shibata Tokyo Institute of Technology, Tokyo 152-8550    J.-G. Shiu Department of Physics, National Taiwan University, Taipei 10617    B. Shwartz Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    F. Simon Max-Planck-Institut für Physik, 80805 München Excellence Cluster Universe, Technische Universität München, 85748 Garching    Y.-S. Sohn Yonsei University, Seoul 120-749    A. Sokolov Institute for High Energy Physics, Protvino 142281    E. Solovieva Institute for Theoretical and Experimental Physics, Moscow 117218    S. Stanič University of Nova Gorica, 5000 Nova Gorica    M. Starič J. Stefan Institute, 1000 Ljubljana    J. Stypula H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342    M. Sumihama Gifu University, Gifu 501-1193    T. Sumiyoshi Tokyo Metropolitan University, Tokyo 192-0397    U. Tamponi INFN - Sezione di Torino, 10125 Torino University of Torino, 10124 Torino    Y. Teramoto Osaka City University, Osaka 558-8585    K. Trabelsi High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    M. Uchida Tokyo Institute of Technology, Tokyo 152-8550    S. Uehara High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    T. Uglov Institute for Theoretical and Experimental Physics, Moscow 117218 Moscow Institute of Physics and Technology, Moscow Region 141700    S. Uno High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193    P. Urquijo School of Physics, University of Melbourne, Victoria 3010    Y. Usov Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    C. Van Hulse University of the Basque Country UPV/EHU, 48080 Bilbao    P. Vanhoefer Max-Planck-Institut für Physik, 80805 München    G. Varner University of Hawaii, Honolulu, Hawaii 96822    A. Vinokurova Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    A. Vossen Indiana University, Bloomington, Indiana 47408    M. N. Wagner Justus-Liebig-Universität Gießen, 35392 Gießen    C. H. Wang National United University, Miao Li 36003    M.-Z. Wang Department of Physics, National Taiwan University, Taipei 10617    X. L. Wang CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061    M. Watanabe Niigata University, Niigata 950-2181    Y. Watanabe Kanagawa University, Yokohama 221-8686    K. M. Williams CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061    E. Won Korea University, Seoul 136-713    J. Yamaoka Pacific Northwest National Laboratory, Richland, Washington 99352    J. Yelton University of Florida, Gainesville, Florida 32611    C. Z. Yuan Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049    Y. Yusa Niigata University, Niigata 950-2181    Z. P. Zhang University of Science and Technology of China, Hefei 230026    V. Zhilich Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    V. Zhulanov Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 Novosibirsk State University, Novosibirsk 630090    A. Zupanc Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana J. Stefan Institute, 1000 Ljubljana
Abstract

We measure the decay using data collected at the resonance with the Belle detector at the KEKB collider. The data sample used corresponds to an integrated luminosity of 121.4 fb. We measure a branching fraction with a significance of 5.1 standard deviations. This measurement constitutes the first observation of this decay.

pacs:
13.25.Hw, 14.40.Nd

The Belle Collaboration

The two-body decays , where is either a pion or kaon, have now all been observed PDG (). In contrast, the neutral-daughter decays have yet to be observed. The decay  charge-conjugate () is of particular interest because the branching fraction is predicted to be relatively large. In the standard model (SM), the decay proceeds mainly via a loop (or “penguin”) transition as shown in Fig. 1, and the branching fraction is predicted to be in the range  SM-branching (). The presence of non-SM particles or couplings could enhance this value Chang:2013hba (). It has been pointed out that asymmetries in decays are promising observables in which to search for new physics susy ().

Figure 1: Loop diagram for decays.

The current upper limit on the branching fraction, at 90% confidence level, was set by the Belle Collaboration using of data recorded at the resonance Peng:2010ze (). Here, we update this result using the full data set of recorded at the . The analysis presented here uses improved tracking, reconstruction, and continuum suppression algorithms. The data set corresponds to pairs Oswald:2015dma () produced in three decay channels: , or , and . The latter two channels dominate, with production fractions of and Esen:2012yz ().

The Belle detector is a large-solid-angle magnetic spectrometer consisting of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov counters, a barrel-like arrangement of time-of-flight scintillation counters, and an electromagnetic calorimeter comprising CsI(Tl) crystals. These detector components are located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside the coil is instrumented to detect mesons and to identify muons. The detector is described in detail elsewhere Belle (); svd2 (). The origin of the coordinate system is defined as the position of the nominal interaction point (IP). The axis is aligned with the direction opposite the beam and is parallel to the direction of the magnetic field within the solenoid. The axis is horizontal and points towards the outside of the storage ring; the axis points vertically upward.

Candidate mesons are reconstructed via the decay using a neural network (NN) technique Feindt:2006pm (). The NN uses the following information: the momentum in the laboratory frame; the distance along between the two track helices at their closest approach; the flight length in the - plane; the angle between the momentum and the vector joining the decay vertex to the IP; the angle between the pion momentum and the laboratory-frame direction in the rest frame; the distance-of-closest-approach in the - plane between the IP and the two pion helices; and the pion hit information in the SVD and CDC. The selection efficiency is 87% over the momentum range of interest. We also require that the invariant mass be within 12 MeV/ (about 3.5 in resolution) of the nominal mass PDG ().

To identify candidates, we define two variables: the beam-energy-constrained mass ; and the energy difference , where is the beam energy and and are the energy and momentum, respectively, of the candidate. These quantities are evaluated in the center-of-mass frame. We require that events satisfy  GeV/ and .

To suppress background arising from continuum production, we use a second NN Feindt:2006pm () that distinguishes jetlike continuum events from more spherical events. This NN uses the following input variables, which characterize the event topology: the cosine of the angle between the thrust axis Brandt:1964sa () of the candidate and the thrust axis of the rest of the event; the cosine of the angle between the thrust axis and the axis; a set of 16 modified Fox-Wolfram moments SFW (); and the ratio of the second to zeroth (unmodified) Fox-Wolfram moments. All quantities are evaluated in the center-of-mass frame. The NN is trained using Monte Carlo (MC) simulated signal events and background events. The MC samples are obtained using EvtGen Lange:2001uf () for event generation and Geant3 geant3 () for modeling the detector response. The NN has a single output variable () that ranges from for backgroundlike events to for signal-like events. We require , which rejects approximately 85% of background while retaining 83% of signal decays. We subsequently translate to a new variable {linenomath}

(1)

where and is the maximum value of obtained from a large sample of signal MC decays. The distribution of is well modeled by a Gaussian function.

After applying all selection criteria, approximately 1.0% of events have multiple candidates. For these events, we retain the candidate having the smallest value of obtained from the deviations of the reconstructed masses from their nominal values PDG (). According to MC simulation, this criterion selects the correct candidate % of the time.

We measure the signal yield by performing an unbinned extended maximum likelihood fit to the variables , , and . The likelihood function is defined as {linenomath}

(2)

where is the yield of component ; is the probability density function (PDF) of component for event ; runs over the two event categories (signal and background); and runs over all events in the sample (). Backgrounds arising from other and non- decays were studied using MC simulation and found to be negligible. As correlations among the variables , , and are found to be small, the three-dimensional PDFs are factorized into the product of separate one-dimensional PDFs.

The signal PDF is defined as {linenomath}

where , , and are the PDFs for signal arising from , and decays. The and PDFs are modeled with Gaussian functions, and the PDFs are each modeled with a sum of two Gaussian functions having a common mean. All parameters of the signal PDF are fixed to the corresponding MC values. The peak positions for and are adjusted according to small data-MC differences observed in a control sample of decays Esen:2012yz (). As this control sample has only modest statistics, the resolutions for , , and , and the peak position for , are adjusted for data-MC differences using a high statistics sample of decays. For background, the , , and PDFs are modeled with an ARGUS function Albrecht:1990am (), a first-order Chebyshev polynomial, and a Gaussian function, respectively. All parameters of the background PDFs except for the end point of the ARGUS function are floated in the fit.

The results of the fit are signal events and continuum background events. Projections of the fit are shown in Fig. 2. The branching fraction is calculated via {linenomath}

(4)

where is the fitted signal yield; is the number of events; is the branching fraction for  PDG (); and is the signal efficiency as determined from MC simulation. The efficiency is corrected by a factor for each reconstructed , to account for a small difference in reconstruction efficiency between data and simulation. This correction is estimated from a high statistics sample of decays. The factor 0.50 accounts for the 50% probability for (since is even). Inserting these values gives , where the error is statistical.

Figure 2: Projections of the 3D fit to the real data: (a) in and ; (b) in and ; and (c) in and . The points with error bars are data, the (green) dashed curves show the signal, (magenta) dotted curves show the continuum background, and (blue) solid curves show the total. The values of these fit projections are 0.30, 0.43, and 0.26, respectively, which indicate that the fit gives a good description of the data. The three peaks in arise from , and decays.

The systematic uncertainty on arises from several sources, as listed in Table 1. The uncertainties due to the fixed parameters in the PDF shape are estimated by varying the parameters individually according to their statistical uncertainties. For each variation the branching fraction is recalculated, and the difference with the nominal branching fraction is taken as the systematic uncertainty associated with that parameter. We add together all uncertainties in quadrature to obtain the overall uncertainty due to fixed parameters. The uncertainties due to errors in the calibration factors and the fractions are evaluated in a similar manner. To test the stability of our fitting procedure, we generate and fit a large ensemble of MC pseudoexperiments. By comparing the mean of the fitted yields with the input value, a bias of is found. We attribute this bias to our neglecting small correlations among the fitted observables. An 0.9% systematic uncertainty is assigned due to the selection; this is obtained by comparing the selection efficiencies in MC simulationand data for the control sample. We assign a 2.0% systematic uncertainty for each reconstructed ; this is determined using a sample. The uncertainty on due to the MC sample size is 0.2%. The total of the above systematic uncertainties is calculated as their sum in quadrature. In addition, there is a 10.1% uncertainty due to the number of pairs. As this large uncertainty does not arise from our analysis, we quote it separately.

Source Uncertainty (%)
PDF parametrization 0.2
Calibration factor
Fit bias
reconstruction 4.0
selection 0.9
MC sample size 0.2
0.1
Total (without )
10.1
Table 1: Systematic uncertainties on . Those listed in the upper section are associated with fitting for the signal yields and are included in the signal significance.

The signal significance is calculated as , where is the likelihood value when the signal yield is fixed to zero, and is the likelihood value of the nominal fit. We include systematic uncertainties in the significance by convolving the likelihood function with a Gaussian function whose width is equal to that part of the systematic uncertainty that affects the signal yield. We obtain a signal significance of 5.1 standard deviations; thus, our measurement constitutes the first observation of this decay.

Figure 3: Background subtracted sPlot distributions of for the (a) higher momentum and (b) lower momentum candidates.

Figure 3 shows the background-subtracted sPlot splot () distributions of , where the selection is removed for the pair being plotted. No contribution is observed. We check this quantitatively by performing our signal fit for events in the mass sidebands of each [ and ]. The extracted signal yields, and for the higher momentum and lower momentum , respectively, are consistent with zero. We calculate the expected number of events in our signal sample using MC simulation and the measured branching fraction,  Aaij:2013uta (); the result is 0.001.

In summary, we report the first observation of the decay . The branching fraction is measured to be {linenomath}

where the first uncertainty is statistical, the second is systematic, and the third reflects the uncertainty due to the total number of pairs. This value is in good agreement with the SM predictions SM-branching (), and it implies that the Belle II experiment Abe:2010gxa () will reconstruct over 1000 of these decays. Such a sample would allow for a much higher sensitivity search for new physics in this penguin-dominated decay.

ACKNOWLEDGMENTS

We thank the KEKB group for the excellent operation of the accelerator; the KEK cryogenics group for the efficient operation of the solenoid; and the KEK computer group, the National Institute of Informatics, and the PNNL/EMSL computing group for valuable computing and SINET4 network support. We acknowledge support from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of Japan, the Japan Society for the Promotion of Science (JSPS), and the Tau-Lepton Physics Research Center of Nagoya University; the Australian Research Council; Austrian Science Fund under Grants No. P 22742-N16 and P 26794-N20; the National Natural Science Foundation of China under Contracts No. 10575109, No. 10775142, No. 10875115, No. 11175187, and No. 11475187; the Chinese Academy of Science Center for Excellence in Particle Physics; the Ministry of Education, Youth and Sports of the Czech Republic under Contract No. LG14034; the Carl Zeiss Foundation, the Deutsche Forschungsgemeinschaft and the VolkswagenStiftung; the Department of Science and Technology of India; the Istituto Nazionale di Fisica Nucleare of Italy; the WCU program of the Ministry of Education, National Research Foundation (NRF) of Korea Grants No. 2011-0029457, No. 2012-0008143, No. 2012R1A1A2008330, No. 2013R1A1A3007772, No. 2014R1A2A2A01005286, No. 2014R1A2A2A01002734, No. 2015R1A2A2A01003280 , No. 2015H1A2A1033649; the Basic Research Lab program under NRF Grant No. KRF-2011-0020333, Center for Korean J-PARC Users, No. NRF-2013K1A3A7A06056592; the Brain Korea 21-Plus program and Radiation Science Research Institute; the Polish Ministry of Science and Higher Education and the National Science Center; the Ministry of Education and Science of the Russian Federation and the Russian Foundation for Basic Research; the Slovenian Research Agency; the Basque Foundation for Science (IKERBASQUE) and the Euskal Herriko Unibertsitatea (UPV/EHU) under program UFI 11/55 (Spain); the Swiss National Science Foundation; the National Science Council and the Ministry of Education of Taiwan; and the U.S. Department of Energy and the National Science Foundation. This work is supported by a Grant-in-Aid from MEXT for Science Research in a Priority Area (“New Development of Flavor Physics”) and from JSPS for Creative Scientific Research (“Evolution of Tau-lepton Physics”).

References

  • (1) K. A. Olive et al. (Particle Data Group), “Review of particle physics,” Chin. Phys. C 38, 090001 (2014).
  • (2) Unless stated otherwise, charge-conjugate modes are implicitly included.
  • (3) C. H. Chen, “Analysis of decays in the PQCD,” Phys. Lett. B 520, 33 (2001); A. R. Williamson and J. Zupan, “Two body decays with isosinglet final states in SCET,” Phys. Rev. D 74, 014003 (2006); A. Ali, G. Kramer, Y. Li, C. D. Lu, Y. L. Shen, W. Wang and Y. M. Wang, “Charmless non-leptonic decays to , and final states in the pQCD approach,” Phys. Rev. D 76, 074018 (2007); C. K. Chua, “Rescattering effects in charmless decays,” Phys. Rev. D 78, 076002 (2008); K. Wang and G. Zhu, “Flavor dependence of annihilation parameters in QCD factorization,” Phys. Rev. D 88, 014043 (2013); J. J. Wang, D. T. Lin, W. Sun, Z. J. Ji, S. Cheng and Z. J. Xiao, “ decays and effects of the next-to-leading order contribution,” Phys. Rev. D 89, 074046 (2014); Q. Chang, J. Sun, Y. Yang and X. Li, “A combined fit on the annihilation corrections in decays within QCDF,” Phys. Lett. B 740, 56 (2015); H. Y. Cheng, C. W. Chiang and A. L. Kuo, “Updating decays in the framework of flavor symmetry,” Phys. Rev. D 91, 014011 (2015).
  • (4) Q. Chang, X. Q. Li and Y. D. Yang, “A comprehensive analysis of hadronic transitions in a family non-universal model,” J. Phys. G 41, 105002 (2014).
  • (5) S. Baek, D. London, J. Matias and J. Virto, “ and decays within supersymmetry,” J. High Energy Phys. 12, (2006) 019; A. Hayakawa, Y. Shimizu, M. Tanimoto and K. Yamamoto, “Searching for the squark flavor mixing in violations of and decays,” Prog. Theor. Exp. Phys. 2014, 023B04 (2014).
  • (6) C.-C. Peng et al. (Belle Collaboration), “Search for decays at the resonance,” Phys. Rev. D 82, 072007 (2010).
  • (7) C. Oswald et al. (Belle Collaboration), “Semi-inclusive studies of semileptonic decays at Belle,” Phys. Rev. D 92, 072013 (2015)
  • (8) S. Esen et al. (Belle Collaboration), “Precise measurement of the branching fractions for and first measurement of the polarization using collisions,” Phys. Rev. D 87, 031101(R) (2013).
  • (9) A. Abashian et al. (Belle Collaboration), “The Belle detector,” Nucl. Instrum. Methods Phys. Res., Sect. A 479, 117 (2002); also see the detector section in J.Brodzicka et al., “Physics achievements from the Belle experiment,” Prog. Theor. Exp. Phys. 2012, 04D001 (2012).
  • (10) Z. Natkaniec et al. (Belle SVD2 Group), “Status of the Belle silicon vertex detector,” Nucl. Instrum. Methods Phys. Res., Sect. A 560, 1 (2006).
  • (11) M. Feindt and U. Kerzel, “The NEUROBAYES neural network package,” Nucl. Instrum. Methods Phys. Res., Sect. A 559, 190 (2006).
  • (12) S. Brandt, C. Peyrou, R. Sosnowski and A. Wroblewski, “The principal axis of jets. An attempt to analyze high-energy collisions as two-body processes,” Phys. Lett. 12, 57 (1964).
  • (13) G. C. Fox and S. Wolfram, “Observables for the Analysis of Event Shapes in Annihilation and Other Processes,” Phys. Rev. Lett. 41, 1581 (1978); The modified moments used in this Letter are described in S. H. Lee et al. (Belle Collaboration), “Evidence for ,” Phys. Rev. Lett. 91, 261801 (2003).
  • (14) D. J. Lange, “The EvtGen particle decay simulation package,” Nucl. Instrum. Methods Phys. Res., Sect. A 462, 152 (2001).
  • (15) R. Brun et al., GEANT 3.21, CERN Report DD/EE/84-1, 1984.
  • (16) H. Albrecht et al. (ARGUS Collaboration), “Search for hadronic decays,” Phys. Lett. B 241, 278 (1990).
  • (17) M. Pivk and F. R. Le Diberder, “sPlot: A statistical tool to unfold data distributions”, Nucl. Instrum. Methods Phys. Res., Sect. A 555, 356 (2005).
  • (18) R. Aaij et al. (LHCb Collaboration), “Study of decays with first observation of and ,” J. High Energy Phys. 10, (2013) 143.
  • (19) T. Abe et al. (Belle II Collaboration), “Belle II Technical Design Report,” arXiv:1011.0352.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
221183
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description