We calculate the lowest order hadronic cross sections for producing co
lored heavy quark-antiquark pairs in L=S=0 and L=S=1 configurations. S
uch Q (Q) over bar[S-1(0)(8)] and Q (Q) over bar[P-3(J)(8)] states had
ronize into psi(Q) quarkonia at the same order in the NRQCD velocity e
xpansion as previously considered Q (Q) over bar[S-3(1)(8)] pairs. The
ir contributions to prompt psi and Upsilon production at the Fermilab
Tevatron bring the shapes of theoretical transverse momentum distribut
ions into line with recent CDF measurements. We find that the best fit
values for the linear combinations of Q (Q) over bar[S-1(0)(8)] and Q
(Q) over bar[P-3(J)(8)] long distance matrix elements which can be ex
tracted from the data are generally consistent with NRQCD scaling rule
s.