#************************************************************ #* MadGraph5_aMC@NLO/MadEvent * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.3.2 2018-06-22 * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* * #************************************************************ #* * #* Command File for MadEvent * #* * #* run as ./bin/madevent.py filename * #* * #************************************************************ generate_events run_01 Traceback (most recent call last): File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1501, in onecmd return self.onecmd_orig(line, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events self.run_generate_events(switch_mode, args) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/common_run_interface.py", line 6711, in new_fct original_fct(obj, *args, **opts) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/madevent_interface.py", line 2562, in run_generate_events self.exec_cmd('decay_events -from_cards', postcmd=False) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd stop = Cmd.onecmd_orig(current_interface, line, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/common_run_interface.py", line 3653, in do_decay_events madspin_cmd.import_command_file(path) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1643, in import_command_file self.exec_cmd(line, precmd=True) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd stop = Cmd.onecmd_orig(current_interface, line, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/madgraph/various/misc.py", line 100, in f_with_no_logger out = f(self, *args, **opt) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/MadSpin/interface_madspin.py", line 573, in do_launch return self.run_bridge(line) File "/nfs/scratch/fynu/mlucente/mg5amcnlo/MadSpin/interface_madspin.py", line 990, in run_bridge nb_decay = len(evt_decayfile[particle.pdg]) KeyError: 13 Run Options ----------- stdout_level : 10 (user set) MadEvent Options ---------------- automatic_html_opening : False (user set) notification_center : True cluster_temp_path : None cluster_memory : None (user set) cluster_size : 150 (user set) cluster_queue : None nb_core : 5 (user set) cluster_time : 5 (user set) run_mode : 1 (user set) Configuration Options --------------------- text_editor : None cluster_local_path : /cvmfs/cp3.uclouvain.be/madgraph/ (user set) cluster_status_update : (900, 60) (user set) pythia8_path : /nfs/scratch/fynu/mlucente/mg5amcnlo/HEPTools/pythia8 (user set) hwpp_path : None (user set) pythia-pgs_path : None (user set) td_path : None (user set) delphes_path : /nfs/scratch/fynu/mlucente/mg5amcnlo/Delphes (user set) thepeg_path : None (user set) cluster_type : slurm (user set) madanalysis5_path : None (user set) cluster_nb_retry : 1 eps_viewer : None web_browser : None syscalc_path : None (user set) madanalysis_path : None (user set) lhapdf : /cvmfs/cp3.uclouvain.be/lhapdf/lhapdf-6.1.5_amd64_gcc44/bin/lhapdf-config (user set) f2py_compiler : None hepmc_path : None (user set) cluster_retry_wait : 300 (user set) fortran_compiler : None auto_update : 7 (user set) exrootanalysis_path : /nfs/scratch/fynu/mlucente/mg5amcnlo/ExRootAnalysis (user set) timeout : 60 cpp_compiler : None #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.3.2 2018-06-22 * #* BZR mg5amcnlo 279 * #* * #* 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_HeavyN_NLO 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~ define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors define j = p generate p p > mu- n1 output pp_crash ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ################################### ## INFORMATION FOR LOOP ################################### BLOCK LOOP # 1 9.118800e+01 # mu_r ################################### ## INFORMATION FOR MASS ################################### BLOCK MASS # 6 1.733000e+02 # mt 23 9.118760e+01 # mz 25 1.257000e+02 # mh 9900012 5.000000e+01 # mn1 9900014 1.000000e+06 # mn2 9900016 1.000000e+06 # mn3 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 24 7.995123e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2))) 9000002 9.118760e+01 # ghz : mz 9000003 7.995123e+01 # ghwp : mw 9000004 7.995123e+01 # ghwm : mw ################################### ## INFORMATION FOR NUMIXING ################################### BLOCK NUMIXING # 1 0.000000e+00 # ven1 2 0.000000e+00 # ven2 3 0.000000e+00 # ven3 4 1.000000e-05 # vmun1 5 0.000000e+00 # vmun2 6 0.000000e+00 # vmun3 7 0.000000e+00 # vtan1 8 0.000000e+00 # vtan2 9 0.000000e+00 # vtan3 ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.279400e+02 # aewm1 2 1.174560e-05 # gf 3 1.184000e-01 # as ################################### ## INFORMATION FOR YUKAWA ################################### BLOCK YUKAWA # 6 1.733000e+02 # ymt ################################### ## INFORMATION FOR DECAY ################################### 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 82 ################################### BLOCK QNUMBERS 82 # 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) ################################### ## INFORMATION FOR QNUMBERS 9900012 ################################### BLOCK QNUMBERS 9900012 # n1 1 0 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 9900014 ################################### BLOCK QNUMBERS 9900014 # n2 1 0 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 9900016 ################################### BLOCK QNUMBERS 9900016 # n3 1 0 # 3 times electric charge 2 2 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) # #************************* # Decay widths * #************************* # # PDG Width DECAY 1 0.000000e+00 # # PDG Width DECAY 2 0.000000e+00 # # PDG Width DECAY 3 0.000000e+00 # # PDG Width DECAY 4 0.000000e+00 # # PDG Width DECAY 5 0.000000e+00 # # PDG Width DECAY 6 1.350000e+00 # # PDG Width DECAY 11 0.000000e+00 # # PDG Width DECAY 12 0.000000e+00 # # PDG Width DECAY 13 0.000000e+00 # # PDG Width DECAY 14 0.000000e+00 # # PDG Width DECAY 15 0.000000e+00 # # PDG Width DECAY 16 0.000000e+00 # # PDG Width DECAY 21 0.000000e+00 # # PDG Width DECAY 22 0.000000e+00 # # PDG Width DECAY 23 2.495200e+00 # # PDG Width DECAY 24 2.085000e+00 # # PDG Width DECAY 25 4.170000e-03 # # PDG Width DECAY 9000002 2.495200e+00 # # PDG Width DECAY 9000003 2.085000e+00 # # PDG Width DECAY 9000004 2.085000e+00 # # PDG Width DECAY 9900012 2.237257e-14 # BR NDA ID1 ID2 ... 1.276474e-01 3 -13 -2 1 # 2.85580039182e-15 1.276474e-01 3 -13 -4 3 # 2.85580039182e-15 1.262349e-01 3 -1 2 13 # 2.82419913669e-15 1.262349e-01 3 -3 4 13 # 2.82419913669e-15 4.308848e-02 3 -11 12 13 # 9.64000034994e-16 4.308848e-02 3 -15 13 16 # 9.64000034994e-16 4.271078e-02 3 -13 -12 11 # 9.55549915305e-16 4.271078e-02 3 -16 -13 15 # 9.55549915305e-16 2.566178e-02 3 -13 13 14 # 5.74119969375e-16 2.565821e-02 3 -14 -13 13 # 5.740400993e-16 2.199271e-02 3 -1 1 14 # 4.92033443965e-16 2.199271e-02 3 -3 3 14 # 4.92033443965e-16 2.199271e-02 3 -5 5 14 # 4.92033443965e-16 2.197632e-02 3 -14 -1 1 # 4.91666757542e-16 2.197632e-02 3 -14 -3 3 # 4.91666757542e-16 2.197632e-02 3 -14 -5 5 # 4.91666757542e-16 1.698687e-02 3 -2 2 14 # 3.80039938156e-16 1.698687e-02 3 -4 4 14 # 3.80039938156e-16 1.695156e-02 3 -14 -2 2 # 3.79249962709e-16 1.695156e-02 3 -14 -4 4 # 3.79249962709e-16 1.002612e-02 3 -14 -12 12 # 2.24310071528e-16 1.002612e-02 3 -16 -14 16 # 2.24310071528e-16 9.841292e-03 3 -12 12 14 # 2.2017499416e-16 9.841292e-03 3 -16 14 16 # 2.2017499416e-16 5.001661e-03 3 -14 14 14 # 1.11900010839e-16 4.988028e-03 3 -11 11 14 # 1.11595005592e-16 4.988028e-03 3 -15 14 15 # 1.11595005592e-16 4.969702e-03 3 -14 -11 11 # 1.11185005874e-16 4.969702e-03 3 -15 -14 15 # 1.11185005874e-16 4.880977e-03 3 -14 -14 14 # 1.09199999601e-16 # # PDG Width DECAY 9900014 1.500000e+00 # # PDG Width DECAY 9900016 1.230000e+01 #********************************************************************* # MadGraph5_aMC@NLO * # * # run_card.dat MadEvent * # * # This file is used to set the parameters of the run. * # * # Some notation/conventions: * # * # Lines starting with a '# ' are info or comments * # * # mind the format: value = variable ! comment * # * # To display more options, you can type the command: * # update full_run_card * #********************************************************************* # #******************* # Running parameters #******************* # #********************************************************************* # Tag name for the run (one word) * #********************************************************************* tag_1 = run_tag ! name of the run #********************************************************************* # Number of events and rnd seed * # Warning: Do not generate more than 1M events in a single run * # If you want to run Pythia, avoid more than 50k events in a run. * #********************************************************************* 10000 = nevents ! Number of unweighted events requested 0 = iseed ! rnd seed (0=assigned automatically=default)) #********************************************************************* # Collider type and energy * # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * # 3=photon from electron * #********************************************************************* 1 = lpp1 ! beam 1 type 1 = lpp2 ! beam 2 type 7000.0 = ebeam1 ! beam 1 total energy in GeV 7000.0 = ebeam2 ! beam 2 total energy in GeV # To see polarised beam options: type "update beam_pol" #********************************************************************* # PDF CHOICE: this automatically fixes also alpha_s and its evol. * #********************************************************************* nn23lo1 = pdlabel ! PDF set 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number # To see heavy ion options: type "update ion_pdf" #********************************************************************* # Renormalization and factorization scales * #********************************************************************* False = fixed_ren_scale ! if .true. use fixed ren scale False = fixed_fac_scale ! if .true. use fixed fac scale 91.188 = scale ! fixed ren scale 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices 1.0 = scalefact ! scale factor for event-by-event scales #********************************************************************* # Type and output format #********************************************************************* False = gridpack !True = setting up the grid pack 1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written) 3.0 = lhe_version ! Change the way clustering information pass to shower. True = clusinfo ! include clustering tag in output average = event_norm ! average/sum. Normalization of the weight in the LHEF #********************************************************************* # Matching parameter (MLM only) #********************************************************************* 0 = ickkw ! 0 no matching, 1 MLM 1.0 = alpsfact ! scale factor for QCD emission vx False = chcluster ! cluster only according to channel diag 4 = asrwgtflavor ! highest quark flavor for a_s reweight False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 ! (turn off for VBF and single top processes) 0.0 = xqcut ! minimum kt jet measure between partons #********************************************************************* # #********************************************************************* # handling of the helicities: # 0: sum over all helicities # 1: importance sampling over helicities #********************************************************************* 0 = nhel ! using helicities importance sampling or not. #********************************************************************* # Generation bias, check the wiki page below for more information: * # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' * #********************************************************************* None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-] {} = bias_parameters ! Specifies the parameters of the module. # #******************************* # Parton level cuts definition * #******************************* # # #********************************************************************* # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay #********************************************************************* 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) #********************************************************************* # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) #********************************************************************* False = cut_decays ! Cut decay products #********************************************************************* # Standard Cuts * #********************************************************************* # Minimum and maximum pt's (for max, -1 means no cut) * #********************************************************************* 20.0 = ptj ! minimum pt for the jets 0.0 = ptb ! minimum pt for the b 10.0 = pta ! minimum pt for the photons 10.0 = ptl ! minimum pt for the charged leptons 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) -1.0 = ptjmax ! maximum pt for the jets -1.0 = ptbmax ! maximum pt for the b -1.0 = ptamax ! maximum pt for the photons -1.0 = ptlmax ! maximum pt for the charged leptons -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) #********************************************************************* # Minimum and maximum E's (in the center of mass frame) * #********************************************************************* 0.0 = ej ! minimum E for the jets 0.0 = eb ! minimum E for the b 0.0 = ea ! minimum E for the photons 0.0 = el ! minimum E for the charged leptons -1.0 = ejmax ! maximum E for the jets -1.0 = ebmax ! maximum E for the b -1.0 = eamax ! maximum E for the photons -1.0 = elmax ! maximum E for the charged leptons {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) #********************************************************************* # Maximum and minimum absolute rapidity (for max, -1 means no cut) * #********************************************************************* 5.0 = etaj ! max rap for the jets -1.0 = etab ! max rap for the b 2.5 = etaa ! max rap for the photons 2.5 = etal ! max rap for the charged leptons 0.0 = etajmin ! min rap for the jets 0.0 = etabmin ! min rap for the b 0.0 = etaamin ! min rap for the photons 0.0 = etalmin ! main rap for the charged leptons {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5}) #********************************************************************* # Minimum and maximum DeltaR distance * #********************************************************************* 0.4 = drjj ! min distance between jets 0.0 = drbb ! min distance between b's 0.4 = drll ! min distance between leptons 0.4 = draa ! min distance between gammas 0.0 = drbj ! min distance between b and jet 0.4 = draj ! min distance between gamma and jet 0.4 = drjl ! min distance between jet and lepton 0.0 = drab ! min distance between gamma and b 0.0 = drbl ! min distance between b and lepton 0.4 = dral ! min distance between gamma and lepton -1.0 = drjjmax ! max distance between jets -1.0 = drbbmax ! max distance between b's -1.0 = drllmax ! max distance between leptons -1.0 = draamax ! max distance between gammas -1.0 = drbjmax ! max distance between b and jet -1.0 = drajmax ! max distance between gamma and jet -1.0 = drjlmax ! max distance between jet and lepton -1.0 = drabmax ! max distance between gamma and b -1.0 = drblmax ! max distance between b and lepton -1.0 = dralmax ! maxdistance between gamma and lepton #********************************************************************* # Minimum and maximum invariant mass for pairs * # WARNING: for four lepton final state mmll cut require to have * # different lepton masses for each flavor! * #********************************************************************* 0.0 = mmjj ! min invariant mass of a jet pair 0.0 = mmbb ! min invariant mass of a b pair 0.0 = mmaa ! min invariant mass of gamma gamma pair 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair -1.0 = mmjjmax ! max invariant mass of a jet pair -1.0 = mmbbmax ! max invariant mass of a b pair -1.0 = mmaamax ! max invariant mass of gamma gamma pair -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250}) {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only ! to pairs of particle/antiparticle and not to pairs of the same pdg codes. #********************************************************************* # Minimum and maximum invariant mass for all letpons * #********************************************************************* 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) #********************************************************************* # Minimum and maximum pt for 4-momenta sum of leptons * #********************************************************************* 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) #********************************************************************* # Inclusive cuts * #********************************************************************* 0.0 = ptheavy ! minimum pt for at least one heavy final state 0.0 = xptj ! minimum pt for at least one jet 0.0 = xptb ! minimum pt for at least one b 0.0 = xpta ! minimum pt for at least one photon 0.0 = xptl ! minimum pt for at least one charged lepton #********************************************************************* # Control the pt's of the jets sorted by pt * #********************************************************************* 0.0 = ptj1min ! minimum pt for the leading jet in pt 0.0 = ptj2min ! minimum pt for the second jet in pt 0.0 = ptj3min ! minimum pt for the third jet in pt 0.0 = ptj4min ! minimum pt for the fourth jet in pt -1.0 = ptj1max ! maximum pt for the leading jet in pt -1.0 = ptj2max ! maximum pt for the second jet in pt -1.0 = ptj3max ! maximum pt for the third jet in pt -1.0 = ptj4max ! maximum pt for the fourth jet in pt 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts #********************************************************************* # Control the pt's of leptons sorted by pt * #********************************************************************* 0.0 = ptl1min ! minimum pt for the leading lepton in pt 0.0 = ptl2min ! minimum pt for the second lepton in pt 0.0 = ptl3min ! minimum pt for the third lepton in pt 0.0 = ptl4min ! minimum pt for the fourth lepton in pt -1.0 = ptl1max ! maximum pt for the leading lepton in pt -1.0 = ptl2max ! maximum pt for the second lepton in pt -1.0 = ptl3max ! maximum pt for the third lepton in pt -1.0 = ptl4max ! maximum pt for the fourth lepton in pt #********************************************************************* # Control the Ht(k)=Sum of k leading jets * #********************************************************************* 0.0 = htjmin ! minimum jet HT=Sum(jet pt) -1.0 = htjmax ! maximum jet HT=Sum(jet pt) 0.0 = ihtmin !inclusive Ht for all partons (including b) -1.0 = ihtmax !inclusive Ht for all partons (including b) 0.0 = ht2min ! minimum Ht for the two leading jets 0.0 = ht3min ! minimum Ht for the three leading jets 0.0 = ht4min ! minimum Ht for the four leading jets -1.0 = ht2max ! maximum Ht for the two leading jets -1.0 = ht3max ! maximum Ht for the three leading jets -1.0 = ht4max ! maximum Ht for the four leading jets #*********************************************************************** # Photon-isolation cuts, according to hep-ph/9801442 * # When ptgmin=0, all the other parameters are ignored * # When ptgmin>0, pta and draj are not going to be used * #*********************************************************************** 0.0 = ptgmin ! Min photon transverse momentum 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) #********************************************************************* # WBF cuts * #********************************************************************* 0.0 = xetamin ! minimum rapidity for two jets in the WBF case 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case #*********************************************************************** # Turn on either the ktdurham or ptlund cut to activate * # CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] * #*********************************************************************** -1.0 = ktdurham 0.4 = dparameter -1.0 = ptlund 1, 2, 3, 4, 5, 6, 21, 82 = pdgs_for_merging_cut ! PDGs for two cuts above #********************************************************************* # maximal pdg code for quark to be considered as a light jet * # (otherwise b cuts are applied) * #********************************************************************* 4 = maxjetflavor ! Maximum jet pdg code #********************************************************************* # #********************************************************************* # Store info for systematics studies * # WARNING: Do not use for interference type of computation * #********************************************************************* False = use_syst ! Enable systematics studies # systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++] ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule # Syscalc is deprecated but to see the associate options type'update syscalc'#********************************************************************* # Additional hidden parameters #********************************************************************* ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments # Choose the argment to pass to the systematics command. like --mur=0.25,1,4. Look at the help of the systematics function for more details.