In principle, one could escape the problem with the IR divergences altogether calculating the matrix element in question in Euclidian space, continuing analytically to Minkowki space and taking the imaginary part.