Not yet reviewed by Pith; the record is open.

In the present work we study production of $J/\psi + \chi_c$ and $J/\psi + J/\psi$ in LHC. The first process is forbidden at the leading order in gluon fusion due to the $C$-parity conservation, and the first non-vanishing contribution is given by the process with additional emission of real gluon. Considering the direct production of $J/\psi + J/\psi$, where the leading order is allowed, we have found that the contribution from the higher order process with real gluon emission is comparable and even more significant than the leading order. Moreover, account of this higher order effect dramatically changes kinematical distributions. Through the present paper we study in details different channels of paired $J/\psi + J/\psi$ production: direct $J/\psi + J/\psi$ production, feed-down from $J/\psi + \chi_c$ channel, double parton scattering. We also try to find kinematical distributions that are most suitable to separate these different channels.