  In this talk, we discuss the speed-up of numerical calculations of jet
observables by replacing the usual sum over all helicity amplitudes with an
integral over a parametrisation for the parton polarisations called random
polarisations. Random polarisations are a linear combination of helicity
eigenstates multiplied by a phase factor depending on a so-called helicity
angle. Instead of a summation over discrete helicities, random polarisations
require an integration over the helicity angle. By combining this integral with
the final-state phase space integral, we only have to evaluate one squared
amplitude per phase space point instead of $2^n$ helicity amplitudes, where $n$
is the total number of particles in the process. While the technique itself has
been known since 1998, so far there has been no way of using it with dipole
subtraction, which is probably the most-used method for dealing with infrared
divergences in NLO calculations. After giving detailed reasons for this
statement, we propose a solution to this problem in terms of extending the
existing subtraction method by a new term.
