We calculate the transverse momentum p_T distribution for production of the
Upsilon states in hadronic reactions. For small pT(≤MΥ), we
resum to all orders in the strong coupling alpha_s the process-independent
large logarithmic contributions that arise from initial-state gluon showers. We
demonstrate that the p_T distribution at low p_T is dominated by the region of
small impact parameter b and that it may be computed reliably in perturbation
theory. We express the cross section at large p_T by the alpha_s^3 lowest-order
non-vanishing perturbative contribution. Our results are consistent with data
from the Fermilab Tevatron collider.Comment: 24 pages latex; 9 postscript files of figures. Presentation improved;
new figure and references added; conclusions unaltered. Version to be
published in Physical Review