  The factorization theorem for $q_T$ spectra in Drell-Yan processes, boson
production and semi-inclusive deep inelastic scattering allows for the
determination of the non-perturbative parts of transverse momentum dependent
parton distribution functions. Here we discuss the fit of Drell-Yan and
$Z$-production data using the transverse momentum dependent formalism and the
resummation of the evolution kernel. We find a good theoretical stability of
the results and a final $\chi^2/{\rm points}\lesssim 1$. We show how the fixing
of the non-perturbative pieces of the evolution can be used to make predictions
at present and future colliders.
