The NRQCD factorization approach for quarkonium production at hadron collid
ers is reviewed. The prediction of J/psi and psi (2S) transverse polarizati
on at large transverse momentum is confronted with recent experimental data
, and potential shortcomings of the theoretical analysis are discussed.