IOError : [Errno 2] No such file or directory:

Hi,

Can you look at the log file of the following directory:
> /home/lihuang/tools/packages/INSTALLATIONS/MG5_aMC_v2_6_2/result/sumrule/backgorund_100Tev/pp_hjjj/SubProcesses/P1_gg_hgbbx/G14/

It should explain what is the reason of the crash

Cheers,

Olivier

> On 1 Jun 2018, at 20:27, Li <email address hidden> wrote:
>
> Question #669848 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/669848
>
> Description changed to:
> Hi,
>
> I met a error when I use:
>
> import model heft
> define jb = j b b~
> generate p p > h jb jb jb
>
>
> Then, ( The other setting like ptj seems irrelevant to this problem )
>
> Error when reading /home/lihuang/tools/packages/INSTALLATIONS/MG5_aMC_v2_6_2/result/sumrule/backgorund_100Tev/pp_hjjj/SubProcesses/P1_gg_hgbbx/G14/results.dat
> Command "survey" interrupted with error:
> IOError : [Errno 2] No such file or directory: '/home/lihuang/tools/packages/INSTALLATIONS/MG5_aMC_v2_6_2/result/sumrule/backgorund_100Tev/pp_hjjj/SubProcesses/P1_gg_hgbbx/G14/results.dat'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/home/lihuang/tools/packages/INSTALLATIONS/MG5_aMC_v2_6_2/result/sumrule/backgorund_100Tev/pp_hjjj/run_01_tag_1_debug.log'.
> Please attach this file to your report.
>
> I tried with version 2.6.1 and 2.6.2 and different PCs. Is that a bug
> and is there any way to avoid it?
>
> Best,
> Li
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

    Find two errors but I don't understand them.

 Process in group number 1
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23_lo_as_0130_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23lo1: 0.13000000000000000
 Define smin to 400.00000000000000
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

  External Params
  ---------------------------------

 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.2000000000000002
 mdl_ymt = 164.50000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_MP = 125.00010000000000
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_WH = 6.3823389999999999E-003
 mdl_WH1 = 6.3823389999999999E-003
  Internal Params
  ---------------------------------

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__4 = 244140625.00000000
 mdl_MT__exp__4 = 895745041.00000000
 mdl_MH__exp__2 = 15625.000000000000
 mdl_MT__exp__2 = 29929.000000000000
 mdl_MH__exp__12 = 1.4551915228366852E+025
 mdl_MH__exp__10 = 9.3132257461547853E+020
 mdl_MH__exp__8 = 59604644775390624.
 mdl_MH__exp__6 = 3814697265625.0000
 mdl_MT__exp__6 = 26808753332089.000
 mdl_aEW = 7.5467711139788835E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.6872153846781555E-002
 mdl_ee = 0.30795376724436879
 mdl_MW__exp__2 = 6467.2159543705357
 mdl_sw2 = 0.22224648578577766
 mdl_cw = 0.88190334743339216
 mdl_sqrt__sw2 = 0.47143025548407230
 mdl_sw = 0.47143025548407230
 mdl_g1 = 0.34919219678733299
 mdl_gw = 0.65323293034757990
 mdl_v = 246.21845810181637
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_MW__exp__12 = 7.3165135181362573E+022
 mdl_MW__exp__10 = 1.1313235200058173E+019
 mdl_MW__exp__8 = 1749320771082757.2
 mdl_MW__exp__6 = 270490545456.51419
 mdl_MW__exp__4 = 41824882.200464800
 mdl_AH = 3.1606819475643554E-005
 mdl_v__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yb = 2.4123686777011714E-002
 mdl_yt = 0.94484439876629200
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_sw__exp__2 = 0.22224648578577769
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_GH = -5.2486599934244080E-005
 mdl_Gphi = -7.9843646768522597E-005
  Couplings of heft
  ---------------------------------

        GC_10 = -0.12177E+01 0.00000E+00
        GC_11 = 0.00000E+00 0.12177E+01
        GC_12 = 0.00000E+00 0.14828E+01
        GC_13 = 0.00000E+00 0.52487E-04
        GC_14 = 0.63914E-04 0.00000E+00
        GC_15 = -0.00000E+00 -0.77829E-04
        GC_74 = -0.00000E+00 -0.17058E-01

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23lo1
 alpha_s(Mz)= 0.1300 running at 2 loops.
 alpha_s(Mz)= 0.1300 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 1000 5 3
Enter desired fractional accuracy:
 Desired fractional accuracy: 0.10000000000000001
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Explicitly summing over helicities
Enter Configuration Number:
Running Configuration Number: 14
 Not subdividing B.W.
 Attempting mappinvarients 1 6
 Completed mapping 6
 about to integrate 10 1000 5 3 10 1
 Using non-zero grid deformation.
 10 dimensions 1000 events 10 invarients 5 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 24
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 13000.00
 Mapping Graph 14 to config 14
Setting grid 1 0.52284E-06 1
Setting grid 2 0.52284E-06 1
Setting grid 3 0.10688E-03 1
 Transforming s_hat 1/s 9 1.4106130177514792E-004 400.00000000000000 168999999.99999997
  14 1 2 3 4 5 6 7 8 9 10
 Masses: 0.000E+00 0.000E+00 0.125E+03 0.000E+00 0.470E+01 0.470E+01
Using random seed offsets 14 : 2
  with seed 24
 Ranmar initialization seeds 31072 9399

 ********************************************
 * You are using the DiscreteSampler module *
 * part of the MG5_aMC framework *
 * Author: Valentin Hirschi *
 ********************************************

  Particle 3 4 5 6
      Et > 0.0 20.0 0.0 0.0
       E > 0.0 0.0 0.0 0.0
     Eta < -1.0 5.0 -1.0 -1.0
   xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.0 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 0.0 0.0
s min # 4> 0.0 0.0 0.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
 alpha_s for scale 276.49796875828514 is 0.10976171725002294
 Added good helicity 1 0.73363657441193242 in event 1 local: 1
 Added good helicity 2 1.1715040797258971 in event 1 local: 1
 Added good helicity 3 3.2234632098845744 in event 1 local: 1
 Added good helicity 4 1.3206935296574118 in event 1 local: 1
 Added good helicity 5 0.28570467333990651 in event 1 local: 1
 Added good helicity 6 1.3444869877660892 in event 1 local: 1
 Added good helicity 7 1.2819594797879881 in event 1 local: 1
 Added good helicity 8 0.54460517257579055 in event 1 local: 1
 Added good helicity 9 0.27797584837047701 in event 1 local: 1
 Added good helicity 10 0.37716791638585012 in event 1 local: 1
 Added good helicity 11 3.1196265049066407 in event 1 local: 1
 Added good helicity 12 3.0547796870709661E-005 in event 1 local: 1
 Added good helicity 13 2.3980646118471671E-002 in event 1 local: 1
 Added good helicity 14 0.77633170582150579 in event 1 local: 1
 Added good helicity 15 1.5187533486073008 in event 1 local: 1
 Added good helicity 16 7.9774843292364011E-005 in event 1 local: 1
 Added good helicity 17 7.9774843292363971E-005 in event 1 local: 1
 Added good helicity 18 1.5187533486073008 in event 1 local: 1
 Added good helicity 19 0.77633170582150579 in event 1 local: 1
 Added good helicity 20 2.3980646118471671E-002 in event 1 local: 1
 Added good helicity 21 3.0547796870709675E-005 in event 1 local: 1
 Added good helicity 22 3.1196265049066390 in event 1 local: 1
 Added good helicity 23 0.37716791638585012 in event 1 local: 1
 Added good helicity 24 0.27797584837047712 in event 1 local: 1
 Added good helicity 25 0.54460517257579055 in event 1 local: 1
 Added good helicity 26 1.2819594797879881 in event 1 local: 1
 Added good helicity 27 1.3444869877660888 in event 1 local: 1
 Added good helicity 28 0.28570467333990651 in event 1 local: 1
 Added good helicity 29 1.3206935296574118 in event 1 local: 1
 Added good helicity 30 3.2234632098845744 in event 1 local: 1
 Added good helicity 31 1.1715040797258967 in event 1 local: 1
 Added good helicity 32 0.73363657441193242 in event 1 local: 1
 Error: Failed despite same graph: 14
 Have jets (>0) 0 0 0 1 0 0
 Should be

ls status:
input_app.txt
run1_app.log

Best,
Li

Hi Olivier,

    I find that for version 3.0.0, it said:

 P1_qq_hgbbx SubProcesses doesn't have available phase-space.
            Please check mass spectrum.

    This information really confused me. If P1_gg_hgbbx subprocess has phase-space, why P1_qq_hgbbx does not?

Best,
Li

Thanks Olivier Mattelaer, that solved my question.

Hi,

I have put a patch on the other thread and it seems that your process is not working
if you do not put some decent cut on your b quark.
If you put set maxjetfavour to 5 then it works nicely.

I will invetigate this issue later but I doubt to have time on Friday, so this will be for next Wednesday.

Cheers,

Olivier

> On 1 Jun 2018, at 21:12, Li <email address hidden> wrote:
>
> Question #669848 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/669848
>
> Status: Answered => Open
>
> Li is still having a problem:
> Hi Olivier,
>
> Find two errors but I don't understand them.
>
> Process in group number 1
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> ****************************************
>
> NNPDFDriver version 1.0.3
> Grid: NNPDF23_lo_as_0130_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23lo1: 0.13000000000000000
> Define smin to 400.00000000000000
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.2000000000000002
> mdl_ymt = 164.50000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MZ = 91.188000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_MP = 125.00010000000000
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WH = 6.3823389999999999E-003
> mdl_WH1 = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__4 = 244140625.00000000
> mdl_MT__exp__4 = 895745041.00000000
> mdl_MH__exp__2 = 15625.000000000000
> mdl_MT__exp__2 = 29929.000000000000
> mdl_MH__exp__12 = 1.4551915228366852E+025
> mdl_MH__exp__10 = 9.3132257461547853E+020
> mdl_MH__exp__8 = 59604644775390624.
> mdl_MH__exp__6 = 3814697265625.0000
> mdl_MT__exp__6 = 26808753332089.000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_v = 246.21845810181637
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_MW__exp__12 = 7.3165135181362573E+022
> mdl_MW__exp__10 = 1.1313235200058173E+019
> mdl_MW__exp__8 = 1749320771082757.2
> mdl_MW__exp__6 = 270490545456.51419
> mdl_MW__exp__4 = 41824882.200464800
> mdl_AH = 3.1606819475643554E-005
> mdl_v__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.4123686777011714E-002
> mdl_yt = 0.94484439876629200
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_gw__exp__2 = 0.42671326129048615
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_sw__exp__2 = 0.22224648578577769
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> mdl_GH = -5.2486599934244080E-005
> mdl_Gphi = -7.9843646768522597E-005
> Couplings of heft
> ---------------------------------
>
> GC_10 = -0.12177E+01 0.00000E+00
> GC_11 = 0.00000E+00 0.12177E+01
> GC_12 = 0.00000E+00 0.14828E+01
> GC_13 = 0.00000E+00 0.52487E-04
> GC_14 = 0.63914E-04 0.00000E+00
> GC_15 = -0.00000E+00 -0.77829E-04
> GC_74 = -0.00000E+00 -0.17058E-01
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = nn23lo1
> alpha_s(Mz)= 0.1300 running at 2 loops.
> alpha_s(Mz)= 0.1300 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> getting user params
> Enter number of events and max and min iterations:
> Number of events and iterations 1000 5 3
> Enter desired fractional accuracy:
> Desired fractional accuracy: 0.10000000000000001
> Enter 0 for fixed, 2 for adjustable grid:
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Explicitly summing over helicities
> Enter Configuration Number:
> Running Configuration Number: 14
> Not subdividing B.W.
> Attempting mappinvarients 1 6
> Completed mapping 6
> about to integrate 10 1000 5 3 10 1
> Using non-zero grid deformation.
> 10 dimensions 1000 events 10 invarients 5 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 24
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 13000.00
> Mapping Graph 14 to config 14
> Setting grid 1 0.52284E-06 1
> Setting grid 2 0.52284E-06 1
> Setting grid 3 0.10688E-03 1
> Transforming s_hat 1/s 9 1.4106130177514792E-004 400.00000000000000 168999999.99999997
> 14 1 2 3 4 5 6 7 8 9 10
> Masses: 0.000E+00 0.000E+00 0.125E+03 0.000E+00 0.470E+01 0.470E+01
> Using random seed offsets 14 : 2
> with seed 24
> Ranmar initialization seeds 31072 9399
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Particle 3 4 5 6
> Et > 0.0 20.0 0.0 0.0
> E > 0.0 0.0 0.0 0.0
> Eta < -1.0 5.0 -1.0 -1.0
> xqcut: 0.0 0.0 0.0 0.0
> d R # 3 > -0.0 0.0 0.0 0.0
> d R # 4 > -0.0 -0.0 0.0 0.0
> d R # 5 > -0.0 -0.0 -0.0 0.0
> s min # 3> 0.0 0.0 0.0 0.0
> s min # 4> 0.0 0.0 0.0 0.0
> s min # 5> 0.0 0.0 0.0 0.0
> xqcutij # 3> 0.0 0.0 0.0 0.0
> xqcutij # 4> 0.0 0.0 0.0 0.0
> xqcutij # 5> 0.0 0.0 0.0 0.0
> alpha_s for scale 276.49796875828514 is 0.10976171725002294
> Added good helicity 1 0.73363657441193242 in event 1 local: 1
> Added good helicity 2 1.1715040797258971 in event 1 local: 1
> Added good helicity 3 3.2234632098845744 in event 1 local: 1
> Added good helicity 4 1.3206935296574118 in event 1 local: 1
> Added good helicity 5 0.28570467333990651 in event 1 local: 1
> Added good helicity 6 1.3444869877660892 in event 1 local: 1
> Added good helicity 7 1.2819594797879881 in event 1 local: 1
> Added good helicity 8 0.54460517257579055 in event 1 local: 1
> Added good helicity 9 0.27797584837047701 in event 1 local: 1
> Added good helicity 10 0.37716791638585012 in event 1 local: 1
> Added good helicity 11 3.1196265049066407 in event 1 local: 1
> Added good helicity 12 3.0547796870709661E-005 in event 1 local: 1
> Added good helicity 13 2.3980646118471671E-002 in event 1 local: 1
> Added good helicity 14 0.77633170582150579 in event 1 local: 1
> Added good helicity 15 1.5187533486073008 in event 1 local: 1
> Added good helicity 16 7.9774843292364011E-005 in event 1 local: 1
> Added good helicity 17 7.9774843292363971E-005 in event 1 local: 1
> Added good helicity 18 1.5187533486073008 in event 1 local: 1
> Added good helicity 19 0.77633170582150579 in event 1 local: 1
> Added good helicity 20 2.3980646118471671E-002 in event 1 local: 1
> Added good helicity 21 3.0547796870709675E-005 in event 1 local: 1
> Added good helicity 22 3.1196265049066390 in event 1 local: 1
> Added good helicity 23 0.37716791638585012 in event 1 local: 1
> Added good helicity 24 0.27797584837047712 in event 1 local: 1
> Added good helicity 25 0.54460517257579055 in event 1 local: 1
> Added good helicity 26 1.2819594797879881 in event 1 local: 1
> Added good helicity 27 1.3444869877660888 in event 1 local: 1
> Added good helicity 28 0.28570467333990651 in event 1 local: 1
> Added good helicity 29 1.3206935296574118 in event 1 local: 1
> Added good helicity 30 3.2234632098845744 in event 1 local: 1
> Added good helicity 31 1.1715040797258967 in event 1 local: 1
> Added good helicity 32 0.73363657441193242 in event 1 local: 1
> Error: Failed despite same graph: 14
> Have jets (>0) 0 0 0 1 0 0
> Should be
>
> ls status:
> input_app.txt
> run1_app.log
>
>
> Best,
> Li
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

I finally decide to put this in higher priority and re-schedule my program since I would need to run a large validation set of run for the changes of the other bugs anyway, so better to combine those validation with the one that are needed for this one.

So here is the patch that I'm setting inside the development branch and which solve the issue that you face.
(Now I have to check if they are no side effect)

Cheers,

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-08 09:36:16 +0000
@@ -1,3 +1,20 @@
+c for cross-checking change in this file.
+c here is a minimal list of process that we have to
+c test
+c
+c SM
+c ----
+c p p > t t~ (up to 2jet)
+c p p > w+ (up to 3 jet)
+c p p > z t t~ j j (no MLM needed)
+c
+c
+c HEFT
+c ----
+c p p > h j b b~
+c q q > a a g q q
+c g g > h g q q
+
       double precision function gamma(q0)
 c**************************************************
 c calculates the branching probability
@@ -522,7 +539,8 @@
       common/to_stot/stot,m1,m2

 C local variables
- integer i, j, idi, idj, k
+ integer i, j, idi, idj, k,m
+ integer get_color
       real*8 PI
       parameter( PI = 3.14159265358979323846d0 )
       integer iforest(2,-max_branch:-1,lmaxconfigs)
@@ -545,6 +563,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 +765,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
+ ! 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
@@ -791,7 +819,57 @@
            goodjet(imocl(n))=
      $ (isjet(ipdgcl(imocl(n),igraphs(1),iproc)).and.
      $ goodjet(idacl(n,1)).and.goodjet(idacl(n,2)))
- endif
+
+c check case with g > g g
+c where the hardest gluon is not goodjet but the other is.
+c in that case change ipart(1,) of the mother gluon
+c pure QCD jet
+c need to take care of the following case:
+c tttttttttt
+c gggggggggggggggg
+c ggggggggg tttttttttt
+c gggggggg
+c
+c in that case the up gluon can be tag as the hardest one
+c but this one is also lead to no QCD one.
+c so in that case we have to change ipart(1) to the sofest gluon
+ pdgm = ipdgcl(imocl(n),igraphs(1),iproc)
+ pdgid1 = ipdgcl(idacl(n,1),igraphs(1),iproc)
+ pdgid2 = ipdgcl(idacl(n,2),igraphs(1),iproc)
+ if (is_octet(pdgm).and.is_octet(pdgid1).and.is_octet(pdgid2))then
+c write(*,*) 'pure QCD vertex (2)'
+c write(*,*) pdgm , '>', pdgid1,' ', pdgid2
+c write(*,*) 'ipart', ipart(1,imocl(n)), ipart(2,imocl(n))
+c write(*,*) 'id', imocl(n), idacl(n,1),idacl(n,2)
+c write(*,*) 'ipart',ipart(1,imocl(n)),'/', ipart(2,imocl(n)), ipart(1,idacl(n,1)),'/', ipart(2,idacl(n,1)),
+c $ ipart(1,idacl(n,2)),'/', ipart(2,idacl(n,2))
+c write(*,*) 'googjet', goodjet(idacl(n,1)),goodjet(idacl(n,2))
+ if (ipart(1,imocl(n)).eq.ipart(1, idacl(n,1))) then
+ if (.not.goodjet(idacl(n,1)).and.goodjet(idacl(n,2))) then
+c write(*,*) 'ggg with hard jet set a QED the second jet lead to', ipart(1,idacl(n,2)), ipart(2,idacl(n,2))
+ do m =n_max_cl,n,-1
+ if(ipart(1,m).eq.ipart(1,imocl(n)).and.ipart(2,m).eq.ipart(2,imocl(n)))then
+ ipart(1,m) = ipart(1,idacl(n,2))
+ ipart(2,m) = ipart(2,idacl(n,2))
+ endif
+ enddo
+ endif
+ else
+ if (.not.goodjet(idacl(n,2)).and.goodjet(idacl(n,1))) then
+c write(*,*) 'ggg with hard jet set a QED the second jet lead to', ipart(1,idacl(n,1)), ipart(2,idacl(n,1))
+ do m =n_max_cl,n,-1
+ if(ipart(1,m).eq.ipart(1,imocl(n)).and.ipart(2,m).eq.ipart(2,imocl(n)))then
+ ipart(1,m) = ipart(1,idacl(n,1))
+ ipart(2,m) = ipart(2,idacl(n,1))
+ endif
+ enddo
+ endif
+ endif
+ endif
+
+
+
+ endif
       enddo

       if (btest(mlevel,4))then