def run(input_file):
import DDG4
from DDG4 import OutputLevel as Output
kernel = DDG4.Kernel()
lcdd = kernel.lcdd()
gen = DDG4.GeneratorAction(kernel, "Geant4InputAction/Input")
kernel.generatorAction().adopt(gen)
gen.Input = "Geant4EventReaderHepMC|/home/frankm/SW/data/hepmc_geant4.dat"
gen.Input = "Geant4EventReaderHepMC|/home/frankm/SW/data/"
gen.Input = "Geant4EventReaderHepMC|" + input_file
gen.OutputLevel = Output.DEBUG
gen.HaveAbort = False
prim_vtx = DDG4.Geant4Vertex()
prim_vtx.x = 0.0
prim_vtx.y = 0.0
prim_vtx.z = 0.0
parts = gen.new_particles()
ret = 1
while ret:
try:
ret = gen.readParticles(0, prim_vtx, parts)
except Exception, X:
print '\nException: readParticles:', str(X)
ret = None
if ret:
for p in parts:
print 'ID:%5d PDG-id:%8d Charge:%1d Mass:%8.3g Momentum:(%8.2g,%8.2g,%8.2g) '\
'Vertex:(%8.2g,%8.2g,%8.2g) NDau:%2d Status:%08X'%\
(p.id,p.pdgID,int(p.charge3()),p.mass,p.psx,p.psy,p.psz,\
p.vsx,p.vsy,p.vsz,len(p.daughters),p.status)
parts.clear()
else:
print '*** End of recordset or read failure.....'
print 132 * '*'
