read.collections = [
"MCParticle:MCParticle",
#"SimCalorimeterHit:EcalBarrelSiliconCollection",
"CalorimeterHit:ECALBarrel",
"CalorimeterHit:ECALEndcap",
"CalorimeterHit:ECALOther",
######### HCAL will effect the reco efficiency close to gap region ######
"CalorimeterHit:HCALBarrel",
"CalorimeterHit:HCALEndcap",
"CalorimeterHit:HCALOther",
#"TrackerHit:VXDTrackerHits",
#"TrackerHit:SITTrackerHits",
"TrackerHit:SITSpacePoints",
"TrackerHit:TPCTrackerHits",
#"TrackerHit:SETTrackerHits",
"TrackerHit:SETSpacePoints",
#"TrackerHit:FTDStripTrackerHits",
"TrackerHit:FTDSpacePoints",
#"TrackerHit:FTDPixelTrackerHits",
"Track:ClupatraTrackSegments",
"Track:ClupatraTracks",
"Track:ForwardTracks",
"Track:SiTracks",
"Track:SubsetTracks",
"Track:MarlinTrkTracks",
"Vertex:KinkVertices",
"Vertex:ProngVertices",
"Vertex:V0Vertices",
"ReconstructedParticle:KinkRecoParticles",
"ReconstructedParticle:ProngRecoParticles",
"ReconstructedParticle:V0RecoParticles"
]
from Configurables import k4DataSvc
dsvc = k4DataSvc("EventDataSvc")
# read LCIO files
from Configurables import LCIOInput
read = LCIOInput("read")
read.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"
"/cefs/data/RecData/CEPC250/CEPC_v1/higgs/E250.Pnnh.whizard195/E250-CDR_ws.Pnnh.eU.pU.01232_rec.slcio"
]
read.collections = [
"SimTrackerHit:COILCollection",
"SimCalorimeterHit:EcalBarrelSiliconCollection", "MCParticle:MCParticle",
"SimTrackerHit:TPCCollection", "SimTrackerHit:VXDCollection",
"Cluster:ArborCharged", "Cluster:ArborNeutral",
"ReconstructedParticle:ArborChargedCore",
"ReconstructedParticle:ArborNeutralCore",
"MCRecoParticleAssociation:RecoMCTruthLink"
]
# write PODIO file
from Configurables import PodioOutput
write = PodioOutput("write")
write.filename = "test.root"
write.outputCommands = ["keep *"]
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=[read, write],
EvtSel='NONE',
from Configurables import k4DataSvc
dsvc = k4DataSvc("EventDataSvc")
# read LCIO files
from Configurables import LCIOInput
read = LCIOInput("read")
read.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"
]
read.collections = [
"SimTrackerHit:COILCollection",
#"SimCalorimeterHit:EcalBarrelSiliconCollection",
"MCParticle:MCParticle",
"SimTrackerHit:TPCCollection",
"SimTrackerHit:VXDCollection",
"MCRecoParticleAssociation:RecoMCTruthLink"
]
# write PODIO file
from Configurables import PodioOutput
write = PodioOutput("write")
write.filename = "test.root"
write.outputCommands = ["keep *"]
# ApplicationMgr
from Configurables import ApplicationMgr
ApplicationMgr(TopAlg=[read, write],
# read LCIO files
from Configurables import LCIOInput
lcioinput = LCIOInput("LCIOInput")
import glob
# inputlists = ["/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__00001.slcio"]
inputlists = ["/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__00003.slcio"]
# inputlists = glob.glob("/cefs/higgs/yudan/CEPC240/Reco_tpc_1800/qqh/Reco_qqh__*.slcio")
lcioinput.inputs = inputlists
lcioinput.collections = [
"MCParticle:MCParticle",
"CalorimeterHit:ECALBarrel",
"CalorimeterHit:ECALEndcap",
"CalorimeterHit:HCALBarrel",
"CalorimeterHit:HCALEndcap",
"CalorimeterHit:HCALOther",
"ReconstructedParticle:AncientPFOs",
"ReconstructedParticle:ArborLICHPFOs"
]
from Configurables import TotalInvMass
total_inv_mass = TotalInvMass("TotalInvMass")
# # write PODIO file
# from Configurables import PodioOutput
# write = PodioOutput("write")
# write.filename = "test.root"
# write.outputCommands = ["keep *"]
# ApplicationMgr
read = LCIOInput("read")
read.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"
#"/junofs/users/wxfang/CEPC/CEPCOFF/doReco/reco_output/nnh_aa.e0.p0.00010_000000_rec.slcio"
"/cefs/higgs/wxfang/cepc/Pandora/CaloDigi/gamma/Digi_sim_0.slcio"
]
read.collections = {
#"COILCollection" : "SimTrackerHit",
#"EcalBarrelSiliconCollection" : "SimCalorimeterHit",
"MCParticle": "MCParticle",
"ECALBarrel": "CalorimeterHit",
"ECALEndcap": "CalorimeterHit",
"ECALOther": "CalorimeterHit",
"HCALBarrel": "CalorimeterHit",
"HCALEndcap": "CalorimeterHit",
"HCALOther": "CalorimeterHit",
"MUON": "CalorimeterHit",
"LCAL": "CalorimeterHit",
"LHCAL": "CalorimeterHit",
"BCAL": "CalorimeterHit",
#"MarlinTrkTracks" : "Track"
#"TPCCollection" : "SimTrackerHit",
#"VXDCollection" : "SimTrackerHit"
}
##############################################################################
from Configurables import GearSvc
gearSvc = GearSvc("GearSvc")
gearSvc.GearXMLFile = "/junofs/users/wxfang/CEPC/CEPCOFF/doSim/fullDet/GearOutput.xml"
##############################################################################
from Configurables import PandoraPFAlg