  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.
