# based on example from Herwig 7.0.4 distribution: # src/LHC-mb.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.nLEMX9aO7U/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: cd /Herwig/MatrixElements/ insert SimpleQCD:MatrixElements[0] MEMinBias # MPI model settings set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0 cd /Herwig/Generators # Change to have no pT cuts for MinBias set LHCGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts # Make particles with c*tau > 10 mm stable: set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm set /Herwig/Decays/DecayHandler:LifeTimeOption Average # tune 'UE-MMHT' parameters: ------------------- #%tuneFile% # 9 Jan 2018: adapted/copied from # https://herwig.hepforge.org/tutorials/mpi/tunes.html#herwig-7-1-tunes # by P. Skands # MPI energy extrapolation settings set /Herwig/UnderlyingEvent/MPIHandler:EnergyExtrapolation Power set /Herwig/UnderlyingEvent/MPIHandler:ReferenceScale 7000.*GeV set /Herwig/UnderlyingEvent/MPIHandler:Power 0.366 set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 4.39*GeV # Energy-independent MPI parameters # Colour reconnection settings set /Herwig/Hadronization/ColourReconnector:ColourReconnection Yes set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.4276 # Colour Disrupt settings set /Herwig/Partons/RemnantDecayer:colourDisrupt 0.798 # inverse hadron radius set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 2.30 # --------------------------------------------- # Run generator cd /Herwig/Generators # Fix to get correct cross section create Herwig::MPIXSecReweighter MPIXSecReweighter insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter run TVT LHCGenerator