We consider semi-inclusive deep inelastic scattering (SIDIS) and Drell-Yan events within transverse momentum dependent (TMD) factorization. Based on the simultaneous fit of multiple data points, we extract the unpolarized TMD distributions and the non-perturbative evolution kernel. The high quality of the fit confirms a complete universality of TMD non-perturbative distributions. The extraction is supplemented by phenomenological analyses of various parts of the TMD factorization, such as sensitivity to non-perturbative parameterizations, perturbative orders, collinear distributions, correlations between parameters, and others.