58 t_ =
new TTree(
"Ntuple",
"Ntuple");
61 t_->Branch(
"subdetId", &
subdetId,
"subdetId/I", bufsize);
63 t_->Branch(
"layer", &
layer,
"layer/I", bufsize);
64 t_->Branch(
"ladder", &
ladder,
"ladder/I", bufsize);
65 t_->Branch(
"mod", &
mod,
"mod/I", bufsize);
66 t_->Branch(
"side", &
side,
"side/I", bufsize);
67 t_->Branch(
"disk", &
disk,
"disk/I", bufsize);
68 t_->Branch(
"blade", &
blade,
"blade/I", bufsize);
69 t_->Branch(
"panel", &
panel,
"panel/I", bufsize);
70 t_->Branch(
"plaq", &
plaq,
"plaq/I", bufsize);
72 t_->Branch(
"rechitx", &
rechitx,
"rechitx/F", bufsize);
73 t_->Branch(
"rechity", &
rechity,
"rechity/F", bufsize);
74 t_->Branch(
"rechitz", &
rechitz,
"rechitz/F", bufsize);
75 t_->Branch(
"rechiterrx", &
rechiterrx,
"rechiterrx/F", bufsize);
76 t_->Branch(
"rechiterry", &
rechiterry,
"rechiterry/F", bufsize);
77 t_->Branch(
"rechitresx", &
rechitresx,
"rechitresx/F", bufsize);
78 t_->Branch(
"rechitresy", &
rechitresy,
"rechitresy/F", bufsize);
79 t_->Branch(
"rechitpullx", &
rechitpullx,
"rechitpullx/F", bufsize);
80 t_->Branch(
"rechitpully", &
rechitpully,
"rechitpully/F", bufsize);
82 t_->Branch(
"npix", &
npix,
"npix/I", bufsize);
83 t_->Branch(
"nxpix", &
nxpix,
"nxpix/I", bufsize);
84 t_->Branch(
"nypix", &
nypix,
"nypix/I", bufsize);
85 t_->Branch(
"charge", &
charge,
"charge/F", bufsize);
87 t_->Branch(
"alpha", &
alpha,
"alpha/F", bufsize);
88 t_->Branch(
"beta", &
beta,
"beta/F", bufsize);
90 t_->Branch(
"phi", &
phi,
"phi/F", bufsize);
91 t_->Branch(
"eta", &
eta,
"eta/F", bufsize);
93 t_->Branch(
"half", &
half,
"half/I", bufsize);
94 t_->Branch(
"flipped", &
flipped,
"flipped/I", bufsize);
96 t_->Branch(
"simhitx", &
simhitx,
"simhitx/F", bufsize);
97 t_->Branch(
"simhity", &
simhity,
"simhity/F", bufsize);
99 t_->Branch(
"nsimhit", &
nsimhit,
"nsimhit/I", bufsize);
100 t_->Branch(
"pidhit", &
pidhit,
"pidhit/I", bufsize);
102 t_->Branch(
"evt", &
evt,
"evt/I", bufsize);
103 t_->Branch(
"run", &
run,
"run/I", bufsize);
110 trackerHitAssociatorConfig_(ps, consumesCollector()),
144 float pullxl = -10.0;
153 float pullyl = -10.0;
156 float barrel_alphal = 80.0;
157 float barrel_alphah = 100.0;
158 float barrel_betal = 10.0;
159 float barrel_betah = 170.0;
165 float forward_p1_alphal = 100.0;
166 float forward_p1_alphah = 115.0;
167 float forward_p2_alphal = 65.0;
168 float forward_p2_alphah = 80.0;
169 float forward_neg_betal = 67.0;
170 float forward_neg_betah = 73.0;
171 float forward_pos_betal = 107.0;
172 float forward_pos_betah = 113.0;
181 float pull_barrel_alphal = 80.0;
182 float pull_barrel_alphah = 100.0;
183 float pull_barrel_betal = 10.0;
184 float pull_barrel_betah = 170.0;
185 float pull_barrel_phil = -180.0;
186 float pull_barrel_phih = 180.0;
187 float pull_barrel_etal = -2.4;
188 float pull_barrel_etah = 2.4;
190 float pull_forward_p1_alphal = 100.0;
191 float pull_forward_p1_alphah = 112.0;
192 float pull_forward_p2_alphal = 68.0;
193 float pull_forward_p2_alphah = 80.0;
194 float pull_forward_neg_betal = 68.0;
195 float pull_forward_neg_betah = 72.0;
196 float pull_forward_pos_betal = 108.0;
197 float pull_forward_pos_betah = 112.0;
198 float pull_forward_phil = -180.0;
199 float pull_forward_phih = 180.0;
200 float pull_forward_neg_etal = -2.4;
201 float pull_forward_neg_etah = -1.4;
202 float pull_forward_pos_etal = 1.5;
203 float pull_forward_pos_etah = 2.5;
212 float barrel_chargel = 0.0;
213 float barrel_chargeh = 250000.0;
214 float forward_chargel = 0.0;
215 float forward_chargeh = 100000.0;
217 ibooker.
setCurrentFolder(
"Tracking/TrackingRecHits/Pixel/Histograms_per_ring-layer_or_disk-plaquette");
220 for (
int i = 0;
i < 3;
i++)
224 sprintf(chisto,
"meResxBarrelLayer_%d",
i + 1);
226 sprintf(chisto,
"meResyBarrelLayer_%d",
i + 1);
228 sprintf(chisto,
"mePullxBarrelLayer_%d",
i + 1);
230 sprintf(chisto,
"mePullyBarrelLayer_%d",
i + 1);
233 sprintf(chisto,
"meResXvsAlphaBarrelFlippedLaddersLayer_%d",
i + 1);
235 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
236 sprintf(chisto,
"meResYvsAlphaBarrelFlippedLaddersLayer_%d",
i + 1);
238 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
239 sprintf(chisto,
"meResXvsBetaBarrelFlippedLaddersLayer_%d",
i + 1);
241 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
242 sprintf(chisto,
"meResYvsBetaBarrelFlippedLaddersLayer_%d",
i + 1);
244 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
246 sprintf(chisto,
"meResXvsAlphaBarrelNonFlippedLaddersLayer_%d",
i + 1);
248 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
249 sprintf(chisto,
"meResYvsAlphaBarrelNonFlippedLaddersLayer_%d",
i + 1);
251 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
252 sprintf(chisto,
"meResXvsBetaBarrelNonFlippedLaddersLayer_%d",
i + 1);
254 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
255 sprintf(chisto,
"meResYvsBetaBarrelNonFlippedLaddersLayer_%d",
i + 1);
257 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
259 for (
int j = 0;
j < 8;
j++)
261 sprintf(chisto,
"mePosxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
263 sprintf(chisto,
"mePosyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
265 sprintf(chisto,
"meErrxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
267 sprintf(chisto,
"meErryBarrelLayerModule_%d_%d",
i + 1,
j + 1);
269 sprintf(chisto,
"meResxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
271 sprintf(chisto,
"meResyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
273 sprintf(chisto,
"mePullxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
275 sprintf(chisto,
"mePullyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
277 sprintf(chisto,
"meNpixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
279 sprintf(chisto,
"meNxpixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
281 sprintf(chisto,
"meNypixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
283 sprintf(chisto,
"meChargeBarrelLayerModule_%d_%d",
i + 1,
j + 1);
286 sprintf(chisto,
"meResXvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
288 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
289 sprintf(chisto,
"meResYvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
291 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
292 sprintf(chisto,
"meResXvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
294 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
295 sprintf(chisto,
"meResYvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
297 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
299 sprintf(chisto,
"mePullXvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
301 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
302 sprintf(chisto,
"mePullYvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
304 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
305 sprintf(chisto,
"mePullXvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
307 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
308 sprintf(chisto,
"mePullYvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
310 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
311 sprintf(chisto,
"mePullXvsPhiBarrelLayerModule_%d_%d",
i + 1,
j + 1);
313 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
314 sprintf(chisto,
"mePullYvsPhiBarrelLayerModule_%d_%d",
i + 1,
j + 1);
316 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
317 sprintf(chisto,
"mePullXvsEtaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
319 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
320 sprintf(chisto,
"mePullYvsEtaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
322 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
333 for (
int i = 0;
i < 2;
i++)
334 for (
int j = 0;
j < 4;
j++)
338 sprintf(chisto,
"mePosxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
340 sprintf(chisto,
"mePosyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
342 sprintf(chisto,
"meErrxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
344 sprintf(chisto,
"meErryZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
346 sprintf(chisto,
"meResxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
348 sprintf(chisto,
"meResyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
350 sprintf(chisto,
"mePullxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
352 sprintf(chisto,
"mePullyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
354 sprintf(chisto,
"meNpixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
356 sprintf(chisto,
"meNxpixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
358 sprintf(chisto,
"meNypixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
360 sprintf(chisto,
"meChargeZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
363 sprintf(chisto,
"meResXvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
365 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
366 sprintf(chisto,
"meResYvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
368 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
369 sprintf(chisto,
"meResXvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
371 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
372 sprintf(chisto,
"meResYvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
374 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
376 sprintf(chisto,
"mePullXvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
378 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
379 sprintf(chisto,
"mePullYvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
381 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
382 sprintf(chisto,
"mePullXvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
384 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
385 sprintf(chisto,
"mePullYvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
387 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
388 sprintf(chisto,
"mePullXvsPhiZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
390 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
391 sprintf(chisto,
"mePullYvsPhiZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
393 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
394 sprintf(chisto,
"mePullXvsEtaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
396 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
397 sprintf(chisto,
"mePullYvsEtaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
399 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
401 sprintf(chisto,
"mePosxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
403 sprintf(chisto,
"mePosyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
405 sprintf(chisto,
"meErrxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
407 sprintf(chisto,
"meErryZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
409 sprintf(chisto,
"meResxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
411 sprintf(chisto,
"meResyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
413 sprintf(chisto,
"mePullxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
415 sprintf(chisto,
"mePullyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
417 sprintf(chisto,
"meNpixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
419 sprintf(chisto,
"meNxpixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
421 sprintf(chisto,
"meNypixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
423 sprintf(chisto,
"meChargeZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
425 sprintf(chisto,
"meResXvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
427 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
428 sprintf(chisto,
"meResYvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
430 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
431 sprintf(chisto,
"meResXvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
433 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
434 sprintf(chisto,
"meResYvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
436 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
438 sprintf(chisto,
"mePullXvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
440 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
441 sprintf(chisto,
"mePullYvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
443 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
444 sprintf(chisto,
"mePullXvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
446 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
447 sprintf(chisto,
"mePullYvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
449 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
450 sprintf(chisto,
"mePullXvsPhiZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
452 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
453 sprintf(chisto,
"mePullYvsPhiZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
455 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
456 sprintf(chisto,
"mePullXvsEtaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
458 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
459 sprintf(chisto,
"mePullYvsEtaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
461 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
466 sprintf(chisto,
"mePosxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
468 sprintf(chisto,
"mePosyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
470 sprintf(chisto,
"meErrxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
472 sprintf(chisto,
"meErryZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
474 sprintf(chisto,
"meResxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
476 sprintf(chisto,
"meResyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
478 sprintf(chisto,
"mePullxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
480 sprintf(chisto,
"mePullyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
482 sprintf(chisto,
"meNpixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
484 sprintf(chisto,
"meNxpixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
486 sprintf(chisto,
"meNypixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
488 sprintf(chisto,
"meChargeZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
490 sprintf(chisto,
"meResXvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
492 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
493 sprintf(chisto,
"meResYvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
495 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
496 sprintf(chisto,
"meResXvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
498 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
499 sprintf(chisto,
"meResYvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
501 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
503 sprintf(chisto,
"mePullXvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
505 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
506 sprintf(chisto,
"mePullYvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
508 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
509 sprintf(chisto,
"mePullXvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
511 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
512 sprintf(chisto,
"mePullYvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
514 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
515 sprintf(chisto,
"mePullXvsPhiZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
517 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
518 sprintf(chisto,
"mePullYvsPhiZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
520 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
521 sprintf(chisto,
"mePullXvsEtaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
523 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
524 sprintf(chisto,
"mePullYvsEtaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
526 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
528 sprintf(chisto,
"mePosxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
530 sprintf(chisto,
"mePosyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
532 sprintf(chisto,
"meErrxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
534 sprintf(chisto,
"meErryZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
536 sprintf(chisto,
"meResxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
538 sprintf(chisto,
"meResyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
540 sprintf(chisto,
"mePullxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
542 sprintf(chisto,
"mePullyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
544 sprintf(chisto,
"meNpixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
546 sprintf(chisto,
"meNxpixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
548 sprintf(chisto,
"meNypixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
550 sprintf(chisto,
"meChargeZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
552 sprintf(chisto,
"meResXvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
554 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
555 sprintf(chisto,
"meResYvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
557 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
558 sprintf(chisto,
"meResXvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
560 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
561 sprintf(chisto,
"meResYvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
563 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
565 sprintf(chisto,
"mePullXvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
567 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
568 sprintf(chisto,
"mePullYvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
570 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
571 sprintf(chisto,
"mePullXvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
573 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
574 sprintf(chisto,
"mePullYvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
576 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
577 sprintf(chisto,
"mePullXvsPhiZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
579 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
580 sprintf(chisto,
"mePullYvsPhiZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
582 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
583 sprintf(chisto,
"mePullXvsEtaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
585 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
586 sprintf(chisto,
"mePullYvsEtaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
588 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
595 sprintf(chisto,
"mePosxBarrel");
597 sprintf(chisto,
"mePosyBarrel");
599 sprintf(chisto,
"meErrxBarrel");
601 sprintf(chisto,
"meErryBarrel");
603 sprintf(chisto,
"meResxBarrel");
605 sprintf(chisto,
"meResyBarrel");
607 sprintf(chisto,
"mePullxBarrel");
609 sprintf(chisto,
"mePullyBarrel");
611 sprintf(chisto,
"meNpixBarrel");
613 sprintf(chisto,
"meNxpixBarrel");
615 sprintf(chisto,
"meNypixBarrel");
617 sprintf(chisto,
"meChargeBarrel");
619 sprintf(chisto,
"meResXvsAlphaBarrel");
621 sprintf(chisto,
"meResYvsAlphaBarrel");
623 sprintf(chisto,
"meResXvsBetaBarrel");
625 sprintf(chisto,
"meResYvsBetaBarrel");
628 sprintf(chisto,
"mePullXvsAlphaBarrel");
630 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
631 sprintf(chisto,
"mePullYvsAlphaBarrel");
633 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
634 sprintf(chisto,
"mePullXvsBetaBarrel");
636 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
637 sprintf(chisto,
"mePullYvsBetaBarrel");
639 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
640 sprintf(chisto,
"mePullXvsPhiBarrel");
642 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
643 sprintf(chisto,
"mePullYvsPhiBarrel");
645 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
646 sprintf(chisto,
"mePullXvsEtaBarrel");
648 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
649 sprintf(chisto,
"mePullYvsEtaBarrel");
651 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
653 sprintf(chisto,
"mePosxBarrelHalfModule");
655 sprintf(chisto,
"mePosxBarrelFullModule");
657 sprintf(chisto,
"mePosxBarrelFlippedLadders");
659 sprintf(chisto,
"mePosxBarrelNonFlippedLadders");
661 sprintf(chisto,
"mePosyBarrelHalfModule");
663 sprintf(chisto,
"mePosyBarrelFullModule");
665 sprintf(chisto,
"mePosyBarrelFlippedLadders");
667 sprintf(chisto,
"mePosyBarrelNonFlippedLadders");
670 sprintf(chisto,
"meResXvsAlphaBarrelFlippedLadders");
672 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
673 sprintf(chisto,
"meResYvsAlphaBarrelFlippedLadders");
675 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
676 sprintf(chisto,
"meResXvsBetaBarrelFlippedLadders");
678 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
679 sprintf(chisto,
"meResYvsBetaBarrelFlippedLadders");
681 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
683 sprintf(chisto,
"mePullXvsAlphaBarrelFlippedLadders");
685 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
686 sprintf(chisto,
"mePullYvsAlphaBarrelFlippedLadders");
688 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
689 sprintf(chisto,
"mePullXvsBetaBarrelFlippedLadders");
691 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
692 sprintf(chisto,
"mePullYvsBetaBarrelFlippedLadders");
694 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
695 sprintf(chisto,
"mePullXvsPhiBarrelFlippedLadders");
697 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
698 sprintf(chisto,
"mePullYvsPhiBarrelFlippedLadders");
700 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
701 sprintf(chisto,
"mePullXvsEtaBarrelFlippedLadders");
703 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
704 sprintf(chisto,
"mePullYvsEtaBarrelFlippedLadders");
706 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
708 sprintf(chisto,
"meWPullXvsAlphaBarrelFlippedLadders");
710 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
711 sprintf(chisto,
"meWPullYvsAlphaBarrelFlippedLadders");
713 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
714 sprintf(chisto,
"meWPullXvsBetaBarrelFlippedLadders");
716 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
717 sprintf(chisto,
"meWPullYvsBetaBarrelFlippedLadders");
719 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
721 sprintf(chisto,
"meResXvsAlphaBarrelNonFlippedLadders");
723 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
724 sprintf(chisto,
"meResYvsAlphaBarrelNonFlippedLadders");
726 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
727 sprintf(chisto,
"meResXvsBetaBarrelNonFlippedLadders");
729 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
730 sprintf(chisto,
"meResYvsBetaBarrelNonFlippedLadders");
732 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
734 sprintf(chisto,
"mePullXvsAlphaBarrelNonFlippedLadders");
736 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
737 sprintf(chisto,
"mePullYvsAlphaBarrelNonFlippedLadders");
739 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
740 sprintf(chisto,
"mePullXvsBetaBarrelNonFlippedLadders");
742 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
743 sprintf(chisto,
"mePullYvsBetaBarrelNonFlippedLadders");
745 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
746 sprintf(chisto,
"mePullXvsPhiBarrelNonFlippedLadders");
748 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
749 sprintf(chisto,
"mePullYvsPhiBarrelNonFlippedLadders");
751 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
752 sprintf(chisto,
"mePullXvsEtaBarrelNonFlippedLadders");
754 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
755 sprintf(chisto,
"mePullYvsEtaBarrelNonFlippedLadders");
757 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
759 sprintf(chisto,
"meWPullXvsAlphaBarrelNonFlippedLadders");
761 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
762 sprintf(chisto,
"meWPullYvsAlphaBarrelNonFlippedLadders");
764 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
765 sprintf(chisto,
"meWPullXvsBetaBarrelNonFlippedLadders");
767 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
768 sprintf(chisto,
"meWPullYvsBetaBarrelNonFlippedLadders");
770 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
772 sprintf(chisto,
"mePosxZmPanel1");
774 sprintf(chisto,
"mePosyZmPanel1");
776 sprintf(chisto,
"meErrxZmPanel1");
778 sprintf(chisto,
"meErryZmPanel1");
780 sprintf(chisto,
"meResxZmPanel1");
782 sprintf(chisto,
"meResyZmPanel1");
784 sprintf(chisto,
"mePullxZmPanel1");
786 sprintf(chisto,
"mePullyZmPanel1");
788 sprintf(chisto,
"meNpixZmPanel1");
790 sprintf(chisto,
"meNxpixZmPanel1");
792 sprintf(chisto,
"meNypixZmPanel1");
794 sprintf(chisto,
"meChargeZmPanel1");
796 sprintf(chisto,
"meResXvsAlphaZmPanel1");
798 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
799 sprintf(chisto,
"meResYvsAlphaZmPanel1");
801 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
802 sprintf(chisto,
"meResXvsBetaZmPanel1");
804 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
805 sprintf(chisto,
"meResYvsBetaZmPanel1");
807 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
809 sprintf(chisto,
"mePullXvsAlphaZmPanel1");
811 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
812 sprintf(chisto,
"mePullYvsAlphaZmPanel1");
814 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
815 sprintf(chisto,
"mePullXvsBetaZmPanel1");
817 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
818 sprintf(chisto,
"mePullYvsBetaZmPanel1");
820 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
821 sprintf(chisto,
"mePullXvsPhiZmPanel1");
823 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
824 sprintf(chisto,
"mePullYvsPhiZmPanel1");
826 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
827 sprintf(chisto,
"mePullXvsEtaZmPanel1");
829 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
830 sprintf(chisto,
"mePullYvsEtaZmPanel1");
832 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
834 sprintf(chisto,
"meWPullXvsAlphaZmPanel1");
836 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
837 sprintf(chisto,
"meWPullYvsAlphaZmPanel1");
839 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
840 sprintf(chisto,
"meWPullXvsBetaZmPanel1");
842 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
843 sprintf(chisto,
"meWPullYvsBetaZmPanel1");
845 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
847 sprintf(chisto,
"mePosxZpPanel1");
849 sprintf(chisto,
"mePosyZpPanel1");
851 sprintf(chisto,
"meErrxZpPanel1");
853 sprintf(chisto,
"meErryZpPanel1");
855 sprintf(chisto,
"meResxZpPanel1");
857 sprintf(chisto,
"meResyZpPanel1");
859 sprintf(chisto,
"mePullxZpPanel1");
861 sprintf(chisto,
"mePullyZpPanel1");
863 sprintf(chisto,
"meNpixZpPanel1");
865 sprintf(chisto,
"meNxpixZpPanel1");
867 sprintf(chisto,
"meNypixZpPanel1");
869 sprintf(chisto,
"meChargeZpPanel1");
871 sprintf(chisto,
"meResXvsAlphaZpPanel1");
873 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
874 sprintf(chisto,
"meResYvsAlphaZpPanel1");
876 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
877 sprintf(chisto,
"meResXvsBetaZpPanel1");
879 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
880 sprintf(chisto,
"meResYvsBetaZpPanel1");
882 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
884 sprintf(chisto,
"mePullXvsAlphaZpPanel1");
886 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
887 sprintf(chisto,
"mePullYvsAlphaZpPanel1");
889 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
890 sprintf(chisto,
"mePullXvsBetaZpPanel1");
892 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
893 sprintf(chisto,
"mePullYvsBetaZpPanel1");
895 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
896 sprintf(chisto,
"mePullXvsPhiZpPanel1");
898 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
899 sprintf(chisto,
"mePullYvsPhiZpPanel1");
901 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
902 sprintf(chisto,
"mePullXvsEtaZpPanel1");
904 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
905 sprintf(chisto,
"mePullYvsEtaZpPanel1");
907 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
909 sprintf(chisto,
"meWPullXvsAlphaZpPanel1");
911 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
912 sprintf(chisto,
"meWPullYvsAlphaZpPanel1");
914 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
915 sprintf(chisto,
"meWPullXvsBetaZpPanel1");
917 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
918 sprintf(chisto,
"meWPullYvsBetaZpPanel1");
920 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
922 sprintf(chisto,
"mePosxZmPanel2");
924 sprintf(chisto,
"mePosyZmPanel2");
926 sprintf(chisto,
"meErrxZmPanel2");
928 sprintf(chisto,
"meErryZmPanel2");
930 sprintf(chisto,
"meResxZmPanel2");
932 sprintf(chisto,
"meResyZmPanel2");
934 sprintf(chisto,
"mePullxZmPanel2");
936 sprintf(chisto,
"mePullyZmPanel2");
938 sprintf(chisto,
"meNpixZmPanel2");
940 sprintf(chisto,
"meNxpixZmPanel2");
942 sprintf(chisto,
"meNypixZmPanel2");
944 sprintf(chisto,
"meChargeZmPanel2");
946 sprintf(chisto,
"meResXvsAlphaZmPanel2");
948 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
949 sprintf(chisto,
"meResYvsAlphaZmPanel2");
951 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
952 sprintf(chisto,
"meResXvsBetaZmPanel2");
954 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
955 sprintf(chisto,
"meResYvsBetaZmPanel2");
957 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
959 sprintf(chisto,
"mePullXvsAlphaZmPanel2");
961 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
962 sprintf(chisto,
"mePullYvsAlphaZmPanel2");
964 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
965 sprintf(chisto,
"mePullXvsBetaZmPanel2");
967 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
968 sprintf(chisto,
"mePullYvsBetaZmPanel2");
970 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
971 sprintf(chisto,
"mePullXvsPhiZmPanel2");
973 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
974 sprintf(chisto,
"mePullYvsPhiZmPanel2");
976 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
977 sprintf(chisto,
"mePullXvsEtaZmPanel2");
979 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
980 sprintf(chisto,
"mePullYvsEtaZmPanel2");
982 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
984 sprintf(chisto,
"meWPullXvsAlphaZmPanel2");
986 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
987 sprintf(chisto,
"meWPullYvsAlphaZmPanel2");
989 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
990 sprintf(chisto,
"meWPullXvsBetaZmPanel2");
992 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
993 sprintf(chisto,
"meWPullYvsBetaZmPanel2");
995 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
997 sprintf(chisto,
"mePosxZpPanel2");
999 sprintf(chisto,
"mePosyZpPanel2");
1001 sprintf(chisto,
"meErrxZpPanel2");
1003 sprintf(chisto,
"meErryZpPanel2");
1005 sprintf(chisto,
"meResxZpPanel2");
1007 sprintf(chisto,
"meResyZpPanel2");
1009 sprintf(chisto,
"mePullxZpPanel2");
1011 sprintf(chisto,
"mePullyZpPanel2");
1013 sprintf(chisto,
"meNpixZpPanel2");
1015 sprintf(chisto,
"meNxpixZpPanel2");
1017 sprintf(chisto,
"meNypixZpPanel2");
1019 sprintf(chisto,
"meChargeZpPanel2");
1021 sprintf(chisto,
"meResXvsAlphaZpPanel2");
1023 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
1024 sprintf(chisto,
"meResYvsAlphaZpPanel2");
1026 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
1027 sprintf(chisto,
"meResXvsBetaZpPanel2");
1029 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
1030 sprintf(chisto,
"meResYvsBetaZpPanel2");
1032 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
1034 sprintf(chisto,
"mePullXvsAlphaZpPanel2");
1036 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
1037 sprintf(chisto,
"mePullYvsAlphaZpPanel2");
1039 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
1040 sprintf(chisto,
"mePullXvsBetaZpPanel2");
1042 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
1043 sprintf(chisto,
"mePullYvsBetaZpPanel2");
1045 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
1046 sprintf(chisto,
"mePullXvsPhiZpPanel2");
1048 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
1049 sprintf(chisto,
"mePullYvsPhiZpPanel2");
1051 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
1052 sprintf(chisto,
"mePullXvsEtaZpPanel2");
1054 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
1055 sprintf(chisto,
"mePullYvsEtaZpPanel2");
1057 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
1059 sprintf(chisto,
"meWPullXvsAlphaZpPanel2");
1061 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
1062 sprintf(chisto,
"meWPullYvsAlphaZpPanel2");
1064 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
1065 sprintf(chisto,
"meWPullXvsBetaZpPanel2");
1067 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
1068 sprintf(chisto,
"meWPullYvsBetaZpPanel2");
1070 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
1073 sprintf(chisto,
"mePosxBarrel_all_hits");
1075 sprintf(chisto,
"mePosyBarrel_all_hits");
1078 sprintf(chisto,
"mePosxZmPanel1_all_hits");
1080 sprintf(chisto,
"mePosyZmPanel1_all_hits");
1082 sprintf(chisto,
"mePosxZmPanel2_all_hits");
1084 sprintf(chisto,
"mePosyZmPanel2_all_hits");
1087 sprintf(chisto,
"mePosxZpPanel1_all_hits");
1089 sprintf(chisto,
"mePosyZpPanel1_all_hits");
1091 sprintf(chisto,
"mePosxZpPanel2_all_hits");
1093 sprintf(chisto,
"mePosyZpPanel2_all_hits");
1114 evt =
e.id().event();
1120 float radtodeg = 180.0 /
math_pi;
1126 float mindist = 999999.9;
1144 for (
const auto& it :
tracker->dets()) {
1145 DetId detId = it->geographicalId();
1147 unsigned int subid = detId.
subdetId();
1148 if (!((subid == 1) || (subid == 2)))
1162 for (; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
1167 detId = it->geographicalId();
1173 side = tTopo->pxfSide(detId);
1174 disk = tTopo->pxfDisk(detId);
1175 blade = tTopo->pxfBlade(detId);
1176 panel = tTopo->pxfPanel(detId);
1177 plaq = tTopo->pxfModule(detId);
1183 }
else if (
panel == 2) {
1188 <<
"..............................................Wrong panel number !";
1190 else if (
side == 2) {
1194 }
else if (
panel == 2) {
1199 <<
"..............................................Wrong panel number !";
1203 <<
".......................................................Wrong side !";
1208 <<
"Pixel rechit collection but we are not in the pixel detector" << (
int)detId.
subdetId();
1217 reco::TrackCollection::const_iterator tciter;
1224 for (tciter =
tracks->begin(); tciter !=
tracks->end(); tciter++) {
1225 phi = tciter->momentum().phi() /
math_pi * 180.0;
1226 eta = tciter->momentum().eta();
1233 const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
1274 error = (*it)->localPositionError();
1282 npix = (*matchedhit).cluster()->size();
1283 nxpix = (*matchedhit).cluster()->sizeX();
1284 nypix = (*matchedhit).cluster()->sizeY();
1285 charge = (*matchedhit).cluster()->charge();
1295 float distx, disty, dist;
1296 bool found_hit_from_generated_particle =
false;
1298 int n_assoc_muon = 0;
1300 std::vector<PSimHit>::const_iterator closestit =
matched.begin();
1301 for (std::vector<PSimHit>::const_iterator
m =
matched.begin();
m <
matched.end();
m++) {
1303 int pid = (*m).particleType();
1308 float simhitx = 0.5 * ((*m).entryPoint().x() + (*m).exitPoint().x());
1309 float simhity = 0.5 * ((*m).entryPoint().y() + (*m).exitPoint().y());
1313 dist =
sqrt(distx * distx + disty * disty);
1315 if (dist < mindist) {
1320 found_hit_from_generated_particle =
true;
1326 if (
checkType_ && !found_hit_from_generated_particle)
1329 if (n_assoc_muon > 1) {
1331 <<
" ----- This is not good: n_assoc_muon = " << n_assoc_muon;
1335 pidhit = (*closestit).particleType();
1337 simhitx = 0.5 * ((*closestit).entryPoint().x() + (*closestit).exitPoint().x());
1338 simhity = 0.5 * ((*closestit).entryPoint().y() + (*closestit).exitPoint().y());
1345 float simhitpx = (*closestit).momentumAtEntry().x();
1346 float simhitpy = (*closestit).momentumAtEntry().y();
1347 float simhitpz = (*closestit).momentumAtEntry().z();
1352 beta = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
1353 alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;
1355 detId = (*it)->geographicalId();
1390 if (tmp_nrows == 80) {
1394 }
else if (tmp_nrows == 160) {
1400 <<
"-------------------------------------------------- Wrong module size !!!";
1451 layer = tTopo->pxbLayer(detId);
1452 ladder = tTopo->pxbLadder(detId);
1453 mod = tTopo->pxbModule(detId);
1498 side = tTopo->pxfSide(detId);
1499 disk = tTopo->pxfDisk(detId);
1500 blade = tTopo->pxfBlade(detId);
1501 panel = tTopo->pxfPanel(detId);
1502 plaq = tTopo->pxfModule(detId);
1561 }
else if (
panel == 2) {
1619 <<
"..............................................Wrong panel number !";
1621 else if (
side == 2) {
1678 }
else if (
panel == 2) {
1736 <<
"..............................................Wrong panel number !";
1740 <<
".......................................................Wrong side !";
1745 <<
"Pixel rechit but we are not in the pixel detector" << (
int)detId.
subdetId();