QCD Loops: Error #6 fsk_Sij Infinity

Asked by Douglas Tuckler

Hi!

I'm having some trouble generating events for a 2HDM that I've been implementing in MadGraph by modifying the 2HDM model files already in MG. In particular, I am generating c s~ > h+ [QCD]. When I launch the even I get the following errors:

Error detected in "launch auto "
write debug file /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/run_02_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : An error occurred during the collection of results.
 Please check the .log files inside the directories which failed:
 /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/SubProcesses/P0_csx_hp/GF1/log.txt

Now, checking the log file, I see the following:

===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = -1.0000000000000000
  > PrecisionVirtualAtRunTime = -1.0000000000000000
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 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: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
 *****************************************************
 * 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
 mdl_tanbeta = 50.000000000000000
 mdl_Varx = 5.0000000000000003E-002
 mdl_CKM1x1 = 0.97198499999999999
 mdl_CKM1x2 = -0.23504200000000000
 mdl_CKM1x3 = 1.1680299999999999E-002
 mdl_CKM2x1 = 0.23482000000000000
 mdl_CKM2x2 = 0.97082700000000000
 mdl_CKM2x3 = -4.8518199999999997E-002
 mdl_CKM3x1 = -1.0269800000000001E-002
 mdl_CKM3x2 = 4.7433200000000002E-002
 mdl_CKM3x3 = 0.99882199999999999
 mdl_l2 = 0.50000000000000000
 mdl_l3 = 1.0000000000000000
 mdl_lR7 = 0.10000000000000001
 mdl_lI7 = 0.20000000000000001
 mdl_mixh = 0.29999999999999999
 mdl_mixh2 = 0.10000000000000001
 mdl_mixh3 = 0.20000000000000001
 MU_R = 91.188000000000002
 mdl_MD = 4.7999999999999996E-003
 mdl_MU = 2.3000000000000000E-003
 mdl_MS = 9.5000000000000001E-002
 mdl_MC = 1.2749999999999999
 mdl_MB = 4.1799999999999997
 mdl_MT = 176.69999999999999
 mdl_Me = 5.1099999999999995E-004
 mdl_MMU = 0.10565800000000000
 mdl_MTA = 1.7768200000000001
 mdl_MZ = 91.187600000000003
 mdl_mh1 = 125.02000000000000
 mdl_mh2 = 500.00000000000000
 mdl_mh3 = 500.00000000000000
 mdl_mhc = 500.00000000000000
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymdo = 1.8240600000000000E-003
 mdl_ymup = 1.1673800000000000E-003
 mdl_yms = 4.7142200000000002E-002
 mdl_ymc = 0.55251399999999995
 mdl_ymb = 2.5662900000000000
 mdl_ymt = 151.45400000000001
 mdl_yme = 5.0851700000000000E-004
 mdl_ymm = 9.9399500000000002E-002
 mdl_ymtau = 1.8887400000000001
 mdl_GDI1x1 = 0.0000000000000000
 mdl_GDI1x2 = 0.0000000000000000
 mdl_GDI1x3 = 0.0000000000000000
 mdl_GDI2x1 = 0.0000000000000000
 mdl_GDI2x2 = 0.0000000000000000
 mdl_GDI2x3 = 0.0000000000000000
 mdl_GDI3x1 = 0.0000000000000000
 mdl_GDI3x2 = 0.0000000000000000
 mdl_GDI3x3 = 0.0000000000000000
 mdl_GLI1x1 = 0.0000000000000000
 mdl_GLI1x2 = 0.0000000000000000
 mdl_GLI1x3 = 0.0000000000000000
 mdl_GLI2x1 = 0.0000000000000000
 mdl_GLI2x2 = 0.0000000000000000
 mdl_GLI2x3 = 0.0000000000000000
 mdl_GLI3x1 = 0.0000000000000000
 mdl_GLI3x2 = 0.0000000000000000
 mdl_GLI3x3 = 0.0000000000000000
 mdl_GUI1x1 = 0.0000000000000000
 mdl_GUI1x2 = 0.0000000000000000
 mdl_GUI1x3 = 0.0000000000000000
 mdl_GUI2x1 = 0.0000000000000000
 mdl_GUI2x2 = 0.0000000000000000
 mdl_GUI2x3 = 0.0000000000000000
 mdl_GUI3x1 = 0.0000000000000000
 mdl_GUI3x2 = 0.0000000000000000
 mdl_GUI3x3 = 0.0000000000000000
 mdl_WT = 1.5083359999999999
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_Wh1 = 1.0000000000000000
 mdl_Wh2 = 1.0000000000000000
 mdl_Wh3 = 1.0000000000000000
 mdl_whc = 1.0000000000000000
  Internal Params
  ---------------------------------

 mdl_atan__tanbeta = 1.5507989928217460
 mdl_thb = 1.5507989928217460
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_l7 = ( 0.10000000000000001 , 0.20000000000000001 )
 mdl_yup = 0.0000000000000000
 mdl_yc = 0.0000000000000000
 mdl_ydo = 0.0000000000000000
 mdl_ys = 0.0000000000000000
 mdl_ye = 0.0000000000000000
 mdl_ym = 0.0000000000000000
 mdl_cos__mixh = 0.95533648912560598
 mdl_cos__mixh2 = 0.99500416527802571
 mdl_TH1x1 = 0.95056378592206325
 mdl_sin__mixh = 0.29552020666133955
 mdl_TH1x2 = 0.29404383655185579
 mdl_sin__mixh2 = 9.9833416646828155E-002
 mdl_TH1x3 = -9.9833416646828155E-002
 mdl_cos__mixh3 = 0.98006657784124163
 mdl_sin__mixh3 = 0.19866933079506122
 mdl_TH2x1 = -0.27068148839190342
 mdl_TH2x2 = 0.94215466351136834
 mdl_TH2x3 = 0.19767681165408385
 mdl_TH3x1 = 0.15218416716418803
 mdl_TH3x2 = -0.16088136066569614
 mdl_TH3x3 = 0.97517032720181585
 mdl_alpha = 3.0002666026849492E-002
 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_mhc__exp__2 = 250000.00000000000
 mdl_mh1__exp__2 = 15630.000399999999
 mdl_mh2__exp__2 = 250000.00000000000
 mdl_mh3__exp__2 = 250000.00000000000
 mdl_conjg__CKM1x1 = ( 0.97198499999999999 , -0.0000000000000000 )
 mdl_conjg__CKM1x2 = (-0.23504200000000000 , -0.0000000000000000 )
 mdl_conjg__CKM1x3 = ( 1.1680299999999999E-002, -0.0000000000000000 )
 mdl_conjg__CKM2x1 = ( 0.23482000000000000 , -0.0000000000000000 )
 mdl_conjg__CKM2x2 = ( 0.97082700000000000 , -0.0000000000000000 )
 mdl_conjg__CKM2x3 = ( -4.8518199999999997E-002, -0.0000000000000000 )
 mdl_conjg__CKM3x1 = ( -1.0269800000000001E-002, -0.0000000000000000 )
 mdl_conjg__CKM3x2 = ( 4.7433200000000002E-002, -0.0000000000000000 )
 mdl_conjg__CKM3x3 = ( 0.99882199999999999 , -0.0000000000000000 )
 mdl_I7a11 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I7a12 = ( -0.0000000000000000 , 0.0000000000000000 )
 mdl_I7a13 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I7a21 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I7a22 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I7a23 = ( -0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a11 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a12 = ( -0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a21 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a22 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a31 = ( -0.0000000000000000 , 0.0000000000000000 )
 mdl_I8a32 = ( 0.0000000000000000 , 0.0000000000000000 )
 mdl_TH1x1__exp__2 = 0.90357151110648615
 mdl_TH2x1__exp__2 = 7.3268468158056155E-002
 mdl_TH3x1__exp__2 = 2.3160020735457528E-002
 mdl_TH1x2__exp__2 = 8.6461777814134486E-002
 mdl_TH2x2__exp__2 = 0.88765540997621972
 mdl_TH3x2__exp__2 = 2.5882812209645799E-002
 mdl_TH1x3__exp__2 = 9.9667110793791851E-003
 mdl_TH2x3__exp__2 = 3.9076121865724144E-002
 mdl_TH3x3__exp__2 = 0.95095716705489652
 mdl_cos__alpha = 0.99954995377662637
 mdl_cos__thb = 1.9996001199600179E-002
 mdl_sin__alpha = 2.9998165029616445E-002
 mdl_tan__thb = 49.999999999999901
 mdl_conjg__l7 = ( 0.10000000000000001 ,-0.20000000000000001 )
 mdl_TH2x1__exp__3 = -1.9832418013217421E-002
 mdl_TH3x1__exp__3 = 3.5245884671309294E-003
 mdl_TH1x1__exp__4 = 0.81644147568325887
 mdl_TH1x1__exp__3 = 0.85890235644870105
 mdl_TH2x1__exp__4 = 5.3682684262280882E-003
 mdl_TH3x1__exp__4 = 5.3638656046682272E-004
 mdl_TH2x2__exp__3 = 0.83630868410019099
 mdl_TH3x2__exp__3 = -4.1640620461425096E-003
 mdl_TH1x2__exp__3 = 2.5423552863562231E-002
 mdl_TH1x2__exp__4 = 7.4756390227807583E-003
 mdl_TH2x2__exp__4 = 0.78793212686005076
 mdl_TH3x2__exp__4 = 6.6991996787978960E-004
 mdl_TH2x3__exp__3 = 7.7244431822227795E-003
 mdl_TH3x3__exp__3 = 0.92734521175183526
 mdl_TH1x3__exp__3 = -9.9501081978622063E-004
 mdl_TH1x3__exp__4 = 9.9335329739819796E-005
 mdl_TH2x3__exp__4 = 1.5269433000649245E-003
 mdl_TH3x3__exp__4 = 0.90431953357307437
 mdl_MB__exp__2 = 17.472399999999997
 mdl_MC__exp__2 = 1.6256249999999999
 mdl_MD__exp__2 = 2.3039999999999996E-005
 mdl_MS__exp__2 = 9.0250000000000000E-003
 mdl_MT__exp__2 = 31222.889999999996
 mdl_MU__exp__2 = 5.2900000000000002E-006
 mdl_tanbeta__exp__2 = 2500.0000000000000
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824660036055974
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9763498719121
 mdl_sw2 = 0.23369336794341478
 mdl_cw = 0.87538941737753784
 mdl_sqrt__sw2 = 0.48341841911889827
 mdl_sw = 0.48341841911889827
 mdl_g1 = 0.35807035569216145
 mdl_gw = 0.64840516548963911
 mdl_vev = 246.21845810181625
 mdl_vev__exp__2 = 60623.529110035844
 mdl_mu2 = 219688.23544498207
 mdl_GDR1x1 = 1.4971115143756042E-008
 mdl_GDR1x2 = 1.5933719228415516E-004
 mdl_GDR1x3 = 8.0407679564786592E-003
 mdl_GDR2x1 = 6.8584643214418930E-006
 mdl_GDR2x2 = 7.2994435811177131E-004
 mdl_GDR2x3 = 3.6829482231482816E-002
 mdl_GDR3x1 = 1.3429483370028448E-004
 mdl_GDR3x2 = 1.4292968395082640E-002
 mdl_GDR3x3 = 0.72115292847721524
 mdl_GLR1x1 = 9.8859913542204794E-013
 mdl_GLR1x2 = 4.1734582824645412E-008
 mdl_GLR1x3 = 7.0915233972442430E-008
 mdl_GLR2x1 = 4.2050614608894160E-008
 mdl_GLR2x2 = 1.7861829144117122E-003
 mdl_GLR2x3 = 3.0143686184889551E-002
 mdl_GLR3x1 = 7.0915233972442430E-008
 mdl_GLR3x2 = 3.0143686184889551E-002
 mdl_GLR3x3 = 0.50869625199166457
 mdl_GUR1x1 = 1.1595866957817220E-014
 mdl_GUR1x2 = 2.6056044896929652E-009
 mdl_GUR1x3 = 7.0901156192962249E-007
 mdl_GUR2x1 = 2.6056562039849330E-009
 mdl_GUR2x2 = 5.8549772159655089E-004
 mdl_GUR2x3 = 0.15931909228196639
 mdl_GUR3x1 = 7.0902592701072475E-007
 mdl_GUR3x2 = 0.15931909228196639
 mdl_GUR3x3 = 43.352378258376959
 mdl_l1 = 0.31530616795676347
 mdl_l4 = -0.37279164105377238
 mdl_lI5 = -0.15925344377464387
 mdl_lI6 = 0.55925714512245139
 mdl_lR5 = -9.6859395805504928E-002
 mdl_lR6 = 0.99471948532090637
 mdl_yb = 1.4743037280067919E-002
 mdl_yt = 0.87008559758071258
 mdl_ytau = 1.0850591411085842E-002
 mdl_mu1 = -19114.972651700700
 mdl_GD1x1 = ( 1.4971115143756042E-008, 0.0000000000000000 )
 mdl_GD1x2 = ( 1.5933719228415516E-004, 0.0000000000000000 )
 mdl_GD1x3 = ( 8.0407679564786592E-003, 0.0000000000000000 )
 mdl_GD2x1 = ( 6.8584643214418930E-006, 0.0000000000000000 )
 mdl_GD2x2 = ( 7.2994435811177131E-004, 0.0000000000000000 )
 mdl_GD2x3 = ( 3.6829482231482816E-002, 0.0000000000000000 )
 mdl_GD3x1 = ( 1.3429483370028448E-004, 0.0000000000000000 )
 mdl_GD3x2 = ( 1.4292968395082640E-002, 0.0000000000000000 )
 mdl_GD3x3 = ( 0.72115292847721524 , 0.0000000000000000 )
 mdl_GL1x1 = ( 9.8859913542204794E-013, 0.0000000000000000 )
 mdl_GL1x2 = ( 4.1734582824645412E-008, 0.0000000000000000 )
 mdl_GL1x3 = ( 7.0915233972442430E-008, 0.0000000000000000 )
 mdl_GL2x1 = ( 4.2050614608894160E-008, 0.0000000000000000 )
 mdl_GL2x2 = ( 1.7861829144117122E-003, 0.0000000000000000 )
 mdl_GL2x3 = ( 3.0143686184889551E-002, 0.0000000000000000 )
 mdl_GL3x1 = ( 7.0915233972442430E-008, 0.0000000000000000 )
 mdl_GL3x2 = ( 3.0143686184889551E-002, 0.0000000000000000 )
 mdl_GL3x3 = ( 0.50869625199166457 , 0.0000000000000000 )
 mdl_GU1x1 = ( 1.1595866957817220E-014, 0.0000000000000000 )
 mdl_GU1x2 = ( 2.6056044896929652E-009, 0.0000000000000000 )
 mdl_GU1x3 = ( 7.0901156192962249E-007, 0.0000000000000000 )
 mdl_GU2x1 = ( 2.6056562039849330E-009, 0.0000000000000000 )
 mdl_GU2x2 = ( 5.8549772159655089E-004, 0.0000000000000000 )
 mdl_GU2x3 = ( 0.15931909228196639 , 0.0000000000000000 )
 mdl_GU3x1 = ( 7.0902592701072475E-007, 0.0000000000000000 )
 mdl_GU3x2 = ( 0.15931909228196639 , 0.0000000000000000 )
 mdl_GU3x3 = ( 43.352378258376959 , 0.0000000000000000 )
 mdl_l5 = ( -9.6859395805504928E-002,-0.15925344377464387 )
 mdl_l6 = ( 0.99471948532090637 , 0.55925714512245139 )
 mdl_mu3 = ( -30151.702837335921 , -16952.070908663238 )
 mdl_conjg__GU1x1 = ( 1.1595866957817220E-014, -0.0000000000000000 )
 mdl_conjg__GU2x1 = ( 2.6056562039849330E-009, -0.0000000000000000 )
 mdl_conjg__GU3x1 = ( 7.0902592701072475E-007, -0.0000000000000000 )
 mdl_I1a11 = ( -6.6696830043862544E-009, 0.0000000000000000 )
 mdl_I1a12 = ( 3.6161007271115428E-008, 0.0000000000000000 )
 mdl_I1a13 = ( 7.0806427285531305E-007, 0.0000000000000000 )
 mdl_conjg__GU1x2 = ( 2.6056044896929652E-009, -0.0000000000000000 )
 mdl_conjg__GU2x2 = ( 5.8549772159655089E-004, -0.0000000000000000 )
 mdl_conjg__GU3x2 = ( 0.15931909228196639 , -0.0000000000000000 )
 mdl_I1a21 = ( -1.4986861063235567E-003, 0.0000000000000000 )
 mdl_I1a22 = ( 8.1254307521668925E-003, 0.0000000000000000 )
 mdl_I1a23 = ( 0.15910300712613651 , 0.0000000000000000 )
 mdl_conjg__GU1x3 = ( 7.0901156192962249E-007, -0.0000000000000000 )
 mdl_conjg__GU2x3 = ( 0.15931909228196639 , -0.0000000000000000 )
 mdl_conjg__GU3x3 = ( 43.352378258376959 , -0.0000000000000000 )
 mdl_I1a31 = (-0.40780825583962538 , 0.0000000000000000 )
 mdl_I1a32 = ( 2.2110131381605749 , 0.0000000000000000 )
 mdl_I1a33 = ( 43.293579289486907 , 0.0000000000000000 )
 mdl_I2a11 = ( -2.8871525617908834E-008, 0.0000000000000000 )
 mdl_I2a12 = ( 1.5025193776809135E-004, 0.0000000000000000 )
 mdl_I2a13 = ( 7.5823132300181419E-003, 0.0000000000000000 )
 mdl_I2a21 = ( 1.4615425861338342E-007, 0.0000000000000000 )
 mdl_I2a22 = ( 5.2596151658443422E-005, 0.0000000000000000 )
 mdl_I2a23 = ( 2.6541468634408635E-003, 0.0000000000000000 )
 mdl_I2a31 = ( 1.3446179954567907E-004, 0.0000000000000000 )
 mdl_I2a32 = ( 1.4309118513943101E-002, 0.0000000000000000 )
 mdl_I2a33 = ( 0.72196777344529195 , 0.0000000000000000 )
 mdl_conjg__GD1x1 = ( 1.4971115143756042E-008, -0.0000000000000000 )
 mdl_conjg__GD2x1 = ( 6.8584643214418930E-006, -0.0000000000000000 )
 mdl_conjg__GD3x1 = ( 1.3429483370028448E-004, -0.0000000000000000 )
 mdl_I3a11 = ( -2.8871525617908834E-008, 0.0000000000000000 )
 mdl_I3a12 = ( 1.4615425861338342E-007, 0.0000000000000000 )
 mdl_I3a13 = ( 1.3446179954567907E-004, 0.0000000000000000 )
 mdl_conjg__GD1x2 = ( 1.5933719228415516E-004, -0.0000000000000000 )
 mdl_conjg__GD2x2 = ( 7.2994435811177131E-004, -0.0000000000000000 )
 mdl_conjg__GD3x2 = ( 1.4292968395082640E-002, -0.0000000000000000 )
 mdl_I3a21 = ( 1.5025193776809135E-004, 0.0000000000000000 )
 mdl_I3a22 = ( 5.2596151658443422E-005, 0.0000000000000000 )
 mdl_I3a23 = ( 1.4309118513943101E-002, 0.0000000000000000 )
 mdl_conjg__GD1x3 = ( 8.0407679564786592E-003, -0.0000000000000000 )
 mdl_conjg__GD2x3 = ( 3.6829482231482816E-002, -0.0000000000000000 )
 mdl_conjg__GD3x3 = ( 0.72115292847721524 , -0.0000000000000000 )
 mdl_I3a31 = ( 7.5823132300181419E-003, 0.0000000000000000 )
 mdl_I3a32 = ( 2.6541468634408635E-003, 0.0000000000000000 )
 mdl_I3a33 = ( 0.72196777344529195 , 0.0000000000000000 )
 mdl_I4a11 = ( -6.6696830043862544E-009, 0.0000000000000000 )
 mdl_I4a12 = ( -1.4986861063235567E-003, 0.0000000000000000 )
 mdl_I4a13 = (-0.40780825583962538 , 0.0000000000000000 )
 mdl_I4a21 = ( 3.6161007271115428E-008, 0.0000000000000000 )
 mdl_I4a22 = ( 8.1254307521668925E-003, 0.0000000000000000 )
 mdl_I4a23 = ( 2.2110131381605749 , 0.0000000000000000 )
 mdl_I4a31 = ( 7.0806427285531305E-007, 0.0000000000000000 )
 mdl_I4a32 = ( 0.15910300712613651 , 0.0000000000000000 )
 mdl_I4a33 = ( 43.293579289486907 , 0.0000000000000000 )
 mdl_I5a11 = ( -5.0917031664077543E-004, 0.0000000000000000 )
 mdl_I5a21 = ( -1.2300948446075495E-004, 0.0000000000000000 )
 mdl_I5a31 = ( 5.3797921962143821E-006, 0.0000000000000000 )
 mdl_I5a12 = ( 3.1821431657991263E-003, 0.0000000000000000 )
 mdl_I5a22 = ( -1.3143653062955847E-002, 0.0000000000000000 )
 mdl_I5a32 = ( -6.4217983684610889E-004, 0.0000000000000000 )
 mdl_I5a13 = ( -8.6084334025725920E-003, 0.0000000000000000 )
 mdl_I5a23 = ( 3.5758130656977784E-002, 0.0000000000000000 )
 mdl_I5a33 = (-0.73613628657006780 , 0.0000000000000000 )
 mdl_I6a11 = ( -3.2586386645182085E-004, 0.0000000000000000 )
 mdl_I6a21 = ( -3.7259992707120132E-002, 0.0000000000000000 )
 mdl_I6a31 = ( 0.44669092424890250 , 0.0000000000000000 )
 mdl_I6a12 = ( 7.8799256057005909E-005, 0.0000000000000000 )
 mdl_I6a22 = (-0.15404568154277881 , 0.0000000000000000 )
 mdl_I6a32 = ( -2.0631346226881773 , 0.0000000000000000 )
 mdl_I6a13 = ( -3.9158914173749626E-006, 0.0000000000000000 )
 mdl_I6a23 = ( 7.6986107578681370E-003, 0.0000000000000000 )
 mdl_I6a33 = ( -43.444343837283810 , 0.0000000000000000 )
 mdl_I7a31 = ( -8.9356050700344029E-003, 0.0000000000000000 )
 mdl_I7a32 = ( 4.1270944167165459E-002, 0.0000000000000000 )
 mdl_I7a33 = ( 0.86906063674676248 , 0.0000000000000000 )
 mdl_I8a13 = ( 1.7220309834237730E-004, 0.0000000000000000 )
 mdl_I8a23 = ( -7.1530563136179125E-004, 0.0000000000000000 )
 mdl_I8a33 = ( 1.4725669982151999E-002, 0.0000000000000000 )
 mdl_gw__exp__2 = 0.42042925863364627
 mdl_sw__exp__2 = 0.23369336794341478
 mdl_cw__exp__2 = 0.76630663205658511
 mdl_g1__exp__2 = 0.12821437962551102
 mdl_conjg__GL1x1 = ( 9.8859913542204794E-013, -0.0000000000000000 )
 mdl_conjg__GL1x2 = ( 4.1734582824645412E-008, -0.0000000000000000 )
 mdl_conjg__GL1x3 = ( 7.0915233972442430E-008, -0.0000000000000000 )
 mdl_conjg__GL2x1 = ( 4.2050614608894160E-008, -0.0000000000000000 )
 mdl_conjg__GL2x2 = ( 1.7861829144117122E-003, -0.0000000000000000 )
 mdl_conjg__GL2x3 = ( 3.0143686184889551E-002, -0.0000000000000000 )
 mdl_conjg__GL3x1 = ( 7.0915233972442430E-008, -0.0000000000000000 )
 mdl_conjg__GL3x2 = ( 3.0143686184889551E-002, -0.0000000000000000 )
 mdl_conjg__GL3x3 = ( 0.50869625199166457 , -0.0000000000000000 )
 mdl_conjg__l5 = ( -9.6859395805504928E-002, 0.15925344377464387 )
 mdl_conjg__l6 = ( 0.99471948532090637 ,-0.55925714512245139 )
 mdl_yb__exp__2 = 2.1735714824147246E-004
 mdl_yt__exp__2 = 0.75704894711738568
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = ( 1.4828317324943823 , 0.0000000000000000 )
 mdl_G__exp__3 = ( 1.8056676068262196 , 0.0000000000000000 )
 mdl_G__exp__4 = ( 2.1987899468922913 , 0.0000000000000000 )
  Couplings of F2HDM_NLO_NEW
  ---------------------------------

 UVGC_702_188 = 0.00000E+00 0.10540E-03
 UVGC_702_189 = 0.00000E+00 0.16084E-03
 UVGC_713_215 = 0.00000E+00 -0.13834E-03
 UVGC_713_216 = 0.00000E+00 -0.21111E-03
 R2GC_702_222 = -0.00000E+00 -0.28467E-04
 R2GC_713_233 = 0.00000E+00 0.37364E-04
 UVGC_702_188_1e 0.00000E+00 0.00000E+00
 UVGC_702_189_1e 0.00000E+00 0.00000E+00
 UVGC_702_190_1e -0.00000E+00 -0.28467E-04
 UVGC_713_215_1e 0.00000E+00 0.00000E+00
 UVGC_713_216_1e 0.00000E+00 0.00000E+00
 UVGC_713_217_1e 0.00000E+00 0.37364E-04
         GC_2 = 0.00000E+00 0.12177E+01
        GC_43 = 0.00000E+00 -0.14592E+00
        GC_52 = 0.00000E+00 -0.13091E-01

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

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

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 1
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 Not subdividing B.W.
 about to integrate 4 -1 12 1
 imode is 0
Note: The following floating-point exceptions are signalling: IEEE_DIVIDE_BY_ZERO
 ------- iteration 1
 Update # PS points (even): 320 --> 243
Using random seed offsets: 1 , 1 , 0
  with seed 34
 Ranmar initialization seeds 13168 9409
 Total number of FKS directories is 4
 FKS process map (sum= 3 ) :
           1 --> 2 : 1 3
           2 --> 2 : 2 4
 ================================
 process combination map (specified per FKS dir):
  1 map 1
  1 inv. map 1
  2 map 1
  2 inv. map 1
  3 map 1
  3 inv. map 1
  4 map 1
  4 inv. map 1
 ================================
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 1. Lower bound for tau is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 2. Lower bound for tau is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 3. Lower bound for tau is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 4. Lower bound for tau is 0.14793E-02 0.50000E+03 0.13000E+05
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.14793E-02 0.50000E+03 0.13000E+05
 Mass shell violation [nocms]
 j= 1
 mass= 1.2749999999999999
 mass computed= 0.0000000000000000
0.25000000D+03 0.00000000D+00 0.00000000D+00 0.25000000D+03
 bpower is 0.0000000000000000
 Mass shell violation [nocms]
 j= 1
 mass= 1.2749999999999999
 mass computed= 0.0000000000000000
0.25000000D+03 0.00000000D+00 0.00000000D+00 0.25000000D+03
 Scale values (may change event by event):
 muR, muR_reference: 0.250672D+03 0.250672D+03 1.00
 muF1, muF1_reference: 0.250672D+03 0.250672D+03 1.00
 muF2, muF2_reference: 0.250672D+03 0.250672D+03 1.00
 QES, QES_reference: 0.250672D+03 0.250672D+03 1.00

 muR_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF1_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF2_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 QES_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state

 alpha_s= 0.10317716300178234
 Mass shell violation [nocms]
 j= 1
 mass= 1.2749999999999999
 mass computed= 0.0000000000000000
0.25000000D+03 0.00000000D+00 0.00000000D+00 0.25000000D+03
 Error #6 in fks_Sij Infinity
Time in seconds: 0

So, in particular the error "Error #y in fks_Sij Infinity" seems to be the problem. I'm not sure why this error is generated, and I'm hoping you could provide me with some suggestions to fix this error.

Thanks!

-Douglas.

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said :
#1

Dear Douglas,
thanks for reporting this error.
Are you using a model where the charm and strange quarks are massless?
If not, can you try to generate a restriction card for the model with MC=0 and regenerate the process folder?
Let us know,

Cheers,

Marco

Revision history for this message
Douglas Tuckler (dtuck001) said :
#2

Hi Marco,

Thanks for your reply. I am using a model where the charm and strange quarks are not massless, but I have generated a run where they are massless. In particular I run the following:

generate c s~ > h+ [QCD]
output Runs/F2HDM_NLO/cs_Hpl_QCD
launch Runs/F2HDM_NLO/cs_Hpl_QCD

Then I use the following settings:

The following switches determine which operations are executed:
 1 Perturbative order of the calculation: order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=OFF
 3 Shower the generated events: shower=ON
 4 Decay particles with the MadSpin module: madspin=OFF
 5 Add weights to the events based on changing model parameters: reweight=OFF
 6 Run MadAnalysis5 on the events generated: madanalysis5=OFF

Finally, I set the mass of the charm quark to zero using the ‘set’ command. After running, I get the following:

INFO: Result for test_ME:
INFO: Output of the failing test:
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 Sum over helicity (0), or random helicity (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: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 4 j= 1
 with PDGs: i= 21 j= 4

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 0 , 1 , 0
  with seed 32
 Ranmar initialization seeds 10729 9407
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 1. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 2. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 3. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 4. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04

 Soft tests done for (Born) config 1
 Failures: 100
     Soft test 1 FAILED. Fraction of failures: 1.00

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 1 PASSED. Fraction of failures: 0.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 2
 FKS partons are: i= 4 j= 2
 with PDGs: i= 21 j= -3

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 100
     Soft test 2 FAILED. Fraction of failures: 1.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 3
 FKS partons are: i= 4 j= 1
 with PDGs: i= -4 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 3 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 3 PASSED. Fraction of failures: 0.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 4 j= 2
 with PDGs: i= 3 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 31
     Soft test 4 FAILED. Fraction of failures: 0.31

 No collinear test for massive j_fks
Error detected in "launch auto "
write debug file /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Some tests failed, run cannot continue.
Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

If I choose not to shower the final states, I receive the same error.

If I set both the charm and quark mass equal to zero, I get the following error:

INFO: Result for check_poles:
Error detected in "launch auto "
write debug file /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/run_03_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Poles do not cancel, run cannot continue

-Douglas.

On Mar 27, 2017, at 3:27 AM, marco zaro <<email address hidden><mailto:<email address hidden>>> wrote:

Your question #591230 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/591230

   Status: Open => Needs information

marco zaro requested more information:
Dear Douglas,
thanks for reporting this error.
Are you using a model where the charm and strange quarks are massless?
If not, can you try to generate a restriction card for the model with MC=0 and regenerate the process folder?
Let us know,

Cheers,

Marco

--
To answer this request for more information, you can either reply to
this email or enter your reply at the following page:
https://answers.launchpad.net/mg5amcnlo/+question/591230

You received this question notification because you asked the question.

Revision history for this message
marco zaro (marco-zaro) said :
#3

Hi Douglas,
be aware that generating a process with a model with massive particles and then setting the mass of these particles to zero is extremely dangerous and may lead to the code crashing horribly. This is because if a particle is massive it will not have collinear singularities, while if it is massless it does have. The structure of the code is generated depending on the model.
So I would suggest to create a restriction of the model with massless charm and strange quarks (this can be done while keeping their yukawas different from zero) and use this restriction to generate the new process.
This way the code should not crash, hopefully.

Please let me know if you need help on doing this..

Cheers,

Marco
On 27 Mar 2017, at 16:43, Douglas Tuckler <email address hidden> wrote:

> Question #591230 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/591230
>
> Status: Needs information => Open
>
> Douglas Tuckler gave more information on the question:
> Hi Marco,
>
> Thanks for your reply. I am using a model where the charm and strange
> quarks are not massless, but I have generated a run where they are
> massless. In particular I run the following:
>
> generate c s~ > h+ [QCD]
> output Runs/F2HDM_NLO/cs_Hpl_QCD
> launch Runs/F2HDM_NLO/cs_Hpl_QCD
>
> Then I use the following settings:
>
> The following switches determine which operations are executed:
> 1 Perturbative order of the calculation: order=NLO
> 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=OFF
> 3 Shower the generated events: shower=ON
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weights to the events based on changing model parameters: reweight=OFF
> 6 Run MadAnalysis5 on the events generated: madanalysis5=OFF
>
>
> Finally, I set the mass of the charm quark to zero using the ‘set’ command. After running, I get the following:
>
> INFO: Result for test_ME:
> INFO: Output of the failing test:
> Enter xi_i, y_ij to be used in coll/soft tests
> Enter -2 to generate them randomly
> Enter number of tests for soft and collinear limits
> Sum over helicity (0), or random helicity (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: NNPDF23nlo_as_0119_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23nlo: 0.11899999999999999
> Give FKS configuration number ("0" loops over all)
>
> =================================================
>
> NEW FKS CONFIGURATION:
> FKS configuration number is 1
> FKS partons are: i= 4 j= 1
> with PDGs: i= 21 j= 4
>
>
> Enter graph number (iconfig), '0' loops over all graphs
> Using random seed offsets: 0 , 1 , 0
> with seed 32
> Ranmar initialization seeds 10729 9407
> nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 1. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 2. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 3. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 4. Lower bound for tau is 0.25000E+00 0.50000E+03 0.10000E+04
> nFKSprocess: 4. Lower bound for tau is (taking resonances into account) 0.25000E+00 0.50000E+03 0.10000E+04
>
>
> Soft tests done for (Born) config 1
> Failures: 100
> Soft test 1 FAILED. Fraction of failures: 1.00
>
>
> Collinear tests done for (Born) config 1
> Failures: 0
> Collinear test 1 PASSED. Fraction of failures: 0.00
>
> =================================================
>
> NEW FKS CONFIGURATION:
> FKS configuration number is 2
> FKS partons are: i= 4 j= 2
> with PDGs: i= 21 j= -3
>
>
> Enter graph number (iconfig), '0' loops over all graphs
>
>
> Soft tests done for (Born) config 1
> Failures: 100
> Soft test 2 FAILED. Fraction of failures: 1.00
>
>
> No collinear test for massive j_fks
>
> =================================================
>
> NEW FKS CONFIGURATION:
> FKS configuration number is 3
> FKS partons are: i= 4 j= 1
> with PDGs: i= -4 j= 21
>
>
> Enter graph number (iconfig), '0' loops over all graphs
>
>
> Soft tests done for (Born) config 1
> Failures: 0
> Soft test 3 PASSED. Fraction of failures: 0.00
>
>
> Collinear tests done for (Born) config 1
> Failures: 0
> Collinear test 3 PASSED. Fraction of failures: 0.00
>
> =================================================
>
> NEW FKS CONFIGURATION:
> FKS configuration number is 4
> FKS partons are: i= 4 j= 2
> with PDGs: i= 3 j= 21
>
>
> Enter graph number (iconfig), '0' loops over all graphs
>
>
> Soft tests done for (Born) config 1
> Failures: 31
> Soft test 4 FAILED. Fraction of failures: 0.31
>
>
> No collinear test for massive j_fks
> Error detected in "launch auto "
> write debug file /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/run_01_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
> aMCatNLOError : Some tests failed, run cannot continue.
> Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.
>
> If I choose not to shower the final states, I receive the same error.
>
> If I set both the charm and quark mass equal to zero, I get the
> following error:
>
> INFO: Result for check_poles:
> Error detected in "launch auto "
> write debug file /Users/dtuckler/Desktop/Douglas/Physics/MG5_aMC_v2_5_2/Runs/F2HDM_NLO/cs_Hpl_QCD/run_03_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
> aMCatNLOError : Poles do not cancel, run cannot continue
>
> -Douglas.
>
>
> On Mar 27, 2017, at 3:27 AM, marco zaro
> <<email address hidden><mailto:<email address hidden>>>
> wrote:
>
> Your question #591230 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/591230
>
> Status: Open => Needs information
>
> marco zaro requested more information:
> Dear Douglas,
> thanks for reporting this error.
> Are you using a model where the charm and strange quarks are massless?
> If not, can you try to generate a restriction card for the model with MC=0 and regenerate the process folder?
> Let us know,
>
> Cheers,
>
> Marco
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https://answers.launchpad.net/mg5amcnlo/+question/591230
>
> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#4

Olivier Mattelaer suggests this article as an answer to your question:
FAQ #2312: “FR Model much slower than build-in MG model. Why and how to fix?”.

Can you help with this problem?

Provide an answer of your own, or ask Douglas Tuckler for more information if necessary.

To post a message you must log in.