# 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 0 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.U1tuBubDVd/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 ## W+jet #cd /Herwig/MatrixElements #insert SubProcess:MatrixElements[0] MEWJet # W+jet cd /Herwig/MatrixElements insert SubProcess:MatrixElements[0] MEqq2W2ff # Set cuts # Use this for lepton pair invariant mass in a certain window set /Herwig/Cuts/MassCut:MinM 10*GeV # min inv. mass of lep. pair set /Herwig/Cuts/MassCut:MaxM 7000*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