# based on example from Herwig++ 2.4.2 distribution: # share/Herwig++/TVT.in # rungen.sh uses the following command to insert PPCollider.in snippet #%beamsetup71% read snippets/PPCollider.in # Run options: cd /Herwig/Generators set EventGenerator:NumberOfEvents 100000 set EventGenerator:RandomNumberGenerator:Seed 0 set EventGenerator:DebugLevel 9 set EventGenerator:PrintEvent 1 set EventGenerator:MaxErrors 100000 # redirect all log output to stdout set EventGenerator:UseStdout true # do output to a HepMC file cd /Herwig/Generators insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile set /Herwig/Analysis/HepMCFile:PrintEvent 1000000 set /Herwig/Analysis/HepMCFile:Format GenEvent set /Herwig/Analysis/HepMCFile:Filename /shared/tmp/tmp.CJKTG839lg/generator.hepmc # set /Herwig/Analysis/HepMCFile:Units GeV_mm # Beam parameters: set EventGenerator:EventHandler:LuminosityFunction:Energy 7000 set EventGenerator:EventHandler:BeamA /Herwig/Particles/p+ set EventGenerator:EventHandler:BeamB /Herwig/Particles/p+ set EventGenerator:MaxErrors -1 # Drell-Yan Z/gamma process: cd /Herwig/MatrixElements insert SubProcess:MatrixElements[0] MEqq2gZ2ff # Select only electron channel #set MEqq2gZ2ff:Process Electron # Set cuts # Use this for hard leading-jets in a certain pT window set /Herwig/Cuts/JetKtCut:MinKT 0*GeV # minimum jet pT set /Herwig/Cuts/JetKtCut:MaxKT 7000*GeV # maximum jet pT # Use this for lepton pair invariant mass in a certain window set /Herwig/Cuts/MassCut:MinM 50*GeV # min inv. mass of lep. pair set /Herwig/Cuts/MassCut:MaxM 130*GeV # max inv. mass of lep. pair # Make particles with c*tau > 10 mm stable: set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm set /Herwig/Decays/DecayHandler:LifeTimeOption Average # tune 'default' parameters: ------------------- #%tuneFile% # --------------------------------------------- set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE -1 # Run generator cd /Herwig/Generators run TVT EventGenerator