# Number of events to be processed (default is 10)
theApp.EvtMax = 5
#--------------------------------------------------------------
# Algorithms Private Options
#--------------------------------------------------------------
from AthenaServices.AthenaServicesConf import AtRndmGenSvc
ServiceMgr += AtRndmGenSvc()
ServiceMgr.AtRndmGenSvc.Seeds = [
"HERWIG 390020611 821000366", "HERWIG_INIT 820021 2347532"
]
from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()
from Herwig_i.Herwig_iConf import Herwig
job += Herwig()
job.Herwig.HerwigCommand = [
"iproc 11500", "modpdf 10042", "autpdf HWLHAPDF", "msflag 1", "jmbug 0",
"ptmin 10."
]
job.Herwig.UpdateHadronStatCode = True
from TruthExamples.TruthExamplesConf import DumpMC
job += DumpMC()
#==============================================================
#
# End of job options file
#
###############################################################
#--------------------------------------------------------------
# Private Application Configuration options
#--------------------------------------------------------------
# Load "user algorithm"
#top algorithms to be run, and the libraries that house them
########
# Read the McEventCollection
#
if not 'MCEVENTKEY' in dir():
MCEVENTKEY = "GEN_EVENT"
pass
if 'DUMP' in dir():
from TruthExamples.TruthExamplesConf import DumpMC
topSequence += DumpMC("ReadGenEvent", McEventKey=MCEVENTKEY)
########
# Dump the McEventCollection into a dumb ASCII file
#
if not 'OUTPUT' in dir():
OUTPUT = "hepmc.ascii"
pass
from McParticleTools.McParticleToolsConf import HepMcWriterTool
from McParticleAlgs.McParticleAlgsConf import GenEventAsciiWriter
topSequence += GenEventAsciiWriter(
McWriter=HepMcWriterTool(Output=OUTPUT, McEvents=MCEVENTKEY))
#--------------------------------------------------------------
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
"zLimits -4000. 4000.", # required for shape=cylinder
# "ptLimits 0. 0.", # lower and upper limits. The upper limit is optional
# "phiLimits 0. 0.", # lower and upper limits. The upper limit is optional
# "etaLimits 0. 0.", # lower and upper limits. The upper limit is optional
"rLimits 650. 1950."
] # required for shape=cylinder
# "weightLimits 0. 0."] # lower and upper limits. The upper limit is optional
topSequence += BeamHaloGenerator
print BeamHaloGenerator
# Required when monitoring plots are enabled.
theApp.HistogramPersistency = "ROOT"
if not hasattr(ServiceMgr, 'THistSvc'):
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output += [
"BeamHalo DATAFILE='BeamHalo.root' OPT='RECREATE'"
]
# For debugging
from TruthExamples.TruthExamplesConf import DumpMC
topSequence += DumpMC()
from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream
Stream1 = AthenaPoolOutputStream("StreamEVGEN")
Stream1.OutputFile = "evgen.pool.root"
Stream1.ItemList += ['EventInfo#*', 'McEventCollection#*']
# configuring the Athena application for a 'generator' job import AthenaCommon.AtlasUnixGeneratorJob # make sure we are loading the ParticleProperty service from PartPropSvc.PartPropSvcConf import PartPropSvc svcMgr += PartPropSvc() #-------------------------------------------------------------- # Private Application Configuration options #-------------------------------------------------------------- from AthenaCommon.AlgSequence import AlgSequence topAlg = AlgSequence() from Hijing_i.Hijing_iConf import Hijing topAlg += Hijing() from TruthExamples.TruthExamplesConf import DumpMC topAlg += DumpMC() theApp.ExtSvc += ["AtRndmGenSvc"] # Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) svcMgr.MessageSvc.OutputLevel = 4 svcMgr.MessageSvc.defaultLimit = 100000 #-------------------------------------------------------------- # Event related parameters #-------------------------------------------------------------- # Number of events to be processed (default is 10) theApp.EvtMax = 5 # Set run number (default 0 causes problems) svcMgr.EventSelector.RunNumber = 12345 #--------------------------------------------------------------