Higher Order QCD Corrections to Tagged Production Processes
  We extend the phase space slicing method to allow for heavy quarks and fragmentation functions. The method can be used to calculate differential cross section in which any particular particle (massive or massless) is tagged. 
Stephane Keller (Fermilab) and Eric Laenen (CERN)