  We examine the effect of the Z-boson transverse momentum distribution
measured at the Run-1 of the Tevatron on the nonperturbative function of the
Collins-Soper-Sterman (CSS) formalism, which resums large logarithmic terms
from multiple soft gluon emission in hadron collisions. The inclusion of the
Tevatron Run-1 Z-boson data strongly favors a Gaussian form of the CSS
nonperturbative function, when combined with the other low energy Drell-Yan
data in a global fit.
