File indexing completed on 2024-04-06 12:05:51
0001 #include <string>
0002 #include <iomanip>
0003
0004 #include <iostream>
0005 #include <fstream>
0006 #include "TH1.h"
0007 #include "TH2.h"
0008 #include "TCanvas.h"
0009 #include "TROOT.h"
0010 #include <TMath.h>
0011 #include "TStyle.h"
0012 #include "TSystem.h"
0013 #include "TLegend.h"
0014 #include "TText.h"
0015 #include "TAxis.h"
0016 #include "TFile.h"
0017 #include "TLine.h"
0018 #include "TGraph.h"
0019 #include <TPaveText.h>
0020 #include <TPostScript.h>
0021 #include <TClass.h>
0022
0023
0024 using namespace std;
0025
0026
0027
0028
0029
0030 int main(int argc, char *argv[]) {
0031 std::string dirnm = "Analyzer";
0032
0033 printf("reco: gROOT Reset \n");
0034 gROOT->Reset();
0035 gROOT->SetStyle("Plain");
0036 gStyle->SetOptStat(0);
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047 gStyle->SetOptFit(0010);
0048
0049
0050
0051
0052
0053
0054
0055
0056 gStyle->SetStatX(0.91);
0057 gStyle->SetStatY(0.75);
0058 gStyle->SetStatW(0.20);
0059 gStyle->SetStatH(0.10);
0060
0061
0062
0063
0064
0065 Float_t LeftOffset = 0.12;
0066 Float_t TopOffset = 0.22;
0067
0068 gStyle->SetLineWidth(1);
0069 gStyle->SetErrorX(0);
0070
0071
0072 gStyle->SetOptTitle(0);
0073
0074 gStyle->SetTitleColor(1);
0075
0076 gStyle->SetTitleX(0.15);
0077 gStyle->SetTitleH(0.15);
0078
0079 gStyle->SetTitleW(0.60);
0080 gStyle->SetTitleFont(42);
0081 gStyle->SetTitleFontSize(0.07);
0082
0083 gStyle->SetPalette(1);
0084
0085
0086
0087
0088
0089
0090
0091
0092
0093
0094
0095
0096 gStyle->SetPadTopMargin(TopOffset);
0097 gStyle->SetPadBottomMargin(LeftOffset);
0098 gStyle->SetPadRightMargin(TopOffset);
0099 gStyle->SetPadLeftMargin(LeftOffset);
0100
0101
0102
0103
0104
0105
0106
0107
0108
0109
0110
0111 TFile *hfile1 = new TFile("Global_362596.root", "READ");
0112
0113
0114 TPostScript psfile("ztsmean.ps", 111);
0115
0116
0117
0118 TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
0119
0120 hfile1->ls();
0121 TDirectory *dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
0122
0123
0124
0125
0126
0127
0128 c1->Clear();
0129 c1->Divide(2, 2);
0130
0131 TH2F *twod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HB");
0132 TH2F *twod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
0133 twod1->Sumw2();
0134 twod0->Sumw2();
0135
0136 TH2F *Ceff = (TH2F *)twod1->Clone("Ceff");
0137 Ceff->Divide(twod1, twod0, 1, 1, "B");
0138 Ceff->Sumw2();
0139
0140 c1->cd(1);
0141 TH2F *twop1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HB");
0142 TH2F *twop0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
0143 twop1->Sumw2();
0144 twop0->Sumw2();
0145
0146 TH2F *Cefz225 = (TH2F *)twop1->Clone("Cefz225");
0147 Cefz225->Divide(twop1, twop0, 1, 1, "B");
0148 Cefz225->Sumw2();
0149
0150 gPad->SetGridy();
0151 gPad->SetGridx();
0152 gPad->SetLogz();
0153 Cefz225->SetMarkerStyle(20);
0154 Cefz225->SetMarkerSize(0.4);
0155
0156 Cefz225->GetZaxis()->SetLabelSize(0.08);
0157 Cefz225->SetXTitle("#eta \b");
0158 Cefz225->SetYTitle("#phi \b");
0159 Cefz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth1 \b");
0160 Cefz225->SetMarkerColor(2);
0161 Cefz225->SetLineColor(2);
0162 Cefz225->SetMaximum(1.000);
0163 Cefz225->SetMinimum(0.0001);
0164 Cefz225->Draw("COLZ");
0165
0166 c1->cd(2);
0167 TH1F *aaaaaa1 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HB");
0168 gPad->SetLogy();
0169 aaaaaa1->SetMarkerStyle(20);
0170 aaaaaa1->SetMarkerSize(0.8);
0171 aaaaaa1->GetYaxis()->SetLabelSize(0.04);
0172 aaaaaa1->SetXTitle("TSmeanA in each event & cell HB \b");
0173 aaaaaa1->SetMarkerColor(2);
0174 aaaaaa1->SetLineColor(2);
0175 aaaaaa1->Draw("");
0176
0177 c1->cd(3);
0178
0179 TH2F *Diffe_Depth1_HB = (TH2F *)Ceff->Clone("Diffe_Depth1_HB");
0180 for (int i = 1; i <= Ceff->GetXaxis()->GetNbins(); i++) {
0181 for (int j = 1; j <= Ceff->GetYaxis()->GetNbins(); j++) {
0182 double ccc1 = Ceff->GetBinContent(i, j);
0183 Diffe_Depth1_HB->SetBinContent(i, j, 0.);
0184
0185 if (ccc1 < 2.5 || ccc1 > 4.5)
0186 Diffe_Depth1_HB->SetBinContent(i, j, ccc1);
0187 }
0188 }
0189 gPad->SetGridy();
0190 gPad->SetGridx();
0191 gPad->SetLogz();
0192 Diffe_Depth1_HB->SetMarkerStyle(20);
0193 Diffe_Depth1_HB->SetMarkerSize(0.4);
0194 Diffe_Depth1_HB->GetZaxis()->SetLabelSize(0.08);
0195
0196 Diffe_Depth1_HB->SetXTitle("#eta \b");
0197 Diffe_Depth1_HB->SetYTitle("#phi \b");
0198 Diffe_Depth1_HB->SetZTitle("<TSmeanA> out 2.5-4.5- HB Depth1 \b");
0199 Diffe_Depth1_HB->SetMarkerColor(2);
0200 Diffe_Depth1_HB->SetLineColor(2);
0201 Diffe_Depth1_HB->Draw("COLZ");
0202
0203 c1->cd(4);
0204 TH1F *diffTSmeanA_Depth1_HB = new TH1F("diffTSmeanA_Depth1_HB", "", 100, 1.0, 6.0);
0205 for (int i = 1; i <= Ceff->GetXaxis()->GetNbins(); i++) {
0206 for (int j = 1; j <= Ceff->GetYaxis()->GetNbins(); j++) {
0207 if (Ceff->GetBinContent(i, j) != 0) {
0208 double ccc1 = Ceff->GetBinContent(i, j);
0209 diffTSmeanA_Depth1_HB->Fill(ccc1);
0210 }
0211 }
0212 }
0213 gPad->SetLogy();
0214 diffTSmeanA_Depth1_HB->SetMarkerStyle(20);
0215 diffTSmeanA_Depth1_HB->SetMarkerSize(0.4);
0216 diffTSmeanA_Depth1_HB->GetYaxis()->SetLabelSize(0.04);
0217 diffTSmeanA_Depth1_HB->SetXTitle("<TSmeanA> in each cell \b");
0218 diffTSmeanA_Depth1_HB->SetMarkerColor(2);
0219 diffTSmeanA_Depth1_HB->SetLineColor(2);
0220 diffTSmeanA_Depth1_HB->Draw("");
0221
0222 c1->Update();
0223
0224
0225
0226
0227
0228
0229 c1->Clear();
0230 c1->Divide(2, 2);
0231
0232 TH2F *awod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HB");
0233 TH2F *awod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
0234 awod1->Sumw2();
0235 awod0->Sumw2();
0236
0237 TH2F *C2ff = (TH2F *)awod1->Clone("C2ff");
0238 C2ff->Divide(awod1, awod0, 1, 1, "B");
0239 C2ff->Sumw2();
0240
0241 c1->cd(1);
0242 TH2F *bwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HB");
0243 TH2F *bwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
0244 bwod1->Sumw2();
0245 bwod0->Sumw2();
0246
0247 TH2F *C2fz225 = (TH2F *)bwod1->Clone("C2fz225");
0248 C2fz225->Divide(bwod1, bwod0, 1, 1, "B");
0249 C2fz225->Sumw2();
0250
0251 gPad->SetGridy();
0252 gPad->SetGridx();
0253 gPad->SetLogz();
0254 C2fz225->SetMarkerStyle(20);
0255 C2fz225->SetMarkerSize(0.4);
0256 C2fz225->GetZaxis()->SetLabelSize(0.08);
0257 C2fz225->SetXTitle("#eta \b");
0258 C2fz225->SetYTitle("#phi \b");
0259 C2fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth2 \b");
0260 C2fz225->SetMarkerColor(2);
0261 C2fz225->SetLineColor(2);
0262 C2fz225->SetMaximum(1.000);
0263 C2fz225->SetMinimum(0.0001);
0264 C2fz225->Draw("COLZ");
0265
0266 c1->cd(2);
0267 TH1F *aaaaaa2 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HB");
0268 gPad->SetLogy();
0269 aaaaaa2->SetMarkerStyle(20);
0270 aaaaaa2->SetMarkerSize(0.8);
0271 aaaaaa2->GetYaxis()->SetLabelSize(0.04);
0272 aaaaaa2->SetXTitle("TSmeanA in each event & cell HB \b");
0273 aaaaaa2->SetMarkerColor(2);
0274 aaaaaa2->SetLineColor(2);
0275 aaaaaa2->Draw("");
0276
0277 c1->cd(3);
0278
0279 TH2F *Diffe_Depth2_HB = (TH2F *)C2ff->Clone("Diffe_Depth2_HB");
0280 for (int i = 1; i <= C2ff->GetXaxis()->GetNbins(); i++) {
0281 for (int j = 1; j <= C2ff->GetYaxis()->GetNbins(); j++) {
0282 double ccc1 = C2ff->GetBinContent(i, j);
0283 Diffe_Depth2_HB->SetBinContent(i, j, 0.);
0284 if (ccc1 < 2.4 || ccc1 > 4.4)
0285 Diffe_Depth2_HB->SetBinContent(i, j, ccc1);
0286 }
0287 }
0288 gPad->SetGridy();
0289 gPad->SetGridx();
0290 gPad->SetLogz();
0291 Diffe_Depth2_HB->SetMarkerStyle(20);
0292 Diffe_Depth2_HB->SetMarkerSize(0.4);
0293 Diffe_Depth2_HB->GetZaxis()->SetLabelSize(0.08);
0294
0295 Diffe_Depth2_HB->SetXTitle("#eta \b");
0296 Diffe_Depth2_HB->SetYTitle("#phi \b");
0297 Diffe_Depth2_HB->SetZTitle("<TSmeanA> out 2.4-4.4 - HB Depth2 \b");
0298 Diffe_Depth2_HB->SetMarkerColor(2);
0299 Diffe_Depth2_HB->SetLineColor(2);
0300 Diffe_Depth2_HB->Draw("COLZ");
0301
0302 c1->cd(4);
0303 TH1F *diffTSmeanA_Depth2_HB = new TH1F("diffTSmeanA_Depth2_HB", "", 100, 1.0, 6.0);
0304 for (int i = 1; i <= C2ff->GetXaxis()->GetNbins(); i++) {
0305 for (int j = 1; j <= C2ff->GetYaxis()->GetNbins(); j++) {
0306 if (C2ff->GetBinContent(i, j) != 0) {
0307 double ccc1 = C2ff->GetBinContent(i, j);
0308 diffTSmeanA_Depth2_HB->Fill(ccc1);
0309 }
0310 }
0311 }
0312 gPad->SetLogy();
0313 diffTSmeanA_Depth2_HB->SetMarkerStyle(20);
0314 diffTSmeanA_Depth2_HB->SetMarkerSize(0.4);
0315 diffTSmeanA_Depth2_HB->GetYaxis()->SetLabelSize(0.04);
0316 diffTSmeanA_Depth2_HB->SetXTitle("<TSmeanA> in each cell \b");
0317 diffTSmeanA_Depth2_HB->SetMarkerColor(2);
0318 diffTSmeanA_Depth2_HB->SetLineColor(2);
0319 diffTSmeanA_Depth2_HB->Draw("");
0320
0321 c1->Update();
0322
0323
0324
0325
0326
0327
0328 c1->Clear();
0329 c1->Divide(2, 2);
0330
0331 TH2F *cwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HE");
0332 TH2F *cwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
0333 cwod1->Sumw2();
0334 cwod0->Sumw2();
0335
0336 TH2F *C3ff = (TH2F *)cwod1->Clone("C3ff");
0337 C3ff->Divide(cwod1, cwod0, 1, 1, "B");
0338 C3ff->Sumw2();
0339
0340 c1->cd(1);
0341 TH2F *dwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HE");
0342 TH2F *dwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
0343 dwod1->Sumw2();
0344 dwod0->Sumw2();
0345
0346 TH2F *C3fz225 = (TH2F *)dwod1->Clone("C3fz225");
0347 C3fz225->Divide(dwod1, dwod0, 1, 1, "B");
0348 C3fz225->Sumw2();
0349
0350 gPad->SetGridy();
0351 gPad->SetGridx();
0352 gPad->SetLogz();
0353 C3fz225->SetMarkerStyle(20);
0354 C3fz225->SetMarkerSize(0.4);
0355 C3fz225->GetZaxis()->SetLabelSize(0.08);
0356 C3fz225->SetXTitle("#eta \b");
0357 C3fz225->SetYTitle("#phi \b");
0358 C3fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth1 \b");
0359 C3fz225->SetMarkerColor(2);
0360 C3fz225->SetLineColor(2);
0361 C3fz225->SetMaximum(1.000);
0362 C3fz225->SetMinimum(0.0001);
0363 C3fz225->Draw("COLZ");
0364
0365 c1->cd(2);
0366 TH1F *aaaaaa3 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
0367 gPad->SetLogy();
0368 aaaaaa3->SetMarkerStyle(20);
0369 aaaaaa3->SetMarkerSize(0.8);
0370 aaaaaa3->GetYaxis()->SetLabelSize(0.04);
0371 aaaaaa3->SetXTitle("TSmeanA in each event & cell HE \b");
0372 aaaaaa3->SetMarkerColor(2);
0373 aaaaaa3->SetLineColor(2);
0374 aaaaaa3->Draw("");
0375
0376 c1->cd(3);
0377
0378 TH2F *Diffe_Depth1_HE = (TH2F *)C3ff->Clone("Diffe_Depth1_HE");
0379 for (int i = 1; i <= C3ff->GetXaxis()->GetNbins(); i++) {
0380 for (int j = 1; j <= C3ff->GetYaxis()->GetNbins(); j++) {
0381 double ccc1 = C3ff->GetBinContent(i, j);
0382 Diffe_Depth1_HE->SetBinContent(i, j, 0.);
0383 if (ccc1 < 2.4 || ccc1 > 4.6)
0384 Diffe_Depth1_HE->SetBinContent(i, j, ccc1);
0385 }
0386 }
0387 gPad->SetGridy();
0388 gPad->SetGridx();
0389 gPad->SetLogz();
0390 Diffe_Depth1_HE->SetMarkerStyle(20);
0391 Diffe_Depth1_HE->SetMarkerSize(0.4);
0392 Diffe_Depth1_HE->GetZaxis()->SetLabelSize(0.08);
0393
0394 Diffe_Depth1_HE->SetXTitle("#eta \b");
0395 Diffe_Depth1_HE->SetYTitle("#phi \b");
0396 Diffe_Depth1_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth1 \b");
0397 Diffe_Depth1_HE->SetMarkerColor(2);
0398 Diffe_Depth1_HE->SetLineColor(2);
0399 Diffe_Depth1_HE->Draw("COLZ");
0400
0401 c1->cd(4);
0402 TH1F *diffTSmeanA_Depth1_HE = new TH1F("diffTSmeanA_Depth1_HE", "", 100, 0.0, 7.0);
0403 for (int i = 1; i <= C3ff->GetXaxis()->GetNbins(); i++) {
0404 for (int j = 1; j <= C3ff->GetYaxis()->GetNbins(); j++) {
0405 if (C3ff->GetBinContent(i, j) != 0) {
0406 double ccc1 = C3ff->GetBinContent(i, j);
0407 diffTSmeanA_Depth1_HE->Fill(ccc1);
0408 }
0409 }
0410 }
0411 gPad->SetLogy();
0412 diffTSmeanA_Depth1_HE->SetMarkerStyle(20);
0413 diffTSmeanA_Depth1_HE->SetMarkerSize(0.4);
0414 diffTSmeanA_Depth1_HE->GetYaxis()->SetLabelSize(0.04);
0415 diffTSmeanA_Depth1_HE->SetXTitle("<TSmeanA> in each cell \b");
0416 diffTSmeanA_Depth1_HE->SetMarkerColor(2);
0417 diffTSmeanA_Depth1_HE->SetLineColor(2);
0418 diffTSmeanA_Depth1_HE->Draw("");
0419
0420 c1->Update();
0421
0422
0423
0424
0425
0426
0427 c1->Clear();
0428 c1->Divide(2, 2);
0429
0430 TH2F *ewod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HE");
0431 TH2F *ewod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
0432 ewod1->Sumw2();
0433 ewod0->Sumw2();
0434
0435 TH2F *C4ff = (TH2F *)ewod1->Clone("C4ff");
0436 C4ff->Divide(ewod1, ewod0, 1, 1, "B");
0437 C4ff->Sumw2();
0438
0439 c1->cd(1);
0440 TH2F *fwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HE");
0441 TH2F *fwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
0442 fwod1->Sumw2();
0443 fwod0->Sumw2();
0444
0445 TH2F *C4fz225 = (TH2F *)fwod1->Clone("C4fz225");
0446 C4fz225->Divide(fwod1, fwod0, 1, 1, "B");
0447 C4fz225->Sumw2();
0448
0449 gPad->SetGridy();
0450 gPad->SetGridx();
0451 gPad->SetLogz();
0452 C4fz225->SetMarkerStyle(20);
0453 C4fz225->SetMarkerSize(0.4);
0454 C4fz225->GetZaxis()->SetLabelSize(0.08);
0455 C4fz225->SetXTitle("#eta \b");
0456 C4fz225->SetYTitle("#phi \b");
0457 C4fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth2 \b");
0458 C4fz225->SetMarkerColor(2);
0459 C4fz225->SetLineColor(2);
0460 C4fz225->SetMaximum(1.000);
0461 C4fz225->SetMinimum(0.0001);
0462 C4fz225->Draw("COLZ");
0463
0464 c1->cd(2);
0465 TH1F *aaaaaa4 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
0466 gPad->SetLogy();
0467 aaaaaa4->SetMarkerStyle(20);
0468 aaaaaa4->SetMarkerSize(0.8);
0469 aaaaaa4->GetYaxis()->SetLabelSize(0.04);
0470 aaaaaa4->SetXTitle("TSmeanA in each event & cell HE \b");
0471 aaaaaa4->SetMarkerColor(2);
0472 aaaaaa4->SetLineColor(2);
0473 aaaaaa4->Draw("");
0474
0475 c1->cd(3);
0476
0477 TH2F *Diffe_Depth2_HE = (TH2F *)C4ff->Clone("Diffe_Depth2_HE");
0478 for (int i = 1; i <= C4ff->GetXaxis()->GetNbins(); i++) {
0479 for (int j = 1; j <= C4ff->GetYaxis()->GetNbins(); j++) {
0480 double ccc1 = C4ff->GetBinContent(i, j);
0481 Diffe_Depth2_HE->SetBinContent(i, j, 0.);
0482 if (ccc1 < 2.4 || ccc1 > 4.6)
0483 Diffe_Depth2_HE->SetBinContent(i, j, ccc1);
0484 }
0485 }
0486 gPad->SetGridy();
0487 gPad->SetGridx();
0488 gPad->SetLogz();
0489 Diffe_Depth2_HE->SetMarkerStyle(20);
0490 Diffe_Depth2_HE->SetMarkerSize(0.4);
0491 Diffe_Depth2_HE->GetZaxis()->SetLabelSize(0.08);
0492
0493 Diffe_Depth2_HE->SetXTitle("#eta \b");
0494 Diffe_Depth2_HE->SetYTitle("#phi \b");
0495 Diffe_Depth2_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth2 \b");
0496 Diffe_Depth2_HE->SetMarkerColor(2);
0497 Diffe_Depth2_HE->SetLineColor(2);
0498 Diffe_Depth2_HE->Draw("COLZ");
0499
0500 c1->cd(4);
0501 TH1F *diffTSmeanA_Depth2_HE = new TH1F("diffTSmeanA_Depth2_HE", "", 100, 0.0, 7.0);
0502 for (int i = 1; i <= C4ff->GetXaxis()->GetNbins(); i++) {
0503 for (int j = 1; j <= C4ff->GetYaxis()->GetNbins(); j++) {
0504 if (C4ff->GetBinContent(i, j) != 0) {
0505 double ccc1 = C4ff->GetBinContent(i, j);
0506 diffTSmeanA_Depth2_HE->Fill(ccc1);
0507 }
0508 }
0509 }
0510 gPad->SetLogy();
0511 diffTSmeanA_Depth2_HE->SetMarkerStyle(20);
0512 diffTSmeanA_Depth2_HE->SetMarkerSize(0.4);
0513 diffTSmeanA_Depth2_HE->GetYaxis()->SetLabelSize(0.04);
0514 diffTSmeanA_Depth2_HE->SetXTitle("<TSmeanA> in each cell \b");
0515 diffTSmeanA_Depth2_HE->SetMarkerColor(2);
0516 diffTSmeanA_Depth2_HE->SetLineColor(2);
0517 diffTSmeanA_Depth2_HE->Draw("");
0518
0519 c1->Update();
0520
0521
0522
0523
0524
0525
0526 c1->Clear();
0527 c1->Divide(2, 2);
0528
0529 TH2F *gwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HE");
0530 TH2F *gwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
0531 gwod1->Sumw2();
0532 gwod0->Sumw2();
0533
0534 TH2F *C5ff = (TH2F *)gwod1->Clone("C5ff");
0535 C5ff->Divide(gwod1, gwod0, 1, 1, "B");
0536 C5ff->Sumw2();
0537
0538 c1->cd(1);
0539 TH2F *jwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA225_HE");
0540 TH2F *jwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
0541 jwod1->Sumw2();
0542 jwod0->Sumw2();
0543
0544 TH2F *C5fz225 = (TH2F *)jwod1->Clone("C5fz225");
0545 C5fz225->Divide(jwod1, jwod0, 1, 1, "B");
0546 C5fz225->Sumw2();
0547
0548 gPad->SetGridy();
0549 gPad->SetGridx();
0550 gPad->SetLogz();
0551 C5fz225->SetMarkerStyle(20);
0552 C5fz225->SetMarkerSize(0.4);
0553 C5fz225->GetZaxis()->SetLabelSize(0.08);
0554 C5fz225->SetXTitle("#eta \b");
0555 C5fz225->SetYTitle("#phi \b");
0556 C5fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth3 \b");
0557 C5fz225->SetMarkerColor(2);
0558 C5fz225->SetLineColor(2);
0559 C5fz225->SetMaximum(1.000);
0560 C5fz225->SetMinimum(0.0001);
0561 C5fz225->Draw("COLZ");
0562
0563 c1->cd(2);
0564 TH1F *aaaaaa5 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
0565 gPad->SetLogy();
0566 aaaaaa5->SetMarkerStyle(20);
0567 aaaaaa5->SetMarkerSize(0.8);
0568 aaaaaa5->GetYaxis()->SetLabelSize(0.04);
0569 aaaaaa5->SetXTitle("TSmeanA in each event & cell HE \b");
0570 aaaaaa5->SetMarkerColor(2);
0571 aaaaaa5->SetLineColor(2);
0572 aaaaaa5->Draw("");
0573
0574 c1->cd(3);
0575
0576 TH2F *Diffe_Depth3_HE = (TH2F *)C5ff->Clone("Diffe_Depth3_HE");
0577 for (int i = 1; i <= C5ff->GetXaxis()->GetNbins(); i++) {
0578 for (int j = 1; j <= C5ff->GetYaxis()->GetNbins(); j++) {
0579 double ccc1 = C5ff->GetBinContent(i, j);
0580 Diffe_Depth3_HE->SetBinContent(i, j, 0.);
0581 if (ccc1 < 2.4 || ccc1 > 4.6)
0582 Diffe_Depth3_HE->SetBinContent(i, j, ccc1);
0583 }
0584 }
0585 gPad->SetGridy();
0586 gPad->SetGridx();
0587 gPad->SetLogz();
0588 Diffe_Depth3_HE->SetMarkerStyle(20);
0589 Diffe_Depth3_HE->SetMarkerSize(0.4);
0590 Diffe_Depth3_HE->GetZaxis()->SetLabelSize(0.08);
0591
0592 Diffe_Depth3_HE->SetXTitle("#eta \b");
0593 Diffe_Depth3_HE->SetYTitle("#phi \b");
0594 Diffe_Depth3_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth3 \b");
0595 Diffe_Depth3_HE->SetMarkerColor(2);
0596 Diffe_Depth3_HE->SetLineColor(2);
0597 Diffe_Depth3_HE->Draw("COLZ");
0598
0599 c1->cd(4);
0600 TH1F *diffTSmeanA_Depth3_HE = new TH1F("diffTSmeanA_Depth3_HE", "", 100, 0.0, 7.0);
0601 for (int i = 1; i <= C5ff->GetXaxis()->GetNbins(); i++) {
0602 for (int j = 1; j <= C5ff->GetYaxis()->GetNbins(); j++) {
0603 if (C5ff->GetBinContent(i, j) != 0) {
0604 double ccc1 = C5ff->GetBinContent(i, j);
0605 diffTSmeanA_Depth3_HE->Fill(ccc1);
0606 }
0607 }
0608 }
0609 gPad->SetLogy();
0610 diffTSmeanA_Depth3_HE->SetMarkerStyle(20);
0611 diffTSmeanA_Depth3_HE->SetMarkerSize(0.4);
0612 diffTSmeanA_Depth3_HE->GetYaxis()->SetLabelSize(0.04);
0613 diffTSmeanA_Depth3_HE->SetXTitle("<TSmeanA> in each cell \b");
0614 diffTSmeanA_Depth3_HE->SetMarkerColor(2);
0615 diffTSmeanA_Depth3_HE->SetLineColor(2);
0616 diffTSmeanA_Depth3_HE->Draw("");
0617
0618 c1->Update();
0619
0620
0621
0622
0623
0624
0625 c1->Clear();
0626 c1->Divide(2, 2);
0627
0628 TH2F *iwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HF");
0629 TH2F *iwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
0630 iwod1->Sumw2();
0631 iwod0->Sumw2();
0632
0633 TH2F *C6ff = (TH2F *)iwod1->Clone("C6ff");
0634 C6ff->Divide(iwod1, iwod0, 1, 1, "B");
0635 C6ff->Sumw2();
0636
0637 c1->cd(1);
0638 TH2F *kwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HF");
0639 TH2F *kwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
0640 kwod1->Sumw2();
0641 kwod0->Sumw2();
0642
0643 TH2F *C6fz225 = (TH2F *)kwod1->Clone("C6fz225");
0644 C6fz225->Divide(kwod1, kwod0, 1, 1, "B");
0645 C6fz225->Sumw2();
0646
0647 gPad->SetGridy();
0648 gPad->SetGridx();
0649 gPad->SetLogz();
0650 C6fz225->SetMarkerStyle(20);
0651 C6fz225->SetMarkerSize(0.4);
0652
0653 C6fz225->GetZaxis()->SetLabelSize(0.08);
0654 C6fz225->SetXTitle("#eta \b");
0655 C6fz225->SetYTitle("#phi \b");
0656 C6fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth1 \b");
0657 C6fz225->SetMarkerColor(2);
0658 C6fz225->SetLineColor(2);
0659 C6fz225->SetMaximum(1.000);
0660 C6fz225->SetMinimum(0.0001);
0661 C6fz225->Draw("COLZ");
0662
0663 c1->cd(2);
0664 TH1F *aaaaaa6 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HF");
0665 gPad->SetLogy();
0666 aaaaaa6->SetMarkerStyle(20);
0667 aaaaaa6->SetMarkerSize(0.8);
0668 aaaaaa6->GetYaxis()->SetLabelSize(0.04);
0669 aaaaaa6->SetXTitle("TSmeanA in each event & cell HF \b");
0670 aaaaaa6->SetMarkerColor(2);
0671 aaaaaa6->SetLineColor(2);
0672 aaaaaa6->Draw("");
0673
0674 c1->cd(3);
0675
0676 TH2F *Diffe_Depth1_HF = (TH2F *)C6ff->Clone("Diffe_Depth1_HF");
0677 for (int i = 1; i <= C6ff->GetXaxis()->GetNbins(); i++) {
0678 for (int j = 1; j <= C6ff->GetYaxis()->GetNbins(); j++) {
0679 double ccc1 = C6ff->GetBinContent(i, j);
0680 Diffe_Depth1_HF->SetBinContent(i, j, 0.);
0681 if (ccc1 < 0.5 || ccc1 > 1.5)
0682 Diffe_Depth1_HF->SetBinContent(i, j, ccc1);
0683 }
0684 }
0685 gPad->SetGridy();
0686 gPad->SetGridx();
0687 gPad->SetLogz();
0688 Diffe_Depth1_HF->SetMarkerStyle(20);
0689 Diffe_Depth1_HF->SetMarkerSize(0.4);
0690 Diffe_Depth1_HF->GetZaxis()->SetLabelSize(0.08);
0691
0692 Diffe_Depth1_HF->SetXTitle("#eta \b");
0693 Diffe_Depth1_HF->SetYTitle("#phi \b");
0694 Diffe_Depth1_HF->SetZTitle("<TSmeanA> out 0.5-1.5 - HF Depth1 \b");
0695 Diffe_Depth1_HF->SetMarkerColor(2);
0696 Diffe_Depth1_HF->SetLineColor(2);
0697 Diffe_Depth1_HF->Draw("COLZ");
0698
0699 c1->cd(4);
0700 TH1F *diffTSmeanA_Depth1_HF = new TH1F("diffTSmeanA_Depth1_HF", "", 100, 0.0, 2.0);
0701 for (int i = 1; i <= C6ff->GetXaxis()->GetNbins(); i++) {
0702 for (int j = 1; j <= C6ff->GetYaxis()->GetNbins(); j++) {
0703 if (C6ff->GetBinContent(i, j) != 0) {
0704 double ccc1 = C6ff->GetBinContent(i, j);
0705 diffTSmeanA_Depth1_HF->Fill(ccc1);
0706 }
0707 }
0708 }
0709 gPad->SetLogy();
0710 diffTSmeanA_Depth1_HF->SetMarkerStyle(20);
0711 diffTSmeanA_Depth1_HF->SetMarkerSize(0.4);
0712 diffTSmeanA_Depth1_HF->GetYaxis()->SetLabelSize(0.04);
0713 diffTSmeanA_Depth1_HF->SetXTitle("<TSmeanA> in each cell \b");
0714 diffTSmeanA_Depth1_HF->SetMarkerColor(2);
0715 diffTSmeanA_Depth1_HF->SetLineColor(2);
0716 diffTSmeanA_Depth1_HF->Draw("");
0717
0718 c1->Update();
0719
0720
0721
0722
0723
0724
0725 c1->Clear();
0726 c1->Divide(2, 2);
0727
0728 TH2F *lwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HF");
0729 TH2F *lwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
0730 lwod1->Sumw2();
0731 lwod0->Sumw2();
0732
0733 TH2F *C7ff = (TH2F *)lwod1->Clone("C7ff");
0734 C7ff->Divide(lwod1, lwod0, 1, 1, "B");
0735 C7ff->Sumw2();
0736
0737 c1->cd(1);
0738 TH2F *mwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HF");
0739 TH2F *mwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
0740 mwod1->Sumw2();
0741 mwod0->Sumw2();
0742
0743 TH2F *C7fz225 = (TH2F *)mwod1->Clone("C7fz225");
0744 C7fz225->Divide(mwod1, mwod0, 1, 1, "B");
0745 C7fz225->Sumw2();
0746
0747 gPad->SetGridy();
0748 gPad->SetGridx();
0749 gPad->SetLogz();
0750 C7fz225->SetMarkerStyle(20);
0751 C7fz225->SetMarkerSize(0.4);
0752 C7fz225->GetZaxis()->SetLabelSize(0.08);
0753 C7fz225->SetXTitle("#eta \b");
0754 C7fz225->SetYTitle("#phi \b");
0755 C7fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth2 \b");
0756 C7fz225->SetMarkerColor(2);
0757 C7fz225->SetLineColor(2);
0758 C7fz225->SetMaximum(1.000);
0759 C7fz225->SetMinimum(0.0001);
0760 C7fz225->Draw("COLZ");
0761
0762 c1->cd(2);
0763 TH1F *aaaaaa7 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HF");
0764 gPad->SetLogy();
0765 aaaaaa7->SetMarkerStyle(20);
0766 aaaaaa7->SetMarkerSize(0.8);
0767 aaaaaa7->GetYaxis()->SetLabelSize(0.04);
0768 aaaaaa7->SetXTitle("TSmeanA in each event & cell HF \b");
0769 aaaaaa7->SetMarkerColor(2);
0770 aaaaaa7->SetLineColor(2);
0771 aaaaaa7->Draw("");
0772
0773 c1->cd(3);
0774
0775 TH2F *Diffe_Depth2_HF = (TH2F *)C7ff->Clone("Diffe_Depth2_HF");
0776 for (int i = 1; i <= C7ff->GetXaxis()->GetNbins(); i++) {
0777 for (int j = 1; j <= C7ff->GetYaxis()->GetNbins(); j++) {
0778 double ccc1 = C7ff->GetBinContent(i, j);
0779 Diffe_Depth2_HF->SetBinContent(i, j, 0.);
0780 if (ccc1 < 0.5 || ccc1 > 1.5)
0781 Diffe_Depth2_HF->SetBinContent(i, j, ccc1);
0782 }
0783 }
0784 gPad->SetGridy();
0785 gPad->SetGridx();
0786 gPad->SetLogz();
0787 Diffe_Depth2_HF->SetMarkerStyle(20);
0788 Diffe_Depth2_HF->SetMarkerSize(0.4);
0789 Diffe_Depth2_HF->GetZaxis()->SetLabelSize(0.08);
0790
0791 Diffe_Depth2_HF->SetXTitle("#eta \b");
0792 Diffe_Depth2_HF->SetYTitle("#phi \b");
0793 Diffe_Depth2_HF->SetZTitle("<TSmeanA> out 0.5-1.5 -HF Depth2 \b");
0794 Diffe_Depth2_HF->SetMarkerColor(2);
0795 Diffe_Depth2_HF->SetLineColor(2);
0796 Diffe_Depth2_HF->Draw("COLZ");
0797
0798 c1->cd(4);
0799 TH1F *diffTSmeanA_Depth2_HF = new TH1F("diffTSmeanA_Depth2_HF", "", 100, 0.0, 2.0);
0800 for (int i = 1; i <= C7ff->GetXaxis()->GetNbins(); i++) {
0801 for (int j = 1; j <= C7ff->GetYaxis()->GetNbins(); j++) {
0802 if (C7ff->GetBinContent(i, j) != 0) {
0803 double ccc1 = C7ff->GetBinContent(i, j);
0804 diffTSmeanA_Depth2_HF->Fill(ccc1);
0805 }
0806 }
0807 }
0808 gPad->SetLogy();
0809 diffTSmeanA_Depth2_HF->SetMarkerStyle(20);
0810 diffTSmeanA_Depth2_HF->SetMarkerSize(0.4);
0811 diffTSmeanA_Depth2_HF->GetYaxis()->SetLabelSize(0.04);
0812 diffTSmeanA_Depth2_HF->SetXTitle("<TSmeanA> in each cell \b");
0813 diffTSmeanA_Depth2_HF->SetMarkerColor(2);
0814 diffTSmeanA_Depth2_HF->SetLineColor(2);
0815 diffTSmeanA_Depth2_HF->Draw("");
0816
0817 c1->Update();
0818
0819
0820
0821
0822
0823
0824 c1->Clear();
0825 c1->Divide(2, 2);
0826
0827 TH2F *nwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HO");
0828 TH2F *nwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
0829 nwod1->Sumw2();
0830 nwod0->Sumw2();
0831
0832 TH2F *C8ff = (TH2F *)nwod1->Clone("C8ff");
0833 C8ff->Divide(nwod1, nwod0, 1, 1, "B");
0834 C8ff->Sumw2();
0835
0836 c1->cd(1);
0837 TH2F *owod1 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HO");
0838 TH2F *owod0 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
0839 owod1->Sumw2();
0840 owod0->Sumw2();
0841
0842 TH2F *C8fz225 = (TH2F *)owod1->Clone("C8fz225");
0843 C8fz225->Divide(owod1, owod0, 1, 1, "B");
0844 C8fz225->Sumw2();
0845
0846 gPad->SetGridy();
0847 gPad->SetGridx();
0848 gPad->SetLogz();
0849 C8fz225->SetMarkerStyle(20);
0850 C8fz225->SetMarkerSize(0.4);
0851 C8fz225->GetZaxis()->SetLabelSize(0.08);
0852 C8fz225->SetXTitle("#eta \b");
0853 C8fz225->SetYTitle("#phi \b");
0854 C8fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5.5 - HO Depth4 \b");
0855 C8fz225->SetMarkerColor(2);
0856 C8fz225->SetLineColor(2);
0857 C8fz225->SetMaximum(1.000);
0858 C8fz225->SetMinimum(0.000005);
0859 C8fz225->Draw("COLZ");
0860
0861 c1->cd(2);
0862 TH1F *aaaaaa8 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HO");
0863 gPad->SetLogy();
0864 aaaaaa8->SetMarkerStyle(20);
0865 aaaaaa8->SetMarkerSize(0.8);
0866 aaaaaa8->GetYaxis()->SetLabelSize(0.04);
0867 aaaaaa8->SetXTitle("TSmeanA in each event & cell HO \b");
0868 aaaaaa8->SetMarkerColor(2);
0869 aaaaaa8->SetLineColor(2);
0870 aaaaaa8->Draw("");
0871
0872 c1->cd(3);
0873
0874 TH2F *Diffe_Depth4_HO = (TH2F *)C8ff->Clone("Diffe_Depth4_HO");
0875 for (int i = 1; i <= C8ff->GetXaxis()->GetNbins(); i++) {
0876 for (int j = 1; j <= C8ff->GetYaxis()->GetNbins(); j++) {
0877 double ccc1 = C8ff->GetBinContent(i, j);
0878 Diffe_Depth4_HO->SetBinContent(i, j, 0.);
0879 if (ccc1 < 4.0 || ccc1 > 5.0)
0880 Diffe_Depth4_HO->SetBinContent(i, j, ccc1);
0881 }
0882 }
0883 gPad->SetGridy();
0884 gPad->SetGridx();
0885 gPad->SetLogz();
0886 Diffe_Depth4_HO->SetMarkerStyle(20);
0887 Diffe_Depth4_HO->SetMarkerSize(0.4);
0888 Diffe_Depth4_HO->GetZaxis()->SetLabelSize(0.08);
0889
0890 Diffe_Depth4_HO->SetXTitle("#eta \b");
0891 Diffe_Depth4_HO->SetYTitle("#phi \b");
0892 Diffe_Depth4_HO->SetZTitle("<TSmeanA> out 4.0-5.0 - HO Depth4 \b");
0893 Diffe_Depth4_HO->SetMarkerColor(2);
0894 Diffe_Depth4_HO->SetLineColor(2);
0895 Diffe_Depth4_HO->Draw("COLZ");
0896
0897 c1->cd(4);
0898 TH1F *diffTSmeanA_Depth4_HO = new TH1F("diffTSmeanA_Depth4_HO", "", 100, 4.0, 5.0);
0899 for (int i = 1; i <= C8ff->GetXaxis()->GetNbins(); i++) {
0900 for (int j = 1; j <= C8ff->GetYaxis()->GetNbins(); j++) {
0901 if (C8ff->GetBinContent(i, j) != 0) {
0902 double ccc1 = C8ff->GetBinContent(i, j);
0903 diffTSmeanA_Depth4_HO->Fill(ccc1);
0904 }
0905 }
0906 }
0907 gPad->SetLogy();
0908 diffTSmeanA_Depth4_HO->SetMarkerStyle(20);
0909 diffTSmeanA_Depth4_HO->SetMarkerSize(0.4);
0910 diffTSmeanA_Depth4_HO->GetYaxis()->SetLabelSize(0.04);
0911 diffTSmeanA_Depth4_HO->SetXTitle("<TSmeanA> in each cell \b");
0912 diffTSmeanA_Depth4_HO->SetMarkerColor(2);
0913 diffTSmeanA_Depth4_HO->SetLineColor(2);
0914 diffTSmeanA_Depth4_HO->Draw("");
0915
0916 c1->Update();
0917
0918
0919
0920
0921
0922
0923
0924
0925 psfile.Close();
0926
0927 hfile1->Close();
0928
0929 hfile1->Close();
0930
0931
0932
0933 gSystem->Exit(0);
0934
0935 }