Newer
Older
#!/usr/bin/env python
from Gaudi.Configuration import *
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"
#"/junofs/users/wxfang/CEPC/CEPCOFF/doReco/reco_output/nnh_aa.e0.p0.00010_000000_rec.slcio"
"/junofs/users/wxfang/MyGit/tmp/fork_update_pandora/CEPCSW/Digi_sim_0.slcio"
"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 GearSvc
gearSvc = GearSvc("GearSvc")
gearSvc.GearXMLFile = "../Detector/DetCEPCv4/compact/FullDetGear.xml"
##############################################################################
from Configurables import PandoraPFAlg
pandoralg = PandoraPFAlg("PandoraPFAlg")
pandoralg.collections = [
"MCParticle:MCParticle",
"CalorimeterHit:ECALBarrel",
"CalorimeterHit:ECALEndcap",
"CalorimeterHit:ECALOther" ,
"CalorimeterHit:HCALBarrel",
"CalorimeterHit:HCALEndcap",
"CalorimeterHit:HCALOther" ,
"CalorimeterHit:MUON",
"CalorimeterHit:LCAL",
"CalorimeterHit:LHCAL",
"CalorimeterHit:BCAL",
"Vertex:KinkVertices",
"Vertex:ProngVertices",
"Vertex:SplitVertices",
"Vertex:V0Vertices",
"Track:MarlinTrkTracks",
"MCRecoCaloAssociation:RecoCaloAssociation_ECALBarrel"
]
pandoralg.WriteClusterCollection = "PandoraClusters"
pandoralg.WriteReconstructedParticleCollection = "PandoraPFOs"
pandoralg.WriteVertexCollection = "PandoraPFANewStartVertices"
pandoralg.AnaOutput = "Pandora_Ana.root"
pandoralg.PandoraSettingsDefault_xml = "../Reconstruction/PFA/Pandora/PandoraSettingsDefault.xml"
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
#### Do not chage the collection name, only add or delete ###############
pandoralg.TrackCollections = ["MarlinTrkTracks"]
pandoralg.ECalCaloHitCollections= ["ECALBarrel", "ECALEndcap", "ECALOther"]
pandoralg.HCalCaloHitCollections= ["HCALBarrel", "HCALEndcap", "HCALOther"]
pandoralg.LCalCaloHitCollections= ["LCAL"]
pandoralg.LHCalCaloHitCollections= ["LHCAL"]
pandoralg.MuonCaloHitCollections= ["MUON"]
pandoralg.MCParticleCollections = ["MCParticle"]
pandoralg.RelCaloHitCollections = ["RecoCaloAssociation_ECALBarrel", "RecoCaloAssociation_ECALEndcap", "RecoCaloAssociation_ECALOther", "RecoCaloAssociation_HCALBarrel", "RecoCaloAssociation_HCALEndcap", "RecoCaloAssociation_HCALOther", "RecoCaloAssociation_LCAL", "RecoCaloAssociation_LHCAL", "RecoCaloAssociation_MUON"]
pandoralg.RelTrackCollections = ["MarlinTrkTracksMCTruthLink"]
pandoralg.KinkVertexCollections = ["KinkVertices"]
pandoralg.ProngVertexCollections= ["ProngVertices"]
pandoralg.SplitVertexCollections= ["SplitVertices"]
pandoralg.V0VertexCollections = ["V0Vertices"]
pandoralg.ECalToMipCalibration = 160.0
pandoralg.HCalToMipCalibration = 34.8
pandoralg.ECalMipThreshold = 0.5
pandoralg.HCalMipThreshold = 0.3
pandoralg.ECalToEMGeVCalibration= 0.9 #for G2CD Digi, 1.007 for NewLDCaloDigi
pandoralg.HCalToEMGeVCalibration= 1.007
pandoralg.ECalToHadGeVCalibrationBarrel= 1.12 #very small effect
pandoralg.ECalToHadGeVCalibrationEndCap= 1.12
pandoralg.HCalToHadGeVCalibration= 1.07
pandoralg.MuonToMipCalibration= 10.0
pandoralg.DigitalMuonHits= 0
pandoralg.MaxHCalHitHadronicEnergy = 1.0
pandoralg.UseOldTrackStateCalculation= 0
pandoralg.AbsorberRadLengthECal= 0.2854
pandoralg.AbsorberIntLengthECal= 0.0101
pandoralg.AbsorberRadLengthHCal= 0.0569
pandoralg.AbsorberIntLengthHCal= 0.006
pandoralg.AbsorberRadLengthOther= 0.0569
pandoralg.AbsorberIntLengthOther= 0.006
##############################################################################
# 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, pandoralg, write],
TopAlg = [read, pandoralg],
EvtSel = 'NONE',