jobOptions_FlukaBeamHalo.py

Go to the documentation of this file.

import AthenaCommon.AtlasUnixGeneratorJob
 
 
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
 
 
from AthenaCommon.AppMgr import ServiceMgr
 
 
from AthenaCommon.AppMgr import theApp
 
#theApp.EvtMax = 2835459 # All of the events
theApp.EvtMax = 3000000   
 
 
 
ServiceMgr.MessageSvc.OutputLevel = INFO
 
from BeamHaloGenerator.BeamHaloGeneratorConf import BeamHaloGeneratorAlg
BeamHaloGenerator = BeamHaloGeneratorAlg()
BeamHaloGenerator.inputType="FLUKA-RB"
BeamHaloGenerator.inputFile="/pcdisk/pcgeneve01/wbell/data/beamhalo/beamgas/beam-gas_IR1_1GeV_cutoff_D1_field"
BeamHaloGenerator.interfacePlane = -22600.0
BeamHaloGenerator.flipProbability = 0.0 
BeamHaloGenerator.randomStream = "BeamHalo"
BeamHaloGenerator.doMonitoringPlots = True
 
# The generator settings determine if the event is accepted.
#   * If the allowedPdgId list is not given all particles are accepted.
#   * Limits are in the form of (lower limit, upper limit)
#   * If a limit is not given it is disabled.
#   * If a limit value is -1 then it is disabled.
#   * All limits are checked against |value|
#   * r = sqrt(x^2 + y^2)
BeamHaloGenerator.generatorSettings = []#
 
 
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#*' ]