  The phenomenology of heavy quarkonia production in hadron collisions is
reviewed. The theoretical predictions are compared to data. Commonly used
production models are shown to fail in explaining all the experimental
findings. The shortcomings of these models are analysed and possible
improvements are discussed.
