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);
108 : trackerHitAssociatorConfig_(ps, consumesCollector()),
dbe_(nullptr), tfile_(nullptr), t_(nullptr) {
139 float pullxl = -10.0;
148 float pullyl = -10.0;
151 float barrel_alphal = 80.0;
152 float barrel_alphah = 100.0;
153 float barrel_betal = 10.0;
154 float barrel_betah = 170.0;
160 float forward_p1_alphal = 100.0;
161 float forward_p1_alphah = 115.0;
162 float forward_p2_alphal = 65.0;
163 float forward_p2_alphah = 80.0;
164 float forward_neg_betal = 67.0;
165 float forward_neg_betah = 73.0;
166 float forward_pos_betal = 107.0;
167 float forward_pos_betah = 113.0;
176 float pull_barrel_alphal = 80.0;
177 float pull_barrel_alphah = 100.0;
178 float pull_barrel_betal = 10.0;
179 float pull_barrel_betah = 170.0;
180 float pull_barrel_phil = -180.0;
181 float pull_barrel_phih = 180.0;
182 float pull_barrel_etal = -2.4;
183 float pull_barrel_etah = 2.4;
185 float pull_forward_p1_alphal = 100.0;
186 float pull_forward_p1_alphah = 112.0;
187 float pull_forward_p2_alphal = 68.0;
188 float pull_forward_p2_alphah = 80.0;
189 float pull_forward_neg_betal = 68.0;
190 float pull_forward_neg_betah = 72.0;
191 float pull_forward_pos_betal = 108.0;
192 float pull_forward_pos_betah = 112.0;
193 float pull_forward_phil = -180.0;
194 float pull_forward_phih = 180.0;
195 float pull_forward_neg_etal = -2.4;
196 float pull_forward_neg_etah = -1.4;
197 float pull_forward_pos_etal = 1.5;
198 float pull_forward_pos_etah = 2.5;
207 float barrel_chargel = 0.0;
208 float barrel_chargeh = 250000.0;
209 float forward_chargel = 0.0;
210 float forward_chargeh = 100000.0;
212 ibooker.
setCurrentFolder(
"Tracking/TrackingRecHits/Pixel/Histograms_per_ring-layer_or_disk-plaquette");
215 for (
int i = 0;
i < 3;
i++)
219 sprintf(chisto,
"meResxBarrelLayer_%d",
i + 1);
221 sprintf(chisto,
"meResyBarrelLayer_%d",
i + 1);
223 sprintf(chisto,
"mePullxBarrelLayer_%d",
i + 1);
225 sprintf(chisto,
"mePullyBarrelLayer_%d",
i + 1);
228 sprintf(chisto,
"meResXvsAlphaBarrelFlippedLaddersLayer_%d",
i + 1);
230 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
231 sprintf(chisto,
"meResYvsAlphaBarrelFlippedLaddersLayer_%d",
i + 1);
233 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
234 sprintf(chisto,
"meResXvsBetaBarrelFlippedLaddersLayer_%d",
i + 1);
236 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
237 sprintf(chisto,
"meResYvsBetaBarrelFlippedLaddersLayer_%d",
i + 1);
239 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
241 sprintf(chisto,
"meResXvsAlphaBarrelNonFlippedLaddersLayer_%d",
i + 1);
243 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
244 sprintf(chisto,
"meResYvsAlphaBarrelNonFlippedLaddersLayer_%d",
i + 1);
246 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
247 sprintf(chisto,
"meResXvsBetaBarrelNonFlippedLaddersLayer_%d",
i + 1);
249 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
250 sprintf(chisto,
"meResYvsBetaBarrelNonFlippedLaddersLayer_%d",
i + 1);
252 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
254 for (
int j = 0;
j < 8;
j++)
256 sprintf(chisto,
"mePosxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
258 sprintf(chisto,
"mePosyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
260 sprintf(chisto,
"meErrxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
262 sprintf(chisto,
"meErryBarrelLayerModule_%d_%d",
i + 1,
j + 1);
264 sprintf(chisto,
"meResxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
266 sprintf(chisto,
"meResyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
268 sprintf(chisto,
"mePullxBarrelLayerModule_%d_%d",
i + 1,
j + 1);
270 sprintf(chisto,
"mePullyBarrelLayerModule_%d_%d",
i + 1,
j + 1);
272 sprintf(chisto,
"meNpixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
274 sprintf(chisto,
"meNxpixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
276 sprintf(chisto,
"meNypixBarrelLayerModule_%d_%d",
i + 1,
j + 1);
278 sprintf(chisto,
"meChargeBarrelLayerModule_%d_%d",
i + 1,
j + 1);
281 sprintf(chisto,
"meResXvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
283 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
284 sprintf(chisto,
"meResYvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
286 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
287 sprintf(chisto,
"meResXvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
289 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
290 sprintf(chisto,
"meResYvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
292 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
294 sprintf(chisto,
"mePullXvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
296 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
297 sprintf(chisto,
"mePullYvsAlphaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
299 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
300 sprintf(chisto,
"mePullXvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
302 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
303 sprintf(chisto,
"mePullYvsBetaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
305 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
306 sprintf(chisto,
"mePullXvsPhiBarrelLayerModule_%d_%d",
i + 1,
j + 1);
308 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
309 sprintf(chisto,
"mePullYvsPhiBarrelLayerModule_%d_%d",
i + 1,
j + 1);
311 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
312 sprintf(chisto,
"mePullXvsEtaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
314 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
315 sprintf(chisto,
"mePullYvsEtaBarrelLayerModule_%d_%d",
i + 1,
j + 1);
317 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
328 for (
int i = 0;
i < 2;
i++)
329 for (
int j = 0;
j < 4;
j++)
333 sprintf(chisto,
"mePosxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
335 sprintf(chisto,
"mePosyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
337 sprintf(chisto,
"meErrxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
339 sprintf(chisto,
"meErryZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
341 sprintf(chisto,
"meResxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
343 sprintf(chisto,
"meResyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
345 sprintf(chisto,
"mePullxZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
347 sprintf(chisto,
"mePullyZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
349 sprintf(chisto,
"meNpixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
351 sprintf(chisto,
"meNxpixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
353 sprintf(chisto,
"meNypixZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
355 sprintf(chisto,
"meChargeZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
358 sprintf(chisto,
"meResXvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
360 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
361 sprintf(chisto,
"meResYvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
363 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
364 sprintf(chisto,
"meResXvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
366 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
367 sprintf(chisto,
"meResYvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
369 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
371 sprintf(chisto,
"mePullXvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
373 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
374 sprintf(chisto,
"mePullYvsAlphaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
376 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
377 sprintf(chisto,
"mePullXvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
379 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
380 sprintf(chisto,
"mePullYvsBetaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
382 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
383 sprintf(chisto,
"mePullXvsPhiZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
385 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
386 sprintf(chisto,
"mePullYvsPhiZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
388 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
389 sprintf(chisto,
"mePullXvsEtaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
391 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
392 sprintf(chisto,
"mePullYvsEtaZmPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
394 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
396 sprintf(chisto,
"mePosxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
398 sprintf(chisto,
"mePosyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
400 sprintf(chisto,
"meErrxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
402 sprintf(chisto,
"meErryZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
404 sprintf(chisto,
"meResxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
406 sprintf(chisto,
"meResyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
408 sprintf(chisto,
"mePullxZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
410 sprintf(chisto,
"mePullyZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
412 sprintf(chisto,
"meNpixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
414 sprintf(chisto,
"meNxpixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
416 sprintf(chisto,
"meNypixZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
418 sprintf(chisto,
"meChargeZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
420 sprintf(chisto,
"meResXvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
422 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
423 sprintf(chisto,
"meResYvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
425 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
426 sprintf(chisto,
"meResXvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
428 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
429 sprintf(chisto,
"meResYvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
431 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
433 sprintf(chisto,
"mePullXvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
435 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
436 sprintf(chisto,
"mePullYvsAlphaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
438 chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
439 sprintf(chisto,
"mePullXvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
441 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
442 sprintf(chisto,
"mePullYvsBetaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
444 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
445 sprintf(chisto,
"mePullXvsPhiZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
447 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
448 sprintf(chisto,
"mePullYvsPhiZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
450 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
451 sprintf(chisto,
"mePullXvsEtaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
453 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
454 sprintf(chisto,
"mePullYvsEtaZpPanel1DiskPlaq_%d_%d",
i + 1,
j + 1);
456 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
461 sprintf(chisto,
"mePosxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
463 sprintf(chisto,
"mePosyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
465 sprintf(chisto,
"meErrxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
467 sprintf(chisto,
"meErryZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
469 sprintf(chisto,
"meResxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
471 sprintf(chisto,
"meResyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
473 sprintf(chisto,
"mePullxZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
475 sprintf(chisto,
"mePullyZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
477 sprintf(chisto,
"meNpixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
479 sprintf(chisto,
"meNxpixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
481 sprintf(chisto,
"meNypixZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
483 sprintf(chisto,
"meChargeZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
485 sprintf(chisto,
"meResXvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
487 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
488 sprintf(chisto,
"meResYvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
490 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
491 sprintf(chisto,
"meResXvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
493 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
494 sprintf(chisto,
"meResYvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
496 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
498 sprintf(chisto,
"mePullXvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
500 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
501 sprintf(chisto,
"mePullYvsAlphaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
503 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
504 sprintf(chisto,
"mePullXvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
506 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
507 sprintf(chisto,
"mePullYvsBetaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
509 chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
510 sprintf(chisto,
"mePullXvsPhiZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
512 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
513 sprintf(chisto,
"mePullYvsPhiZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
515 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
516 sprintf(chisto,
"mePullXvsEtaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
518 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
519 sprintf(chisto,
"mePullYvsEtaZmPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
521 chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
523 sprintf(chisto,
"mePosxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
525 sprintf(chisto,
"mePosyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
527 sprintf(chisto,
"meErrxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
529 sprintf(chisto,
"meErryZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
531 sprintf(chisto,
"meResxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
533 sprintf(chisto,
"meResyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
535 sprintf(chisto,
"mePullxZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
537 sprintf(chisto,
"mePullyZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
539 sprintf(chisto,
"meNpixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
541 sprintf(chisto,
"meNxpixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
543 sprintf(chisto,
"meNypixZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
545 sprintf(chisto,
"meChargeZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
547 sprintf(chisto,
"meResXvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
549 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
550 sprintf(chisto,
"meResYvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
552 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
553 sprintf(chisto,
"meResXvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
555 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
556 sprintf(chisto,
"meResYvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
558 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
560 sprintf(chisto,
"mePullXvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
562 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
563 sprintf(chisto,
"mePullYvsAlphaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
565 chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
566 sprintf(chisto,
"mePullXvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
568 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
569 sprintf(chisto,
"mePullYvsBetaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
571 chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
572 sprintf(chisto,
"mePullXvsPhiZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
574 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
575 sprintf(chisto,
"mePullYvsPhiZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
577 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
578 sprintf(chisto,
"mePullXvsEtaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
580 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
581 sprintf(chisto,
"mePullYvsEtaZpPanel2DiskPlaq_%d_%d",
i + 1,
j + 1);
583 chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
590 sprintf(chisto,
"mePosxBarrel");
592 sprintf(chisto,
"mePosyBarrel");
594 sprintf(chisto,
"meErrxBarrel");
596 sprintf(chisto,
"meErryBarrel");
598 sprintf(chisto,
"meResxBarrel");
600 sprintf(chisto,
"meResyBarrel");
602 sprintf(chisto,
"mePullxBarrel");
604 sprintf(chisto,
"mePullyBarrel");
606 sprintf(chisto,
"meNpixBarrel");
608 sprintf(chisto,
"meNxpixBarrel");
610 sprintf(chisto,
"meNypixBarrel");
612 sprintf(chisto,
"meChargeBarrel");
614 sprintf(chisto,
"meResXvsAlphaBarrel");
616 sprintf(chisto,
"meResYvsAlphaBarrel");
618 sprintf(chisto,
"meResXvsBetaBarrel");
620 sprintf(chisto,
"meResYvsBetaBarrel");
623 sprintf(chisto,
"mePullXvsAlphaBarrel");
625 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
626 sprintf(chisto,
"mePullYvsAlphaBarrel");
628 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
629 sprintf(chisto,
"mePullXvsBetaBarrel");
631 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
632 sprintf(chisto,
"mePullYvsBetaBarrel");
634 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
635 sprintf(chisto,
"mePullXvsPhiBarrel");
637 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
638 sprintf(chisto,
"mePullYvsPhiBarrel");
640 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
641 sprintf(chisto,
"mePullXvsEtaBarrel");
643 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
644 sprintf(chisto,
"mePullYvsEtaBarrel");
646 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
648 sprintf(chisto,
"mePosxBarrelHalfModule");
650 sprintf(chisto,
"mePosxBarrelFullModule");
652 sprintf(chisto,
"mePosxBarrelFlippedLadders");
654 sprintf(chisto,
"mePosxBarrelNonFlippedLadders");
656 sprintf(chisto,
"mePosyBarrelHalfModule");
658 sprintf(chisto,
"mePosyBarrelFullModule");
660 sprintf(chisto,
"mePosyBarrelFlippedLadders");
662 sprintf(chisto,
"mePosyBarrelNonFlippedLadders");
665 sprintf(chisto,
"meResXvsAlphaBarrelFlippedLadders");
667 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
668 sprintf(chisto,
"meResYvsAlphaBarrelFlippedLadders");
670 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
671 sprintf(chisto,
"meResXvsBetaBarrelFlippedLadders");
673 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
674 sprintf(chisto,
"meResYvsBetaBarrelFlippedLadders");
676 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
678 sprintf(chisto,
"mePullXvsAlphaBarrelFlippedLadders");
680 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
681 sprintf(chisto,
"mePullYvsAlphaBarrelFlippedLadders");
683 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
684 sprintf(chisto,
"mePullXvsBetaBarrelFlippedLadders");
686 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
687 sprintf(chisto,
"mePullYvsBetaBarrelFlippedLadders");
689 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
690 sprintf(chisto,
"mePullXvsPhiBarrelFlippedLadders");
692 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
693 sprintf(chisto,
"mePullYvsPhiBarrelFlippedLadders");
695 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
696 sprintf(chisto,
"mePullXvsEtaBarrelFlippedLadders");
698 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
699 sprintf(chisto,
"mePullYvsEtaBarrelFlippedLadders");
701 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
703 sprintf(chisto,
"meWPullXvsAlphaBarrelFlippedLadders");
705 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
706 sprintf(chisto,
"meWPullYvsAlphaBarrelFlippedLadders");
708 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
709 sprintf(chisto,
"meWPullXvsBetaBarrelFlippedLadders");
711 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
712 sprintf(chisto,
"meWPullYvsBetaBarrelFlippedLadders");
714 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
716 sprintf(chisto,
"meResXvsAlphaBarrelNonFlippedLadders");
718 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh,
"");
719 sprintf(chisto,
"meResYvsAlphaBarrelNonFlippedLadders");
721 ibooker.
bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh,
"");
722 sprintf(chisto,
"meResXvsBetaBarrelNonFlippedLadders");
724 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh,
"");
725 sprintf(chisto,
"meResYvsBetaBarrelNonFlippedLadders");
727 ibooker.
bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh,
"");
729 sprintf(chisto,
"mePullXvsAlphaBarrelNonFlippedLadders");
731 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
732 sprintf(chisto,
"mePullYvsAlphaBarrelNonFlippedLadders");
734 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
735 sprintf(chisto,
"mePullXvsBetaBarrelNonFlippedLadders");
737 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
738 sprintf(chisto,
"mePullYvsBetaBarrelNonFlippedLadders");
740 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
741 sprintf(chisto,
"mePullXvsPhiBarrelNonFlippedLadders");
743 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh,
"");
744 sprintf(chisto,
"mePullYvsPhiBarrelNonFlippedLadders");
746 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh,
"");
747 sprintf(chisto,
"mePullXvsEtaBarrelNonFlippedLadders");
749 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh,
"");
750 sprintf(chisto,
"mePullYvsEtaBarrelNonFlippedLadders");
752 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh,
"");
754 sprintf(chisto,
"meWPullXvsAlphaBarrelNonFlippedLadders");
756 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh,
"");
757 sprintf(chisto,
"meWPullYvsAlphaBarrelNonFlippedLadders");
759 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh,
"");
760 sprintf(chisto,
"meWPullXvsBetaBarrelNonFlippedLadders");
762 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh,
"");
763 sprintf(chisto,
"meWPullYvsBetaBarrelNonFlippedLadders");
765 ibooker.
bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh,
"");
767 sprintf(chisto,
"mePosxZmPanel1");
769 sprintf(chisto,
"mePosyZmPanel1");
771 sprintf(chisto,
"meErrxZmPanel1");
773 sprintf(chisto,
"meErryZmPanel1");
775 sprintf(chisto,
"meResxZmPanel1");
777 sprintf(chisto,
"meResyZmPanel1");
779 sprintf(chisto,
"mePullxZmPanel1");
781 sprintf(chisto,
"mePullyZmPanel1");
783 sprintf(chisto,
"meNpixZmPanel1");
785 sprintf(chisto,
"meNxpixZmPanel1");
787 sprintf(chisto,
"meNypixZmPanel1");
789 sprintf(chisto,
"meChargeZmPanel1");
791 sprintf(chisto,
"meResXvsAlphaZmPanel1");
793 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
794 sprintf(chisto,
"meResYvsAlphaZmPanel1");
796 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
797 sprintf(chisto,
"meResXvsBetaZmPanel1");
799 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
800 sprintf(chisto,
"meResYvsBetaZmPanel1");
802 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
804 sprintf(chisto,
"mePullXvsAlphaZmPanel1");
806 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
807 sprintf(chisto,
"mePullYvsAlphaZmPanel1");
809 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
810 sprintf(chisto,
"mePullXvsBetaZmPanel1");
812 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
813 sprintf(chisto,
"mePullYvsBetaZmPanel1");
815 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
816 sprintf(chisto,
"mePullXvsPhiZmPanel1");
818 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
819 sprintf(chisto,
"mePullYvsPhiZmPanel1");
821 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
822 sprintf(chisto,
"mePullXvsEtaZmPanel1");
824 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
825 sprintf(chisto,
"mePullYvsEtaZmPanel1");
827 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
829 sprintf(chisto,
"meWPullXvsAlphaZmPanel1");
831 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
832 sprintf(chisto,
"meWPullYvsAlphaZmPanel1");
834 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
835 sprintf(chisto,
"meWPullXvsBetaZmPanel1");
837 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
838 sprintf(chisto,
"meWPullYvsBetaZmPanel1");
840 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
842 sprintf(chisto,
"mePosxZpPanel1");
844 sprintf(chisto,
"mePosyZpPanel1");
846 sprintf(chisto,
"meErrxZpPanel1");
848 sprintf(chisto,
"meErryZpPanel1");
850 sprintf(chisto,
"meResxZpPanel1");
852 sprintf(chisto,
"meResyZpPanel1");
854 sprintf(chisto,
"mePullxZpPanel1");
856 sprintf(chisto,
"mePullyZpPanel1");
858 sprintf(chisto,
"meNpixZpPanel1");
860 sprintf(chisto,
"meNxpixZpPanel1");
862 sprintf(chisto,
"meNypixZpPanel1");
864 sprintf(chisto,
"meChargeZpPanel1");
866 sprintf(chisto,
"meResXvsAlphaZpPanel1");
868 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh,
"");
869 sprintf(chisto,
"meResYvsAlphaZpPanel1");
871 ibooker.
bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh,
"");
872 sprintf(chisto,
"meResXvsBetaZpPanel1");
874 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
875 sprintf(chisto,
"meResYvsBetaZpPanel1");
877 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
879 sprintf(chisto,
"mePullXvsAlphaZpPanel1");
881 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
882 sprintf(chisto,
"mePullYvsAlphaZpPanel1");
884 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
885 sprintf(chisto,
"mePullXvsBetaZpPanel1");
887 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
888 sprintf(chisto,
"mePullYvsBetaZpPanel1");
890 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
891 sprintf(chisto,
"mePullXvsPhiZpPanel1");
893 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
894 sprintf(chisto,
"mePullYvsPhiZpPanel1");
896 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
897 sprintf(chisto,
"mePullXvsEtaZpPanel1");
899 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
900 sprintf(chisto,
"mePullYvsEtaZpPanel1");
902 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
904 sprintf(chisto,
"meWPullXvsAlphaZpPanel1");
906 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh,
"");
907 sprintf(chisto,
"meWPullYvsAlphaZpPanel1");
909 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh,
"");
910 sprintf(chisto,
"meWPullXvsBetaZpPanel1");
912 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
913 sprintf(chisto,
"meWPullYvsBetaZpPanel1");
915 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
917 sprintf(chisto,
"mePosxZmPanel2");
919 sprintf(chisto,
"mePosyZmPanel2");
921 sprintf(chisto,
"meErrxZmPanel2");
923 sprintf(chisto,
"meErryZmPanel2");
925 sprintf(chisto,
"meResxZmPanel2");
927 sprintf(chisto,
"meResyZmPanel2");
929 sprintf(chisto,
"mePullxZmPanel2");
931 sprintf(chisto,
"mePullyZmPanel2");
933 sprintf(chisto,
"meNpixZmPanel2");
935 sprintf(chisto,
"meNxpixZmPanel2");
937 sprintf(chisto,
"meNypixZmPanel2");
939 sprintf(chisto,
"meChargeZmPanel2");
941 sprintf(chisto,
"meResXvsAlphaZmPanel2");
943 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
944 sprintf(chisto,
"meResYvsAlphaZmPanel2");
946 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
947 sprintf(chisto,
"meResXvsBetaZmPanel2");
949 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh,
"");
950 sprintf(chisto,
"meResYvsBetaZmPanel2");
952 ibooker.
bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh,
"");
954 sprintf(chisto,
"mePullXvsAlphaZmPanel2");
956 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
957 sprintf(chisto,
"mePullYvsAlphaZmPanel2");
959 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
960 sprintf(chisto,
"mePullXvsBetaZmPanel2");
962 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
963 sprintf(chisto,
"mePullYvsBetaZmPanel2");
965 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
966 sprintf(chisto,
"mePullXvsPhiZmPanel2");
968 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
969 sprintf(chisto,
"mePullYvsPhiZmPanel2");
971 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
972 sprintf(chisto,
"mePullXvsEtaZmPanel2");
974 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh,
"");
975 sprintf(chisto,
"mePullYvsEtaZmPanel2");
977 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh,
"");
979 sprintf(chisto,
"meWPullXvsAlphaZmPanel2");
981 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
982 sprintf(chisto,
"meWPullYvsAlphaZmPanel2");
984 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
985 sprintf(chisto,
"meWPullXvsBetaZmPanel2");
987 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh,
"");
988 sprintf(chisto,
"meWPullYvsBetaZmPanel2");
990 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh,
"");
992 sprintf(chisto,
"mePosxZpPanel2");
994 sprintf(chisto,
"mePosyZpPanel2");
996 sprintf(chisto,
"meErrxZpPanel2");
998 sprintf(chisto,
"meErryZpPanel2");
1000 sprintf(chisto,
"meResxZpPanel2");
1002 sprintf(chisto,
"meResyZpPanel2");
1004 sprintf(chisto,
"mePullxZpPanel2");
1006 sprintf(chisto,
"mePullyZpPanel2");
1008 sprintf(chisto,
"meNpixZpPanel2");
1010 sprintf(chisto,
"meNxpixZpPanel2");
1012 sprintf(chisto,
"meNypixZpPanel2");
1014 sprintf(chisto,
"meChargeZpPanel2");
1016 sprintf(chisto,
"meResXvsAlphaZpPanel2");
1018 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh,
"");
1019 sprintf(chisto,
"meResYvsAlphaZpPanel2");
1021 ibooker.
bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh,
"");
1022 sprintf(chisto,
"meResXvsBetaZpPanel2");
1024 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh,
"");
1025 sprintf(chisto,
"meResYvsBetaZpPanel2");
1027 ibooker.
bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh,
"");
1029 sprintf(chisto,
"mePullXvsAlphaZpPanel2");
1031 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
1032 sprintf(chisto,
"mePullYvsAlphaZpPanel2");
1034 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
1035 sprintf(chisto,
"mePullXvsBetaZpPanel2");
1037 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
1038 sprintf(chisto,
"mePullYvsBetaZpPanel2");
1040 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
1041 sprintf(chisto,
"mePullXvsPhiZpPanel2");
1043 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh,
"");
1044 sprintf(chisto,
"mePullYvsPhiZpPanel2");
1046 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh,
"");
1047 sprintf(chisto,
"mePullXvsEtaZpPanel2");
1049 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh,
"");
1050 sprintf(chisto,
"mePullYvsEtaZpPanel2");
1052 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh,
"");
1054 sprintf(chisto,
"meWPullXvsAlphaZpPanel2");
1056 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh,
"");
1057 sprintf(chisto,
"meWPullYvsAlphaZpPanel2");
1059 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh,
"");
1060 sprintf(chisto,
"meWPullXvsBetaZpPanel2");
1062 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh,
"");
1063 sprintf(chisto,
"meWPullYvsBetaZpPanel2");
1065 ibooker.
bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh,
"");
1068 sprintf(chisto,
"mePosxBarrel_all_hits");
1070 sprintf(chisto,
"mePosyBarrel_all_hits");
1073 sprintf(chisto,
"mePosxZmPanel1_all_hits");
1075 sprintf(chisto,
"mePosyZmPanel1_all_hits");
1077 sprintf(chisto,
"mePosxZmPanel2_all_hits");
1079 sprintf(chisto,
"mePosyZmPanel2_all_hits");
1082 sprintf(chisto,
"mePosxZpPanel1_all_hits");
1084 sprintf(chisto,
"mePosyZpPanel1_all_hits");
1086 sprintf(chisto,
"mePosxZpPanel2_all_hits");
1088 sprintf(chisto,
"mePosyZpPanel2_all_hits");
1110 evt =
e.id().event();
1116 float radtodeg = 180.0 /
math_pi;
1122 float mindist = 999999.9;
1142 for (TrackerGeometry::DetContainer::const_iterator it = pDD->
dets().begin(); it != pDD->
dets().end(); it++) {
1143 DetId detId = ((*it)->geographicalId());
1145 unsigned int subid = detId.
subdetId();
1146 if (!((subid == 1) || (subid == 2)))
1160 for (; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
1165 detId = (*it)->geographicalId();
1181 }
else if (
panel == 2) {
1186 <<
"..............................................Wrong panel number !";
1188 else if (
side == 2) {
1192 }
else if (
panel == 2) {
1197 <<
"..............................................Wrong panel number !";
1201 <<
".......................................................Wrong side !";
1206 <<
"Pixel rechit collection but we are not in the pixel detector" << (
int)detId.
subdetId();
1215 reco::TrackCollection::const_iterator tciter;
1222 for (tciter =
tracks->begin(); tciter !=
tracks->end(); tciter++) {
1223 phi = tciter->momentum().phi() /
math_pi * 180.0;
1224 eta = tciter->momentum().eta();
1231 const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
1272 error = (*it)->localPositionError();
1280 npix = (*matchedhit).cluster()->size();
1281 nxpix = (*matchedhit).cluster()->sizeX();
1282 nypix = (*matchedhit).cluster()->sizeY();
1283 charge = (*matchedhit).cluster()->charge();
1293 float distx, disty, dist;
1294 bool found_hit_from_generated_particle =
false;
1296 int n_assoc_muon = 0;
1298 std::vector<PSimHit>::const_iterator closestit =
matched.begin();
1299 for (std::vector<PSimHit>::const_iterator
m =
matched.begin();
m <
matched.end();
m++) {
1301 int pid = (*m).particleType();
1306 float simhitx = 0.5 * ((*m).entryPoint().x() + (*m).exitPoint().x());
1307 float simhity = 0.5 * ((*m).entryPoint().y() + (*m).exitPoint().y());
1311 dist =
sqrt(distx * distx + disty * disty);
1313 if (dist < mindist) {
1318 found_hit_from_generated_particle =
true;
1324 if (
checkType_ && !found_hit_from_generated_particle)
1327 if (n_assoc_muon > 1) {
1329 <<
" ----- This is not good: n_assoc_muon = " << n_assoc_muon;
1333 pidhit = (*closestit).particleType();
1335 simhitx = 0.5 * ((*closestit).entryPoint().x() + (*closestit).exitPoint().x());
1336 simhity = 0.5 * ((*closestit).entryPoint().y() + (*closestit).exitPoint().y());
1343 float simhitpx = (*closestit).momentumAtEntry().x();
1344 float simhitpy = (*closestit).momentumAtEntry().y();
1345 float simhitpz = (*closestit).momentumAtEntry().z();
1350 beta = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
1351 alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;
1353 detId = (*it)->geographicalId();
1388 if (tmp_nrows == 80) {
1392 }
else if (tmp_nrows == 160) {
1398 <<
"-------------------------------------------------- Wrong module size !!!";
1559 }
else if (
panel == 2) {
1617 <<
"..............................................Wrong panel number !";
1619 else if (
side == 2) {
1676 }
else if (
panel == 2) {
1734 <<
"..............................................Wrong panel number !";
1738 <<
".......................................................Wrong side !";
1743 <<
"Pixel rechit but we are not in the pixel detector" << (
int)detId.
subdetId();