# based on example from Herwig++ 2.4.2 distribution: # share/Herwig++/TVT.in # Run options: cd /Herwig/Generators set LHCGenerator:NumberOfEvents 100000 set LHCGenerator:RandomNumberGenerator:Seed 0 set LHCGenerator:DebugLevel 0 set LHCGenerator:PrintEvent 1 set LHCGenerator:MaxErrors 1000000 # redirect all log output to stdout set LHCGenerator:UseStdout true # do output to a HepMC file cd /Herwig/Generators insert LHCGenerator: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.nvCEU9dwCI/generator.hepmc # set /Herwig/Analysis/HepMCFile:Units GeV_mm # Beam parameters: set LHCGenerator:EventHandler:LuminosityFunction:Energy 900 set LHCGenerator:EventHandler:BeamA /Herwig/Particles/p+ set LHCGenerator:EventHandler:BeamB /Herwig/Particles/p+ set LHCGenerator:MaxErrors -1 # MinBias process: # reference: # Herwig++ Physics and Manual (http://arxiv.org/abs/0803.0883) # chapter B.2 Examples / Setup for minimum bias runs cd /Herwig/MatrixElements insert SimpleQCD:MatrixElements[0] MEMinBias cd /Herwig/Cuts set JetKtCut:MinKT 0.0*GeV set QCDCuts:MHatMin 0.0*GeV set QCDCuts:X1Min 0.01 set QCDCuts:X2Min 0.01 set /Herwig/UnderlyingEvent/MPIHandler:softInt Yes set /Herwig/UnderlyingEvent/MPIHandler:twoComp Yes # 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 900*GeV # maximum jet pT # Use this for hard leading-jets in a certain mHat window #set /Herwig/Cuts/QCDCuts:MHatMin 0*GeV # minimum jet mHat #set /Herwig/Cuts/QCDCuts:MHatMax 900*GeV # maximum jet mHat # 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% # Default tune (no additional parameters) # --------------------------------------------- set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0 # Run generator cd /Herwig/Generators #%xsecfix% # Get the cross section right # (as advised e.g. in https://herwig.hepforge.org/trac/browser/Tests/Rivet/SppS-900-UE.in) create Herwig::MPIXSecReweighter MPIXSecReweighter insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter run TVT LHCGenerator