File indexing completed on 2024-09-27 07:03:48
0001 from Gaudi.Configuration import *
0002 from GaudiKernel import SystemOfUnits as units
0003
0004 from GaudiKernel.DataObjectHandleBase import DataObjectHandleBase
0005 from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
0006
0007 geo_service = GeoSvc("GeoSvc")
0008 podioevent = EICDataSvc("EventDataSvc", inputs=["rich_test.root"], OutputLevel=DEBUG)
0009
0010 from Configurables import PodioInput
0011 from Configurables import Jug__Digi__PhotoMultiplierDigi as PhotoMultiplierDigi
0012
0013 qe_data = [(1.0, 0.25), (7.5, 0.25),]
0014 podioinput = PodioInput("PodioReader", collections=["MCParticles", "ForwardRICHHits"], OutputLevel=DEBUG)
0015 pmtdigi = PhotoMultiplierDigi(inputHitCollection="ForwardRICHHits", outputHitCollection="DigiForwardRICHHits",
0016 quantumEfficiency=[(a*units.eV, b) for a, b in qe_data])
0017
0018 out = PodioOutput("out", filename="digi_rich_test.root")
0019 out.outputCommands = ["keep *"]
0020
0021 ApplicationMgr(
0022 TopAlg = [podioinput, pmtdigi, out],
0023 EvtSel = 'NONE',
0024 EvtMax = 100,
0025 ExtSvc = [podioevent],
0026 OutputLevel=DEBUG
0027 )
0028