General subtraction method for numerical calculation of one-loop QCD matrix elements

We present a subtraction scheme for eliminating the ultraviolet, soft, and collinear divergences in the numerical calculation of an arbitrary one-loop QCD amplitude with an arbitrary number of external legs. The subtractions consist of local counter terms in the space of the four-dimensional loop momentum. The ultraviolet subtraction terms reproduce MSbar renormalization. The key point in the method for the soft and collinear subtractions is that, although the subtraction terms are defined graph-by-graph and the matrix element is also calculated graph-by-graph, the sum over graphs of the integral of each the subtraction term can be evaluated analytically and provides the well known simple pole structure that arises from subtractions from real emission graphs, but with the opposite sign.