NNLOCAL: completely local subtractions for color-singlet production in hadron collisions

We present NNLOCAL, a proof-of-concept parton-level Monte Carlo program implementing the extension of the completely local subtraction scheme CoLoRFulNNLO to the case of color-singlet production in hadron collisions. We have built general local subtraction terms that regularize all single and double unresolved infrared singularities in real radiation phase space. The subtractions are then integrated fully analytically to the required order in the parameter of dimensional regularization. Combining the integrated counterterms with the virtual contributions we demonstrate the cancellation of all infrared poles explicitly. We validate our procedure by computing the fully differential cross section for the production of a Higgs boson at the LHC in an effective field theory with gluons only. Our code provides the first public implementation of a completely local analytic subtraction scheme at next-to-next-to-leading order accuracy.