Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 10:52:17

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 //inline void HERE(const char *msg) { std::cout << msg << std::endl; }
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);  //  no statistics _or_
0037                           //                  gStyle->SetOptStat(11111111);
0038                           //gStyle->SetOptStat(1101);// mame mean and rms
0039                           //    gStyle->SetOptStat(0101);// name and entries
0040                           //    gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0041                           //       gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0042                           //                gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0043   //    gStyle->SetOptStat(101110);                                          // entries, meam, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0044   //
0045   //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
0046   //    gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
0047   gStyle->SetOptFit(0010);  // constant, mean and sigma only !!
0048   //    gStyle->SetOptFit(00011);// constant, mean and sigma only !!
0049   // gStyle->SetOptFit(1101);
0050   //       gStyle->SetOptFit(1011);
0051   //
0052   //gStyle->SetStatX(0.98);
0053   //gStyle->SetStatY(0.99);
0054   //gStyle->SetStatW(0.30);
0055   //gStyle->SetStatH(0.25);
0056   gStyle->SetStatX(0.91);
0057   gStyle->SetStatY(0.75);
0058   gStyle->SetStatW(0.20);
0059   gStyle->SetStatH(0.10);
0060   // gStyle->SetStatH(0.35);
0061   //
0062 
0063   //        Float_t LeftOffset = 0.12;
0064   //Float_t TopOffset = 0.12;
0065   Float_t LeftOffset = 0.12;
0066   Float_t TopOffset = 0.22;
0067 
0068   gStyle->SetLineWidth(1);
0069   gStyle->SetErrorX(0);
0070 
0071   //---=[ Titles,Labels ]=-----------
0072   gStyle->SetOptTitle(0);  // title on/off
0073   //      gStyle->SetTitleColor(0);           // title color
0074   gStyle->SetTitleColor(1);  // title color
0075   //      gStyle->SetTitleX(0.35);            // title x-position
0076   gStyle->SetTitleX(0.15);  // title x-position
0077   gStyle->SetTitleH(0.15);  // title height
0078   //      gStyle->SetTitleW(0.53);            // title width
0079   gStyle->SetTitleW(0.60);         // title width
0080   gStyle->SetTitleFont(42);        // title font
0081   gStyle->SetTitleFontSize(0.07);  // title font size
0082 
0083   gStyle->SetPalette(1);
0084   //---=[ Histogram style ]=----------
0085   //        gStyle->SetHistFillColor(38);
0086   //            gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
0087   //    gStyle->SetFrameFillColor(5);// jeltyi
0088   //    gStyle->SetFrameFillColor(17);// seryi
0089   //    gStyle->SetFrameFillColor(18);// svetlo seryi
0090   //    gStyle->SetFrameFillColor(20);// svetlo kori4nev
0091   //        gStyle->SetFrameFillColor(33);// sine-seryi
0092   //    gStyle->SetFrameFillColor(40);// fiolet-seryi
0093   //    gStyle->SetFrameFillColor(23);// sv.kor
0094 
0095   //---=[ Pad style ]=----------------
0096   gStyle->SetPadTopMargin(TopOffset);
0097   gStyle->SetPadBottomMargin(LeftOffset);
0098   gStyle->SetPadRightMargin(TopOffset);
0099   gStyle->SetPadLeftMargin(LeftOffset);
0100   //---=[ SetCanvasDef           ]=----------------
0101   //======================================================================
0102   //
0103   // Connect the input file and get the 2-d histogram in memory
0104   //======================================================================
0105   //  TBrowser *b = new TBrowser
0106 
0107   //    TFile *hfile1= new TFile("/afs/cern.ch/user/z/zhokin/public/html/RMT/lastMarch2014/cmt/Global_229713.root", "READ");
0108   //    TFile *hfile1= new TFile("test.root", "READ");
0109   //
0110 
0111   TFile *hfile1 = new TFile("Global_362596.root", "READ");
0112   //TFile *hfile1 = new TFile("Global_362365.root", "READ");
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   //=============================================================================================== 1
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   //    if(twod0->IsA()->InheritsFrom("TH2F")){
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   //   if(twop0->IsA()->InheritsFrom("TH2F")){
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   //    Cefz225->GetYaxis()->SetLabelSize(0.04);
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       //      Diffe_Depth1_HB->SetBinContent(i,j,ccc1);
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   //    Diffe_Depth1_HB->SetTitle("any Error, HB Depth1 \n");
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   //=============================================================================================== 2
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   //   if(awod0->IsA()->InheritsFrom("TH2F")){
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   //   if(bwod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth2_HB->SetTitle("any Error, HB Depth2 \n");
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   //=============================================================================================== 3
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   //  if(cwod0->IsA()->InheritsFrom("TH2F")){
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   //   if(dwod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth1_HE->SetTitle("any Error, HE Depth1 \n");
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   //=============================================================================================== 4
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   //  if(ewod0->IsA()->InheritsFrom("TH2F")){
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   //  if(fwod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth2_HE->SetTitle("any Error, HE Depth2 \n");
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   //=============================================================================================== 5
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   //  if(gwod0->IsA()->InheritsFrom("TH2F")){
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   // if(jwod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth3_HE->SetTitle("any Error, HE Depth3 \n");
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   //=============================================================================================== 6
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   // if(iwod0->IsA()->InheritsFrom("TH2F")){
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   //  if(kwod0->IsA()->InheritsFrom("TH2F")){
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   //    C6fz225->GetYaxis()->SetLabelSize(0.04);
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   //    Diffe_Depth1_HF->SetTitle("any Error, HF Depth1 \n");
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   //=============================================================================================== 7
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   //  if(lwod0->IsA()->InheritsFrom("TH2F")){
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   //  if(mwod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth2_HF->SetTitle("any Error, HF Depth2 \n");
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   //=============================================================================================== 8
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   //  if(nwod0->IsA()->InheritsFrom("TH2F")){
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   //   if(owod0->IsA()->InheritsFrom("TH2F")){
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   //    Diffe_Depth4_HO->SetTitle("any Error, HO Depth4 \n");
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   //======================================================================================== end
0920   //======================================================================
0921   //======================================================================
0922   // close and delete all possible things:
0923 
0924   //   psfile->Close();
0925   psfile.Close();
0926 
0927   hfile1->Close();
0928   //    hfile1->Delete();
0929   hfile1->Close();
0930   //    hfile1->Delete();
0931 
0932   //  Exit Root
0933   gSystem->Exit(0);
0934   //======================================================================
0935 }