2.5.5 #*********************************************************************** # 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. * #*********************************************************************** # #******************* # 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 ! average or sum #*********************************************************************** # 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!!!! * #*********************************************************************** PYTHIA8 = 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 * #*********************************************************************** False = 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.4 = jetradius ! The radius parameter for the jet algorithm 20.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 * #*********************************************************************** 130.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) #*********************************************************************** # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* #*********************************************************************** 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) #*********************************************************************** ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ################################### ## INFORMATION FOR CKMBLOCK ################################### BLOCK CKMBLOCK # 1 2.277360e-01 # cabi ################################### ## INFORMATION FOR DMINPUTS ################################### BLOCK DMINPUTS # 1 0.000000e+00 # gvxc 2 0.000000e+00 # gvxd 3 1.000000e+00 # gaxd 4 0.000000e+00 # gvd11 5 0.000000e+00 # gvu11 6 0.000000e+00 # gvd22 7 0.000000e+00 # gvu22 8 0.000000e+00 # gvd33 9 0.000000e+00 # gvu33 10 0.000000e+00 # gvl11 11 0.000000e+00 # gvl22 12 0.000000e+00 # gvl33 13 2.500000e-01 # gad11 14 2.500000e-01 # gau11 15 2.500000e-01 # gad22 16 2.500000e-01 # gau22 17 2.500000e-01 # gad33 18 2.500000e-01 # gau33 19 0.000000e+00 # gal11 20 0.000000e+00 # gal22 21 0.000000e+00 # gal33 22 0.000000e+00 # gnu11 23 0.000000e+00 # gnu22 24 0.000000e+00 # gnu33 25 0.000000e+00 # gvu31 26 2.500000e-01 # gau31 27 0.000000e+00 # gvd31 28 2.500000e-01 # gad31 29 0.000000e+00 # gvh ################################### ## INFORMATION FOR LOOP ################################### BLOCK LOOP # 1 9.118800e+01 # mu_r ################################### ## INFORMATION FOR MASS ################################### BLOCK MASS # 6 1.720000e+02 # mt 15 1.777000e+00 # mta 23 9.118760e+01 # mz 25 1.250000e+02 # mh 5000001 8.000000e+02 # my1 5000511 1.000000e+01 # mxr 5000512 1.000000e+01 # mxc 5000521 1.000000e+01 # mxd 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 16 0.000000e+00 # vt : 0.0 21 0.000000e+00 # g : 0.0 22 0.000000e+00 # a : 0.0 24 7.982436e+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.982436e+01 # ghwp : mw 9000004 7.982436e+01 # ghwm : mw ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.279000e+02 # aewm1 2 1.166370e-05 # gf 3 1.184000e-01 # as ################################### ## INFORMATION FOR YUKAWA ################################### BLOCK YUKAWA # 6 1.720000e+02 # ymt 15 1.777000e+00 # ymtau ################################### ## INFORMATION FOR DECAY ################################### BLOCK QNUMBERS 5000511 # xr 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) ################################### ## INFORMATION FOR QNUMBERS 5000512 ################################### BLOCK QNUMBERS 5000512 # xc 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 5000521 ################################### BLOCK QNUMBERS 5000521 # xd 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 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 5000001 ################################### BLOCK QNUMBERS 5000001 # y1 1 0 # 3 times electric charge 2 3 # 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.508336e+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.070000e-03 # # PDG Width DECAY 5000001 4.501948e+01 # BR NDA ID1 ID2 ... 4.709243e-01 2 5000521 -5000521 # 21.2007678199 8.838115e-02 2 5 -5 # 3.9788735773 8.838115e-02 2 4 -4 # 3.9788735773 8.838115e-02 2 3 -3 # 3.9788735773 8.838115e-02 2 2 -2 # 3.9788735773 8.838115e-02 2 1 -1 # 3.9788735773 8.716993e-02 2 6 -6 # 3.92434477843 # # PDG Width DECAY 5000511 0.000000e+00 # # PDG Width DECAY 5000512 0.000000e+00 # # PDG Width DECAY 5000521 0.000000e+00 # # PDG Width DECAY 9000002 2.495200e+00 # # PDG Width DECAY 9000003 2.085000e+00 # # PDG Width DECAY 9000004 2.085000e+00 #*********************************************************************** # MadGraph5_aMC@NLO * # * # shower_card.dat aMC@NLO * # * # This file is used to set the parameters for the shower. * # * # Some notation/conventions: * # * # Lines starting with a hash (#) are info or comments * # * # mind the format: variable = value # comment * #*********************************************************************** # #*********************************************************************** # Shower settings * #*********************************************************************** # Number of events, jobs, errors, and random seeds * #*********************************************************************** nevents = -1 # N evts to shower (< 0 = all) nsplit_jobs = 1 # N jobs to run in parallel (< 100!!) combine_td = T # combine the topdrawer/HwU files if nsplit_jobs>1 maxprint = 2 # N evts to print in the log maxerrs = 0.1 # max fraction of errors rnd_seed = 0 # 1st random seed (0 = default) rnd_seed2 = 0 # 2nd random seed (0 = default) !ONLY FOR HWERIG6! #*********************************************************************** # PDFs and non-perturbative modelling * #*********************************************************************** pdfcode = 0 # 0 = internal, 1 = same as NLO, other = lhaglue ue_enabled = T # underlying event hadronize = T # hadronisation on/off !IGNORED BY HERWIG6! lambda_5 = -1 # Lambda_5 (< 0 = default) !IGNORED BY PYTHIA8! #*********************************************************************** # Stable or unstable particles * #*********************************************************************** b_stable = F # set B hadrons stable pi_stable = T # set pi0's stable wp_stable = F # set w+'s stable wm_stable = F # set w-'s stable z_stable = F # set z0's stable h_stable = F # set Higgs' stable tap_stable = F # set tau+'s stable tam_stable = F # set tau-'s stable mup_stable = F # set mu+'s stable mum_stable = F # set mu-'s stable #*********************************************************************** # Mass of the b quark * #*********************************************************************** b_mass = -1 # if < 0 = read from SubProcesses/MCmasses_*.inc #*********************************************************************** # Special settings * #*********************************************************************** is_4lep = F # T if 4-lepton production !ONLY FOR PYTHIA6! is_bbar = F # T if bb~ production !ONLY FOR HERWIG6! #*********************************************************************** # FxFx merging parameters !ONLY FOR PYTHIA8! #*********************************************************************** Qcut = -1.0 # Merging scale njmax = 0 # Maximal multiplicity in the merging #*********************************************************************** # Decay channels * #*********************************************************************** # Syntax for HERWIG6 * # DM_I = M > D1 D2 @ BR @ ME * # corresponding to call to HWMODK(M,BR,ME,D1,D2) * # I < 100, M is the decaying resonance, D1, D2, ... are the decay * # products (up to five), BR is the branching ratio and ME is the type * # of matrix element to be used in the decay. * # BR's are correctly understood only if they add up to 1, and only if * # no more than three modes are required for a given resonance. * # WARNING: the order of decay products in > 2-body decays IS RELEVANT. * # * # Syntax for PYTHIA6 * # DM_I = M > D1 D2 @ BR @ ME * # WARNING: turning hadronisation off disables top decays * # WARNING: 1 -> n decays (with n > 2) are handled through a sequence * # of 1 -> 2 decays. * # WARNING: entries BR and ME are ignored * # * # Syntax for HERWIG++ * # DM_I = M > D1 D2 @ BR @ ME * # WARNING: entries BR and ME are ignored * # * # Syntax for PYTHIA8 * # DM_I = M:onIfAny = D1 D2 * # or similar, according to the offical PYTHIA8 decay syntax, see * # the online PYTHIA8 manual * # WARNING: 1 -> n decays (with n > 2) are handled through a sequence * # of 1 -> 2 decays. * # * # Examples * # Z -> e+ e- or mu+ mu- with BR = 0.5 each, HERWIG6 * # DM_1 = 23 > -11 11 @ 0.5d0 @ 100 # DM_2 = 23 > -13 13 @ 0.5d0 @ 100 # H -> ta+ ta- with BR = 1, HERWIG6 or HERWIG++ * # DM_3 = 25 > -15 15 @ 1.0d0 @ 0 # t -> ve e+ b with BR = 1, HERWIG6 or HERWIG++ * # DM_4 = 6 > 12 -11 5 @ 1d0 @ 100 # t -> ve e+ b with BR = 1, PYTHIA6 * # DM_5 = 6 > 24 5 @ 1d0 @ 100 # DM_6 = 24 > 12 -11 @ 1d0 @ 100 # W+ -> ve e+, W- -> vm~ mu-, PYTHIA8 * # DM_1 = 24:onMode = off # DM_2 = 24:onPosIfAny = 11 12 # DM_3 = 24:onNegIfAny = 13 14 # W+ -> ve e+, W- -> ve~ e- and vm~ mu-, PYTHIA8 * # DM_1 = 24:onMode = off # DM_2 = 24:onIfAny = 11 12 # DM_3 = 24:onNegIfAny = 13 14 #*********************************************************************** # Extra libraries/analyses * #*********************************************************************** # The following lines need to be changed if the user does not want to * # create a StdHEP/HepMC file, but to directly run an own analysis (to * # be placed in HWAnalyzer or analogous MCatNLO subfolders). * # Please use files in those folders as examples. * # Remember that if your analysis uses hbook or is in the HwU format, * # you must also add to hbook.o or HwU.o to the ANALYSE list as well. * #*********************************************************************** extralibs = pythia8 boost_iostreams z dl stdc++ # Extra-libraries (not LHAPDF) # Default: "stdhep Fmcfio" # PYTHIA > 8.200 may require library dl extrapaths = ../lib /home/krolb/MG5_aMC_v2_5_5/HEPTools/pythia8/lib /usr/lib /home/krolb/MG5_aMC_v2_5_5/HEPTools/zlib/lib # Path to the extra-libraries # Default: "../lib" INCLUDEPATHS = # Path to header files needed by c++ # Dir names separated by white spaces ANALYSE = # User's analysis and histogramming # routines (please use .o as extension # and use spaces to separate files) fixed_order = OFF reweight = OFF madspin = OFF shower = ON madanalysis5 = Not installed order = NLO