The Transverse Momentum Dependent (TMD) Parton Branching (PB) method is a Monte Carlo (MC) approach to obtain QCD high energy collider predictions grounded in ideas originating from the TMD factorization. It provides an evolution equation for the TMD parton distribution functions (TMDs) and a framework to use those within TMD MC generators.

This work focuses on the structure of the PB Sudakov form factor. The Sudakov form factor is factorized in the perturbative and non-perturbative regions by introducing an intermediate separation scale motivated by angular ordering. The logarithmic order of the perturbative low-qt resummation achieved so far by the PB Sudakov is discussed by comparing it to the Collins-Soper-Sterman (CSS) method and is increased up to next-to-next-to-leading logarithm (NNLL) with the use of physical (effective) coupling. A non-perturbative Sudakov form factor provides a term analogous to Collins-Soper (CS) kernel. The effects of different evolution scenarios, including or not the non-perturbative Sudakov contribution, on a numerical extraction of the CS kernel are investigated.