
1 Introduction
Charmonium and bottomonium production provides an ideal case study for the understand-
ing of hadron formation in quantum chromodynamics (QCD) [1]. Its theoretical description
is based on the generally agreed assumption that the charm and beauty quarks (the heav-
iest ones capable of forming bound states) are heavy enough to allow the factorization of
short- and long-distance effects. Within the non-relativistic QCD (NRQCD) framework [2],
in particular, perturbative QCD computations provide the production cross sections of the
QQ pre-resonance (the “short-distance coefficients”, SDCs), while the non-perturbative
evolution of the QQ state to the observed meson (the hadronization step) is described by
phenomenological parameters (the “long-distance matrix elements”, LDMEs), determined
from fits to experimental data. Other theoretical approaches have been considered, such as
the colour-singlet model (CSM) [3,4] and the colour-evaporation model (CEM) [5,6]. These
theoretical models differ in the choice and classification of the allowed pre-resonance con-
figurations. The NRQCD approach foresees the contribution of all possible spin, S, orbital
angular momentum, L, total angular momentum, J, and colour (c= 1,8) configurations,
QQ(2S+1L[c]
J), organized in an expansion in powers of the relative QQ velocity, v < 1, so
that only a small number of leading and sub-leading terms remain quantitatively important.
Instead, the CSM considers that the final-state hadron can only result from a colour-neutral
(singlet) pre-resonance having the same quantum numbers and the CEM is built upon the
assumption that one universal hadronization factor per quarkonium state (independent of
the S, L, J configuration) multiplies the short-distance QQ production cross section.
The fundamental question that all models address is: how are the observable kinematic
properties of the produced quarkonium meson related to the quantum state of the unob-
servable QQ pre-resonance? The answers are different because, among other factors, the
several contributing short-distance processes are scaled by different long-distance weights.
The observable polarization of the quarkonium state provides particularly significant in-
formation regarding the hadronization model, given that it directly reflects the mixture
of S, L, J configurations (and polarizations) of the contributing pre-resonance states. The
polarizations of five vector quarkonia (J/ψ,ψ(2S), Υ(1S), Υ(2S) and Υ(3S)) have recently
been measured at relatively high transverse momentum, pT, both at the Tevatron [7] and
at the LHC [8–12]. These measurements, showing no significant signs of polarization, have
been addressed in many studies, including analyses based on the NRQCD [13–20] and
CEM [21] approaches.
In this paper we devote our attention to low-pTquarkonium hadro-production, a kine-
matical domain complementary to that explored at the LHC. We start by considering the
polarization measurements reported by several fixed-target experiments, at CERN, DESY
and Fermilab, using proton or pion beams, in a broad energy range, colliding on targets
made of several materials. The question we address here is: can this multitude of low-pT
quarkonium polarization measurements be interpreted in a consistent physical picture? At
first sight, we may think that it is very challenging to see coherent patterns emerging from
a collection of results obtained in such a diverse set of kinematical conditions, affected by
several difficulties in the detection and analysis techniques, and reported using three dif-
ferent polarization frames. Nevertheless, a careful look at the experimental results allows
us to see that, while most data points fluctuate around the unpolarized condition, there
are some tendencies towards strong polarizations in certain kinematical regions. These
qualitative patterns motivate us to consider a simple physical interpretation of low-pT
quarkonium production, as a superposition of two 2-to-1 processes: gluon-gluon fusion
2