# MG/ME version : 4.4.57 # madgraph version : 4.4.5 # template version : 2.4.30 # helas version : 3.8 # model version : topBSM_1.4 #********************************************************************* # MadGraph/MadEvent * # http://madgraph.hep.uiuc.edu * # * # proc_card.dat * # * # This file is used to generate the code for a specific process. * # Some notation/conventions: * # * # 0. Do not modify the TAGS and their order. * # 1. hash/pound is a comment. * # 2. The number after the @ is used as an identifier for the class * # of processes. It can be any positive integer. * # 3. The number of lines for the max couplings depends on how many * # different classes of couplings are present in the model * # In the SM these are just two: QED (which include EW) and QCD * # 4. Write "end_coup" after the couplings list, * # to tell MG that the couplings input is over. * # 5. Write "done" after the proc list to * # to tell MG that the proc input is over. * # 6. Some model names available at present are: * # sm = Standard Model * # smckm = Standard Model with Cabibbo matrix * # mssm = Minimal Supersymmetric Standard Model * # 2hdm = Generic Two Higgs Doublet model * # heft = Higgs EFT (+Standard Model) * # usrmod = User Model * # 7. Don't leave spaces between the particles name in the * # definition of the multiparticles. * #********************************************************************* #********************************************************************* # Process(es) requested : mg2 input * #********************************************************************* # Begin PROCESS # This is TAG. Do not modify this line pp>o1>tt~ @1 # First Process QCD=99 # Max QCD couplings QED=99 # Max QED couplings QS1=99 end_coup # End the couplings input done # this tells MG there are no more procs # End PROCESS # This is TAG. Do not modify this line #********************************************************************* # Model information * #********************************************************************* # Begin MODEL # This is TAG. Do not modify this line topBSM # End MODEL # This is TAG. Do not modify this line #********************************************************************* # Start multiparticle definitions * #********************************************************************* # Begin MULTIPARTICLES # This is TAG. Do not modify this line P uu~cc~dd~ss~g J uu~cc~dd~ss~g L+ e+mu+ L- e-mu- vl vevm vl~ ve~vm~ # End MULTIPARTICLES # This is TAG. Do not modify this line #********************************************************************* # MadGraph/MadEvent * # http://madgraph.hep.uiuc.edu * # * # param_card.dat * # * # Les Houches friendly file for the SM parameters of MadGraph * # Spectrum and decay widths produced by SMCalc * # * # Please note the following IMPORTANT issues: * # * # 0. REFRAIN from editing this file by hand! Some of the parame- * # ters are not independent * # (such as G_Fermi, alpha_em, sin(theta_W),MZ,MW) and serious * # problems might be encountered (such as violation of unitarity * # or gauge invariance). Always use a calculator. * # * # 1. alpha_S(MZ) has been used in the calculation of the parameters * # but, for consistency, it will be reset by madgraph to the * # value expected IF the pdfs for collisions with hadrons are * # used. This value is KEPT by madgraph when no pdf are used * # lpp(i)=0 . * # * # 2. Values of the charm and bottom kinematic (pole) masses are * # those used in the matrix elements and phase space UNLESS they * # are set to ZERO from the start in the model (particles.dat) * # This happens, for example, when using 5-flavor QCD where * # charm and bottom are treated as partons in the initial state * # and a zero mass might be hardwired in the model definition. * # * ##******************************************************************** Block MODSEL # Select Model 0 1 # 0 1 = SM Block SMINPUTS # Standard Model inputs 1 1.27900000E+02 # alpha_em(MZ)(-1) SM MSbar 2 1.16600000E-05 # G_Fermi 3 1.18000000E-01 # alpha_s(MZ) SM MSbar 4 9.15445065E+01 # Z mass (as input parameter) Block MGSMPARAM # Standard Model parameters for MadGraph 1 2.31200000E-01 # sin(theta_W)^2 2 8.02673592E+01 # W mass (as input parameter) Block MGYUKAWA # Yukawa masses m/v=y/sqrt(2) # PDG YMASS 5 4.20000000E+00 # mbottom for the Yukawa y_b 4 1.25000000E+00 # mcharm for the Yukawa y_c 6 1.74300000E+02 # mtop for the Yukawa y_t 15 1.77700000E+00 # mtau for the Yukawa y_ta Block MGCKM # CKM elements for MadGraph 1 1 9.75000000E-01 # Vud for Cabibbo matrix #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) IF ANY # (see below for masses and decay tables) # These blocks are automatically created by the MadGraph # qnumbers.pl script from the particles.dat model file #=========================================================== BLOCK QNUMBERS 6000047 # s1 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) BLOCK QNUMBERS 6000048 # o1 1 0 # 3 times electric charge 2 3 # number of spin states (2S+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # Particle/Antiparticle distinction (0=own anti) # END of QNUMBERS blocks #=========================================================== Block MASS # Mass spectrum (kinematic masses) # PDG Mass 4 1.40000000E+00 # charm pole mass 5 4.20000000E+00 # bottom pole mass 6 1.74300000E+02 # top pole mass 15 1.77700000E+00 # tau mass 23 9.15445065E+01 # Z mass 24 8.02673592E+01 # W mass 25 1.00000000E+02 # H mass 6000045 4.00000000e+02 # S0MASS 6000046 4.00000000e+02 # O0MASS 6000047 2.00000000e+03 # S1MASS 6000048 0.22500000e+04 # O1MASS 6000049 9.00000000e+02 # S2MASS !should be equal to mstring below 6000050 6.00000000e+02 # g1mass 6000051 1.00000000e+02 # g2mass !not used... Calculated in couplings.f 6000052 1.00000000e+02 # g3mass !not used... Calculated in couplings.f 6000053 1.00000000e+02 # g4mass !not used... Calculated in couplings.f 6000054 1.00000000e+02 # g5mass !not used... Calculated in couplings.f 6000055 1.00000000e+02 # g6mass !not used... Calculated in couplings.f 6000056 1.00000000e+02 # g7mass !not used... Calculated in couplings.f 6000057 1.00000000e+02 # g8mass !not used... Calculated in couplings.f 6000058 1.00000000e+02 # g9mass !not used... Calculated in couplings.f 6000059 1.00000000e+02 # g0mass !not used... Calculated in couplings.f # PDG Width DECAY 6000045 1.00000000e+00 # S0WIDTH !not used... Calculated in couplings.f DECAY 6000046 1.00000000e+00 # O0WIDTH !not used... Calculated in couplings.f DECAY 6000047 1.00000000e+00 # S1WIDTH DECAY 6000048 0.34425000e+03 # O1WIDTH DECAY 6000049 1.00000000e+00 # S2WIDTH !result is independent of this value DECAY 6000050 1.00000000e+00 # g1width !not used... Calculated in couplings.f DECAY 6000051 1.00000000e+00 # g2width !not used... Calculated in couplings.f DECAY 6000052 1.00000000e+00 # g3width !not used... Calculated in couplings.f DECAY 6000053 1.00000000e+00 # g4width !not used... Calculated in couplings.f DECAY 6000054 1.00000000e+00 # g5width !not used... Calculated in couplings.f DECAY 6000055 1.00000000e+00 # g6width !not used... Calculated in couplings.f DECAY 6000056 1.00000000e+00 # g7width !not used... Calculated in couplings.f DECAY 6000057 1.00000000e+00 # g8width !not used... Calculated in couplings.f DECAY 6000058 1.00000000e+00 # g9width !not used... Calculated in couplings.f DECAY 6000059 1.00000000e+00 # g0width !not used... Calculated in couplings.f DECAY 6 1.38104231E+00 # top width DECAY 23 2.44639985E+00 # Z width DECAY 24 2.03535570E+00 # W width DECAY 25 4.27608700E-03 # H width# # BR NDA ID1 ID2 7.18385415E-02 2 4 -4 # BR( H -> c cbar ) 8.03404251E-01 2 5 -5 # BR( H -> b bbar ) 0.00000000E+00 2 6 -6 # BR( H -> t tbar ) 4.83591503E-02 2 15 -15 # BR( H -> tau- tau+) 2.69892403E-05 2 23 23 # BR( H -> Z Z^(*)) 3.67605190E-03 2 24 -24 # BR( H -> W W^(*)) 2.43358656E-02 2 21 21 # BR( H -> g g ) BLOCK MGUSER 1 1.00000000e+00 # s0scalarf ,spin-0 singlet scalar mult.fac. 2 0.00000000e+00 # s0axialf ,spin-0 singlet pseudo-scalar mult.fac. 3 1.00000000e+00 # o0scalarf ,spin-0 octet scalar mult.fac. 4 0.00000000e+00 # o0axialf ,spin-0 octet pseudo-scalar mult.fac. 5 0.20000000e+00 # s1quleft ,spin-1 singlet left up quark mult.fac. 6 0.20000000e+00 # s1quright ,spin-1 singlet right up quark mult.fac. 7 0.20000000e+00 # s1qdleft ,spin-1 singlet left down quark mult.fac. 8 0.20000000e+00 # s1qdright ,spin-1 singlet right down quark mult.fac. 9 1.00000000e+00 # s1qcleft ,spin-1 singlet left charm quark mult.fac. 10 1.00000000e+00 # s1qcright ,spin-1 singlet right charm quark mult.fac. 11 1.00000000e+00 # s1qsleft ,spin-1 singlet left strange quark mult.fac. 12 1.00000000e+00 # s1qsright ,spin-1 singlet right strange quark mult.fac. 13 1.00000000e+00 # s1qtleft ,spin-1 singlet left top quark mult.fac. 14 1.00000000e+00 # s1qtright ,spin-1 singlet right top quark mult.fac. 15 1.00000000e+00 # s1qbleft ,spin-1 singlet left bottom quark mult.fac. 16 1.00000000e+00 # s1qbright ,spin-1 singlet right bottom quark mult.fac. 17 1.00000000e+00 # s1eleft ,spin-1 singlet left electron mult.fac. 18 1.00000000e+00 # s1eright ,spin-1 singlet right electron mult.fac. 19 1.00000000e+00 # s1muleft ,spin-1 singlet left muon mult.fac. 20 1.00000000e+00 # s1muright ,spin-1 singlet right muon mult.fac. 21 1.00000000e+00 # s1taleft ,spin-1 singlet left tau mult.fac. 22 1.00000000e+00 # s1taright ,spin-1 singlet right tau mult.fac. 23 1.00000000e+00 # s1ne ,spin-1 singlet electron-neutrino mult.fac. 24 1.00000000e+00 # s1nm ,spin-1 singlet muon-neutrino mult.fac. 25 1.00000000e+00 # s1nt ,spin-1 singlet tau-neutrino mult.fac. 26 0.20000000e+00 # o1quleft ,spin-1 octet left up quark mult.fac. 27 0.20000000e+00 # o1quright ,spin-1 octet right up quark mult.fac. 28 0.20000000e+00 # o1qdleft ,spin-1 octet left down quark mult.fac. 29 0.20000000e+00 # o1qdright ,spin-1 octet right down quark mult.fac. 30 0.20000000e+00 # o1qcleft ,spin-1 octet left charm quark mult.fac. 31 0.20000000e+00 # o1qcright ,spin-1 octet right charm quark mult.fac. 32 0.20000000e+00 # o1qsleft ,spin-1 octet left strange quark mult.fac. 33 0.20000000e+00 # o1qsright ,spin-1 octet right strange quark mult.fac. 34 1.00000000e+00 # o1qtleft ,spin-1 octet left top quark mult.fac. 35 4.00000000e+00 # o1qtright ,spin-1 octet right top quark mult.fac. 36 1.00000000e+00 # o1qbleft ,spin-1 octet left bottom quark mult.fac. 37 0.20000000e+00 # o1qbright ,spin-1 octet right bottom quark mult.fac. 38 0.07000000e+00 # kapMpl ,spin-2 RS kappa/M_planck (reduced plack mass) 39 3.00000000e+00 # addn ,spin-2 ADD # extra dim. (only addn=3 for now) 40 9.00000000e+02 # mstring ,spin-2 ADD cut-off scale #********************************************************************* # MadGraph/MadEvent * # http://madgraph.hep.uiuc.edu * # * # run_card.dat * # * # 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 * #********************************************************************* # #******************* # Running parameters #******************* # #********************************************************************* # Tag name for the run (one word) * #********************************************************************* 'fermi' = run_tag ! name of the run #********************************************************************* # Run to generate the grid pack * #********************************************************************* .false. = gridpack !True = setting up the grid pack #********************************************************************* # Number of events and rnd seed * #********************************************************************* 10000 = nevents ! Number of unweighted events requested 1 = 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) 3500 = ebeam1 ! beam 1 energy in GeV 3500 = ebeam2 ! beam 2 energy in GeV #********************************************************************* # Beam polarization from -100 (left-handed) to 100 (right-handed) * #********************************************************************* 0 = polbeam1 ! beam polarization for beam 1 0 = polbeam2 ! beam polarization for beam 2 #********************************************************************* # PDF CHOICE: this automatically fixes also alpha_s and its evol. * #********************************************************************* 'cteq6_l1' = pdlabel ! PDF set #********************************************************************* # Renormalization and factorization scales * #********************************************************************* T = fixed_ren_scale ! if .true. use fixed ren scale T = fixed_fac_scale ! if .true. use fixed fac scale 2250. = scale ! fixed ren scale 2250. = dsqrt_q2fact1 ! fixed fact scale for pdf1 2250. = dsqrt_q2fact2 ! fixed fact scale for pdf2 1 = scalefact ! scale factor for event-by-event scales #********************************************************************* # Matching - Warning! ickkw > 0 is still beta #********************************************************************* 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching 1 = highestmult ! for ickkw=2, highest mult group 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE 1 = alpsfact ! scale factor for QCD emission vx F = chcluster ! cluster only according to channel diag F = pdfwgt ! for ickkw=1, perform pdf reweighting #********************************************************************* # #********************************** # BW cutoff (M+/-bwcutoff*Gamma) #********************************** 15 = bwcutoff #******************* # Standard Cuts #******************* # #********************************************************************* # Minimum and maximum pt's * #********************************************************************* 20 = ptj ! minimum pt for the jets 0 = ptb ! minimum pt for the b 10 = pta ! minimum pt for the photons 10 = ptl ! minimum pt for the charged leptons 0 = misset ! minimum missing Et (sum of neutrino's momenta) 0 = ptheavy ! minimum pt for one heavy final state 1.0 = ptonium ! minimum pt for the quarkonium states 1d5 = ptjmax ! maximum pt for the jets 1d5 = ptbmax ! maximum pt for the b 1d5 = ptamax ! maximum pt for the photons 1d5 = ptlmax ! maximum pt for the charged leptons 1d5 = missetmax ! maximum missing Et (sum of neutrino's momenta) #********************************************************************* # Minimum and maximum E's (in the lab frame) * #********************************************************************* 0 = ej ! minimum E for the jets 0 = eb ! minimum E for the b 0 = ea ! minimum E for the photons 0 = el ! minimum E for the charged leptons 1d5 = ejmax ! maximum E for the jets 1d5 = ebmax ! maximum E for the b 1d5 = eamax ! maximum E for the photons 1d5 = elmax ! maximum E for the charged leptons #********************************************************************* # Maximum and minimum rapidity * #********************************************************************* 1d2 = etaj ! max rap for the jets 1d2 = etab ! max rap for the b 1d2 = etaa ! max rap for the photons 1d2 = etal ! max rap for the charged leptons 1d2 = etaonium ! max rap for the quarkonium states 0d0 = etajmin ! min rap for the jets 0d0 = etabmin ! min rap for the b 0d0 = etaamin ! min rap for the photons 0d0 = etalmin ! main rap for the charged leptons #********************************************************************* # Minimum and maximum DeltaR distance * #********************************************************************* 0.0 = drjj ! min distance between jets 0 = drbb ! min distance between b's 0.0 = drll ! min distance between leptons 0.0 = draa ! min distance between gammas 0 = drbj ! min distance between b and jet 0.0 = draj ! min distance between gamma and jet 0.0 = drjl ! min distance between jet and lepton 0 = drab ! min distance between gamma and b 0 = drbl ! min distance between b and lepton 0.0 = dral ! min distance between gamma and lepton 1d2 = drjjmax ! max distance between jets 1d2 = drbbmax ! max distance between b's 1d2 = drllmax ! max distance between leptons 1d2 = draamax ! max distance between gammas 1d2 = drbjmax ! max distance between b and jet 1d2 = drajmax ! max distance between gamma and jet 1d2 = drjlmax ! max distance between jet and lepton 1d2 = drabmax ! max distance between gamma and b 1d2 = drblmax ! max distance between b and lepton 1d2 = dralmax ! maxdistance between gamma and lepton #********************************************************************* # Minimum and maximum invariant mass for pairs * #********************************************************************* 0 = mmjj ! min invariant mass of a jet pair 0 = mmbb ! min invariant mass of a b pair 0 = mmaa ! min invariant mass of gamma gamma pair 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair 1d5 = mmjjmax ! max invariant mass of a jet pair 1d5 = mmbbmax ! max invariant mass of a b pair 1d5 = mmaamax ! max invariant mass of gamma gamma pair 1d5 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair #********************************************************************* # Minimum and maximum invariant mass for all letpons * #********************************************************************* 0 = mmnl ! min invariant mass for all letpons (l+- and vl) 1d5 = mmnlmax ! max invariant mass for all letpons (l+- and vl) #********************************************************************* # Minimum and maximum pt for 4-momenta sum of all letpons * #********************************************************************* 0 = ptllmin ! min invariant mass for all letpons (l+- and vl) 1d5 = ptllmax ! max invariant mass for all letpons (l+- and vl) #********************************************************************* # Inclusive cuts * #********************************************************************* 0 = xptj ! minimum pt for at least one jet 0 = xptb ! minimum pt for at least one b 0 = xpta ! minimum pt for at least one photon 0 = xptl ! minimum pt for at least one charged lepton #********************************************************************* # Control the pt's of the jets sorted by pt * #********************************************************************* 0 = ptj1min ! minimum pt for the leading jet in pt 0 = ptj2min ! minimum pt for the second jet in pt 0 = ptj3min ! minimum pt for the third jet in pt 0 = ptj4min ! minimum pt for the fourth jet in pt 1d5 = ptj1max ! maximum pt for the leading jet in pt 1d5 = ptj2max ! maximum pt for the second jet in pt 1d5 = ptj3max ! maximum pt for the third jet in pt 1d5 = ptj4max ! maximum pt for the fourth jet in pt 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts #********************************************************************* # Control the Ht(k)=Sum of k leading jets * #********************************************************************* 0 = htjmin ! minimum jet HT=Sum(jet pt) 1d5 = htjmax ! maximum jet HT=Sum(jet pt) 0 = ht2min ! minimum Ht for the two leading jets 0 = ht3min ! minimum Ht for the three leading jets 0 = ht4min ! minimum Ht for the four leading jets 1d5 = ht2max ! maximum Ht for the two leading jets 1d5 = ht3max ! maximum Ht for the three leading jets 1d5 = ht4max ! maximum Ht for the four leading jets #********************************************************************* # WBF cuts * #********************************************************************* 0 = xetamin ! minimum rapidity for two jets in the WBF case 0 = deltaeta ! minimum rapidity for two jets in the WBF case #********************************************************************* # maximal pdg code for quark to be considered as a jet * # otherwise b cuts are applied * #********************************************************************* 4 = maxjetflavor #********************************************************************* # Jet measure cuts * #********************************************************************* 0 = xqcut ! minimum kt jet measure between partons #********************************************************************* # Number of Events : 10000 # Integrated weight (pb) : .31173E-01 # Truncated wgt (pb) : .27412E-04 # Unit wgt : .31173E-05