  Collinear and transverse momentum dependent (TMD) parton densities are
obtained from fits to precision measurements of deep inelastic scattering (DIS)
cross sections at HERA. The parton densities are evolved by DGLAP evolution
with next-to-leading-order (NLO) splitting functions using the parton branching
method, allowing one to determine simultaneously collinear and TMD densities
for all flavors over a wide range in $x$, $\mu^2$ and $k_t$, relevant for
predictions at the LHC. The DIS cross section is computed from the parton
densities using perturbative NLO coefficient functions.
  Parton densities satisfying angular ordering conditions are presented. Two
sets of parton densities are obtained, differing in the renormalization scale
choice for the argument in the strong coupling alpha_s. This is taken to be
either the evolution scale $\mu$ or the transverse momentum $q_t$. While both
choices yield similarly good $\chi^2$ values for the fit to DIS measurements,
especially the gluon density turns out to differ between the two sets.
  The TMD densities are used to predict the transverse momentum spectrum of
Z-bosons at the LHC.
