19 #include <TPaveText.h> 20 #include <TPostScript.h> 33 printf(
"reco: gROOT Reset \n");
35 gROOT->SetStyle(
"Plain");
36 gStyle->SetOptStat(0);
47 gStyle->SetOptFit(0010);
56 gStyle->SetStatX(0.91);
57 gStyle->SetStatY(0.75);
58 gStyle->SetStatW(0.20);
59 gStyle->SetStatH(0.10);
65 Float_t LeftOffset = 0.12;
66 Float_t TopOffset = 0.22;
68 gStyle->SetLineWidth(1);
72 gStyle->SetOptTitle(0);
74 gStyle->SetTitleColor(1);
76 gStyle->SetTitleX(0.15);
77 gStyle->SetTitleH(0.15);
79 gStyle->SetTitleW(0.60);
80 gStyle->SetTitleFont(42);
81 gStyle->SetTitleFontSize(0.07);
83 gStyle->SetPalette(1);
96 gStyle->SetPadTopMargin(TopOffset);
97 gStyle->SetPadBottomMargin(LeftOffset);
98 gStyle->SetPadRightMargin(TopOffset);
99 gStyle->SetPadLeftMargin(LeftOffset);
111 TFile *hfile1 =
new TFile(
"Global_362596.root",
"READ");
114 TPostScript psfile(
"ztsmean.ps", 111);
118 TCanvas *
c1 =
new TCanvas(
"c1",
"Hcal4test", 200, 10, 700, 900);
121 TDirectory *
dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
131 TH2F *twod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA_HB");
132 TH2F *twod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HB");
136 TH2F *Ceff = (TH2F *)twod1->Clone(
"Ceff");
137 Ceff->Divide(twod1, twod0, 1, 1,
"B");
141 TH2F *twop1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA225_HB");
142 TH2F *twop0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HB");
146 TH2F *Cefz225 = (TH2F *)twop1->Clone(
"Cefz225");
147 Cefz225->Divide(twop1, twop0, 1, 1,
"B");
153 Cefz225->SetMarkerStyle(20);
154 Cefz225->SetMarkerSize(0.4);
156 Cefz225->GetZaxis()->SetLabelSize(0.08);
157 Cefz225->SetXTitle(
"#eta \b");
158 Cefz225->SetYTitle(
"#phi \b");
159 Cefz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth1 \b");
160 Cefz225->SetMarkerColor(2);
161 Cefz225->SetLineColor(2);
162 Cefz225->SetMaximum(1.000);
163 Cefz225->SetMinimum(0.0001);
164 Cefz225->Draw(
"COLZ");
167 TH1F *aaaaaa1 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HB");
169 aaaaaa1->SetMarkerStyle(20);
170 aaaaaa1->SetMarkerSize(0.8);
171 aaaaaa1->GetYaxis()->SetLabelSize(0.04);
172 aaaaaa1->SetXTitle(
"TSmeanA in each event & cell HB \b");
173 aaaaaa1->SetMarkerColor(2);
174 aaaaaa1->SetLineColor(2);
179 TH2F *Diffe_Depth1_HB = (TH2F *)Ceff->Clone(
"Diffe_Depth1_HB");
180 for (
int i = 1;
i <= Ceff->GetXaxis()->GetNbins();
i++) {
181 for (
int j = 1;
j <= Ceff->GetYaxis()->GetNbins();
j++) {
182 double ccc1 = Ceff->GetBinContent(
i,
j);
183 Diffe_Depth1_HB->SetBinContent(
i,
j, 0.);
185 if (ccc1 < 2.5 || ccc1 > 4.5)
186 Diffe_Depth1_HB->SetBinContent(
i,
j, ccc1);
192 Diffe_Depth1_HB->SetMarkerStyle(20);
193 Diffe_Depth1_HB->SetMarkerSize(0.4);
194 Diffe_Depth1_HB->GetZaxis()->SetLabelSize(0.08);
196 Diffe_Depth1_HB->SetXTitle(
"#eta \b");
197 Diffe_Depth1_HB->SetYTitle(
"#phi \b");
198 Diffe_Depth1_HB->SetZTitle(
"<TSmeanA> out 2.5-4.5- HB Depth1 \b");
199 Diffe_Depth1_HB->SetMarkerColor(2);
200 Diffe_Depth1_HB->SetLineColor(2);
201 Diffe_Depth1_HB->Draw(
"COLZ");
204 TH1F *diffTSmeanA_Depth1_HB =
new TH1F(
"diffTSmeanA_Depth1_HB",
"", 100, 1.0, 6.0);
205 for (
int i = 1;
i <= Ceff->GetXaxis()->GetNbins();
i++) {
206 for (
int j = 1;
j <= Ceff->GetYaxis()->GetNbins();
j++) {
207 if (Ceff->GetBinContent(
i,
j) != 0) {
208 double ccc1 = Ceff->GetBinContent(
i,
j);
209 diffTSmeanA_Depth1_HB->Fill(ccc1);
214 diffTSmeanA_Depth1_HB->SetMarkerStyle(20);
215 diffTSmeanA_Depth1_HB->SetMarkerSize(0.4);
216 diffTSmeanA_Depth1_HB->GetYaxis()->SetLabelSize(0.04);
217 diffTSmeanA_Depth1_HB->SetXTitle(
"<TSmeanA> in each cell \b");
218 diffTSmeanA_Depth1_HB->SetMarkerColor(2);
219 diffTSmeanA_Depth1_HB->SetLineColor(2);
220 diffTSmeanA_Depth1_HB->Draw(
"");
232 TH2F *awod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA_HB");
233 TH2F *awod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HB");
237 TH2F *C2ff = (TH2F *)awod1->Clone(
"C2ff");
238 C2ff->Divide(awod1, awod0, 1, 1,
"B");
242 TH2F *bwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA225_HB");
243 TH2F *bwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HB");
247 TH2F *C2fz225 = (TH2F *)bwod1->Clone(
"C2fz225");
248 C2fz225->Divide(bwod1, bwod0, 1, 1,
"B");
254 C2fz225->SetMarkerStyle(20);
255 C2fz225->SetMarkerSize(0.4);
256 C2fz225->GetZaxis()->SetLabelSize(0.08);
257 C2fz225->SetXTitle(
"#eta \b");
258 C2fz225->SetYTitle(
"#phi \b");
259 C2fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth2 \b");
260 C2fz225->SetMarkerColor(2);
261 C2fz225->SetLineColor(2);
262 C2fz225->SetMaximum(1.000);
263 C2fz225->SetMinimum(0.0001);
264 C2fz225->Draw(
"COLZ");
267 TH1F *aaaaaa2 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HB");
269 aaaaaa2->SetMarkerStyle(20);
270 aaaaaa2->SetMarkerSize(0.8);
271 aaaaaa2->GetYaxis()->SetLabelSize(0.04);
272 aaaaaa2->SetXTitle(
"TSmeanA in each event & cell HB \b");
273 aaaaaa2->SetMarkerColor(2);
274 aaaaaa2->SetLineColor(2);
279 TH2F *Diffe_Depth2_HB = (TH2F *)C2ff->Clone(
"Diffe_Depth2_HB");
280 for (
int i = 1;
i <= C2ff->GetXaxis()->GetNbins();
i++) {
281 for (
int j = 1;
j <= C2ff->GetYaxis()->GetNbins();
j++) {
282 double ccc1 = C2ff->GetBinContent(
i,
j);
283 Diffe_Depth2_HB->SetBinContent(
i,
j, 0.);
284 if (ccc1 < 2.4 || ccc1 > 4.4)
285 Diffe_Depth2_HB->SetBinContent(
i,
j, ccc1);
291 Diffe_Depth2_HB->SetMarkerStyle(20);
292 Diffe_Depth2_HB->SetMarkerSize(0.4);
293 Diffe_Depth2_HB->GetZaxis()->SetLabelSize(0.08);
295 Diffe_Depth2_HB->SetXTitle(
"#eta \b");
296 Diffe_Depth2_HB->SetYTitle(
"#phi \b");
297 Diffe_Depth2_HB->SetZTitle(
"<TSmeanA> out 2.4-4.4 - HB Depth2 \b");
298 Diffe_Depth2_HB->SetMarkerColor(2);
299 Diffe_Depth2_HB->SetLineColor(2);
300 Diffe_Depth2_HB->Draw(
"COLZ");
303 TH1F *diffTSmeanA_Depth2_HB =
new TH1F(
"diffTSmeanA_Depth2_HB",
"", 100, 1.0, 6.0);
304 for (
int i = 1;
i <= C2ff->GetXaxis()->GetNbins();
i++) {
305 for (
int j = 1;
j <= C2ff->GetYaxis()->GetNbins();
j++) {
306 if (C2ff->GetBinContent(
i,
j) != 0) {
307 double ccc1 = C2ff->GetBinContent(
i,
j);
308 diffTSmeanA_Depth2_HB->Fill(ccc1);
313 diffTSmeanA_Depth2_HB->SetMarkerStyle(20);
314 diffTSmeanA_Depth2_HB->SetMarkerSize(0.4);
315 diffTSmeanA_Depth2_HB->GetYaxis()->SetLabelSize(0.04);
316 diffTSmeanA_Depth2_HB->SetXTitle(
"<TSmeanA> in each cell \b");
317 diffTSmeanA_Depth2_HB->SetMarkerColor(2);
318 diffTSmeanA_Depth2_HB->SetLineColor(2);
319 diffTSmeanA_Depth2_HB->Draw(
"");
331 TH2F *cwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA_HE");
332 TH2F *cwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HE");
336 TH2F *C3ff = (TH2F *)cwod1->Clone(
"C3ff");
337 C3ff->Divide(cwod1, cwod0, 1, 1,
"B");
341 TH2F *dwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA225_HE");
342 TH2F *dwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HE");
346 TH2F *C3fz225 = (TH2F *)dwod1->Clone(
"C3fz225");
347 C3fz225->Divide(dwod1, dwod0, 1, 1,
"B");
353 C3fz225->SetMarkerStyle(20);
354 C3fz225->SetMarkerSize(0.4);
355 C3fz225->GetZaxis()->SetLabelSize(0.08);
356 C3fz225->SetXTitle(
"#eta \b");
357 C3fz225->SetYTitle(
"#phi \b");
358 C3fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth1 \b");
359 C3fz225->SetMarkerColor(2);
360 C3fz225->SetLineColor(2);
361 C3fz225->SetMaximum(1.000);
362 C3fz225->SetMinimum(0.0001);
363 C3fz225->Draw(
"COLZ");
366 TH1F *aaaaaa3 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HE");
368 aaaaaa3->SetMarkerStyle(20);
369 aaaaaa3->SetMarkerSize(0.8);
370 aaaaaa3->GetYaxis()->SetLabelSize(0.04);
371 aaaaaa3->SetXTitle(
"TSmeanA in each event & cell HE \b");
372 aaaaaa3->SetMarkerColor(2);
373 aaaaaa3->SetLineColor(2);
378 TH2F *Diffe_Depth1_HE = (TH2F *)C3ff->Clone(
"Diffe_Depth1_HE");
379 for (
int i = 1;
i <= C3ff->GetXaxis()->GetNbins();
i++) {
380 for (
int j = 1;
j <= C3ff->GetYaxis()->GetNbins();
j++) {
381 double ccc1 = C3ff->GetBinContent(
i,
j);
382 Diffe_Depth1_HE->SetBinContent(
i,
j, 0.);
383 if (ccc1 < 2.4 || ccc1 > 4.6)
384 Diffe_Depth1_HE->SetBinContent(
i,
j, ccc1);
390 Diffe_Depth1_HE->SetMarkerStyle(20);
391 Diffe_Depth1_HE->SetMarkerSize(0.4);
392 Diffe_Depth1_HE->GetZaxis()->SetLabelSize(0.08);
394 Diffe_Depth1_HE->SetXTitle(
"#eta \b");
395 Diffe_Depth1_HE->SetYTitle(
"#phi \b");
396 Diffe_Depth1_HE->SetZTitle(
"<TSmeanA> out 2.4-4.6 - HE Depth1 \b");
397 Diffe_Depth1_HE->SetMarkerColor(2);
398 Diffe_Depth1_HE->SetLineColor(2);
399 Diffe_Depth1_HE->Draw(
"COLZ");
402 TH1F *diffTSmeanA_Depth1_HE =
new TH1F(
"diffTSmeanA_Depth1_HE",
"", 100, 0.0, 7.0);
403 for (
int i = 1;
i <= C3ff->GetXaxis()->GetNbins();
i++) {
404 for (
int j = 1;
j <= C3ff->GetYaxis()->GetNbins();
j++) {
405 if (C3ff->GetBinContent(
i,
j) != 0) {
406 double ccc1 = C3ff->GetBinContent(
i,
j);
407 diffTSmeanA_Depth1_HE->Fill(ccc1);
412 diffTSmeanA_Depth1_HE->SetMarkerStyle(20);
413 diffTSmeanA_Depth1_HE->SetMarkerSize(0.4);
414 diffTSmeanA_Depth1_HE->GetYaxis()->SetLabelSize(0.04);
415 diffTSmeanA_Depth1_HE->SetXTitle(
"<TSmeanA> in each cell \b");
416 diffTSmeanA_Depth1_HE->SetMarkerColor(2);
417 diffTSmeanA_Depth1_HE->SetLineColor(2);
418 diffTSmeanA_Depth1_HE->Draw(
"");
430 TH2F *ewod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA_HE");
431 TH2F *ewod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HE");
435 TH2F *C4ff = (TH2F *)ewod1->Clone(
"C4ff");
436 C4ff->Divide(ewod1, ewod0, 1, 1,
"B");
440 TH2F *fwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA225_HE");
441 TH2F *fwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HE");
445 TH2F *C4fz225 = (TH2F *)fwod1->Clone(
"C4fz225");
446 C4fz225->Divide(fwod1, fwod0, 1, 1,
"B");
452 C4fz225->SetMarkerStyle(20);
453 C4fz225->SetMarkerSize(0.4);
454 C4fz225->GetZaxis()->SetLabelSize(0.08);
455 C4fz225->SetXTitle(
"#eta \b");
456 C4fz225->SetYTitle(
"#phi \b");
457 C4fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth2 \b");
458 C4fz225->SetMarkerColor(2);
459 C4fz225->SetLineColor(2);
460 C4fz225->SetMaximum(1.000);
461 C4fz225->SetMinimum(0.0001);
462 C4fz225->Draw(
"COLZ");
465 TH1F *aaaaaa4 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HE");
467 aaaaaa4->SetMarkerStyle(20);
468 aaaaaa4->SetMarkerSize(0.8);
469 aaaaaa4->GetYaxis()->SetLabelSize(0.04);
470 aaaaaa4->SetXTitle(
"TSmeanA in each event & cell HE \b");
471 aaaaaa4->SetMarkerColor(2);
472 aaaaaa4->SetLineColor(2);
477 TH2F *Diffe_Depth2_HE = (TH2F *)C4ff->Clone(
"Diffe_Depth2_HE");
478 for (
int i = 1;
i <= C4ff->GetXaxis()->GetNbins();
i++) {
479 for (
int j = 1;
j <= C4ff->GetYaxis()->GetNbins();
j++) {
480 double ccc1 = C4ff->GetBinContent(
i,
j);
481 Diffe_Depth2_HE->SetBinContent(
i,
j, 0.);
482 if (ccc1 < 2.4 || ccc1 > 4.6)
483 Diffe_Depth2_HE->SetBinContent(
i,
j, ccc1);
489 Diffe_Depth2_HE->SetMarkerStyle(20);
490 Diffe_Depth2_HE->SetMarkerSize(0.4);
491 Diffe_Depth2_HE->GetZaxis()->SetLabelSize(0.08);
493 Diffe_Depth2_HE->SetXTitle(
"#eta \b");
494 Diffe_Depth2_HE->SetYTitle(
"#phi \b");
495 Diffe_Depth2_HE->SetZTitle(
"<TSmeanA> out 2.4-4.6 - HE Depth2 \b");
496 Diffe_Depth2_HE->SetMarkerColor(2);
497 Diffe_Depth2_HE->SetLineColor(2);
498 Diffe_Depth2_HE->Draw(
"COLZ");
501 TH1F *diffTSmeanA_Depth2_HE =
new TH1F(
"diffTSmeanA_Depth2_HE",
"", 100, 0.0, 7.0);
502 for (
int i = 1;
i <= C4ff->GetXaxis()->GetNbins();
i++) {
503 for (
int j = 1;
j <= C4ff->GetYaxis()->GetNbins();
j++) {
504 if (C4ff->GetBinContent(
i,
j) != 0) {
505 double ccc1 = C4ff->GetBinContent(
i,
j);
506 diffTSmeanA_Depth2_HE->Fill(ccc1);
511 diffTSmeanA_Depth2_HE->SetMarkerStyle(20);
512 diffTSmeanA_Depth2_HE->SetMarkerSize(0.4);
513 diffTSmeanA_Depth2_HE->GetYaxis()->SetLabelSize(0.04);
514 diffTSmeanA_Depth2_HE->SetXTitle(
"<TSmeanA> in each cell \b");
515 diffTSmeanA_Depth2_HE->SetMarkerColor(2);
516 diffTSmeanA_Depth2_HE->SetLineColor(2);
517 diffTSmeanA_Depth2_HE->Draw(
"");
529 TH2F *gwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth3TSmeanA_HE");
530 TH2F *gwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth3_HE");
534 TH2F *C5ff = (TH2F *)gwod1->Clone(
"C5ff");
535 C5ff->Divide(gwod1, gwod0, 1, 1,
"B");
539 TH2F *jwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth3TSmeanA225_HE");
540 TH2F *jwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth3_HE");
544 TH2F *C5fz225 = (TH2F *)jwod1->Clone(
"C5fz225");
545 C5fz225->Divide(jwod1, jwod0, 1, 1,
"B");
551 C5fz225->SetMarkerStyle(20);
552 C5fz225->SetMarkerSize(0.4);
553 C5fz225->GetZaxis()->SetLabelSize(0.08);
554 C5fz225->SetXTitle(
"#eta \b");
555 C5fz225->SetYTitle(
"#phi \b");
556 C5fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth3 \b");
557 C5fz225->SetMarkerColor(2);
558 C5fz225->SetLineColor(2);
559 C5fz225->SetMaximum(1.000);
560 C5fz225->SetMinimum(0.0001);
561 C5fz225->Draw(
"COLZ");
564 TH1F *aaaaaa5 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HE");
566 aaaaaa5->SetMarkerStyle(20);
567 aaaaaa5->SetMarkerSize(0.8);
568 aaaaaa5->GetYaxis()->SetLabelSize(0.04);
569 aaaaaa5->SetXTitle(
"TSmeanA in each event & cell HE \b");
570 aaaaaa5->SetMarkerColor(2);
571 aaaaaa5->SetLineColor(2);
576 TH2F *Diffe_Depth3_HE = (TH2F *)C5ff->Clone(
"Diffe_Depth3_HE");
577 for (
int i = 1;
i <= C5ff->GetXaxis()->GetNbins();
i++) {
578 for (
int j = 1;
j <= C5ff->GetYaxis()->GetNbins();
j++) {
579 double ccc1 = C5ff->GetBinContent(
i,
j);
580 Diffe_Depth3_HE->SetBinContent(
i,
j, 0.);
581 if (ccc1 < 2.4 || ccc1 > 4.6)
582 Diffe_Depth3_HE->SetBinContent(
i,
j, ccc1);
588 Diffe_Depth3_HE->SetMarkerStyle(20);
589 Diffe_Depth3_HE->SetMarkerSize(0.4);
590 Diffe_Depth3_HE->GetZaxis()->SetLabelSize(0.08);
592 Diffe_Depth3_HE->SetXTitle(
"#eta \b");
593 Diffe_Depth3_HE->SetYTitle(
"#phi \b");
594 Diffe_Depth3_HE->SetZTitle(
"<TSmeanA> out 2.4-4.6 - HE Depth3 \b");
595 Diffe_Depth3_HE->SetMarkerColor(2);
596 Diffe_Depth3_HE->SetLineColor(2);
597 Diffe_Depth3_HE->Draw(
"COLZ");
600 TH1F *diffTSmeanA_Depth3_HE =
new TH1F(
"diffTSmeanA_Depth3_HE",
"", 100, 0.0, 7.0);
601 for (
int i = 1;
i <= C5ff->GetXaxis()->GetNbins();
i++) {
602 for (
int j = 1;
j <= C5ff->GetYaxis()->GetNbins();
j++) {
603 if (C5ff->GetBinContent(
i,
j) != 0) {
604 double ccc1 = C5ff->GetBinContent(
i,
j);
605 diffTSmeanA_Depth3_HE->Fill(ccc1);
610 diffTSmeanA_Depth3_HE->SetMarkerStyle(20);
611 diffTSmeanA_Depth3_HE->SetMarkerSize(0.4);
612 diffTSmeanA_Depth3_HE->GetYaxis()->SetLabelSize(0.04);
613 diffTSmeanA_Depth3_HE->SetXTitle(
"<TSmeanA> in each cell \b");
614 diffTSmeanA_Depth3_HE->SetMarkerColor(2);
615 diffTSmeanA_Depth3_HE->SetLineColor(2);
616 diffTSmeanA_Depth3_HE->Draw(
"");
628 TH2F *iwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA_HF");
629 TH2F *iwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HF");
633 TH2F *C6ff = (TH2F *)iwod1->Clone(
"C6ff");
634 C6ff->Divide(iwod1, iwod0, 1, 1,
"B");
638 TH2F *kwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1TSmeanA225_HF");
639 TH2F *kwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth1_HF");
643 TH2F *C6fz225 = (TH2F *)kwod1->Clone(
"C6fz225");
644 C6fz225->Divide(kwod1, kwod0, 1, 1,
"B");
650 C6fz225->SetMarkerStyle(20);
651 C6fz225->SetMarkerSize(0.4);
653 C6fz225->GetZaxis()->SetLabelSize(0.08);
654 C6fz225->SetXTitle(
"#eta \b");
655 C6fz225->SetYTitle(
"#phi \b");
656 C6fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth1 \b");
657 C6fz225->SetMarkerColor(2);
658 C6fz225->SetLineColor(2);
659 C6fz225->SetMaximum(1.000);
660 C6fz225->SetMinimum(0.0001);
661 C6fz225->Draw(
"COLZ");
664 TH1F *aaaaaa6 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HF");
666 aaaaaa6->SetMarkerStyle(20);
667 aaaaaa6->SetMarkerSize(0.8);
668 aaaaaa6->GetYaxis()->SetLabelSize(0.04);
669 aaaaaa6->SetXTitle(
"TSmeanA in each event & cell HF \b");
670 aaaaaa6->SetMarkerColor(2);
671 aaaaaa6->SetLineColor(2);
676 TH2F *Diffe_Depth1_HF = (TH2F *)C6ff->Clone(
"Diffe_Depth1_HF");
677 for (
int i = 1;
i <= C6ff->GetXaxis()->GetNbins();
i++) {
678 for (
int j = 1;
j <= C6ff->GetYaxis()->GetNbins();
j++) {
679 double ccc1 = C6ff->GetBinContent(
i,
j);
680 Diffe_Depth1_HF->SetBinContent(
i,
j, 0.);
681 if (ccc1 < 0.5 || ccc1 > 1.5)
682 Diffe_Depth1_HF->SetBinContent(
i,
j, ccc1);
688 Diffe_Depth1_HF->SetMarkerStyle(20);
689 Diffe_Depth1_HF->SetMarkerSize(0.4);
690 Diffe_Depth1_HF->GetZaxis()->SetLabelSize(0.08);
692 Diffe_Depth1_HF->SetXTitle(
"#eta \b");
693 Diffe_Depth1_HF->SetYTitle(
"#phi \b");
694 Diffe_Depth1_HF->SetZTitle(
"<TSmeanA> out 0.5-1.5 - HF Depth1 \b");
695 Diffe_Depth1_HF->SetMarkerColor(2);
696 Diffe_Depth1_HF->SetLineColor(2);
697 Diffe_Depth1_HF->Draw(
"COLZ");
700 TH1F *diffTSmeanA_Depth1_HF =
new TH1F(
"diffTSmeanA_Depth1_HF",
"", 100, 0.0, 2.0);
701 for (
int i = 1;
i <= C6ff->GetXaxis()->GetNbins();
i++) {
702 for (
int j = 1;
j <= C6ff->GetYaxis()->GetNbins();
j++) {
703 if (C6ff->GetBinContent(
i,
j) != 0) {
704 double ccc1 = C6ff->GetBinContent(
i,
j);
705 diffTSmeanA_Depth1_HF->Fill(ccc1);
710 diffTSmeanA_Depth1_HF->SetMarkerStyle(20);
711 diffTSmeanA_Depth1_HF->SetMarkerSize(0.4);
712 diffTSmeanA_Depth1_HF->GetYaxis()->SetLabelSize(0.04);
713 diffTSmeanA_Depth1_HF->SetXTitle(
"<TSmeanA> in each cell \b");
714 diffTSmeanA_Depth1_HF->SetMarkerColor(2);
715 diffTSmeanA_Depth1_HF->SetLineColor(2);
716 diffTSmeanA_Depth1_HF->Draw(
"");
728 TH2F *lwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA_HF");
729 TH2F *lwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HF");
733 TH2F *C7ff = (TH2F *)lwod1->Clone(
"C7ff");
734 C7ff->Divide(lwod1, lwod0, 1, 1,
"B");
738 TH2F *mwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2TSmeanA225_HF");
739 TH2F *mwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth2_HF");
743 TH2F *C7fz225 = (TH2F *)mwod1->Clone(
"C7fz225");
744 C7fz225->Divide(mwod1, mwod0, 1, 1,
"B");
750 C7fz225->SetMarkerStyle(20);
751 C7fz225->SetMarkerSize(0.4);
752 C7fz225->GetZaxis()->SetLabelSize(0.08);
753 C7fz225->SetXTitle(
"#eta \b");
754 C7fz225->SetYTitle(
"#phi \b");
755 C7fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth2 \b");
756 C7fz225->SetMarkerColor(2);
757 C7fz225->SetLineColor(2);
758 C7fz225->SetMaximum(1.000);
759 C7fz225->SetMinimum(0.0001);
760 C7fz225->Draw(
"COLZ");
763 TH1F *aaaaaa7 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HF");
765 aaaaaa7->SetMarkerStyle(20);
766 aaaaaa7->SetMarkerSize(0.8);
767 aaaaaa7->GetYaxis()->SetLabelSize(0.04);
768 aaaaaa7->SetXTitle(
"TSmeanA in each event & cell HF \b");
769 aaaaaa7->SetMarkerColor(2);
770 aaaaaa7->SetLineColor(2);
775 TH2F *Diffe_Depth2_HF = (TH2F *)C7ff->Clone(
"Diffe_Depth2_HF");
776 for (
int i = 1;
i <= C7ff->GetXaxis()->GetNbins();
i++) {
777 for (
int j = 1;
j <= C7ff->GetYaxis()->GetNbins();
j++) {
778 double ccc1 = C7ff->GetBinContent(
i,
j);
779 Diffe_Depth2_HF->SetBinContent(
i,
j, 0.);
780 if (ccc1 < 0.5 || ccc1 > 1.5)
781 Diffe_Depth2_HF->SetBinContent(
i,
j, ccc1);
787 Diffe_Depth2_HF->SetMarkerStyle(20);
788 Diffe_Depth2_HF->SetMarkerSize(0.4);
789 Diffe_Depth2_HF->GetZaxis()->SetLabelSize(0.08);
791 Diffe_Depth2_HF->SetXTitle(
"#eta \b");
792 Diffe_Depth2_HF->SetYTitle(
"#phi \b");
793 Diffe_Depth2_HF->SetZTitle(
"<TSmeanA> out 0.5-1.5 -HF Depth2 \b");
794 Diffe_Depth2_HF->SetMarkerColor(2);
795 Diffe_Depth2_HF->SetLineColor(2);
796 Diffe_Depth2_HF->Draw(
"COLZ");
799 TH1F *diffTSmeanA_Depth2_HF =
new TH1F(
"diffTSmeanA_Depth2_HF",
"", 100, 0.0, 2.0);
800 for (
int i = 1;
i <= C7ff->GetXaxis()->GetNbins();
i++) {
801 for (
int j = 1;
j <= C7ff->GetYaxis()->GetNbins();
j++) {
802 if (C7ff->GetBinContent(
i,
j) != 0) {
803 double ccc1 = C7ff->GetBinContent(
i,
j);
804 diffTSmeanA_Depth2_HF->Fill(ccc1);
809 diffTSmeanA_Depth2_HF->SetMarkerStyle(20);
810 diffTSmeanA_Depth2_HF->SetMarkerSize(0.4);
811 diffTSmeanA_Depth2_HF->GetYaxis()->SetLabelSize(0.04);
812 diffTSmeanA_Depth2_HF->SetXTitle(
"<TSmeanA> in each cell \b");
813 diffTSmeanA_Depth2_HF->SetMarkerColor(2);
814 diffTSmeanA_Depth2_HF->SetLineColor(2);
815 diffTSmeanA_Depth2_HF->Draw(
"");
827 TH2F *nwod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth4TSmeanA_HO");
828 TH2F *nwod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth4_HO");
832 TH2F *C8ff = (TH2F *)nwod1->Clone(
"C8ff");
833 C8ff->Divide(nwod1, nwod0, 1, 1,
"B");
837 TH2F *owod1 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth4TSmeanA225_HO");
838 TH2F *owod0 = (TH2F *)
dir->FindObjectAny(
"h_mapDepth4_HO");
842 TH2F *C8fz225 = (TH2F *)owod1->Clone(
"C8fz225");
843 C8fz225->Divide(owod1, owod0, 1, 1,
"B");
849 C8fz225->SetMarkerStyle(20);
850 C8fz225->SetMarkerSize(0.4);
851 C8fz225->GetZaxis()->SetLabelSize(0.08);
852 C8fz225->SetXTitle(
"#eta \b");
853 C8fz225->SetYTitle(
"#phi \b");
854 C8fz225->SetZTitle(
"Rate for TSmeanA in each event & cell out 2.3-5.5 - HO Depth4 \b");
855 C8fz225->SetMarkerColor(2);
856 C8fz225->SetLineColor(2);
857 C8fz225->SetMaximum(1.000);
858 C8fz225->SetMinimum(0.000005);
859 C8fz225->Draw(
"COLZ");
862 TH1F *aaaaaa8 = (TH1F *)
dir->FindObjectAny(
"h_TSmeanA_HO");
864 aaaaaa8->SetMarkerStyle(20);
865 aaaaaa8->SetMarkerSize(0.8);
866 aaaaaa8->GetYaxis()->SetLabelSize(0.04);
867 aaaaaa8->SetXTitle(
"TSmeanA in each event & cell HO \b");
868 aaaaaa8->SetMarkerColor(2);
869 aaaaaa8->SetLineColor(2);
874 TH2F *Diffe_Depth4_HO = (TH2F *)C8ff->Clone(
"Diffe_Depth4_HO");
875 for (
int i = 1;
i <= C8ff->GetXaxis()->GetNbins();
i++) {
876 for (
int j = 1;
j <= C8ff->GetYaxis()->GetNbins();
j++) {
877 double ccc1 = C8ff->GetBinContent(
i,
j);
878 Diffe_Depth4_HO->SetBinContent(
i,
j, 0.);
879 if (ccc1 < 4.0 || ccc1 > 5.0)
880 Diffe_Depth4_HO->SetBinContent(
i,
j, ccc1);
886 Diffe_Depth4_HO->SetMarkerStyle(20);
887 Diffe_Depth4_HO->SetMarkerSize(0.4);
888 Diffe_Depth4_HO->GetZaxis()->SetLabelSize(0.08);
890 Diffe_Depth4_HO->SetXTitle(
"#eta \b");
891 Diffe_Depth4_HO->SetYTitle(
"#phi \b");
892 Diffe_Depth4_HO->SetZTitle(
"<TSmeanA> out 4.0-5.0 - HO Depth4 \b");
893 Diffe_Depth4_HO->SetMarkerColor(2);
894 Diffe_Depth4_HO->SetLineColor(2);
895 Diffe_Depth4_HO->Draw(
"COLZ");
898 TH1F *diffTSmeanA_Depth4_HO =
new TH1F(
"diffTSmeanA_Depth4_HO",
"", 100, 4.0, 5.0);
899 for (
int i = 1;
i <= C8ff->GetXaxis()->GetNbins();
i++) {
900 for (
int j = 1;
j <= C8ff->GetYaxis()->GetNbins();
j++) {
901 if (C8ff->GetBinContent(
i,
j) != 0) {
902 double ccc1 = C8ff->GetBinContent(
i,
j);
903 diffTSmeanA_Depth4_HO->Fill(ccc1);
908 diffTSmeanA_Depth4_HO->SetMarkerStyle(20);
909 diffTSmeanA_Depth4_HO->SetMarkerSize(0.4);
910 diffTSmeanA_Depth4_HO->GetYaxis()->SetLabelSize(0.04);
911 diffTSmeanA_Depth4_HO->SetXTitle(
"<TSmeanA> in each cell \b");
912 diffTSmeanA_Depth4_HO->SetMarkerColor(2);
913 diffTSmeanA_Depth4_HO->SetLineColor(2);
914 diffTSmeanA_Depth4_HO->Draw(
"");
int main(int argc, char *argv[])