6 @file SetupEgammaMonitoring.py
9 @brief Example python configuration for the Run III AthenaMonitoring package
11 from AthenaConfiguration.ComponentFactory
import CompFactory
12 import TrigEgammaMonitoring.TrigEgammaMonitCategory
as egammaConf
16 Function to configure the egamma histograms per region (BARREL, CRACK, ENDCAP, FORWARD)
17 The list is different for electrons, forward electrons and photons
19 Regions = [
'BARREL',
'CRACK',
'ENDCAP']
20 for n
in range(len(Regions)):
21 hname = thename +
"in" + Regions[n]
22 htitle = title +
" " + theparttype +
" " + Regions[n]
23 thefinalcut = thecut+Regions[n]
24 thegroupe.defineHistogram(hname,title=htitle, path=path, xbins=xbins,xmin=xmin,xmax=xmax,type = thetype, cutmask = thefinalcut)
28 Function to configure the egamma histograms per region (BARREL, CRACK, ENDCAP, FORWARD)
29 The list is different for electrons, forward electrons and photons
31 Regions = [
'ENDCAP',
'FORWARD']
32 for n
in range(len(Regions)):
33 hname = thename +
"in" + Regions[n]
34 htitle = title +
" " + theparttype +
" " + Regions[n]
35 thefinalcut = thecut+Regions[n]
36 thegroupe.defineHistogram(hname,title=htitle, path=path, xbins=xbins,xmin=xmin,xmax=xmax,type = thetype, cutmask = thefinalcut)
40 Function to configure the egamma histograms
41 The list is different for electrons and photons
42 usage example : BookHistograms(GroupElectronLHTight,"LHTight","Electron","WithTrigger")
45 if particletype
in [
"Electron",
"Photon"] :
46 if particletype
in [
"Electron"] :
49 if particletype
in [
"Photon"] :
54 htitle =
"Number of " + prefix +
"s (" + reconame +
" " + withTrigger+
") ; N" + tlabel +
" ; N_{event}"
55 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
57 if particletype
in [
"Photon"] :
59 htitle =
"Number of converted photons (" + reconame +
" " + withTrigger+
") ; N_{Conv. #gamma} ; N_{event}"
60 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
63 htitle =
"Number of unconverted photons (" + reconame +
" " + withTrigger+
") ; N_{Unconv. #gamma} ; N_{event}"
64 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
67 htitle= particletype +
" transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Et"+ tlabel +
" ; N" + tlabel
68 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=2.,xmax=102000.0)
70 if particletype
in [
"Photon"] :
73 htitle=
"Converted #gamma transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Et" + tlabel +
" ; N" + tlabel
74 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=2.,xmax=102000.0, cutmask =
'is_pt_gt_4gevandconv')
77 htitle=
"Unconverted #gamma transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Et" + tlabel +
" ; N" + tlabel
78 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=2.,xmax=102000.0, cutmask =
'is_pt_gt_4gevandunconv')
81 htitle= particletype +
" #eta" +
" (" + reconame +
" " + withTrigger +
")"+
" ; #eta" + tlabel +
" ; N" + tlabel
82 groupe.defineHistogram(hname,title=htitle, path=
'',type =
"TH1D", xbins=64,xmin=-3.2,xmax=3.2)
84 if particletype
in [
"Photon"] :
86 htitle=
"Converted photons #eta" +
" (" + reconame +
" " + withTrigger +
")"+
" ; #eta" + tlabel +
" ; N" + tlabel
87 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2,cutmask =
'is_pt_gt_4gevandconv')
90 htitle=
"Unconverted photons #eta" +
" (" + reconame +
" " + withTrigger +
")"+
" ; #eta" + tlabel +
" ; N" + tlabel
91 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2,cutmask =
'is_pt_gt_4gevandunconv')
94 htitle= particletype +
" #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #phi" + tlabel +
" ; N" + tlabel
95 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2)
97 if particletype
in [
"Photon"] :
99 htitle=
"Converted photon #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #phi" + tlabel +
" ; N" + tlabel
100 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2,cutmask =
'is_pt_gt_4gevandconv')
103 htitle=
"Unconverted photon #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #phi" + tlabel +
" ; N" + tlabel
104 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2, cutmask =
'is_pt_gt_4gevandunconv')
106 hname=
"Eta,Phi;Eta_Phi_with_Pt_gt_2.5GeV"
107 htitle= particletype +
" #eta,#phi map (candidates with Pt>2.5GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
108 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask =
'is_pt_gt_2_5gev')
110 if particletype
in [
"Photon"] :
111 hname=
"Eta,Phi;Eta_Phi_Conv_with_Pt_gt_2.5GeV"
112 htitle=
"(#eta,#phi) map of Conv. #gamma with Pt>2.5GeV" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
113 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2, cutmask =
'is_pt_gt_2_5gevandconv')
115 hname=
"Eta,Phi;Eta_Phi_Unconv_with_Pt.gt.2.5GeV"
116 htitle=
"(#eta,#phi) map of Unconv. #gamma with Pt>2.5GeV" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
117 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2, cutmask =
'is_pt_gt_2_5gevandunconv')
119 hname=
"Eta,Phi;Eta_Phi_distribution_with_Pt.gt.4GeV"
120 htitle= particletype +
" #eta,#phi map (candidates with Pt>4GeV)" +
" (" + reconame +
" " + prefix +
" " + withTrigger +
") ; #eta ; #phi"
121 groupe.defineHistogram(hname,title=htitle,path=
'',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_4gev')
123 if particletype
in [
"Photon"] :
124 hname=
"Eta,Phi;Eta_Phi_Conv_with_Pt.gt.4GeV"
125 htitle= particletype +
" #eta,#phi map (candidates with Pt>4GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
126 groupe.defineHistogram(hname,title=htitle,path=
'',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_4gevandconv')
128 hname=
"Eta,Phi;Eta_Phi_Unconv_with_Pt.gt.4GeV"
129 htitle= particletype +
" #eta,#phi map (candidates with Pt>4GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
130 groupe.defineHistogram(hname,title=htitle,path=
'',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_4gevandunconv')
132 hname=
"Eta,Phi;Eta_Phi_with_Pt.gt.20GeV"
133 htitle= particletype +
" #eta,#phi map (candidates with Pt>20GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi "
134 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_20gev')
136 if particletype
in [
"Photon"] :
137 hname=
"Eta,Phi;Eta_Phi_distribution_of_Conv._photons_with_Pt.gt.20GeV"
138 htitle= particletype +
" #eta,#phi map (candidates with Pt>20GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
139 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_20gevandconv')
141 hname=
"Eta,Phi;Eta_Phi_distribution_of_Unconv._photons_with_Pt.gt.20GeV"
142 htitle= particletype +
" #eta,#phi map (candidates with Pt>20GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
143 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_20gevandunconv')
145 hname=
"TopoEtCone40"
146 htitle = particletype +
" Topocluster Isolation Energy" +
" (" + reconame +
" " + withTrigger +
") ; TopoEtCone40 ; N" + tlabel
147 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=-10000.,xmax=40000.)
150 htitle = particletype +
" Track Isolation Pt" +
" (" + reconame +
" " + withTrigger +
") ; PtCone20 ; N" + tlabel
151 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-10000.,xmax=40000.)
154 htitle =
"Time associated with " + particletype +
" cluster [ns]" +
" (" + reconame +
" " + withTrigger +
") ; Time [ns] ; N" + tlabel
155 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=90,xmin=-30.,xmax=60.)
157 if particletype
in [
"Photon"] :
159 htitle =
"Photon Conversion radius [mm] " +
" (" + reconame +
" " + withTrigger +
") ; Rconv ; N_{#gamma}"
160 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=0.,xmax=800.)
164 htitle= particletype +
" transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
") ; Et" + tlabel +
" ; N" + tlabel
165 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=100,xmin=2.,xmax=102000.0,thecut=
'is_pt_gt_2_5gev')
168 htitle= particletype +
" #eta" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #eta" + tlabel +
" ; N" + tlabel
169 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=64,xmin=-3.2,xmax=3.2)
172 htitle= particletype +
" #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #eta" + tlabel +
" ; N" + tlabel
173 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=64,xmin=-3.2,xmax=3.2)
176 htitle =
"Time associated with " + particletype +
" cluster [ns]" +
" (" + reconame +
" " + withTrigger +
") ; Time [ns] ; N" + tlabel
177 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=90,xmin=-30.,xmax=60.)
179 if particletype
in [
"Photon"] :
181 htitle =
"Photon Conversion radius [mm] " +
" (" + reconame +
" " + withTrigger +
") ; Rconv [mm] ; N_{#gamma} "
182 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=100,xmin=0.,xmax=800.)
185 htitle =
"Photon Conversion Type " +
" (" + reconame +
" " + withTrigger +
") ; ConvType ; N_{#gamma}"
186 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=6,xmin=-0.5,xmax=5.5)
188 hname=
"ConvTrkMatch1"
189 htitle =
"Photon Conversion Trak1 Match " +
" (" + reconame +
" " + withTrigger +
") ; ConTrkMatch1 ; N_{#gamma}"
190 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=4,xmin=-0.5,xmax=0.)
192 hname=
"ConvTrkMatch2"
193 htitle =
"Photon Conversion Trak1 Match " +
" (" + reconame +
" " + withTrigger +
") ; ConTrkMatch2 ; N_{#gamma}"
194 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=4,xmin=0.,xmax=4.)
199 htitle = particletype +
" energy leakage in 1st hadronic sampling " +
" (" + reconame +
" " + withTrigger +
")" +
"; Ehad1 (MeV) ; N" + tlabel
200 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-1000.,xmax=10000.)
202 if particletype
in [
"Electron"] :
204 htitle = particletype +
" matched track E over P " +
" (" + reconame +
" " + withTrigger +
")" +
"; EoverP ; N" + tlabel
205 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=5.)
208 htitle = particletype +
" core energy in EM calorimeter " +
" (" + reconame +
" " + withTrigger +
")" +
"; E (MeV) ; N" + tlabel
209 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-5000.,xmax=250000.)
212 htitle = particletype +
" fractional energy in Presampler " +
" (" + reconame +
" " + withTrigger +
")" +
"; F0 ; N" + tlabel
213 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
216 htitle = particletype +
" fractional energy in 1st sampling " +
" (" + reconame +
" " + withTrigger +
")" +
"; F1 ; N" + tlabel
217 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
220 htitle = particletype +
" fractional energy in 2nd sampling " +
" (" + reconame +
" " + withTrigger +
")" +
"; F2 ; N" + tlabel
221 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
224 htitle = particletype +
" fractional energy in 3rd sampling " +
" (" + reconame +
" " + withTrigger +
")" +
"; F3 ; N" + tlabel
225 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
228 htitle = particletype +
" uncor. energy fraction in 3x3/3x7 cells in em sampling 2 " +
" (" + reconame +
" " + withTrigger +
")" +
"; R 3x3/3x7 ; N" + tlabel
229 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=2.)
232 htitle = particletype +
" uncor. energy fraction in 3x7/7x7 cells in em sampling 2 " +
" (" + reconame +
" " + withTrigger +
")" +
"; R 3x7/7x7 ; N" + tlabel
233 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=2.)
237 if particletype
in [
"Electron"] :
238 hname=
"NOfBLayerHits"
239 htitle = particletype +
" NOfBLayerHits (" + reconame +
" " + withTrigger +
") ; N_{BlayerHits} ; N_{e}"
240 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=6,xmin=-0.5,xmax=5.5)
242 hname=
"NOfPixelHits"
243 htitle = particletype +
" NOfPixelHits (" + reconame +
" " + withTrigger +
") ; N_{PixelHits} ; N_{e}"
244 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=6,xmin=-0.5,xmax=5.5)
247 htitle = particletype +
" NOfSCTHits (" + reconame +
" " + withTrigger +
") ; N_{SCTHits} ; N_{e}"
248 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=25.5)
251 htitle = particletype +
" NOfTRTHits (" + reconame +
" " + withTrigger +
") ; N_{TRTHits} ; N_{e}"
252 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=50.5)
254 hname=
"NOfTRTHighThresholdHits"
255 htitle = particletype +
" NOfTRTHighThresholdHits (" + reconame +
" " + withTrigger +
") ; N_{TRT HighThres. Hits} ; N_{e}"
256 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=50.5)
259 htitle = particletype +
" track match #Delta #eta (1st sampling) " +
" (" + reconame +
" " + withTrigger +
") ; #Delta #eta ; N_{e} "
260 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=50,xmin=-0.05,xmax=0.05)
263 htitle = particletype +
" track match #Delta #Phi (2st sampling) " +
" (" + reconame +
" " + withTrigger +
") ; #Delta #phi ; N_{e} "
264 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=50,xmin=-0.15,xmax=0.15)
267 htitle = particletype +
" track d0 " +
" (" + reconame +
" " + withTrigger +
") ; d0 ; N_{e}"
268 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=100,xmin=-5.,xmax=5.)
273 htitle=
"Number of " + particletype +
"s per LB (" + reconame +
" " + withTrigger +
") ; LB ; N" + tlabel
274 groupe.defineHistogram(hname,title=htitle,path=
'byLB',type=
'TH1F',xbins=2000,xmin=-0.5,xmax=1999.5)
276 if particletype
in [
"Photon"] :
277 hname=
"LBEvoNPhotonsUnconv;Number_of_Unconv._photons_vs_LB"
278 htitle =
"Unconverted photon number versus LB " +
" (" + reconame +
" " + withTrigger +
") ; LB ; N_{Unconv. #gamma}"
279 groupe.defineHistogram(hname,title=htitle, path=
'byLB',xbins=2000,xmin=-0.5,xmax=1999.)
281 hname=
"LBEvoNPhotonsConv;Number_of_Conv._photons_vs_LB"
282 htitle =
"Converted photon number versus LB " +
" (" + reconame +
" " + withTrigger +
") ; LB ; N_{Conv. #gamma}"
283 groupe.defineHistogram(hname,title=htitle, path=
'byLB',xbins=2000,xmin=-0.5,xmax=1999.)
285 hname=
"is_pt_gt_4gevandconv,LB;Conversion_fraction_vs_LB"
286 htitle =
"Converted photon fraction versus LB " +
" (" + reconame +
" " + withTrigger +
") ; LB ; Conv. #gamma fraction"
287 groupe.defineHistogram(hname,title=htitle, path=
'byLB',type =
'TEfficiency', xbins=2000,xmin=-0.5,xmax=1999.)
293 Function to configure the egamma histograms for a Tag and Probe resonance decayong to ee.
294 usage example : BookHistogramsTnP(GroupElectronLHTight,"LHTight","Electron","WithTrigger")
297 tlabel =
"_{" + TnPType +
"}"
298 particletype =
'Electron'
300 hname =
"LB;Number_of_" + TnPType +
"_candidates_vs_LB"
301 htitle =
"Number of " + TnPType +
" candidates vs LB ; LB ; N" + tlabel
302 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=2000,xmin=-0.5,xmax=1999.5)
304 hname =
"MassZ;" + TnPType +
"_candidate_mass"
305 htitle = TnPType +
" candidate mass ; M_{ee} [MeV/ c^{2}]; N" + tlabel
306 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=100,xmin=MassLowerCut,xmax=MassUpperCut)
308 hname =
"MassZ_BARREL;" + TnPType +
"_candidate_mass_distribution_Lead_in_Barrel"
309 htitle = TnPType +
" candidate mass in Barrel ; M_{ee} [MeV/ c^{2}]; N" + tlabel
310 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=100,xmin=MassLowerCut,xmax=MassUpperCut)
312 hname =
"MassZ_ENDCAP;" + TnPType +
"_candidate_mass_distribution_Lead_in_EndCap"
313 htitle = TnPType +
" candidate mass in EndCap ; M_{ee} [MeV/ c^{2}]; N" + tlabel
314 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=100,xmin=MassLowerCut,xmax=MassUpperCut)
316 hname =
"MassZ_CRACK;" + TnPType +
"_candidate_mass_distribution_Lead_in_Crack"
317 htitle = TnPType +
" candidate mass in Crack Regionl ; M_{ee} [MeV/ c^{2}]; N" + tlabel
318 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=100,xmin=MassLowerCut,xmax=MassUpperCut)
322 hname =
"EtaZ;Number_of_"+ TnPType +
"_candidates_vs_eta_leading_e"
323 htitle =
"Number of "+ TnPType +
" candidates vs #eta of leading e" +
" ; #eta_{e_{lead}} ; N" + tlabel
324 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=64,xmin=-3.2,xmax=3.2)
326 hname =
"PhiZ;Number_of_"+ TnPType +
"_candidates_vs_phi_leading_e"
327 htitle =
"Number of "+ TnPType +
" candidates vs #phi of e" +
" ; #phi_{e_{lead}} ; N" + tlabel
328 groupe.defineHistogram(hname,title=htitle, path=
'TnPCandidate',xbins=64,xmin=-3.2,xmax=3.2)
332 hname =
"is_Tight,Etprobe;TnP_ID_efficiency_vs_etprobe"
333 htitle =
"Tag & Probe ID efficiency vs Et_{e} [MeV] ; Et_{e} [MeV] ; LHTight ID Eff. "
334 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=100,xmin=0.,xmax=250000.)
336 hname =
"is_Tight,Etaprobe;TnP_ID_efficiency_vs_etaprobe"
337 htitle =
"Tag & Probe ID efficiency vs #eta ; #eta _{e} ; LHTight ID Eff. "
338 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
340 hname =
"is_Tight,Phiprobe;TnP_ID_efficiency_vs_phiprobe"
341 htitle =
"Tag & Probe ID efficiency vs #eta ; #phi _{e} ; LHTight ID Eff. "
342 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
344 hname =
"is_Iso,Etprobe;TnP_Calo_Iso_efficiency_vs_etprobe"
345 htitle =
"Tag & Probe Calo. Isolation efficiency vs Et _{e} [MeV] ; Et_{e} [MeV] ; Calo Iso Eff. "
346 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=100,xmin=0.,xmax=250000.)
348 hname =
"is_Iso,Etaprobe;TnP_Calo_Iso_efficiency_vs_etaprobe"
349 htitle =
"Tag & Probe Calo. Isolation efficiency vs #eta_{e} ; #eta_{e} ; Calo. Iso Eff. "
350 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
352 hname =
"is_Iso,Phiprobe;TnP_Calo_Iso_efficiency_vs_phiprobe"
353 htitle =
"Tag & Probe Calo. Isolation efficiency vs #phi_{e} ; #phi_{e}; Calo Iso. Eff. "
354 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
356 hname =
"is_IsoandTight,Etprobe;TnP_Calo_Iso_and_LHTight_efficiency_vs_etprobe"
357 htitle =
"Tag & Probe (Isolation & LHTight ID) efficiency vs Et_{e} [MeV] ; Et_{e} [MeV] ; Calo Iso & LHTight ID Eff. "
358 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=100,xmin=0.,xmax=250000.)
360 hname =
"is_IsoandTight,Etaprobe;TnP_Calo_Iso_and_LHTight_efficiency_vs_etaprobe"
361 htitle =
"Tag & Probe (Isolation & LHTight ID) efficiency vs #eta ; #eta _{e} ; Calo. Iso & LHTight ID Eff. "
362 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
364 hname =
"is_IsoandTight,Phiprobe;TnP_Calo_Iso_and_LHTight_ID_efficiency_vs_phiprobe"
365 htitle =
"Tag & Probe (Isolation & LHTight ID) efficiency vs #phi_{e} ; #phi_{e} ; Calo Iso & LHTight ID Eff. "
366 groupe.defineHistogram(hname,title=htitle, path=
'Efficiencies',type =
'TEfficiency', xbins=64,xmin=-3.2,xmax=3.2)
371 hname =
"Etprobe;Etprobe_distribution"
372 htitle =
"Number of "+ TnPType +
" candidates vs of leading e ; Et_{probe} ; N_{e_{probe}}"
373 groupe.defineHistogram(hname,title=htitle, path=
'ElectronProbes',xbins=100,xmin=-1000.0,xmax=200000.0)
375 hname =
"Etaprobe;Etaprobe_distribution"
376 htitle =
"#eta distribution of probe candidates ; #eta_{e_{probe}} ; N_{e_{probe}}"
377 groupe.defineHistogram(hname,title=htitle,path=
'ElectronProbes', xbins=64,xmin=-3.2,xmax=3.2)
379 hname =
"Phiprobe;Phiprobe_distribution"
380 htitle =
"#phi distribution of probe candidates ; #phi_{e_{probe}} ; N_{e_{probe}}"
381 groupe.defineHistogram(hname,title=htitle,path=
'ElectronProbes', xbins=64,xmin=-3.2,xmax=3.2)
383 hname =
"Etaprobe,Phiprobe;Eta_Phi_map_of_electron_probes_for_T_and_P_"+ TnPType
384 htitle=
"Concerted photons #eta,#phi map (candidates with Pt>4GeV)" +
" ; #eta_{e} ; #phi_{e}"
385 groupe.defineHistogram(hname,title=htitle,path=
'ElectronProbes',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2)
387 hname =
"Timeprobe;Time_of_electron_probe"
388 htitle =
"Time distribution of probe candidates ; Time_{e_{probe}} [ns] ; N_{e_{probe}}"
389 groupe.defineHistogram(hname,title=htitle,path=
'ElectronProbes',xbins=90,xmin=-30.,xmax=60.)
391 hname =
"TopoEtCone40probe;TopoEtCone40_of_electron_probe"
392 htitle =
"Electron probe Topocluster Isolation Energy ; TopoEtCone40 [MeV] ; N_{e_{probe}} "
393 groupe.defineHistogram(hname,title=htitle, path=
'ElectronProbes',xbins=100,xmin=-10000.,xmax=40000.)
395 hname =
"PtCone20probe;PtCone20_of_electron_probe"
396 htitle =
"Electron probe Track Isolation Pt [MeV]; PtCone20 [MeV] ; N_{e_{probe}} "
397 groupe.defineHistogram(hname,title=htitle, path=
'ElectronProbes',xbins=64,xmin=-10000.,xmax=40000.)
401 hname =
"NOfBLayerHitsProbe"
402 htitle =
"Electron probe NOfBLayerHits ; N_{BlayerHits} ; N_{e_{probe}}"
403 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=6,xmin=-0.5,xmax=5.5)
405 hname =
"NOfPixelHitsProbe"
406 htitle =
"Electron probe NOfPixelHits ; N_{PixelHits} ; N_{e_{probe}}"
407 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=6,xmin=-0.5,xmax=5.5)
409 hname =
"NOfSCTHitsProbe"
410 htitle =
"Electron probe NOfSCTHits ; N_{SCTHits} ; N_{e_{probe}}"
411 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=25.5)
413 hname =
"NOfTRTHitsProbe"
414 htitle =
"Electron probe NOfTRTHits ; N_{TRTHits} ; N_{e_{probe}}"
415 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=50.5)
417 hname =
"NOfTRTHighThresholdHitsProbe"
418 htitle =
"Electron probe NOfTRTHighThresholdHits ; N_{TRT HighThres. Hits} ; N_{e_{probe}}"
419 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=26,xmin=-0.5,xmax=50.5)
421 hname =
"DeltaEta1Probe"
422 htitle =
"Electron probe track match #Delta #eta (1st sampling) ; #Delta #eta ; N_{e_{probe}} "
423 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=50,xmin=-0.05,xmax=0.05)
425 hname =
"DeltaPhi2Probe"
426 htitle =
"Electron probe track match #Delta #Phi (2st sampling) ; #Delta #phi ; N_{e_{probe}} "
427 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=50,xmin=-0.15,xmax=0.15)
429 hname =
"Trackd0Probe"
430 htitle =
"Electron probe track d0 ; d0 ; N_{e_{probe}}"
431 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Tracks',xbins=100,xmin=-5.,xmax=5.)
436 htitle =
"Electron probe energy leakage in 1st hadronic sampling ; Ehad1 (MeV) ; N_{e_{probe}} "
437 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-1000.,xmax=10000.)
439 hname =
"EoverPProbe"
440 htitle =
"Electron probe matched track E over P ; EoverP ; N_{e_{probe}} "
441 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=5.)
443 hname =
"CoreEMProbe"
444 htitle =
"Electron probe core energy in EM calorimeter ; E (MeV) ; N_{e_{probe}} "
445 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-5000.,xmax=250000.)
448 htitle =
"Electron probe fractional energy in Presampler ; F0 ; N_{e_{probe}} "
449 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
452 htitle =
"Electron probe fractional energy in 1st sampling ; F1 ; N_{e_{probe}} "
453 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
456 htitle =
"Electron probe fractional energy in 2nd sampling ; F2 ; N_{e_{probe}} "
457 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
460 htitle =
"Electron probe fractional energy in 3rd sampling ; F3 ; N_{e_{probe}} "
461 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=-0.2,xmax=1.)
463 hname =
"Re233e237Probe"
464 htitle =
"Electron probe uncor. energy fraction in 3x3/3x7 cells in em sampling ; R 3x3/3x7 ; N_{e_{probe}} "
465 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=2.)
467 hname =
"Re237e277Probe"
468 htitle =
"Electron probe uncor. energy fraction in 3x7/7x7 cells in em sampling 2 ; R 3x7/7x7 ; N_{e_{probe}} "
469 BookHistogramsPerRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'ID',xbins=50,xmin=0.,xmax=2.)
471 hname =
"LBEvoNProbe"
472 htitle=
"Number of " + TnPType +
" electron probe per LB ; LB ; N_{e_{probe}} "
473 groupe.defineHistogram(hname,title=htitle,path=
'byLB',type=
'TH1F',xbins=2000,xmin=-0.5,xmax=1999.5)
480 Function to configure the forward electrons histograms
481 usage example : BookHistograms(GroupElectronLHTight,"LHTight","Electron","WithTrigger")
484 prefix =
"fwdElectron"
485 tlabel =
"_{fwd. electron}"
488 htitle =
"Number of " + prefix +
"s (" + reconame +
" " + withTrigger+
") ; N" + tlabel +
" ; N_{event}"
489 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
492 htitle =
"Number of " + prefix +
"s in ENDCAP (" + reconame +
" " + withTrigger+
") ; N" + tlabel +
" ; N_{event}"
493 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
496 htitle =
"Number of " + prefix +
"s in FORWARD (" + reconame +
" " + withTrigger+
") ; N" + tlabel +
" ; N_{event}"
497 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=20,xmin=-0.5,xmax=19.5)
500 htitle= particletype +
" transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Et"+ tlabel +
" ; N" + tlabel
501 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=2.,xmax=102000.0)
504 htitle= particletype +
" #eta" +
" (" + reconame +
" " + withTrigger +
")"+
" ; #eta" + tlabel +
" ; N" + tlabel
505 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2)
508 htitle= particletype +
" #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #phi" + tlabel +
" ; N" + tlabel
509 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=64,xmin=-3.2,xmax=3.2)
511 hname=
"Eta,Phi;Eta_Phi_distribution_Pt_gt_2.5GeV"
512 htitle= particletype +
" #eta,#phi map (candidates with Pt>2.5GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi"
513 groupe.defineHistogram(hname,title=htitle,path=
'Expert',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask =
'is_pt_gt_2_5gev')
515 hname=
"Eta,Phi;Eta_Phi_distributionf_Pt_gt_10GeV"
516 htitle= particletype +
" #eta,#phi map (candidates with Pt>10GeV)" +
" (" + reconame +
" " + withTrigger +
") ; #eta ; #phi "
517 groupe.defineHistogram(hname,title=htitle,path=
'',type=
'TH2F',xbins=64,xmin=-3.2,xmax=3.2,ybins=64,ymin=-3.2,ymax=3.2,cutmask=
'is_pt_gt_10gev')
521 hname=
"EnergyDensity"
522 htitle=
"Fwd electron 1st Moment Energy Density" +
" (" + reconame +
" " + withTrigger +
")" +
" ; firstENGdensity ; N" + tlabel
523 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=200,xmin=0.,xmax=2.0)
526 htitle=
"Fwd electron fraction of most energetic cell distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; lateral moment ; N" + tlabel
527 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=50,xmin=0.,xmax=1.0)
530 htitle=
"Fwd electron lateral moment distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; lateral moment ; N" + tlabel
531 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=0.,xmax=1.0)
533 hname=
"Longitudinal"
534 htitle=
"Fwd electron longitudinal moment distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; longitudinal moment ; N" + tlabel
535 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=100,xmin=0.,xmax=1.0)
537 hname=
"SecondLambda"
538 htitle=
"Fwd electron lambda second moment distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Second#Lambda ; N" + tlabel
539 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=500,xmin=0.,xmax=10000.0)
542 htitle=
"Fwd electron SecondR lateral moment distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; SecondR ; N" + tlabel
543 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=500,xmin=0.,xmax=20000.0)
545 hname=
"CenterLambda"
546 htitle=
"Fwd Electron shower center from calo front face distribution" +
" (" + reconame +
" " + withTrigger +
")" +
" ; Center lambda ; N" + tlabel
547 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=500,xmin=0.,xmax=2000.0)
550 htitle =
"Time associated with " + particletype +
" cluster [ns]" +
" (" + reconame +
" " + withTrigger +
") ; Time [ns] ; N" + tlabel
551 groupe.defineHistogram(hname,title=htitle, path=
'',xbins=90,xmin=-30.,xmax=60.)
557 htitle= particletype +
" transverse energy [MeV]" +
" (" + reconame +
" " + withTrigger +
") ; Et" + tlabel +
" ; N" + tlabel
558 BookHistogramsPerForwardRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=100,xmin=2.,xmax=102000.0,thecut=
'is_pt_gt_2_5gev')
561 htitle= particletype +
" #eta" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #eta" + tlabel +
" ; N" + tlabel
562 BookHistogramsPerForwardRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=64,xmin=-3.2,xmax=3.2)
565 htitle= particletype +
" #phi" +
" (" + reconame +
" " + withTrigger +
")" +
" ; #eta" + tlabel +
" ; N" + tlabel
566 BookHistogramsPerForwardRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=64,xmin=-3.2,xmax=3.2)
569 htitle =
"Time associated with " + particletype +
" cluster [ns]" +
" (" + reconame +
" " + withTrigger +
") ; Time [ns] ; N" + tlabel
570 BookHistogramsPerForwardRegions(thegroupe = groupe, theparttype = particletype, thename = hname, title=htitle, path=
'Expert',xbins=90,xmin=-30.,xmax=60.)
575 '''Function to configures some algorithms in the monitoring system.'''
579 from AthenaMonitoring
import AthMonitorCfgHelper
580 helper = AthMonitorCfgHelper(inputFlags,
'ElectronAthMonitorCfg')
584 SpareElectronMonitoringGroups = [ ]
586 MonitorElectronAlgorithm=CompFactory.MonitorElectronAlgorithm
587 elLHTightMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elLHTightMonAlg')
588 elLHTightTrigMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elLHTightTrigMonAlg')
590 elLHLooseMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elLHLooseMonAlg')
591 elLHLooseTrigMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elLHLooseTrigMonAlg')
593 if (
"CBTight" in SpareElectronMonitoringGroups) :
594 elCBTightMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elCBTightMonAlg')
595 if (
"CBTightTrig" in SpareElectronMonitoringGroups) :
596 elCBTightTrigMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elCBTightTrigMonAlg')
597 if (
"CBLoose" in SpareElectronMonitoringGroups) :
598 elCBLooseMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elCBLooseMonAlg')
599 if (
"CBLooseTrig" in SpareElectronMonitoringGroups) :
600 elCBLooseTrigMonAlg = helper.addAlgorithm(MonitorElectronAlgorithm,
'elCBLooseTrigMonAlg')
605 elLHTightMonAlg.ParticleContainerName =
"Electrons"
606 elLHTightMonAlg.RecoName =
"LHTight"
607 elLHTightMonAlg.ParticlePrefix =
"electron"
608 elLHTightMonAlg.WithTrigger =
""
610 elLHTightTrigMonAlg.ParticleContainerName =
"Electrons"
611 elLHTightTrigMonAlg.RecoName =
"LHTight"
612 elLHTightTrigMonAlg.ParticlePrefix =
"electron"
613 elLHTightTrigMonAlg.WithTrigger =
"WithTrigger"
615 elLHTightTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
617 elLHLooseMonAlg.ParticleContainerName =
"Electrons"
618 elLHLooseMonAlg.RecoName =
"LHLoose"
619 elLHLooseMonAlg.ParticlePrefix =
"electron"
621 elLHLooseTrigMonAlg.ParticleContainerName =
"Electrons"
622 elLHLooseTrigMonAlg.RecoName =
"LHLoose"
623 elLHLooseTrigMonAlg.ParticlePrefix =
"electron"
624 elLHLooseTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
626 if (
"CBTight" in SpareElectronMonitoringGroups) :
627 elCBTightMonAlg.ParticleContainerName =
"Electrons"
628 elCBTightMonAlg.RecoName =
"Tight"
629 elCBTightMonAlg.ParticlePrefix =
"electron"
631 if (
"CBTightTrig" in SpareElectronMonitoringGroups) :
632 elCBTightTrigMonAlg.ParticleContainerName =
"Electrons"
633 elCBTightTrigMonAlg.RecoName =
"Tight"
634 elCBTightTrigMonAlg.ParticlePrefix =
"electron"
635 elCBTightTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
637 if (
"CBLoose" in SpareElectronMonitoringGroups) :
638 elCBLooseMonAlg.ParticleContainerName =
"Electrons"
639 elCBLooseMonAlg.RecoName =
"Loose"
640 elCBLooseMonAlg.ParticlePrefix =
"electron"
642 if (
"CBLooseTrig" in SpareElectronMonitoringGroups) :
643 elCBLooseTrigMonAlg.ParticleContainerName =
"Electrons"
644 elCBLooseTrigMonAlg.RecoName =
"Loose"
645 elCBLooseTrigMonAlg.ParticlePrefix =
"electron"
646 elCBLooseTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
650 GroupElectronLHTight = helper.addGroup(elLHTightMonAlg,
'MonitorElectron',
'egamma/LHTightElectrons/')
651 GroupElectronLHTightTriggered = helper.addGroup(elLHTightTrigMonAlg,
'MonitorElectron',
'egamma/LHTightElectronsWithTrigger/')
652 GroupElectronLHLoose = helper.addGroup(elLHLooseMonAlg,
'MonitorElectron',
'egamma/LHLooseElectrons/')
653 GroupElectronLHLooseTriggered = helper.addGroup(elLHLooseTrigMonAlg,
'MonitorElectron',
'egamma/LHLooseElectronsWithTrigger/')
654 if (
"CBTight" in SpareElectronMonitoringGroups) :
655 GroupElectronCBTight = helper.addGroup(elCBTightMonAlg,
'MonitorElectron',
'egamma/CBTightElectrons/')
656 if (
"CBTightTrig" in SpareElectronMonitoringGroups) :
657 GroupElectronCBTightTriggered = helper.addGroup(elCBTightTrigMonAlg,
'MonitorElectron',
'egamma/CBTightTriggeredElectrons/')
658 if (
"CBLoose" in SpareElectronMonitoringGroups) :
659 GroupElectronCBLoose = helper.addGroup(elCBLooseMonAlg,
'MonitorElectron',
'egamma/CBLooseElectrons/')
660 if (
"CBLooseTrig" in SpareElectronMonitoringGroups) :
661 GroupElectronCBLooseTriggered = helper.addGroup(elCBLooseTrigMonAlg,
'MonitorElectron',
'egamma/CBLooseTriggeredElectrons/')
668 BookHistograms(GroupElectronLHTightTriggered,
"LHTightTrig",
"Electron",
"WithTrigger")
671 BookHistograms(GroupElectronLHLooseTriggered,
"LHLooseTrig",
"Electron",
"WithTrigger")
673 if (
"CBTight" in SpareElectronMonitoringGroups) :
675 if (
"CBTightTrig" in SpareElectronMonitoringGroups) :
676 BookHistograms(GroupElectronCBTightTriggered,
"CBTightTrig",
"Electron",
"WithTrigger")
678 if (
"CBLoose" in SpareElectronMonitoringGroups) :
680 if (
"CBLooseTrig" in SpareElectronMonitoringGroups) :
681 BookHistograms(GroupElectronCBLooseTriggered,
"CBLooseTrig",
"Electron",
"WithTrigger")
684 return helper.result()
688 '''Function to configures some algorithms in the monitoring system.'''
692 from AthenaMonitoring
import AthMonitorCfgHelper
693 helper = AthMonitorCfgHelper(inputFlags,
'PhotonAthMonitorCfg')
697 MonitorPhotonAlgorithm=CompFactory.MonitorPhotonAlgorithm
698 phCBTightMonAlg = helper.addAlgorithm(MonitorPhotonAlgorithm,
'phCBTightMonAlg')
699 phCBTightTrigMonAlg = helper.addAlgorithm(MonitorPhotonAlgorithm,
'phCBTightTrigMonAlg')
701 phCBLooseMonAlg = helper.addAlgorithm(MonitorPhotonAlgorithm,
'phCBLooseMonAlg')
702 phCBLooseTrigMonAlg = helper.addAlgorithm(MonitorPhotonAlgorithm,
'phCBLooseTrigMonAlg')
707 phCBTightMonAlg.ParticleContainerName =
"Photons"
708 phCBTightMonAlg.RecoName =
"Tight"
709 phCBTightMonAlg.ParticlePrefix =
"photon"
711 phCBTightTrigMonAlg.ParticleContainerName =
"Photons"
712 phCBTightTrigMonAlg.RecoName =
"Tight"
713 phCBTightTrigMonAlg.ParticlePrefix =
"photon"
714 phCBTightTrigMonAlg.WithTrigger =
"WithTrigger"
716 phCBTightTrigMonAlg.TriggerChain = egammaConf.primary_double_pho[0]
718 phCBLooseMonAlg.ParticleContainerName =
"Photons"
719 phCBLooseMonAlg.RecoName =
"Loose"
720 phCBLooseMonAlg.ParticlePrefix =
"photon"
722 phCBLooseTrigMonAlg.ParticleContainerName =
"Photons"
723 phCBLooseTrigMonAlg.RecoName =
"Loose"
724 phCBLooseTrigMonAlg.ParticlePrefix =
"photon"
725 phCBLooseTrigMonAlg.WithTrigger =
"WithTrigger"
727 phCBLooseTrigMonAlg.TriggerChain = egammaConf.primary_double_pho[0]
732 GroupPhotonCBTight = helper.addGroup(phCBTightMonAlg,
'MonitorPhoton',
'egamma/CBTightPhotons/')
733 GroupPhotonCBTightTriggered = helper.addGroup(phCBTightTrigMonAlg,
'MonitorPhoton',
'egamma/CBTightPhotonsWithTrigger/')
734 GroupPhotonCBLoose = helper.addGroup(phCBLooseMonAlg,
'MonitorPhoton',
'egamma/CBLoosePhotons/')
735 GroupPhotonCBLooseTriggered = helper.addGroup(phCBLooseTrigMonAlg,
'MonitorPhoton',
'egamma/CBLoosePhotonsWithTrigger/')
742 BookHistograms(GroupPhotonCBTightTriggered,
"CBTightTrig",
"Photon")
745 BookHistograms(GroupPhotonCBLooseTriggered,
"CBLooseTrig",
"Photon")
748 return helper.result()
752 '''Function to configures some algorithms in the monitoring system for Tag and Probe candidates (Z, J/Psi).'''
756 from AthenaMonitoring
import AthMonitorCfgHelper
757 helper = AthMonitorCfgHelper(inputFlags,
'TnPPAthMonitorCfg')
762 MonitorTnPAlgorithm = CompFactory.MonitorTnPAlgorithm
763 ZeeMonAlg = helper.addAlgorithm(MonitorTnPAlgorithm,
'TnPZeeMonAlg')
764 JPsiMonAlg = helper.addAlgorithm(MonitorTnPAlgorithm,
'TnPJpsiMonAlg')
770 ZeeMonAlg.ParticleContainerName =
"Electrons"
771 ZeeMonAlg.RecoName =
"LHLoose"
772 ZeeMonAlg.ParticlePrefix =
"electron"
773 ZeeMonAlg.MassPeak = 91188.
774 ZeeMonAlg.ElectronEtCut = 15000.
775 ZeeMonAlg.MassLowerCut = 70000.
776 ZeeMonAlg.MassUpperCut = 110000.
777 ZeeMonAlg.TnPType =
"Z"
781 for el
in egammaConf.monitoring_Zee:
783 chain = chain +
", " + el
788 JPsiMonAlg.ParticleContainerName =
"Electrons"
789 JPsiMonAlg.RecoName =
"LHLoose"
790 JPsiMonAlg.ParticlePrefix =
"electron"
791 JPsiMonAlg.MassPeak = 3097.
792 JPsiMonAlg.ElectronEtCut = 3000.
793 JPsiMonAlg.MassLowerCut = 2000.
794 JPsiMonAlg.MassUpperCut = 5000.
795 JPsiMonAlg.TnPType =
"JPsi"
799 for el
in egammaConf.monitoring_Jpsiee:
801 chain = chain +
", " + el
804 JPsiMonAlg.TriggerChain = chain
808 GroupZee = helper.addGroup(ZeeMonAlg,
'MonitorTnP',
'egamma/TnPZ/')
809 GroupJPsi = helper.addGroup(JPsiMonAlg,
'MonitorTnP',
'egamma/TnPJPsi/')
815 BookTnPHistograms(GroupZee, ZeeMonAlg.RecoName, ZeeMonAlg.TnPType, ZeeMonAlg.MassLowerCut, ZeeMonAlg.MassUpperCut)
817 BookTnPHistograms(GroupJPsi, JPsiMonAlg.RecoName, JPsiMonAlg.TnPType, JPsiMonAlg.MassLowerCut, JPsiMonAlg.MassUpperCut)
820 return helper.result()
824 '''Function to configures some algorithms in the monitoring system.'''
828 from AthenaMonitoring
import AthMonitorCfgHelper
829 helper = AthMonitorCfgHelper(inputFlags,
'FwdElectronAthMonitorCfg')
833 SpareFwdElectronMonitoringGroups = [
"CBFwdTight" ]
835 MonitorFwdElectronAlgorithm=CompFactory.MonitorFwdElectronAlgorithm
838 if (
"CBFwdTight" in SpareFwdElectronMonitoringGroups) :
839 fwdelCBTightMonAlg = helper.addAlgorithm(MonitorFwdElectronAlgorithm,
'fwdelCBTightMonAlg')
840 if (
"CBfwdTightTrig" in SpareFwdElectronMonitoringGroups) :
841 fwdelCBTightTrigMonAlg = helper.addAlgorithm(MonitorFwdElectronAlgorithm,
'fwdelCBTightTrigMonAlg')
842 if (
"CBLoose" in SpareFwdElectronMonitoringGroups) :
843 fwdelCBLooseMonAlg = helper.addAlgorithm(MonitorFwdElectronAlgorithm,
'fwdelCBLooseMonAlg')
844 if (
"CBLooseTrig" in SpareFwdElectronMonitoringGroups) :
845 fwdelCBLooseTrigMonAlg = helper.addAlgorithm(MonitorFwdElectronAlgorithm,
'fwdelCBLooseTrigMonAlg')
871 if (
"CBFwdTight" in SpareFwdElectronMonitoringGroups) :
872 fwdelCBTightMonAlg.ParticleContainerName =
"ForwardElectrons"
873 fwdelCBTightMonAlg.RecoName =
"Tight"
874 fwdelCBTightMonAlg.ParticlePrefix =
"FwdElectron"
876 if (
"CBFwdTightTrig" in SpareFwdElectronMonitoringGroups) :
877 fwdelCBTightTrigMonAlg.ParticleContainerName =
"ForwardElectrons"
878 fwdelCBTightTrigMonAlg.RecoName =
"Tight"
879 fwdelCBTightTrigMonAlg.ParticlePrefix =
"FwdElectron"
881 fwdelCBTightTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
883 if (
"CBFwdLoose" in SpareFwdElectronMonitoringGroups) :
884 fwdelCBLooseMonAlg.ParticleContainerName =
"ForwardElectrons"
885 fwdelCBLooseMonAlg.RecoName =
"Loose"
886 fwdelCBLooseMonAlg.ParticlePrefix =
"FwdElectron"
888 if (
"CBFwdLooseTrig" in SpareFwdElectronMonitoringGroups) :
889 fwdelCBLooseTrigMonAlg.ParticleContainerName =
"ForwardElectrons"
890 fwdelCBLooseTrigMonAlg.RecoName =
"Loose"
891 fwdelCBLooseTrigMonAlg.ParticlePrefix =
"FwdElectron"
893 fwdelCBLooseTrigMonAlg.TriggerChain = egammaConf.primary_single_ele[0]
901 if (
"CBFwdTight" in SpareFwdElectronMonitoringGroups) :
902 GroupFwdElectronCBTight = helper.addGroup(fwdelCBTightMonAlg,
'MonitorFwdElectron',
'egamma/CBTightFwdElectrons/')
903 if (
"CBFwdTightTrig" in SpareFwdElectronMonitoringGroups) :
904 GroupFwdElectronCBTightTriggered = helper.addGroup(fwdelCBTightTrigMonAlg,
'MonitorFwdElectron',
'egamma/CBTightTriggeredFwdElectrons/')
905 if (
"CBFwdLoose" in SpareFwdElectronMonitoringGroups) :
906 GroupFwdElectronCBLoose = helper.addGroup(fwdelCBLooseMonAlg,
'MonitorFwdElectron',
'egamma/CBLooseFwdElectrons/')
907 if (
"CBFwdLooseTrig" in SpareFwdElectronMonitoringGroups) :
908 GroupFwdElectronCBLooseTriggered = helper.addGroup(fwdelCBLooseTrigMonAlg,
'MonitorFwdElectron',
'egamma/CBLooseTriggeredFwdElectrons/')
920 if (
"CBFwdTight" in SpareFwdElectronMonitoringGroups) :
922 if (
"CBFwdTightTrig" in SpareFwdElectronMonitoringGroups) :
925 if (
"CBFwdLoose" in SpareFwdElectronMonitoringGroups) :
927 if (
"CBFwdLooseTrig" in SpareFwdElectronMonitoringGroups) :
931 return helper.result()
934 if __name__==
'__main__':
936 from AthenaCommon.Logging
import log
941 from AthenaConfiguration.AllConfigFlags
import initConfigFlags
948 flags.Input.Files = [
'/eos/atlas/atlascerngroupdisk/trig-egam/ForKyle/AOD/valid/user.eegidiop/user.eegidiop.35311896.EXT0._000012.AOD.root']
949 flags.Input.isMC =
True
953 flags.Output.HISTFileName =
'data16_13TeV.00311321.physics_Main.merge.HIST.f1156_h347._0002.1.root'
958 from AthenaConfiguration.MainServicesConfig
import MainServicesCfg
959 from AthenaPoolCnvSvc.PoolReadConfig
import PoolReadCfg
964 cfg.merge(MonitorElectronAcc)
967 cfg.merge(MonitorPhotonAcc)
970 cfg.merge(MonitorTnPAcc)
973 cfg.merge(MonitorFwdElectronAcc)
984 cfg.merge(MonitorFwdElectronAcc)
986 cfg.printConfig(withDetails=
False)