MadGraph5_aMC@NLO

S -> g g g end with error

Asked by Sarif Khan

Hi,

I am getting the following error if I include 3 gluon decay modes in the final state. I tried but could not solve it after looking at the other similar issues at the launchpad. Can you please let me know if I need to change something?

##########

Error when reading /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/G1/results.dat
Command "generate_events run_01" interrupted with error:
Exception : Reported error: End code 2.0
  Full associated log:
  Process in group number 1
  No PDF is used, alpha_s(MZ) from param_card is used
  Value of alpha_s from param_card: 0.13000000000000003
  The default order of alpha_s running is fixed to 2
  Warning! ptj set to xqcut= 30.000000000000000 to improve integration efficiency
  Note that this might affect non-radiated jets,
  e.g. from decays. Use cut_decays=F in run_card.
  Warning! mmjj set to xqcut= 30.000000000000000 to improve integration efficiency
  Note that this might affect non-radiated jets,
  e.g. from decays. Use cut_decays=F in run_card.
  Define smin to 3600.0000000000000
  *****************************************************
  * MadGraph/MadEvent *
  * -------------------------------- *
  * http://madgraph.hep.uiuc.edu *
  * http://madgraph.phys.ucl.ac.be *
  * http://madgraph.roma2.infn.it *
  * -------------------------------- *
  * *
  * PARAMETER AND COUPLING VALUES *
  * *
  *****************************************************

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

  mdl_cabi = 0.22773599999999999
  aEWM1 = 127.90000000000001
  mdl_Gf = 1.1663700000000000E-005
  aS = 0.13000000000000003
  mdl_ymdo = 5.0400000000000002E-003
  mdl_ymup = 2.5500000000000002E-003
  mdl_yms = 0.10100000000000001
  mdl_ymc = 0.69999999999999996
  mdl_ymb = 3.1499999999999999
  mdl_ymt = 172.00000000000000
  mdl_yme = 5.1099999999999995E-004
  mdl_ymm = 0.10566000000000000
  mdl_ymtau = 1.7769999999999999
  mdl_xpow = 100.00000000000000
  mdl_yDM5 = 4.1099999999999998E-020
  mdl_yDM52 = 2.0550000000000001E-019
  mdl_yDM53 = 4.1100000000000002E-019
  mdl_yDM54 = 1.0270000000000000E-019
  mdl_th = 0.78500000000000003
  mdl_MZ = 91.187600000000003
  mdl_Me = 5.1099999999999995E-004
  mdl_MMU = 0.10566000000000000
  mdl_MTA = 1.7769999999999999
  mdl_MU = 2.5500000000000002E-003
  mdl_MC = 0.69999999999999996
  mdl_MT = 172.00000000000000
  mdl_MD = 5.0400000000000002E-003
  mdl_MS = 0.10100000000000001
  mdl_MB = 3.1499999999999999
  mdl_MH = 125.00000000000000
  mdl_WZ = 2.4952000000000001
  mdl_WW = 2.0850000000000000
  mdl_WT = 1.5083359999999999
  mdl_WH = 4.0699999999999998E-003
  mdl_Wms = 4.0699999999999998E-003
   Internal Params
   ---------------------------------

  mdl_Mms = 100.00000000000000
  mdl_cos__cabi = 0.97418004031982097
  mdl_CKM1x1 = (0.97418004031982097,0.0000000000000000)
  mdl_sin__cabi = 0.22577256042856930
  mdl_CKM1x2 = (0.22577256042856930,0.0000000000000000)
  mdl_CKM1x3 = (0.0000000000000000,0.0000000000000000)
  mdl_CKM2x1 = (-0.22577256042856930,0.0000000000000000)
  mdl_CKM2x2 = (0.97418004031982097,0.0000000000000000)
  mdl_CKM2x3 = (0.0000000000000000,0.0000000000000000)
  mdl_CKM3x1 = (0.0000000000000000,0.0000000000000000)
  mdl_CKM3x2 = (0.0000000000000000,0.0000000000000000)
  mdl_CKM3x3 = (1.0000000000000000,0.0000000000000000)
  mdl_MZ__exp__2 = 8315.1783937600012
  mdl_MZ__exp__4 = 69142191.720053151
  mdl_sqrt__2 = 1.4142135623730951
  mdl_MH__exp__2 = 15625.000000000000
  mdl_complexi = (0.0000000000000000,1.0000000000000000)
  mdl_conjg__CKM1x1 = (0.97418004031982097,-0.0000000000000000)
  mdl_conjg__CKM1x2 = (0.22577256042856930,-0.0000000000000000)
  mdl_conjg__CKM1x3 = (0.0000000000000000,-0.0000000000000000)
  mdl_conjg__CKM2x1 = (-0.22577256042856930,-0.0000000000000000)
  mdl_conjg__CKM2x2 = (0.97418004031982097,-0.0000000000000000)
  mdl_conjg__CKM2x3 = (0.0000000000000000,-0.0000000000000000)
 At line 1235 of file reweight.f
 Fortran runtime error: Index '32681' of dimension 1 of array 'ptclus' above upper bound of 3

 Error termination. Backtrace:
 #0 0x7fa9f28238a2 in ???
 #1 0x7fa9f2824399 in ???
 #2 0x7fa9f2824929 in ???
 #3 0x41dc4f in setclscales_
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/reweight.f:1235
 #4 0x415c4b in passcuts_
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/cuts.f:1236
 #5 0x4252c6 in dsigproc_
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/auto_dsig.f:505
 #6 0x426459 in dsig_
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/auto_dsig.f:242
  mdl_conjg__CKM3x1 = (0.0000000000000000,-0.0000000000000000)
  mdl_conjg__CKM3x2 = (0.0000000000000000,-0.0000000000000000)
  mdl_conjg__CKM3x3 = (1.0000000000000000,-0.0000000000000000)
  mdl_aEW = 7.8186082877247844E-003
  mdl_MW = 79.824359746197842
  mdl_sqrt__aEW = 8.8422894590285753E-002
  mdl_ee = 0.31345100004952897
  mdl_MW__exp__2 = 6371.9284088904105
  mdl_sw2 = 0.23369913342182447
  mdl_cw = 0.87538612427783857
  mdl_sqrt__sw2 = 0.48342438232036300
  mdl_sw = 0.48342438232036300
  mdl_g1 = 0.35807170271074895
  mdl_gw = 0.64839716719502682
  mdl_vev = 246.22056907348590
  mdl_vev__exp__2 = 60624.568634871241
  mdl_lam = 0.12886689630821144
  mdl_yb = 1.8092609964465226E-002
  mdl_yc = 4.0205799921033836E-003
  mdl_ydo = 2.8948175943144366E-005
  mdl_ye = 2.9350233942354699E-006
  mdl_ym = 6.0687783137949077E-004
  mdl_ys = 5.8011225600348830E-004
  mdl_yt = 0.98791394091683138
  mdl_ytau = 1.0206529494239589E-002
  mdl_yup = 1.4646398542662327E-005
  mdl_muH = 88.388347648318444
  mdl_ee__exp__2 = 9.8251529432049817E-002
  mdl_sw__exp__2 = 0.23369913342182450
  mdl_cw__exp__2 = 0.76630086657817542
  mdl_muH__exp__2 = 7812.5000000000009
   Internal Params evaluated point by point
   ----------------------------------------

  mdl_sqrt__aS = 0.36055512754639896
  mdl_G__exp__2 = 1.6336281798666927
   Couplings of StandardModelDM5_UFO
   ---------------------------------

          GC_6 = -0.12781E+01 0.00000E+00
         GC_73 = -0.52531E-18 0.00000E+00
         GC_70 = -0.00000E+00 -0.41100E-18

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

  Running at e e machine @ 100.00000000000000 GeV
  PDF set = none
  alpha_s(Mz)= 0.1300 running at 2 loops.
  alpha_s(Mz)= 0.1300 running at 2 loops.
  Renormalization scale fixed @ 100.00000000000000
  Factorization scales fixed @ 91.188000000000002 91.188000000000002

  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: 1
  Not subdividing B.W.
  Attempting mappinvarients 1 3
  Determine nb_t
  T-channel found: 0
  Completed mapping 3
  about to integrate 2 1000 5 3 2 1
  Using non-zero grid deformation.
   2 dimensions 1000 events 2 invarients 5 iterations 1 config(s), (0.99)
  Using h-tuple random number sequence.
  Error opening grid
  Using Uniform Grid! 16
  Using uniform alpha 1.0000000000000000
  Grid defined OK
  Masses: 0.100E+03 0.000E+00 0.000E+00
  Set CM energy to 100.00
  Mapping Graph 1 to config 1
  Determine nb_t
  T-channel found: 0
 Using random seed offsets 1 : 2
   with seed 27
  Ranmar initialization seeds 4635 9403

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

   Particle 2 3
       Et > 30.0 30.0
        E > 0.0 0.0
      Eta < -1.0 -1.0
    xqcut: 0.0 30.0
 d R # 2 > -0.0 0.0
 s min # 2> 0.0 900.0
 xqcutij # 2> 0.0 0.0
 #7 0x437991 in sample_full_
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/Source/dsample.f:137
 #8 0x40448c in driver
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/driver.f:184
 #9 0x404618 in main
  at /scratch/khan24/package/MG5_aMC_v3_4_1/PROC_StandardModelDM5_UFO_2/SubProcesses/P1_S_gg/driver.f:204

##########

Best regards,
Sarif Khan

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said : 		#1

So looks like you do include two gluon and three gluon.
The issue is that due QCD renormalization you do have a double counting between the two if done naively.

If your three gluon does not contain any "S>gg" + one QCD radiation, then you do not have double counting and you can just remove the part of the code that tries to handle the double counting (so set ickkw to 0)

If you do have double counting, then I do not think that our MLM algorithm is going to work here.
Depending of what you want to do, you can
1) rely on the parton-shower to undo the running of alpha_s and to get the third gluon
2) going to NLO accuracy
3) try CKKW-L
4) implement yourself MLM in this case.

Cheers,

Olivier

Revision history for this message
Sarif Khan (sarifkhan) said : 		#2

Thanks a lot, Olivier, it solved the problem.

Revision history for this message
Sarif Khan (sarifkhan) said : 		#3

Thanks Olivier Mattelaer, that solved my question.