Non-relativistic quantum chromodynamics (NRQCD) [1]} is generally considered as the best
theory approach to study quarkonium production.
It is important to compare its predictions to the experimental measurements
recently made available by the LHC experiments, which already reach rather high
quarkonium transverse momentum, \(p_{\rm T} \approx 100\)  GeV,
where the calculations are expected to be more reliable.
If discrepancies between the theory predictions and the measured patterns are found,
it is crucial to investigate if those differences stem from problems in the conceptual foundations
of the theory or from approximations and inaccuracies of the fixed-order perturbative calculations
available at present.
