def initialize(self, reinit=False):
Helper.service(theApp._svcloc, 'HistogramDataSvc', True)
from CTB_G4Sim.CtbPhotonESpectrum import Photon_Energy_spectrum
p = Photon_Energy_spectrum(SimFlags.RunNumber())
p.build_efunction(theApp)
h2 = p.build_spectrum(theApp)
return True
# ParticleGenerator with an energy = 180 Gev - photon-energy
SimFlags.ParticlePDG.set_Value('22')
#
# - Option: we reduce the beam spot for all the photon runs
#----------
if SimFlags.ParticleGeneratorOrders().has_key('vertY:'):
SimFlags.ParticleGeneratorOrders()['vertY:'] = ' flat -5. 5.'
if SimFlags.ParticleGeneratorOrders().has_key('vertZ:'):
SimFlags.ParticleGeneratorOrders()['vertZ:'] = ' flat -5. 5.'
# - Option: run with a given energy distribution for the photon :
#----------
# In this case you have to customize the ParticleGenerator (see
# in the generator section)
SimFlags.ParticleGeneratorOrders()['energy[0]:']=' histogram '+\
repr(SimFlags.RunNumber())
SimFlags.ParticleGeneratorOrders()['pdgcode[0]:'] = ' 22'
SimFlags.ParticleGeneratorOrders(
)['energy[1]:'] = ' function espectrum energy[0]'
SimFlags.ParticleGeneratorOrders()['pdgcode[1]:'] = ' 11'
#
#--- Output printout level -----------------------------------
#output threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL)
MessageSvc = Service("MessageSvc")
MessageSvc.OutputLevel = 4
#==============================================================
#==============================================================
# Job configuration
# ***>> Do not add flags or simulation options below this line