Dear experts, I am using MG5 3.3.2 with the SMEFTatNLO model to generate the tWa process at NLO in QCD. The restrict_NLO.dat card that I am using is pasted at the bottom. The generation of the pp->tWa process does not converge, as during generation there are jobs stuck indefinitely. To factor out that this could be due to the Diagram Removal I am applying, I have tested the same with generating the tqa process, but also in this case the generation does not converge. A noticeable difference is that everything works well in case I try to generate the tta process. I paste here also the run_banners for tWa and tqa. For tqa, I have tried setting nlo_mixed_expansion at both True and False and both giving the coupling orders linear or squared in the proc_card.dat. I am not sure what it should be of the acknowledged_v3.1_syntax. Is there something I am not considering with the generation of these processes? Thank you for your help. Best, Michel restrict_NLO.dat ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO ##################### ###################################################################### ################################### ## INFORMATION FOR SMINPUTS ################################### Block SMINPUTS 2 1.166370e-05 # Gf 3 1.184000e-01 # aS ################################### ## INFORMATION FOR MASS ################################### Block MASS 6 1.720000e+02 # MT 23 9.118760e+01 # MZ 24 7.982440e+01 # MW 25 1.250000e+02 # MH ## Not dependent paramater. ## Those values should be edited following analytical the ## analytical expression. Some generator could simply ignore ## those values and use the analytical expression 22 0.000000 # a : 0.0 21 0.000000 # g : 0.0 9000001 0.000000 # ghA : 0.0 9000005 0.000000 # ghG : 0.0 12 0.000000 # ve : 0.0 14 0.000000 # vm : 0.0 16 0.000000 # vt : 0.0 11 0.000000 # e- : 0.0 13 0.000000 # mu- : 0.0 15 0.000000 # ta- : 0.0 2 0.000000 # u : 0.0 4 0.000000 # c : 0.0 1 0.000000 # d : 0.0 3 0.000000 # s : 0.0 5 0.000000 # b : 0.0 9000002 91.187600 # ghZ : MZ 9000003 79.824400 # ghWp : MW 9000004 79.824400 # ghWm : MW 250 91.187600 # G0 : MZ 251 79.824400 # G+ : MW ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.508336e+00 DECAY 23 2.495200e+00 DECAY 24 0. DECAY 25 0. ## Not dependent paramater. ## Those values should be edited following analytical the ## analytical expression. Some generator could simply ignore ## those values and use the analytical expression DECAY 22 0. # a : 0.0 DECAY 21 0. # g : 0.0 DECAY 9000001 0. # ghA : 0.0 DECAY 9000005 0. # ghG : 0.0 DECAY 12 0. # ve : 0.0 DECAY 14 0. # vm : 0.0 DECAY 16 0. # vt : 0.0 DECAY 11 0. # e- : 0.0 DECAY 13 0. # mu- : 0.0 DECAY 15 0. # ta- : 0.0 DECAY 2 0. # u : 0.0 DECAY 4 0. # c : 0.0 DECAY 1 0. # d : 0.0 DECAY 3 0. # s : 0.0 DECAY 5 0. # b : 0.0 DECAY 9000002 0. # ghZ : WZ DECAY 9000003 0. # ghWp : WW DECAY 9000004 0. # ghWm : WW DECAY 250 0. # G0 : WZ DECAY 251 0. # G+ : WW ################################### ## INFORMATION FOR DIM6 ################################### Block DIM6 1 1.000000e+03 # Lambda 2 0.0 # cpDC 3 0.0 # cpWB 4 0.0 # cdp 5 0.0 # cp 6 0.0 # cWWW 7 0.0 # cG 8 0.0 # cpG 9 0.0 # cpW 10 0.0 # cpBB ################################### ## INFORMATION FOR DIM62F ################################### Block DIM62F 1 0.0 # cpl1 2 0.0 # cpl2 3 0.0 # cpl3 4 0.0 # c3pl1 5 0.0 # c3pl2 6 0.0 # c3pl3 7 0.0 # cpe 8 0.0 # cpmu 9 0.0 # cpta 10 0.0 # cpqMi 11 0.0 # cpq3i 12 0.0 # cpQ3 13 0.0 # cpQM 14 0.0 # cpu 15 0.0 # cpt 16 0.0 # cpd 19 0.0 # ctp 22 0.049973 # ctZ 23 0.761901 # ctW 24 0.0 # ctG ################################### ## INFORMATION FOR DIM64F ################################### Block DIM64F 1 0.0 # cQq83 2 0.0 # cQq81 3 0.0 # cQu8 4 0.0 # ctq8 6 0.0 # cQd8 7 0.0 # ctu8 8 0.0 # ctd8 10 0.0 # cQq13 11 0.0 # cQq11 12 0.0 # cQu1 13 0.0 # ctq1 14 0.0 # cQd1 16 0.0 # ctu1 17 0.0 # ctd1 19 0.0 # cQQ8 20 0.0 # cQQ1 21 0.0 # cQt1 23 0.0 # ctt1 25 0.0 # cQt8 ################################### ## INFORMATION FOR DIM64F2L ################################### Block DIM64F2L 1 0.0 # cQlM1 2 0.0 # cQlM2 3 0.0 # cQl31 4 0.0 # cQl32 5 0.0 # cQe1 6 0.0 # cQe2 7 0.0 # ctl1 8 0.0 # ctl2 9 0.0 # cte1 10 0.0 # cte2 13 0.0 # cQlM3 14 0.0 # cQl33 15 0.0 # cQe3 16 0.0 # ctl3 17 0.0 # cte3 19 0.0 # ctlS3 20 0.0 # ctlT3 21 0.0 # cblS3 ################################### ## INFORMATION FOR DIM64F4L ################################### Block DIM64F4L 1 0.0 # cll1111 2 0.0 # cll2222 3 0.0 # cll3333 4 0.0 # cll1122 5 0.0 # cll1133 6 0.0 # cll2233 7 0.0 # cll1221 8 0.0 # cll1331 9 0.0 # cll2332 ################################### ## INFORMATION FOR LOOP ################################### Block LOOP 1 9.118800e+01 # MU_R ################################### ## INFORMATION FOR RENOR ################################### Block Renor 1 5.550000e+02 # mueft ################################### ## INFORMATION FOR YUKAWA ################################### Block YUKAWA 6 1.720000e+02 # ymt #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) #=========================================================== Block QNUMBERS 9000001 # ghA 1 0 # 3 times electric charge 2 -1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000002 # ghZ 1 0 # 3 times electric charge 2 -1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000003 # ghWp 1 3 # 3 times electric charge 2 -1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000004 # ghWm 1 -3 # 3 times electric charge 2 -1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000005 # ghG 1 0 # 3 times electric charge 2 -1 # number of spin states (2S+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 250 # G0 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 251 # G+ 1 3 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) PROC_CARD tWa #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 3.3.2 2022-03-18 * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* * #************************************************************ #* * #* Command File for MadGraph5_aMC@NLO * #* * #* run as ./bin/mg5_aMC filename * #* * #************************************************************ set group_subprocesses Auto set ignore_six_quark_processes False set low_mem_multicore_nlo_generation False set complex_mass_scheme False set include_lepton_initiated_processes False set gauge unitary set loop_optimized_output True set loop_color_flows False set max_npoint_for_channel 0 set default_unset_couplings 99 set max_t_for_channel 99 set zerowidth_tchannel True set nlo_mixed_expansion True 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 SMEFTatNLO-NLO define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors define j = p define w = w+ w- define tt = t t~ generate p p > tt w a QED=2 QCD=1 NP=2 [QCD] output twa_SMEFT_DR1 run_banner tWa

#5.3.3.2 ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL #### ###################################################################### ## ## ## Width set on Auto will be computed following the information ## ## present in the decay.py files of the model. ## ## See arXiv:1402.1178 for more details. ## ## ## ###################################################################### ################################### ## INFORMATION FOR DIM6 ################################### Block dim6 1 1.000000e+03 # Lambda ################################### ## INFORMATION FOR DIM62F ################################### Block dim62f # 22 4.997300e-02 # ctZ # 23 7.619010e-01 # ctW 22 0.0 # ctZ 23 0.0 # ctW ################################### ## INFORMATION FOR MASS ################################### Block mass 6 1.720000e+02 # MT 23 9.118760e+01 # MZ 24 7.982440e+01 # MW 25 1.250000e+02 # MH ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. 1 0.000000e+00 # d : 0.0 2 0.000000e+00 # u : 0.0 3 0.000000e+00 # s : 0.0 4 0.000000e+00 # c : 0.0 5 0.000000e+00 # b : 0.0 11 0.000000e+00 # e- : 0.0 12 0.000000e+00 # ve : 0.0 13 0.000000e+00 # mu- : 0.0 14 0.000000e+00 # vm : 0.0 15 0.000000e+00 # ta- : 0.0 16 0.000000e+00 # vt : 0.0 21 0.000000e+00 # g : 0.0 22 0.000000e+00 # a : 0.0 9000002 9.118760e+01 # ghz : MZ 9000003 7.982440e+01 # ghwp : MW 9000004 7.982440e+01 # ghwm : MW ################################### ## INFORMATION FOR RENOR ################################### Block renor 1 5.550000e+02 # mueft ################################### ## INFORMATION FOR SMINPUTS ################################### Block sminputs 2 1.166370e-05 # Gf 3 1.184000e-01 # aS (Note that Parameter not used if you use a PDF set) ################################### ## INFORMATION FOR YUKAWA ################################### Block yukawa 6 1.720000e+02 # ymt ################################### ## INFORMATION FOR DECAY ################################### #DECAY 6 1.508336e+00 # WT #DECAY 23 2.495200e+00 # WZ DECAY 6 0.0 DECAY 23 0.0 ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. DECAY 1 0.000000e+00 # d : 0.0 DECAY 2 0.000000e+00 # u : 0.0 DECAY 3 0.000000e+00 # s : 0.0 DECAY 4 0.000000e+00 # c : 0.0 DECAY 5 0.000000e+00 # b : 0.0 DECAY 11 0.000000e+00 # e- : 0.0 DECAY 12 0.000000e+00 # ve : 0.0 DECAY 13 0.000000e+00 # mu- : 0.0 DECAY 14 0.000000e+00 # vm : 0.0 DECAY 15 0.000000e+00 # ta- : 0.0 DECAY 16 0.000000e+00 # vt : 0.0 DECAY 21 0.000000e+00 # g : 0.0 DECAY 22 0.000000e+00 # a : 0.0 DECAY 24 0.000000e+00 # w+ : 0.0 DECAY 25 0.000000e+00 # h : 0.0 #DECAY 9000002 2.495200e+00 # ghz : WZ DECAY 9000002 0.000000e+00 # ghz : WZ #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) #=========================================================== Block QNUMBERS 9000001 # gha 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000002 # ghz 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000003 # ghwp 1 3 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000004 # ghwm 1 -3 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000005 # ghg 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) order = NLO fixed_order = ON shower = OFF madspin = OFF reweight = OFF madanalysis = OFF runshower = False /FixedOrderAnalysis/ ( is the name of the exported # process directory). See the # /FixedOrderAnalysis/analysis_*_template.f file for details # on how to write your own analysis. # ####################################################################### # # Analysis format. # Can either be 'topdrawer', 'root', 'HwU', 'LHE' or 'none'. # When choosing HwU, it comes with a GnuPlot wrapper. When choosing # topdrawer, the histogramming package 'dbook.f' is included in the # code, while when choosing root the 'rbook_fe8.f' and 'rbook_be8.cc' # are included. If 'none' is chosen, all the other entries below have # to be set empty. FO_ANALYSIS_FORMAT = HwU # # # Needed extra-libraries (FastJet is already linked): FO_EXTRALIBS = # # (Absolute) path to the extra libraries. Directory names should be # separated by white spaces. FO_EXTRAPATHS = # # (Absolute) path to the dirs containing header files needed by the # libraries (e.g. C++ header files): FO_INCLUDEPATHS = # # User's analysis (to be put in the /FixedOrderAnalysis/ # directory). Please use .o as extension and white spaces to separate # files. FO_ANALYSE = analysis_HwU_twa.o # # ## When linking with root, the following settings are a working ## example on lxplus (CERN) as of July 2014. When using this, comment ## out the lines above and replace with the physical ## path to root, ## e.g. /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.11/x86_64-slc6-gcc46-dbg/root/ #FO_ANALYSIS_FORMAT = root #FO_EXTRALIBS = Core Cint Hist Matrix MathCore RIO dl Thread #FO_EXTRAPATHS = /lib #FO_INCLUDEPATHS = /include #FO_ANALYSE = analysis_root_template.o ]]>

PROC_CARD tqa #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 3.3.2 2022-03-18 * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* * #************************************************************ #* * #* Command File for MadGraph5_aMC@NLO * #* * #* run as ./bin/mg5_aMC filename * #* * #************************************************************ set group_subprocesses Auto set ignore_six_quark_processes False set low_mem_multicore_nlo_generation False set complex_mass_scheme False set include_lepton_initiated_processes False set gauge unitary set loop_optimized_output True set loop_color_flows False set max_npoint_for_channel 0 set default_unset_couplings 99 set max_t_for_channel 99 set zerowidth_tchannel True set nlo_mixed_expansion True 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 SMEFTatNLO-NLO define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors define j = p generate p p > t j a NP^2=4 QED^2=9 QCD^2=0 [all=QCD] output tqa_SMEFT_sq_ord run_banner tqa

#5.3.3.2 ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL #### ###################################################################### ## ## ## Width set on Auto will be computed following the information ## ## present in the decay.py files of the model. ## ## See arXiv:1402.1178 for more details. ## ## ## ###################################################################### ################################### ## INFORMATION FOR DIM6 ################################### Block dim6 1 1.000000e+03 # Lambda ################################### ## INFORMATION FOR DIM62F ################################### Block dim62f 22 4.997300e-02 # ctZ 23 7.619010e-01 # ctW ################################### ## INFORMATION FOR MASS ################################### Block mass 6 1.720000e+02 # MT 23 9.118760e+01 # MZ 24 7.982440e+01 # MW 25 1.250000e+02 # MH ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. 1 0.000000e+00 # d : 0.0 2 0.000000e+00 # u : 0.0 3 0.000000e+00 # s : 0.0 4 0.000000e+00 # c : 0.0 5 0.000000e+00 # b : 0.0 11 0.000000e+00 # e- : 0.0 12 0.000000e+00 # ve : 0.0 13 0.000000e+00 # mu- : 0.0 14 0.000000e+00 # vm : 0.0 15 0.000000e+00 # ta- : 0.0 16 0.000000e+00 # vt : 0.0 21 0.000000e+00 # g : 0.0 22 0.000000e+00 # a : 0.0 9000002 9.118760e+01 # ghz : MZ 9000003 7.982440e+01 # ghwp : MW 9000004 7.982440e+01 # ghwm : MW ################################### ## INFORMATION FOR RENOR ################################### Block renor 1 5.550000e+02 # mueft ################################### ## INFORMATION FOR SMINPUTS ################################### Block sminputs 2 1.166370e-05 # Gf 3 1.184000e-01 # aS (Note that Parameter not used if you use a PDF set) ################################### ## INFORMATION FOR YUKAWA ################################### Block yukawa 6 1.720000e+02 # ymt ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.508336e+00 # WT DECAY 23 2.495200e+00 # WZ ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. DECAY 1 0.000000e+00 # d : 0.0 DECAY 2 0.000000e+00 # u : 0.0 DECAY 3 0.000000e+00 # s : 0.0 DECAY 4 0.000000e+00 # c : 0.0 DECAY 5 0.000000e+00 # b : 0.0 DECAY 11 0.000000e+00 # e- : 0.0 DECAY 12 0.000000e+00 # ve : 0.0 DECAY 13 0.000000e+00 # mu- : 0.0 DECAY 14 0.000000e+00 # vm : 0.0 DECAY 15 0.000000e+00 # ta- : 0.0 DECAY 16 0.000000e+00 # vt : 0.0 DECAY 21 0.000000e+00 # g : 0.0 DECAY 22 0.000000e+00 # a : 0.0 DECAY 24 0.000000e+00 # w+ : 0.0 DECAY 25 0.000000e+00 # h : 0.0 DECAY 9000002 2.495200e+00 # ghz : WZ #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) #=========================================================== Block QNUMBERS 9000001 # gha 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000002 # ghz 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000003 # ghwp 1 3 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000004 # ghwm 1 -3 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) Block QNUMBERS 9000005 # ghg 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) order = NLO fixed_order = ON shower = OFF madspin = OFF reweight = OFF madanalysis = OFF runshower = False /FixedOrderAnalysis/ ( is the name of the exported # process directory). See the # /FixedOrderAnalysis/analysis_*_template.f file for details # on how to write your own analysis. # ####################################################################### # # Analysis format. # Can either be 'topdrawer', 'root', 'HwU', 'LHE' or 'none'. # When choosing HwU, it comes with a GnuPlot wrapper. When choosing # topdrawer, the histogramming package 'dbook.f' is included in the # code, while when choosing root the 'rbook_fe8.f' and 'rbook_be8.cc' # are included. If 'none' is chosen, all the other entries below have # to be set empty. FO_ANALYSIS_FORMAT = HwU # # # Needed extra-libraries (FastJet is already linked): FO_EXTRALIBS = # # (Absolute) path to the extra libraries. Directory names should be # separated by white spaces. FO_EXTRAPATHS = # # (Absolute) path to the dirs containing header files needed by the # libraries (e.g. C++ header files): FO_INCLUDEPATHS = # # User's analysis (to be put in the /FixedOrderAnalysis/ # directory). Please use .o as extension and white spaces to separate # files. FO_ANALYSE = analysis_HwU_template.o # # ## When linking with root, the following settings are a working ## example on lxplus (CERN) as of July 2014. When using this, comment ## out the lines above and replace with the physical ## path to root, ## e.g. /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.11/x86_64-slc6-gcc46-dbg/root/ #FO_ANALYSIS_FORMAT = root #FO_EXTRALIBS = Core Cint Hist Matrix MathCore RIO dl Thread #FO_EXTRAPATHS = /lib #FO_INCLUDEPATHS = /include #FO_ANALYSE = analysis_root_template.o ]]>

When	Who	Action	New State
2023-06-29	marco zaro	Request for more information	Needs information
2023-06-29	Olivier Mattelaer	Answer	Answered
2023-06-29	Michele Mormile	Reopen	Open

MadGraph5_aMC@NLO

History of question #707151