CMS 3D CMS Logo
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 layerscheme = { 'lastLayerEEzm': 0, 'lastLayerFHzm': 0, 'maxlayerzm': 0, 'lastLayerEEzp': 0, 'lastLayerFHzp': 0, 'maxlayerzp': 0 }
25 
26 #Let's take the relevant values of layerscheme from the dqm file.
27 theDQMfile = "DQM_V0001_R000000001__Global__CMSSW_X_Y_Z__RECO.root"
28 if not os.path.isfile(theDQMfile):
29  print("Error: file", theDQMfile, "not found, exit")
30  sys.exit(0)
31 
32 #Take general info from the first file is sufficient.
33 thefile = TFile( theDQMfile )
34 GeneralInfoDirectory = 'DQMData/Run 1/HGCAL/Run summary/HGCalValidator/GeneralInfo'
35 
36 if not gDirectory.GetDirectory( GeneralInfoDirectory ):
37  print("Error: GeneralInfo directory not found in DQM file, exit")
38  sys.exit(0)
39 
40 keys = gDirectory.GetDirectory( GeneralInfoDirectory ).GetListOfKeys()
41 key = keys[0]
42 layvalue = 0
43 while key:
44  obj = key.ReadObj()
45  for laykey in layerscheme.keys():
46  if laykey in obj.GetName():
47  layvalue = obj.GetName()[len("<"+laykey+">i="):-len("</"+laykey+">")]
48  layerscheme[laykey] = layvalue
49  #print(layvalue)
50  key = keys.After(key)
51 
52 thefile.Close()
53 
54 print(layerscheme)
55 #TODO: Anticipating the fine/coarse layer information in CMSSW we overwrite values from DQM file
56 #For now values returned for
57 # 'lastLayerFHzp': '104', 'lastLayerFHzm': '52'
58 #are not the one expected. Will come back to this when there will be info in CMSSW to put in DQM file.
59 layerscheme = { 'lastLayerEEzm': 28, 'lastLayerFHzm': 40, 'maxlayerzm': 52, 'lastLayerEEzp': 80, 'lastLayerFHzp': 92, 'maxlayerzp': 104 }
60 print(layerscheme)
61 
62 lastLayerEEzm = layerscheme['lastLayerEEzm'] # last layer of EE -z
63 lastLayerFHzm = layerscheme['lastLayerFHzm'] # last layer of FH -z
64 maxlayerzm = layerscheme['maxlayerzm'] # last layer of BH -z
65 lastLayerEEzp = layerscheme['lastLayerEEzp'] # last layer of EE +z
66 lastLayerFHzp = layerscheme['lastLayerFHzp'] # last layer of FH +z
67 maxlayerzp = layerscheme['maxlayerzp'] # last layer of BH +z
68 
69 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
70 _legend_common = {"legendDx": -0.3,
71  "legendDy": -0.05,
72  "legendDw": 0.1}
73 
74 _SelectedCaloParticles = PlotGroup("SelectedCaloParticles", [
75  Plot("num_caloparticle_eta", xtitle="", **_common),
76  Plot("caloparticle_energy", xtitle="", **_common),
77  Plot("caloparticle_pt", xtitle="", **_common),
78  Plot("caloparticle_phi", xtitle="", **_common),
79  Plot("Eta vs Zorigin", xtitle="", **_common),
80  ])
81 
82 #Need to adjust the statbox to see better the plot
83 _common = {"stat": True, "drawStyle": "hist", "statx": 0.38, "staty": 0.68 }
84 _num_reco_cluster_eta = PlotGroup("num_reco_cluster_eta", [
85  Plot("num_reco_cluster_eta", xtitle="", **_common),
86 ],ncols=1)
87 #Back to normal
88 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
89 
90 _mixedhitsclusters = PlotGroup("mixedhitsclusters", [
91  Plot("mixedhitscluster_zminus", xtitle="", **_common),
92  Plot("mixedhitscluster_zplus", xtitle="", **_common),
93 ],ncols=2)
94 
95 #Just to prevent the stabox covering the plot
96 _common = {"stat": True, "drawStyle": "hist", "statx": 0.45, "staty": 0.65 }
97 
98 _energyclustered = PlotGroup("energyclustered", [
99  Plot("energyclustered_zminus", xtitle="", **_common),
100  Plot("energyclustered_zplus", xtitle="", **_common),
101 ],ncols=2)
102 
103 #Coming back to the usual box definition
104 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
105 
106 _longdepthbarycentre = PlotGroup("longdepthbarycentre", [
107  Plot("longdepthbarycentre_zminus", xtitle="", **_common),
108  Plot("longdepthbarycentre_zplus", xtitle="", **_common),
109 ],ncols=2)
110 
111 _common_layerperthickness = {}
112 _common_layerperthickness.update(_common)
113 _common_layerperthickness['xmin'] = 0.
114 _common_layerperthickness['xmax'] = 100
115 
116 _totclusternum_thick = PlotGroup("totclusternum_thick", [
117  Plot("totclusternum_thick_120", xtitle="", **_common_layerperthickness),
118  Plot("totclusternum_thick_200", xtitle="", **_common_layerperthickness),
119  Plot("totclusternum_thick_300", xtitle="", **_common_layerperthickness),
120  Plot("totclusternum_thick_-1", xtitle="", **_common_layerperthickness),
121  Plot("mixedhitscluster", xtitle="", **_common_layerperthickness),
122 ])
123 
124 #We will plot the density in logy scale.
125 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65}
126 
127 _cellsenedens_thick = PlotGroup("cellsenedens_thick", [
128  Plot("cellsenedens_thick_120", xtitle="", **_common),
129  Plot("cellsenedens_thick_200", xtitle="", **_common),
130  Plot("cellsenedens_thick_300", xtitle="", **_common),
131  Plot("cellsenedens_thick_-1", xtitle="", **_common),
132 ])
133 
134 #Coming back to the usual box definition
135 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
136 
137 
138 #--------------------------------------------------------------------------------------------
139 # z-
140 #--------------------------------------------------------------------------------------------
141 _totclusternum_layer_EE_zminus = PlotGroup("totclusternum_layer_EE_zminus", [
142  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
143 ], ncols=4)
144 
145 _totclusternum_layer_FH_zminus = PlotGroup("totclusternum_layer_FH_zminus", [
146  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
147 ], ncols=4)
148 
149 _totclusternum_layer_BH_zminus = PlotGroup("totclusternum_layer_BH_zminus", [
150  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
151 ], ncols=4)
152 
153 _energyclustered_perlayer_EE_zminus = PlotGroup("energyclustered_perlayer_EE_zminus", [
154  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
155 ], ncols=4)
156 
157 _energyclustered_perlayer_FH_zminus = PlotGroup("energyclustered_perlayer_FH_zminus", [
158  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
159 ], ncols=4)
160 
161 _energyclustered_perlayer_BH_zminus = PlotGroup("energyclustered_perlayer_BH_zminus", [
162  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
163 ], ncols=4)
164 
165 #----------------------------------------------------------------------------------------------------------------
166 #120 um
167 _common_cells = {}
168 _common_cells.update(_common)
169 _common_cells["xmin"] = 0
170 _common_cells["xmax"] = 50
171 _common_cells["ymin"] = 0.1
172 _common_cells["ymax"] = 10000
173 _common_cells["ylog"] = True
174 _cellsnum_perthick_perlayer_120_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_120_EE_zminus", [
175  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
176 ], ncols=4)
177 
178 _cellsnum_perthick_perlayer_120_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_120_FH_zminus", [
179  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
180 ], ncols=4)
181 
182 _cellsnum_perthick_perlayer_120_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_120_BH_zminus", [
183  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
184 ], ncols=4)
185 
186 #200 um
187 _cellsnum_perthick_perlayer_200_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_200_EE_zminus", [
188  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
189 ], ncols=4)
190 
191 _cellsnum_perthick_perlayer_200_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_200_FH_zminus", [
192  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
193 ], ncols=4)
194 
195 _cellsnum_perthick_perlayer_200_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_200_BH_zminus", [
196  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
197 ], ncols=4)
198 
199 #300 um
200 _cellsnum_perthick_perlayer_300_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_300_EE_zminus", [
201  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
202 ], ncols=4)
203 
204 _cellsnum_perthick_perlayer_300_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_300_FH_zminus", [
205  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
206 ], ncols=4)
207 
208 _cellsnum_perthick_perlayer_300_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_300_BH_zminus", [
209  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
210 ], ncols=4)
211 
212 #scint um
213 _cellsnum_perthick_perlayer_scint_EE_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_EE_zminus", [
214  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm)
215 ], ncols=4)
216 
217 _cellsnum_perthick_perlayer_scint_FH_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_FH_zminus", [
218  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzm,lastLayerFHzm)
219 ], ncols=4)
220 
221 _cellsnum_perthick_perlayer_scint_BH_zminus = PlotGroup("cellsnum_perthick_perlayer_-1_BH_zminus", [
222  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzm,maxlayerzm)
223 ], ncols=4)
224 
225 #----------------------------------------------------------------------------------------------------------------
226 #120 um
227 _common_distance = {}
228 _common_distance.update(_common)
229 _common_distance.update(_legend_common)
230 _common_distance["xmax"] = 150
231 _common_distance["stat"] = False
232 _common_distance["ymin"] = 1e-3
233 _common_distance["ymax"] = 10000
234 _common_distance["ylog"] = True
235 
236 _distancetomaxcell_perthickperlayer_120_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_EE_zminus", [
237  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
238 ], ncols=4)
239 
240 _distancetomaxcell_perthickperlayer_120_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_FH_zminus", [
241  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
242 ], ncols=4)
243 
244 _distancetomaxcell_perthickperlayer_120_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_120_BH_zminus", [
245  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
246 ], ncols=4)
247 
248 #200 um
249 _distancetomaxcell_perthickperlayer_200_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_EE_zminus", [
250  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
251 ], ncols=4)
252 
253 _distancetomaxcell_perthickperlayer_200_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_FH_zminus", [
254  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
255 ], ncols=4)
256 
257 _distancetomaxcell_perthickperlayer_200_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_200_BH_zminus", [
258  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
259 ], ncols=4)
260 
261 #300 um
262 _distancetomaxcell_perthickperlayer_300_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_EE_zminus", [
263  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
264 ], ncols=4)
265 
266 _distancetomaxcell_perthickperlayer_300_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_FH_zminus", [
267  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
268 ], ncols=4)
269 
270 _distancetomaxcell_perthickperlayer_300_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_300_BH_zminus", [
271  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
272 ], ncols=4)
273 
274 #scint um
275 _distancetomaxcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_EE_zminus", [
276  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
277 ], ncols=4)
278 
279 _distancetomaxcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_FH_zminus", [
280  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
281 ], ncols=4)
282 
283 _distancetomaxcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_-1_BH_zminus", [
284  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
285 ], ncols=4)
286 
287 #----------------------------------------------------------------------------------------------------------------
288 #120 um
289 _distancebetseedandmaxcell_perthickperlayer_120_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_EE_zminus", [
290  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
291 ], ncols=4)
292 
293 _distancebetseedandmaxcell_perthickperlayer_120_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_FH_zminus", [
294  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
295 ], ncols=4)
296 
297 _distancebetseedandmaxcell_perthickperlayer_120_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_BH_zminus", [
298  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
299 ], ncols=4)
300 
301 #200 um
302 _distancebetseedandmaxcell_perthickperlayer_200_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_EE_zminus", [
303  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
304 ], ncols=4)
305 
306 _distancebetseedandmaxcell_perthickperlayer_200_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_FH_zminus", [
307  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
308 ], ncols=4)
309 
310 _distancebetseedandmaxcell_perthickperlayer_200_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_BH_zminus", [
311  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
312 ], ncols=4)
313 
314 #300 um
315 _distancebetseedandmaxcell_perthickperlayer_300_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_EE_zminus", [
316  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
317 ], ncols=4)
318 
319 _distancebetseedandmaxcell_perthickperlayer_300_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_FH_zminus", [
320  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
321 ], ncols=4)
322 
323 _distancebetseedandmaxcell_perthickperlayer_300_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_BH_zminus", [
324  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
325 ], ncols=4)
326 
327 #scint um
328 _distancebetseedandmaxcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_EE_zminus", [
329  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
330 ], ncols=4)
331 
332 _distancebetseedandmaxcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_FH_zminus", [
333  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
334 ], ncols=4)
335 
336 _distancebetseedandmaxcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_BH_zminus", [
337  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
338 ], ncols=4)
339 
340 #----------------------------------------------------------------------------------------------------------------
341 #120 um
342 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus", [
343  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
344 ], ncols=4)
345 
346 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus", [
347  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
348 ], ncols=4)
349 
350 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus", [
351  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
352 ], ncols=4)
353 
354 #200 um
355 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus", [
356  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
357 ], ncols=4)
358 
359 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus", [
360  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
361 ], ncols=4)
362 
363 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus", [
364  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
365 ], ncols=4)
366 
367 #300 um
368 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus", [
369  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
370 ], ncols=4)
371 
372 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus", [
373  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
374 ], ncols=4)
375 
376 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus", [
377  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
378 ], ncols=4)
379 
380 #scint um
381 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_EE_zminus", [
382  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
383 ], ncols=4)
384 
385 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_FH_zminus", [
386  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
387 ], ncols=4)
388 
389 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zminus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_BH_zminus", [
390  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
391 ], ncols=4)
392 
393 #----------------------------------------------------------------------------------------------------------------
394 #120 um
395 _distancetoseedcell_perthickperlayer_120_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_EE_zminus", [
396  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
397 ], ncols=4)
398 
399 _distancetoseedcell_perthickperlayer_120_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_FH_zminus", [
400  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
401 ], ncols=4)
402 
403 _distancetoseedcell_perthickperlayer_120_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_120_BH_zminus", [
404  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
405 ], ncols=4)
406 
407 #200 um
408 _distancetoseedcell_perthickperlayer_200_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_EE_zminus", [
409  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
410 ], ncols=4)
411 
412 _distancetoseedcell_perthickperlayer_200_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_FH_zminus", [
413  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
414 ], ncols=4)
415 
416 _distancetoseedcell_perthickperlayer_200_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_200_BH_zminus", [
417  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
418 ], ncols=4)
419 
420 #300 um
421 _distancetoseedcell_perthickperlayer_300_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_EE_zminus", [
422  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
423 ], ncols=4)
424 
425 _distancetoseedcell_perthickperlayer_300_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_FH_zminus", [
426  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
427 ], ncols=4)
428 
429 _distancetoseedcell_perthickperlayer_300_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_300_BH_zminus", [
430  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
431 ], ncols=4)
432 
433 #scint um
434 _distancetoseedcell_perthickperlayer_scint_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_EE_zminus", [
435  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
436 ], ncols=4)
437 
438 _distancetoseedcell_perthickperlayer_scint_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_FH_zminus", [
439  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
440 ], ncols=4)
441 
442 _distancetoseedcell_perthickperlayer_scint_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_-1_BH_zminus", [
443  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
444 ], ncols=4)
445 
446 #=====================================================================================================================
447 #----------------------------------------------------------------------------------------------------------------
448 #We need points for the weighted plots
449 _common = {"stat": True, "drawStyle": "EP", "staty": 0.65 }
450 #120 um
451 _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus", [
452  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
453 ], ncols=4)
454 
455 _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus", [
456  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
457 ], ncols=4)
458 
459 _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus", [
460  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
461 ], ncols=4)
462 
463 #200 um
464 _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus", [
465  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
466 ], ncols=4)
467 
468 _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus", [
469  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
470 ], ncols=4)
471 
472 _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus", [
473  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
474 ], ncols=4)
475 
476 #300 um
477 _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus", [
478  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
479 ], ncols=4)
480 
481 _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus", [
482  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
483 ], ncols=4)
484 
485 _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus", [
486  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
487 ], ncols=4)
488 #scint um
489 _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_EE_zminus", [
490  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
491 ], ncols=4)
492 
493 _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_FH_zminus", [
494  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
495 ], ncols=4)
496 
497 _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zminus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_BH_zminus", [
498  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
499 ], ncols=4)
500 
501 
502 #----------------------------------------------------------------------------------------------------------------
503 #120 um
504 _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus", [
505  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
506 ], ncols=4)
507 
508 _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus", [
509  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
510 ], ncols=4)
511 
512 _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus", [
513  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
514 ], ncols=4)
515 
516 #200 um
517 _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus", [
518  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
519 ], ncols=4)
520 
521 _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus", [
522  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
523 ], ncols=4)
524 
525 _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus", [
526  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
527 ], ncols=4)
528 
529 #300 um
530 _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus", [
531  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
532 ], ncols=4)
533 
534 _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus", [
535  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
536 ], ncols=4)
537 
538 _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus", [
539  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
540 ], ncols=4)
541 
542 #scint um
543 _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_EE_zminus", [
544  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm)
545 ], ncols=4)
546 
547 _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_FH_zminus", [
548  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzm,lastLayerFHzm)
549 ], ncols=4)
550 
551 _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zminus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_BH_zminus", [
552  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzm,maxlayerzm)
553 ], ncols=4)
554 
555 #Coming back to the usual definition
556 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
557 
558 #--------------------------------------------------------------------------------------------
559 # z+
560 #--------------------------------------------------------------------------------------------
561 _totclusternum_layer_EE_zplus = PlotGroup("totclusternum_layer_EE_zplus", [
562  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
563 ], ncols=4)
564 
565 _totclusternum_layer_FH_zplus = PlotGroup("totclusternum_layer_FH_zplus", [
566  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
567 ], ncols=4)
568 
569 _totclusternum_layer_BH_zplus = PlotGroup("totclusternum_layer_BH_zplus", [
570  Plot("totclusternum_layer_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
571 ], ncols=4)
572 
573 _energyclustered_perlayer_EE_zplus = PlotGroup("energyclustered_perlayer_EE_zplus", [
574  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
575 ], ncols=4)
576 
577 _energyclustered_perlayer_FH_zplus = PlotGroup("energyclustered_perlayer_FH_zplus", [
578  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
579 ], ncols=4)
580 
581 _energyclustered_perlayer_BH_zplus = PlotGroup("energyclustered_perlayer_BH_zplus", [
582  Plot("energyclustered_perlayer{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
583 ], ncols=4)
584 
585 #----------------------------------------------------------------------------------------------------------------
586 #120 um
587 _cellsnum_perthick_perlayer_120_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_120_EE_zplus", [
588  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
589 ], ncols=4)
590 
591 _cellsnum_perthick_perlayer_120_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_120_FH_zplus", [
592  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
593 ], ncols=4)
594 _cellsnum_perthick_perlayer_120_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_120_BH_zplus", [
595  Plot("cellsnum_perthick_perlayer_120_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
596 ], ncols=4)
597 
598 #200 um
599 _cellsnum_perthick_perlayer_200_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_200_EE_zplus", [
600  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
601 ], ncols=4)
602 
603 _cellsnum_perthick_perlayer_200_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_200_FH_zplus", [
604  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
605 ], ncols=4)
606 
607 _cellsnum_perthick_perlayer_200_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_200_BH_zplus", [
608  Plot("cellsnum_perthick_perlayer_200_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
609 ], ncols=4)
610 #300 um
611 _cellsnum_perthick_perlayer_300_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_300_EE_zplus", [
612  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
613 ], ncols=4)
614 
615 _cellsnum_perthick_perlayer_300_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_300_FH_zplus", [
616  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
617 ], ncols=4)
618 _cellsnum_perthick_perlayer_300_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_300_BH_zplus", [
619  Plot("cellsnum_perthick_perlayer_300_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
620 ], ncols=4)
621 
622 #scint um
623 _cellsnum_perthick_perlayer_scint_EE_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_EE_zplus", [
624  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(maxlayerzm,lastLayerEEzp)
625 ], ncols=4)
626 
627 _cellsnum_perthick_perlayer_scint_FH_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_FH_zplus", [
628  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerEEzp,lastLayerFHzp)
629 ], ncols=4)
630 
631 _cellsnum_perthick_perlayer_scint_BH_zplus = PlotGroup("cellsnum_perthick_perlayer_-1_BH_zplus", [
632  Plot("cellsnum_perthick_perlayer_-1_{:02d}".format(i), xtitle="", **_common_cells) for i in range(lastLayerFHzp,maxlayerzp)
633 ], ncols=4)
634 
635 #----------------------------------------------------------------------------------------------------------------
636 #120 um
637 _common_distance = {}
638 _common_distance.update(_common)
639 _common_distance.update(_legend_common)
640 _common_distance["xmax"] = 150
641 _common_distance["stat"] = False
642 _common_distance["ymin"] = 1e-3
643 _common_distance["ymax"] = 10000
644 _common_distance["ylog"] = True
645 
646 _distancetomaxcell_perthickperlayer_120_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_EE_zplus", [
647  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
648 ], ncols=4)
649 
650 _distancetomaxcell_perthickperlayer_120_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_FH_zplus", [
651  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
652 ], ncols=4)
653 
654 _distancetomaxcell_perthickperlayer_120_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_120_BH_zplus", [
655  Plot("distancetomaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
656 ], ncols=4)
657 
658 #200 um
659 _distancetomaxcell_perthickperlayer_200_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_EE_zplus", [
660  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
661 ], ncols=4)
662 
663 _distancetomaxcell_perthickperlayer_200_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_FH_zplus", [
664  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
665 ], ncols=4)
666 
667 _distancetomaxcell_perthickperlayer_200_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_200_BH_zplus", [
668  Plot("distancetomaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
669 ], ncols=4)
670 
671 #300 um
672 _distancetomaxcell_perthickperlayer_300_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_EE_zplus", [
673  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
674 ], ncols=4)
675 
676 _distancetomaxcell_perthickperlayer_300_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_FH_zplus", [
677  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
678 ], ncols=4)
679 
680 _distancetomaxcell_perthickperlayer_300_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_300_BH_zplus", [
681  Plot("distancetomaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
682 ], ncols=4)
683 
684 #scint um
685 _distancetomaxcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_EE_zplus", [
686  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
687 ], ncols=4)
688 
689 _distancetomaxcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_FH_zplus", [
690  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
691 ], ncols=4)
692 
693 _distancetomaxcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_-1_BH_zplus", [
694  Plot("distancetomaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
695 ], ncols=4)
696 
697 #----------------------------------------------------------------------------------------------------------------
698 #120 um
699 _distancebetseedandmaxcell_perthickperlayer_120_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_EE_zplus", [
700  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
701 ], ncols=4)
702 
703 _distancebetseedandmaxcell_perthickperlayer_120_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_FH_zplus", [
704  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
705 ], ncols=4)
706 
707 _distancebetseedandmaxcell_perthickperlayer_120_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_120_BH_zplus", [
708  Plot("distancebetseedandmaxcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
709 ], ncols=4)
710 
711 #200 um
712 _distancebetseedandmaxcell_perthickperlayer_200_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_EE_zplus", [
713  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
714 ], ncols=4)
715 
716 _distancebetseedandmaxcell_perthickperlayer_200_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_FH_zplus", [
717  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
718 ], ncols=4)
719 
720 _distancebetseedandmaxcell_perthickperlayer_200_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_200_BH_zplus", [
721  Plot("distancebetseedandmaxcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
722 ], ncols=4)
723 
724 #300 um
725 _distancebetseedandmaxcell_perthickperlayer_300_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_EE_zplus", [
726  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
727 ], ncols=4)
728 
729 _distancebetseedandmaxcell_perthickperlayer_300_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_FH_zplus", [
730  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
731 ], ncols=4)
732 
733 _distancebetseedandmaxcell_perthickperlayer_300_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_300_BH_zplus", [
734  Plot("distancebetseedandmaxcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
735 ], ncols=4)
736 
737 #scint um
738 _distancebetseedandmaxcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_EE_zplus", [
739  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
740 ], ncols=4)
741 
742 _distancebetseedandmaxcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_FH_zplus", [
743  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
744 ], ncols=4)
745 
746 _distancebetseedandmaxcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancebetseedandmaxcell_perthickperlayer_-1_BH_zplus", [
747  Plot("distancebetseedandmaxcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
748 ], ncols=4)
749 
750 #----------------------------------------------------------------------------------------------------------------
751 #120 um
752 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus", [
753  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
754 ], ncols=4)
755 
756 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus", [
757  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
758 ], ncols=4)
759 
760 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus", [
761  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
762 ], ncols=4)
763 
764 #200 um
765 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus", [
766  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
767 ], ncols=4)
768 
769 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus", [
770  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
771 ], ncols=4)
772 
773 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus", [
774  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
775 ], ncols=4)
776 
777 #300 um
778 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus", [
779  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
780 ], ncols=4)
781 
782 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus", [
783  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
784 ], ncols=4)
785 
786 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus", [
787  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
788 ], ncols=4)
789 
790 #scint um
791 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_EE_zplus", [
792  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
793 ], ncols=4)
794 
795 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_FH_zplus", [
796  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
797 ], ncols=4)
798 
799 _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zplus = PlotGroup("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_BH_zplus", [
800  Plot("distancebetseedandmaxcellvsclusterenergy_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
801 ], ncols=4)
802 
803 
804 #----------------------------------------------------------------------------------------------------------------
805 #120 um
806 _distancetoseedcell_perthickperlayer_120_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_EE_zplus", [
807  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
808 ], ncols=4)
809 
810 _distancetoseedcell_perthickperlayer_120_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_FH_zplus", [
811  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
812 ], ncols=4)
813 
814 _distancetoseedcell_perthickperlayer_120_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_120_BH_zplus", [
815  Plot("distancetoseedcell_perthickperlayer_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
816 ], ncols=4)
817 
818 #200 um
819 _distancetoseedcell_perthickperlayer_200_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_EE_zplus", [
820  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
821 ], ncols=4)
822 
823 _distancetoseedcell_perthickperlayer_200_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_FH_zplus", [
824  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
825 ], ncols=4)
826 
827 _distancetoseedcell_perthickperlayer_200_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_200_BH_zplus", [
828  Plot("distancetoseedcell_perthickperlayer_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
829 ], ncols=4)
830 
831 #300 um
832 _distancetoseedcell_perthickperlayer_300_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_EE_zplus", [
833  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
834 ], ncols=4)
835 
836 _distancetoseedcell_perthickperlayer_300_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_FH_zplus", [
837  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
838 ], ncols=4)
839 
840 _distancetoseedcell_perthickperlayer_300_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_300_BH_zplus", [
841  Plot("distancetoseedcell_perthickperlayer_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
842 ], ncols=4)
843 
844 #scint um
845 _distancetoseedcell_perthickperlayer_scint_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_EE_zplus", [
846  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
847 ], ncols=4)
848 
849 _distancetoseedcell_perthickperlayer_scint_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_FH_zplus", [
850  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
851 ], ncols=4)
852 
853 _distancetoseedcell_perthickperlayer_scint_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_-1_BH_zplus", [
854  Plot("distancetoseedcell_perthickperlayer_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
855 ], ncols=4)
856 
857 #=====================================================================================================================
858 #----------------------------------------------------------------------------------------------------------------
859 #We need points for the weighted plots
860 _common = {"stat": True, "drawStyle": "EP", "staty": 0.65 }
861 
862 #120 um
863 _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus", [
864  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
865 ], ncols=4)
866 
867 _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus", [
868  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
869 ], ncols=4)
870 
871 _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus", [
872  Plot("distancetomaxcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
873 ], ncols=4)
874 
875 #200 um
876 _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus", [
877  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
878 ], ncols=4)
879 _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus", [
880  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
881 ], ncols=4)
882 
883 _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus", [
884  Plot("distancetomaxcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
885 ], ncols=4)
886 
887 #300 um
888 _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus", [
889  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
890 ], ncols=4)
891 
892 _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus", [
893  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
894 ], ncols=4)
895 
896 _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus", [
897  Plot("distancetomaxcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
898 ], ncols=4)
899 
900 #scint um
901 _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_EE_zplus", [
902  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
903 ], ncols=4)
904 
905 _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_FH_zplus", [
906  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
907 ], ncols=4)
908 
909 _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zplus = PlotGroup("distancetomaxcell_perthickperlayer_eneweighted_-1_BH_zplus", [
910  Plot("distancetomaxcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
911 ], ncols=4)
912 
913 #----------------------------------------------------------------------------------------------------------------
914 #120 um
915 _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus", [
916  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
917 ], ncols=4)
918 
919 _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus", [
920  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
921 ], ncols=4)
922 
923 _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus", [
924  Plot("distancetoseedcell_perthickperlayer_eneweighted_120_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
925 ], ncols=4)
926 
927 #200 um
928 _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus", [
929  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
930 ], ncols=4)
931 
932 _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus", [
933  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
934 ], ncols=4)
935 
936 _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus", [
937  Plot("distancetoseedcell_perthickperlayer_eneweighted_200_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
938 ], ncols=4)
939 
940 #300 um
941 _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus", [
942  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
943 ], ncols=4)
944 
945 _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus", [
946  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
947 ], ncols=4)
948 
949 _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus", [
950  Plot("distancetoseedcell_perthickperlayer_eneweighted_300_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
951 ], ncols=4)
952 
953 #scint um
954 _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_EE_zplus", [
955  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(maxlayerzm,lastLayerEEzp)
956 ], ncols=4)
957 
958 _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_FH_zplus", [
959  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerEEzp,lastLayerFHzp)
960 ], ncols=4)
961 
962 _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zplus = PlotGroup("distancetoseedcell_perthickperlayer_eneweighted_-1_BH_zplus", [
963  Plot("distancetoseedcell_perthickperlayer_eneweighted_-1_{:02d}".format(i), xtitle="", **_common) for i in range(lastLayerFHzp,maxlayerzp)
964 ], ncols=4)
965 #Just in case we add some plots below to be on the safe side.
966 _common = {"stat": True, "drawStyle": "hist", "staty": 0.65 }
967 
968 #--------------------------------------------------------------------------------------------
969 # z-
970 #--------------------------------------------------------------------------------------------
971 
972 _common_score = {"title": "Score CaloParticle to LayerClusters in z-",
973  "stat": False,
974  "ymin": 0.1,
975  "ymax": 1000,
976  "xmin": 0,
977  "xmax": 1,
978  "drawStyle": "hist",
979  "lineWidth": 1,
980  "ylog": True
981  }
982 _common_score.update(_legend_common)
983 _score_caloparticle_to_layerclusters_zminus = PlotGroup("score_caloparticle_to_layercluster_zminus", [
984  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)
985  ], ncols=8 )
986 
987 _common_score = {"title": "Score LayerCluster to CaloParticles in z-",
988  "stat": False,
989  "ymin": 0.1,
990  "ymax": 1000,
991  "xmin": 0,
992  "xmax": 1,
993  "drawStyle": "hist",
994  "lineWidth": 1,
995  "ylog": True
996  }
997 _common_score.update(_legend_common)
998 _score_layercluster_to_caloparticles_zminus = PlotGroup("score_layercluster_to_caloparticle_zminus", [
999  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)
1000  ], ncols=8 )
1001 
1002 _common_shared= {"title": "Shared Energy CaloParticle To Layer Cluster in z-",
1003  "stat": False,
1004  "legend": False,
1005  }
1006 _common_shared.update(_legend_common)
1007 _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)]
1008 _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)])
1009 _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)])
1010 _sharedEnergy_caloparticle_to_layercluster_zminus = PlotGroup("sharedEnergy_caloparticle_to_layercluster_zminus", _shared_plots_zminus, ncols=8)
1011 
1012 _common_shared= {"title": "Shared Energy Layer Cluster To CaloParticle in z-",
1013  "stat": False,
1014  "legend": False,
1015  }
1016 _common_shared.update(_legend_common)
1017 _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)]
1018 _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)])
1019 _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)])
1020 _sharedEnergy_layercluster_to_caloparticle_zminus = PlotGroup("sharedEnergy_layercluster_to_caloparticle_zminus", _shared_plots2_zminus, ncols=8)
1021 
1022 
1023 _common_assoc = {#"title": "Cell Association Table in z-",
1024  "stat": False,
1025  "legend": False,
1026  "xbinlabels": ["", "TN(pur)", "FN(ineff.)", "FP(fake)", "TP(eff)"],
1027  "xbinlabeloption": "h",
1028  "drawStyle": "hist",
1029  "ymin": 0.1,
1030  "ymax": 10000,
1031  "ylog": True}
1032 _common_assoc.update(_legend_common)
1033 _cell_association_table_zminus = PlotGroup("cellAssociation_table_zminus", [
1034  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)
1035  ], ncols=8 )
1036 
1037 _bin_count = 0
1038 _xbinlabels = [ "Layer {:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1039 _common_eff = {"stat": False, "legend": False, "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1040 _effplots_zminus = [Plot("effic_eta_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(0,maxlayerzm)]
1041 _effplots_zminus.extend([Plot("effic_phi_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(0,maxlayerzm)])
1042 _common_eff["xmin"] = 0.
1043 _bin_count += maxlayerzm
1044 _common_eff["xmax"] =_bin_count
1045 _effplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Efficiencies in z-", **_common_eff)])
1046 _efficiencies_zminus = PlotGroup("Efficiencies_zminus", _effplots_zminus, ncols=8)
1047 
1048 
1049 _common_dup = {"stat": False, "legend": False, "title": "Global Duplicates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1050 _dupplots_zminus = [Plot("duplicate_eta_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(0,maxlayerzm)]
1051 _dupplots_zminus.extend([Plot("duplicate_phi_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(0,maxlayerzm)])
1052 _common_dup["xmin"] = _bin_count+maxlayerzm+1
1053 _bin_count += maxlayerzp
1054 _common_dup["xmax"] = _bin_count
1055 _dupplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Duplicates in z-", **_common_dup)])
1056 _duplicates_zminus = PlotGroup("Duplicates_zminus", _dupplots_zminus, ncols=8)
1057 
1058 _common_fake = {"stat": False, "legend": False, "title": "Global Fake Rates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1059 _fakeplots_zminus = [Plot("fake_eta_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(0,maxlayerzm)]
1060 _fakeplots_zminus.extend([Plot("fake_phi_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(0,maxlayerzm)])
1061 _common_fake["xmin"] = _bin_count+maxlayerzm+1
1062 _bin_count += maxlayerzp
1063 _common_fake["xmax"] = _bin_count
1064 _common_fake["xbinlabels"] = [ "Layer {:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1065 _common_fake["xbinlabelsize"] = 10.
1066 _fakeplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global Fake Rate in z-", **_common_fake)])
1067 _fakes_zminus = PlotGroup("FakeRate_zminus", _fakeplots_zminus, ncols=8)
1068 
1069 _common_merge = {"stat": False, "legend": False, "title": "Global Merge Rates in z-", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1070 _mergeplots_zminus = [Plot("merge_eta_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(0,maxlayerzm)]
1071 _mergeplots_zminus.extend([Plot("merge_phi_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(0,maxlayerzm)])
1072 _common_merge["xmin"] = _bin_count+maxlayerzm+1
1073 _bin_count += maxlayerzp
1074 _common_merge["xmax"] = _bin_count
1075 _common_merge["xbinlabels"] = [ "Layer {:02d}".format(i+1) for i in range(0,maxlayerzm) ]
1076 _common_merge["xbinlabelsize"] = 10.
1077 _mergeplots_zminus.extend([Plot("globalEfficiencies", xtitle="Global merge Rate in z-", **_common_merge)])
1078 _merges_zminus = PlotGroup("MergeRate_zminus", _mergeplots_zminus, ncols=8)
1079 
1080 
1081 _common_energy_score = dict(removeEmptyBins=True, xbinlabelsize=10, xbinlabeloption="d", ncols=4)
1082 _energyscore_cp2lc_zminus = []
1083 for i in range(0, maxlayerzm):
1084  _energyscore_cp2lc_zminus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_caloparticle2layer_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginLeft=0.1, adjustMarginRight=0.1, **_common_energy_score)))
1085 
1086 _energyscore_lc2cp_zminus = []
1087 for i in range(0, maxlayerzm):
1088  _energyscore_lc2cp_zminus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_layer2caloparticle_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginLeft=0.1, adjustMarginRight=0.1, **_common_energy_score)))
1089 #_energyclustered =
1090 
1091 #--------------------------------------------------------------------------------------------
1092 # z+
1093 #--------------------------------------------------------------------------------------------
1094 _common_score = {"title": "Score CaloParticle to LayerClusters in z+",
1095  "stat": False,
1096  "ymin": 0.1,
1097  "ymax": 1000,
1098  "xmin": 0,
1099  "xmax": 1,
1100  "drawStyle": "hist",
1101  "lineWidth": 1,
1102  "ylog": True
1103  }
1104 _common_score.update(_legend_common)
1105 _score_caloparticle_to_layerclusters_zplus = PlotGroup("score_caloparticle_to_layercluster_zplus", [
1106  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)
1107  ], ncols=8 )
1108 
1109 _common_score = {"title": "Score LayerCluster to CaloParticles in z+",
1110  "stat": False,
1111  "ymin": 0.1,
1112  "ymax": 1000,
1113  "xmin": 0,
1114  "xmax": 1,
1115  "drawStyle": "hist",
1116  "lineWidth": 1,
1117  "ylog": True
1118  }
1119 _common_score.update(_legend_common)
1120 _score_layercluster_to_caloparticles_zplus = PlotGroup("score_layercluster_to_caloparticle_zplus", [
1121  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)
1122  ], ncols=8 )
1123 
1124 _common_shared= {"title": "Shared Energy CaloParticle To Layer Cluster in z+",
1125  "stat": False,
1126  "legend": False,
1127  }
1128 _common_shared.update(_legend_common)
1129 _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)]
1130 _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)])
1131 _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)])
1132 _sharedEnergy_caloparticle_to_layercluster_zplus = PlotGroup("sharedEnergy_caloparticle_to_layercluster_zplus", _shared_plots_zplus, ncols=8)
1133 
1134 _common_shared= {"title": "Shared Energy Layer Cluster To CaloParticle in z+",
1135  "stat": False,
1136  "legend": False,
1137  }
1138 _common_shared.update(_legend_common)
1139 _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)]
1140 _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)])
1141 _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)])
1142 _sharedEnergy_layercluster_to_caloparticle_zplus = PlotGroup("sharedEnergy_layercluster_to_caloparticle_zplus", _shared_plots2_zplus, ncols=8)
1143 
1144 
1145 _common_assoc = {#"title": "Cell Association Table in z+",
1146  "stat": False,
1147  "legend": False,
1148  "xbinlabels": ["", "TN(pur)", "FN(ineff.)", "FP(fake)", "TP(eff)"],
1149  "xbinlabeloption": "h",
1150  "drawStyle": "hist",
1151  "ymin": 0.1,
1152  "ymax": 10000,
1153  "ylog": True}
1154 _common_assoc.update(_legend_common)
1155 _cell_association_table_zplus = PlotGroup("cellAssociation_table_zplus", [
1156  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)
1157  ], ncols=8 )
1158 
1159 _bin_count = 0
1160 _common_eff = {"stat": False, "legend": False, "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1161 _effplots_zplus = [Plot("effic_eta_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(maxlayerzm,maxlayerzp)]
1162 _effplots_zplus.extend([Plot("effic_phi_layer{:02d}".format(i), xtitle="", **_common_eff) for i in range(maxlayerzm,maxlayerzp)])
1163 _common_eff["xmin"] = maxlayerzm
1164 _bin_count += maxlayerzp
1165 _common_eff["xmax"] =_bin_count
1166 _effplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Efficiencies in z+", **_common_eff)])
1167 _efficiencies_zplus = PlotGroup("Efficiencies_zplus", _effplots_zplus, ncols=8)
1168 
1169 
1170 _common_dup = {"stat": False, "legend": False, "title": "Global Duplicates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1171 _dupplots_zplus = [Plot("duplicate_eta_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(maxlayerzm,maxlayerzp)]
1172 _dupplots_zplus.extend([Plot("duplicate_phi_layer{:02d}".format(i), xtitle="", **_common_dup) for i in range(maxlayerzm,maxlayerzp)])
1173 _common_dup["xmin"] = _bin_count+maxlayerzm+1
1174 _bin_count += maxlayerzp
1175 _common_dup["xmax"] = _bin_count
1176 _dupplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Duplicates in z+", **_common_dup)])
1177 _duplicates_zplus = PlotGroup("Duplicates_zplus", _dupplots_zplus, ncols=8)
1178 
1179 _common_fake = {"stat": False, "legend": False, "title": "Global Fake Rates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1180 _fakeplots_zplus = [Plot("fake_eta_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(maxlayerzm,maxlayerzp)]
1181 _fakeplots_zplus.extend([Plot("fake_phi_layer{:02d}".format(i), xtitle="", **_common_fake) for i in range(maxlayerzm,maxlayerzp)])
1182 _common_fake["xmin"] = _bin_count+maxlayerzm+1
1183 _bin_count += maxlayerzp
1184 _common_fake["xmax"] = _bin_count
1185 _fakeplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global Fake Rate in z+", **_common_fake)])
1186 _fakes_zplus = PlotGroup("FakeRate_zplus", _fakeplots_zplus, ncols=8)
1187 
1188 _common_merge = {"stat": False, "legend": False, "title": "Global Merge Rates in z+", "xbinlabels": _xbinlabels, "xbinlabelsize": 12, "xbinlabeloptions": "v"}
1189 _mergeplots_zplus = [Plot("merge_eta_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(maxlayerzm,maxlayerzp)]
1190 _mergeplots_zplus.extend([Plot("merge_phi_layer{:02d}".format(i), xtitle="", **_common_merge) for i in range(maxlayerzm,maxlayerzp)])
1191 _common_merge["xmin"] = _bin_count+maxlayerzm+1
1192 _bin_count += maxlayerzp
1193 _common_merge["xmax"] = _bin_count
1194 _mergeplots_zplus.extend([Plot("globalEfficiencies", xtitle="Global merge Rate in z+", **_common_merge)])
1195 _merges_zplus = PlotGroup("MergeRate_zplus", _mergeplots_zplus, ncols=8)
1196 
1197 
1198 _common_energy_score = dict(removeEmptyBins=False, xbinlabelsize=10, xbinlabeloption="d", ncols=4)
1199 _energyscore_cp2lc_zplus = []
1200 for i in range(maxlayerzm,maxlayerzp):
1201  _energyscore_cp2lc_zplus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_caloparticle2layer_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginLeft=0.1, adjustMarginRight=0.1, **_common_energy_score)))
1202 
1203 _energyscore_lc2cp_zplus = []
1204 for i in range(maxlayerzm,maxlayerzp):
1205  _energyscore_lc2cp_zplus.append(PlotOnSideGroup("Energy_vs_Score_Layer{:02d}".format(i), Plot("Energy_vs_Score_layer2caloparticle_perlayer{:02d}".format(i), drawStyle="COLZ", adjustMarginLeft=0.1, adjustMarginRight=0.1, **_common_energy_score)))
1206 #_energyclustered =
1207 
1208 #=================================================================================================
1209 hgcalLayerClustersPlotter = Plotter()
1210 #We follow Chris categories in folders
1211 # [A] calculated "energy density" for cells in a) 120um, b) 200um, c) 300um, d) scint
1212 # (one entry per rechit, in the appropriate histo)
1213 hgcalLayerClustersPlotter.append("CellsEnergyDensityPerThickness", [
1214  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1215  ], PlotFolder(
1216  _cellsenedens_thick,
1217  loopSubFolders=False,
1218  purpose=PlotPurpose.Timing, page="CellsEnergyDensityPerThickness"
1219  ))
1220 
1221 # [B] number of layer clusters per event in a) 120um, b) 200um, c) 300um, d) scint
1222 # (one entry per event in each of the four histos)
1223 hgcalLayerClustersPlotter.append("TotalNumberofLayerClustersPerThickness", [
1224  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1225  ], PlotFolder(
1226  _totclusternum_thick,
1227  loopSubFolders=False,
1228  purpose=PlotPurpose.Timing, page="TotalNumberofLayerClustersPerThickness"
1229  ))
1230 
1231 # [C] number of layer clusters per layer (one entry per event in each histo)
1232 # z-
1233 hgcalLayerClustersPlotter.append("NumberofLayerClustersPerLayer_zminus", [
1234  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1235  ], PlotFolder(
1236  _totclusternum_layer_EE_zminus,
1237  _totclusternum_layer_FH_zminus,
1238  _totclusternum_layer_BH_zminus,
1239  loopSubFolders=False,
1240  purpose=PlotPurpose.Timing, page="NumberofLayerClustersPerLayer_zminus"
1241  ))
1242 
1243 # z+
1244 hgcalLayerClustersPlotter.append("NumberofLayerClustersPerLayer_zplus", [
1245  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1246  ], PlotFolder(
1247  _totclusternum_layer_EE_zplus,
1248  _totclusternum_layer_FH_zplus,
1249  _totclusternum_layer_BH_zplus,
1250  loopSubFolders=False,
1251  purpose=PlotPurpose.Timing, page="NumberofLayerClustersPerLayer_zplus"
1252  ))
1253 
1254 # [D] For each layer cluster:
1255 # number of cells in layer cluster, by layer - separate histos in each layer for 120um Si, 200/300um Si, Scint
1256 # NB: not all combinations exist; e.g. no 120um Si in layers with scint.
1257 # (One entry in the appropriate histo per layer cluster).
1258 # z-
1259 hgcalLayerClustersPlotter.append("CellsNumberPerLayerPerThickness_zminus", [
1260  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1261  ], PlotFolder(
1262  _cellsnum_perthick_perlayer_120_EE_zminus,
1263  _cellsnum_perthick_perlayer_120_FH_zminus,
1264  _cellsnum_perthick_perlayer_120_BH_zminus,
1265  _cellsnum_perthick_perlayer_200_EE_zminus,
1266  _cellsnum_perthick_perlayer_200_FH_zminus,
1267  _cellsnum_perthick_perlayer_200_BH_zminus,
1268  _cellsnum_perthick_perlayer_300_EE_zminus,
1269  _cellsnum_perthick_perlayer_300_FH_zminus,
1270  _cellsnum_perthick_perlayer_300_BH_zminus,
1271  _cellsnum_perthick_perlayer_scint_EE_zminus,
1272  _cellsnum_perthick_perlayer_scint_FH_zminus,
1273  _cellsnum_perthick_perlayer_scint_BH_zminus,
1274  loopSubFolders=False,
1275  purpose=PlotPurpose.Timing, page="CellsNumberPerLayerPerThickness_zminus"
1276  ))
1277 
1278 # z+
1279 hgcalLayerClustersPlotter.append("CellsNumberPerLayerPerThickness_zplus", [
1280  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1281  ], PlotFolder(
1282  _cellsnum_perthick_perlayer_120_EE_zplus,
1283  _cellsnum_perthick_perlayer_120_FH_zplus,
1284  _cellsnum_perthick_perlayer_120_BH_zplus,
1285  _cellsnum_perthick_perlayer_200_EE_zplus,
1286  _cellsnum_perthick_perlayer_200_FH_zplus,
1287  _cellsnum_perthick_perlayer_200_BH_zplus,
1288  _cellsnum_perthick_perlayer_300_EE_zplus,
1289  _cellsnum_perthick_perlayer_300_FH_zplus,
1290  _cellsnum_perthick_perlayer_300_BH_zplus,
1291  _cellsnum_perthick_perlayer_scint_EE_zplus,
1292  _cellsnum_perthick_perlayer_scint_FH_zplus,
1293  _cellsnum_perthick_perlayer_scint_BH_zplus,
1294  loopSubFolders=False,
1295  purpose=PlotPurpose.Timing, page="CellsNumberPerLayerPerThickness_zplus"
1296  ))
1297 
1298 # [E] For each layer cluster:
1299 # distance of cells from a) seed cell, b) max cell; and c), d): same with entries weighted by cell energy
1300 # separate histos in each layer for 120um Si, 200/300um Si, Scint
1301 # NB: not all combinations exist; e.g. no 120um Si in layers with scint.
1302 # (One entry in each of the four appropriate histos per cell in a layer cluster)
1303 # z-
1304 hgcalLayerClustersPlotter.append("CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zminus", [
1305  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1306  ], PlotFolder(
1307  _distancetomaxcell_perthickperlayer_120_EE_zminus,
1308  _distancetomaxcell_perthickperlayer_120_FH_zminus,
1309  _distancetomaxcell_perthickperlayer_120_BH_zminus,
1310  _distancetomaxcell_perthickperlayer_200_EE_zminus,
1311  _distancetomaxcell_perthickperlayer_200_FH_zminus,
1312  _distancetomaxcell_perthickperlayer_200_BH_zminus,
1313  _distancetomaxcell_perthickperlayer_300_EE_zminus,
1314  _distancetomaxcell_perthickperlayer_300_FH_zminus,
1315  _distancetomaxcell_perthickperlayer_300_BH_zminus,
1316  _distancetomaxcell_perthickperlayer_scint_EE_zminus,
1317  _distancetomaxcell_perthickperlayer_scint_FH_zminus,
1318  _distancetomaxcell_perthickperlayer_scint_BH_zminus,
1319  _distancetoseedcell_perthickperlayer_120_EE_zminus,
1320  _distancetoseedcell_perthickperlayer_120_FH_zminus,
1321  _distancetoseedcell_perthickperlayer_120_BH_zminus,
1322  _distancetoseedcell_perthickperlayer_200_EE_zminus,
1323  _distancetoseedcell_perthickperlayer_200_FH_zminus,
1324  _distancetoseedcell_perthickperlayer_200_BH_zminus,
1325  _distancetoseedcell_perthickperlayer_300_EE_zminus,
1326  _distancetoseedcell_perthickperlayer_300_FH_zminus,
1327  _distancetoseedcell_perthickperlayer_300_BH_zminus,
1328  _distancetoseedcell_perthickperlayer_scint_EE_zminus,
1329  _distancetoseedcell_perthickperlayer_scint_FH_zminus,
1330  _distancetoseedcell_perthickperlayer_scint_BH_zminus,
1331  _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zminus,
1332  _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zminus,
1333  _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zminus,
1334  _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zminus,
1335  _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zminus,
1336  _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zminus,
1337  _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zminus,
1338  _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zminus,
1339  _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zminus,
1340  _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zminus,
1341  _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zminus,
1342  _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zminus,
1343  _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zminus,
1344  _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zminus,
1345  _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zminus,
1346  _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zminus,
1347  _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zminus,
1348  _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zminus,
1349  _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zminus,
1350  _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zminus,
1351  _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zminus,
1352  _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zminus,
1353  _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zminus,
1354  _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zminus,
1355  _distancebetseedandmaxcell_perthickperlayer_120_EE_zminus,
1356  _distancebetseedandmaxcell_perthickperlayer_120_FH_zminus,
1357  _distancebetseedandmaxcell_perthickperlayer_120_BH_zminus,
1358  _distancebetseedandmaxcell_perthickperlayer_200_EE_zminus,
1359  _distancebetseedandmaxcell_perthickperlayer_200_FH_zminus,
1360  _distancebetseedandmaxcell_perthickperlayer_200_BH_zminus,
1361  _distancebetseedandmaxcell_perthickperlayer_300_EE_zminus,
1362  _distancebetseedandmaxcell_perthickperlayer_300_FH_zminus,
1363  _distancebetseedandmaxcell_perthickperlayer_300_BH_zminus,
1364  _distancebetseedandmaxcell_perthickperlayer_scint_EE_zminus,
1365  _distancebetseedandmaxcell_perthickperlayer_scint_FH_zminus,
1366  _distancebetseedandmaxcell_perthickperlayer_scint_BH_zminus,
1367  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zminus,
1368  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zminus,
1369  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zminus,
1370  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zminus,
1371  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zminus,
1372  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zminus,
1373  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zminus,
1374  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zminus,
1375  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zminus,
1376  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zminus,
1377  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zminus,
1378  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zminus,
1379  loopSubFolders=False,
1380  purpose=PlotPurpose.Timing, page="CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zminus"
1381  ))
1382 
1383 # z+
1384 hgcalLayerClustersPlotter.append("CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zplus", [
1385  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1386  ], PlotFolder(
1387  _distancetomaxcell_perthickperlayer_120_EE_zplus,
1388  _distancetomaxcell_perthickperlayer_120_FH_zplus,
1389  _distancetomaxcell_perthickperlayer_120_BH_zplus,
1390  _distancetomaxcell_perthickperlayer_200_EE_zplus,
1391  _distancetomaxcell_perthickperlayer_200_FH_zplus,
1392  _distancetomaxcell_perthickperlayer_200_BH_zplus,
1393  _distancetomaxcell_perthickperlayer_300_EE_zplus,
1394  _distancetomaxcell_perthickperlayer_300_FH_zplus,
1395  _distancetomaxcell_perthickperlayer_300_BH_zplus,
1396  _distancetomaxcell_perthickperlayer_scint_EE_zplus,
1397  _distancetomaxcell_perthickperlayer_scint_FH_zplus,
1398  _distancetomaxcell_perthickperlayer_scint_BH_zplus,
1399  _distancetoseedcell_perthickperlayer_120_EE_zplus,
1400  _distancetoseedcell_perthickperlayer_120_FH_zplus,
1401  _distancetoseedcell_perthickperlayer_120_BH_zplus,
1402  _distancetoseedcell_perthickperlayer_200_EE_zplus,
1403  _distancetoseedcell_perthickperlayer_200_FH_zplus,
1404  _distancetoseedcell_perthickperlayer_200_BH_zplus,
1405  _distancetoseedcell_perthickperlayer_300_EE_zplus,
1406  _distancetoseedcell_perthickperlayer_300_FH_zplus,
1407  _distancetoseedcell_perthickperlayer_300_BH_zplus,
1408  _distancetoseedcell_perthickperlayer_scint_EE_zplus,
1409  _distancetoseedcell_perthickperlayer_scint_FH_zplus,
1410  _distancetoseedcell_perthickperlayer_scint_BH_zplus,
1411  _distancetomaxcell_perthickperlayer_eneweighted_120_EE_zplus,
1412  _distancetomaxcell_perthickperlayer_eneweighted_120_FH_zplus,
1413  _distancetomaxcell_perthickperlayer_eneweighted_120_BH_zplus,
1414  _distancetomaxcell_perthickperlayer_eneweighted_200_EE_zplus,
1415  _distancetomaxcell_perthickperlayer_eneweighted_200_FH_zplus,
1416  _distancetomaxcell_perthickperlayer_eneweighted_200_BH_zplus,
1417  _distancetomaxcell_perthickperlayer_eneweighted_300_EE_zplus,
1418  _distancetomaxcell_perthickperlayer_eneweighted_300_FH_zplus,
1419  _distancetomaxcell_perthickperlayer_eneweighted_300_BH_zplus,
1420  _distancetomaxcell_perthickperlayer_eneweighted_scint_EE_zplus,
1421  _distancetomaxcell_perthickperlayer_eneweighted_scint_FH_zplus,
1422  _distancetomaxcell_perthickperlayer_eneweighted_scint_BH_zplus,
1423  _distancetoseedcell_perthickperlayer_eneweighted_120_EE_zplus,
1424  _distancetoseedcell_perthickperlayer_eneweighted_120_FH_zplus,
1425  _distancetoseedcell_perthickperlayer_eneweighted_120_BH_zplus,
1426  _distancetoseedcell_perthickperlayer_eneweighted_200_EE_zplus,
1427  _distancetoseedcell_perthickperlayer_eneweighted_200_FH_zplus,
1428  _distancetoseedcell_perthickperlayer_eneweighted_200_BH_zplus,
1429  _distancetoseedcell_perthickperlayer_eneweighted_300_EE_zplus,
1430  _distancetoseedcell_perthickperlayer_eneweighted_300_FH_zplus,
1431  _distancetoseedcell_perthickperlayer_eneweighted_300_BH_zplus,
1432  _distancetoseedcell_perthickperlayer_eneweighted_scint_EE_zplus,
1433  _distancetoseedcell_perthickperlayer_eneweighted_scint_FH_zplus,
1434  _distancetoseedcell_perthickperlayer_eneweighted_scint_BH_zplus,
1435  _distancebetseedandmaxcell_perthickperlayer_120_EE_zplus,
1436  _distancebetseedandmaxcell_perthickperlayer_120_FH_zplus,
1437  _distancebetseedandmaxcell_perthickperlayer_120_BH_zplus,
1438  _distancebetseedandmaxcell_perthickperlayer_200_EE_zplus,
1439  _distancebetseedandmaxcell_perthickperlayer_200_FH_zplus,
1440  _distancebetseedandmaxcell_perthickperlayer_200_BH_zplus,
1441  _distancebetseedandmaxcell_perthickperlayer_300_EE_zplus,
1442  _distancebetseedandmaxcell_perthickperlayer_300_FH_zplus,
1443  _distancebetseedandmaxcell_perthickperlayer_300_BH_zplus,
1444  _distancebetseedandmaxcell_perthickperlayer_scint_EE_zplus,
1445  _distancebetseedandmaxcell_perthickperlayer_scint_FH_zplus,
1446  _distancebetseedandmaxcell_perthickperlayer_scint_BH_zplus,
1447  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_EE_zplus,
1448  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_FH_zplus,
1449  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_120_BH_zplus,
1450  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_EE_zplus,
1451  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_FH_zplus,
1452  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_200_BH_zplus,
1453  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_EE_zplus,
1454  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_FH_zplus,
1455  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_300_BH_zplus,
1456  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_EE_zplus,
1457  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_FH_zplus,
1458  _distancebetseedandmaxcellvsclusterenergy_perthickperlayer_scint_BH_zplus,
1459  loopSubFolders=False,
1460  purpose=PlotPurpose.Timing, page="CellsDistanceToSeedAndMaxCellPerLayerPerThickness_zplus"
1461  ))
1462 
1463 # [F] Looking at the fraction of true energy that has been clustered; by layer and overall
1464 # z-
1465 hgcalLayerClustersPlotter.append("EnergyClusteredByLayerAndOverall_zminus", [
1466  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1467  ], PlotFolder(
1468  _energyclustered_perlayer_EE_zminus,
1469  _energyclustered_perlayer_FH_zminus,
1470  _energyclustered_perlayer_BH_zminus,
1471  loopSubFolders=False,
1472  purpose=PlotPurpose.Timing, page="EnergyClusteredByLayerAndOverall_zminus"
1473  ))
1474 # z+
1475 hgcalLayerClustersPlotter.append("EnergyClusteredByLayerAndOverall_zplus", [
1476  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1477  ], PlotFolder(
1478  _energyclustered_perlayer_EE_zplus,
1479  _energyclustered_perlayer_FH_zplus,
1480  _energyclustered_perlayer_BH_zplus,
1481  loopSubFolders=False,
1482  purpose=PlotPurpose.Timing, page="EnergyClusteredByLayerAndOverall_zplus"
1483  ))
1484 
1485 # [G] Miscellaneous plots:
1486 # longdepthbarycentre: The longitudinal depth barycentre. One entry per event.
1487 # mixedhitscluster: Number of clusters per event with hits in different thicknesses.
1488 # num_reco_cluster_eta: Number of reco clusters vs eta
1489 
1490 hgcalLayerClustersPlotter.append("Miscellaneous", [
1491  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1492  ], PlotFolder(
1493  _num_reco_cluster_eta,
1494  _energyclustered,
1495  _mixedhitsclusters,
1496  _longdepthbarycentre,
1497  loopSubFolders=False,
1498  purpose=PlotPurpose.Timing, page="Miscellaneous"
1499  ))
1500 
1501 # [H] SelectedCaloParticles plots
1502 hgcalLayerClustersPlotter.append("SelectedCaloParticles_Photons", [
1503  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/SelectedCaloParticles/22",
1504  ], PlotFolder(
1505  _SelectedCaloParticles,
1506  loopSubFolders=False,
1507  purpose=PlotPurpose.Timing, page="SelectedCaloParticles_Photons"
1508  ))
1509 
1510 # [I] Score of CaloParticles wrt Layer Clusters
1511 # z-
1512 hgcalLayerClustersPlotter.append("ScoreCaloParticlesToLayerClusters_zminus", [
1513  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1514  ], PlotFolder(
1515  _score_caloparticle_to_layerclusters_zminus,
1516  loopSubFolders=False,
1517  purpose=PlotPurpose.Timing, page="ScoreCaloParticlesToLayerClusters_zminus"))
1518 
1519 # z+
1520 hgcalLayerClustersPlotter.append("ScoreCaloParticlesToLayerClusters_zplus", [
1521  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1522  ], PlotFolder(
1523  _score_caloparticle_to_layerclusters_zplus,
1524  loopSubFolders=False,
1525  purpose=PlotPurpose.Timing, page="ScoreCaloParticlesToLayerClusters_zplus"))
1526 
1527 # [J] Score of LayerClusters wrt CaloParticles
1528 # z-
1529 hgcalLayerClustersPlotter.append("ScoreLayerClustersToCaloParticles_zminus", [
1530  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1531  ], PlotFolder(
1532  _score_layercluster_to_caloparticles_zminus,
1533  loopSubFolders=False,
1534  purpose=PlotPurpose.Timing, page="ScoreLayerClustersToCaloParticles_zminus"))
1535 
1536 # z+
1537 hgcalLayerClustersPlotter.append("ScoreLayerClustersToCaloParticles_zplus", [
1538  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1539  ], PlotFolder(
1540  _score_layercluster_to_caloparticles_zplus,
1541  loopSubFolders=False,
1542  purpose=PlotPurpose.Timing, page="ScoreLayerClustersToCaloParticles_zplus"))
1543 
1544 # [K] Shared Energy between CaloParticle and LayerClusters
1545 # z-
1546 hgcalLayerClustersPlotter.append("SharedEnergy_zminus", [
1547  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1548  ], PlotFolder(
1549  _sharedEnergy_caloparticle_to_layercluster_zminus,
1550  loopSubFolders=False,
1551  purpose=PlotPurpose.Timing, page="SharedEnergyCaloParticleToLayerCluster_zminus"))
1552 
1553 # z+
1554 hgcalLayerClustersPlotter.append("SharedEnergy_zplus", [
1555  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1556  ], PlotFolder(
1557  _sharedEnergy_caloparticle_to_layercluster_zplus,
1558  loopSubFolders=False,
1559  purpose=PlotPurpose.Timing, page="SharedEnergyCaloParticleToLayerCluster_zplus"))
1560 
1561 # [K2] Shared Energy between LayerClusters and CaloParticle
1562 # z-
1563 hgcalLayerClustersPlotter.append("SharedEnergy_zminus", [
1564  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1565  ], PlotFolder(
1566  _sharedEnergy_layercluster_to_caloparticle_zminus,
1567  loopSubFolders=False,
1568  purpose=PlotPurpose.Timing, page="SharedEnergyLayerClusterToCaloParticle_zminus"))
1569 
1570 # z+
1571 hgcalLayerClustersPlotter.append("SharedEnergy_zplus", [
1572  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1573  ], PlotFolder(
1574  _sharedEnergy_layercluster_to_caloparticle_zplus,
1575  loopSubFolders=False,
1576  purpose=PlotPurpose.Timing, page="SharedEnergyLayerClusterToCaloParticle_zplus"))
1577 
1578 # [L] Cell Association per Layer
1579 # z-
1580 hgcalLayerClustersPlotter.append("CellAssociation_zminus", [
1581  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1582  ], PlotFolder(
1583  _cell_association_table_zminus,
1584  loopSubFolders=False,
1585  purpose=PlotPurpose.Timing, page="CellAssociation_zminus"))
1586 
1587 # z+
1588 hgcalLayerClustersPlotter.append("CellAssociation_zplus", [
1589  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1590  ], PlotFolder(
1591  _cell_association_table_zplus,
1592  loopSubFolders=False,
1593  purpose=PlotPurpose.Timing, page="CellAssociation_zplus"))
1594 
1595 # [M] Efficiency Plots
1596 # z-
1597 hgcalLayerClustersPlotter.append("Efficiencies_zminus", [
1598  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1599  ], PlotFolder(
1600  _efficiencies_zminus,
1601  loopSubFolders=False,
1602  purpose=PlotPurpose.Timing, page="Efficiencies_zminus"))
1603 
1604 # z+
1605 hgcalLayerClustersPlotter.append("Efficiencies_zplus", [
1606  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1607  ], PlotFolder(
1608  _efficiencies_zplus,
1609  loopSubFolders=False,
1610  purpose=PlotPurpose.Timing, page="Efficiencies_zplus"))
1611 
1612 # [L] Duplicate Plots
1613 # z-
1614 hgcalLayerClustersPlotter.append("Duplicates_zminus", [
1615  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1616  ], PlotFolder(
1617  _duplicates_zminus,
1618  loopSubFolders=False,
1619  purpose=PlotPurpose.Timing, page="Duplicates_zminus"))
1620 
1621 # z+
1622 hgcalLayerClustersPlotter.append("Duplicates_zplus", [
1623  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1624  ], PlotFolder(
1625  _duplicates_zplus,
1626  loopSubFolders=False,
1627  purpose=PlotPurpose.Timing, page="Duplicates_zplus"))
1628 
1629 # [M] Fake Rate Plots
1630 # z-
1631 hgcalLayerClustersPlotter.append("FakeRate_zminus", [
1632  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1633  ], PlotFolder(
1634  _fakes_zminus,
1635  loopSubFolders=False,
1636  purpose=PlotPurpose.Timing, page="Fakes_zminus"))
1637 
1638 # z+
1639 hgcalLayerClustersPlotter.append("FakeRate_zplus", [
1640  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1641  ], PlotFolder(
1642  _fakes_zplus,
1643  loopSubFolders=False,
1644  purpose=PlotPurpose.Timing, page="Fakes_zplus"))
1645 
1646 # [N] Merge Rate Plots
1647 # z-
1648 hgcalLayerClustersPlotter.append("MergeRate_zminus", [
1649  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1650  ], PlotFolder(
1651  _merges_zminus,
1652  loopSubFolders=False,
1653  purpose=PlotPurpose.Timing, page="Merges_zminus"))
1654 
1655 # z+
1656 hgcalLayerClustersPlotter.append("MergeRate_zplus", [
1657  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1658  ], PlotFolder(
1659  _merges_zplus,
1660  loopSubFolders=False,
1661  purpose=PlotPurpose.Timing, page="Merges_zplus"))
1662 
1663 # [O] Energy vs Score 2D plots CP to LC
1664 # z-
1665 for i,item in enumerate(_energyscore_cp2lc_zminus, start=1):
1666  hgcalLayerClustersPlotter.append("Energy_vs_Score_CP2LC_zminus", [
1667  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1668  ], PlotFolder(
1669  item,
1670  loopSubFolders=False,
1671  purpose=PlotPurpose.Timing, page="Energy_vs_Score_CP2LC_zminus"))
1672 
1673 # z+
1674 for i,item in enumerate(_energyscore_cp2lc_zplus, start=1):
1675  hgcalLayerClustersPlotter.append("Energy_vs_Score_CP2LC_zplus", [
1676  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1677  ], PlotFolder(
1678  item,
1679  loopSubFolders=False,
1680  purpose=PlotPurpose.Timing, page="Energy_vs_Score_CP2LC_zplus"))
1681 
1682 # [P] Energy vs Score 2D plots LC to CP
1683 # z-
1684 for i,item in enumerate(_energyscore_lc2cp_zminus, start=1):
1685  hgcalLayerClustersPlotter.append("Energy_vs_Score_LC2CP_zminus", [
1686  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1687  ], PlotFolder(
1688  item,
1689  loopSubFolders=False,
1690  purpose=PlotPurpose.Timing, page="Energy_vs_Score_LC2CP_zminus"))
1691 
1692 # z+
1693 for i,item in enumerate(_energyscore_lc2cp_zplus, start=1):
1694  hgcalLayerClustersPlotter.append("Energy_vs_Score_LC2CP_zplus", [
1695  "DQMData/Run 1/HGCAL/Run summary/HGCalValidator/hgcalLayerClusters",
1696  ], PlotFolder(
1697  item,
1698  loopSubFolders=False,
1699  purpose=PlotPurpose.Timing, page="Energy_vs_Score_LC2CP_zplus"))
1700 
edm::print
S & print(S &os, JobReport::InputFile const &f)
Definition: JobReport.cc:66
Generated for CMSSW Reference Manual by doxygen 1.8.11