CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
hgcalPlots.py
Go to the documentation of this file.
1 from __future__ import print_function
2 import os
3 import sys
4 import copy
5 import collections
6 
7 import six
8 import ROOT
9 from ROOT import TFile
10 from ROOT import gDirectory
11 ROOT.gROOT.SetBatch(True)
12 ROOT.PyConfig.IgnoreCommandLineOptions = True
13 
14 from Validation.RecoTrack.plotting.plotting import Plot, PlotGroup, PlotFolder, Plotter, PlotOnSideGroup
15 from Validation.RecoTrack.plotting.html import PlotPurpose
16 import Validation.RecoTrack.plotting.plotting as plotting
17 import Validation.RecoTrack.plotting.validation as validation
18 import Validation.RecoTrack.plotting.html as html
19 
20 #To be able to spot any issues both in -z and +z a layer id was introduced
21 #that spans from 0 to 103 for hgcal_v9 geometry. The mapping for hgcal_v9 is:
22 #-z: 0->51
23 #+z: 52->103
24 #while for V10 is:
25 #-z: 0->49
26 #+z: 50->99
27 '''
28 layerscheme = { 'lastLayerEEzm': 0, 'lastLayerFHzm': 0, 'maxlayerzm': 0, 'lastLayerEEzp': 0, 'lastLayerFHzp': 0, 'maxlayerzp': 0 }
29 
30 #Let's take the relevant values of layerscheme from the dqm file.
31 theDQMfile = "DQM_V0001_R000000001__Global__CMSSW_X_Y_Z__RECO.root"
32 if not os.path.isfile(theDQMfile):
33  print("Error: file", theDQMfile, "not found, exit")
34  sys.exit(0)
35 
36 
37 #Take general info from the first file is sufficient.
38 thefile = TFile( theDQMfile )
39 GeneralInfoDirectory = 'DQMData/Run 1/HGCAL/Run summary/HGCalValidator/GeneralInfo'
40 
41 if not gDirectory.GetDirectory( GeneralInfoDirectory ):
42  print("Error: GeneralInfo directory not found in DQM file, exit")
43  sys.exit(0)
44 
45 keys = gDirectory.GetDirectory( GeneralInfoDirectory ).GetListOfKeys()
46 key = keys[0]
47 layvalue = 0
48 while key:
49  obj = key.ReadObj()
50  for laykey in layerscheme.keys():
51  if laykey in obj.GetName():
52  layvalue = obj.GetName()[len("<"+laykey+">i="):-len("</"+laykey+">")]
53  layerscheme[laykey] = layvalue
54  #print(layvalue)
55  key = keys.After(key)
56 
57 thefile.Close()
58 
59 print(layerscheme)
60 #TODO: Anticipating the fine/coarse layer information in CMSSW we overwrite values from DQM file
61 #For now values returned for
62 # 'lastLayerFHzp': '104', 'lastLayerFHzm': '52'
63 #are not the one expected. Will come back to this when there will be info in CMSSW to put in DQM file.
64 #For V9:
65 #layerscheme = { 'lastLayerEEzm': 28, 'lastLayerFHzm': 40, 'maxlayerzm': 52, 'lastLayerEEzp': 80, 'lastLayerFHzp': 92, 'maxlayerzp': 104 }
66 #For V10:
67 '''
68 layerscheme = { 'lastLayerEEzm': 28, 'lastLayerFHzm': 40, 'maxlayerzm': 50, 'lastLayerEEzp': 78, 'lastLayerFHzp': 90, 'maxlayerzp': 100 }
69 #print(layerscheme)
70 
71 lastLayerEEzm = layerscheme['lastLayerEEzm'] # last layer of EE -z
72 lastLayerFHzm = layerscheme['lastLayerFHzm'] # last layer of FH -z
73 maxlayerzm = layerscheme['maxlayerzm'] # last layer of BH -z
74 lastLayerEEzp = layerscheme['lastLayerEEzp'] # last layer of EE +z
75 lastLayerFHzp = layerscheme['lastLayerFHzp'] # last layer of FH +z
76 maxlayerzp = layerscheme['maxlayerzp'] # last layer of BH +z
77 
78 hitlayerscheme = { 'EE_min': 1,'EE_max': 28, 'HESilicon_min': 1, 'HESilicon_max': 22, 'HEScintillator_min': 9 , 'HEScintillator_max': 22 }
79 #print(hitlayerscheme)
80 
81 EE_min = hitlayerscheme['EE_min']
82 EE_max = hitlayerscheme['EE_max']
83 HESilicon_min = hitlayerscheme['HESilicon_min']
84 HESilicon_max = hitlayerscheme['HESilicon_max']
85 HEScintillator_min = hitlayerscheme['HEScintillator_min']
86 HEScintillator_max = hitlayerscheme['HEScintillator_max']
87 
88 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
89 _legend_common = {"legendDx": -0.3,
90  "legendDy": -0.05,
91  "legendDw": 0.1}
92 
93 _SelectedCaloParticles = PlotGroup("SelectedCaloParticles", [
94  Plot("num_caloparticle_eta", xtitle="", **_common),
95  Plot("caloparticle_energy", xtitle="", **_common),
96  Plot("caloparticle_pt", xtitle="", **_common),
97  Plot("caloparticle_phi", xtitle="", **_common),
98  Plot("Eta vs Zorigin", xtitle="", **_common),
99  ])
100 
101 #Need to adjust the statbox to see better the plot
102 _common = {"stat": True, "drawStyle": "hist", "statx": 0.38, "staty": 0.68 }
103 _num_reco_cluster_eta = PlotGroup("num_reco_cluster_eta", [
104  Plot("num_reco_cluster_eta", xtitle="", **_common),
105 ],ncols=1)
106 #Back to normal
107 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
108 
109 _mixedhitsclusters = PlotGroup("mixedhitsclusters", [
110  Plot("mixedhitscluster_zminus", xtitle="", **_common),
111  Plot("mixedhitscluster_zplus", xtitle="", **_common),
112 ],ncols=2)
113 
114 #Just to prevent the stabox covering the plot
115 _common = {"stat": True, "drawStyle": "hist", "statx": 0.45, "staty": 0.65 }
116 
117 _energyclustered = PlotGroup("energyclustered", [
118  Plot("energyclustered_zminus", xtitle="", **_common),
119  Plot("energyclustered_zplus", xtitle="", **_common),
120 ],ncols=2)
121 
122 #Coming back to the usual box definition
123 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
124 
125 _longdepthbarycentre = PlotGroup("longdepthbarycentre", [
126  Plot("longdepthbarycentre_zminus", xtitle="", **_common),
127  Plot("longdepthbarycentre_zplus", xtitle="", **_common),
128 ],ncols=2)
129 
130 _common_layerperthickness = {}
131 _common_layerperthickness.update(_common)
132 _common_layerperthickness['xmin'] = 0.
133 _common_layerperthickness['xmax'] = 100
134 
135 _totclusternum_thick = PlotGroup("totclusternum_thick", [
136  Plot("totclusternum_thick_120", xtitle="", **_common_layerperthickness),
137  Plot("totclusternum_thick_200", xtitle="", **_common_layerperthickness),
138  Plot("totclusternum_thick_300", xtitle="", **_common_layerperthickness),
139  Plot("totclusternum_thick_-1", xtitle="", **_common_layerperthickness),
140  Plot("mixedhitscluster", xtitle="", **_common_layerperthickness),
141 ])
142 
143 #We will plot the density in logy scale.
144 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ylog": True}
145 
146 _cellsenedens_thick = PlotGroup("cellsenedens_thick", [
147  Plot("cellsenedens_thick_120", xtitle="", **_common),
148  Plot("cellsenedens_thick_200", xtitle="", **_common),
149  Plot("cellsenedens_thick_300", xtitle="", **_common),
150  Plot("cellsenedens_thick_-1", xtitle="", **_common),
151 ])
152 
153 #Coming back to the usual box definition
154 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
155 
156 
157 #--------------------------------------------------------------------------------------------
158 # z-
159 #--------------------------------------------------------------------------------------------
160 _totclusternum_layer_EE_zminus = PlotGroup("totclusternum_layer_EE_zminus", [
161  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
162 ], ncols=4)
163 
164 _totclusternum_layer_FH_zminus = PlotGroup("totclusternum_layer_FH_zminus", [
165  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
166 ], ncols=4)
167 
168 _totclusternum_layer_BH_zminus = PlotGroup("totclusternum_layer_BH_zminus", [
169  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
170 ], ncols=4)
171 
172 _energyclustered_perlayer_EE_zminus = PlotGroup("energyclustered_perlayer_EE_zminus", [
173  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
174 ], ncols=4)
175 
176 _energyclustered_perlayer_FH_zminus = PlotGroup("energyclustered_perlayer_FH_zminus", [
177  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
178 ], ncols=4)
179 
180 _energyclustered_perlayer_BH_zminus = PlotGroup("energyclustered_perlayer_BH_zminus", [
181  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
182 ], ncols=4)
183 
184 #----------------------------------------------------------------------------------------------------------------
185 #120 um
186 _common_cells = {}
187 _common_cells.update(_common)
188 _common_cells["xmin"] = 0
189 _common_cells["xmax"] = 50
190 _common_cells["ymin"] = 0.1
191 _common_cells["ymax"] = 10000
192 _common_cells["ylog"] = True
193 _cellsnum_perthick_perlayer_120_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_120_EE_zminus", [
194  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
195 ], ncols=4)
196 
197 _cellsnum_perthick_perlayer_120_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_120_FH_zminus", [
198  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
199 ], ncols=4)
200 
201 _cellsnum_perthick_perlayer_120_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_120_BH_zminus", [
202  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
203 ], ncols=4)
204 
205 #200 um
206 _cellsnum_perthick_perlayer_200_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_200_EE_zminus", [
207  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
208 ], ncols=4)
209 
210 _cellsnum_perthick_perlayer_200_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_200_FH_zminus", [
211  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
212 ], ncols=4)
213 
214 _cellsnum_perthick_perlayer_200_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_200_BH_zminus", [
215  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
216 ], ncols=4)
217 
218 #300 um
219 _cellsnum_perthick_perlayer_300_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_300_EE_zminus", [
220  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
221 ], ncols=4)
222 
223 _cellsnum_perthick_perlayer_300_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_300_FH_zminus", [
224  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
225 ], ncols=4)
226 
227 _cellsnum_perthick_perlayer_300_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_300_BH_zminus", [
228  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
229 ], ncols=4)
230 
231 #scint um
232 _cellsnum_perthick_perlayer_scint_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_EE_zminus", [
233  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
234 ], ncols=4)
235 
236 _cellsnum_perthick_perlayer_scint_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_FH_zminus", [
237  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
238 ], ncols=4)
239 
240 _cellsnum_perthick_perlayer_scint_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_BH_zminus", [
241  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
242 ], ncols=4)
243 
244 #----------------------------------------------------------------------------------------------------------------
245 #120 um
246 _common_distance = {}
247 _common_distance.update(_common)
248 _common_distance.update(_legend_common)
249 _common_distance["xmax"] = 150
250 _common_distance["stat"] = False
251 _common_distance["ymin"] = 1e-3
252 _common_distance["ymax"] = 10000
253 _common_distance["ylog"] = True
254 
255 _distancetomaxcell_perthickperlayer_120_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_EE_zminus", [
256  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
257 ], ncols=4)
258 
259 _distancetomaxcell_perthickperlayer_120_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_FH_zminus", [
260  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
261 ], ncols=4)
262 
263 _distancetomaxcell_perthickperlayer_120_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_BH_zminus", [
264  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
265 ], ncols=4)
266 
267 #200 um
268 _distancetomaxcell_perthickperlayer_200_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_EE_zminus", [
269  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
270 ], ncols=4)
271 
272 _distancetomaxcell_perthickperlayer_200_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_FH_zminus", [
273  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
274 ], ncols=4)
275 
276 _distancetomaxcell_perthickperlayer_200_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_BH_zminus", [
277  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
278 ], ncols=4)
279 
280 #300 um
281 _distancetomaxcell_perthickperlayer_300_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_EE_zminus", [
282  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
283 ], ncols=4)
284 
285 _distancetomaxcell_perthickperlayer_300_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_FH_zminus", [
286  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
287 ], ncols=4)
288 
289 _distancetomaxcell_perthickperlayer_300_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_BH_zminus", [
290  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
291 ], ncols=4)
292 
293 #scint um
294 _distancetomaxcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_EE_zminus", [
295  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
296 ], ncols=4)
297 
298 _distancetomaxcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_FH_zminus", [
299  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
300 ], ncols=4)
301 
302 _distancetomaxcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_BH_zminus", [
303  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
304 ], ncols=4)
305 
306 #----------------------------------------------------------------------------------------------------------------
307 #120 um
308 _distancebetseedandmaxcell_perthickperlayer_120_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_EE_zminus", [
309  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
310 ], ncols=4)
311 
312 _distancebetseedandmaxcell_perthickperlayer_120_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_FH_zminus", [
313  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
314 ], ncols=4)
315 
316 _distancebetseedandmaxcell_perthickperlayer_120_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_BH_zminus", [
317  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
318 ], ncols=4)
319 
320 #200 um
321 _distancebetseedandmaxcell_perthickperlayer_200_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_EE_zminus", [
322  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
323 ], ncols=4)
324 
325 _distancebetseedandmaxcell_perthickperlayer_200_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_FH_zminus", [
326  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
327 ], ncols=4)
328 
329 _distancebetseedandmaxcell_perthickperlayer_200_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_BH_zminus", [
330  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
331 ], ncols=4)
332 
333 #300 um
334 _distancebetseedandmaxcell_perthickperlayer_300_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_EE_zminus", [
335  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
336 ], ncols=4)
337 
338 _distancebetseedandmaxcell_perthickperlayer_300_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_FH_zminus", [
339  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
340 ], ncols=4)
341 
342 _distancebetseedandmaxcell_perthickperlayer_300_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_BH_zminus", [
343  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
344 ], ncols=4)
345 
346 #scint um
347 _distancebetseedandmaxcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_EE_zminus", [
348  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
349 ], ncols=4)
350 
351 _distancebetseedandmaxcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_FH_zminus", [
352  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
353 ], ncols=4)
354 
355 _distancebetseedandmaxcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_BH_zminus", [
356  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
357 ], ncols=4)
358 
359 #----------------------------------------------------------------------------------------------------------------
360 #120 um
361 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus", [
362  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
363 ], ncols=4)
364 
365 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus", [
366  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
367 ], ncols=4)
368 
369 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus", [
370  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
371 ], ncols=4)
372 
373 #200 um
374 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus", [
375  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
376 ], ncols=4)
377 
378 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus", [
379  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
380 ], ncols=4)
381 
382 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus", [
383  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
384 ], ncols=4)
385 
386 #300 um
387 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus", [
388  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
389 ], ncols=4)
390 
391 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus", [
392  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
393 ], ncols=4)
394 
395 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus", [
396  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
397 ], ncols=4)
398 
399 #scint um
400 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_EE_zminus", [
401  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
402 ], ncols=4)
403 
404 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_FH_zminus", [
405  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
406 ], ncols=4)
407 
408 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_BH_zminus", [
409  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
410 ], ncols=4)
411 
412 #----------------------------------------------------------------------------------------------------------------
413 #120 um
414 _distancetoseedcell_perthickperlayer_120_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_EE_zminus", [
415  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
416 ], ncols=4)
417 
418 _distancetoseedcell_perthickperlayer_120_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_FH_zminus", [
419  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
420 ], ncols=4)
421 
422 _distancetoseedcell_perthickperlayer_120_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_BH_zminus", [
423  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
424 ], ncols=4)
425 
426 #200 um
427 _distancetoseedcell_perthickperlayer_200_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_EE_zminus", [
428  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
429 ], ncols=4)
430 
431 _distancetoseedcell_perthickperlayer_200_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_FH_zminus", [
432  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
433 ], ncols=4)
434 
435 _distancetoseedcell_perthickperlayer_200_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_BH_zminus", [
436  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
437 ], ncols=4)
438 
439 #300 um
440 _distancetoseedcell_perthickperlayer_300_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_EE_zminus", [
441  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
442 ], ncols=4)
443 
444 _distancetoseedcell_perthickperlayer_300_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_FH_zminus", [
445  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
446 ], ncols=4)
447 
448 _distancetoseedcell_perthickperlayer_300_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_BH_zminus", [
449  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
450 ], ncols=4)
451 
452 #scint um
453 _distancetoseedcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_EE_zminus", [
454  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
455 ], ncols=4)
456 
457 _distancetoseedcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_FH_zminus", [
458  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
459 ], ncols=4)
460 
461 _distancetoseedcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_BH_zminus", [
462  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
463 ], ncols=4)
464 
465 #=====================================================================================================================
466 #----------------------------------------------------------------------------------------------------------------
467 #We need points for the weighted plots
468 _common = {"stat": True, "drawStyle": "EP", "staty": 0.65 }
469 #120 um
470 _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus", [
471  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
472 ], ncols=4)
473 
474 _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus", [
475  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
476 ], ncols=4)
477 
478 _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus", [
479  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
480 ], ncols=4)
481 
482 #200 um
483 _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus", [
484  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
485 ], ncols=4)
486 
487 _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus", [
488  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
489 ], ncols=4)
490 
491 _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus", [
492  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
493 ], ncols=4)
494 
495 #300 um
496 _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus", [
497  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
498 ], ncols=4)
499 
500 _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus", [
501  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
502 ], ncols=4)
503 
504 _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus", [
505  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
506 ], ncols=4)
507 #scint um
508 _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_EE_zminus", [
509  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
510 ], ncols=4)
511 
512 _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_FH_zminus", [
513  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
514 ], ncols=4)
515 
516 _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_BH_zminus", [
517  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
518 ], ncols=4)
519 
520 
521 #----------------------------------------------------------------------------------------------------------------
522 #120 um
523 _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus", [
524  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
525 ], ncols=4)
526 
527 _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus", [
528  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
529 ], ncols=4)
530 
531 _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus", [
532  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
533 ], ncols=4)
534 
535 #200 um
536 _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus", [
537  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
538 ], ncols=4)
539 
540 _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus", [
541  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
542 ], ncols=4)
543 
544 _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus", [
545  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
546 ], ncols=4)
547 
548 #300 um
549 _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus", [
550  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
551 ], ncols=4)
552 
553 _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus", [
554  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
555 ], ncols=4)
556 
557 _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus", [
558  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
559 ], ncols=4)
560 
561 #scint um
562 _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_EE_zminus", [
563  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
564 ], ncols=4)
565 
566 _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_FH_zminus", [
567  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
568 ], ncols=4)
569 
570 _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_BH_zminus", [
571  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
572 ], ncols=4)
573 
574 #Coming back to the usual definition
575 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
576 
577 #--------------------------------------------------------------------------------------------
578 # z+
579 #--------------------------------------------------------------------------------------------
580 _totclusternum_layer_EE_zplus = PlotGroup("totclusternum_layer_EE_zplus", [
581  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
582 ], ncols=4)
583 
584 _totclusternum_layer_FH_zplus = PlotGroup("totclusternum_layer_FH_zplus", [
585  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
586 ], ncols=4)
587 
588 _totclusternum_layer_BH_zplus = PlotGroup("totclusternum_layer_BH_zplus", [
589  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
590 ], ncols=4)
591 
592 _energyclustered_perlayer_EE_zplus = PlotGroup("energyclustered_perlayer_EE_zplus", [
593  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
594 ], ncols=4)
595 
596 _energyclustered_perlayer_FH_zplus = PlotGroup("energyclustered_perlayer_FH_zplus", [
597  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
598 ], ncols=4)
599 
600 _energyclustered_perlayer_BH_zplus = PlotGroup("energyclustered_perlayer_BH_zplus", [
601  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
602 ], ncols=4)
603 
604 #----------------------------------------------------------------------------------------------------------------
605 #120 um
606 _cellsnum_perthick_perlayer_120_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_120_EE_zplus", [
607  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
608 ], ncols=4)
609 
610 _cellsnum_perthick_perlayer_120_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_120_FH_zplus", [
611  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
612 ], ncols=4)
613 _cellsnum_perthick_perlayer_120_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_120_BH_zplus", [
614  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
615 ], ncols=4)
616 
617 #200 um
618 _cellsnum_perthick_perlayer_200_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_200_EE_zplus", [
619  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
620 ], ncols=4)
621 
622 _cellsnum_perthick_perlayer_200_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_200_FH_zplus", [
623  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
624 ], ncols=4)
625 
626 _cellsnum_perthick_perlayer_200_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_200_BH_zplus", [
627  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
628 ], ncols=4)
629 #300 um
630 _cellsnum_perthick_perlayer_300_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_300_EE_zplus", [
631  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
632 ], ncols=4)
633 
634 _cellsnum_perthick_perlayer_300_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_300_FH_zplus", [
635  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
636 ], ncols=4)
637 _cellsnum_perthick_perlayer_300_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_300_BH_zplus", [
638  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
639 ], ncols=4)
640 
641 #scint um
642 _cellsnum_perthick_perlayer_scint_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_EE_zplus", [
643  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
644 ], ncols=4)
645 
646 _cellsnum_perthick_perlayer_scint_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_FH_zplus", [
647  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
648 ], ncols=4)
649 
650 _cellsnum_perthick_perlayer_scint_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_BH_zplus", [
651  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
652 ], ncols=4)
653 
654 #----------------------------------------------------------------------------------------------------------------
655 #120 um
656 _common_distance = {}
657 _common_distance.update(_common)
658 _common_distance.update(_legend_common)
659 _common_distance["xmax"] = 150
660 _common_distance["stat"] = False
661 _common_distance["ymin"] = 1e-3
662 _common_distance["ymax"] = 10000
663 _common_distance["ylog"] = True
664 
665 _distancetomaxcell_perthickperlayer_120_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_EE_zplus", [
666  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
667 ], ncols=4)
668 
669 _distancetomaxcell_perthickperlayer_120_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_FH_zplus", [
670  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
671 ], ncols=4)
672 
673 _distancetomaxcell_perthickperlayer_120_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_BH_zplus", [
674  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
675 ], ncols=4)
676 
677 #200 um
678 _distancetomaxcell_perthickperlayer_200_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_EE_zplus", [
679  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
680 ], ncols=4)
681 
682 _distancetomaxcell_perthickperlayer_200_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_FH_zplus", [
683  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
684 ], ncols=4)
685 
686 _distancetomaxcell_perthickperlayer_200_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_BH_zplus", [
687  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
688 ], ncols=4)
689 
690 #300 um
691 _distancetomaxcell_perthickperlayer_300_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_EE_zplus", [
692  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
693 ], ncols=4)
694 
695 _distancetomaxcell_perthickperlayer_300_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_FH_zplus", [
696  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
697 ], ncols=4)
698 
699 _distancetomaxcell_perthickperlayer_300_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_BH_zplus", [
700  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
701 ], ncols=4)
702 
703 #scint um
704 _distancetomaxcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_EE_zplus", [
705  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
706 ], ncols=4)
707 
708 _distancetomaxcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_FH_zplus", [
709  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
710 ], ncols=4)
711 
712 _distancetomaxcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_BH_zplus", [
713  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
714 ], ncols=4)
715 
716 #----------------------------------------------------------------------------------------------------------------
717 #120 um
718 _distancebetseedandmaxcell_perthickperlayer_120_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_EE_zplus", [
719  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
720 ], ncols=4)
721 
722 _distancebetseedandmaxcell_perthickperlayer_120_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_FH_zplus", [
723  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
724 ], ncols=4)
725 
726 _distancebetseedandmaxcell_perthickperlayer_120_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_BH_zplus", [
727  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
728 ], ncols=4)
729 
730 #200 um
731 _distancebetseedandmaxcell_perthickperlayer_200_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_EE_zplus", [
732  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
733 ], ncols=4)
734 
735 _distancebetseedandmaxcell_perthickperlayer_200_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_FH_zplus", [
736  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
737 ], ncols=4)
738 
739 _distancebetseedandmaxcell_perthickperlayer_200_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_BH_zplus", [
740  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
741 ], ncols=4)
742 
743 #300 um
744 _distancebetseedandmaxcell_perthickperlayer_300_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_EE_zplus", [
745  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
746 ], ncols=4)
747 
748 _distancebetseedandmaxcell_perthickperlayer_300_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_FH_zplus", [
749  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
750 ], ncols=4)
751 
752 _distancebetseedandmaxcell_perthickperlayer_300_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_BH_zplus", [
753  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
754 ], ncols=4)
755 
756 #scint um
757 _distancebetseedandmaxcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_EE_zplus", [
758  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
759 ], ncols=4)
760 
761 _distancebetseedandmaxcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_FH_zplus", [
762  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
763 ], ncols=4)
764 
765 _distancebetseedandmaxcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_BH_zplus", [
766  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
767 ], ncols=4)
768 
769 #----------------------------------------------------------------------------------------------------------------
770 #120 um
771 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus", [
772  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
773 ], ncols=4)
774 
775 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus", [
776  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
777 ], ncols=4)
778 
779 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus", [
780  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
781 ], ncols=4)
782 
783 #200 um
784 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus", [
785  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
786 ], ncols=4)
787 
788 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus", [
789  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
790 ], ncols=4)
791 
792 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus", [
793  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
794 ], ncols=4)
795 
796 #300 um
797 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus", [
798  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
799 ], ncols=4)
800 
801 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus", [
802  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
803 ], ncols=4)
804 
805 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus", [
806  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
807 ], ncols=4)
808 
809 #scint um
810 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_EE_zplus", [
811  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
812 ], ncols=4)
813 
814 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_FH_zplus", [
815  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
816 ], ncols=4)
817 
818 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_BH_zplus", [
819  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
820 ], ncols=4)
821 
822 
823 #----------------------------------------------------------------------------------------------------------------
824 #120 um
825 _distancetoseedcell_perthickperlayer_120_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_EE_zplus", [
826  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
827 ], ncols=4)
828 
829 _distancetoseedcell_perthickperlayer_120_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_FH_zplus", [
830  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
831 ], ncols=4)
832 
833 _distancetoseedcell_perthickperlayer_120_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_BH_zplus", [
834  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
835 ], ncols=4)
836 
837 #200 um
838 _distancetoseedcell_perthickperlayer_200_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_EE_zplus", [
839  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
840 ], ncols=4)
841 
842 _distancetoseedcell_perthickperlayer_200_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_FH_zplus", [
843  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
844 ], ncols=4)
845 
846 _distancetoseedcell_perthickperlayer_200_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_BH_zplus", [
847  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
848 ], ncols=4)
849 
850 #300 um
851 _distancetoseedcell_perthickperlayer_300_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_EE_zplus", [
852  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
853 ], ncols=4)
854 
855 _distancetoseedcell_perthickperlayer_300_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_FH_zplus", [
856  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
857 ], ncols=4)
858 
859 _distancetoseedcell_perthickperlayer_300_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_BH_zplus", [
860  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
861 ], ncols=4)
862 
863 #scint um
864 _distancetoseedcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_EE_zplus", [
865  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
866 ], ncols=4)
867 
868 _distancetoseedcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_FH_zplus", [
869  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
870 ], ncols=4)
871 
872 _distancetoseedcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_BH_zplus", [
873  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
874 ], ncols=4)
875 
876 #=====================================================================================================================
877 #----------------------------------------------------------------------------------------------------------------
878 #We need points for the weighted plots
879 _common = {"stat": True, "drawStyle": "EP", "staty": 0.65 }
880 
881 #120 um
882 _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus", [
883  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
884 ], ncols=4)
885 
886 _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus", [
887  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
888 ], ncols=4)
889 
890 _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus", [
891  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
892 ], ncols=4)
893 
894 #200 um
895 _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus", [
896  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
897 ], ncols=4)
898 _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus", [
899  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
900 ], ncols=4)
901 
902 _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus", [
903  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
904 ], ncols=4)
905 
906 #300 um
907 _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus", [
908  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
909 ], ncols=4)
910 
911 _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus", [
912  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
913 ], ncols=4)
914 
915 _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus", [
916  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
917 ], ncols=4)
918 
919 #scint um
920 _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_EE_zplus", [
921  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
922 ], ncols=4)
923 
924 _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_FH_zplus", [
925  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
926 ], ncols=4)
927 
928 _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_BH_zplus", [
929  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
930 ], ncols=4)
931 
932 #----------------------------------------------------------------------------------------------------------------
933 #120 um
934 _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus", [
935  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
936 ], ncols=4)
937 
938 _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus", [
939  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
940 ], ncols=4)
941 
942 _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus", [
943  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
944 ], ncols=4)
945 
946 #200 um
947 _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus", [
948  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
949 ], ncols=4)
950 
951 _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus", [
952  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
953 ], ncols=4)
954 
955 _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus", [
956  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
957 ], ncols=4)
958 
959 #300 um
960 _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus", [
961  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
962 ], ncols=4)
963 
964 _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus", [
965  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
966 ], ncols=4)
967 
968 _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus", [
969  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
970 ], ncols=4)
971 
972 #scint um
973 _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_EE_zplus", [
974  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
975 ], ncols=4)
976 
977 _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_FH_zplus", [
978  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
979 ], ncols=4)
980 
981 _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_BH_zplus", [
982  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
983 ], ncols=4)
984 #Just in case we add some plots below to be on the safe side.
985 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
986 
987 #--------------------------------------------------------------------------------------------
988 # z-
989 #--------------------------------------------------------------------------------------------
990 
991 _common_score = {"title": "Score CaloParticle to LayerClusters in z-",
992  "stat": False,
993  "ymin": 0.1,
994  "ymax": 1000,
995  "xmin": 0,
996  "xmax": 1,
997  "drawStyle": "hist",
998  "lineWidth": 1,
999  "ylog": True
1000  }
1001 _common_score.update(_legend_common)
1002 _score_caloparticle_to_layerclusters_zminus = PlotGroup("score_caloparticle_to_layercluster_zminus", [
1003  Plot("Score_caloparticle2layercl_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_score) for i in range(0,maxlayerzm)
1004  ], ncols=8 )
1005 
1006 _common_score = {"title": "Score LayerCluster to CaloParticles in z-",
1007  "stat": False,
1008  "ymin": 0.1,
1009  "ymax": 1000,
1010  "xmin": 0,
1011  "xmax": 1,
1012  "drawStyle": "hist",
1013  "lineWidth": 1,
1014  "ylog": True
1015  }
1016 _common_score.update(_legend_common)
1017 _score_layercluster_to_caloparticles_zminus = PlotGroup("score_layercluster_to_caloparticle_zminus", [
1018  Plot("Score_layercl2caloparticle_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_score) for i in range(0,maxlayerzm)
1019  ], ncols=8 )
1020 
1021 _common_shared= {"title": "Shared Energy CaloParticle To Layer Cluster in z-",
1022  "stat": False,
1023  "legend": False,
1024  }
1025 _common_shared.update(_legend_common)
1026 _shared_plots_zminus = [Plot("SharedEnergy_caloparticle2layercl_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)]
1027 _shared_plots_zminus.extend([Plot("SharedEnergy_caloparticle2layercl_vs_eta_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)])
1028 _shared_plots_zminus.extend([Plot("SharedEnergy_caloparticle2layercl_vs_phi_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)])
1029 _sharedEnergy_caloparticle_to_layercluster_zminus = PlotGroup("sharedEnergy_caloparticle_to_layercluster_zminus", _shared_plots_zminus, ncols=8)
1030 
1031 _common_shared= {"title": "Shared Energy Layer Cluster To CaloParticle in z-",
1032  "stat": False,
1033  "legend": False,
1034  }
1035 _common_shared.update(_legend_common)
1036 _shared_plots2_zminus = [Plot("SharedEnergy_layercluster2caloparticle_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)]
1037 _common_shared= {"title": "Shared Energy Layer Cluster To Best CaloParticle in z-",
1038  "stat": False,
1039  "legend": False,
1040  "ymin": 0,
1041  "ymax": 1
1042  }
1043 _common_shared.update(_legend_common)
1044 _shared_plots2_zminus.extend([Plot("SharedEnergy_layercl2caloparticle_vs_eta_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)])
1045 _shared_plots2_zminus.extend([Plot("SharedEnergy_layercl2caloparticle_vs_phi_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(0,maxlayerzm)])
1046 _sharedEnergy_layercluster_to_caloparticle_zminus = PlotGroup("sharedEnergy_layercluster_to_caloparticle_zminus", _shared_plots2_zminus, ncols=8)
1047 
1048 
1049 _common_assoc = {#"title": "Cell Association Table in z-",
1050  "stat": False,
1051  "legend": False,
1052  "xbinlabels": ["", "TN(pur)", "FN(ineff.)", "FP(fake)", "TP(eff)"],
1053  "xbinlabeloption": "h",
1054  "drawStyle": "hist",
1055  "ymin": 0.1,
1056  "ymax": 10000,
1057  "ylog": True}
1058 _common_assoc.update(_legend_common)
1059 _cell_association_table_zminus = PlotGroup("cellAssociation_table_zminus", [
1060  Plot("cellAssociation_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_assoc) for i in range(0,maxlayerzm)
1061  ], ncols=8 )
1062 
1063 _bin_count = 0
1064 _xbinlabels = [ "L{:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1065 _common_eff = {"stat": False, "legend": False}
1066 _effplots_zminus = [Plot("effic_eta_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(0,maxlayerzm)]
1067 _effplots_zminus.extend([Plot("effic_phi_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(0,maxlayerzm)])
1068 _common_eff = {"stat": False, "legend": False, "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1069 _common_eff["xmin"] = _bin_count
1070 _common_eff["xmax"] = maxlayerzm
1071 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1072 _effplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Efficiencies in z-", **_common_eff)])
1073 _efficiencies_zminus = PlotGroup("Efficiencies_zminus", _effplots_zminus, ncols=8)
1074 
1075 _common_dup = {"stat": False, "legend": False}
1076 _dupplots_zminus = [Plot("duplicate_eta_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(0,maxlayerzm)]
1077 _dupplots_zminus.extend([Plot("duplicate_phi_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(0,maxlayerzm)])
1078 _common_dup = {"stat": False, "legend": False, "title": "Global Duplicates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1079 _common_dup["xmin"] = _bin_count
1080 _common_dup["xmax"] = _common_dup["xmin"] + maxlayerzm
1081 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1082 _dupplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Duplicates in z-", **_common_dup)])
1083 _duplicates_zminus = PlotGroup("Duplicates_zminus", _dupplots_zminus, ncols=8)
1084 
1085 _common_fake = {"stat": False, "legend": False}
1086 _fakeplots_zminus = [Plot("fake_eta_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(0,maxlayerzm)]
1087 _fakeplots_zminus.extend([Plot("fake_phi_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(0,maxlayerzm)])
1088 _common_fake = {"stat": False, "legend": False, "title": "Global Fake Rates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1089 _common_fake["xmin"] = _bin_count
1090 _common_fake["xmax"] = _common_fake["xmin"] + maxlayerzm
1091 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1092 _common_fake["xbinlabels"] = [ "L{:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1093 _common_fake["xbinlabelsize"] = 10.
1094 _fakeplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Fake Rate in z-", **_common_fake)])
1095 _fakes_zminus = PlotGroup("FakeRate_zminus", _fakeplots_zminus, ncols=8)
1096 
1097 _common_merge = {"stat": False, "legend": False}
1098 _mergeplots_zminus = [Plot("merge_eta_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(0,maxlayerzm)]
1099 _mergeplots_zminus.extend([Plot("merge_phi_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(0,maxlayerzm)])
1100 _common_merge = {"stat": False, "legend": False, "title": "Global Merge Rates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1101 _common_merge["xmin"] = _bin_count
1102 _common_merge["xmax"] = _common_merge["xmin"] + maxlayerzm
1103 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1104 _common_merge["xbinlabels"] = [ "L{:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1105 _common_merge["xbinlabelsize"] = 10.
1106 _mergeplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global merge Rate in z-", **_common_merge)])
1107 _merges_zminus = PlotGroup("MergeRate_zminus", _mergeplots_zminus, ncols=8)
1108 
1109 
1110 _common_energy_score = dict(removeEmptyBins=False, xbinlabelsize=10,
1111  stat=True,
1112  xbinlabeloption="d",
1113  ncols=1,
1114  ylog=True,
1115  xlog=True,
1116  xmin=0.001,
1117  xmax=1.,
1118  ymin=0.01,
1119  ymax=1.)
1120 _energyscore_cp2lc_zminus = []
1121 for i in range(0, maxlayerzm):
1122  _energyscore_cp2lc_zminus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_caloparticle2layer_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1))
1123 
1124 _energyscore_lc2cp_zminus = []
1125 _common_energy_score["xlog"]=False
1126 _common_energy_score["ylog"]=False
1127 _common_energy_score["xmin"]=-0.1
1128 for i in range(0, maxlayerzm):
1129  _energyscore_lc2cp_zminus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_layer2caloparticle_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1))
1130 #_energyclustered =
1131 
1132 #--------------------------------------------------------------------------------------------
1133 # z+
1134 #--------------------------------------------------------------------------------------------
1135 _common_score = {"title": "Score CaloParticle to LayerClusters in z+",
1136  "stat": False,
1137  "ymin": 0.1,
1138  "ymax": 1000,
1139  "xmin": 0,
1140  "xmax": 1,
1141  "drawStyle": "hist",
1142  "lineWidth": 1,
1143  "ylog": True
1144  }
1145 _common_score.update(_legend_common)
1146 _score_caloparticle_to_layerclusters_zplus = PlotGroup("score_caloparticle_to_layercluster_zplus", [
1147  Plot("Score_caloparticle2layercl_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_score) for i in range(maxlayerzm,maxlayerzp)
1148  ], ncols=8 )
1149 
1150 _common_score = {"title": "Score LayerCluster to CaloParticles in z+",
1151  "stat": False,
1152  "ymin": 0.1,
1153  "ymax": 1000,
1154  "xmin": 0,
1155  "xmax": 1,
1156  "drawStyle": "hist",
1157  "lineWidth": 1,
1158  "ylog": True
1159  }
1160 _common_score.update(_legend_common)
1161 _score_layercluster_to_caloparticles_zplus = PlotGroup("score_layercluster_to_caloparticle_zplus", [
1162  Plot("Score_layercl2caloparticle_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_score) for i in range(maxlayerzm,maxlayerzp)
1163  ], ncols=8 )
1164 
1165 _common_shared= {"title": "Shared Energy CaloParticle To Layer Cluster in z+",
1166  "stat": False,
1167  "legend": False,
1168  }
1169 _common_shared.update(_legend_common)
1170 _shared_plots_zplus = [Plot("SharedEnergy_caloparticle2layercl_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)]
1171 _shared_plots_zplus.extend([Plot("SharedEnergy_caloparticle2layercl_vs_eta_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)])
1172 _shared_plots_zplus.extend([Plot("SharedEnergy_caloparticle2layercl_vs_phi_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)])
1173 _sharedEnergy_caloparticle_to_layercluster_zplus = PlotGroup("sharedEnergy_caloparticle_to_layercluster_zplus", _shared_plots_zplus, ncols=8)
1174 
1175 _common_shared= {"title": "Shared Energy Layer Cluster To CaloParticle in z+",
1176  "stat": False,
1177  "legend": False,
1178  }
1179 _common_shared.update(_legend_common)
1180 _shared_plots2_zplus = [Plot("SharedEnergy_layercluster2caloparticle_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)]
1181 _common_shared= {"title": "Shared Energy Layer Cluster To Best CaloParticle in z+",
1182  "stat": False,
1183  "legend": False,
1184  "ymin": 0,
1185  "ymax": 1,
1186  }
1187 _common_shared.update(_legend_common)
1188 _shared_plots2_zplus.extend([Plot("SharedEnergy_layercl2caloparticle_vs_eta_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)])
1189 _shared_plots2_zplus.extend([Plot("SharedEnergy_layercl2caloparticle_vs_phi_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_shared) for i in range(maxlayerzm,maxlayerzp)])
1190 _sharedEnergy_layercluster_to_caloparticle_zplus = PlotGroup("sharedEnergy_layercluster_to_caloparticle_zplus", _shared_plots2_zplus, ncols=8)
1191 
1192 
1193 _common_assoc = {#"title": "Cell Association Table in z+",
1194  "stat": False,
1195  "legend": False,
1196  "xbinlabels": ["", "TN(pur)", "FN(ineff.)", "FP(fake)", "TP(eff)"],
1197  "xbinlabeloption": "h",
1198  "drawStyle": "hist",
1199  "ymin": 0.1,
1200  "ymax": 10000,
1201  "ylog": True}
1202 _common_assoc.update(_legend_common)
1203 _cell_association_table_zplus = PlotGroup("cellAssociation_table_zplus", [
1204  Plot("cellAssociation_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_assoc) for i in range(maxlayerzm,maxlayerzp)
1205  ], ncols=8 )
1206 
1207 
1208 _bin_count = 50
1209 _common_eff = {"stat": False, "legend": False}
1210 _effplots_zplus = [Plot("effic_eta_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(maxlayerzm,maxlayerzp)]
1211 _effplots_zplus.extend([Plot("effic_phi_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(maxlayerzm,maxlayerzp)])
1212 _common_eff = {"stat": False, "legend": False, "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1213 _common_eff["xmin"] = _bin_count
1214 _common_eff["xmax"] = _common_eff["xmin"] + maxlayerzm
1215 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1216 _effplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Efficiencies in z+", **_common_eff)])
1217 _efficiencies_zplus = PlotGroup("Efficiencies_zplus", _effplots_zplus, ncols=8)
1218 
1219 
1220 _common_dup = {"stat": False, "legend": False}
1221 _dupplots_zplus = [Plot("duplicate_eta_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(maxlayerzm,maxlayerzp)]
1222 _dupplots_zplus.extend([Plot("duplicate_phi_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(maxlayerzm,maxlayerzp)])
1223 _common_dup = {"stat": False, "legend": False, "title": "Global Duplicates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1224 _common_dup["xmin"] = _bin_count
1225 _common_dup["xmax"] = _common_dup["xmin"] + maxlayerzm
1226 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1227 _dupplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Duplicates in z+", **_common_dup)])
1228 _duplicates_zplus = PlotGroup("Duplicates_zplus", _dupplots_zplus, ncols=8)
1229 
1230 _common_fake = {"stat": False, "legend": False}
1231 _fakeplots_zplus = [Plot("fake_eta_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(maxlayerzm,maxlayerzp)]
1232 _fakeplots_zplus.extend([Plot("fake_phi_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(maxlayerzm,maxlayerzp)])
1233 _common_fake = {"stat": False, "legend": False, "title": "Global Fake Rates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1234 _common_fake["xmin"] = _bin_count
1235 _common_fake["xmax"] = _common_fake["xmin"] + maxlayerzm
1236 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1237 _fakeplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Fake Rate in z+", **_common_fake)])
1238 _fakes_zplus = PlotGroup("FakeRate_zplus", _fakeplots_zplus, ncols=8)
1239 
1240 _common_merge = {"stat": False, "legend": False}
1241 _mergeplots_zplus = [Plot("merge_eta_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(maxlayerzm,maxlayerzp)]
1242 _mergeplots_zplus.extend([Plot("merge_phi_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(maxlayerzm,maxlayerzp)])
1243 _common_merge = {"stat": False, "legend": False, "title": "Global Merge Rates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1244 _common_merge["xmin"] = _bin_count
1245 _common_merge["xmax"] = _common_merge["xmin"] + maxlayerzm
1246 _bin_count += 4*maxlayerzm # 2 for the eta{-,+} and 2 for phi{+,-}
1247 _mergeplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global merge Rate in z+", **_common_merge)])
1248 _merges_zplus = PlotGroup("MergeRate_zplus", _mergeplots_zplus, ncols=8)
1249 
1250 
1251 _common_energy_score = dict(removeEmptyBins=False, xbinlabelsize=10,
1252  stat=True,
1253  xbinlabeloption="d",
1254  ncols=1,
1255  ylog=True,
1256  xlog=True,
1257  xmin=0.001,
1258  xmax=1.,
1259  ymin=0.01,
1260  ymax=1.)
1261 _energyscore_cp2lc_zplus = []
1262 for i in range(maxlayerzm,maxlayerzp):
1263  _energyscore_cp2lc_zplus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_caloparticle2layer_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1))
1264 
1265 _common_energy_score["xlog"]=False
1266 _common_energy_score["ylog"]=False
1267 _common_energy_score["xmin"]=-0.1
1268 _energyscore_lc2cp_zplus = []
1269 for i in range(maxlayerzm,maxlayerzp):
1270  _energyscore_lc2cp_zplus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_layer2caloparticle_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1))
1271 #_energyclustered =
1272 
1273 #--------------------------------------------------------------------------------------------
1274 # MULTICLUSTERS
1275 #--------------------------------------------------------------------------------------------
1276 _common_score = {#"title": "Score CaloParticle to MultiClusters",
1277  "stat": False,
1278  "ymin": 0.1,
1279  "ymax": 100000,
1280  "xmin": 0,
1281  "xmax": 1.0,
1282  "drawStyle": "hist",
1283  "lineWidth": 1,
1284  "ylog": True
1285  }
1286 _common_score.update(_legend_common)
1287 _score_caloparticle_to_multiclusters = PlotGroup("score_caloparticle_to_multicluster", [
1288  Plot("Score_caloparticle2multicl", **_common_score)
1289  ], ncols=1)
1290 
1291 _common_score = {#"title": "Score MultiCluster to CaloParticles",
1292  "stat": False,
1293  "ymin": 0.1,
1294  "ymax": 100000,
1295  "xmin": 0,
1296  "xmax": 1.0,
1297  "drawStyle": "hist",
1298  "lineWidth": 1,
1299  "ylog": True
1300  }
1301 _common_score.update(_legend_common)
1302 _score_multicluster_to_caloparticles = PlotGroup("score_multicluster_to_caloparticle", [
1303  Plot("Score_multicl2caloparticle", **_common_score)
1304  ])
1305 
1306 _common_shared= {"title": "Shared Energy CaloParticle To Multi Cluster ",
1307  "stat": False,
1308  "legend": True,
1309  "xmin": 0,
1310  "xmax": 2.0,
1311  }
1312 _common_shared.update(_legend_common)
1313 _shared_plots = [ Plot("SharedEnergy_caloparticle2multicl", **_common_shared) ]
1314 _common_shared["xmin"] = -4.0
1315 _common_shared["xmax"] = 4.0
1316 _shared_plots.extend([Plot("SharedEnergy_caloparticle2multicl_vs_eta", **_common_shared)])
1317 _shared_plots.extend([Plot("SharedEnergy_caloparticle2multicl_vs_phi", **_common_shared)])
1318 _sharedEnergy_caloparticle_to_multicluster = PlotGroup("sharedEnergy_caloparticle_to_multicluster", _shared_plots, ncols=3)
1319 
1320 _common_shared= {"title": "Shared Energy Multi Cluster To CaloParticle ",
1321  "stat": False,
1322  "legend": True,
1323  "xmin": 0,
1324  "xmax": 2.0,
1325  }
1326 _common_shared.update(_legend_common)
1327 _shared_plots2 = [Plot("SharedEnergy_multicluster2caloparticle", **_common_shared)]
1328 _common_shared["xmin"] = -4.0
1329 _common_shared["xmax"] = 4.0
1330 _shared_plots2.extend([Plot("SharedEnergy_multicl2caloparticle_vs_eta", **_common_shared)])
1331 _shared_plots2.extend([Plot("SharedEnergy_multicl2caloparticle_vs_phi", **_common_shared)])
1332 _sharedEnergy_multicluster_to_caloparticle = PlotGroup("sharedEnergy_multicluster_to_caloparticle", _shared_plots2, ncols=3)
1333 
1334 
1335 _common_assoc = {#"title": "Cell Association Table",
1336  "stat": False,
1337  "legend": False,
1338  "xbinlabels": ["", "TN(pur)", "FN(ineff.)", "FP(fake)", "TP(eff)"],
1339  "xbinlabeloption": "h",
1340  "drawStyle": "hist",
1341  "ymin": 0.1,
1342  "ymax": 10000000,
1343  "ylog": True}
1344 _common_assoc.update(_legend_common)
1345 _cell_association_table = PlotGroup("cellAssociation_table", [
1346  Plot("cellAssociation_perlayer{:02d}".format(i), xtitle="Layer {:02d} in z-".format(i%maxlayerzm+1) if (i<maxlayerzm) else "Layer {:02d} in z+".format(i%maxlayerzm+1), **_common_assoc) for i in range(0,maxlayerzm)
1347  ], ncols=8 )
1348 
1349 _common_eff = {"stat": False, "legend": False}
1350 _effplots = [Plot("effic_eta", xtitle="", **_common_eff)]
1351 _effplots.extend([Plot("effic_phi", xtitle="", **_common_eff)])
1352 _effplots.extend([Plot("globalEfficiencies", xtitle="", **_common_eff)])
1353 _efficiencies = PlotGroup("Efficiencies", _effplots, ncols=3)
1354 
1355 
1356 _common_dup = {"stat": False, "legend": False}
1357 _dupplots = [Plot("duplicate_eta", xtitle="", **_common_dup)]
1358 _dupplots.extend([Plot("duplicate_phi", xtitle="", **_common_dup)])
1359 _dupplots.extend([Plot("globalEfficiencies", xtitle="", **_common_dup)])
1360 _duplicates = PlotGroup("Duplicates", _dupplots, ncols=3)
1361 
1362 _common_fake = {"stat": False, "legend": False}
1363 _fakeplots = [Plot("fake_eta", xtitle="", **_common_fake)]
1364 _fakeplots.extend([Plot("fake_phi", xtitle="", **_common_fake)])
1365 _fakeplots.extend([Plot("globalEfficiencies", xtitle="", **_common_fake)])
1366 _fakes = PlotGroup("FakeRate", _fakeplots, ncols=3)
1367 
1368 _common_merge = {"stat": False, "legend": False}
1369 _mergeplots = [Plot("merge_eta", xtitle="", **_common_merge)]
1370 _mergeplots.extend([Plot("merge_phi", xtitle="", **_common_merge)])
1371 _mergeplots.extend([Plot("globalEfficiencies", xtitle="", **_common_merge)])
1372 _merges = PlotGroup("MergeRate", _mergeplots, ncols=3)
1373 
1374 _common_energy_score = dict(removeEmptyBins=True, xbinlabelsize=10, xbinlabeloption="d")
1375 _common_energy_score["ymax"] = 1.
1376 _common_energy_score["xmax"] = 1.0
1377 _energyscore_cp2mcl = PlotOnSideGroup("_energyscore_cp2mcl", Plot("Energy_vs_Score_caloparticle2multi", drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1)
1378 _common_energy_score["ymax"] = 1.
1379 _common_energy_score["xmax"] = 1.0
1380 _energyscore_mcl2cp = PlotOnSideGroup("_energyscore_mcl2cp", Plot("Energy_vs_Score_multi2caloparticle", drawStyle="COLZ", adjustMarginRight=0.1, **_common_energy_score), ncols=1)
1381 
1382 #Coming back to the usual box definition
1383 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
1384 
1385 _totmulticlusternum = PlotGroup("totmulticlusternum", [
1386  Plot("totmulticlusternum", xtitle="", **_common)
1387 ],ncols=1)
1388 
1389 _multicluster_firstlayer = PlotGroup("multicluster_firstlayer", [
1390  Plot("multicluster_firstlayer", xtitle="Layer number", **_common)
1391 ],ncols=1)
1392 
1393 _multicluster_lastlayer = PlotGroup("multicluster_lastlayer", [
1394  Plot("multicluster_lastlayer", xtitle="Layer number", **_common)
1395 ],ncols=1)
1396 
1397 _multicluster_layersnum = PlotGroup("multicluster_layersnum", [
1398  Plot("multicluster_layersnum", xtitle="", **_common)
1399 ],ncols=1)
1400 
1401 _common["xmax"] = 50
1402 _clusternum_in_multicluster = PlotGroup("clusternum_in_multicluster",[
1403  Plot("clusternum_in_multicluster", xtitle="", **_common)
1404 ],ncols=1)
1405 
1406 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
1407 _common = {"stat": True, "drawStyle": "pcolz", "staty": 0.65}
1408 
1409 _clusternum_in_multicluster_vs_layer = PlotGroup("clusternum_in_multicluster_vs_layer",[
1410  Plot("clusternum_in_multicluster_vs_layer", xtitle="Layer number", ytitle = "<2d Layer Clusters in Multicluster>", **_common)
1411 ],ncols=1)
1412 
1413 _common["scale"] = 100.
1414 #, ztitle = "% of clusters" normalizeToUnitArea=True
1415 _multiplicity_numberOfEventsHistogram = "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA/multiplicity_numberOfEventsHistogram"
1416 _multiplicity_zminus_numberOfEventsHistogram = "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA/multiplicity_zminus_numberOfEventsHistogram"
1417 _multiplicity_zplus_numberOfEventsHistogram = "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA/multiplicity_zplus_numberOfEventsHistogram"
1418 _multiplicityOfLCinMCL = PlotGroup("multiplicityOfLCinMCL",[
1419  Plot("multiplicityOfLCinMCL", xtitle="Layer Cluster multiplicity in Multiclusters", ytitle = "Cluster size (n_{hit}) ", drawCommand = "colz text45", normalizeToNumberOfEvents = True, **_common)
1420 ],ncols=1)
1421 
1422 _multiplicityOfLCinMCL_vs_layercluster_zminus = PlotGroup("multiplicityOfLCinMCL_vs_layercluster_zminus",[
1423  Plot("multiplicityOfLCinMCL_vs_layercluster_zminus", xtitle="Layer Cluster multiplicity in Multiclusters", ytitle = "Layer Number ", drawCommand = "colz text45", normalizeToNumberOfEvents = True, **_common)
1424 ],ncols=1)
1425 
1426 _multiplicityOfLCinMCL_vs_layercluster_zplus = PlotGroup("multiplicityOfLCinMCL_vs_layercluster_zplus",[
1427  Plot("multiplicityOfLCinMCL_vs_layercluster_zplus", xtitle="Layer Cluster multiplicity in Multiclusters", ytitle = "Layer Number ", drawCommand = "colz text45", normalizeToNumberOfEvents = True, **_common)
1428 ],ncols=1)
1429 
1430 _multiplicityOfLCinMCL_vs_layerclusterenergy = PlotGroup("multiplicityOfLCinMCL_vs_layerclusterenergy",[
1431  Plot("multiplicityOfLCinMCL_vs_layerclusterenergy", xtitle="Layer Cluster multiplicity in Multiclusters", ytitle = "Cluster Energy (GeV) ", drawCommand = "colz text45", normalizeToNumberOfEvents = True, **_common)
1432 ],ncols=1)
1433 
1434 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
1435 #--------------------------------------------------------------------------------------------
1436 # z-
1437 #--------------------------------------------------------------------------------------------
1438 _clusternum_in_multicluster_perlayer_zminus_EE = PlotGroup("clusternum_in_multicluster_perlayer_zminus_EE", [
1439  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
1440 ], ncols=4)
1441 
1442 _clusternum_in_multicluster_perlayer_zminus_FH = PlotGroup("clusternum_in_multicluster_perlayer_zminus_FH", [
1443  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
1444 ], ncols=4)
1445 
1446 _clusternum_in_multicluster_perlayer_zminus_BH = PlotGroup("clusternum_in_multicluster_perlayer_zminus_BH", [
1447  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
1448 ], ncols=4)
1449 
1450 #--------------------------------------------------------------------------------------------
1451 # z+
1452 #--------------------------------------------------------------------------------------------
1453 _clusternum_in_multicluster_perlayer_zplus_EE = PlotGroup("clusternum_in_multicluster_perlayer_zplus_EE", [
1454  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
1455 ], ncols=4)
1456 
1457 _clusternum_in_multicluster_perlayer_zplus_FH = PlotGroup("clusternum_in_multicluster_perlayer_zplus_FH", [
1458  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
1459 ], ncols=4)
1460 
1461 _clusternum_in_multicluster_perlayer_zplus_BH = PlotGroup("clusternum_in_multicluster_perlayer_zplus_BH", [
1462  Plot("clusternum_in_multicluster_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
1463 ], ncols=4)
1464 
1465 #Coming back to the usual box definition
1466 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
1467 
1468 #Some multiclusters quantities
1469 _multicluster_pt = PlotGroup("multicluster_pt", [
1470  Plot("multicluster_pt", xtitle="", **_common)
1471 ],ncols=1)
1472 
1473 _multicluster_eta = PlotGroup("multicluster_eta", [
1474  Plot("multicluster_eta", xtitle="", **_common)
1475 ],ncols=1)
1476 
1477 _multicluster_phi = PlotGroup("multicluster_phi", [
1478  Plot("multicluster_phi", xtitle="", **_common)
1479 ],ncols=1)
1480 
1481 _multicluster_energy = PlotGroup("multicluster_energy", [
1482  Plot("multicluster_energy", xtitle="", **_common)
1483 ],ncols=1)
1484 
1485 _multicluster_x = PlotGroup("multicluster_x", [
1486  Plot("multicluster_x", xtitle="", **_common)
1487 ],ncols=1)
1488 
1489 _multicluster_y = PlotGroup("multicluster_y", [
1490  Plot("multicluster_y", xtitle="", **_common)
1491 ],ncols=1)
1492 
1493 _multicluster_z = PlotGroup("multicluster_z", [
1494  Plot("multicluster_z", xtitle="", **_common)
1495 ],ncols=1)
1496 
1497 #--------------------------------------------------------------------------------------------
1498 # SIMHITS, DIGIS, RECHITS
1499 #--------------------------------------------------------------------------------------------
1500 
1501 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ymin": 0.1, "ylog": True}
1502 
1503 _HitValidation = PlotGroup("HitValidation", [
1504  Plot("heeEnSim", title="SimHits_EE_Energy", **_common),
1505  Plot("hebEnSim", title="SimHits_HE_Silicon_Energy", **_common),
1506  Plot("hefEnSim", title="SimHits_HE_Scintillator_Energy", **_common),
1507  ])
1508 
1509 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
1510 
1511 _SimHits_Occupancy_EE_zplus = PlotGroup("SimHits_Occupancy_EE_zplus", [
1512  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_Occupancy_EE_zplus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1513  ], ncols=4)
1514 
1515 _SimHits_Occupancy_HE_Silicon_zplus = PlotGroup("SimHits_Occupancy_HE_Silicon_zplus", [
1516  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_HE_Occupancy_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1517  ], ncols=4)
1518 
1519 _SimHits_Occupancy_HE_Scintillator_zplus = PlotGroup("SimHits_Occupancy_HE_Scintillator_zplus", [
1520  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_Occupancy_HE_Scintillator_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1521  ], ncols=4)
1522 
1523 _SimHits_Occupancy_EE_zminus = PlotGroup("SimHits_Occupancy_EE_zminus", [
1524  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_Occupancy_EE_zminus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1525  ], ncols=4)
1526 
1527 _SimHits_Occupancy_HE_Silicon_zminus = PlotGroup("SimHits_Occupancy_HE_Silicon_zminus", [
1528  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_Occupancy_HE_Silicon_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1529  ], ncols=4)
1530 
1531 _SimHits_Occupancy_HE_Scintillator_zminus = PlotGroup("SimHits_Occupancy_HE_Scintillator_zminus", [
1532  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="SimHits_Occupancy_HE_Scintillator_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1533  ], ncols=4)
1534 
1535 _common_etaphi = dict(removeEmptyBins=False, xbinlabelsize=10, xbinlabeloption="d", ymin=None)
1536 
1537 _SimHits_EtaPhi_EE_zplus=[]
1538 for i in range(EE_min,EE_max+1):
1539  _SimHits_EtaPhi_EE_zplus.append(PlotOnSideGroup("EE_EtaPhi_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="SimHits_EtaPhi_EE_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1540 
1541 _SimHits_EtaPhi_HE_Silicon_zplus=[]
1542 for i in range(HESilicon_min,HESilicon_max+1):
1543  _SimHits_EtaPhi_HE_Silicon_zplus.append(PlotOnSideGroup("HE_Silicon_EtaPhi_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="SimHits_EtaPhi_HE_Silicon_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1544 
1545 _SimHits_EtaPhi_HE_Scintillator_zplus=[]
1546 for i in range(HEScintillator_min,HEScintillator_max+1):
1547  _SimHits_EtaPhi_HE_Scintillator_zplus.append(PlotOnSideGroup("HE_Scintillator_EtaPhi_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="SimHits_EtaPhi_HE_Scintillator_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1548 
1549 _SimHits_EtaPhi_EE_zminus=[]
1550 for i in range(EE_min,EE_max+1):
1551  _SimHits_EtaPhi_EE_zminus.append(PlotOnSideGroup("EE_EtaPhi_zminus_layer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="SimHits_EtaPhi_EE_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1552 
1553 _SimHits_EtaPhi_HE_Silicon_zminus=[]
1554 for i in range(HESilicon_min,HESilicon_max+1):
1555  _SimHits_EtaPhi_HE_Silicon_zminus.append(PlotOnSideGroup("HE_Silicon_EtaPhi_zminus_layer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="SimHits_EtaPhi_HE_Silicon_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1556 
1557 _SimHits_EtaPhi_HE_Scintillator_zminus=[]
1558 for i in range(HEScintillator_min,HEScintillator_max+1):
1559  _SimHits_EtaPhi_HE_Scintillator_zminus.append(PlotOnSideGroup("HE_Scintillator_EtaPhi_zminuslayer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="SimHits_EtaPhi_HE_Scintillator_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1560 
1561 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ymin": 0.1, "ylog": True}
1562 
1563 _SimHits_Energy_EE_0 = PlotGroup("SimHits_Energy_Time_0_EE", [
1564  Plot("energy_time_0_layer_{:02d}".format(i), title="SimHits_Energy_Time_0_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1565  ], ncols=4)
1566 
1567 _SimHits_Energy_HE_Silicon_0 = PlotGroup("SimHits_Energy_Time_0_HE_Silicon", [
1568  Plot("energy_time_0_layer_{:02d}".format(i), title="SimHits_Energy_Time_0_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1569  ], ncols=4)
1570 
1571 _SimHits_Energy_HE_Scintillator_0 = PlotGroup("SimHits_Energy_Time_0_HE_Scintillator", [
1572  Plot("energy_time_0_layer_{:02d}".format(i), title="SimHits_Energy_Time_0_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1573  ], ncols=4)
1574 
1575 _SimHits_Energy_EE_1 = PlotGroup("SimHits_Energy_Time_1_EE", [
1576  Plot("energy_time_1_layer_{:02d}".format(i), title="SimHits_Energy_Time_1_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1577  ], ncols=4)
1578 
1579 _SimHits_Energy_HE_Silicon_1 = PlotGroup("SimHits_Energy_Time_1_HE_Silicon", [
1580  Plot("energy_time_1_layer_{:02d}".format(i), title="SimHits_Energy_Time_1_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1581  ], ncols=4)
1582 
1583 _SimHits_Energy_HE_Scintillator_1 = PlotGroup("SimHits_Energy_Time_1_HE_Scintillator", [
1584  Plot("energy_time_1_layer_{:02d}".format(i), title="SimHits_Energy_Time_1_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1585  ], ncols=4)
1586 
1587 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
1588 
1589 _RecHits_Occupancy_EE_zplus = PlotGroup("RecHits_Occupancy_EE_zplus", [
1590  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="RecHits_Occupancy_EE_zplus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1591  ], ncols=4)
1592 
1593 _RecHits_Occupancy_HE_Silicon_zplus = PlotGroup("RecHits_Occupancy_HE_Silicon_zplus", [
1594  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="RecHits_Occupancy_HE_Silicon_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1595  ], ncols=4)
1596 
1597 _RecHits_Occupancy_HE_Scintillator_zplus = PlotGroup("RecHits_Occupancy_HE_Scintillator_zplus", [
1598  Plot("HitOccupancy_Plus_layer_{:02d}".format(i), title="RecHits_Occupancy_HE_Scintillator_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1599  ], ncols=4)
1600 
1601 _RecHits_Occupancy_EE_zminus = PlotGroup("RecHits_Occupancy_EE_zminus", [
1602  Plot("HitOccupancy_Minus_layer_{:02d}".format(i), title="RecHits_Occupancy_EE_zminus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1603  ], ncols=4)
1604 
1605 _RecHits_Occupancy_HE_Silicon_zminus = PlotGroup("RecHits_Occupancy_HE_Silicon_zminus", [
1606  Plot("HitOccupancy_Minus_layer_{:02d}".format(i), title="RecHits_Occupancy_HE_Silicon_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1607  ], ncols=4)
1608 
1609 _RecHits_Occupancy_HE_Scintillator_zminus = PlotGroup("RecHits_Occupancy_HE_Scintillator_zminus", [
1610  Plot("HitOccupancy_Minus_layer_{:02d}".format(i), title="RecHits_Occupancy_HE_Scintillator_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1611  ], ncols=4)
1612 
1613 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ymin": 0.1, "ylog": True}
1614 
1615 _RecHits_Energy_EE = PlotGroup("RecHits_Energy_EE", [
1616  Plot("energy_layer_{:02d}".format(i), title="RecHits_Energy_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1617  ], ncols=4)
1618 
1619 _RecHits_Energy_HE_Silicon = PlotGroup("RecHits_Energy_HE_Silicon", [
1620  Plot("energy_layer_{:02d}".format(i), title="RecHits_Energy_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1621  ], ncols=4)
1622 
1623 _RecHits_Energy_HE_Scintillator = PlotGroup("RecHits_Energy_HE_Scintillator", [
1624  Plot("energy_layer_{:02d}".format(i), title="RecHits_Energy_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1625  ], ncols=4)
1626 
1627 _RecHits_EtaPhi_EE_zplus=[]
1628 for i in range(EE_min,EE_max+1):
1629  _RecHits_EtaPhi_EE_zplus.append(PlotOnSideGroup("RecHits_EtaPhi_EE_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="RecHits_EtaPhi_EE_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1630 
1631 _RecHits_EtaPhi_HE_Silicon_zplus=[]
1632 for i in range(HESilicon_min,HESilicon_max+1):
1633  _RecHits_EtaPhi_HE_Silicon_zplus.append(PlotOnSideGroup("RecHits_EtaPhi_HE_Silicon_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="RecHits_EtaPhi_HE_Silicon_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1634 
1635 _RecHits_EtaPhi_HE_Scintillator_zplus=[]
1636 for i in range(HEScintillator_min,HEScintillator_max+1):
1637  _RecHits_EtaPhi_HE_Scintillator_zplus.append(PlotOnSideGroup("RecHits_EtaPhi_HE_Scintillator_zplus_layer_{:02d}".format(i), Plot("EtaPhi_Plus_layer_{:02d}".format(i), title="RecHits_EtaPhi_HE_Scintillator_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1638 
1639 _RecHits_EtaPhi_EE_zminus=[]
1640 for i in range(EE_min,EE_max+1):
1641  _RecHits_EtaPhi_EE_zminus.append(PlotOnSideGroup("RecHits_EtaPhi_EE_zminus_layer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="RecHits_EtaPhi_EE_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1642 
1643 _RecHits_EtaPhi_HE_Silicon_zminus=[]
1644 for i in range(HESilicon_min,HESilicon_max+1):
1645  _RecHits_EtaPhi_HE_Silicon_zminus.append(PlotOnSideGroup("RecHits_EtaPhi_HE_Silicon_zminus_layer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="RecHits_EtaPhi_HE_Silicon_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1646 
1647 _RecHits_EtaPhi_HE_Scintillator_zminus=[]
1648 for i in range(HEScintillator_min,HEScintillator_max+1):
1649  _RecHits_EtaPhi_HE_Scintillator_zminus.append(PlotOnSideGroup("RecHits_EtaPhi_HE_Scintillator_zminus_layer_{:02d}".format(i), Plot("EtaPhi_Minus_layer_{:02d}".format(i), title="RecHits_EtaPhi_HE_Scintillator_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_etaphi), ncols=1))
1650 
1651 _DigiHits_ADC_EE = PlotGroup("DigiHits_ADC_EE", [
1652  Plot("ADC_layer_{:02d}".format(i), title="DigiHits_ADC_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1653  ], ncols=4)
1654 
1655 _DigiHits_ADC_HE_Silicon = PlotGroup("DigiHits_ADC_HE_Silicon", [
1656  Plot("ADC_layer_{:02d}".format(i), title="DigiHits_ADC_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1657  ], ncols=4)
1658 
1659 _DigiHits_ADC_HE_Scintillator = PlotGroup("DigiHits_ADC_HE_Scintillator", [
1660  Plot("ADC_layer_{:02d}".format(i), title="DigiHits_ADC_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1661  ], ncols=4)
1662 
1663 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
1664 
1665 _DigiHits_Occupancy_EE_zplus = PlotGroup("DigiHits_Occupancy_EE_zplus", [
1666  Plot("DigiOccupancy_Plus_layer_{:02d}".format(i), title="DigiHits_Occupancy_EE_zplus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1667  ], ncols=4)
1668 
1669 _DigiHits_Occupancy_HE_Silicon_zplus = PlotGroup("DigiHits_Occupancy_HE_Silicon_zplus", [
1670  Plot("DigiOccupancy_Plus_layer_{:02d}".format(i), title="DigiHits_Occupancy_HE_Silicon_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1671  ], ncols=4)
1672 
1673 _DigiHits_Occupancy_HE_Scintillator_zplus = PlotGroup("DigiHits_Occupancy_HE_Scintillator_zplus", [
1674  Plot("DigiOccupancy_Plus_layer_{:02d}".format(i), title="DigiHits_Occupancy_HE_Scintillator_zplus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1675  ], ncols=4)
1676 
1677 _DigiHits_Occupancy_EE_zminus = PlotGroup("DigiHits_Occupancy_EE_zminus", [
1678  Plot("DigiOccupancy_Minus_layer_{:02d}".format(i), title="DigiHits_Occupancy_EE_zminus", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1679  ], ncols=4)
1680 
1681 _DigiHits_Occupancy_HE_Silicon_zminus = PlotGroup("DigiHits_Occupancy_HE_Silicon_zminus", [
1682  Plot("DigiOccupancy_Minus_layer_{:02d}".format(i), title="DigiHits_Occupancy_HE_Silicon_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1683  ], ncols=4)
1684 
1685 _DigiHits_Occupancy_HE_Scintillator_zminus = PlotGroup("DigiHits_Occupancy_HE_Scintillator_zminus", [
1686  Plot("DigiOccupancy_Minus_layer_{:02d}".format(i), title="DigiHits_Occupancy_HE_Scintillator_zminus", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1687  ], ncols=4)
1688 
1689 _common_XY = dict(removeEmptyBins=True, xbinlabelsize=10, xbinlabeloption="d", ymin=None)
1690 
1691 _DigiHits_Occupancy_XY_EE_zplus=[]
1692 for i in range(EE_min,EE_max+1):
1693  _DigiHits_Occupancy_XY_EE_zplus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_EE_zplus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_EE_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1694 
1695 _DigiHits_Occupancy_XY_HE_Silicon_zplus=[]
1696 for i in range(HESilicon_min,HESilicon_max+1):
1697  _DigiHits_Occupancy_XY_HE_Silicon_zplus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_HE_Silicon_zplus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_HE_Silicon_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1698 
1699 _DigiHits_Occupancy_XY_HE_Scintillator_zplus=[]
1700 for i in range(HEScintillator_min,HEScintillator_max+1):
1701  _DigiHits_Occupancy_XY_HE_Scintillator_zplus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_HE_Scintillator_zplus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_HE_Scintillator_zplus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1702 
1703 _DigiHits_Occupancy_XY_EE_zminus=[]
1704 for i in range(EE_min,EE_max+1):
1705  _DigiHits_Occupancy_XY_EE_zminus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_EE_zminus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_EE_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1706 
1707 _DigiHits_Occupancy_XY_HE_Silicon_zminus=[]
1708 for i in range(HESilicon_min,HESilicon_max+1):
1709  _DigiHits_Occupancy_XY_HE_Silicon_zminus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_HE_Silicon_zminus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_HE_Silicon_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1710 
1711 _DigiHits_Occupancy_XY_HE_Scintillator_zminus=[]
1712 for i in range(HEScintillator_min,HEScintillator_max+1):
1713  _DigiHits_Occupancy_XY_HE_Scintillator_zminus.append(PlotOnSideGroup("DigiHits_Occupancy_XY_HE_Scintillator_zminus_layer_{:02d}".format(i), Plot("DigiOccupancy_XY_layer_{:02d}".format(i), title="DigiHits_Occupancy_XY_HE_Scintillator_zminus", xtitle="Layer {}".format(i), drawStyle="COLZ", adjustMarginRight=0.1, **_common_XY), ncols=1))
1714 
1715 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ymin": 0.1, "ylog": True}
1716 
1717 _DigiHits_TOA_EE = PlotGroup("DigiHits_TOA_EE", [
1718  Plot("TOA_layer_{:02d}".format(i), title="DigiHits_TOA_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1719  ], ncols=4)
1720 
1721 _DigiHits_TOA_HE_Silicon = PlotGroup("DigiHits_TOA_HE_Silicon", [
1722  Plot("TOA_layer_{:02d}".format(i), title="DigiHits_TOA_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1723  ], ncols=4)
1724 
1725 _DigiHits_TOA_HE_Scintillator = PlotGroup("DigiHits_TOA_HE_Scintillator", [
1726  Plot("TOA_layer_{:02d}".format(i), title="DigiHits_TOA_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1727  ], ncols=4)
1728 
1729 _DigiHits_TOT_EE = PlotGroup("DigiHits_TOT_EE", [
1730  Plot("TOT_layer_{:02d}".format(i), title="DigiHits_TOT_EE", xtitle="Layer {}".format(i), **_common) for i in range(EE_min,EE_max+1)
1731  ], ncols=4)
1732 
1733 _DigiHits_TOT_HE_Silicon = PlotGroup("DigiHits_TOT_HE_Silicon", [
1734  Plot("TOT_layer_{:02d}".format(i), title="DigiHits_TOT_HE_Silicon", xtitle="Layer {}".format(i), **_common) for i in range(HESilicon_min,HESilicon_max+1)
1735  ], ncols=4)
1736 
1737 _DigiHits_TOT_HE_Scintillator = PlotGroup("DigiHits_TOT_HE_Scintillator", [
1738  Plot("TOT_layer_{:02d}".format(i), title="DigiHits_TOT_HE_Scintillator", xtitle="Layer {}".format(i), **_common) for i in range(HEScintillator_min,HEScintillator_max+1)
1739  ], ncols=4)
1740 
1741 #===================================================================================================================
1742 #Plot definition for HitCalibration
1743 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65, "ymin": 0.1, "ylog": False}
1744 
1745 _LayerOccupancy = PlotGroup("LayerOccupancy", [
1746  Plot("LayerOccupancy", title="LayerOccupancy", **_common)], ncols=1)
1747 
1748 _ReconstructableEnergyOverCPenergy = PlotGroup("ReconstructableEnergyOverCPenergy", [
1749  Plot("h_EoP_CPene_100_calib_fraction", title="EoP_CPene_100_calib_fraction", **_common),
1750  Plot("h_EoP_CPene_200_calib_fraction", title="EoP_CPene_200_calib_fraction", **_common),
1751  Plot("h_EoP_CPene_300_calib_fraction", title="EoP_CPene_300_calib_fraction", **_common),
1752 ])
1753 
1754 _ParticleFlowClusterHGCalFromMultiCl_Closest_EoverCPenergy = PlotGroup("ParticleFlowClusterHGCalFromMultiCl", [
1755  Plot("hgcal_EoP_CPene_100_calib_fraction", title="hgcal_EoP_CPene_100_calib_fraction", **_common),
1756  Plot("hgcal_EoP_CPene_200_calib_fraction", title="hgcal_EoP_CPene_200_calib_fraction", **_common),
1757  Plot("hgcal_EoP_CPene_300_calib_fraction", title="hgcal_EoP_CPene_300_calib_fraction", **_common),
1758 ])
1759 
1760 _EcalDrivenGsfElectronsFromMultiCl_Closest_EoverCPenergy = PlotGroup("EcalDrivenGsfElectronsFromMultiCl", [
1761  Plot("hgcal_ele_EoP_CPene_100_calib_fraction", title="hgcal_ele_EoP_CPene_100_calib_fraction", **_common),
1762  Plot("hgcal_ele_EoP_CPene_200_calib_fraction", title="hgcal_ele_EoP_CPene_200_calib_fraction", **_common),
1763  Plot("hgcal_ele_EoP_CPene_300_calib_fraction", title="hgcal_ele_EoP_CPene_300_calib_fraction", **_common),
1764 ])
1765 
1766 _PhotonsFromMultiCl_Closest_EoverCPenergy = PlotGroup("PhotonsFromMultiCl", [
1767  Plot("hgcal_photon_EoP_CPene_100_calib_fraction", title="hgcal_photon_EoP_CPene_100_calib_fraction", **_common),
1768  Plot("hgcal_photon_EoP_CPene_200_calib_fraction", title="hgcal_photon_EoP_CPene_200_calib_fraction", **_common),
1769  Plot("hgcal_photon_EoP_CPene_300_calib_fraction", title="hgcal_photon_EoP_CPene_300_calib_fraction", **_common),
1770 ])
1771 
1772 #=================================================================================================
1773 hgcalLayerClustersPlotter = Plotter()
1774 #We follow Chris categories in folders
1775 # [A] calculated "energy density" for cells in a) 120um, b) 200um, c) 300um, d) scint
1776 # (one entry per rechit, in the appropriate histo)
1777 hgcalLayerClustersPlotter.append("CellsEnergyDensityPerThickness", [
1778  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1779  ], PlotFolder(
1780  _cellsenedens_thick,
1781  loopSubFolders=False,
1782  purpose=PlotPurpose.Timing, page="CellsEnergyDensityPerThickness"
1783  ))
1784 
1785 # [B] number of layer clusters per event in a) 120um, b) 200um, c) 300um, d) scint
1786 # (one entry per event in each of the four histos)
1787 hgcalLayerClustersPlotter.append("TotalNumberofLayerClustersPerThickness", [
1788  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1789  ], PlotFolder(
1790  _totclusternum_thick,
1791  loopSubFolders=False,
1792  purpose=PlotPurpose.Timing, page="TotalNumberofLayerClustersPerThickness"
1793  ))
1794 
1795 # [C] number of layer clusters per layer (one entry per event in each histo)
1796 # z-
1797 hgcalLayerClustersPlotter.append("NumberofLayerClustersPerLayer_zminus", [
1798  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1799  ], PlotFolder(
1800  _totclusternum_layer_EE_zminus,
1801  _totclusternum_layer_FH_zminus,
1802  _totclusternum_layer_BH_zminus,
1803  loopSubFolders=False,
1804  purpose=PlotPurpose.Timing, page="NumberofLayerClustersPerLayer_zminus"
1805  ))
1806 
1807 # z+
1808 hgcalLayerClustersPlotter.append("NumberofLayerClustersPerLayer_zplus", [
1809  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1810  ], PlotFolder(
1811  _totclusternum_layer_EE_zplus,
1812  _totclusternum_layer_FH_zplus,
1813  _totclusternum_layer_BH_zplus,
1814  loopSubFolders=False,
1815  purpose=PlotPurpose.Timing, page="NumberofLayerClustersPerLayer_zplus"
1816  ))
1817 
1818 # [D] For each layer cluster:
1819 # number of cells in layer cluster, by layer - separate histos in each layer for 120um Si, 200/300um Si, Scint
1820 # NB: not all combinations exist; e.g. no 120um Si in layers with scint.
1821 # (One entry in the appropriate histo per layer cluster).
1822 # z-
1823 hgcalLayerClustersPlotter.append("CellsNumberPerLayerPerThickness_zminus", [
1824  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1825  ], PlotFolder(
1826  _cellsnum_perthick_perlayer_120_EE_zminus,
1827  _cellsnum_perthick_perlayer_120_FH_zminus,
1828  _cellsnum_perthick_perlayer_120_BH_zminus,
1829  _cellsnum_perthick_perlayer_200_EE_zminus,
1830  _cellsnum_perthick_perlayer_200_FH_zminus,
1831  _cellsnum_perthick_perlayer_200_BH_zminus,
1832  _cellsnum_perthick_perlayer_300_EE_zminus,
1833  _cellsnum_perthick_perlayer_300_FH_zminus,
1834  _cellsnum_perthick_perlayer_300_BH_zminus,
1835  _cellsnum_perthick_perlayer_scint_EE_zminus,
1836  _cellsnum_perthick_perlayer_scint_FH_zminus,
1837  _cellsnum_perthick_perlayer_scint_BH_zminus,
1838  loopSubFolders=False,
1839  purpose=PlotPurpose.Timing, page="CellsNumberPerLayerPerThickness_zminus"
1840  ))
1841 
1842 # z+
1843 hgcalLayerClustersPlotter.append("CellsNumberPerLayerPerThickness_zplus", [
1844  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1845  ], PlotFolder(
1846  _cellsnum_perthick_perlayer_120_EE_zplus,
1847  _cellsnum_perthick_perlayer_120_FH_zplus,
1848  _cellsnum_perthick_perlayer_120_BH_zplus,
1849  _cellsnum_perthick_perlayer_200_EE_zplus,
1850  _cellsnum_perthick_perlayer_200_FH_zplus,
1851  _cellsnum_perthick_perlayer_200_BH_zplus,
1852  _cellsnum_perthick_perlayer_300_EE_zplus,
1853  _cellsnum_perthick_perlayer_300_FH_zplus,
1854  _cellsnum_perthick_perlayer_300_BH_zplus,
1855  _cellsnum_perthick_perlayer_scint_EE_zplus,
1856  _cellsnum_perthick_perlayer_scint_FH_zplus,
1857  _cellsnum_perthick_perlayer_scint_BH_zplus,
1858  loopSubFolders=False,
1859  purpose=PlotPurpose.Timing, page="CellsNumberPerLayerPerThickness_zplus"
1860  ))
1861 
1862 # [E] For each layer cluster:
1863 # distance of cells from a) seed cell, b) max cell; and c), d): same with entries weighted by cell energy
1864 # separate histos in each layer for 120um Si, 200/300um Si, Scint
1865 # NB: not all combinations exist; e.g. no 120um Si in layers with scint.
1866 # (One entry in each of the four appropriate histos per cell in a layer cluster)
1867 # z-
1868 hgcalLayerClustersPlotter.append("CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zminus", [
1869  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1870  ], PlotFolder(
1871  _distancetomaxcell_perthickperlayer_120_EE_zminus,
1872  _distancetomaxcell_perthickperlayer_120_FH_zminus,
1873  _distancetomaxcell_perthickperlayer_120_BH_zminus,
1874  _distancetomaxcell_perthickperlayer_200_EE_zminus,
1875  _distancetomaxcell_perthickperlayer_200_FH_zminus,
1876  _distancetomaxcell_perthickperlayer_200_BH_zminus,
1877  _distancetomaxcell_perthickperlayer_300_EE_zminus,
1878  _distancetomaxcell_perthickperlayer_300_FH_zminus,
1879  _distancetomaxcell_perthickperlayer_300_BH_zminus,
1880  _distancetomaxcell_perthickperlayer_scint_EE_zminus,
1881  _distancetomaxcell_perthickperlayer_scint_FH_zminus,
1882  _distancetomaxcell_perthickperlayer_scint_BH_zminus,
1883  _distancetoseedcell_perthickperlayer_120_EE_zminus,
1884  _distancetoseedcell_perthickperlayer_120_FH_zminus,
1885  _distancetoseedcell_perthickperlayer_120_BH_zminus,
1886  _distancetoseedcell_perthickperlayer_200_EE_zminus,
1887  _distancetoseedcell_perthickperlayer_200_FH_zminus,
1888  _distancetoseedcell_perthickperlayer_200_BH_zminus,
1889  _distancetoseedcell_perthickperlayer_300_EE_zminus,
1890  _distancetoseedcell_perthickperlayer_300_FH_zminus,
1891  _distancetoseedcell_perthickperlayer_300_BH_zminus,
1892  _distancetoseedcell_perthickperlayer_scint_EE_zminus,
1893  _distancetoseedcell_perthickperlayer_scint_FH_zminus,
1894  _distancetoseedcell_perthickperlayer_scint_BH_zminus,
1895  _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus,
1896  _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus,
1897  _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus,
1898  _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus,
1899  _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus,
1900  _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus,
1901  _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus,
1902  _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus,
1903  _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus,
1904  _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zminus,
1905  _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zminus,
1906  _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zminus,
1907  _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus,
1908  _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus,
1909  _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus,
1910  _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus,
1911  _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus,
1912  _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus,
1913  _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus,
1914  _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus,
1915  _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus,
1916  _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zminus,
1917  _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zminus,
1918  _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zminus,
1919  _distancebetseedandmaxcell_perthickperlayer_120_EE_zminus,
1920  _distancebetseedandmaxcell_perthickperlayer_120_FH_zminus,
1921  _distancebetseedandmaxcell_perthickperlayer_120_BH_zminus,
1922  _distancebetseedandmaxcell_perthickperlayer_200_EE_zminus,
1923  _distancebetseedandmaxcell_perthickperlayer_200_FH_zminus,
1924  _distancebetseedandmaxcell_perthickperlayer_200_BH_zminus,
1925  _distancebetseedandmaxcell_perthickperlayer_300_EE_zminus,
1926  _distancebetseedandmaxcell_perthickperlayer_300_FH_zminus,
1927  _distancebetseedandmaxcell_perthickperlayer_300_BH_zminus,
1928  _distancebetseedandmaxcell_perthickperlayer_scint_EE_zminus,
1929  _distancebetseedandmaxcell_perthickperlayer_scint_FH_zminus,
1930  _distancebetseedandmaxcell_perthickperlayer_scint_BH_zminus,
1931  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus,
1932  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus,
1933  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus,
1934  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus,
1935  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus,
1936  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus,
1937  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus,
1938  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus,
1939  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus,
1940  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zminus,
1941  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zminus,
1942  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zminus,
1943  loopSubFolders=False,
1944  purpose=PlotPurpose.Timing, page="CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zminus"
1945  ))
1946 
1947 # z+
1948 hgcalLayerClustersPlotter.append("CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zplus", [
1949  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1950  ], PlotFolder(
1951  _distancetomaxcell_perthickperlayer_120_EE_zplus,
1952  _distancetomaxcell_perthickperlayer_120_FH_zplus,
1953  _distancetomaxcell_perthickperlayer_120_BH_zplus,
1954  _distancetomaxcell_perthickperlayer_200_EE_zplus,
1955  _distancetomaxcell_perthickperlayer_200_FH_zplus,
1956  _distancetomaxcell_perthickperlayer_200_BH_zplus,
1957  _distancetomaxcell_perthickperlayer_300_EE_zplus,
1958  _distancetomaxcell_perthickperlayer_300_FH_zplus,
1959  _distancetomaxcell_perthickperlayer_300_BH_zplus,
1960  _distancetomaxcell_perthickperlayer_scint_EE_zplus,
1961  _distancetomaxcell_perthickperlayer_scint_FH_zplus,
1962  _distancetomaxcell_perthickperlayer_scint_BH_zplus,
1963  _distancetoseedcell_perthickperlayer_120_EE_zplus,
1964  _distancetoseedcell_perthickperlayer_120_FH_zplus,
1965  _distancetoseedcell_perthickperlayer_120_BH_zplus,
1966  _distancetoseedcell_perthickperlayer_200_EE_zplus,
1967  _distancetoseedcell_perthickperlayer_200_FH_zplus,
1968  _distancetoseedcell_perthickperlayer_200_BH_zplus,
1969  _distancetoseedcell_perthickperlayer_300_EE_zplus,
1970  _distancetoseedcell_perthickperlayer_300_FH_zplus,
1971  _distancetoseedcell_perthickperlayer_300_BH_zplus,
1972  _distancetoseedcell_perthickperlayer_scint_EE_zplus,
1973  _distancetoseedcell_perthickperlayer_scint_FH_zplus,
1974  _distancetoseedcell_perthickperlayer_scint_BH_zplus,
1975  _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus,
1976  _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus,
1977  _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus,
1978  _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus,
1979  _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus,
1980  _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus,
1981  _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus,
1982  _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus,
1983  _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus,
1984  _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zplus,
1985  _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zplus,
1986  _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zplus,
1987  _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus,
1988  _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus,
1989  _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus,
1990  _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus,
1991  _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus,
1992  _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus,
1993  _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus,
1994  _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus,
1995  _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus,
1996  _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zplus,
1997  _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zplus,
1998  _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zplus,
1999  _distancebetseedandmaxcell_perthickperlayer_120_EE_zplus,
2000  _distancebetseedandmaxcell_perthickperlayer_120_FH_zplus,
2001  _distancebetseedandmaxcell_perthickperlayer_120_BH_zplus,
2002  _distancebetseedandmaxcell_perthickperlayer_200_EE_zplus,
2003  _distancebetseedandmaxcell_perthickperlayer_200_FH_zplus,
2004  _distancebetseedandmaxcell_perthickperlayer_200_BH_zplus,
2005  _distancebetseedandmaxcell_perthickperlayer_300_EE_zplus,
2006  _distancebetseedandmaxcell_perthickperlayer_300_FH_zplus,
2007  _distancebetseedandmaxcell_perthickperlayer_300_BH_zplus,
2008  _distancebetseedandmaxcell_perthickperlayer_scint_EE_zplus,
2009  _distancebetseedandmaxcell_perthickperlayer_scint_FH_zplus,
2010  _distancebetseedandmaxcell_perthickperlayer_scint_BH_zplus,
2011  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus,
2012  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus,
2013  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus,
2014  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus,
2015  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus,
2016  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus,
2017  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus,
2018  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus,
2019  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus,
2020  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zplus,
2021  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zplus,
2022  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zplus,
2023  loopSubFolders=False,
2024  purpose=PlotPurpose.Timing, page="CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zplus"
2025  ))
2026 
2027 # [F] Looking at the fraction of true energy that has been clustered; by layer and overall
2028 # z-
2029 hgcalLayerClustersPlotter.append("EnergyClusteredByLayerAndOverall_zminus", [
2030  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2031  ], PlotFolder(
2032  _energyclustered_perlayer_EE_zminus,
2033  _energyclustered_perlayer_FH_zminus,
2034  _energyclustered_perlayer_BH_zminus,
2035  loopSubFolders=False,
2036  purpose=PlotPurpose.Timing, page="EnergyClusteredByLayerAndOverall_zminus"
2037  ))
2038 # z+
2039 hgcalLayerClustersPlotter.append("EnergyClusteredByLayerAndOverall_zplus", [
2040  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2041  ], PlotFolder(
2042  _energyclustered_perlayer_EE_zplus,
2043  _energyclustered_perlayer_FH_zplus,
2044  _energyclustered_perlayer_BH_zplus,
2045  loopSubFolders=False,
2046  purpose=PlotPurpose.Timing, page="EnergyClusteredByLayerAndOverall_zplus"
2047  ))
2048 
2049 # [G] Miscellaneous plots:
2050 # longdepthbarycentre: The longitudinal depth barycentre. One entry per event.
2051 # mixedhitscluster: Number of clusters per event with hits in different thicknesses.
2052 # num_reco_cluster_eta: Number of reco clusters vs eta
2053 
2054 hgcalLayerClustersPlotter.append("Miscellaneous", [
2055  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2056  ], PlotFolder(
2057  _num_reco_cluster_eta,
2058  _energyclustered,
2059  _mixedhitsclusters,
2060  _longdepthbarycentre,
2061  loopSubFolders=False,
2062  purpose=PlotPurpose.Timing, page="Miscellaneous"
2063  ))
2064 
2065 # [H] SelectedCaloParticles plots
2066 hgcalLayerClustersPlotter.append("SelectedCaloParticles_Photons", [
2067  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/SelectedCaloParticles/22",
2068  ], PlotFolder(
2069  _SelectedCaloParticles,
2070  loopSubFolders=False,
2071  purpose=PlotPurpose.Timing, page="SelectedCaloParticles_Photons"
2072  ))
2073 
2074 # [I] Score of CaloParticles wrt Layer Clusters
2075 # z-
2076 hgcalLayerClustersPlotter.append("ScoreCaloParticlesToLayerClusters_zminus", [
2077  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2078  ], PlotFolder(
2079  _score_caloparticle_to_layerclusters_zminus,
2080  loopSubFolders=False,
2081  purpose=PlotPurpose.Timing, page="ScoreCaloParticlesToLayerClusters_zminus"))
2082 
2083 # z+
2084 hgcalLayerClustersPlotter.append("ScoreCaloParticlesToLayerClusters_zplus", [
2085  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2086  ], PlotFolder(
2087  _score_caloparticle_to_layerclusters_zplus,
2088  loopSubFolders=False,
2089  purpose=PlotPurpose.Timing, page="ScoreCaloParticlesToLayerClusters_zplus"))
2090 
2091 # [J] Score of LayerClusters wrt CaloParticles
2092 # z-
2093 hgcalLayerClustersPlotter.append("ScoreLayerClustersToCaloParticles_zminus", [
2094  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2095  ], PlotFolder(
2096  _score_layercluster_to_caloparticles_zminus,
2097  loopSubFolders=False,
2098  purpose=PlotPurpose.Timing, page="ScoreLayerClustersToCaloParticles_zminus"))
2099 
2100 # z+
2101 hgcalLayerClustersPlotter.append("ScoreLayerClustersToCaloParticles_zplus", [
2102  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2103  ], PlotFolder(
2104  _score_layercluster_to_caloparticles_zplus,
2105  loopSubFolders=False,
2106  purpose=PlotPurpose.Timing, page="ScoreLayerClustersToCaloParticles_zplus"))
2107 
2108 # [K] Shared Energy between CaloParticle and LayerClusters
2109 # z-
2110 hgcalLayerClustersPlotter.append("SharedEnergyC2L_zminus", [
2111  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2112  ], PlotFolder(
2113  _sharedEnergy_caloparticle_to_layercluster_zminus,
2114  loopSubFolders=False,
2115  purpose=PlotPurpose.Timing, page="SharedEnergyCaloParticleToLayerCluster_zminus"))
2116 
2117 # z+
2118 hgcalLayerClustersPlotter.append("SharedEnergyC2L_zplus", [
2119  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2120  ], PlotFolder(
2121  _sharedEnergy_caloparticle_to_layercluster_zplus,
2122  loopSubFolders=False,
2123  purpose=PlotPurpose.Timing, page="SharedEnergyCaloParticleToLayerCluster_zplus"))
2124 
2125 # [K2] Shared Energy between LayerClusters and CaloParticle
2126 # z-
2127 hgcalLayerClustersPlotter.append("SharedEnergyL2C_zminus", [
2128  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2129  ], PlotFolder(
2130  _sharedEnergy_layercluster_to_caloparticle_zminus,
2131  loopSubFolders=False,
2132  purpose=PlotPurpose.Timing, page="SharedEnergyLayerClusterToCaloParticle_zminus"))
2133 
2134 # z+
2135 hgcalLayerClustersPlotter.append("SharedEnergyL2C_zplus", [
2136  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2137  ], PlotFolder(
2138  _sharedEnergy_layercluster_to_caloparticle_zplus,
2139  loopSubFolders=False,
2140  purpose=PlotPurpose.Timing, page="SharedEnergyLayerClusterToCaloParticle_zplus"))
2141 
2142 # [L] Cell Association per Layer
2143 # z-
2144 hgcalLayerClustersPlotter.append("CellAssociation_zminus", [
2145  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2146  ], PlotFolder(
2147  _cell_association_table_zminus,
2148  loopSubFolders=False,
2149  purpose=PlotPurpose.Timing, page="CellAssociation_zminus"))
2150 
2151 # z+
2152 hgcalLayerClustersPlotter.append("CellAssociation_zplus", [
2153  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2154  ], PlotFolder(
2155  _cell_association_table_zplus,
2156  loopSubFolders=False,
2157  purpose=PlotPurpose.Timing, page="CellAssociation_zplus"))
2158 
2159 # [M] Efficiency Plots
2160 # z-
2161 hgcalLayerClustersPlotter.append("Efficiencies_zminus", [
2162  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2163  ], PlotFolder(
2164  _efficiencies_zminus,
2165  loopSubFolders=False,
2166  purpose=PlotPurpose.Timing, page="Efficiencies_zminus"))
2167 
2168 # z+
2169 hgcalLayerClustersPlotter.append("Efficiencies_zplus", [
2170  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2171  ], PlotFolder(
2172  _efficiencies_zplus,
2173  loopSubFolders=False,
2174  purpose=PlotPurpose.Timing, page="Efficiencies_zplus"))
2175 
2176 # [L] Duplicate Plots
2177 # z-
2178 hgcalLayerClustersPlotter.append("Duplicates_zminus", [
2179  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2180  ], PlotFolder(
2181  _duplicates_zminus,
2182  loopSubFolders=False,
2183  purpose=PlotPurpose.Timing, page="Duplicates_zminus"))
2184 
2185 # z+
2186 hgcalLayerClustersPlotter.append("Duplicates_zplus", [
2187  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2188  ], PlotFolder(
2189  _duplicates_zplus,
2190  loopSubFolders=False,
2191  purpose=PlotPurpose.Timing, page="Duplicates_zplus"))
2192 
2193 # [M] Fake Rate Plots
2194 # z-
2195 hgcalLayerClustersPlotter.append("FakeRate_zminus", [
2196  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2197  ], PlotFolder(
2198  _fakes_zminus,
2199  loopSubFolders=False,
2200  purpose=PlotPurpose.Timing, page="Fakes_zminus"))
2201 
2202 # z+
2203 hgcalLayerClustersPlotter.append("FakeRate_zplus", [
2204  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2205  ], PlotFolder(
2206  _fakes_zplus,
2207  loopSubFolders=False,
2208  purpose=PlotPurpose.Timing, page="Fakes_zplus"))
2209 
2210 # [N] Merge Rate Plots
2211 # z-
2212 hgcalLayerClustersPlotter.append("MergeRate_zminus", [
2213  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2214  ], PlotFolder(
2215  _merges_zminus,
2216  loopSubFolders=False,
2217  purpose=PlotPurpose.Timing, page="Merges_zminus"))
2218 
2219 # z+
2220 hgcalLayerClustersPlotter.append("MergeRate_zplus", [
2221  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2222  ], PlotFolder(
2223  _merges_zplus,
2224  loopSubFolders=False,
2225  purpose=PlotPurpose.Timing, page="Merges_zplus"))
2226 
2227 # [O] Energy vs Score 2D plots CP to LC
2228 # z-
2229 for i,item in enumerate(_energyscore_cp2lc_zminus, start=1):
2230  hgcalLayerClustersPlotter.append("Energy_vs_Score_CP2LC_zminus", [
2231  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2232  ], PlotFolder(
2233  item,
2234  loopSubFolders=False,
2235  purpose=PlotPurpose.Timing, page="Energy_vs_Score_CP2LC_zminus"))
2236 
2237 # z+
2238 for i,item in enumerate(_energyscore_cp2lc_zplus, start=1):
2239  hgcalLayerClustersPlotter.append("Energy_vs_Score_CP2LC_zplus", [
2240  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2241  ], PlotFolder(
2242  item,
2243  loopSubFolders=False,
2244  purpose=PlotPurpose.Timing, page="Energy_vs_Score_CP2LC_zplus"))
2245 
2246 # [P] Energy vs Score 2D plots LC to CP
2247 # z-
2248 for i,item in enumerate(_energyscore_lc2cp_zminus, start=1):
2249  hgcalLayerClustersPlotter.append("Energy_vs_Score_LC2CP_zminus", [
2250  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2251  ], PlotFolder(
2252  item,
2253  loopSubFolders=False,
2254  purpose=PlotPurpose.Timing, page="Energy_vs_Score_LC2CP_zminus"))
2255 
2256 # z+
2257 for i,item in enumerate(_energyscore_lc2cp_zplus, start=1):
2258  hgcalLayerClustersPlotter.append("Energy_vs_Score_LC2CP_zplus", [
2259  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
2260  ], PlotFolder(
2261  item,
2262  loopSubFolders=False,
2263  purpose=PlotPurpose.Timing, page="Energy_vs_Score_LC2CP_zplus"))
2264 
2265 #=================================================================================================
2266 hgcalMultiClustersPlotter = Plotter()
2267 # [A] Score of CaloParticles wrt Multi Clusters
2268 hgcalMultiClustersPlotter.append("ScoreCaloParticlesToMultiClusters", [
2269  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2270  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2271  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2272  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2273  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2274  ], PlotFolder(
2275  _score_caloparticle_to_multiclusters,
2276  loopSubFolders=False,
2277  purpose=PlotPurpose.Timing, page="ScoreCaloParticlesToMultiClusters"))
2278 
2279 # [B] Score of MultiClusters wrt CaloParticles
2280 hgcalMultiClustersPlotter.append("ScoreMultiClustersToCaloParticles", [
2281  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2282  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2283  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2284  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2285  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2286  ], PlotFolder(
2287  _score_multicluster_to_caloparticles,
2288  loopSubFolders=False,
2289  purpose=PlotPurpose.Timing, page="ScoreMultiClustersToCaloParticles"))
2290 
2291 # [C] Shared Energy between CaloParticle and MultiClusters
2292 hgcalMultiClustersPlotter.append("SharedEnergy_CP2MCL", [
2293  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2294  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2295  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2296  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2297  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2298  ], PlotFolder(
2299  _sharedEnergy_caloparticle_to_multicluster,
2300  loopSubFolders=False,
2301  purpose=PlotPurpose.Timing, page="SharedEnergyCaloParticleToMultiCluster"))
2302 
2303 # [C2] Shared Energy between MultiClusters and CaloParticle
2304 hgcalMultiClustersPlotter.append("SharedEnergy_MCL2CP", [
2305  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2306  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2307  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2308  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2309  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2310  ], PlotFolder(
2311  _sharedEnergy_multicluster_to_caloparticle,
2312  loopSubFolders=False,
2313  purpose=PlotPurpose.Timing, page="SharedEnergyMultiClusterToCaloParticle"))
2314 
2315 # [E] Efficiency Plots
2316 hgcalMultiClustersPlotter.append("Efficiencies", [
2317  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2318  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2319  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2320  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2321  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2322  ], PlotFolder(
2323  _efficiencies,
2324  loopSubFolders=False,
2325  purpose=PlotPurpose.Timing, page="Efficiencies"))
2326 
2327 # [F] Duplicate Plots
2328 hgcalMultiClustersPlotter.append("Duplicates", [
2329  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2330  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2331  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2332  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2333  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2334  ], PlotFolder(
2335  _duplicates,
2336  loopSubFolders=False,
2337  purpose=PlotPurpose.Timing, page="Duplicates"))
2338 
2339 # [G] Fake Rate Plots
2340 hgcalMultiClustersPlotter.append("FakeRate", [
2341  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2342  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2343  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2344  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2345  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2346  ], PlotFolder(
2347  _fakes,
2348  loopSubFolders=False,
2349  purpose=PlotPurpose.Timing, page="Fakes"))
2350 
2351 # [H] Merge Rate Plots
2352 hgcalMultiClustersPlotter.append("MergeRate", [
2353  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2354  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2355  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2356  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2357  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2358  ], PlotFolder(
2359  _merges,
2360  loopSubFolders=False,
2361  purpose=PlotPurpose.Timing, page="Merges"))
2362 
2363 # [I] Energy vs Score 2D plots CP to MCL and MCL to CP
2364 hgcalMultiClustersPlotter.append("Energy_vs_Score_CP2MCL_MCL2CP", [
2365  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2366  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2367  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2368  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2369  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2370  ], PlotFolder(
2371  #_energyscore_cp2mcl_mcl2cp,
2372  _energyscore_cp2mcl,
2373  loopSubFolders=False,
2374  purpose=PlotPurpose.Timing, page="Energy_vs_Score_CP2MCL"))
2375 
2376 # [J] Energy vs Score 2D plots MCL to CP
2377 hgcalMultiClustersPlotter.append("Energy_vs_Score_MCL2CP", [
2378  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2379  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2380  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2381  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2382  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2383  ], PlotFolder(
2384  _energyscore_mcl2cp,
2385  loopSubFolders=False,
2386  purpose=PlotPurpose.Timing, page="Energy_vs_Score_MCL2CP"))
2387 
2388 #[K] Number of multiclusters per event.
2389 hgcalMultiClustersPlotter.append("NumberofMultiClusters", [
2390  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2391  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2392  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2393  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2394  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2395  ], PlotFolder(
2396  _totmulticlusternum,
2397  loopSubFolders=False,
2398  purpose=PlotPurpose.Timing, page="NumberofMultiClusters"
2399  ))
2400 
2401 #[L] total number of layer clusters in multicluster per event and per layer
2402 hgcalMultiClustersPlotter.append("NumberofLayerClustersinMultiClusterPerEventAndPerLayer", [
2403  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2404  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2405  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2406  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2407  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2408  ], PlotFolder(
2409  _clusternum_in_multicluster,
2410  _clusternum_in_multicluster_vs_layer,
2411  _clusternum_in_multicluster_perlayer_zminus_EE,
2412  _clusternum_in_multicluster_perlayer_zminus_FH,
2413  _clusternum_in_multicluster_perlayer_zminus_BH,
2414  _clusternum_in_multicluster_perlayer_zplus_EE,
2415  _clusternum_in_multicluster_perlayer_zplus_FH,
2416  _clusternum_in_multicluster_perlayer_zplus_BH,
2417  loopSubFolders=False,
2418  purpose=PlotPurpose.Timing, page="NumberofLayerClustersinMultiClusterPerEventAndPerLayer"
2419  ))
2420 
2421 #[M] For each multicluster: pt, eta, phi, energy, x, y, z.
2422 hgcalMultiClustersPlotter.append("MultiClustersPtEtaPhiEneXYZ", [
2423  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2424  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2425  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2426  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2427  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2428  ], PlotFolder(
2429  _multicluster_pt,
2430  _multicluster_eta,
2431  _multicluster_phi,
2432  _multicluster_energy,
2433  _multicluster_x,
2434  _multicluster_y,
2435  _multicluster_z,
2436  loopSubFolders=False,
2437  purpose=PlotPurpose.Timing, page="MultiClustersPtEtaPhiEneXYZ"
2438  ))
2439 
2440 #[N] Multicluster first, last, total number of layers
2441 hgcalMultiClustersPlotter.append("NumberofMultiClusters_First_Last_NLayers", [
2442  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2443  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2444  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2445  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2446  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2447  ], PlotFolder(
2448  _multicluster_firstlayer,
2449  _multicluster_lastlayer,
2450  _multicluster_layersnum,
2451  loopSubFolders=False,
2452  purpose=PlotPurpose.Timing, page="NumberofMultiClusters_First_Last_NLayers"
2453  ))
2454 
2455 #[O] Multiplicity of layer clusters in multicluster
2456 hgcalMultiClustersPlotter.append("Multiplicity", [
2457  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2458  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2459  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2460  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2461  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2462  ], PlotFolder(
2463  _multiplicityOfLCinMCL,
2464  _multiplicityOfLCinMCL_vs_layerclusterenergy,
2465  loopSubFolders=False,
2466  purpose=PlotPurpose.Timing, page="Multiplicity",
2467  numberOfEventsHistogram=_multiplicity_numberOfEventsHistogram
2468  ))
2469 
2470 #We append here two PlotFolder because we want the text to be in percent
2471 #and the number of events are different in zplus and zminus
2472 hgcalMultiClustersPlotter.append("Multiplicity", [
2473  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2474  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2475  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2476  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2477  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2478  ], PlotFolder(
2479  _multiplicityOfLCinMCL_vs_layercluster_zminus,
2480  loopSubFolders=False,
2481  purpose=PlotPurpose.Timing, page="Multiplicity",
2482  numberOfEventsHistogram=_multiplicity_zminus_numberOfEventsHistogram
2483  ))
2484 
2485 hgcalMultiClustersPlotter.append("Multiplicity", [
2486  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersMIP_MIPMultiClustersFromTracksterByCA",
2487  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersTrk_TrkMultiClustersFromTracksterByCA",
2488  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersEM_MultiClustersFromTracksterByCA",
2489  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/multiClustersFromTrackstersHAD_MultiClustersFromTracksterByCA",
2490  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalMultiClusters",
2491  ], PlotFolder(
2492  _multiplicityOfLCinMCL_vs_layercluster_zplus,
2493  loopSubFolders=False,
2494  purpose=PlotPurpose.Timing, page="Multiplicity",
2495  numberOfEventsHistogram=_multiplicity_zplus_numberOfEventsHistogram
2496  ))
2497 
2498 #=================================================================================================
2499 # hitValidation
2500 hgcalHitPlotter = Plotter()
2501 
2502 hgcalHitPlotter.append("SimHits_Validation", [
2503  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HitValidation",
2504  ], PlotFolder(
2505  _HitValidation,
2506  loopSubFolders=False,
2507  purpose=PlotPurpose.Timing, page="SimHits_Validation"
2508  ))
2509 
2510 hgcalHitPlotter.append("SimHits_Occupancy_zplus", [
2511  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2512  ], PlotFolder(
2513  _SimHits_Occupancy_EE_zplus,
2514  loopSubFolders=False,
2515  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zplus"
2516  ))
2517 
2518 hgcalHitPlotter.append("SimHits_Occupancy_zplus", [
2519  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2520  ], PlotFolder(
2521  _SimHits_Occupancy_HE_Silicon_zplus,
2522  loopSubFolders=False,
2523  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zplus"
2524  ))
2525 
2526 hgcalHitPlotter.append("SimHits_Occupancy_zplus", [
2527  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive"
2528  ], PlotFolder(
2529  _SimHits_Occupancy_HE_Scintillator_zplus,
2530  loopSubFolders=False,
2531  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zplus"
2532  ))
2533 
2534 hgcalHitPlotter.append("SimHits_Occupancy_zminus", [
2535  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2536  ], PlotFolder(
2537  _SimHits_Occupancy_EE_zminus,
2538  loopSubFolders=False,
2539  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zminus"
2540  ))
2541 
2542 hgcalHitPlotter.append("SimHits_Occupancy_zminus", [
2543  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2544  ], PlotFolder(
2545  _SimHits_Occupancy_HE_Silicon_zminus,
2546  loopSubFolders=False,
2547  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zminus"
2548  ))
2549 
2550 hgcalHitPlotter.append("SimHits_Occupancy_zminus", [
2551  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive"
2552  ], PlotFolder(
2553  _SimHits_Occupancy_HE_Scintillator_zminus,
2554  loopSubFolders=False,
2555  purpose=PlotPurpose.Timing, page="SimHits_Occupancy_zminus"
2556  ))
2557 
2558 for i,item in enumerate(_SimHits_EtaPhi_EE_zplus, start=1):
2559  hgcalHitPlotter.append("SimHits_EtaPhi_zplus", [
2560  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2561  ], PlotFolder(
2562  item,
2563  loopSubFolders=False,
2564  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zplus"
2565  ))
2566 
2567 for i,item in enumerate(_SimHits_EtaPhi_HE_Silicon_zplus, start=1):
2568  hgcalHitPlotter.append("SimHits_EtaPhi_zplus", [
2569  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2570  ], PlotFolder(
2571  item,
2572  loopSubFolders=False,
2573  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zplus"
2574  ))
2575 
2576 for i,item in enumerate(_SimHits_EtaPhi_HE_Scintillator_zplus, start=1):
2577  hgcalHitPlotter.append("SimHits_EtaPhi_zplus", [
2578  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive",
2579  ], PlotFolder(
2580  item,
2581  loopSubFolders=False,
2582  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zplus"))
2583 
2584 for i,item in enumerate(_SimHits_EtaPhi_EE_zminus, start=1):
2585  hgcalHitPlotter.append("SimHits_EtaPhi_zminus", [
2586  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2587  ], PlotFolder(
2588  item,
2589  loopSubFolders=False,
2590  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zminus"
2591  ))
2592 
2593 for i,item in enumerate(_SimHits_EtaPhi_HE_Silicon_zminus, start=1):
2594  hgcalHitPlotter.append("SimHits_EtaPhi_zminus", [
2595  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2596  ], PlotFolder(
2597  item,
2598  loopSubFolders=False,
2599  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zminus"
2600  ))
2601 
2602 for i,item in enumerate(_SimHits_EtaPhi_HE_Scintillator_zminus, start=1):
2603  hgcalHitPlotter.append("SimHits_EtaPhi_zminus", [
2604  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive"
2605  ], PlotFolder(
2606  item,
2607  loopSubFolders=False,
2608  purpose=PlotPurpose.Timing, page="SimHits_EtaPhi_zminus"
2609  ))
2610 
2611 hgcalHitPlotter.append("SimHits_Energy", [
2612  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2613  ], PlotFolder(
2614  _SimHits_Energy_EE_0,
2615  loopSubFolders=False,
2616  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2617  ))
2618 
2619 hgcalHitPlotter.append("SimHits_Energy", [
2620  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2621  ], PlotFolder(
2622  _SimHits_Energy_HE_Silicon_0,
2623  loopSubFolders=False,
2624  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2625  ))
2626 
2627 hgcalHitPlotter.append("SimHits_Energy", [
2628  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive"
2629  ], PlotFolder(
2630  _SimHits_Energy_HE_Scintillator_0,
2631  loopSubFolders=False,
2632  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2633  ))
2634 
2635 hgcalHitPlotter.append("SimHits_Energy", [
2636  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalEESensitive"
2637  ], PlotFolder(
2638  _SimHits_Energy_EE_1,
2639  loopSubFolders=False,
2640  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2641  ))
2642 
2643 hgcalHitPlotter.append("SimHits_Energy", [
2644  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHESiliconSensitive"
2645  ], PlotFolder(
2646  _SimHits_Energy_HE_Silicon_1,
2647  loopSubFolders=False,
2648  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2649  ))
2650 
2651 hgcalHitPlotter.append("SimHits_Energy", [
2652  "DQMData/Run 1/HGCAL/Run summary/HGCalSimHitsV/HGCalHEScintillatorSensitive"
2653  ], PlotFolder(
2654  _SimHits_Energy_HE_Scintillator_1,
2655  loopSubFolders=False,
2656  purpose=PlotPurpose.Timing, page="SimHits_Energy"
2657  ))
2658 
2659 hgcalHitPlotter.append("RecHits_Occupancy_zplus", [
2660  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalEESensitive"
2661  ], PlotFolder(
2662  _RecHits_Occupancy_EE_zplus,
2663  loopSubFolders=False,
2664  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zplus"
2665  ))
2666 
2667 hgcalHitPlotter.append("RecHits_Occupancy_zplus", [
2668  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHESiliconSensitive"
2669  ], PlotFolder(
2670  _RecHits_Occupancy_HE_Silicon_zplus,
2671  loopSubFolders=False,
2672  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zplus"
2673  ))
2674 
2675 hgcalHitPlotter.append("RecHits_Occupancy_zplus", [
2676  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHEScintillatorSensitive"
2677  ], PlotFolder(
2678  _RecHits_Occupancy_HE_Scintillator_zplus,
2679  loopSubFolders=False,
2680  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zplus"
2681  ))
2682 
2683 hgcalHitPlotter.append("RecHits_Occupancy_zminus", [
2684  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalEESensitive"
2685  ], PlotFolder(
2686  _RecHits_Occupancy_EE_zminus,
2687  loopSubFolders=False,
2688  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zminus"
2689  ))
2690 
2691 hgcalHitPlotter.append("RecHits_Occupancy_zminus", [
2692  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHESiliconSensitive"
2693  ], PlotFolder(
2694  _RecHits_Occupancy_HE_Silicon_zminus,
2695  loopSubFolders=False,
2696  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zminus"
2697  ))
2698 
2699 hgcalHitPlotter.append("RecHits_Occupancy_zminus", [
2700  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHEScintillatorSensitive"
2701  ], PlotFolder(
2702  _RecHits_Occupancy_HE_Scintillator_zminus,
2703  loopSubFolders=False,
2704  purpose=PlotPurpose.Timing, page="RecHits_Occupancy_zminus"
2705  ))
2706 
2707 hgcalHitPlotter.append("RecHits_Energy", [
2708  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalEESensitive"
2709  ], PlotFolder(
2710  _RecHits_Energy_EE,
2711  loopSubFolders=False,
2712  purpose=PlotPurpose.Timing, page="RecHits_Energy"
2713  ))
2714 
2715 hgcalHitPlotter.append("RecHits_Energy", [
2716  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHESiliconSensitive"
2717  ], PlotFolder(
2718  _RecHits_Energy_HE_Silicon,
2719  loopSubFolders=False,
2720  purpose=PlotPurpose.Timing, page="RecHits_Energy"
2721  ))
2722 
2723 hgcalHitPlotter.append("RecHits_Energy", [
2724  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHEScintillatorSensitive"
2725  ], PlotFolder(
2726  _RecHits_Energy_HE_Scintillator,
2727  loopSubFolders=False,
2728  purpose=PlotPurpose.Timing, page="RecHits_Energy"
2729  ))
2730 
2731 for i,item in enumerate(_RecHits_EtaPhi_EE_zplus, start=1):
2732  hgcalHitPlotter.append("RecHits_EtaPhi_zplus", [
2733  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalEESensitive"
2734  ], PlotFolder(
2735  item,
2736  loopSubFolders=False,
2737  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zplus"
2738  ))
2739 
2740 for i,item in enumerate(_RecHits_EtaPhi_HE_Silicon_zplus, start=1):
2741  hgcalHitPlotter.append("RecHits_EtaPhi_zplus", [
2742  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHESiliconSensitive"
2743  ], PlotFolder(
2744  item,
2745  loopSubFolders=False,
2746  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zplus"
2747  ))
2748 
2749 for i,item in enumerate(_RecHits_EtaPhi_HE_Scintillator_zplus, start=1):
2750  hgcalHitPlotter.append("RecHits_EtaPhi_zplus", [
2751  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHEScintillatorSensitive",
2752  ], PlotFolder(
2753  item,
2754  loopSubFolders=False,
2755  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zplus"))
2756 
2757 for i,item in enumerate(_RecHits_EtaPhi_EE_zminus, start=1):
2758  hgcalHitPlotter.append("RecHits_EtaPhi_zminus", [
2759  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalEESensitive"
2760  ], PlotFolder(
2761  item,
2762  loopSubFolders=False,
2763  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zminus"
2764  ))
2765 
2766 for i,item in enumerate(_RecHits_EtaPhi_HE_Silicon_zminus, start=1):
2767  hgcalHitPlotter.append("RecHits_EtaPhi_zminus", [
2768  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHESiliconSensitive"
2769  ], PlotFolder(
2770  item,
2771  loopSubFolders=False,
2772  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zminus"
2773  ))
2774 
2775 for i,item in enumerate(_RecHits_EtaPhi_HE_Scintillator_zminus, start=1):
2776  hgcalHitPlotter.append("RecHits_EtaPhi_zminus", [
2777  "DQMData/Run 1/HGCAL/Run summary/HGCalRecHitsV/HGCalHEScintillatorSensitive"
2778  ], PlotFolder(
2779  item,
2780  loopSubFolders=False,
2781  purpose=PlotPurpose.Timing, page="RecHits_EtaPhi_zminus"
2782  ))
2783 
2784 hgcalHitPlotter.append("DigiHits_ADC", [
2785  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2786  ], PlotFolder(
2787  _DigiHits_ADC_EE,
2788  loopSubFolders=False,
2789  purpose=PlotPurpose.Timing, page="DigiHits_ADC"
2790  ))
2791 
2792 hgcalHitPlotter.append("DigiHits_ADC", [
2793  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2794  ], PlotFolder(
2795  _DigiHits_ADC_HE_Silicon,
2796  loopSubFolders=False,
2797  purpose=PlotPurpose.Timing, page="DigiHits_ADC"
2798  ))
2799 
2800 hgcalHitPlotter.append("DigiHits_ADC", [
2801  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive"
2802  ], PlotFolder(
2803  _DigiHits_ADC_HE_Scintillator,
2804  loopSubFolders=False,
2805  purpose=PlotPurpose.Timing, page="DigiHits_ADC"
2806  ))
2807 
2808 hgcalHitPlotter.append("DigiHits_Occupancy_zplus", [
2809  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2810  ], PlotFolder(
2811  _DigiHits_Occupancy_EE_zplus,
2812  loopSubFolders=False,
2813  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zplus"
2814  ))
2815 
2816 hgcalHitPlotter.append("DigiHits_Occupancy_zplus", [
2817  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2818  ], PlotFolder(
2819  _DigiHits_Occupancy_HE_Silicon_zplus,
2820  loopSubFolders=False,
2821  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zplus"
2822  ))
2823 
2824 hgcalHitPlotter.append("DigiHits_Occupancy_zplus", [
2825  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive"
2826  ], PlotFolder(
2827  _DigiHits_Occupancy_HE_Scintillator_zplus,
2828  loopSubFolders=False,
2829  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zplus"
2830  ))
2831 
2832 hgcalHitPlotter.append("DigiHits_Occupancy_zminus", [
2833  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2834  ], PlotFolder(
2835  _DigiHits_Occupancy_EE_zminus,
2836  loopSubFolders=False,
2837  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zminus"
2838  ))
2839 
2840 hgcalHitPlotter.append("DigiHits_Occupancy_zminus", [
2841  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2842  ], PlotFolder(
2843  _DigiHits_Occupancy_HE_Silicon_zminus,
2844  loopSubFolders=False,
2845  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zminus"
2846  ))
2847 
2848 hgcalHitPlotter.append("DigiHits_Occupancy_zminus", [
2849  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive"
2850  ], PlotFolder(
2851  _DigiHits_Occupancy_HE_Scintillator_zminus,
2852  loopSubFolders=False,
2853  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_zminus"
2854  ))
2855 
2856 for i,item in enumerate(_DigiHits_Occupancy_XY_EE_zplus, start=1):
2857  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zplus", [
2858  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2859  ], PlotFolder(
2860  item,
2861  loopSubFolders=False,
2862  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zplus"
2863  ))
2864 
2865 for i,item in enumerate(_DigiHits_Occupancy_XY_HE_Silicon_zplus, start=1):
2866  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zplus", [
2867  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2868  ], PlotFolder(
2869  item,
2870  loopSubFolders=False,
2871  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zplus"
2872  ))
2873 
2874 for i,item in enumerate(_DigiHits_Occupancy_XY_HE_Scintillator_zplus, start=1):
2875  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zplus", [
2876  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive",
2877  ], PlotFolder(
2878  item,
2879  loopSubFolders=False,
2880  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zplus"))
2881 
2882 for i,item in enumerate(_DigiHits_Occupancy_XY_EE_zminus, start=1):
2883  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zminus", [
2884  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2885  ], PlotFolder(
2886  item,
2887  loopSubFolders=False,
2888  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zminus"
2889  ))
2890 
2891 for i,item in enumerate(_DigiHits_Occupancy_XY_HE_Silicon_zminus, start=1):
2892  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zminus", [
2893  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2894  ], PlotFolder(
2895  item,
2896  loopSubFolders=False,
2897  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zminus"
2898  ))
2899 
2900 for i,item in enumerate(_DigiHits_Occupancy_XY_HE_Scintillator_zminus, start=1):
2901  hgcalHitPlotter.append("DigiHits_Occupancy_XY_zminus", [
2902  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive",
2903  ], PlotFolder(
2904  item,
2905  loopSubFolders=False,
2906  purpose=PlotPurpose.Timing, page="DigiHits_Occupancy_XY_zminus"))
2907 
2908 hgcalHitPlotter.append("DigiHits_TOA", [
2909  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2910  ], PlotFolder(
2911  _DigiHits_TOA_EE,
2912  loopSubFolders=False,
2913  purpose=PlotPurpose.Timing, page="DigiHits_TOA"
2914  ))
2915 
2916 hgcalHitPlotter.append("DigiHits_TOA", [
2917  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2918  ], PlotFolder(
2919  _DigiHits_TOA_HE_Silicon,
2920  loopSubFolders=False,
2921  purpose=PlotPurpose.Timing, page="DigiHits_TOA"
2922  ))
2923 
2924 hgcalHitPlotter.append("DigiHits_TOA", [
2925  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive"
2926  ], PlotFolder(
2927  _DigiHits_TOA_HE_Scintillator,
2928  loopSubFolders=False,
2929  purpose=PlotPurpose.Timing, page="DigiHits_TOA"
2930  ))
2931 
2932 hgcalHitPlotter.append("DigiHits_TOT", [
2933  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalEESensitive"
2934  ], PlotFolder(
2935  _DigiHits_TOT_EE,
2936  loopSubFolders=False,
2937  purpose=PlotPurpose.Timing, page="DigiHits_TOT"
2938  ))
2939 
2940 hgcalHitPlotter.append("DigiHits_TOT", [
2941  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHESiliconSensitive"
2942  ], PlotFolder(
2943  _DigiHits_TOT_HE_Silicon,
2944  loopSubFolders=False,
2945  purpose=PlotPurpose.Timing, page="DigiHits_TOT"
2946  ))
2947 
2948 hgcalHitPlotter.append("DigiHits_TOT", [
2949  "DQMData/Run 1/HGCAL/Run summary/HGCalDigisV/HGCalHEScintillatorSensitive"
2950  ], PlotFolder(
2951  _DigiHits_TOT_HE_Scintillator,
2952  loopSubFolders=False,
2953  purpose=PlotPurpose.Timing, page="DigiHits_TOT"
2954  ))
2955 #=================================================================================================
2956 # hitCalibration
2957 hgcalHitCalibPlotter = Plotter()
2958 
2959 hgcalHitCalibPlotter.append("Layer_Occupancy", [
2960  "DQMData/Run 1/HGCalHitCalibration/Run summary",
2961  ], PlotFolder(
2962  _LayerOccupancy,
2963  loopSubFolders=False,
2964  purpose=PlotPurpose.Timing, page="Layer_Occupancy"
2965  ))
2966 hgcalHitCalibPlotter.append("ReconstructableEnergyOverCPenergy", [
2967  "DQMData/Run 1/HGCalHitCalibration/Run summary",
2968  ], PlotFolder(
2969  _ReconstructableEnergyOverCPenergy,
2970  loopSubFolders=False,
2971  purpose=PlotPurpose.Timing, page="ReconstructableEnergyOverCPenergy"
2972  ))
2973 
2974 hgcalHitCalibPlotter.append("ParticleFlowClusterHGCalFromMultiCl_Closest_EoverCPenergy", [
2975  "DQMData/Run 1/HGCalHitCalibration/Run summary",
2976  ], PlotFolder(
2977  _ParticleFlowClusterHGCalFromMultiCl_Closest_EoverCPenergy,
2978  loopSubFolders=False,
2979  purpose=PlotPurpose.Timing, page="ParticleFlowClusterHGCalFromMultiCl_Closest_EoverCPenergy"
2980  ))
2981 
2982 hgcalHitCalibPlotter.append("PhotonsFromMultiCl_Closest_EoverCPenergy", [
2983  "DQMData/Run 1/HGCalHitCalibration/Run summary",
2984  ], PlotFolder(
2985  _PhotonsFromMultiCl_Closest_EoverCPenergy,
2986  loopSubFolders=False,
2987  purpose=PlotPurpose.Timing, page="PhotonsFromMultiCl_Closest_EoverCPenergy"
2988  ))
2989 
2990 hgcalHitCalibPlotter.append("EcalDrivenGsfElectronsFromMultiCl_Closest_EoverCPenergy", [
2991  "DQMData/Run 1/HGCalHitCalibration/Run summary",
2992  ], PlotFolder(
2993  _EcalDrivenGsfElectronsFromMultiCl_Closest_EoverCPenergy,
2994  loopSubFolders=False,
2995  purpose=PlotPurpose.Timing, page="EcalDrivenGsfElectronsFromMultiCl_Closest_EoverCPenergy"
2996  ))