# based on example from Herwig++ 2.4.2 distribution:
#   share/Herwig++/TVT.in

# Run options:
cd /Herwig/Generators
set LHCGenerator:NumberOfEvents 1000
set LHCGenerator:RandomNumberGenerator:Seed 1
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 /tmp/akorneev/tmp.MxWuuj2928/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+


# 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
run TVT LHCGenerator