Phenomenologically relevant inclusive quarkonium production processes: For two-photon scattering [770,1189], e+ e- scattering , photoproduction [770,1191] and hadroproduction [1184,1192?195]. The polarized cross sections have been calculated for photoproduction  and hadroproduction [1181,1183,1197,1198]. In , a global fit of the J/ CO LDMEs to 26 sets of inclusive J/ production yield data from 10 different pp, p, , and e+ e- experiments was done; see the second column of Table 13 for the fit results. This fit describes all data, except perhaps the two-photon scattering at LEP , reasonably well. This fit is overconstrained, and get CPI-203 practically independent of possible low- pT cuts (unless such high pT cuts are chosen that all data except hadroproduction drop out of the fit ). Furthermore, the resulting LDMEs are in accordance with the velocity scaling rules predicted by NRQCD; see Table 5. Thus the fit is in itself already a nontrivial test of the NRQCD factorization conjecture, especially since the high-z photoproduction region can now also be well described, which had been plagued by divergent behavior in the earlier Born analyses [1200,1201]. However, in  it was shown that these CO LDME values lead to predictions of a strong transverse J/ polarization in the hadroproduction helicity frame, which is in contrast to the precise CDF title= fpsyg.2017.00209 Tevatron run II measurement ; see Fig. 33d. On the other hand, in  it was shown that both the measured J/ hadroproduction cross sections and the CDF run II polarization measurement  can, even at the highest measured pT values, be well described when choosing one of the three CO LDME sets listed in columns four through six of Table 13. These LDMEs, however, result in predictions for e+ e- annihilation and photoproduction which are factors four to six above the data; see Fig. 13e . Third, the calculation  is the first NLO polarization analysis to include feed-down contributions. To this end, the CO LDMEs of J/, (2S) and c J were fitted to CDF [1141,1146] and LHCb [1147,1148,1182] unpolarized production data with pT > 7 GeV; see column three of Table 13. These fit results were then used for the predictions of Fig. 33e , taking the (2S) and c J feed-down con-2981 Page 78 of 241 Table 13 Overview of different NLO fits of the CO LDMEs. Analysis  is a global fit to inclusive J/ yield data from 10 different pp, p, ee, and experiments. In , fits to pp yields from CDF [1141,1146] and LHCb [1147,1148,1182] were made. In , three values for their combined fit to CDF J/ yield and polarization [1158,1159] data are given: A default set, and two alternative sets. Analysis  is a fit to the c2 /c1 production ratio measured by CDF Butenschoen, Kniehl :O title= 17470919.2015.1029593 J/ (3 S1 ) /GeV O J/ (1 S0 ) /GeV3  O J/ (3 S1 ) /GeV3  O J/ (3 P0 ) /GeV5  O(2S) (3 S1 ) /GeV3 (2S) (1 S  ) /GeV3 O 0  O(2S) (3 S1 ) /GeV3  O(2S) (3 P0 ) /GeV5  O0 (3 P0 ) /GeV5  O0 (3 S1 ) /GeV3 Eur. Phys. J. C (2014) 74:2981 .