[madgraph5amc] # suffix keyword '[QCD]' means activate NLO QCD corrections generate p p > l+ l- [QCD] add process p p > l+ l- j [QCD] display processes display particles display multiparticles # generate .ps with digrams pictures #display diagrams # '-f' to force rewrite of existing MG5RUN/ subdirectory output MG5RUN -f # '-p' to Stop the run after the parton level file generation launch MG5RUN -p # set which parton-shower code will be used on LesHouches events after MG5 run set parton_shower PYTHIA8 set nevents 100000 set iseed 4 set lpp1 1 set lpp2 1 set ebeam1 3500. set ebeam2 3500. # Min inv. mass of all opp. sign same-flavor lepton pairs set mll_sf 50 # Min jet transverse momentum (run_card.dat, default = 10) #set ptj 10 # FxFx merging set ickkw 3 # FxFx merging scale (shower_card.dat, default = -1), limit Qcut >= ptj set Qcut 20 [pythia8] Main:numberOfEvents = 100000 ! number of events to generate Next:numberShowEvent = 0 ! suppress full listing of first events Init:showChangedParticleData = off ! suppress changes particle data (a lot of output with LHE) # random seed Random:setSeed = on Random:seed = 4 # pythia8 showering # http://home.thep.lu.se/~torbjorn/pythia82html/LesHouchesAccord.html Beams:frameType = 4 Beams:LHEF = /shared/tmp/tmp.5pe9y8ma4o/MG5RUN/Events/run_01/events.lhe # pythia8 merging steering (FxFx) as generated with MG5_aMC_v2_6_4/MG5aMC_PY8_interface BeamRemnants:primordialKT = off Check:epTolErr = 0.0010000 JetMatching:coneRadius = 1.00000 JetMatching:doFxFx = on JetMatching:etaJetMax = 1000.000 JetMatching:merge = on JetMatching:nJetMax = 0 JetMatching:qCut = 20.00000 JetMatching:setMad = off PartonLevel:MPI = off SpaceShower:alphaEMorder = 0 SpaceShower:alphaSvalue = 0.11800 SpaceShower:MEcorrections = off SpaceShower:pTmaxMatch = 1 SpaceShower:QEDshowerByL = off SpaceShower:QEDshowerByQ = off SpaceShower:rapidityOrder = off TimeShower:alphaEMorder = 0 TimeShower:alphaSvalue = 0.11800 TimeShower:dampenBeamRecoil = off TimeShower:globalRecoil = on TimeShower:globalRecoilMode = 2 TimeShower:limitPTmaxGlobal = on TimeShower:MEcorrections = off TimeShower:nMaxGlobalBranch = 1 TimeShower:nMaxGlobalRecoil = 1 TimeShower:QEDshowerByGamma = off TimeShower:QEDshowerByL = off TimeShower:QEDshowerByQ = off TimeShower:weightGluonToQuark = 1