Implementing Transverse Momentum Dependent splitting functions in Parton Branching evolution equations

The Parton Branching (PB) approach describes the evolution of transverse momentum dependent (TMD) parton densities. We propose to extend the PB method by including TMD splitting functions, instead of the DGLAP splitting functions which assume strong ordering in transverse momentum. We present the evolution equations and their numerical solution, which is the first Monte Carlo implementation including TMD splitting functions.