There is another issue that should be fixed. In particular, when requiring e.g.,

p p > j j QED=0 QCD=4 [QED]

the code also includes QCD corrections (i.e., NLO_1), since at the squared order level, they are allowed. To remove them, one needs to set QCD=2.

In my opinion, the QED=0 and QCD=4 should be used to define which LO contributions should be included, and the ones between the square brackets to determine if one should include QCD and/or QED corrections on top of the Born contributions. Hence,

p p > j j QED=0 QCD=4 [QED]
and
p p > j j QED=0 QCD=2 [QED]

Hi,

There is another issue that should be fixed. In particular, when requiring e.g.,

p p > j j QED=0 QCD=4 [QED]

the code also includes QCD corrections (i.e., NLO_1), since at the squared order level, they are allowed. To remove them, one needs to set QCD=2.

In my opinion, the QED=0 and QCD=4 should be used to define which LO contributions should be included, and the ones between the square brackets to determine if one should include QCD and/or QED corrections on top of the Born contributions. Hence,

p p > j j QED=0 QCD=4 [QED]

and

p p > j j QED=0 QCD=2 [QED]

should lead to exactly the same results.

Who knows how to fix this?

Cheers,

Rikkert