We present an 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 (N4LL) (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.