Relativistic effect of $J/\psi$ hadroproduction in large $p_T$ region

By combining NRQCD factorization and collinear factorization, we compute a series of relativistic corrections for $J/\psi$ hadroproduction to all orders in $v^2$ at large $p_T$ limit. The $v^2$ expansion converges well for all channels. We find that the ratio of relativistic correction term to the corresponding leading term is independent of kinematic variables for any channel, which generalizes the proportional relations found in previous works to all orders.