We consider higher-order QCD corrections to the production of colorless high-mass systems (lepton pairs, vector bosons, Higgs bosons, etc.) in hadron collisions. We propose a new formulation of the subtraction method to numerically compute arbitrary infrared-safe observables for this class of processes. To cancel the infrared divergences, we exploit the universal behavior of the associated transverse-momentum ($q_T$) distributions in the small-$q_T$ region. The method is illustrated in general terms up to the next-to-next-to-leading order in QCD perturbation theory. As a first explicit application, we study Higgs-boson production through gluon fusion. Our calculation is implemented in a parton level Monte Carlo program that includes the decay of the Higgs boson into two photons. We present selected numerical results at the CERN Large Hadron Collider.

• Received 5 March 2007

DOI:https://doi.org/10.1103/PhysRevLett.98.222002