Higgs pair production via gluon fusion in the Two-Higgs-Doublet Model

We study the production of Higgs boson pairs via gluon fusion at the LHC in the Two-Higgs-Doublet Model. We present predictions at NLO accuracy in QCD, matched to parton showers through the MC@NLO method. A dedicated reweighting technique is used to improve the NLO calculation upon the infinite top-mass limit. We perform our calculation within the MadGraph5_aMC@NLO framework, along with the 2HDM implementation based on the NLOCT package. The inclusion of the NLO corrections leads to large K-factors and significantly reduced theoretical uncertainties. We examine the seven 2HDM Higgs pair combinations using a number of representative 2HDM scenarios. We show how the model-specific features modify the Higgs pair total rates and distribution shapes, leading to trademark signatures of an extended Higgs sector.