Issue when generating one loop gluon fusion to higgs

Hello,

I am working with a 2HDM trying to generate g g > h2 > h1 h1. I am using the 2HDM_NLO UFO file that is available here http://feynrules.irmp.ucl.ac.be/wiki/2HDM . When I execute the following commands:
import model 2HDM_NLO
generate g g > h2 > h1 h1 [QCD]

I get the following error:

INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=6: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=7: WEIGTHED IS QCD+2*QED
Error detected in "generate g g > h2 > h1 h1 [QCD]"
write debug file MG5_debug
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str :
        Automatic process-order determination lead to negative constraints:
        QED: 6, QCD: -4
        Please specify the coupling orders from the command line.

I also tried generating g g > h using the loop_sm UFO. I get a similar error:
import model loop_sm
MG5_aMC>generate g g > h [QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS QCD+2*QED
Error detected in "generate g g > h [QCD]"
write debug file MG5_debug
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str :
        Automatic process-order determination lead to negative constraints:
        QED: 3, QCD: -2
        Please specify the coupling orders from the command line.

I am using MG5_aMC_v3_4_2.

Is there a way to get madgraph to perform this calculation?

Thank you,
Artemis