$J/\psi$ production in hadron scattering: three-pomeron contribution.

In this paper we discuss the inclusive $J/\psi$ production in proton-proton collisions from fusion of three pomerons. We demonstrate that this mechanism gets dominant contribution from the region which can be theoretically described by CGC/Saturation approach. Numerically, it gives a substantial contribution to the $J/\psi$ production, and is able to describe the experimentally observable shapes of the rapidity, momenta and multiplicity distributions. The latter fact provides a natural explanation of the experimentally observed enhancement of multiplicity distribution in $J/\psi$ production.