“…Since we include parton masses in the dynamics of initial state splittings, we need evolution kernels for parton evolution that have non-trivial mass dependence. This view is in accord with the work of Jadach, Kusina, Placzek, Skrzypek and Slawinska [16,17], who connect parton shower splittings with the parton distribution function kernels with the aim of extending parton shower splittings to next-to-leading order in α s , albeit for massless partons. In this paper, we work only at lowest order, but include masses.…”