We employ a small Color Glass Condensate (CGC)+ Non-Relativistic QCD (NRQCD) formalism to compute production at low in proton-proton collisions at collider energies. Very good agreement is obtained for total cross-sections, rapidity distributions and low momentum distributions. Similar agreement is obtained for production. We observe an overlap region in where our results match smoothly to those obtained in a next-to-leading order (NLO) collinearly factorized NRQCD formalism. The relative contribution of color singlet and color octet contributions can be quantified in the CGC+NRQCD framework, with the former contributing approximately of the total cross-section.