#5.1.5.11 #************************************************************ #* MadGraph 5 * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 1.5.11 2013-06-21 * #* * #* The MadGraph Development Team - Please visit us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* * #************************************************************ #* * #* Command File for MadGraph 5 * #* * #* run as ./bin/mg5 filename * #* * #************************************************************ set group_subprocesses Auto set ignore_six_quark_processes False set gauge unitary set complex_mass_scheme False import model 2HDM_GF define p = g u c d s u~ c~ d~ s~ define j = g u c d s u~ c~ d~ s~ define p = g u c d s u~ c~ d~ s~ b b~ define j = g u c d s u~ c~ d~ s~ b b~ a define l+ = e+ mu+ define l- = e- mu- define vl = ve vm vt define vl~ = ve~ vm~ vt~ # Define multiparticle labels # Specify process(es) to run generate g g > h3, (h3 > z h1, (z > l+ l-, h1 > ta- ta+ )) # Output processes to MadEvent directory output -f #********************************************************************* # MadGraph/MadEvent * # http://madgraph.hep.uiuc.edu * # * # proc_card.dat * #********************************************************************* # * # This Files is generated by MADGRAPH 5 * # * # WARNING: This Files is generated for MADEVENT (compatibility issue)* # This files is NOT a valid MG4 proc_card.dat * # Running this in MG4 will NEVER reproduce the result of MG5* # * #********************************************************************* #********************************************************************* # Process(es) requested : mg2 input * #********************************************************************* # Begin PROCESS # This is TAG. Do not modify this line g g > h3 , (h3 > z h1 , (z > l+ l- , h1 > ta- ta+ )) #Process # Be carefull the coupling are here in MG5 convention 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 2HDM_GF # End MODEL # This is TAG. Do not modify this line #********************************************************************* # Start multiparticle definitions * #********************************************************************* # Begin MULTIPARTICLES # This is TAG. Do not modify this line # End MULTIPARTICLES # This is TAG. Do not modify this line #********************************************************************* # 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) * #********************************************************************* tag_1 = 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 * # 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. * #********************************************************************* 5000 = nevents ! Number of unweighted events requested 48 = 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 4000 = ebeam1 ! beam 1 total energy in GeV 4000 = ebeam2 ! beam 2 total 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. * #********************************************************************* 'cteq6l1' = pdlabel ! PDF set #********************************************************************* # Renormalization and factorization scales * #********************************************************************* F = fixed_ren_scale ! if .true. use fixed ren scale F = fixed_fac_scale ! if .true. use fixed fac scale 91.1880 = scale ! fixed ren scale 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2 1 = scalefact ! scale factor for event-by-event scales #********************************************************************* # Matching - Warning! ickkw > 1 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 T = pdfwgt ! for ickkw=1, perform pdf reweighting 5 = asrwgtflavor ! highest quark flavor for a_s reweight #********************************************************************* # Automatic ptj and mjj cuts if xqcut > 0 # (turn off for VBF and single top processes) #********************************************************** T = auto_ptj_mjj ! Automatic setting of ptj and mjj #********************************************************** # #********************************** # BW cutoff (M+/-bwcutoff*Gamma) #********************************** 15 = bwcutoff ! (M+/-bwcutoff*Gamma) #********************************************************** # Apply pt/E/eta/dr/mij cuts on decay products or not # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) #********************************************************** T = cut_decays ! Cut decay products #************************************************************* # Number of helicities to sum per event (0 = all helicities) # 0 gives more stable result, but longer run time (needed for # long decay chains e.g.). # Use >=2 if most helicities contribute, e.g. pure QCD. #************************************************************* 0 = nhel ! Number of helicities used per event #******************* # Standard Cuts #******************* # #********************************************************************* # Minimum and maximum pt's (for max, -1 means no cut) * #********************************************************************* 1 = ptj ! minimum pt for the jets 1 = ptb ! minimum pt for the b 0 = pta ! minimum pt for the photons 0 = 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 0 = ptonium ! minimum pt for the quarkonium states -1 = ptjmax ! maximum pt for the jets -1 = ptbmax ! maximum pt for the b -1 = ptamax ! maximum pt for the photons -1 = ptlmax ! maximum pt for the charged leptons -1 = 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 -1 = ejmax ! maximum E for the jets -1 = ebmax ! maximum E for the b -1 = eamax ! maximum E for the photons -1 = elmax ! maximum E for the charged leptons #********************************************************************* # Maximum and minimum absolute rapidity (for max, -1 means no cut) * #********************************************************************* 30 = etaj ! max rap for the jets 30 = etab ! max rap for the b 30 = etaa ! max rap for the photons 30 = etal ! max rap for the charged leptons 30 = etaonium ! max rap for the quarkonium states 0 = etajmin ! min rap for the jets 0 = etabmin ! min rap for the b 0 = etaamin ! min rap for the photons 0 = etalmin ! main rap for the charged leptons #********************************************************************* # Minimum and maximum DeltaR distance * #********************************************************************* 0 = drjj ! min distance between jets 0 = drbb ! min distance between b's 0 = drll ! min distance between leptons 0 = draa ! min distance between gammas 0 = drbj ! min distance between b and jet 0 = draj ! min distance between gamma and jet 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 = dral ! min distance between gamma and lepton -1 = drjjmax ! max distance between jets -1 = drbbmax ! max distance between b's -1 = drllmax ! max distance between leptons -1 = draamax ! max distance between gammas -1 = drbjmax ! max distance between b and jet -1 = drajmax ! max distance between gamma and jet -1 = drjlmax ! max distance between jet and lepton -1 = drabmax ! max distance between gamma and b -1 = drblmax ! max distance between b and lepton -1 = 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 -1 = mmjjmax ! max invariant mass of a jet pair -1 = mmbbmax ! max invariant mass of a b pair -1 = mmaamax ! max invariant mass of gamma gamma pair -1 = 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) -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl) #********************************************************************* # Minimum and maximum pt for 4-momenta sum of leptons * #********************************************************************* 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(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 -1 = ptj1max ! maximum pt for the leading jet in pt -1 = ptj2max ! maximum pt for the second jet in pt -1 = ptj3max ! maximum pt for the third jet in pt -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt 0 = ptl2min ! minimum pt for the second lepton in pt 0 = ptl3min ! minimum pt for the third lepton in pt 0 = ptl4min ! minimum pt for the fourth lepton in pt -1 = ptl1max ! maximum pt for the leading lepton in pt -1 = ptl2max ! maximum pt for the second lepton in pt -1 = ptl3max ! maximum pt for the third lepton in pt -1 = ptl4max ! maximum pt for the fourth lepton in pt #********************************************************************* # Control the Ht(k)=Sum of k leading jets * #********************************************************************* 0 = htjmin ! minimum jet HT=Sum(jet pt) -1 = htjmax ! maximum jet HT=Sum(jet pt) 0 = ihtmin !inclusive Ht for all partons (including b) -1 = ihtmax !inclusive Ht for all partons (including b) 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 -1 = ht2max ! maximum Ht for the two leading jets -1 = ht3max ! maximum Ht for the three leading jets -1 = 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 light jet * # (otherwise b cuts are applied) * #********************************************************************* 5 = maxjetflavor ! Maximum jet pdg code #********************************************************************* # Jet measure cuts * #********************************************************************* 10 = xqcut ! minimum kt jet measure between partons #********************************************************************* ###################################################################### ## 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. By default, ## ## this is only 1->2 decay modes. ## ## ## ###################################################################### ################################### ## INFORMATION FOR CKMI ################################### Block ckmi 1 1 0.000000e+00 # CKMI1x1 1 2 0.000000e+00 # CKMI1x2 1 3 0.000000e+00 # CKMI1x3 2 1 0.000000e+00 # CKMI2x1 2 2 0.000000e+00 # CKMI2x2 2 3 0.000000e+00 # CKMI2x3 3 1 0.000000e+00 # CKMI3x1 3 2 0.000000e+00 # CKMI3x2 3 3 0.000000e+00 # CKMI3x3 ################################### ## INFORMATION FOR CKMR ################################### Block ckmr 1 1 9.745891e-01 # CKMR1x1 1 2 2.240000e-01 # CKMR1x2 1 3 0.000000e+00 # CKMR1x3 2 1 -2.240000e-01 # CKMR2x1 2 2 9.745891e-01 # CKMR2x2 2 3 0.000000e+00 # CKMR2x3 3 1 0.000000e+00 # CKMR3x1 3 2 0.000000e+00 # CKMR3x2 3 3 1.000000e+00 # CKMR3x3 ################################### ## INFORMATION FOR HIGGS ################################### Block higgs 1 1.000000e+00 # l1 2 1.000000e+00 # l2 3 1.000000e+00 # l3 4 5.000000e-01 # l4 5 4.000000e-01 # lR5 6 3.000000e-01 # lI5 7 0.000000e+00 # lR6 8 2.000000e-01 # lI6 9 0.000000e+00 # lR7 10 0.000000e+00 # lI7 ################################### ## INFORMATION FOR HIGGSMIX ################################### Block higgsmix 1 1 0.000000e+00 # TH1x1 1 2 7.806441e-01 # TH1x2 1 3 6.249758e-01 # TH1x3 2 1 0.000000e+00 # TH2x1 2 2 -6.249758e-01 # TH2x2 2 3 7.806441e-01 # TH2x3 3 1 1.000000e+00 # TH3x1 3 2 0.000000e+00 # TH3x2 3 3 0.000000e+00 # TH3x3 ################################### ## INFORMATION FOR MASS ################################### Block mass 3 1.050000e-01 # MS 4 1.250000e+00 # MC 5 4.200000e+00 # MB 6 1.743000e+02 # MT 13 1.060000e-01 # MM 15 1.777000e+00 # MTA 23 9.118760e+01 # MZ 25 1.250000e+02 # mh1 35 3.000000e+02 # mh2 36 3.000000e+02 # mh3 37 3.000000e+02 # mhc ## 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.000000 # d : 0.0 2 0.000000 # u : 0.0 11 0.000000 # e- : 0.0 12 0.000000 # ve : 0.0 14 0.000000 # vm : 0.0 16 0.000000 # vt : 0.0 21 0.000000 # g : 0.0 22 0.000000 # a : 0.0 24 79.829013 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2))) ################################### ## INFORMATION FOR SMINPUTS ################################### Block sminputs 1 1.279340e+02 # aEWM1 2 1.166370e-05 # Gf 3 1.172000e-01 # aS ################################### ## INFORMATION FOR YUKAWAD ################################### Block yukawad 1 0.000000e+00 # yukd1 2 0.000000e+00 # yukd2 3 3.000000e+00 # yukd3 ################################### ## INFORMATION FOR YUKAWAL ################################### Block yukawal 1 0.000000e+00 # yukl1 2 0.000000e+00 # yukl2 3 1.777000e+00 # yukl3 ################################### ## INFORMATION FOR YUKAWAU ################################### Block yukawau 1 0.000000e+00 # yuku1 2 6.000000e-01 # yuku2 3 1.750000e+02 # yuku3 ################################### ## INFORMATION FOR YUKAWAGDI ################################### Block yukawagdi 1 1 0.000000e+00 # GDI1x1 1 2 0.000000e+00 # GDI1x2 1 3 0.000000e+00 # GDI1x3 2 1 0.000000e+00 # GDI2x1 2 2 0.000000e+00 # GDI2x2 2 3 0.000000e+00 # GDI2x3 3 1 0.000000e+00 # GDI3x1 3 2 0.000000e+00 # GDI3x2 3 3 0.000000e+00 # GDI3x3 ################################### ## INFORMATION FOR YUKAWAGDR ################################### Block yukawagdr 1 1 0.000000e+00 # GDR1x1 1 2 0.000000e+00 # GDR1x2 1 3 0.000000e+00 # GDR1x3 2 1 0.000000e+00 # GDR2x1 2 2 4.000000e-01 # GDR2x2 2 3 0.000000e+00 # GDR2x3 3 1 0.000000e+00 # GDR3x1 3 2 2.000000e-01 # GDR3x2 3 3 5.000000e+00 # GDR3x3 ################################### ## INFORMATION FOR YUKAWAGLI ################################### Block yukawagli 1 1 0.000000e+00 # GLI1x1 1 2 0.000000e+00 # GLI1x2 1 3 0.000000e+00 # GLI1x3 2 1 0.000000e+00 # GLI2x1 2 2 0.000000e+00 # GLI2x2 2 3 0.000000e+00 # GLI2x3 3 1 0.000000e+00 # GLI3x1 3 2 0.000000e+00 # GLI3x2 3 3 0.000000e+00 # GLI3x3 ################################### ## INFORMATION FOR YUKAWAGLR ################################### Block yukawaglr 1 1 0.000000e+00 # GLR1x1 1 2 0.000000e+00 # GLR1x2 1 3 0.000000e+00 # GLR1x3 2 1 0.000000e+00 # GLR2x1 2 2 1.000000e-01 # GLR2x2 2 3 0.000000e+00 # GLR2x3 3 1 0.000000e+00 # GLR3x1 3 2 5.000000e-01 # GLR3x2 3 3 3.000000e+00 # GLR3x3 ################################### ## INFORMATION FOR YUKAWAGUI ################################### Block yukawagui 1 1 0.000000e+00 # GUI1x1 1 2 0.000000e+00 # GUI1x2 1 3 0.000000e+00 # GUI1x3 2 1 0.000000e+00 # GUI2x1 2 2 0.000000e+00 # GUI2x2 2 3 0.000000e+00 # GUI2x3 3 1 0.000000e+00 # GUI3x1 3 2 0.000000e+00 # GUI3x2 3 3 0.000000e+00 # GUI3x3 ################################### ## INFORMATION FOR YUKAWAGUR ################################### Block yukawagur 1 1 0.000000e+00 # GUR1x1 1 2 0.000000e+00 # GUR1x2 1 3 0.000000e+00 # GUR1x3 2 1 0.000000e+00 # GUR2x1 2 2 2.000000e+00 # GUR2x2 2 3 0.000000e+00 # GUR2x3 3 1 0.000000e+00 # GUR3x1 3 2 1.000000e+00 # GUR3x2 3 3 1.000000e+02 # GUR3x3 ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.514871e+00 # WT DECAY 23 2.412628e+00 # WZ DECAY 24 2.002884e+00 # WW DECAY 25 4.070000e-03 # Wh1 DECAY 35 1.000000e+00 # Wh2 DECAY 36 1.000000e+00 # Wh3 DECAY 37 1.000000e+00 # whc ## 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.000000 # d : 0.0 DECAY 2 0.000000 # u : 0.0 DECAY 3 0.000000 # s : 0.0 DECAY 4 0.000000 # c : 0.0 DECAY 5 0.000000 # b : 0.0 DECAY 11 0.000000 # e- : 0.0 DECAY 12 0.000000 # ve : 0.0 DECAY 13 0.000000 # mu- : 0.0 DECAY 14 0.000000 # vm : 0.0 DECAY 15 0.000000 # ta- : 0.0 DECAY 16 0.000000 # vt : 0.0 DECAY 21 0.000000 # g : 0.0 DECAY 22 0.000000 # a : 0.0 #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) #=========================================================== Block QNUMBERS 37 # h+ 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 35 # h2 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 36 # h3 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) # Number of Events : 5000 # Integrated weight (pb) : 2.8527e-12