Hi, Sorry I'm one day late on this thread. The problem is linked to the dynamical scale choice. So one easy way to fix the problem is to not use the default dynamical scale choice but a simpler one (like sum HT/2/...) Another solution is to use the following patch. This is not yet fully tested (such change can have huge impact on MLM merging and i have to be carefull about that.) but it solve the problem of the directory that you point out. Thanks a lot for such report, Olivier === modified file 'Template/LO/SubProcesses/reweight.f' --- Template/LO/SubProcesses/reweight.f 2018-03-11 12:29:51 +0000 +++ Template/LO/SubProcesses/reweight.f 2018-06-07 19:04:43 +0000 @@ -545,6 +545,7 @@ logical chclusold,fail,increasecode save chclusold integer tmpindex + integer pdgm, pdgid1, pdgid2 logical isqcd,isjet,isparton,cluster,isjetvx,is_octet integer ifsno @@ -746,8 +747,17 @@ c The ishft gives the FS particle corresponding to imocl if(.not.is_octet(ipdgcl(ishft(1,ipart(1,imocl(n))-1),igraphs(1),iproc)))then ! split case for q a > q and for g > g h (with the gluon splitting into quark) - if (ipart(2,imocl(n)).eq.0) then ! q a > q case - iqjets(ipart(1,imocl(n)))=0 +c ! also check for case of three scalar interaction (then do nothing) + pdgm = ipdgcl(imocl(n),igraphs(1),iproc) + pdgid1 = ipdgcl(idacl(n,1),igraphs(1),iproc) + pdgid2 = ipdgcl(idacl(n,2),igraphs(1),iproc) + + if (.not.isqcd(pdgm).and..not.isqcd(pdgid1).and..not.isqcd(pdgid2)) then + ! this is to avoid to do weird stuff for w+ w- z (or h h h) + ! this fix an issue for qq_zttxqq G1594.08 + continue + elseif (ipart(2,imocl(n)).eq.0) then ! q a > q case + iqjets(ipart(1,imocl(n)))=0 else ! octet. want to be sure that both are tagged as jet before removing one ! this prevent that both are removed in case of g > g h , g > q1 q2, q1 > a q1. ! at least one of the two should be kept as jet > On 4 Jun 2018, at 19:11, jackaraz