  We present the extraction of unpolarized transverse momentum dependent parton
distributions functions (TMDPDFs) and Collins-Soper kernel from the fit of
Drell-Yan and weak-vector boson production data. The TMDPDF are parameterized,
as commonly done, using their (large transverse momentum) asymptotic matching
to PDF. The analysis is done at the next-to-next-to-next-to-next-to leading
logarithmic accuracy (N$^4$LL) (performed only approximately because PDF
evolution is known so far at next-to-next-to leading order (NNLO)). The
non-perturbative model used for TMDPDF is flavor dependent to reduce the
colllinear PDF bias. The estimation of uncertainties is done with the replica
method and, for the first time, it includes the propagation of uncertainties
due to the collinear distributions.
