The parton branching method is a formalism for the QCD evolution of transverse momentum dependent parton distribution functions and of a backward parton shower in the Cascade event generator. The parton evolution is iteratively calculated by treating resolvable branchings separately from non-resolvable branchings. With the implication of angular ordering in the scale of the strong coupling, in the soft-gluon resolution scale and in the calculation of transverse momenta in each branching, one performs soft-gluon resummation and avoids splittings close to the Landau pole. With the implementation of this method in a TMD evolution program updfevolv and MC event generator Cascade, predictions for observables in broad kinematic regimes can be made. Recent results show that the Z boson pT spectrum can be predicted accurately. This method has also been compared to the CSS theory for TMD factorization of the Drell Yan cross section. Both analytical and numerical results are presented.