ATLAS Offline Software
Functions
TrigMETMonitorAlgorithm Namespace Reference

Functions

def TrigMETMonConfig (inputFlags)
 

Function Documentation

◆ TrigMETMonConfig()

def TrigMETMonitorAlgorithm.TrigMETMonConfig (   inputFlags) 
Function to configures some algorithms in the monitoring system.

Definition at line 13 of file TrigMETMonitorAlgorithm.py.

13 def TrigMETMonConfig(inputFlags):
14  '''Function to configures some algorithms in the monitoring system.'''
15 
16  # The following class will make a sequence, configure algorithms, and link
17  # them to GenericMonitoringTools
18  from AthenaMonitoring import AthMonitorCfgHelper
19  helper = AthMonitorCfgHelper(inputFlags,'TrigMETAthMonitorCfg')
20 
21 
28 
29  from AthenaConfiguration.ComponentFactory import CompFactory
30  TrigMETMonAlg = helper.addAlgorithm(CompFactory.TrigMETMonitorAlgorithm,'TrigMETMonAlg')
31 
32  # You can actually make multiple instances of the same algorithm and give
33  # them different configurations
34  TrigMETMonChain1Alg = helper.addAlgorithm(CompFactory.TrigMETMonitorAlgorithm,'TrigMETMonChain1Alg')
35  TrigMETMonChain2Alg = helper.addAlgorithm(CompFactory.TrigMETMonitorAlgorithm,'TrigMETMonChain2Alg')
36  TrigMETMonChain3Alg = helper.addAlgorithm(CompFactory.TrigMETMonitorAlgorithm,'TrigMETMonChain3Alg')
37 
38  # # If for some really obscure reason you need to instantiate an algorithm
39  # # yourself, the AddAlgorithm method will still configure the base
40  # # properties and add the algorithm to the monitoring sequence.
41  # helper.AddAlgorithm(myExistingAlg)
42 
43 
44 
45  mt_chains = True
46  if ( inputFlags.Trigger.EDMVersion < 3 ) :
47  mt_chains = False
48 
49 
50 
57  TrigMETMonChain1 = 'HLT_xe80_cell_xe115_tcpufit_L1jXE100'
58  TrigMETMonChain2 = 'HLT_xe65_cell_xe90_pfopufit_L1jXE100'
59  TrigMETMonChain3 = 'HLT_xe65_cell_xe105_nn_L1jXE100'
60  if mt_chains:
61  TrigMETMonChain1Alg.TriggerChain = TrigMETMonChain1
62  TrigMETMonChain2Alg.TriggerChain = TrigMETMonChain2
63  TrigMETMonChain3Alg.TriggerChain = TrigMETMonChain3
64  else:
65  TrigMETMonChain1Alg.TriggerChain = 'HLT_xe110_pufit_xe65_L1XE50'
66  TrigMETMonChain2Alg.TriggerChain = 'HLT_xe110_pufit_xe65_L1XE50'
67  TrigMETMonChain3Alg.TriggerChain = 'HLT_xe110_pufit_xe65_L1XE50'
68 
69 
70 
71  from TrigConfigSvc.TriggerConfigAccess import getHLTMonitoringAccess
72  moniAccess=getHLTMonitoringAccess(inputFlags)
73  metChains=moniAccess.monitoredChains(signatures="metMon",monLevels=["shifter"])
74  metChainsVal=moniAccess.monitoredChains(signatures="metMon",monLevels=["val"])
75  metChainsT0=moniAccess.monitoredChains(signatures="metMon",monLevels=["t0"])
76 
77 
78  if mt_chains: # these are temporary, needs to be changed
79  TrigMETMonAlg.hlt_electron_key = 'HLT_egamma_Electrons'
80  TrigMETMonAlg.hlt_muon_key = 'HLT_MuonsCB_RoI'
81  TrigMETMonAlg.offline_met_key = 'MET_Reference_AntiKt4EMTopo' #this used to be 'MET_EMTopo'
82  TrigMETMonAlg.hlt_pfsum_key = 'HLT_MET_pfsum'
83  TrigMETMonAlg.l1_jFexMet_key = 'L1_jFexMETRoI'
84  TrigMETMonAlg.l1_jFexSumEt_key = 'L1_jFexSumETRoI'
85  TrigMETMonAlg.l1_gFexJwojScalar_key = 'L1_gScalarEJwoj'
86  TrigMETMonAlg.l1_gFexJwojMETComponents_key = 'L1_gMETComponentsJwoj'
87  TrigMETMonAlg.l1_gFexJwojMHTComponents_key = 'L1_gMHTComponentsJwoj'
88  TrigMETMonAlg.l1_gFexJwojMSTComponents_key = 'L1_gMSTComponentsJwoj'
89  TrigMETMonAlg.l1_gFexNCMETScalar_key = 'L1_gScalarENoiseCut'
90  TrigMETMonAlg.l1_gFexNCMETComponents_key = 'L1_gMETComponentsNoiseCut'
91  TrigMETMonAlg.l1_gFexRhoMETScalar_key = 'L1_gScalarERms'
92  TrigMETMonAlg.l1_gFexRhoMETComponents_key = 'L1_gMETComponentsRms'
93  else:
94  TrigMETMonAlg.hlt_electron_key = 'HLT_xAOD__ElectronContainer_egamma_Electrons'
95  TrigMETMonAlg.hlt_muon_key = 'HLT_xAOD__MuonContainer_MuonEFInfo'
96  TrigMETMonAlg.offline_met_key = 'MET_Reference_AntiKt4LCTopo'
97  TrigMETMonAlg.hlt_cell_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET'
98  TrigMETMonAlg.hlt_mht_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_mht'
99  TrigMETMonAlg.hlt_tc_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
100  TrigMETMonAlg.hlt_tc_em_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
101  TrigMETMonAlg.hlt_tcpufit_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl_PUC'
102  TrigMETMonAlg.hlt_trkmht_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
103  TrigMETMonAlg.hlt_pfsum_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
104  TrigMETMonAlg.hlt_pfsum_vssk_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
105  TrigMETMonAlg.hlt_pfsum_cssk_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
106  TrigMETMonAlg.hlt_pfopufit_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
107  TrigMETMonAlg.hlt_cvfpufit_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
108  TrigMETMonAlg.hlt_mhtpufit_pf_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
109  TrigMETMonAlg.hlt_mhtpufit_em_key = 'HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl'
110 
111 
112  L1Chains = ["L1_XE50",
113  "L1_XE55",
114  "L1_jXE100",
115  "L1_gXENC100",
116  "L1_gXERHO100",
117  "L1_gXEJWOJ100"]
118  HLTChains = []
119  HLTChainsVal = []
120  HLTChainsT0 = []
121  if mt_chains == 0:
122  L1Chains = ["L1_XE50"]
123  HLTChains = ["HLT_xe70_mht",
124  "HLT_xe90_mht_L1XE50",
125  "HLT_xe100_mht_L1XE50",
126  "HLT_xe110_mht_L1XE50",
127  "HLT_xe90_pufit_L1XE50",
128  "HLT_xe100_pufit_L1XE50",
129  "HLT_xe100_pufit_L1XE55",
130  "HLT_xe110_pufit_L1XE50",
131  "HLT_xe110_pufit_L1XE55",
132  "HLT_xe110_pufit_xe65_L1XE50",
133  "HLT_xe110_pufit_xe65_L1XE60",
134  "HLT_xe110_pufit_xe70_L1XE50"]
135 
136  if len(metChains) > 0:
137  HLTChains = metChains
138  if len(metChainsVal) > 0:
139  HLTChainsVal = metChainsVal
140  if len(metChainsT0) > 0:
141  HLTChainsT0 = metChainsT0
142 
143 
144  HLTChainEl = ["HLT_e26_lhtight_ivarloose_L1.*","HLT_e28_lhtight_ivarloose_L1.*"]
145  HLTChainMu = ["HLT_mu24_ivarmedium_L1.*"]
146 
147  TrigMETMonAlg.L1Chains = L1Chains
148  TrigMETMonAlg.HLTChains = HLTChains
149  TrigMETMonAlg.HLTChainsVal = HLTChainsVal
150  TrigMETMonAlg.HLTChainsT0 = HLTChainsT0
151  TrigMETMonAlg.HLTChainEl = HLTChainEl
152  TrigMETMonAlg.HLTChainMu = HLTChainMu
153 
154 
155  algsL1 = ["roi"]
156  algsL1Fex = ["jFex",
157  "gFexJwoj",
158  "gFexNC",
159  "gFexRho"]
160  algsHLT = ["cell",
161  "tcpufit",
162  "pfopufit",
163  "met_nn"]
164  algsHLTChain1 = ["cell",
165  "tcpufit"]
166  algsHLTChain2 = ["cell",
167  "pfopufit"]
168  algsHLTChain3 = ["cell",
169  "met_nn"]
170  algsHLTPreSel = ["cell",
171  "tcpufit",
172  "tcpufit_sig30",
173  "pfsum_cssk",
174  "pfsum_vssk",
175  "pfopufit",
176  "pfopufit_sig30",
177  "mhtpufit_pf",
178  "mhtpufit_em",
179  "met_nn"]
180  algsHLT2d = ["cell",
181  "tcpufit",
182  "pfopufit"]
183  algsHLTExpert = ["mht",
184  "tc_em",
185  "cvfpufit",
186  "pfsum",
187  "trkmht",
188  "mhtpufit_pf",
189  "mhtpufit_em",
190  "tcpufit_sig30",
191  "pfopufit_sig30",
192  "pfsum_cssk",
193  "pfsum_vssk"]
194  algsMET2d_tcpufit = ["pfopufit",
195  "pfsum_cssk",
196  "mhtpufit_pf"]
197 
198 
199  TrigMETMonAlg.algsL1 = algsL1
200  TrigMETMonAlg.algsHLT = algsHLT
201  TrigMETMonAlg.algsHLTPreSel = algsHLTPreSel
202  TrigMETMonAlg.algsHLT2d = algsHLT2d
203  TrigMETMonAlg.algsHLTExpert = algsHLTExpert
204  TrigMETMonAlg.algsMET2d_tcpufit = algsMET2d_tcpufit
205  TrigMETMonChain1Alg.algsHLT = algsHLTChain1
206  TrigMETMonChain2Alg.algsHLT = algsHLTChain2
207  TrigMETMonChain3Alg.algsHLT = algsHLTChain3
208 
209 
210  l1MetAlg = 1 #0=legacy,1=jFex,2=gFexJWOJ
211  l1MetCut = 48.0 #50.0,48.0,?
212  TrigMETMonAlg.L1MetAlg = l1MetAlg
213  TrigMETMonAlg.L1MetCut = l1MetCut
214 
215 
216  comp_names = ["PreSamplB", "EMB1", "EMB2", "EMB3", # LAr barrel
217  "PreSamplE", "EME1", "EME2", "EME3", # LAr EM endcap
218  "HEC0", "HEC1", "HEC2", "HEC3", # Hadronic end cap cal.
219  "TileBar0", "TileBar1", "TileBar2", # Tile barrel
220  "TileGap1", "TileGap2", "TileGap3", # Tile gap (ITC & scint)
221  "TileExt0", "TileExt1", "TileExt2", # Tile extended barrel
222  "FCalEM", "FCalHad2", "FCalHad3"] # Forward cal endcap
223  #"Muons" ] # Muons
224  bit_names = [
225  "Processing", # bit 0
226  "ErrBSconv", # bit 1
227  "ErrMuon", # bit 2
228  "ErrFEB", # bit 3
229  "Skipped", # bit 4
230  "CompBigMEtSEtRatio", # bit 5
231  "BadCompEnergy", # bit 6
232  "BadEnergyRatio", # bit 7
233  "spare", # bit 8
234  "BadCellQuality", # bit 9
235  "BadCellEnergy", # bit 10
236  "BadCellTime", # bit 11
237  "NoMuonTrack", # bit 12
238  "spare", # bit 13
239  "Processed", # bit 14
240  "CompError", # bit 15
241  "EMB_A_Missing", # bit 16
242  "EMB_C_Missing", # bit 17
243  "EME_A_Missing", # bit 18
244  "EME_C_Missing", # bit 19
245  "HEC_A_Missing", # bit 20
246  "HEC_C_Missing", # bit 21
247  "FCAL_A_Missing", # bit 22
248  "FCAL_C_Missing", # bit 23
249  "TileB_A_Missing", # bit 24
250  "TileB_C_Missing", # bit 25
251  "TileE_A_Missing", # bit 26
252  "TileE_C_Missing", # bit 27
253  "BadEMfraction", # bit 28
254  "GlobBigMEtSEtRatio", # bit 29
255  "ObjInCrack", # bit 30
256  "GlobError" # bit 31
257  ]
258 
259  TrigMETMonAlg.compNames = comp_names
260  TrigMETMonAlg.bitNames = bit_names
261 
262  electronPtCut = 30.0
263  muonPtCut = 30.0
264  electronEtaCut = 2.5
265  muonEtaCut = 2.5
266  LArNoiseBurstVetoAlgs = ["pfopufit",
267  "cell"]
268  signalLepAlgs = ["pfopufit",
269  "cell",
270  "tcpufit"]
271 
272  TrigMETMonAlg.electronPtCut = electronPtCut
273  TrigMETMonAlg.electronEtaCut = electronEtaCut
274  TrigMETMonAlg.muonPtCut = muonPtCut
275  TrigMETMonAlg.muonEtaCut = muonEtaCut
276  TrigMETMonAlg.signalLepAlgs = signalLepAlgs
277  TrigMETMonAlg.LArNoiseBurstVetoAlgs = LArNoiseBurstVetoAlgs
278 
279 
282 
283  # # Add a tool that doesn't have its own configuration function. In
284  # # this example, no accumulator is returned, so no merge is necessary.
285  # from MyDomainPackage.MyDomainPackageConf import MyDomainTool
286  # expertTrigMETMonAlg.MyDomainTool = MyDomainTool()
287 
288  # Add a generic monitoring tool (a "group" in old language). The returned
289  # object here is the standard GenericMonitoringTool.
290  metGroup = helper.addGroup(TrigMETMonAlg,'TrigMETMonitor','HLT/METMon/')
291 
292  # Add a GMT for the other example monitor algorithm
293  metChain1Group = helper.addGroup(TrigMETMonChain1Alg,'TrigMETMonitor','HLT/METMon/{}'.format(TrigMETMonChain1))
294  metChain2Group = helper.addGroup(TrigMETMonChain2Alg,'TrigMETMonitor','HLT/METMon/{}'.format(TrigMETMonChain2))
295  metChain3Group = helper.addGroup(TrigMETMonChain3Alg,'TrigMETMonitor','HLT/METMon/{}'.format(TrigMETMonChain3))
296 
297 
302  et_bins=205 #Et
303  et_min=-13.5
304  et_max=401.5
305  et_bins_log=20 #Et_log
306  et_min_log=-1.875
307  et_max_log=4.125
308  ec_bins=199 #Ex, Ey, Ez
309  ec_min=-298.5
310  ec_max=298.5
311  ec_bins_log=27 #Ex_log, Ey_log, Ez_log
312  ec_min_log=-4.125
313  ec_max_log=4.125
314  sumet_bins=305 #sumEt
315  sumet_min=-27.0
316  sumet_max=4203.0
317  sumet_bins_log=20 #sumEt_log
318  sumet_min_log=-1.875
319  sumet_max_log=4.125
320 
328  phi_bins=64 # phi, used to be 100
329  phi_bins_2d=32 # phi, used to be 24
330  phi_min=-3.1416
331  phi_max=3.1416
332  eta_bins_2d=24# eta
333  eta_min=-4.8
334  eta_max=4.8
335  eff_bins=42 # efficiency
336  eff_min=-13.5
337  eff_max=401.5
338  # Histograms
339 
340  metGroup.defineHistogram('offline_Ex',title='Offline Missing E_{x};E_{x} [GeV];Events',
341  path='Expert/Offline',xbins=ec_bins,xmin=ec_min,xmax=ec_max)
342  metGroup.defineHistogram('offline_Ey',title='Offline Missing E_{y};E_{y} [GeV];Events',
343  path='Expert/Offline',xbins=ec_bins,xmin=ec_min,xmax=ec_max)
344  metGroup.defineHistogram('offline_Et',title='Offline Missing E_{T};E_{T} [GeV];Events',
345  path='Expert/Offline',xbins=et_bins,xmin=et_min,xmax=et_max)
346  metGroup.defineHistogram('offline_sumEt',title='Offline sumE_{T};sumE_{T} [GeV];Events',
347  path='Expert/Offline',xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
348  metGroup.defineHistogram('offline_NoMu_Ex',title='Offline (No Mu) Missing E_{x};E_{x} [GeV];Events',
349  path='Expert/Offline',xbins=ec_bins,xmin=ec_min,xmax=ec_max)
350  metGroup.defineHistogram('offline_NoMu_Ey',title='Offline (No Mu) Missing E_{y};E_{y} [GeV];Events',
351  path='Expert/Offline',xbins=ec_bins,xmin=ec_min,xmax=ec_max)
352  metGroup.defineHistogram('offline_NoMu_Et',title='Offline (No Mu) Missing E_{T};E_{T} [GeV];Events',
353  path='Expert/Offline',xbins=et_bins,xmin=et_min,xmax=et_max)
354  metGroup.defineHistogram('offline_NoMu_sumEt',title='Offline (No Mu) sumE_{T};sumE_{T} [GeV];Events',
355  path='Expert/Offline',xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
356 
357  metGroup.defineHistogram('hlt_el_mult',title='HLT Electron Multiplicity;Number of electrons;Events',
358  path='Expert/ElMu',xbins=10,xmin=0,xmax=10)
359  metGroup.defineHistogram('hlt_el_pt',title='HLT Electron p_{T};p_{T} [GeV];Events',
360  path='Expert/ElMu',xbins=100,xmin=0,xmax=100)
361  metGroup.defineHistogram('hlt_mu_mult',title='HLT Muon Multiplicity;Number of muons;Events',
362  path='Expert/ElMu',xbins=10,xmin=0,xmax=10)
363  metGroup.defineHistogram('hlt_mu_pt',title='HLT Muon p_{T};p_{T} [GeV];Events',
364  path='Expert/ElMu',xbins=100,xmin=0,xmax=100)
365 
366  metGroup.defineHistogram('hlt_topoclusters_mult',title='HLT Topoclusters Multiplicity;Number of Clusters;Events',
367  path='Expert/Topoclusters',xbins=120,xmin=0,xmax=1200)
368  metGroup.defineHistogram('hlt_topoclusters_pt',title='HLT Topoclusters p_{T};p_{T} [GeV];Events',
369  path='Expert/Topoclusters',xbins=50,xmin=0,xmax=20)
370 
371  metGroup.defineHistogram('hlt_tracks_mult',title='HLT Tracks Multiplicity;Number of Tracks;Events',
372  path='Expert/Tracks',xbins=120,xmin=0,xmax=1200)
373  metGroup.defineHistogram('hlt_tracks_pt',title='HLT Tracks p_{T};p_{T} [GeV];Events',
374  path='Expert/Tracks',xbins=50,xmin=0,xmax=20)
375  metGroup.defineHistogram('hlt_tracks_leading_pt',title='HLT Tracks Leading p_{T};p_{T} [GeV];Events',
376  path='Expert/Tracks',xbins=50,xmin=0,xmax=20)
377  metGroup.defineHistogram('hlt_tracks_vec_sumPt',title='HLT Tracks Vector Sum p_{T};p_{T} [GeV];Events',
378  path='Expert/Tracks',xbins=100,xmin=0,xmax=100)
379  metGroup.defineHistogram('hlt_tracks_sca_sumPt',title='HLT Tracks Scalar Sum p_{T};p_{T} [GeV];Events',
380  path='Expert/Tracks',xbins=100,xmin=0,xmax=1000)
381  metGroup.defineHistogram('hlt_tracks_eta,hlt_tracks_phi;hlt_tracks_eta_phi',
382  type='TH2F',
383  title='HLT Tracks #eta - #phi (p_{T} > 3 GeV);#eta;#phi',
384  path='Expert/Tracks',
385  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
386 
387  metGroup.defineHistogram('hlt_vertex_mult',title='HLT Vertex Multiplicity;Number of Vertexs;Events',
388  path='Expert/Vertex',xbins=55,xmin=-5,xmax=50)
389  metGroup.defineHistogram('hlt_vertex_z',title='HLT Vertex Z;Vertex Z [mm];Events',
390  path='Expert/Vertex',xbins=100,xmin=-200,xmax=200)
391  metGroup.defineHistogram('hlt_vertex_z_diff',title='(HLT-Offline) Vertex Z Diff;Vertex Z [mm];Events',
392  path='Expert/Vertex',xbins=100,xmin=-200,xmax=200)
393  metGroup.defineHistogram('hlt_vertex_mult_mu,act_IPBC;hlt_vertex_mult_mu',
394  type='TProfile',
395  title='Average Vertex Mult. per IPBC;Actual IPBC;Average Vertex Mult.',
396  path='Expert/Vertex',
397  xbins=55, xmin=-5, xmax=55)
398 
399 
400  for alg in algsL1:
401  metGroup.defineHistogram('L1_{}_Ex'.format(alg),
402  title='L1_{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
403  path='Shifter/L1_{}'.format(alg),
404  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
405  metGroup.defineHistogram('L1_{}_Ex_log'.format(alg),
406  title='L1_{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
407  path='Shifter/L1_{}'.format(alg),
408  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
409  metGroup.defineHistogram('L1_{}_Ey'.format(alg),
410  title='L1_{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
411  path='Shifter/L1_{}'.format(alg),
412  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
413  metGroup.defineHistogram('L1_{}_Ey_log'.format(alg),
414  title='L1_{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
415  path='Shifter/L1_{}'.format(alg),
416  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
417  metGroup.defineHistogram('L1_{}_Et'.format(alg),
418  title='L1_{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
419  path='Shifter/L1_{}'.format(alg),
420  xbins=et_bins,xmin=et_min,xmax=et_max)
421  metGroup.defineHistogram('L1_{}_Et_log'.format(alg),
422  title='L1_{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
423  path='Shifter/L1_{}'.format(alg),
424  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
425  metGroup.defineHistogram('L1_{}_sumEt'.format(alg),
426  title='L1_{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
427  path='Shifter/L1_{}'.format(alg),
428  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
429  metGroup.defineHistogram('L1_{}_sumEt_log'.format(alg),
430  title='L1_{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
431  path='Shifter/L1_{}'.format(alg),
432  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
433  metGroup.defineHistogram('L1_{}_phi'.format(alg),
434  title='L1_{} #phi;#phi;Events'.format(alg),
435  path='Shifter/L1_{}'.format(alg),
436  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
437 
438  for alg in algsL1Fex:
439  metGroup.defineHistogram('L1_{}_Ex'.format(alg),
440  title='L1_{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
441  path='Shifter/L1_{}'.format(alg),
442  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
443  metGroup.defineHistogram('L1_{}_HT_Ex'.format(alg),
444  title='L1_{} HT Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
445  path='Shifter/L1_{}'.format(alg),
446  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
447  metGroup.defineHistogram('L1_{}_ST_Ex'.format(alg),
448  title='L1_{} ST Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
449  path='Shifter/L1_{}'.format(alg),
450  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
451  metGroup.defineHistogram('L1_{}_Ex_log'.format(alg),
452  title='L1_{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
453  path='Shifter/L1_{}'.format(alg),
454  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
455  metGroup.defineHistogram('L1_{}_HT_Ex_log'.format(alg),
456  title='L1_{} HT Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
457  path='Shifter/L1_{}'.format(alg),
458  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
459  metGroup.defineHistogram('L1_{}_ST_Ex_log'.format(alg),
460  title='L1_{} ST Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
461  path='Shifter/L1_{}'.format(alg),
462  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
463  metGroup.defineHistogram('L1_{}_Ey'.format(alg),
464  title='L1_{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
465  path='Shifter/L1_{}'.format(alg),
466  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
467  metGroup.defineHistogram('L1_{}_HT_Ey'.format(alg),
468  title='L1_{} HT Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
469  path='Shifter/L1_{}'.format(alg),
470  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
471  metGroup.defineHistogram('L1_{}_ST_Ey'.format(alg),
472  title='L1_{} ST Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
473  path='Shifter/L1_{}'.format(alg),
474  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
475  metGroup.defineHistogram('L1_{}_Ey_log'.format(alg),
476  title='L1_{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
477  path='Shifter/L1_{}'.format(alg),
478  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
479  metGroup.defineHistogram('L1_{}_HT_Ey_log'.format(alg),
480  title='L1_{} HT Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
481  path='Shifter/L1_{}'.format(alg),
482  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
483  metGroup.defineHistogram('L1_{}_ST_Ey_log'.format(alg),
484  title='L1_{} ST Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
485  path='Shifter/L1_{}'.format(alg),
486  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
487  metGroup.defineHistogram('L1_{}_Et'.format(alg),
488  title='L1_{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
489  path='Shifter/L1_{}'.format(alg),
490  xbins=et_bins,xmin=et_min,xmax=et_max)
491  metGroup.defineHistogram('L1_{}_Et_log'.format(alg),
492  title='L1_{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
493  path='Shifter/L1_{}'.format(alg),
494  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
495  metGroup.defineHistogram('L1_{}_sumEt'.format(alg),
496  title='L1_{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
497  path='Shifter/L1_{}'.format(alg),
498  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
499  metGroup.defineHistogram('L1_{}_sumEt_log'.format(alg),
500  title='L1_{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
501  path='Shifter/L1_{}'.format(alg),
502  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
503  metGroup.defineHistogram('L1_{}_phi'.format(alg),
504  title='L1_{} #phi;#phi;Events'.format(alg),
505  path='Shifter/L1_{}'.format(alg),
506  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
507  metGroup.defineHistogram('L1_{}_ST_phi'.format(alg),
508  title='L1_{} ST #phi;#phi;Events'.format(alg),
509  path='Shifter/L1_{}'.format(alg),
510  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
511  metGroup.defineHistogram('L1_{}_HT_phi'.format(alg),
512  title='L1_{} HT #phi;#phi;Events'.format(alg),
513  path='Shifter/L1_{}'.format(alg),
514  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
515 
516  for alg in algsHLT:
517  metGroup.defineHistogram('{}_Ex'.format(alg),
518  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
519  path='Shifter/HLT_{}'.format(alg),
520  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
521  metGroup.defineHistogram('{}_Ex_log'.format(alg),
522  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
523  path='Shifter/HLT_{}'.format(alg),
524  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
525  metGroup.defineHistogram('{}_Ey'.format(alg),
526  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
527  path='Shifter/HLT_{}'.format(alg),
528  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
529  metGroup.defineHistogram('{}_Ey_log'.format(alg),
530  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
531  path='Shifter/HLT_{}'.format(alg),
532  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
533  metGroup.defineHistogram('{}_Et'.format(alg),
534  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
535  path='Shifter/HLT_{}'.format(alg),
536  xbins=et_bins,xmin=et_min,xmax=et_max)
537  metGroup.defineHistogram('{}_Et_log'.format(alg),
538  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
539  path='Shifter/HLT_{}'.format(alg),
540  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
541  metGroup.defineHistogram('{}_sumEt'.format(alg),
542  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
543  path='Shifter/HLT_{}'.format(alg),
544  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
545  metGroup.defineHistogram('{}_sumEt_log'.format(alg),
546  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
547  path='Shifter/HLT_{}'.format(alg),
548  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
549  metGroup.defineHistogram('{}_phi'.format(alg),
550  title='{} #phi;#phi;Events'.format(alg),
551  path='Shifter/HLT_{}'.format(alg),
552  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
553  metGroup.defineHistogram('{0}_phi;{0}_phi_etweight'.format(alg),
554  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
555  weight='{}_Et'.format(alg),
556  path='Shifter/HLT_{}'.format(alg),
557  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
558 
559  for alg in algsHLTPreSel:
560  metGroup.defineHistogram('{}_presel_Et'.format(alg),
561  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
562  path='Shifter/preSel',
563  xbins=et_bins,xmin=et_min,xmax=et_max)
564 
565  for alg in LArNoiseBurstVetoAlgs:
566  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Ex'.format(alg),
567  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
568  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
569  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
570  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Ex_log'.format(alg),
571  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
572  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
573  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
574  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Ey'.format(alg),
575  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
576  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
577  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
578  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Ey_log'.format(alg),
579  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
580  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
581  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
582  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Et'.format(alg),
583  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
584  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
585  xbins=et_bins,xmin=et_min,xmax=et_max)
586  metGroup.defineHistogram('{}_LArNoiseBurstVeto_Et_log'.format(alg),
587  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
588  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
589  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
590  metGroup.defineHistogram('{}_LArNoiseBurstVeto_sumEt'.format(alg),
591  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
592  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
593  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
594  metGroup.defineHistogram('{}_LArNoiseBurstVeto_sumEt_log'.format(alg),
595  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
596  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
597  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
598  metGroup.defineHistogram('{}_LArNoiseBurstVeto_phi'.format(alg),
599  title='{} #phi;#phi;Events'.format(alg),
600  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
601  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
602  metGroup.defineHistogram('{0}_LArNoiseBurstVeto_phi;{0}_phi_etweight'.format(alg),
603  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
604  weight='{}_Et'.format(alg),
605  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
606  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
607  metGroup.defineHistogram('{0}_LArNoiseBurstVeto_eta,{0}_LArNoiseBurstVeto_phi;{0}_LArNoiseBurstVeto_eta_phi'.format(alg),
608  type='TH2F',
609  title='{} #eta - #phi;#eta;#phi'.format(alg),
610  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
611  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
612  metGroup.defineHistogram('{0}_LArNoiseBurstVeto_eta,{0}_LArNoiseBurstVeto_phi;{0}_LArNoiseBurstVeto_eta_phi_etweight'.format(alg),
613  type='TH2F',
614  title='{} #eta - #phi (etweighted);#eta;#phi'.format(alg),
615  weight='{}_LArNoiseBurstVeto_Et'.format(alg),
616  path='Shifter/LArNoiseBurstVetoed/{}'.format(alg),
617  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
618  # for alg in signalLepAlgs:
619  for alg in signalLepAlgs:
620  metGroup.defineHistogram('{}_SigEl_Ex'.format(alg),
621  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
622  path='Shifter/SignalEl/{}'.format(alg),
623  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
624  metGroup.defineHistogram('{}_SigEl_Ex_log'.format(alg),
625  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
626  path='Shifter/SignalEl/{}'.format(alg),
627  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
628  metGroup.defineHistogram('{}_SigEl_Ey'.format(alg),
629  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
630  path='Shifter/SignalEl/{}'.format(alg),
631  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
632  metGroup.defineHistogram('{}_SigEl_Ey_log'.format(alg),
633  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
634  path='Shifter/SignalEl/{}'.format(alg),
635  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
636  metGroup.defineHistogram('{}_SigEl_Et'.format(alg),
637  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
638  path='Shifter/SignalEl/{}'.format(alg),
639  xbins=et_bins,xmin=et_min,xmax=et_max)
640  metGroup.defineHistogram('{}_SigEl_Et_log'.format(alg),
641  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
642  path='Shifter/SignalEl/{}'.format(alg),
643  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
644  metGroup.defineHistogram('{}_SigEl_sumEt'.format(alg),
645  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
646  path='Shifter/SignalEl/{}'.format(alg),
647  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
648  metGroup.defineHistogram('{}_SigEl_sumEt_log'.format(alg),
649  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
650  path='Shifter/SignalEl/{}'.format(alg),
651  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
652  metGroup.defineHistogram('{}_SigEl_phi'.format(alg),
653  title='{} #phi;#phi;Events'.format(alg),
654  path='Shifter/SignalEl/{}'.format(alg),
655  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
656  metGroup.defineHistogram('{0}_SigEl_phi;{0}_SigEl_phi_etweight'.format(alg),
657  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
658  weight='{}_SigEl_Et'.format(alg),
659  path='Shifter/SignalEl/{}'.format(alg),
660  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
661  metGroup.defineHistogram('{0}_SigEl_eta,{0}_SigEl_phi;{0}_SigEl_eta_phi'.format(alg),
662  type='TH2F',
663  title='{} #eta - #phi;#eta;#phi'.format(alg),
664  path='Shifter/SignalEl/{}'.format(alg),
665  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
666  metGroup.defineHistogram('{0}_SigEl_eta,{0}_SigEl_phi;{0}_SigEl_eta_phi_etweight'.format(alg),
667  type='TH2F',
668  title='{} #eta - #phi (etweighted);#eta;#phi'.format(alg),
669  weight='{}_SigEl_Et'.format(alg),
670  path='Shifter/SignalEl/{}'.format(alg),
671  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
672  metGroup.defineHistogram('{}_SigMu_Ex'.format(alg),
673  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
674  path='Shifter/SignalMu/{}'.format(alg),
675  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
676  metGroup.defineHistogram('{}_SigMu_Ex_log'.format(alg),
677  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
678  path='Shifter/SignalMu/{}'.format(alg),
679  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
680  metGroup.defineHistogram('{}_SigMu_Ey'.format(alg),
681  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
682  path='Shifter/SignalMu/{}'.format(alg),
683  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
684  metGroup.defineHistogram('{}_SigMu_Ey_log'.format(alg),
685  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
686  path='Shifter/SignalMu/{}'.format(alg),
687  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
688  metGroup.defineHistogram('{}_SigMu_Et'.format(alg),
689  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
690  path='Shifter/SignalMu/{}'.format(alg),
691  xbins=et_bins,xmin=et_min,xmax=et_max)
692  metGroup.defineHistogram('{}_SigMu_Et_log'.format(alg),
693  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
694  path='Shifter/SignalMu/{}'.format(alg),
695  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
696  metGroup.defineHistogram('{}_SigMu_sumEt'.format(alg),
697  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
698  path='Shifter/SignalMu/{}'.format(alg),
699  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
700  metGroup.defineHistogram('{}_SigMu_sumEt_log'.format(alg),
701  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
702  path='Shifter/SignalMu/{}'.format(alg),
703  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
704  metGroup.defineHistogram('{}_SigMu_phi'.format(alg),
705  title='{} #phi;#phi;Events'.format(alg),
706  path='Shifter/SignalMu/{}'.format(alg),
707  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
708  metGroup.defineHistogram('{0}_SigMu_phi;{0}_SigMu_phi_etweight'.format(alg),
709  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
710  weight='{}_SigMu_Et'.format(alg),
711  path='Shifter/SignalMu/{}'.format(alg),
712  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
713  metGroup.defineHistogram('{0}_SigMu_eta,{0}_SigMu_phi;{0}_SigMu_eta_phi'.format(alg),
714  type='TH2F',
715  title='{} #eta - #phi;#eta;#phi'.format(alg),
716  path='Shifter/SignalMu/{}'.format(alg),
717  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
718  metGroup.defineHistogram('{0}_SigMu_eta,{0}_SigMu_phi;{0}_SigMu_eta_phi_etweight'.format(alg),
719  type='TH2F',
720  title='{} #eta - #phi (etweighted);#eta;#phi'.format(alg),
721  weight='{}_SigMu_Et'.format(alg),
722  path='Shifter/SignalMu/{}'.format(alg),
723  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
724 
725 
726  for alg in algsHLT2d:
727  metGroup.defineHistogram('{0}_eta,{0}_phi;{0}_eta_phi'.format(alg),
728  type='TH2F',
729  title='{} #eta - #phi;#eta;#phi'.format(alg),
730  path='Shifter/HLT_{}'.format(alg),
731  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
732  metGroup.defineHistogram('{a}_eta,{a}_phi;{a}_eta_phi_etweight'.format(a=alg),
733  type='TH2F',
734  title='{} #eta - #phi (etweighted);#eta;#phi'.format(alg),
735  weight='{}_Et'.format(alg),
736  path='Shifter/HLT_{}'.format(alg),
737  xbins=eta_bins_2d,xmin=eta_min,xmax=eta_max,ybins=phi_bins_2d,ymin=phi_min,ymax=phi_max)
738 
739 
740  for chain in L1Chains:
741  metGroup.defineHistogram('offline_NoMu_Et_eff,pass_{};{}_eff'.format(chain,chain),
742  type='TProfile',
743  title='{} efficiency;offline E_{{T}} [GeV];Efficiency'.format(chain),
744  path='Shifter/eff',
745  xbins=eff_bins, xmin=eff_min, xmax=eff_max)
746 
747  for chain in HLTChains:
748  metGroup.defineHistogram('offline_NoMu_Et_eff,pass_{};{}_eff'.format(chain,chain),
749  type='TProfile',
750  title='{} efficiency;offline E_{{T}} [GeV];Efficiency'.format(chain),
751  path='Shifter/eff',
752  xbins=eff_bins, xmin=eff_min, xmax=eff_max)
753 
754  for chain in HLTChainsVal:
755  metGroup.defineHistogram('offline_NoMu_Et_eff,pass_{};{}_eff'.format(chain,chain),
756  type='TProfile',
757  title='{} efficiency;offline E_{{T}} [GeV];Efficiency'.format(chain),
758  path='Expert/eff',
759  xbins=eff_bins, xmin=eff_min, xmax=eff_max)
760  for chain in HLTChainsT0:
761  metGroup.defineHistogram('offline_NoMu_Et_eff,pass_{};{}_eff'.format(chain,chain),
762  type='TProfile',
763  title='{} efficiency;offline E_{{T}} [GeV];Efficiency'.format(chain),
764  path='Expert/eff',
765  xbins=eff_bins, xmin=eff_min, xmax=eff_max)
766 
767 
768  metGroup.defineHistogram('act_IPBC', type='TH1F',title='Actual IPBC;Actual IPBC;Events',
769  path='Shifter/Component',xbins=100,xmin=0,xmax=100)
770 
771  metGroup.defineHistogram('HLT_MET_status',type='TH1F',title='HLT MET Status;;',
772  weight='MET_status',
773  path='Shifter/Component',
774  xbins=len(bit_names),xmin=-0.5,xmax=31.5, xlabels=bit_names)
775  metGroup.defineHistogram('HLT_MET_component,component_Et;compN_compEt',
776  type='TH2F',
777  title='HLT Missing E_{T} VS component;;Missing E_{T} [GeV]',
778  path='Shifter/Component',
779  xbins=len(comp_names),xmin=-0.5,xmax=24.5,ybins=et_bins,ymin=et_min,ymax=et_max,
780  xlabels=comp_names)
781  metGroup.defineHistogram('component,component_status;compN_HLT_MET_status',
782  type='TH2F',
783  title='HLT MET Status VS component;;',
784  weight='component_status_weight',
785  path='Shifter/Component',
786  xbins=len(comp_names),xmin=-0.5,xmax=24.5,ybins=len(comp_names),ymin=-0.5,ymax=31.5,
787  xlabels=comp_names, ylabels=bit_names)
788 
789  for alg in algsHLTExpert:
790  metGroup.defineHistogram('{}_Ex'.format(alg),
791  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
792  path='Expert/HLT_{}'.format(alg),
793  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
794  metGroup.defineHistogram('{}_Ey'.format(alg),
795  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
796  path='Expert/HLT_{}'.format(alg),
797  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
798  metGroup.defineHistogram('{}_Et'.format(alg),
799  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
800  path='Expert/HLT_{}'.format(alg),
801  xbins=et_bins,xmin=et_min,xmax=et_max)
802  metGroup.defineHistogram('{}_sumEt'.format(alg),
803  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
804  path='Expert/HLT_{}'.format(alg),
805  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
806  #2D MET tcpufit vs pfopufit, pfsum_cssk, trkmht
807  for alg in algsMET2d_tcpufit:
808  metGroup.defineHistogram('{}_2D_Et,tcpufit_2D_Et;hlt_tcpufit_Et_{}_Et'.format(alg, alg),
809  type='TH2F',
810  title='HLT tcpufit Missing E_{{T}} vs. HLT {} Missing Et;{} E_{{T}} [GeV];tcpufit E_{{T}} [GeV]'.format(alg, alg),
811  path='Expert/HLT_MET2D',
812  xbins=40,xmin=et_min,xmax=et_max,ybins=40,ymin=et_min,ymax=et_max)
813 
814  for alg in algsHLTChain1:
815  metChain1Group.defineHistogram('{}_Ex'.format(alg),
816  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
817  path='HLT_{}'.format(alg),
818  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
819  metChain1Group.defineHistogram('{}_Ex_log'.format(alg),
820  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
821  path='HLT_{}'.format(alg),
822  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
823  metChain1Group.defineHistogram('{}_Ey'.format(alg),
824  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
825  path='HLT_{}'.format(alg),
826  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
827  metChain1Group.defineHistogram('{}_Ey_log'.format(alg),
828  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
829  path='HLT_{}'.format(alg),
830  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
831  metChain1Group.defineHistogram('{}_Et'.format(alg),
832  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
833  path='HLT_{}'.format(alg),
834  xbins=et_bins,xmin=et_min,xmax=et_max)
835  metChain1Group.defineHistogram('{}_Et_log'.format(alg),
836  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
837  path='HLT_{}'.format(alg),
838  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
839  metChain1Group.defineHistogram('{}_sumEt'.format(alg),
840  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
841  path='HLT_{}'.format(alg),
842  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
843  metChain1Group.defineHistogram('{}_sumEt_log'.format(alg),
844  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
845  path='HLT_{}'.format(alg),
846  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
847  metChain1Group.defineHistogram('{}_phi'.format(alg),
848  title='{} #phi;#phi;Events'.format(alg),
849  path='HLT_{}'.format(alg),
850  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
851  metChain1Group.defineHistogram('{0}_phi;{0}_phi_etweight'.format(alg),
852  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
853  weight='{}_Et'.format(alg),
854  path='HLT_{}'.format(alg),
855  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
856 
857  for alg in algsHLTChain2:
858  metChain2Group.defineHistogram('{}_Ex'.format(alg),
859  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
860  path='HLT_{}'.format(alg),
861  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
862  metChain2Group.defineHistogram('{}_Ex_log'.format(alg),
863  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
864  path='HLT_{}'.format(alg),
865  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
866  metChain2Group.defineHistogram('{}_Ey'.format(alg),
867  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
868  path='HLT_{}'.format(alg),
869  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
870  metChain2Group.defineHistogram('{}_Ey_log'.format(alg),
871  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
872  path='HLT_{}'.format(alg),
873  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
874  metChain2Group.defineHistogram('{}_Et'.format(alg),
875  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
876  path='HLT_{}'.format(alg),
877  xbins=et_bins,xmin=et_min,xmax=et_max)
878  metChain2Group.defineHistogram('{}_Et_log'.format(alg),
879  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
880  path='HLT_{}'.format(alg),
881  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
882  metChain2Group.defineHistogram('{}_sumEt'.format(alg),
883  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
884  path='HLT_{}'.format(alg),
885  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
886  metChain2Group.defineHistogram('{}_sumEt_log'.format(alg),
887  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
888  path='HLT_{}'.format(alg),
889  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
890  metChain2Group.defineHistogram('{}_phi'.format(alg),
891  title='{} #phi;#phi;Events'.format(alg),
892  path='HLT_{}'.format(alg),
893  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
894  metChain2Group.defineHistogram('{0}_phi;{0}_phi_etweight'.format(alg),
895  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
896  weight='{}_Et'.format(alg),
897  path='HLT_{}'.format(alg),
898  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
899 
900  for alg in algsHLTChain3:
901  metChain3Group.defineHistogram('{}_Ex'.format(alg),
902  title='{} Missing E_{{x}};E_{{x}} [GeV];Events'.format(alg),
903  path='HLT_{}'.format(alg),
904  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
905  metChain3Group.defineHistogram('{}_Ex_log'.format(alg),
906  title='{} Missing E_{{x}} log;sgn(E_{{x}}) log(E_{{x}}/GeV);Events'.format(alg),
907  path='HLT_{}'.format(alg),
908  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
909  metChain3Group.defineHistogram('{}_Ey'.format(alg),
910  title='{} Missing E_{{y}};E_{{y}} [GeV];Events'.format(alg),
911  path='HLT_{}'.format(alg),
912  xbins=ec_bins,xmin=ec_min,xmax=ec_max)
913  metChain3Group.defineHistogram('{}_Ey_log'.format(alg),
914  title='{} Missing E_{{y}} log;sgn(E_{{y}}) log(E_{{y}}/GeV);Events'.format(alg),
915  path='HLT_{}'.format(alg),
916  xbins=ec_bins_log,xmin=ec_min_log,xmax=ec_max_log)
917  metChain3Group.defineHistogram('{}_Et'.format(alg),
918  title='{} Missing E_{{T}};E_{{T}} [GeV];Events'.format(alg),
919  path='HLT_{}'.format(alg),
920  xbins=et_bins,xmin=et_min,xmax=et_max)
921  metChain3Group.defineHistogram('{}_Et_log'.format(alg),
922  title='{} Missing E_{{T}} log;log(E_{{T}}/GeV);Events'.format(alg),
923  path='HLT_{}'.format(alg),
924  xbins=et_bins_log,xmin=et_min_log,xmax=et_max_log)
925  metChain3Group.defineHistogram('{}_sumEt'.format(alg),
926  title='{} sumE_{{T}};sumE_{{T}} [GeV];Events'.format(alg),
927  path='HLT_{}'.format(alg),
928  xbins=sumet_bins,xmin=sumet_min,xmax=sumet_max)
929  metChain3Group.defineHistogram('{}_sumEt_log'.format(alg),
930  title='{} sumE_{{T}} log;log(sumE_{{T}}/GeV);Events'.format(alg),
931  path='HLT_{}'.format(alg),
932  xbins=sumet_bins_log,xmin=sumet_min_log,xmax=sumet_max_log)
933  metChain3Group.defineHistogram('{}_phi'.format(alg),
934  title='{} #phi;#phi;Events'.format(alg),
935  path='HLT_{}'.format(alg),
936  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
937  metChain3Group.defineHistogram('{0}_phi;{0}_phi_etweight'.format(alg),
938  title='{} #phi (etweighted);#phi;Et weighted events'.format(alg),
939  weight='{}_Et'.format(alg),
940  path='HLT_{}'.format(alg),
941  xbins=phi_bins,xmin=phi_min,xmax=phi_max)
942 
943 
944 
949  return helper.result()
950 
vtune_athena.format
format
Definition: vtune_athena.py:14
python.TriggerConfigAccess.getHLTMonitoringAccess
HLTMonitoringAccess getHLTMonitoringAccess(flags=None)
Definition: TriggerConfigAccess.py:256
TrigMETMonitorAlgorithm.TrigMETMonConfig
def TrigMETMonConfig(inputFlags)
Definition: TrigMETMonitorAlgorithm.py:13
Generated on Thu Jul 4 2024 22:29:34 for ATLAS Offline Software by   doxygen 1.8.18