Processes $e^+e^- \to c \bar c c \bar c$ and $e^+e^- \to J/\psi +gg$ at $\sqrt{s}=10.59$ GeV

BELLE Collaboration data for $J/\psi$ inclusive production in the processes $e^+e^-\to J/\psi+gg$ and $e^+e^-\to J/\psi +c\bar c$ have been discussed. These data have been compared with the predictions of two pQCD methods: the former method use the information about the $J/\psi$ wave function, the latter one do not use such information and is based on the quark-hadron duality hypothesis. Both these calculation methods lead to essential discrepancy between the theoretic predictions and the experimental data. The production cross section dependence on effective gluon mass has been investigated for the process $e^+e^-\to J/\psi+gg$. The production cross section of doubly charmed baryons $\Xi ^{*}_{cc}$ has been estimated.