We study the production of $J/\psi$ vector mesons as a function charged hadron multiplicity in p+p and p+Pb collisions at LHC energies. employ color glass condensate framework, using running coupling Balitsky-Kovchegov evolved dipole amplitudes, to compute gluon $c\bar{c}$-pair production. use fragmentation functions obtain hadrons, explore two different hadronization schemes for $J/\psi$: non-...