Add to ORKGWe present results from a parton branching solution of the QCD evolution equations which includes the transverse momentum distributions for all parton flavours. The transverse momentum dependent (TMD) parton densities were extracted by fitting inclusive HERA DIS data with two alternative evolution scenarios. In both cases the same ordering variable is used, but as αS argument in the branching we use either the evolution scale or $p_T$. The resulting distributions differ mainly in the gluon component at low scales and/or small $k_T$ region. This was demonstrated by two phenomenological applications from LHC – the $Z$ boson $p_T$ spectrum and the azimuthal decorrelation of two leading jets