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