#!/usr/bin/env python from Gaudi.Configuration import * from Configurables import K4DataSvc dsvc = K4DataSvc("EventDataSvc") from Configurables import PlcioWriteAlg alg = PlcioWriteAlg("PlcioWriteAlg") alg.HeaderCol.Path = "EventHeader" alg.OutputCol.Path = "MCParticle" from Configurables import PodioOutput out = PodioOutput("out") out.filename = "test.root" out.outputCommands = ["keep *"] # ApplicationMgr from Configurables import ApplicationMgr ApplicationMgr(TopAlg=[alg, out], EvtSel='NONE', EvtMax=10, ExtSvc=[dsvc], OutputLevel=DEBUG)
from Gaudi.Configuration import * from Configurables import K4DataSvc podioevent = K4DataSvc("EventDataSvc") podioevent.input = "output_k4test_exampledata.root" from Configurables import CreateExampleEventData producer = CreateExampleEventData() from Configurables import PodioInput inp = PodioInput("InputReader") inp.collections = ["MCParticles", "SimTrackerHit"] from Configurables import ApplicationMgr ApplicationMgr( TopAlg=[inp], EvtSel="NONE", EvtMax=100, ExtSvc=[podioevent], OutputLevel=DEBUG, )
#!/usr/bin/env python from Gaudi.Configuration import * from Configurables import K4DataSvc dsvc = K4DataSvc("EventDataSvc", input="test.root") from Configurables import PodioInput podioinput = PodioInput("PodioReader", collections=[ "EventHeader", "MCParticle", "SimCalorimeterHit" ]) from Configurables import Edm4hepReadAlg alg = Edm4hepReadAlg("Edm4hepReadAlg") #alg.HeaderCol.Path = "EventHeader" #alg.MCParticleCol.Path = "MCParticle" alg.SimCalorimeterHitCol.Path = "SimCalorimeterHit" # ApplicationMgr from Configurables import ApplicationMgr ApplicationMgr( TopAlg = [podioinput, alg], EvtSel = 'NONE', EvtMax = 10, ExtSvc = [dsvc], OutputLevel=DEBUG )
#!/usr/bin/env python import os import sys from Gaudi.Configuration import * ############################################################################## # Event Data Svc ############################################################################## from Configurables import K4DataSvc dsvc = K4DataSvc("EventDataSvc", input="you_input.root") ############################################################################## from Configurables import PodioInput ## set the input collection data podioinput = PodioInput("PodioReader", collections=["MCParticle", "ECALBarrel"]) ############################################################################## from Configurables import GearSvc gearSvc = GearSvc("GearSvc") gearSvc.GearXMLFile = "../Pandora/FullDetGear.xml" ############################################################################## from Configurables import PandoraPFAlg pandoralg = PandoraPFAlg("PandoraPFAlg") ## KEEP same with lcioinput name for the ReadXXX ########### pandoralg.ReadMCParticle = "MCParticle" pandoralg.ReadECALBarrel = "ECALBarrel" pandoralg.ReadECALEndcap = "ECALEndcap" pandoralg.ReadECALOther = "ECALOther" pandoralg.ReadHCALBarrel = "HCALBarrel" pandoralg.ReadHCALEndcap = "HCALEndcap"
geometry_path = os.path.join(os.getenv("DETCEPCV4ROOT"), "compact", geometry_option) if not os.path.exists(geometry_path): print("Can't find the compact geometry file: %s" % geometry_path) sys.exit(-1) from Configurables import GeoSvc geosvc = GeoSvc("GeoSvc") geosvc.compact = geometry_path ############################################################################## # Event Data Svc ############################################################################## from Configurables import K4DataSvc dsvc = K4DataSvc("EventDataSvc", input="test-detsim10.root") ############################################################################## # NTuple Svc ############################################################################## from Configurables import NTupleSvc ntsvc = NTupleSvc("NTupleSvc") ntsvc.Output = ["MyTuples DATAFILE='result.root' OPT='NEW' TYP='ROOT'"] ############################################################################## # DumpAlg ############################################################################## from Configurables import DumpIDAlg alg = DumpIDAlg("DumpAlg")
#!/usr/bin/env python from Gaudi.Configuration import * from Configurables import K4DataSvc #dsvc = K4DataSvc("EventDataSvc", input="detsim_ECAL_gamma_10000evt.root") dsvc = K4DataSvc("EventDataSvc", input="ECALonly_gamma_30degree_10000evt.root") ######################################################################### geometry_option = "CepC_v4-onlyECAL.xml" if not os.getenv("DETCEPCV4ROOT"): print("Can't find the geometry. Please setup envvar DETCEPCV4ROOT.") sys.exit(-1) geometry_path = os.path.join(os.getenv("DETCEPCV4ROOT"), "compact", geometry_option) if not os.path.exists(geometry_path): print("Can't find the compact geometry file: %s" % geometry_path) sys.exit(-1) from Configurables import GeomSvc geosvc = GeomSvc("GeomSvc") geosvc.compact = geometry_path ######################################################################## from Configurables import PodioInput podioinput = PodioInput("PodioReader", collections=[ "MCParticle", "EcalBarrelCollection", "EcalEndcapsCollection" ])