Efficiency improvements for the numerical computation of NLO corrections
  In this paper we discuss techniques, which lead to a significant improvement of the efficiency of the Monte Carlo integration, when one-loop QCD amplitudes are calculated numerically with the help of the subtraction method and contour deformation. The techniques discussed are: holomorphic and non-holomorphic division into sub-channels, optimisation of the integration contour, improvement of the ultraviolet subtraction terms, importance sampling and antithetic variates in loop momentum space, recurrence relations. 
Sebastian Becker, Christian Reuschle and Stefan Weinzierl
hep-ph