We present estimates of transverse single spin asymmetry (TSSA) in p + p ↑ → J/ψ + X within the colour evaporation model of charmonium production in a generalized parton model (GPM) framework, using the recently obtained best fit parameters for the gluon Sivers function (GSF) extracted from PHENIX data on TSSA in p+p ↑ → π 0 +X at midrapidity. We calculate asymmetry at √ s = 200 GeV, and compare the results with PHENIX data on TSSA in the process p + p ↑ → J/ψ+X. We also present estimates for asymmetry at √ s = 115 GeV corresponding to the proposed fixed target experiment AFTER@LHC and at √ s = 500 GeV corresponding to the higher RHIC energy. Finally, we investigate the effect of the transverse momentum dependent (TMD) evolution of the densities involved, on the asymmetry.