Combining Monte Carlo generators with next-to-next-to-leading order calculations: event reweighting for Higgs boson production at the LHC
Giovanna Davatz, Fabian Stoeckli, Charalampos Anastasiou, Guenther Dissertori, Michael Dittmar, Kirill Melnikov, Frank Petriello
TL;DR
This work proposes and tests a pragmatic method to incorporate NNLO QCD corrections into Monte Carlo event generators by bin-by-bin reweighting to NNLO Higgs production distributions in Higgs $p_T$ and rapidity. The approach is validated for $pp \to H+X$ using PYTHIA and MC@NLO, and extended to Higgs decays $H\to\gamma\gamma$ and $H\to W^+W^- \to l^+l^-\nu\bar{\nu}$, showing good agreement with NNLO predictions in regions where fixed-order or resummed calculations are reliable. Key findings include the necessity of NNLO-level reweighting to handle jet-veto observables, the robustness of lepton distributions against reweighting in the $WW$ channel, and the preservation of resummation features near kinematic boundaries. The study suggests that NNLO reweighting provides a practical and accurate way to include higher-order corrections in MC-based Higgs analyses at the LHC, with potential applicability to other processes.
Abstract
We study a phenomenological ansatz for merging next-to-next-to-leading order (NNLO) calculations with Monte Carlo event generators. We reweight them to match bin-integrated NNLO differential distributions. To test this procedure, we study the Higgs boson production cross-section at the LHC, for which a fully differential partonic NNLO calculation is available. We normalize PYTHIA and MC@NLO Monte Carlo events for Higgs production in the gluon fusion channel to reproduce the bin integrated NNLO double differential distribution in the transverse momentum and rapidity of the Higgs boson. These events are used to compute differential distributions for the photons in the pp \to H \to γγdecay channel, and are compared to predictions from fixed-order perturbation theory at NNLO. We find agreement between the reweighted generators and the NNLO result in kinematic regions where we expect a good description using fixed-order perturbation theory. Kinematic boundaries where resummation is required are also modeled correctly using this procedure. We then use these events to compute distributions in the pp \to H \to W^+W^- \to l^+l^- ν\barν channel, for which an accurate description is needed for measurements at the LHC. We find that the final state lepton distributions obtained from PYTHIA are not significantly changed by the reweighting procedure.
