We predict the shape of the transverse momentum p_T spectrum of Upsilon
production. The distribution at low p_T is dominated by the region of small
impact parameter b and may be computed reliably in perturbation theory. We
resum to all orders in the strong coupling alpha_s the process-independent
large logarithmic contributions that arise from initial-state gluon showers in
the small p_T (< M_Upsilon) region. The cross section at large p_T is
represented by the alpha_s^3 lowest-order non-vanishing perturbative
contribution.Comment: 4 pages, 5 figures. Presented at the 2004 Meeting of the APS Division
of Particles and Fields, University of California, Riverside, CA, Aug 26-31,
200