Measurement of Double Charmonium Production in $e^+e^-$ Annihilations at   $\sqrt{s}=10.6$ GeV
  We study $e^+e^-\to J/\psi c\bar{c}$ by measuring the invariant mass distribution recoiling against fully reconstructed $J/\psi$ decays, using 124 fb$^{-1}$ of data collected with a center-of-mass energy of 10.6 GeV with the $BaBar$ detector. We observe signals for $\eta_c(1S)$, $\chi_{c0}$, and $\eta_c(2S)$ in the recoil mass distribution, thus confirming previous measurements. We measure $\sigma(e^+e^-\to J/\psi+c\bar{c}) {\cal B} (c\bar{c} \to >2 \mathrm{charged})$ to be   $17.6\pm2.8\textrm{(stat)}^{+1.5}_{-2.1}\textrm{(syst)}$ fb,   $10.3\pm2.5\textrm{(stat)}^{+1.4}_{-1.8}\textrm{(syst)}$ fb, and $16.4\pm3.7\textrm{(stat)}^{+2.4}_{-3.0}\textrm{(syst)}$ fb with $c\bar{c}=\eta_c(1S)$, $\chi_{c0}$, and $\eta_c(2S)$, respectively. 
The BABAR Collaboration, B. Aubert, et al