At the LHC energies, the underlying observables are of major topic of interest in high multiplicity $p+p$ collisions. Multiple Parton Interactions (MPI) is one of them, in which several interactions occur in a single $p+p$ collision. It is believed that MPI is the main reason behind the high multiplicity in $p+p$ collisions at the LHC. It was believed that MPI has only effect to the soft particle production, but recent ALICE result reveals that it can also affect the hard-particle production. In such case, the self normalized yield of heavy particle like $rm J/psi$ shows an increasing trend with event multiplicity. In the present contribution, we discuss the energy and multiplicity dependence of charmonium production to understand the effects of MPI on charmonium production.