  We present analytic formulas that describe fully-differential production of
color-singlet final states in $q\bar q$ and $gg$ annihilation, including all
the relevant partonic channels, through NNLO QCD. We work within the nested
soft-collinear scheme which allows for fully local subtraction of infrared
divergences. We demonstrate analytic cancellation of soft and collinear poles
and present formulas for finite parts of all integrated subtraction terms.
These results provide an important building block for calculating NNLO QCD
corrections to arbitrary processes at hadron colliders within the nested
soft-collinear subtraction scheme.
