ATLAS Offline Software
Loading...
Searching...
No Matches
tauMonitorAlgorithm Namespace Reference

Functions

 tauMonitoringConfig (inputFlags, tauContainer='TauJets', **kwargs)

Detailed Description

@file tauMonitorAlgorithm.py
@date 2018-01-11
@brief tau python configuration for the Run III AthenaMonitoring package

Function Documentation

◆ tauMonitoringConfig()

tauMonitorAlgorithm.tauMonitoringConfig ( inputFlags,
tauContainer = 'TauJets',
** kwargs )
Function to configures some algorithms in the monitoring system.

Definition at line 12 of file tauMonitorAlgorithm.py.

12def tauMonitoringConfig(inputFlags,tauContainer='TauJets',**kwargs):
13 '''Function to configures some algorithms in the monitoring system.'''
14
15
16 if not (inputFlags.Common.isOnline == 'online' or inputFlags.Input.isMC):
17 kwargs.setdefault('useReadyFilterTool', True)
18 else:
19 kwargs.setdefault('useReadyFilterTool', False)
20
21 from AthenaConfiguration.ComponentFactory import CompFactory
22 from AthenaMonitoring.BadLBFilterToolConfig import LArBadLBFilterToolCfg
23 from AthenaMonitoring.AtlasReadyFilterConfig import AtlasReadyFilterCfg
24
25
26
27 # The following class will make a sequence, configure algorithms, and link
28 # them to GenericMonitoringTools
29 from AthenaMonitoring import AthMonitorCfgHelper
30 cfgHelper = AthMonitorCfgHelper(inputFlags, monName='tauMonitorAlgFamily')
31 cfg = cfgHelper.result()
32
33 if kwargs['useReadyFilterTool'] and 'ReadyFilterTool' not in kwargs:
34 readyFilterTool = cfg.popToolsAndMerge(AtlasReadyFilterCfg(inputFlags))
35 kwargs['ReadyFilterTool'] = readyFilterTool
36
37
38
44
45
46 tauMonitorAlgorithm = CompFactory.tauMonitorAlgorithm
47 tauMonAlgBA = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgBA', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
48 tauMonAlgCR = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgCR', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
49 tauMonAlgEC = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgEC', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
50 tauMonAlgGlobal = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgGlobal', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
51 tauMonAlgTauTrig1 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig1', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
52 tauMonAlgTauTrig2 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig2', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
53 tauMonAlgTauTrig3 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig3', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
54 tauMonAlgTauTrig4 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig4', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
55 tauMonAlgTauTrig5 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig5', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
56 tauMonAlgTauTrig6 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig6', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
57 tauMonAlgTauTrig7 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig7', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
58 tauMonAlgTauTrig8 = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgTauTrig8', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
59 tauMonAlgEleTrig = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgEleTrig', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
60 tauMonAlgJetTrig = cfgHelper.addAlgorithm( tauMonitorAlgorithm, name='tauMonAlgJetTrig', addFilterTools = [LArBadLBFilterToolCfg(inputFlags)])
61
62
63
64
70
71 tauMonAlgTauTrig1.TriggerChain = "HLT_tau25_mediumGNTau_L1eTAU20"
72 tauMonAlgTauTrig2.TriggerChain = "HLT_tau25_mediumGNTau_L1cTAU20M"
73 tauMonAlgTauTrig3.TriggerChain = "HLT_tau160_mediumGNTau_L1eTAU140"
74 tauMonAlgTauTrig4.TriggerChain = "HLT_tau35_mediumGNTau_tau25_mediumGNTau_03dRAB30_L1cTAU30M_2cTAU20M_DR-eTAU30eTAU20-jJ55"
75 tauMonAlgTauTrig5.TriggerChain = "HLT_tau35_mediumGNTau_tau25_mediumGNTau_03dRAB_L1cTAU30M_2cTAU20M_4jJ30p0ETA25"
76 tauMonAlgTauTrig6.TriggerChain = "HLT_tau40_mediumGNTau_tau35_mediumGNTau_03dRAB_L1cTAU35M_2cTAU30M_2jJ55_3jJ50"
77 tauMonAlgTauTrig7.TriggerChain = "HLT_tau80_mediumGNTau_tau35_mediumGNTau_03dRAB30_L1eTAU80_2cTAU30M_DR-eTAU30eTAU20"
78 tauMonAlgTauTrig8.TriggerChain = "HLT_tau80_mediumGNTau_tau60_mediumGNTau_03dRAB_L1eTAU80_2eTAU60"
79
80
81 tauMonAlgEleTrig.TriggerChain = "HLT_e[2-9][0-9]_.*"
82 tauMonAlgJetTrig.TriggerChain = "HLT_j[2-9][0-9]_.*"
83
84 tauMonAlgBA.TauRecContainer = tauContainer
85 tauMonAlgBA.etaMin = -1.
86 tauMonAlgBA.etaMax = 1.3
87 tauMonAlgCR.TauRecContainer = tauContainer
88 tauMonAlgCR.etaMin = 1.3
89 tauMonAlgCR.etaMax = 1.7
90 tauMonAlgEC.TauRecContainer = tauContainer
91 tauMonAlgEC.etaMin = 1.7
92 tauMonAlgEC.etaMax = 3.5
93 tauMonAlgGlobal.TauRecContainer = tauContainer
94 tauMonAlgGlobal.etaMin = -100
95 tauMonAlgGlobal.etaMax = 100
96 tauMonAlgTauTrig1.TauRecContainer = tauContainer
97 tauMonAlgTauTrig1.etaMin = -100
98 tauMonAlgTauTrig1.etaMax = 100
99 tauMonAlgTauTrig2.TauRecContainer = tauContainer
100 tauMonAlgTauTrig2.etaMin = -100
101 tauMonAlgTauTrig2.etaMax = 100
102 tauMonAlgTauTrig3.TauRecContainer = tauContainer
103 tauMonAlgTauTrig3.etaMin = -100
104 tauMonAlgTauTrig3.etaMax = 100
105 tauMonAlgTauTrig4.TauRecContainer = tauContainer
106 tauMonAlgTauTrig4.etaMin = -100
107 tauMonAlgTauTrig4.etaMax = 100
108 tauMonAlgTauTrig5.TauRecContainer = tauContainer
109 tauMonAlgTauTrig5.etaMin = -100
110 tauMonAlgTauTrig5.etaMax = 100
111 tauMonAlgTauTrig6.TauRecContainer = tauContainer
112 tauMonAlgTauTrig6.etaMin = -100
113 tauMonAlgTauTrig6.etaMax = 100
114 tauMonAlgTauTrig7.TauRecContainer = tauContainer
115 tauMonAlgTauTrig7.etaMin = -100
116 tauMonAlgTauTrig7.etaMax = 100
117 tauMonAlgTauTrig8.TauRecContainer = tauContainer
118 tauMonAlgTauTrig8.etaMin = -100
119 tauMonAlgTauTrig8.etaMax = 100
120
121 tauMonAlgEleTrig.TauRecContainer = tauContainer
122 tauMonAlgEleTrig.etaMin = -100
123 tauMonAlgEleTrig.etaMax = 100
124 tauMonAlgJetTrig.TauRecContainer = tauContainer
125 tauMonAlgJetTrig.etaMin = -100
126 tauMonAlgJetTrig.etaMax = 100
127
128 tauMonAlgBA.kinGroupName = 'tauMonKinGroupBA'
129 tauMonAlgCR.kinGroupName = 'tauMonKinGroupCR'
130 tauMonAlgEC.kinGroupName = 'tauMonKinGroupEC'
131 tauMonAlgGlobal.kinGroupName = 'tauMonKinGroupGlobal'
132 tauMonAlgTauTrig1.kinGroupName = 'tauMonKinGroupTauTrig1'
133 tauMonAlgTauTrig2.kinGroupName = 'tauMonKinGroupTauTrig2'
134 tauMonAlgTauTrig3.kinGroupName = 'tauMonKinGroupTauTrig3'
135 tauMonAlgTauTrig4.kinGroupName = 'tauMonKinGroupTauTrig4'
136 tauMonAlgTauTrig5.kinGroupName = 'tauMonKinGroupTauTrig5'
137 tauMonAlgTauTrig6.kinGroupName = 'tauMonKinGroupTauTrig6'
138 tauMonAlgTauTrig7.kinGroupName = 'tauMonKinGroupTauTrig7'
139 tauMonAlgTauTrig8.kinGroupName = 'tauMonKinGroupTauTrig8'
140
141 tauMonAlgEleTrig.kinGroupName = 'tauMonKinGroupEleTrig'
142 tauMonAlgJetTrig.kinGroupName = 'tauMonKinGroupJetTrig'
143
144
147
148 # Add a generic monitoring tool (a "group" in old language). The returned
149 # object here is the standard GenericMonitoringTool.
150 myKinGroupBA = cfgHelper.addGroup(alg=tauMonAlgBA, name='tauMonKinGroupBA', topPath='Tau/TauB/' )
151 myKinGroupCR = cfgHelper.addGroup(alg=tauMonAlgCR, name='tauMonKinGroupCR', topPath='Tau/TauCR/' )
152 myKinGroupEC = cfgHelper.addGroup(alg=tauMonAlgEC, name='tauMonKinGroupEC', topPath='Tau/TauE/' )
153 myKinGroupGlobal = cfgHelper.addGroup(alg=tauMonAlgGlobal, name='tauMonKinGroupGlobal', topPath='Tau/' )
154 myKinGroupTauTrig1 = cfgHelper.addGroup(alg=tauMonAlgTauTrig1, name='tauMonKinGroupTauTrig1', topPath='Tau/Trigger/tauTrigger1/' )
155 myKinGroupTauTrig2 = cfgHelper.addGroup(alg=tauMonAlgTauTrig2, name='tauMonKinGroupTauTrig2', topPath='Tau/Trigger/tauTrigger2/' )
156 myKinGroupTauTrig3 = cfgHelper.addGroup(alg=tauMonAlgTauTrig3, name='tauMonKinGroupTauTrig3', topPath='Tau/Trigger/tauTrigger3/' )
157 myKinGroupTauTrig4 = cfgHelper.addGroup(alg=tauMonAlgTauTrig4, name='tauMonKinGroupTauTrig4', topPath='Tau/Trigger/tauTrigger4/' )
158 myKinGroupTauTrig5 = cfgHelper.addGroup(alg=tauMonAlgTauTrig5, name='tauMonKinGroupTauTrig5', topPath='Tau/Trigger/tauTrigger5/' )
159 myKinGroupTauTrig6 = cfgHelper.addGroup(alg=tauMonAlgTauTrig6, name='tauMonKinGroupTauTrig6', topPath='Tau/Trigger/tauTrigger6/' )
160 myKinGroupTauTrig7 = cfgHelper.addGroup(alg=tauMonAlgTauTrig7, name='tauMonKinGroupTauTrig7', topPath='Tau/Trigger/tauTrigger7/' )
161 myKinGroupTauTrig8 = cfgHelper.addGroup(alg=tauMonAlgTauTrig8, name='tauMonKinGroupTauTrig8', topPath='Tau/Trigger/tauTrigger8/' )
162 myKinGroupEleTrig = cfgHelper.addGroup(alg=tauMonAlgEleTrig, name='tauMonKinGroupEleTrig', topPath='Tau/Trigger/EleTrig/' )
163 myKinGroupJetTrig = cfgHelper.addGroup(alg=tauMonAlgJetTrig, name='tauMonKinGroupJetTrig', topPath='Tau/Trigger/JetTrig/' )
164
165
166 naming= {
167 'BA': "Tau_TauB_",
168 'CR': "Tau_TauCR_",
169 'EC': "Tau_TauE_",
170 'Global': "",
171 'EleTrig': "emTriggered_",
172 'JetTrig': "jetTriggered_",
173 'TauTrig1': "tauTriggered1_",
174 'TauTrig2': "tauTriggered2_",
175 'TauTrig3': "tauTriggered3_",
176 'TauTrig4': "tauTriggered4_",
177 'TauTrig5': "tauTriggered5_",
178 'TauTrig6': "tauTriggered6_",
179 'TauTrig7': "tauTriggered7_",
180 'TauTrig8': "tauTriggered8_",
181
182 }
183
184 def namer(variables,histogramName, folder, algorithm):
185 if folder != "":
186 #This wouldn't have been necessary except the calo folder in the naming convention:
187 #This adds the Tau prefix to the histograms in the calo folder
188 if postfix == "Global" and (folder =="Calo" or folder=="Track" or folder=="Track_allTracks" or folder=="Track_leadTrack" or folder=="Cluster" ):
189 return variables+";"+"Tau_"+naming[algorithm]+folder+"_"+histogramName
190 else:
191 return variables+";"+naming[algorithm]+folder+"_"+histogramName
192 else:
193 return variables+";"+naming[algorithm]+folder+histogramName
194
195 # Add a GMT for the other example monitor algorithm
196 # anotherGroup = cfgHelper.addGroup(anotherTauMonAlg,'tauMonitor')
197
198
200
201 # define which ID is running
202 tauid = "GNTau" if "TMTauJets" in tauContainer else "RNN"
203
204 # Configure histograms
205 for itup in [(myKinGroupBA,'BA'),
206 (myKinGroupCR,'CR'),
207 (myKinGroupEC,'EC'),
208 (myKinGroupGlobal,'Global'),
209 (myKinGroupEleTrig,'EleTrig'),
210 (myKinGroupJetTrig,'JetTrig'),
211 (myKinGroupTauTrig1, 'TauTrig1'),
212 (myKinGroupTauTrig2, 'TauTrig2'),
213 (myKinGroupTauTrig3, 'TauTrig3'),
214 (myKinGroupTauTrig4, 'TauTrig4'),
215 (myKinGroupTauTrig5, 'TauTrig5'),
216 (myKinGroupTauTrig6, 'TauTrig6'),
217 (myKinGroupTauTrig7, 'TauTrig7'),
218 (myKinGroupTauTrig8, 'TauTrig8')]:
219
220 (igroup, postfix) = itup
221
222 PHIMIN = -pi
223 PHIMAX = pi
224
225 folder = ""
226
227 if(postfix =="BA" or postfix =="CR" or postfix=="EC" or postfix == "Global" or postfix.startswith('TauTrig') or postfix=="EleTrig" or postfix =="JetTrig"):
228
229
230 #potentialHigh Pt replacement
231 igroup.defineHistogram(namer('tauEtaEt15,tauPhiEt15','tauPhiVsEta_et15','',postfix), type='TH2F', title='Eta Vs phi (Et>15GeV);Eta;Phi',
232 xbins=30,xmin=-2.55,xmax=2.55,ybins=32,ymin=PHIMIN,ymax=PHIMAX)
233
234 igroup.defineHistogram(namer('nTauCandidates', 'nTauCandidates',"",postfix), title='Number of tau candidates;Number of Taus per Event',
235 xbins=61, xmin=-0.5, xmax=60.5,path=folder)
236
237
238 igroup.defineHistogram(namer('tauCharge','tauCharge','',postfix), title='Charge of tau candidates;Charge;Number of Candidates',
239 xbins=11, xmin=-5.5, xmax=5.5,path=folder)
240
241
242 igroup.defineHistogram(namer(tauid+'JetScore',tauid+'JetScore','',postfix), title= tauid +' Jet Score',
243 xbins=100, xmin=0, xmax=1,path=folder)
244
245 igroup.defineHistogram(namer(tauid+'JetScoreSigTrans',tauid+'JetScoreSigTrans','',postfix), title=tauid+' Jet Score Sig Trans',
246 xbins=48, xmin=0, xmax=1.1,path=folder)
247
248 igroup.defineHistogram(namer('RNNEleScore','RNNEleScore','',postfix), title='RNN Ele Score',
249 xbins=100, xmin=0, xmax=1,path=folder)
250
251 igroup.defineHistogram(namer('RNNEleScoreSigTrans','RNNEleScoreSigTrans','',postfix), title='RNN Ele Score Sig Trans',
252 xbins=48, xmin=0, xmax=1.1,path=folder)
253
254 igroup.defineHistogram(namer('tauEt','tauEt','',postfix), title='Et of tau candidates;Transverse Energy (GeV);Number of Candidates',
255 xbins=60, xmin=0., xmax=300.,path=folder)
256
257 igroup.defineHistogram(namer('tauEtaEXP','tauEtaEXP','',postfix), title='Eta of tau candidates;Eta;Number of Candidates',
258 xbins=51, xmin=-2.55, xmax=2.55,path=folder)
259
260
261
262 igroup.defineHistogram(namer('tauEta','tauEta','',postfix), title='Eta of tau candidates;Eta;Number of Candidates',
263 xbins=51, xmin=-2.55, xmax=2.55,path=folder)
264
265 igroup.defineHistogram(namer('LB,tauEta','tauEtaVsLB','',postfix), type='TH2F', title='Tau Eta Vs Lumiblock;Lumiblock;eta',
266 xbins=1200,xmin=0,xmax=1200,ybins=51,ymin=-2.55,ymax=2.55,path=folder)
267
268 igroup.defineHistogram(namer('tauPhi','tauPhi','',postfix), title='Phi of tau candidates;Phi;Number of Candidates',
269 xbins=65, xmin=PHIMIN-0.098174/2., xmax=PHIMAX+0.098174/2.,path=folder)
270
271 igroup.defineHistogram(namer('tauEta,tauPhi','tauPhiVsEta','',postfix), type='TH2F', title='Phi Vs Eta;Eta;Phi',
272 xbins=30,xmin=-2.55,xmax=2.55,ybins=32,ymin=PHIMIN,ymax=PHIMAX,path=folder)
273
274
275 igroup.defineHistogram(namer('LB,tauPhi','tauPhiVsLB','',postfix), type='TH2F', title='Tau Phi vs Lumiblock;Lumiblock;Phi',
276 xbins=1200,xmin=0,xmax=1200,ybins=80,ymin=PHIMIN,ymax=PHIMAX,path=folder)
277
278 if(postfix =="BA" or postfix =="CR" or postfix=="EC" or postfix.startswith('TauTrig')):
279
280 igroup.defineHistogram(namer('tauPhi'+tauid+'Loose','phi','Identification_'+tauid+'Loose',postfix), title='Phi of tau candidates ('+tauid+'Loose) ;Phi;Number of Candidates',
281 xbins=65, xmin=PHIMIN-0.098174/2., xmax=PHIMAX+0.098174/2., path=folder+"Identification/"+tauid+"Loose" )
282
283 igroup.defineHistogram(namer('tauEta'+tauid+'Loose','eta','Identification_'+tauid+'Loose',postfix), title='Eta of tau candidates ('+tauid+'Loose) ;Eta;Number of Candidates',
284 xbins=51, xmin=-2.55, xmax=2.55, path=folder+"Identification/"+tauid+"Loose")
285
286 igroup.defineHistogram(namer('tauEt'+tauid+'Loose','et','Identification_'+tauid+'Loose',postfix), title='Et of tau candidates;Transverse Energy (GeV);Number of Candidates',
287 xbins=60, xmin=0., xmax=300.,path=folder+"Identification/"+tauid+"Loose")
288
289 igroup.defineHistogram(namer('NumTracks'+tauid+'Loose','NumTracks','Identification_'+tauid+'Loose',postfix), title='Number Of Tracks for Tau Candidates ('+tauid+'Loose);Number Of Tracks;Number Of Candidates',
290 xbins=21, xmin=-0.5, xmax=20.5,path=folder+"Identification/"+tauid+"Loose")
291
292 igroup.defineHistogram(namer('tauPhi'+tauid+'Medium','phi','Identification_'+tauid+'Medium',postfix), title='Phi of tau candidates ('+tauid+'Medium) ;Phi;Number of Candidates',
293 xbins=65, xmin=PHIMIN-0.098174/2., xmax=PHIMAX+0.098174/2., path=folder+"Identification/"+tauid+"Medium" )
294
295 igroup.defineHistogram(namer('tauEta'+tauid+'Medium','eta','Identification_'+tauid+'Medium',postfix), title='Eta of tau candidates ('+tauid+'Medium) ;Eta;Number of Candidates',
296 xbins=51, xmin=-2.55, xmax=2.55, path=folder+"Identification/"+tauid+"Medium")
297
298 igroup.defineHistogram(namer('tauEt'+tauid+'Medium','et','Identification_'+tauid+'Medium',postfix), title='Et of tau candidates;Transverse Energy (GeV);Number of Candidates',
299 xbins=60, xmin=0., xmax=300.,path=folder+"Identification/"+tauid+"Medium")
300
301 igroup.defineHistogram(namer('NumTracks'+tauid+'Medium','NumTracks','Identification_'+tauid+'Medium',postfix), title='Number Of Tracks for Tau Candidates ('+tauid+'Medium);Number Of Tracks;Number Of Candidates',
302 xbins=21, xmin=-0.5, xmax=20.5,path=folder+"Identification/"+tauid+"Medium")
303
304 igroup.defineHistogram(namer('tauPhi'+tauid+'Tight','phi','Identification_'+tauid+'Tight',postfix), title='Phi of tau candidates ('+tauid+'Tight) ;Phi;Number of Candidates',
305 xbins=65, xmin=PHIMIN-0.098174/2., xmax=PHIMAX+0.098174/2., path=folder+"Identification/"+tauid+"Tight" )
306
307 igroup.defineHistogram(namer('tauEta'+tauid+'Tight','eta','Identification_'+tauid+'Tight',postfix), title='Eta of tau candidates ('+tauid+'Tight) ;Eta;Number of Candidates',
308 xbins=51, xmin=-2.55, xmax=2.55, path=folder+"Identification/"+tauid+"Tight")
309
310 igroup.defineHistogram(namer('tauEt'+tauid+'Tight','et','Identification_'+tauid+'Tight',postfix), title='Et of tau candidates;Transverse Energy (GeV);Number of Candidates',
311 xbins=60, xmin=0., xmax=300.,path=folder+"Identification/"+tauid+"Tight")
312
313 igroup.defineHistogram(namer('NumTracks'+tauid+'Tight','NumTracks','Identification_'+tauid+'Tight',postfix), title='Number Of Tracks for Tau Candidates ('+tauid+'Tight);Number Of Tracks;Number Of Candidates',
314 xbins=21, xmin=-0.5, xmax=20.5,path=folder+"Identification/"+tauid+"Tight")
315
316
317
318
319 igroup.defineHistogram(namer('tauPhiEt15'+tauid+'Loose','phi','Identification_'+tauid+'Loose15GeV',postfix), title='Phi of tau candidates (Et>15, '+tauid+'Loose) ;Phi;Number of Candidates',
320 xbins=65, xmin=PHIMIN-0.098174/2., xmax=PHIMAX+0.098174/2., path=folder+"Identification/"+tauid+"Loose15GeV" )
321
322 igroup.defineHistogram(namer('tauEtaEt15'+tauid+'Loose','eta','Identification_'+tauid+'Loose15GeV',postfix), title='Eta of tau candidates (Et>15, '+tauid+'Loose) ;Eta;Number of Candidates',
323 xbins=51, xmin=-2.55, xmax=2.55, path=folder+"Identification/"+tauid+"Loose15GeV")
324
325 igroup.defineHistogram(namer('nClustersEt15'+tauid+'Loose','nCluster','Identification_'+tauid+'Loose15GeV',postfix), title='Number Of CaloTopoClusters (Et>15,'+tauid+'Loose);Number Of Clusters;Number Of Candidates',
326 xbins=40, xmin=0., xmax=40.,path=folder+"Identification/"+tauid+"Loose15GeV" )
327
328 igroup.defineHistogram(namer('NumTracksEt15'+tauid+'Loose','NumTracks','Identification_'+tauid+'Loose15GeV',postfix), title='Number Of Tracks for Tau Candidates (Et>15,'+tauid+'Loose);Number Of Tracks;Number Of Candidates',
329 xbins=21, xmin=-0.5, xmax=20.5,path=folder+"Identification/"+tauid+"Loose15GeV")
330
331 igroup.defineHistogram(namer('tauEtEt15'+tauid+'Loose','et','Identification_'+tauid+'Loose15GeV',postfix), title='Et of tau candidates;Transverse Energy (GeV);Number of Candidates',
332 xbins=60, xmin=0., xmax=300.,path=folder+"Identification/"+tauid+"Loose15GeV")
333
334 igroup.defineHistogram(namer('panModeEt15'+tauid+'Loose','panMode','Identification_'+tauid+'Loose15GeV',postfix), title='tau decay mode from panTau upon Jet'+tauid+'SigMedium;mode',
335 xbins=5, xmin=0., xmax=5., path=folder+"Identification/"+tauid+"Loose15GeV", xlabels=["1p0n","1p1n","1pXn","3p0n","3pXn"])
336
337 igroup.defineHistogram(namer('tauEtaEt15'+tauid+'Loose,tauPhiEt15'+tauid+'Loose','tauPhiVsEta_et15_'+tauid+'Loose','',postfix), type='TH2F', title='Phi vs Eta (Et>15, '+tauid+'Loose) ;Eta;Phi',
338 xbins=30,xmin=-2.55,xmax=2.55,ybins=32,ymin=PHIMIN,ymax=PHIMAX)
339
340
341 igroup.defineHistogram(namer('jetSeedEta','jetSeedEta','Calo',postfix), title='Calorimeter eta of tau candidates;Eta;Numbers of Candidates',path=folder+"Calo",
342 xbins=50, xmin=-2.5, xmax=2.5 )
343
344 igroup.defineHistogram(namer('jetSeedPhi','jetSeedPhi','Calo',postfix), title='Calorimeter phi of tau candidates;Phi;Numbers of Candidates',path=folder+"Calo",
345 xbins=64, xmin=PHIMIN, xmax=PHIMAX )
346
347 igroup.defineHistogram(namer('jetSeedPt','jetSeedPt','Calo',postfix), title='Calorimeter pt of tau candidates;Et;Numbers of Candidates',path=folder+"Calo",
348 xbins=100, xmin=0., xmax=300 )
349
350 igroup.defineHistogram(namer('centFrac,LB','CentFracVsLB','Calo',postfix), type='TH2F', title='Centrality Fraction vs Lumiblock;Centrality Fraction;Lumiblock', path=folder+"Calo",
351 xbins=51,xmin=0,xmax=1.02,ybins=1200,ymin=0.,ymax=1200.)
352
353 igroup.defineHistogram('isolFrac,LB', type='TH2F', title='Isolation Fraction vs Lumiblock;Isolation Fraction;Lumiblock', path=folder+"Calo",
354 xbins=51,xmin=0,xmax=1.02,ybins=1200,ymin=0.,ymax=1200.)
355
356
357
358 igroup.defineHistogram(namer('muonVeto','muonVeto','Identification',postfix), title='Muon Veto',path=folder+"Identification",
359 xbins=2, xmin=-0.5, xmax=1.5, xlabels=["False","True"] )
360
361 igroup.defineHistogram(namer('tau'+tauid+'Loose','tau'+tauid+'Loose','Identification',postfix), title='Identification Flag: tau'+tauid+'Loose',path=folder+"Identification",
362 xbins=2, xmin=-0.5, xmax=1.5 , xlabels=["False","True"])
363
364 igroup.defineHistogram(namer('tau'+tauid+'Medium','tau'+tauid+'Medium','Identification',postfix), title='Identification Flag: tau'+tauid+'Medium',path=folder+"Identification",
365 xbins=2, xmin=-0.5, xmax=1.5 , xlabels=["False","True"])
366
367 igroup.defineHistogram(namer('tau'+tauid+'Tight','tau'+tauid+'Tight','Identification',postfix), title='Identification Flag: tau'+tauid+'Tight',path=folder+"Identification",
368 xbins=2, xmin=-0.5, xmax=1.5, xlabels=["False","True"])
369
370
371 #Substructure Histograms
372 igroup.defineHistogram(namer('panModeSubstructure','panMode','SubStructure',postfix), title='tau decay mode from PanTau upon JetBDTSigMedium;mode',
373 xbins=8, xmin=0., xmax=8., path=folder+"SubStructure")
374
375 igroup.defineHistogram(namer('BDTScoreAsP0','BDTScoreAsP0','SubStructure',postfix), title='BDT score indentifying pion zero ; BDT score',path=folder+"SubStructure",
376 xbins=15, xmin=-0, xmax=1.2 )
377
378 igroup.defineHistogram(namer('EMFracTrk','EMFracTrk','SubStructure',postfix), title='Ratio of pt to shot electromagnetic energy for associated tracks; track pt ratio in EM',path=folder+"SubStructure",
379 xbins=15, xmin=0, xmax=1.5 )
380
381 igroup.defineHistogram(namer('nNeutPFO','nNeutPFO','SubStructure',postfix), title='_NumNeutPFO;Number of neutral ParticleFlow objects ; PFO number',path=folder+"SubStructure",
382 xbins=20, xmin=0, xmax=20. )
383
384 igroup.defineHistogram(namer('nShot','nShot','SubStructure',postfix), title='number of shots ; shot number ',path=folder+"SubStructure",
385 xbins=20, xmin=0, xmax=20. )
386
387 igroup.defineHistogram(namer('panEta','panEta','SubStructure',postfix), title='tau Eta from PanTau ; substructure Eta',path=folder+"SubStructure",
388 xbins=16, xmin=-3.2, xmax=3.2 )
389
390 igroup.defineHistogram(namer('panPhi','panPhi','SubStructure',postfix), title='tau Phi from PanTau ; substructure Phi',path=folder+"SubStructure",
391 xbins=16, xmin=PHIMIN, xmax=PHIMAX )
392
393 igroup.defineHistogram(namer('panPt','panPt','SubStructure',postfix), title='tau Pt from PanTau ; substructure pt (GeV)',path=folder+"SubStructure",
394 xbins=20, xmin=0, xmax=200 )
395
396 igroup.defineHistogram(namer('PSSFrac','PSSFracEffCluster','SubStructure',postfix), title='Energy fraction for PreSampling and sampling layers in effective clusters in shot;Sampling Energy Fraction',path=folder+"SubStructure",
397 xbins=10, xmin=0, xmax=1. )
398
399 igroup.defineHistogram(namer('ptRatioEflowApprox','ptRatioApprox','SubStructure',postfix), title='Ratio of pt to shot total energy for associated tracks; track pt ratio',
400 xbins=20, xmin=0., xmax=2., path=folder+"SubStructure" )
401
402
403 if(postfix =="BA" or postfix =="CR" or postfix=="EC" or postfix.startswith('TauTrig') or postfix=="Global"):
404 #Cluster Histograms
405 igroup.defineHistogram(namer('clusterLogEt','logEt','Cluster',postfix), title='log(cluster E_{T});log(cluster E_{T});Entries',path=folder+"Cluster",
406 xbins=50, xmin=1.5, xmax=6.5 )
407 igroup.defineHistogram(namer('clusterEta','eta','Cluster',postfix), title='cluster eta;cluster eta;Entries',path=folder+"Cluster",
408 xbins=50, xmin=-2.6, xmax=2.6 )
409 igroup.defineHistogram(namer('clusterPhi','phi','Cluster',postfix), title='cluster phi;cluster phi;Entries',path=folder+"Cluster",
410 xbins=32, xmin=-3.2, xmax=3.2 )
411 igroup.defineHistogram(namer('clusterSecondR','secondR','Cluster',postfix), title='cluster SECOND_R;cluster SECOND_R;Entries',path=folder+"Cluster",
412 xbins=50, xmin=-1.0, xmax=6.0 )
413 igroup.defineHistogram(namer('clusterSecondLambda','secondLambda','Cluster',postfix), title='cluster SECOND_LAMBDA;cluster SECOND_LAMBDA;Entries',path=folder+"Cluster",
414 xbins=50, xmin=2.0, xmax=6.5 )
415 igroup.defineHistogram(namer('clusterCenterLambda','centerLambda','Cluster',postfix), title='cluster CENTER_LAMBDA;cluster CENTER_LAMBDA;Entries',path=folder+"Cluster",
416 xbins=50, xmin=1.0, xmax=4.0 )
417 igroup.defineHistogram(namer('clusterFirstEngDens','firstEngDens','Cluster',postfix), title='cluster FIRST_ENG_DENS;cluster FIRST_ENG_DENS;Entries',path=folder+"Cluster",
418 xbins=50, xmin=-6.0, xmax=1.0 )
419 igroup.defineHistogram(namer('clusterEMproba','EMproba','Cluster',postfix), title='cluster EM_PROBABILITY;cluster EM_PROBABILITY;Entries',path=folder+"Cluster",
420 xbins=20, xmin=0.0, xmax=1.0 )
421
422 igroup.defineHistogram(namer('clustersMeanCenterLambda','meanCenterLambda','Cluster',postfix), title='clustersMeanCenterLambda;clustersMeanCenterLambda;Entries',path=folder+"Cluster",
423 xbins=40, xmin=0., xmax=2500.0 )
424 igroup.defineHistogram(namer('clustersMeanFirstEngDens','meanFirstEngDens','Cluster',postfix), title='clustersMeanFirstEngDens;clustersMeanFirstEngDens;Entries',path=folder+"Cluster",
425 xbins=40, xmin=-8.7, xmax=-5.5 )
426 igroup.defineHistogram(namer('clustersMeanEMProbability','meanEMProbability','Cluster',postfix), title='clustersMeanEMProbability;clustersMeanEMProbability;Entries',path=folder+"Cluster",
427 xbins=20, xmin=0.0, xmax=1.0 )
428 igroup.defineHistogram(namer('clustersMeanSecondLambda','meanSecondLambda','Cluster',postfix), title='clustersMeanSecondLambda;clustersMeanSecondLambda;Entries',path=folder+"Cluster",
429 xbins=20, xmin=0., xmax=600000 )
430 igroup.defineHistogram(namer('clustersMeanPresamplerFrac','meanPresamplerFrac','Cluster',postfix), title='clustersMeanPresamplerFrac;clustersMeanPresamplerFrac;Entries',path=folder+"Cluster",
431 xbins=20, xmin=0., xmax=0.2 )
432 igroup.defineHistogram(namer('clustersPFOEngRelDiff','PFOEngRelDiff','Cluster',postfix), title='PFOEngRelDiff;PFOEngRelDiff;Entries',path=folder+"Cluster",
433 xbins=20, xmin=-1.0, xmax=1.0 )
434 #Calo Histograms
435 igroup.defineHistogram(namer('EMRadius','EMRadius','Calo',postfix), title='Uncalibrated EM Radius;EM Radius;Number Of Candidates',path=folder+"Calo",
436 xbins=50, xmin=0., xmax=1.2 )
437
438 igroup.defineHistogram(namer('hadRadius','hadRadius','Calo',postfix), title='Hadronic Radius of tau candidates;Hadronic Radius; Number Of Candidates',path=folder+"Calo",
439 xbins=50, xmin=0., xmax=1. )
440
441 igroup.defineHistogram(namer('isolFrac','isolFrac','Calo',postfix), title='Isolation Fraction;Et Isolation Fraction;Number of Candidates',path=folder+"Calo",
442 xbins=51, xmin=0.0, xmax=1.02 )
443
444 igroup.defineHistogram(namer('etEMAtEMScale','etEMAtEMScale','Calo',postfix), title='EM energy at the EM scale;EM Et (GeV) );Number of Candidates',path=folder+"Calo",
445 xbins=50, xmin=0., xmax=200. )
446
447 igroup.defineHistogram(namer('etHadAtEMScale','etHadAtEMScale','Calo',postfix), title='Hadronic Energy at the EM Scale;Had Et (GeV) );Number of Candidates',path=folder+"Calo",
448 xbins=50, xmin=0., xmax=200. )
449
450 igroup.defineHistogram(namer('centFrac,LB','centFracVsLB','Calo',postfix), type='TH2F', title='Centrality Fraction vs Lumiblock;Centrality Fraction;Lumiblock', path=folder+"Calo",
451 xbins=51,xmin=0,xmax=1.02,ybins=1200,ymin=0.,ymax=1200.)
452
453 igroup.defineHistogram(namer('isolFrac,LB','isolFracVsLB','Calo',postfix), type='TH2F', title='Isolation Fraction vs Lumiblock;Isolation Fraction;Lumiblock', path=folder+"Calo",
454 xbins=51,xmin=0,xmax=1.02,ybins=1200,ymin=0.,ymax=1200.)
455
456 igroup.defineHistogram(namer('centFrac',"centralityFraction","Calo",postfix), title='Centrality Fraction;Centrality Fraction;Number Of Candidates',
457 xbins=70, xmin=0., xmax=1.4, path=folder+"Calo")
458
459 igroup.defineHistogram(namer('dRmax',"dRmax","Calo",postfix), title='Maximum delta R;delta R;Number of Candidates',
460 xbins=42, xmin=0., xmax=0.21, path=folder+"Calo" )
461
462 igroup.defineHistogram(namer('tauEta','eta','Calo',postfix), title='Eta of tau candidates;Eta;Number of Candidates',
463 xbins=51, xmin=-2.55, xmax=2.55,path=folder+"Calo")
464
465 igroup.defineHistogram(namer('tauPhi','phi','Calo',postfix), title='Phi of tau candidates;Phi;Number of Candidates',
466 xbins=65, xmin=PHIMIN, xmax=PHIMAX,path=folder+"Calo")
467
468 igroup.defineHistogram(namer('trkAvgDist','trkAvgDist','Calo',postfix), title='Average Track Distance from Calorimeter Seed;Distance (mm);Number of Candidates',
469 xbins=20, xmin=0., xmax=1., path=folder+"Calo" )
470
471 #Track Histograms!
472
473 #track variables used for track classification and RNN tau ID
474 igroup.defineHistogram(namer('trackLogSeedJetPt','logSeedJetPt','Track_allTracks',postfix), title='log(ptJetSeed);log(ptJetSeed);Entries',path=folder+"Track/allTracks",
475 xbins=50, xmin=10.0, xmax=14.5 )
476 igroup.defineHistogram(namer('trackLogPt','logPt','Track_allTracks',postfix), title='track log(p_{T});track log(p_{T});Entries',path=folder+"Track/allTracks",
477 xbins=50, xmin=6.0, xmax=15.0 )
478 igroup.defineHistogram(namer('trackEta','eta','Track_allTracks',postfix), title='track eta;track eta;Entries',path=folder+"Track/allTracks",
479 xbins=50, xmin=-2.6, xmax=2.6 )
480 igroup.defineHistogram(namer('trackd0TJVA','d0TJVA','Track_allTracks',postfix), title='track tanh(d0TJVA);track tanh(d0TJVA);Entries',path=folder+"Track/allTracks",
481 xbins=50, xmin=-0.1, xmax=0.1 )
482 igroup.defineHistogram(namer('trackZ0SinthetaTJVA','Z0SinthetaTJVA','Track_allTracks',postfix), title='track z0SinthetaTJVA;track z0SinthetaTJVA;Entries',path=folder+"Track/allTracks",
483 xbins=50, xmin=-100.0, xmax=100.0 )
484 igroup.defineHistogram(namer('trackD0SigTJVA','d0SigTJVA','Track_allTracks',postfix), title='track d0SigTJVA;track d0SigTJVA;Entries',path=folder+"Track/allTracks",
485 xbins=50, xmin=-10.0, xmax=10.0 )
486 igroup.defineHistogram(namer('trackZ0sinthetaSigTJVA','Z0sinthetaSigTJVA','Track_allTracks',postfix), title='track z0sinthetaSigTJVA;track z0sinthetaSigTJVA;Entries',path=folder+"Track/allTracks",
487 xbins=50, xmin=-10.0, xmax=10.0 )
488 igroup.defineHistogram(namer('trackCharge','charge','Track_allTracks',postfix), title='track charge;track charge;Entries',path=folder+"Track/allTracks",
489 xbins=8, xmin=-4.0, xmax=4.0 )
490 igroup.defineHistogram(namer('trackqOverP','qOverP','Track_allTracks',postfix), title='track qOverP;track qOverP;Entries',path=folder+"Track/allTracks",
491 xbins=50, xmin=-2.0, xmax=2.0 )
492 igroup.defineHistogram(namer('trackLogRConv','logRConv','Track_allTracks',postfix), title='track log(rConv);track log(rConv);Entries',path=folder+"Track/allTracks",
493 xbins=50, xmin=0.0, xmax=7.0 )
494 igroup.defineHistogram(namer('trackTanhRConvII','tanhRConvII','Track_allTracks',postfix), title='track tanh(rConvII/500);track tanh(rConvII/500);Entries',path=folder+"Track/allTracks",
495 xbins=50, xmin=-0.4, xmax=0.4 )
496 igroup.defineHistogram(namer('trackPtRatioSeedJet','ptRatioSeedJet','Track_allTracks',postfix), title='(track p_{T})/(ptJetSeed);(track p_{T})/(ptJetSeed);Entries',path=folder+"Track/allTracks",
497 xbins=50, xmin=0.0, xmax=0.1 )
498 igroup.defineHistogram(namer('trackdRJetSeedAxis','dRJetSeedAxis','Track_allTracks',postfix), title='#DeltaR(tau,track);#DeltaR(tau,track);Entries',path=folder+"Track/allTracks",
499 xbins=50, xmin=0.0, xmax=0.5 )
500 igroup.defineHistogram(namer('trackNInnermostPixHits','nInnermostPixHits','Track_allTracks',postfix), title='track nInnermostPixHits;track nInnermostPixHits;Entries',path=folder+"Track/allTracks",
501 xbins=5, xmin=0, xmax=5.0 )
502 igroup.defineHistogram(namer('trackNPixHits','nPixHits','Track_allTracks',postfix), title='track nPixHits;track nPixHits;Entries',path=folder+"Track/allTracks",
503 xbins=10, xmin=0.0, xmax=10.0 )
504 #nPixelSharedHits;track
505 igroup.defineHistogram(namer('trackNSiHits','nSiHits','Track_allTracks',postfix), title='track nSiHits;track nSiHits;Entries',path=folder+"Track/allTracks",
506 xbins=20, xmin=0.0, xmax=20.0 )
507 #nSCTSharedHits
508 #nTRTHits
509 igroup.defineHistogram(namer('trackeProbabilityHT','eProbabilityHT','Track_allTracks',postfix), title='track eProbabilityHT;track eProbabilityHT;Entries',path=folder+"Track/allTracks",
510 xbins=50, xmin=0.0, xmax=1.0 )
511 igroup.defineHistogram(namer('trackeProbabilityNN','eProbabilityNN','Track_allTracks',postfix), title='track eProbabilityNN;track eProbabilityNN;Entries',path=folder+"Track/allTracks",
512 xbins=50, xmin=0.0, xmax=1.0 )
513 igroup.defineHistogram(namer('trackeProbabilityHTorNN','eProbabilityHTorNN','Track_allTracks',postfix), title='track eProbabilityHTorNN;track eProbabilityHTorNN;Entries',path=folder+"Track/allTracks",
514 xbins=50, xmin=0.0, xmax=1.0 )
515
516 igroup.defineHistogram(namer('trackIdScoreCharged','idScoreCharged','Track_allTracks',postfix), title='track RNN idScoreCharged;track idScoreCharged;Entries',path=folder+"Track/allTracks",
517 xbins=50, xmin=0.0, xmax=1.0 )
518 igroup.defineHistogram(namer('trackIdScoreIso','idScoreIso','Track_allTracks',postfix), title='track RNN idScoreIso;track idScoreIso;Entries',path=folder+"Track/allTracks",
519 xbins=50, xmin=0.0, xmax=1.0 )
520 igroup.defineHistogram(namer('trackIdScoreConv','idScoreConv','Track_allTracks',postfix), title='track RNN idScoreConv;track idScoreConv;Entries',path=folder+"Track/allTracks",
521 xbins=50, xmin=0.0, xmax=1.0 )
522 igroup.defineHistogram(namer('trackIdScoreFake','idScoreFake','Track_allTracks',postfix), title='track RNN idScoreFake;track idScoreFake;Entries',path=folder+"Track/allTracks",
523 xbins=50, xmin=0.0, xmax=1.0 )
524
525 #original Track Histograms
526 igroup.defineHistogram(namer('d0','d0','Track_leadTrack',postfix), title='Track d0;Transverse Impact Parameter (mm)',path=folder+"Track/leadTrack",
527 xbins=50, xmin=-5.0, xmax=5.0 )
528
529 igroup.defineHistogram(namer('etaTrack','eta','Track_leadTrack',postfix), title='Track Eta;Eta',path=folder+"Track/leadTrack",
530 xbins=51, xmin=-2.55, xmax=2.55 )
531
532 igroup.defineHistogram(namer('leadTrkPt','leadTrkPt','Track_leadTrack',postfix), title='Pt of Leading track;Transverse Momentum (GeV)',path=folder+"Track/leadTrack",
533 xbins=50, xmin=0., xmax=200. )
534
535 igroup.defineHistogram(namer('nHighPtTaus','nHighPtTaus','Track',postfix), title='Number of High Pt tau candidates;Number of Taus per Event',
536 xbins=15, xmin=-0.5, xmax=14.5,path=folder+"Track")
537
538 igroup.defineHistogram(namer('numberOfTRTHighThresholdHits','numberOfTRTHighThresholdHits','Track_leadTrack',postfix), title='Number of TRT High Threshold Hits;Number of High Threshold TRT Hits',path=folder+"Track/leadTrack",
539 xbins=15, xmin=-0.5, xmax=14.5 )
540
541 igroup.defineHistogram(namer('numberOfTRTHighThresholdOutliers','numberOfTRTHighThresholdOutliers','Track_leadTrack',postfix), title="Number of TRT High Threshold Outliers;Number of TRT High Threshold Outliers",path=folder+"Track/leadTrack",
542 xbins=26, xmin=-0.5, xmax=25.5 )
543
544 igroup.defineHistogram(namer('numberOfTRTHits','numberOfTRTHits','Track_leadTrack',postfix), title='Number of TRT Low Threshold Hits;Number of Low Threshold TRT Hits',path=folder+"Track/leadTrack",
545 xbins=101, xmin=-0.5, xmax=100.5 )
546
547 igroup.defineHistogram(namer('numberOfTRTOutliers','numberOfTRTOutliers','Track_leadTrack',postfix), title="Number of TRT Low Threshold Outliers;Number of TRT Low Threshold Outliers",path=folder+"Track/leadTrack",
548 xbins=31, xmin=-0.5, xmax=25.5)
549
550 igroup.defineHistogram(namer('ptTrack','pT','Track_leadTrack',postfix), title='Track pT;Transverse Momentum (GeV)',path=folder+"Track/leadTrack",
551 xbins=60, xmin=0, xmax=300 )
552
553 igroup.defineHistogram(namer('phiTrack','phi','Track_leadTrack',postfix), title='Track Phi;Phi',path=folder+"Track/leadTrack",
554 xbins=64, xmin=PHIMIN, xmax=PHIMAX )
555
556 igroup.defineHistogram(namer('ipZ0SinThetaSigLeadTrk','ipZ0SinThetaSigLeadTrk','Track_leadTrack',postfix), title='Impact Parameter z0 Sine Theta Significance of Leading Track;Z0SinTheta Significance',path=folder+"Track/leadTrack",
557 xbins=50, xmin=-10.0, xmax=10.0 )
558
559 igroup.defineHistogram(namer('numberOfPixelHits','numberOfPixelHits','Track_leadTrack',postfix), title='Number of Pixel Hits;Number of Pixel Hits',path=folder+"Track/leadTrack",
560 xbins=25, xmin=-0.5, xmax=25.5 )
561
562 igroup.defineHistogram(namer('numberOfSCTHits','numberOfSCTHits','Track_leadTrack',postfix), title='Number of SCT Hits;Number of SCT Hits',path=folder+"Track/leadTrack",
563 xbins=26, xmin=0, xmax=26 )
564
565 igroup.defineHistogram(namer('etOverPtLeadTrack','etOverPtLeadTrack','Track_leadTrack',postfix), title='Et over Pt of lead track of tau candidates;Et/Pt;Number of Candidates',
566 xbins=50, xmin=0., xmax=12., path=folder+"Track/leadTrack" )
567
568 igroup.defineHistogram(namer('ipSigLeadTrk','ipSigLeadTrk','Track_leadTrack',postfix), title='Impact Parameter Significance of Leading Track;Transverse Impact Parameter Significance',
569 xbins=100, xmin=-5., xmax=5., path=folder+"Track/leadTrack" )
570
571 igroup.defineHistogram(namer('massTrkSys','massTrkSys','Track',postfix), title='Mass of the Track System;Invariant Mass (GeV)',
572 xbins=30, xmin=-1., xmax=5., path=folder+"Track" )
573
574 igroup.defineHistogram(namer('trFlightPathSig','trFlightPathSig','Track',postfix),title='Track Transverse Flight Path Significance;Transverse Flight Path Significance',
575 xbins=125, xmin=-4.5, xmax=8., path=folder+"Track" )
576
577 igroup.defineHistogram(namer('z0','z0','Track_leadTrack',postfix), title='Track z0;Longitudinal Impact Parameter (mm)',path=folder+"Track/leadTrack",
578 xbins=50, xmin=-35.0, xmax=35.0 )
579
580
581
582
583 if postfix == 'Global' or postfix.startswith('TauTrig') or postfix == "EleTrig" or postfix=="JetTrig":
584 igroup.defineHistogram(namer('LB','nTauPerLB','',postfix), title='Total number of tau candidates per LB;Luminosity Block);Number of Candidates',
585 xbins=1000, xmin=0., xmax=1000. ,path=folder)
586
587 igroup.defineHistogram(namer('nHighPtTauCandidates','nHightPtTauCandidates','',postfix), title='Number of High Pt tau candidates;Number of Taus per Event',
588 xbins=15, xmin=-0.5, xmax=14.5,path=folder)
589
590 igroup.defineHistogram(namer('NumTracks','tauNumTracks','',postfix), title='Number Of Tracks for Tau Candidates;Number Of Tracks;Number Of Candidates',
591 xbins=21, xmin=-0.5, xmax=20.5,path=folder)
592
593 igroup.defineHistogram(namer('nClusters','nCluster','',postfix), title='Number Of CaloTopoClusters;Number Of Clusters;Number Of Candidates',
594 xbins=40, xmin=0., xmax=40. ,path=folder)
595
596 igroup.defineHistogram(namer('tauEta,tauEt','tauEtVsEta','',postfix), type='TH2F', title='Tau Et Vs Eta;#eta;Transverse Energy (Gev)',
597 xbins=40,xmin=-2.55,xmax=2.55 ,ybins=300,ymin=0,ymax=300,path=folder)
598
599 igroup.defineHistogram(namer('tauPhi,tauEt','tauEtVsPhi','',postfix), type='TH2F', title='EtVsPhi;Phi;Transverse Energy (GeV);',
600 xbins=80,xmin=PHIMIN,xmax=PHIMAX,ybins=300,ymin=-0.,ymax=300.,path=folder)
601
602
603
608
609 accumulator = cfgHelper.result()
610 cfg.merge(accumulator)
611 return cfg
612
613 # # Otherwise, merge with result object and return
614 # acc = cfgHelper.result()
615 # result.merge(acc)
616 # return result
617
if(febId1==febId2)