Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
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
#!/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"
"/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
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"
pandoralg.ReadHCALOther = "HCALOther"
pandoralg.ReadMUON = "MUON"
pandoralg.ReadLCAL = "LCAL"
pandoralg.ReadLHCAL = "LHCAL"
pandoralg.ReadBCAL = "BCAL"
pandoralg.ReadKinkVertices = "KinkVertices"
pandoralg.ReadProngVertices = "ProngVertices"
pandoralg.ReadSplitVertices = "SplitVertices"
pandoralg.ReadV0Vertices = "V0Vertices"
pandoralg.ReadTracks = "MarlinTrkTracks"
pandoralg.WriteClusterCollection = "PandoraClusters"
pandoralg.WriteReconstructedParticleCollection = "PandoraPFOs"
pandoralg.WriteVertexCollection = "PandoraPFANewStartVertices"
pandoralg.AnaOutput = "/cefs/higgs/wxfang/cepc/Pandora/Ana/gamma/Ana_gamma_test.root"
pandoralg.PandoraSettingsDefault_xml = "/junofs/users/wxfang/MyGit/MarlinPandora/scripts/PandoraSettingsDefault_wx.xml"
#### 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',
EvtMax = 10,
ExtSvc = [dsvc, gearSvc],
OutputLevel=INFO
)