  We present a first attempt at a global fit of unpolarized quark transverse
momentum dependent distribution and fragmentation functions from available data
on semi-inclusive deep-inelastic scattering, Drell-Yan and $Z$ boson production
processes. This analysis is performed in the low transverse momentum region, at
leading order in perturbative QCD and with the inclusion of energy scale
evolution effects at the next-to-leading logarithmic accuracy.
