Introduction

1. Introduction A major issue for the charmonium particle system are the two c-cbar states below open charm threshold that are still not well established: the n=2singlet S state, the hc(2S), and the n=1 singlet P state, the hc(1P); see Fig. 1. The observation of these states and the determination of their masses would provide useful information about the spin-spin part of the c cbar potential. It is expected [1] that the (2S)-c(2S) mass splitting is smaller than that for the n=1c-J/ system by factors corresponding the the reduced value of the wave-function at zero quark separations and the running of the strong coupling constants(Q2) between Q = M(J/) and M((2S)). The expected mass range for the c(2S) is 3625 < M(c(2S) )<3645 MeV/c2. The Crystal Ball group [2] identified an excess of E = 91 MeV -rays in the reaction (2S)X as evidence for the c(2S) with mass M(c(2S))  3594 MeV/c2. This is well below potential model expectations and has not been confirmed by other experiments [3]. Signal Extraction Figures 3(a) through (y) show the Mbc projections for events with |E|<40 MeV for 40 MeV/c2 bins of M(KsK) with central values from 2840 thru 3800 MeV/c2. An c signal is apparent around M(KsK) = 2980 MeV/c2. A smaller, but distinct signal is evident near 3660 MeV/c2, in the expected c(2S). mass region.We perform simultaneous fits to each Mbc plot and its corresponding E distribution (for events in the same KsK mass bin and in the 5.27<Mbc<5.29 GeV/c2 signal band) simultaneously. Observation of the hc(2S) in exclusive BK KSK-p+ decaysS.K. Choi, Gyeongsang Univ. & S.L. Olsen, Univ. of Hawaii(For the Belle Collaboration) 1. The charmonium level diagram for ccbar bound states below the MD “open charm” threshold (from the PDG tables [4]). All of the levels in this diagram are experimentally well established except for the c(2S) and the (unshown) hc(1P1) state. 3. The Mbc distributions for 40 MeV/c2M(KsK) bins. We use Gaussian functions with MC-determined widths for the Mbc and E signals; the areas of the Mbc and E signals are constrained to be equal. For the Mbc background, we use a smooth function with a phase-space-like shape near the kinematic endpoint; for E we use a 2nd-order polynomial. As an example, Figs. 4(a) & (b)shows the fit to the 3640 MeV/c2KsK mass bin. B-meson decays to final states containing charmonium particles are relatively common and, thus, the large numbers of B meson decays recorded in the Belle detector provide an excellent opportunity to search for exclusive decays of the type BKc(2S). In this experiment we concentrate on the c(2S)KSK+p- decay mode. This a strong decay mode for the c (Br1.8%), has low combinatoric background, is strongly suppressed in (2S) decays, and, since the final state is all charged particles, has good energy resolution in the Belle detector. 4. Simultaneous fits to the Mbc(a) & E (b) plots for the M(KsK)=3640 MeV/c2 bin. Results Figure 5 shows the results of the fits vs M(KsK). In addition to the c signal, the J/ and, possibly, non-resonant BKKsK decays, there is a narrow peak at the high mass side that we identify as the c(2S).The curve in the figure is the result of a fit that includes Breit Wigner shapes to represent the c and c(2S), a Gaussian with a negligibly small  and mass fixed at M(J/), and a 2nd-order polynomial to represent the non-resonant contribution. The functions are convolved with a Gaussian resolution function with a MC-determined width of 15 MeV/c2. The fit results are: The Belle Detector The Belle detector [5] at the KEKB e+e storage ring [6] is a nearly 4 magnetic spectrometer that consists of a three-layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), aerogel Cherenkov counters (ACC), time-of-flight scintillators (TOF) and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a 1.5 T axial magnetic field. The iron flux-return is instrumented to detect KL mesons and identify muons (KLM). cc(2S) Nevts: 104±14 39±11 Mass (MeV /c2) 2979±2 3654±6 Width (MeV/c2) 11±11 15±24, where only the statistical errors from the fit are listed. The c mass is in good agreement with the world average value of 2979.8±1.8 MeV/c2. The c width value agrees, within errors, both Fig. 5 with the PDG [4] value of 13.2±3.5MeV/c2 and a new CLEO result of 26±6 MeV/c2. The sum of observed events in the three bins centered around the hc(2S) is 56, while the integral of the 2nd-order polynomial over the same interval is 21±2. The probability for the latter number of events to fluctuate up to 56 is ~10, which corresponds to a signal significance of greater that 6. The fitted mass of the candidate hc(2S) is substantially above that of the Crystal Ball value and consistent, within errors, with the upper end of potential model expectations. We repeated the analysis with different bin sizes and bin centers shifted by half a bin width and take the maximum variation in mass, namely 8 MeV/c2, as the systematic error. Changes in the resolution width, the order of the polynomial function, and the inclusion of a possible cc1 contribution have negligible effects on the results. The limited statistics and wide bin size preclude an accurate determination of the width. We set a 90% CL upper limit of (c(2S))<55 MeV/c2, from the worst-case result of all the binning combinations that were tried. A MC simulation indicates that the experimental acceptance for the c and c(2S) are very nearly the same. Thus the ratio of product branching fractions for the two states is the ratio of observed event yields. Event Selection We select events with KKsK or KsKsK (inclusion of charge-conjugate modes is implicit). Charged kaons and pions are required to originate near the IP, and identified as such by the combined ACC, TOF and CDC dE/dx measurements. Candidate Ksp+p- decays arep+p- pairs from a common vertex that is displaced from the IP and with an invariant mass |M(p+p-) -–M(K0)|<12 MeV/c2 (3). We identify BKKsK signals by their beam-constrained mass: Mbc = [(Ecm/2)2 – p2]1/2 and cm energy difference: E = E - Ecm/2, where p and E are the total momentum and energy of the KKsK system. 2a 2c 2b Summary In a sample of exclusive BKKsK decays, we observe a narrow peak in the KsK mass spectrum near M((2S)). The significance of the signal is greater than 6. The observed features of this peak lead us to conclude that this is most likely the c(2S).The measured parameters are: c(2S) properties M = 3654 ± 6 (stat) ± 8 (syst) MeV/c2  < 55 MeV/c2 (90% CL) r (KsK) = 0.38 ± 0.12 (stat) ± 0.5 (syst), where r(KsK) is the ratio of c(2S) to c(1S) product branching fractions: r(KsK)=Br(B Kc(2S))Br(c(2S)) KsK). Br(B Kc) Br(cKsK). Acknowledgement This research was made possible by the strong and dedicated efforts of the KEKB team and the extraordinary performance of the KEKB storage ring. • References • See, eg, W.Buchmuller and S-H.H. Tye, Phys. Rev. D24, 132 (1981) • C.Edwards et al. (Crystal Ball), Phys.Rev. Lett. 48, 70 (1982). • T.A.Armstrong et al. (E760), Phys. Rev. D52, 4839 (1985). M. Masuzawa, Northwestern Univ. PhD thesis (1992), unpublished. • D.E.Groom et al. (PDG), Eur. Phys. Journ. C15, 1 (2002) • A.Abashian et al. (Belle), Nucl. Instr. & Meth A479, 117 (2002) • E.Kikutani et al. (KEKB), KEK preprint 2001-157; submitted to Nucl. Instr. & Meth. Background Reduction To suppress e+eq qbar continuum events (q = u, d, s & c), we use a Fisher discriminant of event-shape variables (Fig. 2a) and the cosine of the angle between the B candidate direction and the beam axis (Fig. 2b). Monte Carlo is used to determine signal (Lsig, solid histograms) and continuum (Lcont, dashed histograms) likelihood functions. We select events with Lsig/(Lsig + Lcont)>0.6, a value that optimizes the BKc; cKSK+p- signal. Backgrounds from BD(Ds)X decays are reduced by vetoing candidates with |M(K)-M(D)| or |M(KsK)-M(Ds)|<10 MeV/c2. Also, since the decays c(nS)K*K are suppressed by a momentum barrier, candidates with |M(K)-M(K*)|<50 MeV/c2 are also rejected.