Color-octet quarkonia production. II

We calculate the lowest order hadronic cross sections for producing colored heavy quark-antiquark pairs in $L=S=0\mathrm{}$ and $L=S=1\mathrm{}$ configurations. Such $Q\overline{Q}\left[{}_{}{}^1{}_{}{}^{}S_0^{\mathrm{}\left(8\right)\mathrm{}}{}_{}{}^{}\right]$ and $Q\overline{Q}\left[{}_{}{}^3{}_{}{}^{}P_J^{\mathrm{}\left(8\right)\mathrm{}}{}_{}{}^{}\right]$ states hadronize into ${\psi }_Q$ quarkonia at the same order in the NRQCD velocity expansion as previously considered $Q\overline{Q}\left[{}_{}{}^3{}_{}{}^{}S_1^{\mathrm{}\left(8\right)\mathrm{}}{}_{}{}^{}\right]$ pairs. Their contributions to prompt $\psi$ and $\Upsilon$ production at the Fermilab Tevatron bring the shapes of theoretical transverse momentum distributions into line with recent CDF measurements. We find that the best fit values for the linear combinations of $Q\overline{Q}\left[{}_{}{}^1{}_{}{}^{}S_0^{\mathrm{}\left(8\right)\mathrm{}}{}_{}{}^{}\right]$ and $Q\overline{Q}\left[{}_{}{}^3{}_{}{}^{}P_J^{\mathrm{}\left(8\right)\mathrm{}}{}_{}{}^{}\right]$ long distance matrix elements which can be extracted from the data are generally consistent with NRQCD scaling rules.