Skip to content
Snippets Groups Projects
Commit ad836d98 authored by zoujh@ihep.ac.cn's avatar zoujh@ihep.ac.cn
Browse files

update the example options for reading LCIO data

parent 5cf624dc
No related branches found
No related tags found
No related merge requests found
#!/usr/bin/env python #!/usr/bin/env python
from Gaudi.Configuration import * from Gaudi.Configuration import *
from Configurables import k4DataSvc
dsvc = k4DataSvc("EventDataSvc")
from Configurables import LCIODataSvc, CEPCDataSvc # read LCIO files
svcname = "LCIODataSvc"
rsvc = LCIODataSvc(svcname, inputs = [
"/cefs/data/FullSim/CEPC240/CEPC_v4/higgs/smart_final_states/E240.Pffh_invi.e0.p0.whizard195//ffh_inv.e0.p0.00001_1000_sim.slcio"
])
wsvc = CEPCDataSvc("EventDataSvc")
from Configurables import PlcioReadAlg
alg = PlcioReadAlg("PlcioReadAlg")
alg.InputCol.Path = "MCParticle"
alg.HeaderCol.Path = "EventHeader"
from Configurables import LCIOInput from Configurables import LCIOInput
lcioinput = LCIOInput("LCIOReader", collections=[ read = LCIOInput("read")
"EventHeader", read.inputs = [
"MCParticle", "/cefs/data/FullSim/CEPC240/CEPC_v4/higgs/smart_final_states/E240.Pffh_invi.e0.p0.whizard195/ffh_inv.e0.p0.00001_1000_sim.slcio"
"TPCCollection" ]
#"EventHeader", read.collections = [
#"MCParticle", "MCParticle:MCParticle",
#"COILCollection", "SimCalorimeterHit:EcalBarrelSiliconCollection",
#"EcalBarrelSiliconCollection", #"CalorimeterHit:ECALBarrel",
#"EcalBarrelSiliconPreShowerCollection", #"CalorimeterHit:ECALEndcap",
#"EcalEndcapRingCollection", #"CalorimeterHit:ECALOther" ,
#"EcalEndcapRingPreShowerCollection", ########## HCAL will effect the reco efficiency close to gap region ######
#"EcalEndcapSiliconCollection", #"CalorimeterHit:HCALBarrel",
#"EcalEndcapSiliconPreShowerCollection", #"CalorimeterHit:HCALEndcap",
#"FTD_PIXELCollection", #"CalorimeterHit:HCALOther",
#"FTD_STRIPCollection", ##"TrackerHit:VXDTrackerHits",
#"HcalBarrelCollection", ##"TrackerHit:SITTrackerHits",
#"HcalEndCapRingsCollection", #"TrackerHit:SITSpacePoints",
#"HcalEndCapsCollection", #"TrackerHit:TPCTrackerHits",
#"LumiCalCollection", ##"TrackerHit:SETTrackerHits",
#"MuonBarrelCollection", #"TrackerHit:SETSpacePoints",
#"MuonEndCapCollection", ##"TrackerHit:FTDStripTrackerHits",
#"SETCollection", #"TrackerHit:FTDSpacePoints",
#"SITCollection", ##"TrackerHit:FTDPixelTrackerHits",
#"TPCCollection", #"Track:ClupatraTrackSegments",
#"TPCSpacePointCollection", #"Track:ClupatraTracks",
#"VXDCollection" #"Track:ForwardTracks",
]) #"Track:SiTracks",
lcioinput.DataSvc = svcname #"Track:SubsetTracks",
#"Track:MarlinTrkTracks",
#"Vertex:KinkVertices",
#"Vertex:ProngVertices",
#"Vertex:V0Vertices",
#"ReconstructedParticle:KinkRecoParticles",
#"ReconstructedParticle:ProngRecoParticles",
#"ReconstructedParticle:V0RecoParticles"
]
from Configurables import PodioOutput from Configurables import PodioOutput
plcioout = PodioOutput("PlcioWriter") write = PodioOutput("write")
plcioout.filename = "lcio2plcio.root" write.filename = "lcio2plcio.root"
plcioout.outputCommands = ["keep *"] write.outputCommands = ["keep *"]
# ApplicationMgr # ApplicationMgr
from Configurables import ApplicationMgr from Configurables import ApplicationMgr
ApplicationMgr( TopAlg = [lcioinput, alg, plcioout], ApplicationMgr( TopAlg = [read, write],
EvtSel = 'NONE', EvtSel = 'NONE',
EvtMax = 10, EvtMax = 10,
ExtSvc = [rsvc, wsvc], ExtSvc = [dsvc],
OutputLevel=DEBUG OutputLevel=DEBUG
) )
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment