On stochastically sampling color configurations

Parton shower algorithms are key components of theoretical predictions for high-energy collider physics. Work towards more accurate parton shower algorithms is thus pursued along many different avenues. The systematic treatment of subleading color corrections in parton shower algorithms is however technically challenging and remains elusive. In this article, we present an efficient and numerically stable algorithm to sample color configurations at fixed $N_C=3$, using the correct color factor including subleading corrections with a parton shower. The algorithm is implemented as stand-alone program that can be interfaced to the PYTHIA event generator. Preliminary comparisons to to LEP data are presented.