GaussGen.RunNumber = 1082
LHCbApp().EvtMax = 10000
#LHCbApp().setProp("OutputLevel", 3)
#LHCbApp().setProp("OutputLevel", 4)
LHCbApp().OutputLevel = WARNING
Gauss().Histograms = "NONE"
Gauss().OutputType = "NONE"
Gauss().DatasetName = "GaussDevWHepMC"
#Special settings for beam conditions:
#** ignore multiple pp collisions, force nPV===1 (algorithm will complain if not)
from Configurables import Generation
gaussGen = Generation("Generation")
gaussGen.PileUpTool = "FixedNInteractions"
#** force head-on pp collisions:
# - theoreticians always assume pp collision is in CM frame
# - you may experience some problems in generating event with settings below (issue under investigation)
# - see CorrectCrossingAngles option below to control CM boost at runtime
#Gauss().BeamCrossingAngle = 0.0
#Gauss().BeamBetaStar = 0.0
GenMonitor = GaudiSequencer("GenMonitor")
GenMonitor.Members += ["RivetAnalysisHandler"]
rivet = RivetAnalysisHandler()
rivet.BaseFileName = "testBird"
rivet.RunName = "LHCbMBPy8"
rivet.Analyses = ["MC_GENERIC", "MC_IDENTIFIED"]
#to search Rivet plugins in current directory (when developing/LbRivetPlugins)
# Generate Pile Up
from Gaudi.Configuration import *
from Configurables import Generation, Pythia8Production, Gauss
gen = Generation()
gen.Special.PileUpProductionTool = "Pythia8Production/Pythia8PileUp"
gen.PileUpTool = "FixedLuminosityForRareProcess"
gen.Special.addTool(Pythia8Production, name="Pythia8PileUp")
gen.Special.Pythia8PileUp.Tuning = "LHCbDefault.cmd"
gen.Special.ReinitializePileUpGenerator = False
# AlpGen specific options
from Configurables import AlpGenProductionwqq, AlpGenProductionzbb, AlpGenProductiontop, AlpGenProduction4q
# AlpGen conf dict
try:
AlpGenDict = Generation.AlpGenDict
except:
raise NameError("AlpGenDict not defined")
# File type
try:
alpgen_FileLabel = AlpGenDict["alpgen_FileLabel"]
except:
raise NameError('value for alpgen_FileLabel property not set')
else:
if alpgen_FileLabel == "wbb" or alpgen_FileLabel == "wcc":
gen.Special.ProductionTool = "AlpGenProductionwqq"
gen.Special.addTool(AlpGenProductionwqq)
alpgen = gen.Special.AlpGenProductionwqq