Update of hadronic decays of $J/\psi$ and $\psi(2S)$ though virtual photons

The hadronic decay branching ratios of $J/\psi$ and $\psi(2S)$ through virtual photons, $B(J/\psi, \psi(2S) \rightarrow \gamma^*\rightarrow \text{hadrons})$, are updated by using the latest published measurements of the $R$ value and the branching ratios of $J/\psi, \psi(2S) \rightarrow l^+l^-$. Their respective precision increases by about 4 and 3 times.

The J/ψ(1 3 S 1 ) and ψ(2 3 S 1 ), J PC (1 −− ) states of charmonium have extremely small width. The reason is usually referred to by the phenomenological Okubo-Zweig-Iizuka rule, which states that processes in which the initial quark pairs cannot appear as part of the final state particles are suppressed. Therefore, the J/ψ only decays through virtual photons (vacuum polarization), three gluons, two gluons plus a single photon, or into η c via M1 transition and three photons. The determination of the partial widths could play an important role in understanding the nature of J/ψ decays. The hadronic final state through virtual photons and gluons is indistinguishable. Fortunately, the branching ratios B(J/ψ → γ * → hadrons) can be calculated using the relation from Ref. [1] where assumes that the ratios of the cross sections for virtual photons decaying to hadrons and leptons, denoted as σ h and σ l respectively, are the same for both electron-positron and cc annihilations. This relationship also holds for the decay of ψ(2S).
The best published determinations on B(J/ψ, ψ(2S) → γ * → hadrons) are from 2004 and were determined by Seth [2]: The R value used in this result was 2.28 ± 0.04.
In the studies of B(J/ψ, ψ(2S) → l + l − ), the most recent and the most precise results were published by CLEO c and BESIII, which reached relative precision of 0.6% and 2.1%, respectively. Reasons for these improvements are better detectors, larger data samples, and carefully studied systematic uncertainties. By combining all available results [10], the precisions on B(J/ψ → l + l − ) and B(ψ(2S) → l + l − ) improved by factors of 4 and 2.5, respectively, as shown in Table I.
The R values between 2.0 and 3.73 GeV are used to approximate the ones at J/ψ and ψ(2S) peaks, as shown in Fig. 1. The results of the R scans used in this study are from the MARK-I, the DM 2, the BESII, the KEDR, the BESIII, etc. It can be seen that the precision of R values has improved significantly in the last decades, benefiting from better detectors, increasing data samples, and precise knowledge of initial state radiation correction and systematic uncertainties.
Various combinations of the existing measurements of the R value are tested to determine the R. The results are summarized in Table. II. It can be seen clearly that the BESII, the KEDR, and the BESIII can achieve comparable precision, and all three are better than the one used in Ref. [2]. It should be noted that the BESII published new measurements of the R value in 2006 [6,7] and 2009 [8] after the publication of Ref [2]. The precision reaches 0.4% if all the results are combined, which is going to be used in the following calculations.
The branching ratios B(J/ψ, ψ(2S) → γ * → hadrons) are calculated using equation (1) and the combined fit result on R from Table II  I. Summary of the branching ratios B(J/ψ, ψ(2S) → l + l − ) and their averages in PDG 2002 [11] and 2022 [10]. In the 3rd and 5th columns, the branching ratios B(J/ψ → l + l − ) and ψ(2S) → l + l − ) were assuming lepton universality and neglecting the mass difference between electrons and muons. The relative uncertainty of the branching ratios is also provided.  [9]. The R values of the BESII experiment [6-8, 12, 13] are indicated as green dots. Blue rectangles represent the results of KEDR detector [3][4][5]. MARK-1 [14] and DM 2 [15] results are shown as cyan stars and orange crosses, respectively. The values from PLUTO [16] are denoted by brown diamonds. Magenta triangles stand for the results of the Crystal Ball Collaboration [17]. B(J/ψ, ψ(2S) → l + l − ), the relative uncertainties are reduced significantly. The B(J/ψ → γ * → hadrons) and B(ψ(2S) → γ * → hadrons) are improved by about 4 times and 3 times, respectively. But it should be noted that the correlations among R values within a single experiment are neglected, as well as the correlation between the branching ratios for B(J/ψ → e + e − ) and B(J/ψ → µ + µ − ).