#************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.5 2018-02-03 * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* and * #* http://amcatnlo.cern.ch * #* * #************************************************************ #* * #* Command File for aMCatNLO * #* * #* run as ./bin/aMCatNLO.py filename * #* * #************************************************************ reweight /afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/ttW_SMEFTNLOQED1_v1/Events/run_01/events.lhe.gz_0.lhe -from_cards --multicore=create Traceback (most recent call last): File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd return self.onecmd_orig(line, **opt) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig return func(arg, **opt) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 2142, in do_reweight reweight_cmd.import_command_file(path) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1660, in import_command_file self.exec_cmd(line, precmd=True) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd stop = Cmd.onecmd_orig(current_interface, line, **opt) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig return func(arg, **opt) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/various/misc.py", line 100, in f_with_no_logger out = f(self, *args, **opt) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 546, in do_launch weight = self.calculate_weight(event) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 954, in calculate_weight return self.calculate_nlo_weight(event) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 1028, in calculate_nlo_weight w_origV = self.calculate_matrix_element(cevent, 'V0', scale2=scale2) File "/afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 1155, in calculate_matrix_element orig_order, Pdir, hel_dict = self.id_to_path[tag] KeyError: (((-3, 4), (-6, 6, 24)), 'V') Value of current Options: text_editor : None web_browser : None cluster_temp_path : None timeout : 60 cluster_local_path : None cluster_queue : None madanalysis_path : None lhapdf : lhapdf-config cluster_size : 100 cluster_memory : None pjfry : None cluster_status_update : (600, 30) cluster_time : None f2py_compiler : None ninja : /afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/HEPTools/lib hepmc_path : None pythia8_path : None hwpp_path : None automatic_html_opening : True cluster_retry_wait : 300 stdout_level : None pythia-pgs_path : None mg5_path : /afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5 td_path : None delphes_path : None thepeg_path : None cluster_type : condor madanalysis5_path : None exrootanalysis_path : None fortran_compiler : None nb_core : 64 collier : /afs/crc.nd.edu/user/r/rgoldouz/gridpacks/MadGraph5/MG5_aMC_v2_6_5/HEPTools/lib auto_update : 7 cluster_nb_retry : 1 eps_viewer : None syscalc_path : None golem : None cpp_compiler : None notification_center : True run_mode : 2 #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.5 2018-02-03 * #* * #* 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 default_unset_couplings 99 set group_subprocesses Auto set ignore_six_quark_processes False set loop_optimized_output True set loop_color_flows False set gauge unitary set complex_mass_scheme False set max_npoint_for_channel 0 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_U2_2_U3_3_cG_4F_LO_UFO-NLO define j = p define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 define WW = w+ w- generate p p > t t~ WW QCD=2 QED=1 NP=2 [QCD] output ttW_SMEFTNLOQED1_v1 -nojpeg ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ################################### ## INFORMATION FOR DIM6 ################################### BLOCK DIM6 # 1 1.000000e+03 # lambda 2 2.000000e-05 # cpdc 3 3.000000e-05 # cpwb 4 4.000000e-05 # cdp 5 5.000000e-05 # cp 6 6.000000e-05 # cwww 8 8.000000e-05 # cpg 9 9.000000e-05 # cpw 10 1.000000e-04 # cpbb ################################### ## INFORMATION FOR DIM62F ################################### BLOCK DIM62F # 1 1.000000e-05 # cpl1 2 2.000000e-05 # cpl2 3 3.000000e-05 # cpl3 4 4.000000e-05 # c3pl1 5 5.000000e-05 # c3pl2 6 6.000000e-05 # c3pl3 7 7.000000e-05 # cpe 8 8.000000e-05 # cpmu 9 9.000000e-05 # cpta 10 1.000000e-04 # cpqmi 11 1.100000e-04 # cpq3i 12 1.000000e+01 # cpq3 13 1.000000e+01 # cpqm 14 1.400000e-04 # cpu 15 1.000000e+01 # cpt 16 1.600000e-04 # cpd 19 1.000000e+01 # ctp 22 1.000000e+01 # ctz 23 1.000000e+01 # ctw 24 1.000000e+00 # ctg ################################### ## INFORMATION FOR DIM64F4L ################################### BLOCK DIM64F4L # 1 1.000000e-05 # cll1111 2 2.000000e-05 # cll2222 3 3.000000e-05 # cll3333 4 4.000000e-05 # cll1122 5 5.000000e-05 # cll1133 6 6.000000e-05 # cll2233 7 7.000000e-05 # cll1221 8 8.000000e-05 # cll1331 9 9.000000e-05 # cll2332 ################################### ## INFORMATION FOR LOOP ################################### BLOCK LOOP # 1 9.118800e+01 # mu_r ################################### ## 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 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 # mreno ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.325070e+02 # aewm1 2 1.166370e-05 # gf 3 1.184000e-01 # as ################################### ## INFORMATION FOR YUKAWA ################################### BLOCK YUKAWA # 6 1.720000e+02 # ymt ################################### ## INFORMATION FOR DECAY ################################### DECAY 23 2.495200e+00 # wz DECAY 24 2.085000e+00 # ww DECAY 25 5.753088e-03 # wh 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 6 0.000000e+00 # t : 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 9000002 2.495200e+00 # ghz : wz DECAY 9000003 2.085000e+00 # ghwp : ww DECAY 9000004 2.085000e+00 # ghwm : ww ################################### ## INFORMATION FOR QNUMBERS 9000001 ################################### 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) ################################### ## INFORMATION FOR QNUMBERS 9000002 ################################### 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) ################################### ## INFORMATION FOR QNUMBERS 9000003 ################################### 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) ################################### ## INFORMATION FOR QNUMBERS 9000004 ################################### 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) ################################### ## INFORMATION FOR QNUMBERS 9000005 ################################### 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) #*********************************************************************** # MadGraph5_aMC@NLO * # * # run_card.dat aMC@NLO * # * # This file is used to set the parameters of the run. * # * # Some notation/conventions: * # * # Lines starting with a hash (#) are info or comments * # * # mind the format: value = variable ! comment * # * # Some of the values of variables can be list. These can either be * # comma or space separated. * # * # To display additional parameter, you can use the command: * # update to_full * #*********************************************************************** # #******************* # Running parameters #******************* # #*********************************************************************** # Tag name for the run (one word) * #*********************************************************************** tag_1 = run_tag ! name of the run #*********************************************************************** # Number of LHE events (and their normalization) and the required * # (relative) accuracy on the Xsec. * # These values are ignored for fixed order runs * #*********************************************************************** 10000 = nevents ! Number of unweighted events requested -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents) -1 = nevt_job ! Max number of events per job in event generation. ! (-1= no split). #*********************************************************************** # Normalize the weights of LHE events such that they sum or average to * # the total cross section * #*********************************************************************** average = event_norm ! valid settings: average, sum, bias #*********************************************************************** # Number of points per itegration channel (ignored for aMC@NLO runs) * #*********************************************************************** 0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the ! number of points and iter. below) # These numbers are ignored except if req_acc_FO is equal to -1 5000 = npoints_fo_grid ! number of points to setup grids 4 = niters_fo_grid ! number of iter. to setup grids 10000 = npoints_fo ! number of points to compute Xsec 6 = niters_fo ! number of iter. to compute Xsec #*********************************************************************** # Random number seed * #*********************************************************************** 0 = iseed ! rnd seed (0=assigned automatically=default)) #*********************************************************************** # Collider type and energy * #*********************************************************************** 1 = lpp1 ! beam 1 type (0 = no PDF) 1 = lpp2 ! beam 2 type (0 = no PDF) 6500.0 = ebeam1 ! beam 1 energy in GeV 6500.0 = ebeam2 ! beam 2 energy in GeV #*********************************************************************** # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * #*********************************************************************** nn23nlo = pdlabel ! PDF set 244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only ! numbers for central PDF sets are allowed. Can be a list; ! PDF sets beyond the first are included via reweighting. #*********************************************************************** # Include the NLO Monte Carlo subtr. terms for the following parton * # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * # WARNING: PYTHIA6PT works only for processes without FSR!!!! * #*********************************************************************** HERWIG6 = parton_shower 1.0 = shower_scale_factor ! multiply default shower starting ! scale by this factor #*********************************************************************** # Renormalization and factorization scales * # (Default functional form for the non-fixed scales is the sum of * # the transverse masses divided by two of all final state particles * # and partons. This can be changed in SubProcesses/set_scales.f or via * # dynamical_scale_choice option) * #*********************************************************************** False = fixed_ren_scale ! if .true. use fixed ren scale False = fixed_fac_scale ! if .true. use fixed fac scale 91.118 = mur_ref_fixed ! fixed ren reference scale 91.118 = muf_ref_fixed ! fixed fact reference scale -1 = dynamical_scale_choice ! Choose one (or more) of the predefined ! dynamical choices. Can be a list; scale choices beyond the ! first are included via reweighting 1.0 = mur_over_ref ! ratio of current muR over reference muR 1.0 = muf_over_ref ! ratio of current muF over reference muF #*********************************************************************** # Reweight variables for scale dependence and PDF uncertainty * #*********************************************************************** 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting True = reweight_scale ! Reweight to get scale variation using the ! rw_rscale and rw_fscale factors. Should be a list of ! booleans of equal length to dynamical_scale_choice to ! specify for which choice to include scale dependence. False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a ! list booleans of equal length to lhaid to specify for ! which PDF set to include the uncertainties. #*********************************************************************** # Store reweight information in the LHE file for off-line model- * # parameter reweighting at NLO+PS accuracy * #*********************************************************************** True = store_rwgt_info ! Store info for reweighting in LHE file #*********************************************************************** # ickkw parameter: * # 0: No merging * # 3: FxFx Merging - WARNING! Applies merging only at the hard-event * # level. After showering an MLM-type merging should be applied as * # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * # 4: UNLOPS merging (with pythia8 only). No interface from within * # MG5_aMC available, but available in Pythia8. * # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * #*********************************************************************** 0 = ickkw #*********************************************************************** # #*********************************************************************** # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are * # written in the LHE event file * #*********************************************************************** 15.0 = bwcutoff #*********************************************************************** # Cuts on the jets. Jet clustering is performed by FastJet. * # - When matching to a parton shower, these generation cuts should be * # considerably softer than the analysis cuts. * # - More specific cuts can be specified in SubProcesses/cuts.f * #*********************************************************************** 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) 0.7 = jetradius ! The radius parameter for the jet algorithm 10.0 = ptj ! Min jet transverse momentum -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) #*********************************************************************** # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * # More specific cuts can be specified in SubProcesses/cuts.f * #*********************************************************************** 0.0 = ptl ! Min lepton transverse momentum -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) 0.0 = drll ! Min distance between opposite sign lepton pairs 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs #*********************************************************************** # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, * # all the other parameters are ignored. * # More specific cuts can be specified in SubProcesses/cuts.f * #*********************************************************************** 20.0 = ptgmin ! Min photon transverse momentum -1.0 = etagamma ! Max photon abs(pseudo-rap) 0.4 = r0gamma ! Radius of isolation code 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 True = isoem ! isolate photons from EM energy (photons and leptons) #*********************************************************************** # Cuts associated to MASSIVE particles identified by their PDG codes. * # All cuts are applied to both particles and anti-particles, so use * # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or * # {6:100, 25:200} for multiple particles * #*********************************************************************** {} = pt_min_pdg ! Min pT for a massive particle {} = pt_max_pdg ! Max pT for a massive particle {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles #*********************************************************************** # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* #*********************************************************************** 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) #***********************************************************************