Incorrectly terminated jobs in w +jj

Asked by Matthew Lim

Hi Madgraph team,
I am trying to generate events for the process p p > w+ j j [QCD] with a no b mass model. Everything compiles fine, but after running the error gives
INFO: Return code of the event collection: None
INFO: Output of the event collection:
N of directories: 122
Determining the number of unweighted events per channel
50 jobs did not terminate correctly:
P0_dxg_wpuxg/GF1
P0_dxg_wpuxg/GF2
P0_dxg_wpuxg/GF3
P0_dxg_wpuxg/GF4
P0_dxg_wpuxg/GF5
P0_dxg_wpuxg/GF6
P0_dxg_wpuxg/GF7
P0_dxg_wpuxg/GF8
 et cetera and then this
The results are probably not correct. Please check the relevant log files corresponding to the above jobs for more information.
Found 72 correctly terminated jobs
72 122
Integrated abs(cross-section)
P0_uxu_wpuxd/GF4/res_0: Final result [ABS]: 164.20421054826977 +/- 4.5787611448815042
Integrated cross-section
P0_uxu_wpuxd/GF4/res_0: Final result: 39.930273101102422 +/- 4.4005328720454360

Error detected in "launch "
write debug file /home/malim/Applications/MG5_aMC_v2_2_2/BensingleW/run_02_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : An error occurred during the collection of results.

A typical log file for a subprocess looks like
 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.00000000000000008E-005
  > PrecisionVirtualAtRunTime = 1.00000000000000002E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.00000000000000010E-003
 ===============================================================
 Process in group number 0
 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
WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 MU_R = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.16639000000000003E-005
 aS = 0.11799999999999999
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 0.0000000000000000
 mdl_WH = 6.38233899999999989E-003
  Internal Params
  ---------------------------------

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM22 = 1.0000000000000000
 mdl_I4x33 = 0.0000000000000000
 mdl_I1x33 = 0.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM22 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 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__2 = 15625.000000000000
 mdl_Ncol__exp__2_m_1 = 8.0000000000000000
 mdl_MT__exp__2 = 29929.000000000000
 mdl_Ncol__exp__2_m_1_0 = 8.0000000000000000
 mdl_aEW = 7.54677111397888348E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.68721538467815552E-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_v__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yt = 0.99366614581500623
 mdl_ytau = 1.02066170006547170E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.54305075882732989E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.99366614581500623 , 0.0000000000000000 )
 mdl_Axial_tbGp = (-0.99366614581500623 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_Axial_tbGm = (-0.99366614581500623 , 0.0000000000000000 )
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_ee__exp__2 = 9.48355227599988748E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_yt__exp__2 = 0.98737240933884918
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__3 = 1.8056676068262196
 mdl_G__exp__4 = 2.1987899468922913
  Couplings of loop_sm-no_b_mass
  ---------------------------------

       UV_3Gt = 0.48815E-02 0.00000E+00
      UV_GQQt = 0.00000E+00 -0.48815E-02
   UVWfct_G_2 = 0.40088E-02 0.00000E+00
      R2_sxcW = -0.00000E+00 -0.11566E-01
         GC_4 = -0.12177E+01 0.00000E+00
         GC_5 = 0.00000E+00 0.12177E+01
         GC_6 = 0.00000E+00 0.14828E+01
       R2_3Gq = 0.76230E-02 0.00000E+00
       R2_3Gg = 0.31445E-01 0.00000E+00
       R2_GQQ = -0.00000E+00 -0.30492E-01
       R2_GGq = 0.00000E+00 0.62601E-02
       R2_GGt = -0.00000E+00 -0.11242E+04
     R2_GGg_1 = 0.00000E+00 0.28170E-01
     R2_GGg_2 = -0.00000E+00 -0.18780E-01
       R2_QQq = 0.00000E+00 0.12520E-01
  UV_3Gg_1eps = 0.62890E-01 0.00000E+00
  UV_3Gb_1eps = -0.38115E-02 0.00000E+00
 UV_GQQg_1eps = 0.00000E+00 -0.62890E-01
 UV_GQQq_1eps = 0.00000E+00 0.38115E-02
 UVWfct_G_2_1eps -0.31300E-02 0.00000E+00
        GC_47 = 0.00000E+00 0.46191E+00

 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
 Diagram information for clustering has been set-up for nFKSprocess 5
 Diagram information for clustering has been set-up for nFKSprocess 6
 Diagram information for clustering has been set-up for nFKSprocess 7
 Diagram information for clustering has been set-up for nFKSprocess 8
 Diagram information for clustering has been set-up for nFKSprocess 9
 Diagram information for clustering has been set-up for nFKSprocess 10
 Diagram information for clustering has been set-up for nFKSprocess 11
 Diagram information for clustering has been set-up for nFKSprocess 12
 Diagram information for clustering has been set-up for nFKSprocess 13
 Diagram information for clustering has been set-up for nFKSprocess 14
 Diagram information for clustering has been set-up for nFKSprocess 15
 Diagram information for clustering has been set-up for nFKSprocess 16
 Diagram information for clustering has been set-up for nFKSprocess 17
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.99999999999999989E-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)?
 Summing over 1 helicities/event for virt
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 10 -1 12 1
 imode is 0
 ------- iteration 1
 Update # PS points (even): 800 --> 800
Using random seed offsets: 1 , 2 , 0
  with seed 34
 Ranmar initialization seeds 13168 9410
 Total number of FKS directories is 17
 FKS process map (sum= 3 ) :
           1 --> 9 : 1 5 6 7 8 9 10 11 12
           2 --> 2 : 2 13
           3 --> 1 : 3
           4 --> 5 : 4 14 15 16 17
 ================================
 process combination map (specified per FKS dir):
  1 map 1 2
  1 inv. map 1 2
  2 map 1 2
  2 inv. map 1 2
  3 map 1 2
  3 inv. map 1 2
  4 map 1 2
  4 inv. map 1 2
  5 map 1 2
  5 inv. map 1 2
  6 map 1 2
  6 inv. map 1 2
  7 map 1 2
  7 inv. map 1 2
  8 map 1 2
  8 inv. map 1 2
  9 map 1 1 2 2 2
  9 inv. map 2 5
 10 map 1 1 2 2 2
 10 inv. map 2 5
 Number of unique IPROCs not identical among nFKSprocesses 11 1 2
Time in seconds: 0

This is not as far as I can tell a problem with the machine or the running time - the runs took about 20 minutes and still gave this result. What can I do to fix this?
Many thanks,
Matthew

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Rikkert Frederix Edit question
Solved by:
Matthew Lim
Solved:
Last query:
Last reply:
Revision history for this message
Rikkert Frederix (frederix) said :
#1

Dear Matthew,

What are the exact commands you did to generate the process? Or equivalently, can you copy the contents of the Cards/proc_card_mg5.dat here?

Best regards,
Rikkert

Revision history for this message
Matthew Lim (sengyanglim) said :
#2

Dear Rikkert,
Here are the contents of the proc card:
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_sm-no_b_mass
generate p p > w+ j j [QCD]
output BensingleW

Thanks,
Matthew

Revision history for this message
Rikkert Frederix (frederix) said :
#3

Dear Matthew,

When using a model with a massless b quark you should include the b in the definition of the p and j:

import model loop_sm-no_b_mass
define p = g u c d s u~ c~ d~ s~ b b~
define j = p
generate p p > w+ j j [QCD]
output BensingleW

That should solve your problem.

Best regards,
Rikkert

Revision history for this message
Matthew Lim (sengyanglim) said :
#4

Dear Rikkert,
Thank you - this solved my problem.
Best,
Matthew