In this talk, I would discuss the effect of linearly polarized gluons on transverse momentum and rapidity
spectrum of J/psi and Upsilon in unpolarized proton-proton collision. I will start with a brief introduction
to gluon TMDs and quarkonium models. I will show that quarkonium production spectra can be modulated
due to the presence of linearly polarized gluons inside an unpolarized proton.