21 TFile
f1(
"/afs/cern.ch/cms/tracking/output/rootfiles/seedmultiplicity_Run2011A_express_ge1711_v3_v2.root");
23 "/afs/cern.ch/cms/tracking/output/rootfiles/"
24 "seedmultiplicity_highmult_default_Run2011A_express_ge1711_v3_v2.root");
25 TFile f3(
"/afs/cern.ch/cms/tracking/output/rootfiles/seedmultiplicity_highmult_Run2011A_express_ge1711_v3_v2.root");
31 TH2F* iter3Pixel1 = (TH2F*)castat1.
getObject(
"thTripletsPixel");
32 TH2F* iter3Pixel2 = (TH2F*)castat2.
getObject(
"thTripletsPixel");
33 TH2F* iter3Pixel3 = (TH2F*)castat3.
getObject(
"thTripletsPixel");
35 if (iter3Pixel1 && iter3Pixel2 && iter3Pixel3) {
36 TProfile* iter3Pixel1_prof = iter3Pixel1->ProfileX(
"prof1");
37 TProfile* iter3Pixel2_prof = iter3Pixel2->ProfileX(
"prof2");
38 TProfile* iter3Pixel3_prof = iter3Pixel3->ProfileX(
"prof3");
40 iter3Pixel1_prof->SetMarkerColor(kBlack);
41 iter3Pixel1_prof->SetLineColor(kBlack);
42 iter3Pixel2_prof->SetMarkerColor(kBlue);
43 iter3Pixel2_prof->SetLineColor(kBlue);
44 iter3Pixel3_prof->SetMarkerColor(kRed);
45 iter3Pixel3_prof->SetLineColor(kRed);
47 new TCanvas(
"profile",
"profile");
49 iter3Pixel1_prof->GetXaxis()->SetRangeUser(0, 10000);
50 iter3Pixel1_prof->GetXaxis()->SetTitle(
"clusters");
51 iter3Pixel1_prof->GetYaxis()->SetTitle(
"seeds");
53 TLegend* leg =
new TLegend(.2, .6, .5, .85,
"Seeds vs clusters");
54 leg->AddEntry(iter3Pixel1_prof->DrawCopy(),
"standard RECO",
"l");
55 leg->AddEntry(iter3Pixel2_prof->DrawCopy(
"same"),
"standard RECO (high stat)",
"l");
56 leg->AddEntry(iter3Pixel3_prof->DrawCopy(
"same"),
"iter2 thr = 2M",
"l");
69 TH1F* iter3_2 = (TH1F*)castat2.
getObject(
"thTriplets");
70 TH1F* iter3_3 = (TH1F*)castat3.
getObject(
"thTriplets");
72 if (iter3_2 && iter3_3) {
73 iter3_2->SetLineColor(kBlue);
74 iter3_3->SetLineColor(kRed);
76 new TCanvas(
"iter3",
"iter3");
79 TLegend* legiter3 =
new TLegend(.5, .6, .85, .85,
"iter3 seeds");
80 legiter3->AddEntry(iter3_2->DrawCopy(),
"standard RECO",
"l");
81 legiter3->AddEntry(iter3_3->DrawCopy(
"same"),
"iter2 thr = 2M",
"l");
85 TH1F* iter2_2 = (TH1F*)castat2.
getObject(
"secTriplets");
86 TH1F* iter2_3 = (TH1F*)castat3.
getObject(
"secTriplets");
88 if (iter2_2 && iter2_3) {
89 iter2_2->SetLineColor(kBlue);
90 iter2_3->SetLineColor(kRed);
92 new TCanvas(
"iter2",
"iter2");
95 TLegend* legiter2 =
new TLegend(.5, .6, .85, .85,
"iter2 seeds");
96 legiter2->AddEntry(iter2_2->DrawCopy(),
"standard RECO",
"l");
97 legiter2->AddEntry(iter2_3->DrawCopy(
"same"),
"iter2 thr = 2M",
"l");
104 const char* fullname,
const char*
module,
const char* postfix,
const char* shortname,
const char* outtrunk) {
108 sprintf(modfull,
"%s%s",
module, postfix);
112 sprintf(
dirname,
"%s%s",
"seedmultiplicity_", shortname);
128 gStyle->SetOptStat(111111);
134 TH1F* iter0 = (TH1F*)castat.
getObject(
"newSeedFromTriplets");
139 plotfilename += outtrunk;
141 plotfilename +=
"/iter0_";
143 plotfilename +=
".gif";
144 gPad->Print(plotfilename.c_str());
149 TH1F* iter1 = (TH1F*)castat.
getObject(
"newSeedFromPairs");
154 plotfilename += outtrunk;
156 plotfilename +=
"/iter1_";
158 plotfilename +=
".gif";
159 gPad->Print(plotfilename.c_str());
164 TH1F* iter2 = (TH1F*)castat.
getObject(
"secTriplets");
169 plotfilename += outtrunk;
171 plotfilename +=
"/iter2_";
173 plotfilename +=
".gif";
174 gPad->Print(plotfilename.c_str());
179 TH1F* iter3 = (TH1F*)castat.
getObject(
"thTriplets");
184 plotfilename += outtrunk;
186 plotfilename +=
"/iter3_";
188 plotfilename +=
".gif";
189 gPad->Print(plotfilename.c_str());
194 TH1F* iter3A = (TH1F*)castat.
getObject(
"thTripletsA");
199 plotfilename += outtrunk;
201 plotfilename +=
"/iter3A_";
203 plotfilename +=
".gif";
204 gPad->Print(plotfilename.c_str());
209 TH1F* iter3B = (TH1F*)castat.
getObject(
"thTripletsB");
214 plotfilename += outtrunk;
216 plotfilename +=
"/iter3B_";
218 plotfilename +=
".gif";
219 gPad->Print(plotfilename.c_str());
224 TH1F* iter4 = (TH1F*)castat.
getObject(
"fourthPLSeeds");
229 plotfilename += outtrunk;
231 plotfilename +=
"/iter4_";
233 plotfilename +=
".gif";
234 gPad->Print(plotfilename.c_str());
239 TH1F* iter5 = (TH1F*)castat.
getObject(
"fifthSeeds");
244 plotfilename += outtrunk;
246 plotfilename +=
"/iter5_";
248 plotfilename +=
".gif";
249 gPad->Print(plotfilename.c_str());
254 TH2F* iter0TK = (TH2F*)castat.
getObject(
"newSeedFromTripletsTK");
256 iter0TK->Draw(
"colz");
257 iter0TK->GetYaxis()->SetRangeUser(0, 1000);
258 iter0TK->GetYaxis()->SetTitle(
"seeds");
259 iter0TK->GetXaxis()->SetRangeUser(0, 50000);
260 iter0TK->GetXaxis()->SetTitle(
"clusters");
263 plotfilename += outtrunk;
265 plotfilename +=
"/iter0TK_";
267 plotfilename +=
".gif";
268 gPad->Print(plotfilename.c_str());
273 TH2F* iter1TK = (TH2F*)castat.
getObject(
"newSeedFromPairsTK");
275 iter1TK->Draw(
"colz");
276 iter1TK->GetYaxis()->SetTitle(
"seeds");
277 iter1TK->GetXaxis()->SetRangeUser(0, 50000);
278 iter1TK->GetXaxis()->SetTitle(
"clusters");
281 plotfilename += outtrunk;
283 plotfilename +=
"/iter1TK_";
285 plotfilename +=
".gif";
286 gPad->Print(plotfilename.c_str());
291 TH2F* iter2TK = (TH2F*)castat.
getObject(
"secTripletsTK");
293 iter2TK->Draw(
"colz");
294 iter2TK->GetYaxis()->SetRangeUser(0, 50000);
295 iter2TK->GetYaxis()->SetTitle(
"seeds");
296 iter2TK->GetXaxis()->SetRangeUser(0, 50000);
297 iter2TK->GetXaxis()->SetTitle(
"clusters");
300 plotfilename += outtrunk;
302 plotfilename +=
"/iter2TK_";
304 plotfilename +=
".gif";
305 gPad->Print(plotfilename.c_str());
310 TH2F* iter3TK = (TH2F*)castat.
getObject(
"thTripletsTK");
312 iter3TK->Draw(
"colz");
314 iter3TK->GetYaxis()->SetTitle(
"seeds");
315 iter3TK->GetXaxis()->SetRangeUser(0, 50000);
316 iter3TK->GetXaxis()->SetTitle(
"clusters");
319 plotfilename += outtrunk;
321 plotfilename +=
"/iter3TK_";
323 plotfilename +=
".gif";
324 gPad->Print(plotfilename.c_str());
329 TH2F* iter3ATK = (TH2F*)castat.
getObject(
"thTripletsATK");
331 iter3ATK->Draw(
"colz");
333 iter3ATK->GetYaxis()->SetTitle(
"seeds");
334 iter3ATK->GetXaxis()->SetRangeUser(0, 50000);
335 iter3ATK->GetXaxis()->SetTitle(
"clusters");
338 plotfilename += outtrunk;
340 plotfilename +=
"/iter3ATK_";
342 plotfilename +=
".gif";
343 gPad->Print(plotfilename.c_str());
348 TH2F* iter3BTK = (TH2F*)castat.
getObject(
"thTripletsBTK");
350 iter3BTK->Draw(
"colz");
352 iter3BTK->GetYaxis()->SetTitle(
"seeds");
353 iter3BTK->GetXaxis()->SetRangeUser(0, 50000);
354 iter3BTK->GetXaxis()->SetTitle(
"clusters");
357 plotfilename += outtrunk;
359 plotfilename +=
"/iter3BTK_";
361 plotfilename +=
".gif";
362 gPad->Print(plotfilename.c_str());
367 TH2F* iter4TK = (TH2F*)castat.
getObject(
"fourthPLSeedsTK");
369 iter4TK->Draw(
"colz");
370 iter4TK->GetYaxis()->SetTitle(
"seeds");
371 iter4TK->GetXaxis()->SetRangeUser(0, 50000);
372 iter4TK->GetXaxis()->SetTitle(
"clusters");
375 plotfilename += outtrunk;
377 plotfilename +=
"/iter4TK_";
379 plotfilename +=
".gif";
380 gPad->Print(plotfilename.c_str());
385 TH2F* iter5TK = (TH2F*)castat.
getObject(
"fifthSeedsTK");
387 iter5TK->Draw(
"colz");
388 iter5TK->GetYaxis()->SetRangeUser(0, 30000);
389 iter5TK->GetYaxis()->SetTitle(
"seeds");
390 iter5TK->GetXaxis()->SetRangeUser(0, 50000);
391 iter5TK->GetXaxis()->SetTitle(
"clusters");
394 plotfilename += outtrunk;
396 plotfilename +=
"/iter5TK_";
398 plotfilename +=
".gif";
399 gPad->Print(plotfilename.c_str());
404 TH2F* iter0Pixel = (TH2F*)castat.
getObject(
"newSeedFromTripletsPixel");
406 iter0Pixel->Draw(
"colz");
407 iter0Pixel->GetYaxis()->SetRangeUser(0, 1000);
408 iter0Pixel->GetYaxis()->SetTitle(
"seeds");
409 iter0Pixel->GetXaxis()->SetRangeUser(0, 10000);
410 iter0Pixel->GetXaxis()->SetTitle(
"clusters");
413 plotfilename += outtrunk;
415 plotfilename +=
"/iter0Pixel_";
417 plotfilename +=
".gif";
418 gPad->Print(plotfilename.c_str());
423 TH2F* iter1Pixel = (TH2F*)castat.
getObject(
"newSeedFromPairsPixel");
425 iter1Pixel->Draw(
"colz");
426 iter1Pixel->GetYaxis()->SetTitle(
"seeds");
427 iter1Pixel->GetXaxis()->SetRangeUser(0, 10000);
428 iter1Pixel->GetXaxis()->SetTitle(
"clusters");
431 plotfilename += outtrunk;
433 plotfilename +=
"/iter1Pixel_";
435 plotfilename +=
".gif";
436 gPad->Print(plotfilename.c_str());
441 TH2F* iter2Pixel = (TH2F*)castat.
getObject(
"secTripletsPixel");
443 iter2Pixel->Draw(
"colz");
444 iter2Pixel->GetYaxis()->SetRangeUser(0, 50000);
445 iter2Pixel->GetYaxis()->SetTitle(
"seeds");
446 iter2Pixel->GetXaxis()->SetRangeUser(0, 10000);
447 iter2Pixel->GetXaxis()->SetTitle(
"clusters");
450 plotfilename += outtrunk;
452 plotfilename +=
"/iter2Pixel_";
454 plotfilename +=
".gif";
455 gPad->Print(plotfilename.c_str());
460 TH2F* iter3Pixel = (TH2F*)castat.
getObject(
"thTripletsPixel");
462 iter3Pixel->Draw(
"colz");
464 iter3Pixel->GetYaxis()->SetTitle(
"seeds");
465 iter3Pixel->GetXaxis()->SetRangeUser(0, 10000);
466 iter3Pixel->GetXaxis()->SetTitle(
"clusters");
469 plotfilename += outtrunk;
471 plotfilename +=
"/iter3Pixel_";
473 plotfilename +=
".gif";
474 gPad->Print(plotfilename.c_str());
479 TH2F* iter3APixel = (TH2F*)castat.
getObject(
"thTripletsAPixel");
481 iter3APixel->Draw(
"colz");
483 iter3APixel->GetYaxis()->SetTitle(
"seeds");
484 iter3APixel->GetXaxis()->SetRangeUser(0, 10000);
485 iter3APixel->GetXaxis()->SetTitle(
"clusters");
488 plotfilename += outtrunk;
490 plotfilename +=
"/iter3APixel_";
492 plotfilename +=
".gif";
493 gPad->Print(plotfilename.c_str());
498 TH2F* iter3BPixel = (TH2F*)castat.
getObject(
"thTripletsBPixel");
500 iter3BPixel->Draw(
"colz");
502 iter3BPixel->GetYaxis()->SetTitle(
"seeds");
503 iter3BPixel->GetXaxis()->SetRangeUser(0, 10000);
504 iter3BPixel->GetXaxis()->SetTitle(
"clusters");
507 plotfilename += outtrunk;
509 plotfilename +=
"/iter3BPixel_";
511 plotfilename +=
".gif";
512 gPad->Print(plotfilename.c_str());
517 TH2F* iter4Pixel = (TH2F*)castat.
getObject(
"fourthPLSeedsPixel");
519 iter4Pixel->Draw(
"colz");
520 iter4Pixel->GetYaxis()->SetTitle(
"seeds");
521 iter4Pixel->GetXaxis()->SetRangeUser(0, 10000);
522 iter4Pixel->GetXaxis()->SetTitle(
"clusters");
525 plotfilename += outtrunk;
527 plotfilename +=
"/iter4Pixel_";
529 plotfilename +=
".gif";
530 gPad->Print(plotfilename.c_str());
535 TH2F* iter5Pixel = (TH2F*)castat.
getObject(
"fifthSeedsPixel");
537 iter5Pixel->Draw(
"colz");
538 iter5Pixel->GetYaxis()->SetRangeUser(0, 30000);
539 iter5Pixel->GetYaxis()->SetTitle(
"seeds");
540 iter5Pixel->GetXaxis()->SetRangeUser(0, 10000);
541 iter5Pixel->GetXaxis()->SetTitle(
"clusters");
544 plotfilename += outtrunk;
546 plotfilename +=
"/iter5Pixel_";
548 plotfilename +=
".gif";
549 gPad->Print(plotfilename.c_str());
553 gStyle->SetOptStat(1111);