A non-linear EFT description of $gg\to HH$ at NLO interfaced to POWHEG

We present the implementation of Higgs boson pair production in gluon fusion within a non-linear Effective Field Theory framework containing five anomalous couplings for this process. The code, available within the POWHEG-BOX-V2, includes full NLO QCD corrections with massive top quarks. All five couplings can be modified by the user. We show m$_{hh}$ distributions at seven benchmark points provided by an m$_{hh}$ shape analysis at NLO and showered $ {p}_T^{hh} $ distributions resulting from an interface to PYTHIA-8 and HERWIG-7.2.