Interactions between hard partons and the quark-gluon plasma range from frequent soft interactions to rare hard scatterings. The larger number of soft interactions makes possible an effective stochastic description of parton-plasma interactions in terms of drag and diffusion transport coefficients. In this work, we present a numerical implementation that builds upon this systematic division between soft and hard parton-plasma interactions. We study the dependence of the single parton distribution on the cutoff between soft and hard parton-plasma interactions, both for small and phenomenological values of the strong coupling constant.