A new approach to QCD final-state evolution in processes with massive partons

We present an algorithm for massive parton evolution which is based on the differentially accurate simulation of soft-gluon radiation by means of a non-trivial azimuthal angle dependence of the splitting functions. The kinematics mapping is chosen such as to to reflect the symmetry of the final state in soft-gluon radiation and collinear splitting processes. We compute the counterterms needed for a fully differential NLO matching and discuss the analytic structure of the parton shower in the NLL limit. We implement the new algorithm in the numerical code Alaric and present a first comparison to experimental data.