Masses, fermions and generalized D-dimensional unitarity
  We extend the generalized D-dimensional unitarity method for numerical evaluation of one-loop amplitudes by incorporating massive particles. The issues related to extending the spinor algebra to higher dimensions, treatment of external self-energy diagrams and mass renormalization are discussed within the context of the D-dimensional unitarity method. To validate our approach, we calculate in QCD the one-loop scattering amplitudes of a massive quark pair with up to three additional gluons for arbitrary spin states of the external quarks and gluons. 
R. K. Ellis, W. T. Giele, Z. Kunszt and K. Melnikov