ggxy: NLO QCD corrections to loop induced $gg$ initiated processes

We present the program package ggxy, which in its first version can be used to calculate partonic and hadronic cross sections to Higgs boson pair production at NLO QCD and it offers a direct interface to POWHEG-BOX. The 2-loop virtual amplitudes are implemented using analytical approximations in different kinematic regions, while all other parts of the calculation are exact. This implementation allows to freely modify all input parameter, such as the three Higgs coupling, the masses of the top quark and the Higgs boson, and the top-quark mass renormalization scheme.