arxiv: v2 [hep-ex] 11 Feb 2015

Transcription

1 A Comprehensive Study of the Radiative Decays of J/ψ and ψ(s) to Pseudoscalar Meson Pairs, and Search for Glueballs S. Dobbs, A. Tomaradze, T. Xiao, and Kamal K. Seth Northwestern University, Evanston, Illinois, USA (Dated: September, 7) arxiv:.v [hepex] Feb Using pb of e + e annihilation data taken at s =. GeV, a comprehensive study has been made of the radiative decays of samples of. million J/ψ and. million ψ(s) into pairs of pseudoscalar mesons, π + π, π π, K + K, K SK S, and ηη. Product branching fractions for the radiative decays of J/ψ and ψ(s) to scalar resonances f (7,,7,, and ), and tensor resonances f (7,, and ) have been determined, and are discussed in relation to predicted glueballs. For ψ(s) radiative decays the search for glueballs has been extended to masses between. GeV and. GeV. I. INTRODUCTION Quantum Chromodynamics(QCD) is presently the accepted theory of strong interactions. Its predictions for hadronic phenomena at high energies have received remarkable experimental confirmation, and it is an integral part of the Standard Model of particle physics. In contrast, at lower energies the success of QCD has been comparatively limited. This is mainly due to the fact that while at high energies reliable predictions can be made in perturbative QCD (pqcd) models, at low energies pqcd predictions are unreliable, and lattice calculations with the required precision are in many cases not yet available. In 97, Fritzsch and GellMann [] pointed out that if the quarkgluon field theory indeed yields a correct description of strong interactions, there must exist glue states in the hadron spectrum, and Fritzsch and Minkowski [] presented a detailed discussion of the phenomenology of the spectrum of what they called gluestates. The gluestates, which we here call glueballs, are bosons with baryon number and isospin equal to zero, and they are SU() flavor singlets to the extent to which SU() breaking effect can be neglected []. Ever since 9, when η() was claimed as the first glueball discovered [], numerous experimental searches for glueballs have been reported, and numerous claims and counterclaims abound. Many reviews of QCD exotics, including glueballs, have been published. Two of the latest comprehensive reviews are by Klempt and Zaitsev [] and Ochs []. Fritzsch and Minkowski [] have detailed colorsinglet states of two and three gluons allowed by conservation rules, which are identical to those for states of two and three photons. The twogluon states have charge conjugation, C = +, while the three gluon states can have C = + or. The lowest mass states have S waves between the gluons, and are twogluon states with J PC = ++, ++, and +, and threegluon states with J PC = ++, +,, and. Predictions for glueball spectra have been made in bag models [], QCDbased potential models [7], QCD sum rules [], and lattice gauge calculations [9 ]. Lattice gauge calculations have supplanted most other predictions. The lattice calculations of Refs.,, and were done in the quenched approximation with no quarks, and their predictions are in general agreement with each other. The first unquenched lattice calculations with + flavor seaquarks have been reported by the UKQCD Collaboration []. They obtain results for the masses of the ++, +, and ++ glueballs which are in good agreement with the lattice predictions of Refs. and in the quenched approximation, and conclude that they find no evidence of strong unquenching effects. Thus, the results of Ref., reproduced in Fig., can be considered as representative of lattice predictions so far. The order of the lowest lying glueballs is as anticipated by Fritzsch and Minkowski [], with the prediction for the first scalar ( ++ ) and tensor ( ++ ) glueball masses are G( ++ ) = 7()() MeV, G( ++ ) = 9()() MeV. These are followed by five heavier glueballs with masses less than M(ψ(S)). These are + (()() MeV), + (9()() MeV), + (()()MeV), + (()(7)MeV), and ++ (7()()MeV). Inthepresentpaper,ourglueball search is limited to masses less than. GeV, below the masses of the χ cj resonances of charmonium. Early searches for glueballs were misled by the wishful expectation of pure glueballs. With the realization that the predicted glueballs have masses close to those of normal q q mesons with the same quantum numbers with which they must mix, the more realistic searches aspire to identify states which contain substantial glueball content in their wave functions. Searches for such states have been based mainly on three expectations. These are: (a) for a given J PC, the glueball state is supernumary, i.e. a state in excess of the expected number of q q states, e.g., the nonet of ++ SU() states, (b) a pure glueball state does not decay preferentially to states containing up, down, or strange quarks, i.e., pure glueball decays to pseudoscalar pairs should have equal branching fractions for each available charged state. Thus decays to ππ, KK, ηη, η η and ηη, should have branching fractions proportional to : : : :, and (c) certain decays are particularly favorable for populating glueballs, e.g.,

2 p p annihilation, and radiative decays of vector states of quarkonia, J/ψ,ψ,Υ, to γgg..., where g stands for gluons. Numerous experimental attempts have been made to exploit these expectations, but they have been generally limited in their scope. An extensive review of searches for glueballs in p p annihilations, central production, and hadronic and radiative decays of vector states of charmonium and bottomonium has been made by Klempt and Zaitsev []. Radiative decays of J/ψ have been studied by Crystal Ball [], Mark III [, ], DM [7], and BES [ ]. Among the latest of these are BES measurements of J/ψ γ(π + π,π π ) [] and J/ψ γ(k + K,K S K S ) [] with million directly produced J/ψ in e + e J/ψ. The invariant mass spectra in the decays of these directly produced J/ψ contain large backgrounds due to ISR production of π + π π and KKπ with one missing photon which simulate γπ + π and γk + K decays. More recently, BES III has reported a study of J/ψ γηη with. directly produced J/ψ []. No studies of radiative decays of ψ(s) have been published, and only an unpublished report with million ψ(s) by BES exists []. Inthepresentpaperwepresentresultsofastudyofthe radiative decays of the charmonium resonances,. million J/ψ, and. million ψ(s) resonances populated in e + e annihilation, and decaying into pairs of pseudoscalar, ππ, KK, and ηη. We determine the masses, widths, and product branching fractions of the resonances observed, and discuss their possible relations to glueballs. Since the J/ψ used in the present investigation are produced by the exclusive decay ψ(s) π + π J/ψ, the invariant mass spectra we observe from the J/ψ radiative decays are free from the ISRproduced backgrounds. The invariant mass spectra for ψ(s) decays contain, of course, background contribution from ISR. The paper is organized as follows. In Sec. II we describe the event selection criteria. In Sec. III we describe Monte Carlo event simulations and use them to determine instrumental resolutions and efficiencies. In Sec. IV we present the Dalitz plots for J/ψ and ψ(s) radiative decays to pseudoscalar pairs. In Sec. V.A we present the ππ and KK invariant mass distributions for the radiative decays of J/ψ, and discuss the resonances observed in them. In Sec. V.B we present the similar mass distributions and discussion for the radiative decays of ψ(s). In Sec. V.C we present the results of our search for the putative tensor glueball ξ(). In Sec. VI we summarize and compare our results for J/ψ and ψ(s) radiative decays to pseudoscalar pairs for the observed resonances with masses below. GeV. In Sec. VII.A we present the invariant mass spectra for ψ(s) radiative decays to ππ, KK, and ηη decays in search for potential glueball candidates in the mass range.. GeV, and establish upper limits for any resonance enhancements. We observestrong decaysof χ c ( ++ ) and χ c ( ++ ) resonances of charmonium in our data for ψ(s) radiative decays. Because the decays of these resonances are necessarily r M G M G (GeV) FIG.. Quenched lattice predictions of glueball spectrum by Chen et al. []. The states are grouped by their parity and charge conjugation indicated in the abscissa. TABLE I. Summary of the requirements for the number of reconstructed particles for the reactions investigated in this paper. Decay mode # tracks # showers Other ψ(s) γπ + π γk + K γπ π π γηη η γksk S KS J/ψ γπ + π γk + K γπ π π γηη η γksk S KS mediated by two gluons, just as the decays of twogluon glueballs, it is instructive to examine and compare their characteristics. We do so in Sec. VII.B in terms of the branching fractions measured using the same data as we use in the present paper, and published as Ref. []. In Sec. VIII we discuss systematic uncertainties in our results. In Sec. IX we summarize the results of the present investigation for the decays of known and proposed resonances into pseudoscalar pairs, and their implications for scalar and tensor glueball candidates.

3 Counts / MeV J/ψ γ+(π,η) γ π η DATA Counts / MeV K S J/ψ γ+k S K S DATA M(γγ) (GeV) M(π + π ) (MeV) Counts / MeV ψ(s) γ+(π,η) γ π η DATA Counts / MeV K S ψ(s) γ+k S K S DATA M(γγ) (GeV) M(π + π ) (MeV) FIG.. Two photon invariant mass distributions for radiative decays of (top panel) J/ψ and (bottom panel) ψ(s), with at least photons in the final state. In both panels, the top histograms show invariant mass distribution for all two photon combinations in the selected events. The middle histograms show distributions for all two photon combinations in the events that pass the C energymomentum kinematic fit with χ C <. The shaded solid histograms show the distributions of events which contain two π or two η candidates which best fit π and η masses. Signals from π γγ and η γγ decays are clearly visible in all three histograms. FIG.. Two pion invariant mass distributions for radiative decays of (top panel) J/ψ and (bottom panel) ψ(s) with at least two K S π + π candidates reconstructed. In both panels the top histograms show all K S π + π candidates in the selected events. The dashed histograms show distributions for all K S π + π candidates in the events that pass the C energymomentum kinematic fit with χ C <. The shaded histograms show the distribution for the best two K S candidates in the event, as determined by the kinematic fit. The mass peaks corresponding to the selected K S candidates are clearly visible in all three histograms. II. EVENT SELECTIONS We use pb of e + e annihilation data taken at s =. GeV at the CESR e + e collider at Cornell University. It corresponds to. million ψ(s) produced, and. million J/ψ tagged by the decay ψ(s) π + π J/ψ. We also use pb of e + e annihilation data taken at s =.7 GeV to study backgrounds from offresonance decays. The decay products were analyzed by the CLEOc detector, which has been described elsewhere in detail []. Briefly, it consists of a CsI(Tl) electromagnetic calorimeter, an inner vertex drift chamber, a central drift chamber, and a ring imaging Cherenkov (RICH) detector, all inside a superconducting solenoid magnet providing a nominal. Tesla magnetic field. The acceptance for charged and neutral particles is cos θ <.9. Charged particle resolution is σ p /p GeV/c. Photon resolution is σ E /E GeV, and MeV. The detector response was studied using a GEANTbased Monte Carlo (MC) simulation including radiation corrections []. The particle multiplicity requirements for the various decay modes are given in Table I. Charged particle tracks were accepted in the full range of the drift chamber s acceptance, cosθ <.9 (θ is the polar angle with respect to the positron beam), and were selected by standard CLEO quality requirements. The total charge in the event was required to be zero. Photon candidates were accepted in the Good Barrel and Good Endcap regions of the calorimeter ( cosθ <. and cosθ =..9, respectively), and were required to contain none of the few noisy calorimeter cells, and to have a transverse energy deposition that was consistent with that of an electromagnetic shower. The π and η were reconstructed in their decay to two photons. Pairs of photon candidates were kinematically

4 Counts /. MeV (a) Counts /. J/ψ γ f (7), f (7) π + π Recoil Mass against π + π (MeV) cos θ (f ) Counts /. MeV (b) Counts / J/ψ γ f (7), f (7) K K Recoil Mass against π + π (MeV) cos θ (f ) FIG.. Mass distribution of events recoiling against π + π in ψ(s) π + π J/ψ decay (a) before and (b) after C kinematic fit. The fitted mass of the J/ψ peak is found to be M(J/ψ) = 97. ±.(stat) MeV. FIG.. Distributions of cosθ(f,) for f (7) π + π and f (7) K + K decays. The point show the distributions from data, and the shaded histograms show the results of MC simulations. fitted to the known π and η masses, and the pairs were accepted if they had a mass σ of the particle mass. Invariant mass distributions at different stages of these selections are shown in Fig.. The KS were reconstructed in their decay to two charged pions. Pairs of oppositely charged pions were kinematically fitted to originate from a common vertex, and this common vertex was required to be displaced by σ from the e + e interaction point. Invariant mass distributions at different stages of these selections are shown in Fig.. Charged pions and kaons were identified using information from both de/dx and RICH measurements. For the pion and kaon tracks with momenta > MeV, the variable L i,j was constructed, combining information from de/dx measurements, σ de/dx i,j, and the RICH loglikelihood, L RICH i,j, for the particle hypotheses i and j L i,j = (σ de/dx i ) (σ de/dx j ) +L RICH i L RICH j To identify pions, it was required that L π,k <. For kaons, it was required that L π,k >. For pion and kaon tracks with momenta < MeV, for which RICH information was not available, it was required that the de/dx of the track be consistent within σ of the respective mass hypothesis, and that it be more consistent with its mass hypothesis than with the alternate mass hypothesis. For the γπ + π final state, electron contamination was removed by rejecting events with E CC /p >., and muon contamination was removed by rejecting events whichhavetracksmatchedtoanyhits inthe muonchambers. Since we study resonances with masses <. GeV for J/ψ decays and <. GeV for ψ(s) decays, we safely require E γ > MeV for the radiative photon candidate. We obtain our J/ψ sample from ψ(s) by tagging events with two oppositely charged tracks, which were assumed to be pions, to have a recoil mass in the range M(J/ψ) ± MeV. The mass distribution of events recoiling against π + π in the decay ψ(s) π + π J/ψ is shown in Fig. (a) before C kinematic fit, and Fig. (b) after C fit. The strong peak corresponds to J/ψ events. To reconstructthe full ψ(s) and J/ψ events, we make a C energy/momentum kinematic fit, constraining the sum of the vectorsof all final state particles to be equal to the vector of the e + e interaction. For J/ψ decays, the J/ψ decay products were additionally constrained to have M(J/ψ). In all cases the final fit was required to have χ C <. When constructing the invariant masses from the radiative decays, the vectors resulting from these fits were used.

5 ) γ) (GeV (π M 9 7 J/ψ γ π + π ) γ) (GeV M (π 9 7 J/ψ γ π π 7 9 M (π + γ) (GeV ) M 7 9 (π γ) (GeV ) ) γ) (GeV (K 7 + J/ψ γ K K ) (GeV γ) S, 7 J/ψ γ K K S S M M(K 7 + M (K γ) (GeV ) 7 M(K γ) (GeV ) S, FIG.. Dalitz plots for (top two panels) J/ψ γπ + π and γπ π, and (bottom two panels) J/ψ γk + K and γk SK S. The dashed lines indicate the regions used to reject J/ψ ρπ and K K decays in the γπ + π and γk + K decays, respectively. III. MONTE CARLO SAMPLES AND EFFICIENCIES Monte Carlo samples of simulated events were generated using EvtGen [], and then passed through the standard detector simulation and reconstruction routines, to determine mass resolutions and experimental efficiencies. Monte Carlo samples of k events were generated for each pseudoscalar pair decay chain for the light quark resonances seen in the data and the presumed glueball candidates in the mass region MeV. To investigate background contributions from other decays, such as J/ψ ρπ,k K, Monte Carlo samples of million events were generated. When considering radiative decays to a particular resonance (e.g. ψ γf, ), the expected + αcos θ angular distributions for the corresponding E transitions were used, with α = for decays to a ++ state and α = / for decays to a ++ state. The angular distributions for scalar and tensor mesons decaying into two pseudoscalars were simulated using the helicity amplitude formalism of EvtGen. Fig. illustrates excellent agreement between the data and simulated angular distribution for the wellseparated resonances f (7) and f (7). The efficiencies obtained for decays to individual resonances are shown in Table II. For J/ψ decays to charged final states the efficiencies range from % to %, and for neutral final states they range form % to %. Corresponding efficiencies for ψ(s) decays are approximately a factor two larger. Monte Carlo simulations give instrumental Gaussian resolution widths for masses < MeV as σ = MeV, for π + π, K + K, and KS K S, 9 MeV for π π, and 7 MeV for ηη. IV. DALITZ PLOTS FOR J/ψ AND ψ(s) RADIATIVE DECAYS J/ψ γππ: Dalitz plots for the radiative decays of J/ψ γπ + π andπ π areshowninthetoptwopanelsoffig.. Inthe γπ + π plottherearehorizontalandverticalbandscorresponding to ρ(77). As confirmed by Monte Carlo simu

6 TABLE II. MC determined efficiencies for the reactions ψ(s) γr, R PS,PS, and ψ(s) π + π J/ψ, J/ψ γr, R PS, PS, for selected intermediate resonances R. Resonance Efficiency, (%) X J/ψ γx ψ(s) γx X π + π f (7) 9.. f ()..9 f (7) 7.. f () X π π, π γγ f (7).. f ().. f (7)..7 f () X K + K f (7).. f ().. f (7).. f ().... X KSK S, KS π + π f (7).7. f ().. f (7).. f () X ηη, η γγ f (7) lations, these arise from ρ π decay of J/ψ with one photon from π missing. In subsequent analysis we remove these by rejecting events with M (π ± γ) <. GeV. We do not remove π + π events from ρ decays which populate the diagonal band in the plot in order to be able to take the effect of their tails in subsequent fits to the π + π mass spectrum. Fig. 7(a) shows that the tail of the residual ρ peak does not contribute significantly to the π + π mass spectrum in the M(π + π ) >. GeV. As expected, the ρ(77) contributions are absent in the π π Dalitz plot. The only other feature which is identifiable in these Dalitz plots is the excitation of f (7) as a diagonal band with J = characteristic angular distribution which produces enhancements near the maxima of M (πγ). J/ψ γkk: Dalitz plots for the radiative decays of J/ψ γk + K and KS K S are shown in the middle two panels of Fig.. The bands due to ρ(77) excitation seen in the γπ + π plot are, of course, absent in these plots. Instead, there is evidence for contributions from the decay J/ψ K + K (9) K + K π (+c.c.), with one photon from π missing. We remove these events in the horizontal Counts / MeV Counts / MeV... (a) J/ψ γ π + π DATA J/ψ ρπ MC... M(π + π ) (GeV)... (b) + J/ψ γ K K DATA J/ψ KKπ MC... + M(K K ) (GeV) FIG. 7. (a) Distributions of π + π invariant mass for J/ψ γπ + π decays for data (solid histogram) and J/ψ ρπ MC decays (shaded histogram), illustrating the small residual nonρ background from ρπ decay. (b) Distributions of K + K invariant mass for J/ψ γk + K decays for data (solid histogram) and J/ψ KKπ MC decays (shaded histogram). and vertical bands by rejecting events with M (K ± γ) <. GeV. Monte Carlo simulation, shown in Fig. 7(b), confirms that the residual background, which is too small to be discernible, is completely negligable. Excitation of f () with its characteristic J = angular distribution is identifiable as a diagonal band, followed by the nearly uniform population of f (7). A weakly populated diagonal band corresponding to f () is also discernible. The J/ψ γks K S Dalitz plot is similar to the of J/ψ γk + K plot but has smaller yield. ψ(s) γππ,γkk: The four Dalitz plots forthe ππ and KK decaysofψ(s) are shown in Fig.. As is well known, all hadronic and radiative decays of ψ(s) are suppressed by a factor 7 (the socalled % rule) compared to those of J/ψ. As a result all plots in Fig. have smaller yields than in the corresponding plots in Fig.. Because of the well known anomalously small branching fractions for the decays ψ(s) ρπ and K K, the ρ bands and the K bands seen in J/ψ decays are not seen in the corresponding ψ(s) Dalitz plots. However, there is a strong diag

7 M(πγ ) (GeV) M(πγ ) (GeV) 7 ψ(s) γ π π+ ψ(s) γ K K + M(π+γ ) (GeV) M(πγ ) (GeV) ψ(s) γ K KS S M(K S, γ ) (GeV) M(K γ ) (GeV) ψ(s) γ ππ M(K γ ) (GeV) + M(K γ ) (GeV) S, FIG.. Dalitz plots for (top two panels) ψ(s) γπ + π and γπ π, and (bottom two panels) ψ(s) γk + K and γks KS. onal band for ρ(77) in the ψ(s) γπ + π plot and a narrow band corresponding to φ() in γk + K. These arise from ISR contributions as we show later. The strong diagonal bands due to the population of χ and χ states of charmonium are visible in all four Dalitz plots for ψ(s) radiative decays. V. INVARIANT MASS DISTRIBUTIONS FOR J/ψ AND ψ(s) RADIATIVE DECAYS In order to draw quantitative conclusions from our measurements, we make projections from the Dalitz plots for the invariant masses for the ππ, KK, and ηη final states. An overall presentation of these is made in Fig. 9. Several gross features are apparent in this figure. The first feature is the general similarity in the observed spectra for J/ψ radiative decays (left five panels) and ψ(s) radiative decays (right five panels). The yields for ψ(s) decays are lower. As mentioned earlier, this is mainly the consequence of the wellknown % rule of pqcd which stipulates that the hadronic and radiative decays of ψ(s) have nearly a factor seven smaller widths compared to those of J/ψ, a ratio which follows that for the leptonic decays. The second feature is the expected similarity between the corresponding charged and neutral final states, π + π and π π, and K + K and KS KS. We find that resonance masses determined from charged and neutral spectra are always in agreement within their statistical errors. However, the mass resolution and statistics is better for the charged final states, as expected. This leads to our decision to make resonance mass determinations from the charged decay spectra, and to fix them for the final analysis of the neutral decay spectra. The limited event statistics also does not allow us to make partial wave analysis of the data. One additional observation is that the yields for radiative decays to the ηη final states, though suggestive of the population of f () and f (7), are too small to lead to any meaningful results for the branching fractions, and we do not attempt to analyze the ηη spectra. We analyze the spectra of invariant mass of pseudoscalar pairs for the radiative decays of J/ψ and ψ(s) in terms of noninterfering relativistic BreitWigner resonances. The resonance masses, M, the counts, N, and the background are kept free. It was not found possible

8 Counts/ MeV... J/ψ γ π + π Counts/ MeV... + ψ(s) γ π π... M(π + π ) (GeV)... M(π + π ) (GeV) Counts/ MeV... 7 J/ψ γ π π Counts/ MeV ψ(s) γπ π... M(π π ) (GeV)... M(π π ) (GeV) Counts/ MeV... + J/ψ γ K K Counts/ MeV ψ(s) γ K K Counts/ MeV... + M(K K ) (GeV)... J/ψ γ K K S S... M(K K s) (GeV) s Counts/ MeV... + M(K K ) (GeV)... ψ(s) γ K K S S... M(K K s) (GeV) s Counts/ MeV... J/ψ γηη Counts/ MeV... ψ(s) γηη... M(ηη) (GeV)... M(ηη) (GeV) FIG. 9. Two particle invariant mass distributions for J/ψ,ψ(S) γps,ps (PS π ±,π,k ±,K S,η) decays, for M(PS,PS) <. GeV. to keep resonance widths also free, and they were fixed at their PDG values []. The instrumental resolution widths, as given in Table II, were convolved with the PDG widths to obtain the peak shapes used in the fits. We use the following relativistic BreitWigner lineshape BW(M) [MΓ(M)] [M M ] +[MΓ(M)] () where M and Γ are the nonrelativistic values for the

9 9 Counts/ MeV... J/ψ γ π + π Counts/ MeV... 7 J/ψ γ π π... M(π + π ) (GeV)... M(π π ) (GeV) Counts/ MeV... + J/ψ γ K K Counts/ MeV... J/ψ γ K K S S... + M(K K ) (GeV)... M(K K s ) (GeV) s FIG.. Illustrating fits to the invariant mass distributions for J/ψ γππ and γkk. mass and width of the resonance, respectively. The massdependent resonance width Γ(M) is given by ( )( ) l+ [ ] M p B(p) Γ(M) = Γ M p B(p ) () where Γ is the nonrelativistic width, p and p are the momenta of the resonance s decay products in its centerofmass frame for masses M and M, respectively. B(p) is the L dependent BlattWeisskopf barrier factor, defined as B(p) =, / +(Rp), and / +(Rp) /+(Rp) /9 for resonances with L =,, and, respectively, and R = GeV. The product branching fractions B (J/ψ,ψ(S) γr) B (R PS,PS) are obtained as N R /ǫ R B (J/ψ,ψ(S) γr) B (R PS,PS) = N(J/ψ,ψ(S)) () where N R is the number of events in the resonance R, for which the efficiency is ǫ R as given in Table II, and N(J/ψ,ψ(S)) are the number of J/ψ and ψ(s) in the data. A. Resonances in J/ψ Radiative Decays to Pseudoscalar Pairs The invariant mass spectra for J/ψ radiative decays and fits to the observed resonances are shown in Fig.. Table III lists the fit results for the resonances. We note that in general, the results shown in Table III show that for both ππ and KK the masses determined for charged and neutral decays are consistent within their respective statistical errors, and the branching fractions are consistent with the isospin expectation of factor two between charged and neutral decays. We now discuss the results for individual resonances separately. The scalar f () in J/ψ decay The ππ spectrum: In the π + π spectrum there is a nontrivial yield in a broad bump in the low energytail region of the strong peak of ρ(77). It is tempting to ascribe it to f (). However, such an assignment is unsupportable because we note that there is no evidence for a corresponding enhancement in the π π spectrum which is cleaner in that mass region because of the absence of the ρ peak. The vector ρ(77) in J/ψ decay As mentioned in Sec. IV the prominent peak for ρ(77) arises in the π + π spectrum because of the decay with largebranching fraction B(J/ψ ρ π ) =.(7), ρ π + π, with one photon missing from the π γγ decay. As expected, the peak is missing in the π π spectrum. The tensor resonance f (7) in J/ψ decay The ππ Spectra: This is the strongest excited peak in both the π + π and π π spectra. Fitted sepa

10 TABLE III. Fit results for J/ψ decays to γr, R PS,PS. For each observed resonance R, results are presented for decays to π + π (χ /d.o.f =.), decays to π π (χ /d.o.f. =.), decays to K + K (χ /d.o.f. =.), and decays to K SK S (χ /d.o.f. =.). For fits to π + π and K + K, the resonance masses were kept free. For all fits their widths were fixed at their PDG values. The summed branching fractions are B B (ππ) = B B (π + π )+B B (π π ), and B B (KK) = [B B (K + K )+B B (K SK S)] (/). J/ψ γr R ππ R, M(R), Γ(R), MeV M(MeV) N B B f (7), π + π 9() 7().(9) 7(),.(.9.) π π ().() ππ (7) 7.() f (), π + π 7() ().() (), 9(7) π π ().(7) ππ 7().(9) f (7), π + π 7(7) () 7.9() 7(), () π π () 9.() ππ () 7.() f (), π + π 9() 9(9).() (), 9(9) π π ().(9) ππ 7().() J/ψ γr R KK R, M(R), Γ(R), MeV M(MeV) N B B f (7), K + K () (9).() (), (77) KSK S ().() KK 7(9).9(7) f (), K + K () (7).() (), 7( ) KSK S ().() KK 7(9) 7.9() f (7), K + K 7() (7).() 7(), () KSK S ().() KK 9() 7.() f (), K + K () ().(9) 9(), () KSK S 9(9) 7.() KK 9().(9) rately we obtain for π + π decay the mass M(π + π ) = 9() MeV, and B B (π + π ) =.(9) and for π π decay and B B (π π ) =.(). The total branching fraction for ππ decay is therefore B B (ππ) = 7.(). These results compare with M = ()MeV,B B (π + π ) = 9.(7)()), B B (π π ) =.(9)(), and B B (ππ) =.()(9),obtainedbyBESII[9]fordirectly produced J/ψ. We note that the BES spectra have much larger backgrounds which may account for the substantial difference from our results. The KK Spectra: The f (7) is not strongly excited in KK decays. In the K + K spectrum there is a prominent bump in the low energy tail of the strong peak due to the excitation of f () resonance. It can not be fitted with a single peak corresponding to f (7), and requires including another larger mass resonance, for which the best candidate is the controversial f (7), which we discuss below. Because of the isospin difference of factor two, and the nearly factor two smaller efficiency, the corresponding enhancement the KS K S spectrum has too small statistics to make independent analysis. We do not fit for f (7) KK, but include its small contribution by appropriately scaling down N ππ (f (7)) by the PDG [] ratio B(f (7) KK)/B(f (7) ππ) =.(), to obtain N K+ K (f (7)) = and N K S KS (f (7)) =. We do not, however, report B B (KK) for f (7). BES II reports that there is evidence for f (7) decays into KK, but also does not report any numerical results []. The scalars f (7) and f () in J/ψ decay The ππ Spectra: The resonance f (7) is considered to be the crucial ingredient in the search for scalar glue

11 balls [,]. However, its mere existence is considered questionable. As a result, PDG () estimates the mass of f (7) to be in the range of MeV and its width to be between MeV. In both the π + π and π π spectra a small, relatively narrow, bump is visible in the high energy tail of the large f (7) peak. Since the wellknown resonance f () is expected in the same tail region, it is difficult to assign the observed bump uniquely to either f (7) or f (), This problem was also noted by Belle [7]. Our attempt to make the best possible choice is based on the following exercise. We analyze our π + π spectrum assuming that a relatively narrow resonance with width between and MeV is responsible for the bump we observe in the high energy tail region of f (7). As noted before, we do not have enough statistics to fit the observed structures keeping both the mass and width of this resonance free. We therefore assume five different values of the width, Γ =,,,7, and MeV, and fit the π + π spectrum for the best value of the mass corresponding to each assumed value of the width. The corresponding best fit values are found to be respectively: M (MeV) = (), (), 9(7), (9), and (). These results indicate that the observed narrow bump can be only assigned to f () (for which many recent measurements report M 7 MeV []), and not to f (7). This does not mean that a broad f (7) does not exist, but only that we have no evidence for its existence in our ππ spectra which are dominated by the f (7) resonance. Our attribution of the bump to the excitation of f (), assuming the PDG [] values of the width Γ = 9 MeV [], leads to M(f ()) = 7() MeV, B B (π + π ) =.() and B B (π π ) =.(7). The summed branching fraction is B B (ππ) =.(9). BES II has analyzed their data for radiative decays of directly produced J/ψ into pion pairs, and reported M() = () MeV, B B (π + π ) =.7()(), B B (π π ) =.()(), and B B (ππ) =.()() []. The KK Spectra: Both the K + K and KS K S spectra are dominated by the f () and f (7) resonances. However, there is appreciable yield spread over morethanmevinthe K + K spectrum in the region of the low energy tail of f (). To determine whether this yield should be assigned to the resonance f (7) or f (7), we have done the same exercise as described abovefor f (). We fit the K + K spectrum to determine values of masses in this region with assumed values of width, Γ =,,,, and MeV. The best fit masses are found to vary from to MeV. Allowing the width to be a free parameter results in M = () MeV, Γ = (77) MeV. We therefore assign this enhancement to f (7). The K + K spectrum leads to B B (K + K ) =.(), and the KS K S spectrum leadstob B (KS K S ) =.(). The summed result is B B (KK) =.9(7), taking account of a factor / for the undetected K L K L decay, as noted for all kaon decays in the Table III caption. BESIIalsoobserveskaonyieldsinthisgeneralmass region, but does not quote any numerical results [9]. The tensor f () in J/ψ decay There is no evidence for the excitation of f () in the ππ decay spectra. In the KK spectra f () is prominently excited, with the K + K spectrum leading to M(f ()) = () MeV. The K+ K spectrum leads to B B (K + K ) =.(), and the KS K S spectrum leads to B B (KS K S ) =.(). The summed result is B B (KK) = 7.9(). The result from BES II with large background in the invariant mass spectra is B B (KK) =.()( 9 )( ) []. The scalar f (7) in J/ψ decay The scalar f (7) is prominently excited in both ππ decays and KK decays (Fig. ). The π + π spectrum leads to the mass M(f (7)) = 7(7) MeV and B B (π + π ) = 7.9(), and the π π spectrum leads to B B (π π ) = 9.(), with the sum B B (ππ) = 7.(). BES II has reported B B (π + π ) =.()(7) and B B (π π ) =.()(), and therefore B B (ππ) = 9.7()() [9]. The KK Spectra: The K + K spectrum leads to M(7) = 7() MeV, and B B (K + K ) =.(), the KS K S spectrum leads to B B (KS K S ) =.(). Our summed result is B B (KK) = 7.(). The corresponding result of BES II is B B (KK) = 9.(9)( )( ) []. Of the seven resonances listed in Table III, unambiguous results for both ππ and KK decays are observed for only one resonance, f (7). This is largely due to the fact that all the resonances have widths much larger than the experimental resolution widths, and they strongly overlap. As a result, as described in Sec. IX, only for f (7) we are able to test the flavor independent character attributed to pure glueballs. Both K + K and KS K S spectra are very well defined in their high energy tail region, and show no evidence for enhancement in the M = 79 region which has been occasionally reported. The M = MeV region in J/ψ decay The mass region above f (7) is complicated with several f, and f resonances reported []. In particular, radiative decays of J/ψ have been reported to f (), f (), f (9), f (9), and to the much talked about tensor f (). In the ππ and KK mass spectra in the present investigation a broad enhancement at M MeV is seen in the ππ decays, and another broad enhancement at M MeV is seen in the KK decays. The background is generally much larger than in the lower mass regions. This makes resonance assignments in this mass region rather speculative. With this

12 Counts/ MeV... + ψ(s) γ π π Counts/ MeV 7... ψ(s) γ π π... M(π + π ) (GeV)... M(π π ) (GeV) Counts/ MeV ψ(s) γ K K Counts / MeV... ψ(s) γ K K S S... + M(K K ) (GeV)... M(K K S ) (GeV) S FIG.. Illustrating fits to the invariant mass distributions for ψ(s) γππ and γkk. Counts/ MeV... + ψ(s) γ π π... M(π + π ) (GeV) caveat we fit the ππ spectra in this region with f () with the PDG width of 9 MeV, and the KK spectra with f () with the PDG width of MeV. Fits to the ππ spectra lead to M(f ()) = 9() MeV, B B (π + π ) =.(), B B (π π ) =.(9), and the sum B B (ππ) =.(). Fits tothekk spectraleadtom(f ()) = ()MeV, B B (K + K ) =.(9), and B B (K S K S ) = 7.(), and the total B B (KK) =.(9). Counts/ MeV ψ(s) γ K K... + M(K K ) (GeV) FIG.. Data for π + π and K + K mass distributions with superposed continuum contributions (in blue) estimated by normalizing continuum yields observed in offresonance data taken at s 7 MeV. B. Resonances in ψ(s) Radiative Decays to Pseudoscalar Pairs While searches for glueballs in radiative decays of J/ψ have been reported before, similar searches in ψ(s) radiative decays have not been reported before. The only reports for ψ(s) radiative decays which exist are a BES II publication with a sample of. million ψ(s) [9], and an unpublished report with million ψ(s) []. We present in this paper our results for the radiative decays observed in a sample of. million ψ(s). As pointed out in Sections IV and V, and shown in Fig. 7, the results for ψ(s) decays are very similar to those of J/ψ decays. The event statistics is smaller, due primarily to the % rule, and the backgrounds are relatively larger. The fits to the invariant mass spectra, and the determination of the product branching fractions

13 TABLE IV. Fit results for ψ(s) decays to γr, R PS,PS. For each observed resonance R, results are presented for decays to π + π (χ /d.o.f. =.), decays to π π (χ /d.o.f. =.9), decays to K + K (χ /d.o.f. =.9), and decays to K SK S (χ /d.o.f. =.9. For fits to π + π and K + K, the resonance masses were kept free. For all fits their widths were fixed at their PDG values. The summed branching fractions are B B (ππ) = B B (π + π ) + B B (π π ) and B B (KK) = [B B (K + K )+B B (K SK S)] (/). ψ(s) γr R ππ R, M(R), Γ(R), MeV M(MeV) N B B f (7), π + π 7() 9() 7.(7) 7(),.(.9.) π π ().9() ππ 9(7).9(9) f (), π + π (9) ().() (), 9(7) π π (7).() ππ 7().() f (7), π + π 7() 7().() 7(), () π π ().() ππ 9().() f (), π + π 99(7) ().() (), 9(9) π π 9(7).() ππ 7(9).() ψ(s) γr R KK R, M(R), Γ(R), MeV M(MeV) N B B f (7), K + K () ().9() (), (77) KSK S 7().() KK 7().() f (), K + K (9) ().() (), 7( ) KSK S ().9() KK ().9() f (7), K + K 9() 9().() 7(), () KSK S ().9() KK 7().7() f (), K + K (7) ().() 9(), () KSK S (9).() KK 7().() for the ψ(s) decays was done in the same manner as for J/ψ decays. The results are tabulated in Table IV, and the fits to the chargedand neutraldecaysfor ππ andkk spectra are shown in Fig.. The vectors ρ(77) and φ() in ψ(s) decay In the π + π and K + K invariant mass spectra we note the presence ofthe ρ(77)and φ()peaks asalsoindicated by the diagonal enhancements in the corresponding Dalitz plots of Fig.. These arise due to the production of ρ(77) and φ() vector mesons via initial state radiation. To confirm this, we have analyzed our offψ(s) resonance data at s = 7 MeV with the integrated luminosity of. pb. To compare the π + π and K + K yields observed in these data with those observed in ψ(s) datawe multiply the yields by the relativeluminosity factor L(ψ(S))/L(7) = /. =.7, and the factor (s /s) = (/7) =.. The normalized yields are shown in Fig. superposed on the ψ(s) data. It is seen that the normalized contributions fully account for the yield of ρ(77) in the π + π spectrum, and for the yield of φ() in the K + K spectrum. It is also seen that the background contribution due to this ISR yield is predicted to be negligible for masses > GeV. The tensor f (7) in ψ(s) decay As in J/ψ decay, the f (7) is prominently excited in both π + π and π π decays. On the high energy shoulders in both spectra we again see enhancement which is well fitted with the narrow f () rather than narrow or wide f (7). The fit to f (7) results in M = 7() MeV, B B (π + π ) = 7.(7), B B (π π ) =.9(), and the sum

14 B B (ππ) =.9(9). In the K + K spectrum there is a broad enhancement below MeV which cannot be assigned to f (7). Instead, it is well fitted with a wide resonance with mass around MeV, as discussed below. The scalar f (7) in ψ(s) decay In the J/ψ decay to ππ it was not possible to discern the possible excitation of f (7) because of the strong excitation of f (7). As described in Sec. V.A, the small shoulder seen on the high energy tail of f (7) could be assigned to f () but not to f (7). In the KK spectraforψ(s)the f (7)excitationismuchweaker. This allows us the opportunity to fit the mass region below MeV. The ψ(s) spectra have much smaller statistics than the J/ψ spectra, which makes it difficult to determine the properties of a wide, weakly excited resonance like the f (7). To explore its possible excitation, we have fitted the K + K spectrum with the width of f (7) fixed to Γ =,, MeV. The resulting masses were M = (), (), 9() MeV, i.e., in essential agreement with the determination of M = () MeV and Γ = MeV for J/ψ decay. Forourfinalfit, wethereforeassumeastatewithwidth ofmev. ThefittotheK + K spectrumyieldsamass of M = () MeV. We identify it with f (7). The fit results to B B (K + K ) =.9(), B B (K S K S ) =.(), and the total B B (KK) =.(). This fit has been made excluding the two bins near M MeV, which would correspond to a very narrow state of Γ < MeV, for which no plausible candidate exists. The scalar f () in ψ(s) decay The absence of f () in ππ decay allows us to fit for f () with the results: M(f ()) = (9) MeV, B B (π + π ) =.(), B B (π π ) =.(), and the sum B B (ππ) =.(). The strong excitation of f ()in KK decaysmakes it impossible to say anything about possible excitation of f () in KK decays. The tensor f () in ψ(s) decay Noevidencefortheππ decayoff ()wasfoundinthe radiative decay of J/ψ, and none is found in the radiative decay ψ(s) γππ, either. In K + K decay strong excitation of f () is observed, and the fit to the K + K spectrum leads to M(f ()) = (9) MeV andb B (K + K ) =.(). Thefit tothe KS K S spectrum leads to B B (KS K S ) =.9(). The total branching fraction is B B (KK) =.9(). The scalar f (7) in ψ(s) decay The scalar f (7) is strongly excited in both ππ and KK decays. Fits to the ππ spectra lead to M(f (7)) = 7() MeV, B B (π + π ) =.(), and B B (π π ) =.(), and the sum B B (ππ) =.(). Fits to the KK spectra lead tom(f (7)) = 9()MeV,B B (K + K ) =.(), B B (K S K S ) =.9(), and the total B B (KK) =.7(). The M = MeV region in ψ(s) decay As mentioned for J/ψ decays, in this region backgrounds are large and only a large, wide peak is observed in both ππ and KK spectra. Each of them may be composite of the excitation of more than one resonance. If analyzed as single resonances, fits to the ππ spectra lead to M(f ()) = 99(7) MeV, B B (π + π ) =.(), B B (π π ) =.(), and the sum B B (ππ) =.(). Anisovich and Uman have reported the observation of a resonance with this mass in p p annihilations []. Fits to the KK spectra lead to M(f ()) = (7) MeV, B B (K + K ) =.(), B B (K S K S ) =.(), and the total B B (KK) =.(). Augustin [7] has reported excitation of a resonance with this mass in radiative decay of J/ψ. C. The ξ() Tensor Glueball Candidate While the main interest in glueball searches has been in the scalars, special interest has existed in the claim for a narrow tensor state with mass MeV, and width Γ MeV. This claim has a checkered history. Originally claimed by Mark III [] in the radiative decay of J/ψ, and the decay of ξ() into K + K and KS K S, it has gone through several claimed observations [,, ], and a number of nonobservations [7,, 9 ]. The latest of these consist of its observation by ITEP in π p formation and decay into KS K S [], and its nonobservation by BaBar in the radiative decay of J/ψ into K + K and KS K S [], and its nonobservation by BES III in the radiative decay of J/ψ and decay into ηη []. The present measurement allow us to search for it in the radiative decays of both J/ψ and ψ(s), with its decays into ππ, KK, and ηη pseudoscalar pairs. In the analyses of the invariant mass spectra as described above, we add a resonance with M = MeV, and widths of and MeV, and determine upper limits at 9% confidence level for the product branching fractions for the existence of ξ(). The results are listed in Table V. They can be summarized as upper limits of in J/ψ decays and < in ψ(s) decays. Our limits in J/ψ decays are consistent with. reported by DM [7] for K + K and KS K S decays, and reported by BaBar [] in K + K and KS K S decays. No evidence for its existence was found by BES III in radiative decay of a very large sample of. J/ψ to ηη [].

15 TABLE V. Upper limits at 9% confidence level for the production of the presumed tensor glueball candidate ξ() in the reactions J/ψ,ψ(S) γr, R PS,PS. 9% CL UL Γ(ξ), MeV N B B B B (ref.) J/ψ γ(π + π ) /./../..( 7 ) [] γ(π π ) /./../.9 γ(k + K ) /./..7/..( ) [],.( 9 ) []. (9% CL) [7] γ(ksk S) /./9../..( 7 ) [],.7( ) []. (9% CL) [7] ψ(s) γ(π + π ) /.9/../. γ(π π ) /.7/../. γ(k + K ) /.9/../. γ(ksk S) /.7/..7/. TABLE VI. Summary of results for J/ψ and ψ(s) radiative decays to light quark resonances. For entries with two uncertainties, the first listed uncertainties are statistical and the second are systematic as described in Sec. VIII and listed in Tables VIII and IX. From J/ψ From ψ(s) B (ψ(s)) B (J/ψ) M (MeV) B B M (MeV) B B (%) f (7) ππ 9()() 7.()() 7()().9(9)(9).7(7)() f (7) KK ()().9(7)() ()().()() 7.(7)() f () ππ 7()().(9)() (9)().()().()() f () KK ()() 7.9()(7) 7(9)().9()().(9)() f (7) ππ 7(7)() 7.()() 7()().()() 9.7()() KK 7()() 7.()(9) 9()().7()().7()(7) f () ππ 9()().()(7) 99(7)().()(9) 7.7()() f () KK ()().(9)() (7)().()().()() TABLE VII. Product branching fractions for χ c, PS,PS decays (based on Ref. []). B B χ c χ c χ cj π + π.().() χ cj π π 7.().(7) χ cj ππ.9().() χ cj K + K.().() χ cj KSK S.().9() χ cj KK.9(9).7() χ cj ηη 9.().(7) χ cj η η 9.().() χ cj ηη.().(9) VI. SUMMARY AND COMPARISON OF RESULTS FOR RADIATIVE DECAYS OF J/ψ AND ψ(s) TO RESONANCES WITH MASSES M < MEV In Table VI we compare the results for all resonances observed in the radiative decays of J/ψ and ψ(s) to pairs of charged and neutral pseudoscalars. Several features are notable in Table VI. As mentioned in Sec. V.B, the results for J/ψ and ψ(s) decays are very similar for both charged and neutral decays. The masses are generally in agreement well within their statistical errors. Branching fraction results from BES, when available, are generally not in good agreement with ours. The BES results are based on directly produced J/ψ, and as mentioned before, the BES invariant mass spectra contain large backgrounds due to ISR contributions. The Socalled % Rule: The % rule mentioned earlier is generally invoked for hadronic decays of vec

16 Counts / MeV ψ(s) γ π π tor resonances of charmonium J/ψ, ψ(s),...ψ(ns) because both their leptonic decays via one photon and their hadronic decays via three gluons are proportional to the wave function at origin. Actually, the radiative decays to hadrons, via a photon plus two gluons have the same proportionality to the wave function at origin, so that we expect, for example, the ratio of radiative decays to pseudoscalar pairs Counts / MeV..... M(π + π ) (GeV)..... ψ(s) γπ π B(ψ(S) γps,ps)/b(j/ψ γps,ps) = B(ψ(S) e + e )/B(J/ψ e + e ) =.79 /.97 =.()% In Table VI we list this ratio for the various resonances we observe. The ratio is found to vary between.()% and.()%. This is similar to the variation which has been observed for the pure hadronic decays[] M(π π ) (GeV)..... VII. RESONANCES WITH MASSES > MEV Counts / MeV + ψ(s) γ K K A. Search for Glueballs with Masses in the.. GeV Region Counts / MeV Counts / MeV M(K K ) (GeV)..... ψ(s) γ K K S S..... M(K K S ) (GeV) S..... ψ(s) γηη One of the main reasons for our studying ψ(s) radiative decays in addition to those of J/ψ was to extend our glueball searches to masses in the region.. GeV. As stated earlier, and illustrated in Fig., lattice calculationspredictaj PC = + glueballwithmassm( + ) = ()() MeV, and a J PC = + glueball with mass M( + ) = ()() MeV. These + and + positive charge conjugation states can be excited in radiative decays of ψ(s), and we have searched for them and any other enhancements in the entire region, M =.. GeV, in their decays into pseudoscalar pairs, assuming three different values of their widths, Γ(X) = MeV, MeV and MeV. The results are displayed in Fig.. No notable enhancement is found anywhere in the mass region, and upper limits for their branching fractions at 9% confidence level have been determined. The upper limits for all product branching fractions B B are determined to be all <, except for decay to ηη for which they become as large as because of the very small number of observed counts and small B(η γγ) M(ηη) (GeV) FIG.. Two particle invariant mass distributions for selected ψ(s) γps,ps, (PS π ±,π,k ±,K S,η) decays, for the region M >. GeV. The χ c and χ c states are clearly seen. No evidence for any other states is seen. B. The χ c and χ c resonances of charmonium As shown in the Dalitz plot of Fig., and invariant mass projections in Fig. for ψ(s) radiative decay, we observe strongly populated χ c and χ c states of charmonium in their decay info ππ, KK, and ηη. The branching fractions for these decays were reported by us earlier in the CLEO publication of Ref. []. The B (χ c,c PS,PS) branching fractions reported there, multiplied by B (ψ(s) γχ c ) = 9.% and

17 7 (9% CL) B B + ψ(s) γ π π Counts / MeV (9% CL) B B ψ(s) γ π π Counts / MeV M (GeV) M (GeV) (9% CL) B B + ψ(s) γ K K Counts / MeV (9% CL) B B ψ(s) γ K K S S Counts / MeV M (GeV) M (GeV) (9% CL) B B ψ(s) γ ηη Counts / MeV M (GeV) FIG.. Upper limits at 9% confidence level of B (ψ(s) R) B (R PS,PS) as function of resonance mass, M. The differently colored curves, from top to bottom, correspond to upper limits for assumed widths Γ(R) =, and MeV. The twopseudoscalar invariant mass spectra are shown as the shaded histograms, with the scale of the observed counts on the right ordinate. The thick horizontal lines illustrate the mass ranges predicated in Ref. [] for the glueball candidates with J PC = + and +. B (ψ(s) γχ c ) = 9.% used there, are quoted as B B in Table VII. Ochs [] has emphasized that while it is difficult to study the flavorindependent decay of twogluon glueball candidates to pseudoscalar pairs, the corresponding decays of χ c (J = ) and χ c (J = ), which are mediated by two gluons, provide an excellent opportunity to test the flavor independence ansatz. The values for ππ, KK, andηη decayslisted intablevii leadto theratiosbelow:

18 TABLE VIII. Summary of systematic uncertainties on branching fractions (in units of percent). J/ψ γπ + π /γπ π ψ(s) γπ + π /γπ π f (7) f () f (7) f () f (7) f () f (7) f () Reconstruction & PID / / / / / / / / Kinematic Fit / / / / / /9 / / BW Widths / / / 7/ / / / 7/ Background / / / 7/ / / /9 / Binning / / / / / /9 / / ρπ Rejection / / / / Number of ψ(s) / / / / / / / / Total / / / / / / / / J/ψ γk + K /γksk S ψ(s) γk + K /γksk S f (7) f () f (7) f () f (7) f () f (7) f () Reconstruction & PID / / / / / / / / Kinematic Fit / / / / /7 / / / BW Widths / / / / / 7/ / / Background / / / / / / 7/ / Binning / / / / / / / / K K Rejection / / / / Number of ψ(s) / / / / / / / / Total / /9 9/7 / / / / / Pseudoscalar pairs ππ KK ηη η η ηη Twogluon decays (theoretical) =... χ c decays =.().().().().() χ c decays =.().9().().().() χ c decays (divided by phase space factor p) =.().7(9).().7().() χ c decays (divided by phase space factor p ) =.().().().().() The branching fraction ratios for χ c decays are seen to agree within experimental uncertainties with the theoretical expectations for gluon decays. The χ c ratios also agree with these expectations, except for decays to η η. VIII. SYSTEMATIC UNCERTAINTIES Our estimates of systematic uncertainties in product branching fractions J/ψ γππ, γkk, and ψ(s) γππ, γkk, are listed in Table VIII. We estimate uncertainties due to reconstruction and particle identification as follows. We find a reconstruction uncertainty of % per charged particle, with an additional % and % per identified π ± and K ±, respectively, for the use of de/dx and RICH information. We find a K S π+ π identification uncertainty of %. We also find uncertainties of % per π and % for the radiative photon. We add the contributions for each mode in quadrature to obtain the systematic uncertainties for this source. We estimate the uncertainty in our requirement that the kinematic fit have χ < by analyzing data with the requirements χ < and χ <, and we take the largest variations as the systematic uncertainty due to this source. In resonance fits, our fixed parameters are the Breit Wigner widths of the resonances and the MCdetermined mass resolution. To estimate the uncertainties due to the fixed resonance widths, we refit the data, varying the widths by ±σ of their PDG values [], and we take the largest variation for each resonance as the systematic uncertainty due to this source. We also refit the data with the MCdetermined mass resolution varied by ±%, and find that this has no measurable effect our results, since the resonances in our fits are all much wider than the mass resolutions. The uncertainty due to the backgrounds in the fits is obtained by varying background normalization by ±σ around the central value obtained fromthe fits. Theuncertaintydueto the MeVbinning of the invariant mass spectra is obtained by analyzing the data with MeV and MeV bins and taking the largest