Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:08:26

0001 #include <iostream>
0002 #include <fstream>
0003 #include "TH1.h"
0004 #include "TH2.h"
0005 #include "TCanvas.h"
0006 #include "TROOT.h"
0007 #include <TMath.h>
0008 #include "TStyle.h"
0009 #include "TSystem.h"
0010 #include "TLegend.h"
0011 #include "TText.h"
0012 #include "TAxis.h"
0013 #include "TFile.h"
0014 #include "TLine.h"
0015 #include "TGraph.h"
0016 #include <TPaveText.h>
0017 #include <TPostScript.h>
0018 #include <TClass.h>
0019 
0020 //
0021 using namespace std;
0022 //
0023 //
0024 
0025 //inline void HERE(const char *msg) { std::cout << msg << std::endl; }
0026 
0027 int main(int argc, char *argv[]) {
0028   //======================================================================
0029   printf("reco: gROOT Reset \n");
0030   gROOT->Reset();
0031   gROOT->SetStyle("Plain");
0032   //        gStyle->SetOptStat(0);   //  no statistics _or_
0033   //                  gStyle->SetOptStat(11111111);
0034   //gStyle->SetOptStat(1101);// name mean and rms
0035   //    gStyle->SetOptStat(0101);// name and entries
0036   //       gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0037   gStyle->SetOptStat(1110000);  // und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0038   //    gStyle->SetOptStat(101110);                                          // entries, mean, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0039   //    gStyle->SetOptStat(100000);//  over !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0040   //
0041   //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
0042   //    gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
0043   //    gStyle->SetOptFit(0010);// constant, mean and sigma only !!
0044   //    gStyle->SetOptFit(00011);// constant, mean and sigma only !!
0045   // gStyle->SetOptFit(1101);
0046   //       gStyle->SetOptFit(1011);
0047   //
0048   //gStyle->SetStatX(0.98);
0049   //gStyle->SetStatY(0.99);
0050   //gStyle->SetStatW(0.30);
0051   //gStyle->SetStatH(0.25);
0052   gStyle->SetStatX(0.91);
0053   gStyle->SetStatY(0.75);
0054   gStyle->SetStatW(0.20);
0055   gStyle->SetStatH(0.10);
0056   // gStyle->SetStatH(0.35);
0057   //
0058 
0059   //        Float_t LeftOffset = 0.12;
0060   //Float_t TopOffset = 0.12;
0061   Float_t LeftOffset = 0.12;
0062   Float_t TopOffset = 0.22;
0063 
0064   gStyle->SetLineWidth(1);
0065   gStyle->SetErrorX(0);
0066 
0067   //---=[ Titles,Labels ]=-----------
0068   gStyle->SetOptTitle(0);  // title on/off
0069   //      gStyle->SetTitleColor(0);           // title color
0070   gStyle->SetTitleColor(1);  // title color
0071   //      gStyle->SetTitleX(0.35);            // title x-position
0072   gStyle->SetTitleX(0.15);  // title x-position
0073   gStyle->SetTitleH(0.15);  // title height
0074   //      gStyle->SetTitleW(0.53);            // title width
0075   gStyle->SetTitleW(0.60);         // title width
0076   gStyle->SetTitleFont(42);        // title font
0077   gStyle->SetTitleFontSize(0.07);  // title font size
0078 
0079   gStyle->SetPalette(1);
0080   //---=[ Histogram style ]=----------
0081   //        gStyle->SetHistFillColor(38);
0082   //            gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
0083   //    gStyle->SetFrameFillColor(5);// jeltyi
0084   //    gStyle->SetFrameFillColor(17);// seryi
0085   //    gStyle->SetFrameFillColor(18);// svetlo seryi
0086   //    gStyle->SetFrameFillColor(20);// svetlo kori4nev
0087   //        gStyle->SetFrameFillColor(33);// sine-seryi
0088   //    gStyle->SetFrameFillColor(40);// fiolet-seryi
0089   //    gStyle->SetFrameFillColor(23);// sv.kor
0090 
0091   //---=[ Pad style ]=----------------
0092   gStyle->SetPadTopMargin(TopOffset);
0093   gStyle->SetPadBottomMargin(LeftOffset);
0094   gStyle->SetPadRightMargin(TopOffset);
0095   gStyle->SetPadLeftMargin(LeftOffset);
0096   //---=[ SetCanvasDef           ]=----------------
0097   //======================================================================
0098   //
0099   // Connect the input file and get the 2-d histogram in memory
0100   //======================================================================
0101   //  TBrowser *b = new TBrowser
0102 
0103   //    TFile *hfile1= new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/CMTweb/histos/Global_234556.root", "READ");
0104   //    TFile *hfile1= new TFile("Global_321177_41.root", "READ");
0105   //    TFile *hfile1= new TFile("Global_321177_ls1to600.root", "READ");
0106   //    TFile *hfile1= new TFile("Global_321177_ls1to600.root_no41", "READ");
0107   //    TFile *hfile1= new TFile("Global_325001_ls1to600.root", "READ");
0108   //    TFile *hfile1= new TFile("Global_RBX_325001_40.root", "READ");
0109   //    TFile *hfile1= new TFile("Global_RBX_325001_ls1to600.root", "READ");
0110   ////////////////////////////////////////////////////////////
0111   //    TFile *hfile1= new TFile("Global_321177_41_abortgap.root", "READ");
0112   //    TFile *hfile1= new TFile("Global_321177_ls1to600_abortgap.root", "READ");
0113   //    TFile *hfile1= new TFile("Global_321177_ls1to600_abortgap_no41.root", "READ");
0114   //    TFile *hfile1= new TFile("Global_325001_ls1to600_abortgap.root", "READ");
0115   //    TFile *hfile1= new TFile("Global_321624_1.root", "READ");
0116   //    TFile *hfile1= new TFile("Global_321625.root", "READ");
0117   //    TFile *hfile1= new TFile("Global_321313.root", "READ");
0118   //    TFile *hfile1= new TFile("Global_RBX_325001.root", "READ");
0119 
0120   //    TFile *hfile1= new TFile("Global_RBX_325001test.root", "READ");
0121   //    TFile *hfile1= new TFile("Global_RBX_321177test.root", "READ");
0122   //    TFile *hfile1= new TFile("LED_327785test.root", "READ");
0123 
0124   TFile *hfile1 = new TFile("Global_RBX_325001.root", "READ");
0125   ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
0126 
0127   //    getchar();
0128   //
0129   TPostScript psfile("zadcamplitude.ps", 111);
0130 
0131   //
0132 
0133   TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
0134 
0135   //=============================================================================================== 1
0136   //======================================================================
0137   //======================================================================
0138   //================
0139   //======================================================================
0140   c1->Clear();
0141   c1->Divide(2, 2);
0142 
0143   TH2F *twoddepth1hb1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HB");
0144   TH2F *twoddepth1hb0 = (TH2F *)hfile1->Get("h_mapDepth1_HB");
0145   twoddepth1hb1->Sumw2();
0146   twoddepth1hb0->Sumw2();
0147   //    if(twoddepth1hb0->IsA()->InheritsFrom("TH2F")){
0148   TH2F *Cdepth1hbff = (TH2F *)twoddepth1hb1->Clone("Cdepth1hbff");
0149   Cdepth1hbff->Divide(twoddepth1hb1, twoddepth1hb0, 1, 1, "B");  // average A
0150   Cdepth1hbff->Sumw2();
0151   //    }
0152 
0153   c1->cd(1);
0154   TH2F *twoedepth1hb1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HB");
0155   TH2F *twoedepth1hb0 = (TH2F *)hfile1->Get("h_mapDepth1_HB");
0156   twoedepth1hb1->Sumw2();
0157   twoedepth1hb0->Sumw2();
0158   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0159   TH2F *Cdepth1hbfz225 = (TH2F *)twoedepth1hb1->Clone("Cdepth1hbfz225");
0160   Cdepth1hbfz225->Divide(twoedepth1hb1, twoedepth1hb0, 1, 1, "B");  // just RATE
0161   Cdepth1hbfz225->Sumw2();
0162   //    }
0163   gPad->SetGridy();
0164   gPad->SetGridx();
0165   gPad->SetLogz();
0166   Cdepth1hbfz225->SetMarkerStyle(20);
0167   Cdepth1hbfz225->SetMarkerSize(0.4);
0168   //    Cdepth1hbfz225->GetYaxis()->SetLabelSize(0.04);
0169   Cdepth1hbfz225->GetZaxis()->SetLabelSize(0.08);
0170   Cdepth1hbfz225->SetXTitle("#eta \b");
0171   Cdepth1hbfz225->SetYTitle("#phi \b");
0172   Cdepth1hbfz225->SetZTitle("Rate for ADCAmpl in each event & cell - HB Depth1 \b");
0173   Cdepth1hbfz225->SetMarkerColor(2);
0174   Cdepth1hbfz225->SetLineColor(2);
0175   Cdepth1hbfz225->SetMaximum(1.000);
0176   Cdepth1hbfz225->SetMinimum(0.00001);
0177   //      Cdepth1hbfz225->SetMinimum(0.0001);
0178   Cdepth1hbfz225->Draw("COLZ");
0179 
0180   c1->cd(2);
0181 
0182   //         TH1F *adepth1hb= (TH1F*)hfile1->Get("h_ADCAmpl_HB");
0183   //    TH1F *adepth1hb= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HB");
0184   TH1F *adepth1hb = (TH1F *)hfile1->Get("h_ADCAmplZoom_HB");
0185   gPad->SetLogy();
0186   //      gPad->SetLogx();
0187   adepth1hb->SetMarkerStyle(20);
0188   adepth1hb->SetMarkerSize(0.8);
0189   adepth1hb->GetYaxis()->SetLabelSize(0.04);
0190   adepth1hb->SetXTitle("ADCAmpl in each event & cell HB \b");
0191   adepth1hb->SetMarkerColor(2);
0192   adepth1hb->SetLineColor(2);
0193   adepth1hb->Draw("");
0194 
0195   c1->cd(3);
0196 
0197   ///////////////////////////////////////
0198   TH2F *Diffe_Depth1hb = (TH2F *)Cdepth1hbff->Clone("Diffe_Depth1hb");
0199   for (int i = 1; i <= Cdepth1hbff->GetXaxis()->GetNbins(); i++) {
0200     for (int j = 1; j <= Cdepth1hbff->GetYaxis()->GetNbins(); j++) {
0201       //    TH2F* Cdepth1hbff = (TH2F*)twoddepth1hb1->Clone("Cdepth1hbff");
0202       double ccc1 = Cdepth1hbff->GetBinContent(i, j);
0203       Diffe_Depth1hb->SetBinContent(i, j, 0.);
0204       //    if(ccc1>0.) cout<<" ibin=     "<< i <<" jbin=     "<< j <<" nevents=     "<< ccc1 <<endl;
0205 
0206       //    if(ccc1 <18|| ccc1>30.)  Diffe_Depth1hb->SetBinContent(i,j,ccc1);
0207       //    if(ccc1 <500.)  Diffe_Depth1hb->SetBinContent(i,j,ccc1);
0208       Diffe_Depth1hb->SetBinContent(i, j, ccc1);
0209     }
0210   }
0211   gPad->SetGridy();
0212   gPad->SetGridx();
0213   gPad->SetLogz();
0214   Diffe_Depth1hb->SetMarkerStyle(20);
0215   Diffe_Depth1hb->SetMarkerSize(0.4);
0216   Diffe_Depth1hb->GetZaxis()->SetLabelSize(0.08);
0217   //    Diffe_Depth1hb->SetTitle("any Error, HB Depth1hb \n");
0218   Diffe_Depth1hb->SetXTitle("#eta \b");
0219   Diffe_Depth1hb->SetYTitle("#phi \b");
0220   Diffe_Depth1hb->SetZTitle("<ADCAmpl> - HB Depth1 \b");
0221   Diffe_Depth1hb->SetMarkerColor(2);
0222   Diffe_Depth1hb->SetLineColor(2);
0223   //      Diffe_Depth1hb->SetMaximum(1.000);
0224   //      Diffe_Depth1hb->SetMinimum(0.0000001);
0225   Diffe_Depth1hb->Draw("COLZ");
0226 
0227   c1->cd(4);
0228   //    TH1F* diffADCAmpl_Depth1hb = new TH1F("diffADCAmpl_Depth1hb","", 100, 0.,3000.);
0229   TH1F *diffADCAmpl_Depth1hb = new TH1F("diffADCAmpl_Depth1hb", "", 250, 0., 1000.);
0230   TH2F *Cdepth1hbfw = (TH2F *)Cdepth1hbff->Clone("Cdepth1hbfw");
0231   //    TH2F* Cdepth1hbfw = (TH2F*)twoddepth1hb1->Clone("diffADCAmpl_Depth1hb");
0232   for (int i = 1; i <= Cdepth1hbfw->GetXaxis()->GetNbins(); i++) {
0233     for (int j = 1; j <= Cdepth1hbfw->GetYaxis()->GetNbins(); j++) {
0234       if (Cdepth1hbfw->GetBinContent(i, j) != 0) {
0235         double ccc1 = Cdepth1hbfw->GetBinContent(i, j);
0236         diffADCAmpl_Depth1hb->Fill(ccc1);
0237       }
0238     }
0239   }
0240   gPad->SetLogy();
0241   diffADCAmpl_Depth1hb->SetMarkerStyle(20);
0242   diffADCAmpl_Depth1hb->SetMarkerSize(0.4);
0243   diffADCAmpl_Depth1hb->GetYaxis()->SetLabelSize(0.04);
0244   diffADCAmpl_Depth1hb->SetXTitle("<ADCAmpl> in each cell \b");
0245   diffADCAmpl_Depth1hb->SetMarkerColor(2);
0246   diffADCAmpl_Depth1hb->SetLineColor(2);
0247   diffADCAmpl_Depth1hb->Draw("");
0248 
0249   c1->Update();
0250 
0251   delete twoddepth1hb0;
0252   delete twoddepth1hb1;
0253   delete adepth1hb;
0254   delete Cdepth1hbfz225;
0255   delete Cdepth1hbff;
0256   delete diffADCAmpl_Depth1hb;
0257   delete Diffe_Depth1hb;
0258 
0259   //=============================================================================================== 2
0260   //======================================================================
0261   //======================================================================
0262   //================
0263   //======================================================================
0264   c1->Clear();
0265   c1->Divide(2, 2);
0266 
0267   TH2F *twoddepth2hb1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HB");
0268   TH2F *twoddepth2hb0 = (TH2F *)hfile1->Get("h_mapDepth2_HB");
0269   twoddepth2hb1->Sumw2();
0270   twoddepth2hb0->Sumw2();
0271   //    if(twoddepth2hb0->IsA()->InheritsFrom("TH2F")){
0272   TH2F *Cdepth2hbff = (TH2F *)twoddepth2hb1->Clone("Cdepth2hbff");
0273   Cdepth2hbff->Divide(twoddepth2hb1, twoddepth2hb0, 1, 1, "B");
0274   Cdepth2hbff->Sumw2();
0275   //    }
0276   c1->cd(1);
0277   TH2F *twoedepth2hb1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HB");
0278   TH2F *twoedepth2hb0 = (TH2F *)hfile1->Get("h_mapDepth2_HB");
0279   twoedepth2hb1->Sumw2();
0280   twoedepth2hb0->Sumw2();
0281   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0282   TH2F *Cdepth2hbfz225 = (TH2F *)twoedepth2hb1->Clone("Cdepth2hbfz225");
0283   Cdepth2hbfz225->Divide(twoedepth2hb1, twoedepth2hb0, 1, 1, "B");
0284   Cdepth2hbfz225->Sumw2();
0285   //    }
0286   gPad->SetGridy();
0287   gPad->SetGridx();
0288   gPad->SetLogz();
0289   Cdepth2hbfz225->SetMarkerStyle(20);
0290   Cdepth2hbfz225->SetMarkerSize(0.4);
0291   //    Cdepth2hbfz225->GetYaxis()->SetLabelSize(0.04);
0292   Cdepth2hbfz225->GetZaxis()->SetLabelSize(0.08);
0293   Cdepth2hbfz225->SetXTitle("#eta \b");
0294   Cdepth2hbfz225->SetYTitle("#phi \b");
0295   Cdepth2hbfz225->SetZTitle("Rate for ADCAmpl in each event & cell - HB Depth2 \b");
0296   Cdepth2hbfz225->SetMarkerColor(2);
0297   Cdepth2hbfz225->SetLineColor(2);
0298   Cdepth2hbfz225->SetMaximum(1.000);
0299   Cdepth2hbfz225->SetMinimum(0.00001);
0300   //      Cdepth2hbfz225->SetMinimum(0.0001);
0301   Cdepth2hbfz225->Draw("COLZ");
0302 
0303   c1->cd(2);
0304 
0305   //         TH1F *adepth2hb= (TH1F*)hfile1->Get("h_ADCAmpl_HB");
0306   TH1F *adepth2hb = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HB");
0307   //      TH1F *adepth2hb= (TH1F*)hfile1->Get("h_ADCAmplZoom_HB");
0308   gPad->SetLogy();
0309   //      gPad->SetLogx();
0310   adepth2hb->SetMarkerStyle(20);
0311   adepth2hb->SetMarkerSize(0.8);
0312   adepth2hb->GetYaxis()->SetLabelSize(0.04);
0313   adepth2hb->SetXTitle("ADCAmpl in each event & cell HB \b");
0314   adepth2hb->SetMarkerColor(2);
0315   adepth2hb->SetLineColor(2);
0316   adepth2hb->Draw("");
0317 
0318   c1->cd(3);
0319 
0320   ///////////////////////////////////////
0321   TH2F *Diffe_Depth2hb = (TH2F *)Cdepth2hbff->Clone("Diffe_Depth2hb");
0322   for (int i = 1; i <= Cdepth2hbff->GetXaxis()->GetNbins(); i++) {
0323     for (int j = 1; j <= Cdepth2hbff->GetYaxis()->GetNbins(); j++) {
0324       //    TH2F* Cdepth2hbff = (TH2F*)twoddepth2hb1->Clone("Cdepth2hbff");
0325       double ccc1 = Cdepth2hbff->GetBinContent(i, j);
0326       Diffe_Depth2hb->SetBinContent(i, j, 0.);
0327       if (ccc1 < 500.)
0328         Diffe_Depth2hb->SetBinContent(i, j, ccc1);
0329     }
0330   }
0331   gPad->SetGridy();
0332   gPad->SetGridx();
0333   gPad->SetLogz();
0334   Diffe_Depth2hb->SetMarkerStyle(20);
0335   Diffe_Depth2hb->SetMarkerSize(0.4);
0336   Diffe_Depth2hb->GetZaxis()->SetLabelSize(0.08);
0337   //    Diffe_Depth2hb->SetTitle("any Error, HB Depth2hb \n");
0338   Diffe_Depth2hb->SetXTitle("#eta \b");
0339   Diffe_Depth2hb->SetYTitle("#phi \b");
0340   Diffe_Depth2hb->SetZTitle("<ADCAmpl> smalle 508 - HB Depth2hb \b");
0341   Diffe_Depth2hb->SetMarkerColor(2);
0342   Diffe_Depth2hb->SetLineColor(2);
0343   //      Diffe_Depth2hb->SetMaximum(1.000);
0344   //      Diffe_Depth2hb->SetMinimum(0.0000001);
0345   Diffe_Depth2hb->Draw("COLZ");
0346 
0347   c1->cd(4);
0348   //    TH1F* diffADCAmpl_Depth2hb = new TH1F("diffADCAmpl_Depth2hb","", 100, 0.,3000.);
0349   TH1F *diffADCAmpl_Depth2hb = new TH1F("diffADCAmpl_Depth2hb", "", 250, 0., 1000.);
0350   TH2F *Cdepth2hbfw = (TH2F *)Cdepth2hbff->Clone("diffADCAmpl_Depth2hb");
0351   //      TH2F* Cdepth2hbfw = (TH2F*)twoddepth2hb1->Clone("diffADCAmpl_Depth2hb");
0352   for (int i = 1; i <= Cdepth2hbfw->GetXaxis()->GetNbins(); i++) {
0353     for (int j = 1; j <= Cdepth2hbfw->GetYaxis()->GetNbins(); j++) {
0354       if (Cdepth2hbfw->GetBinContent(i, j) != 0) {
0355         double ccc1 = Cdepth2hbfw->GetBinContent(i, j);
0356         diffADCAmpl_Depth2hb->Fill(ccc1);
0357       }
0358     }
0359   }
0360   gPad->SetLogy();
0361   diffADCAmpl_Depth2hb->SetMarkerStyle(20);
0362   diffADCAmpl_Depth2hb->SetMarkerSize(0.4);
0363   diffADCAmpl_Depth2hb->GetYaxis()->SetLabelSize(0.04);
0364   diffADCAmpl_Depth2hb->SetXTitle("<ADCAmpl> in each cell \b");
0365   diffADCAmpl_Depth2hb->SetMarkerColor(2);
0366   diffADCAmpl_Depth2hb->SetLineColor(2);
0367   diffADCAmpl_Depth2hb->Draw("");
0368 
0369   c1->Update();
0370 
0371   delete twoddepth2hb0;
0372   delete twoddepth2hb1;
0373   delete adepth2hb;
0374   delete Cdepth2hbfz225;
0375   delete Cdepth2hbff;
0376   delete diffADCAmpl_Depth2hb;
0377   delete Diffe_Depth2hb;
0378 
0379   //=============================================================================================== 3
0380   //======================================================================
0381   //======================================================================
0382   //================
0383   //======================================================================
0384   c1->Clear();
0385   c1->Divide(2, 2);
0386 
0387   TH2F *twoddepth1he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HE");
0388   TH2F *twoddepth1he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
0389   twoddepth1he1->Sumw2();
0390   twoddepth1he0->Sumw2();
0391   //    if(twoddepth1he0->IsA()->InheritsFrom("TH2F")){
0392   TH2F *Cdepth1heff = (TH2F *)twoddepth1he1->Clone("Cdepth1heff");
0393   Cdepth1heff->Divide(twoddepth1he1, twoddepth1he0, 1, 1, "B");
0394   Cdepth1heff->Sumw2();
0395   //    }
0396   c1->cd(1);
0397   TH2F *twoedepth1he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HE");
0398   TH2F *twoedepth1he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
0399   twoedepth1he1->Sumw2();
0400   twoedepth1he0->Sumw2();
0401   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0402   TH2F *Cdepth1hefz225 = (TH2F *)twoedepth1he1->Clone("Cdepth1hefz225");
0403   Cdepth1hefz225->Divide(twoedepth1he1, twoedepth1he0, 1, 1, "B");
0404   Cdepth1hefz225->Sumw2();
0405   //    }
0406   gPad->SetGridy();
0407   gPad->SetGridx();
0408   gPad->SetLogz();
0409   Cdepth1hefz225->SetMarkerStyle(20);
0410   Cdepth1hefz225->SetMarkerSize(0.4);
0411   //    Cdepth1hefz225->GetYaxis()->SetLabelSize(0.04);
0412   Cdepth1hefz225->GetZaxis()->SetLabelSize(0.08);
0413   Cdepth1hefz225->SetXTitle("#eta \b");
0414   Cdepth1hefz225->SetYTitle("#phi \b");
0415   Cdepth1hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth1 \b");
0416   Cdepth1hefz225->SetMarkerColor(2);
0417   Cdepth1hefz225->SetLineColor(2);
0418   Cdepth1hefz225->SetMaximum(1.000);
0419   Cdepth1hefz225->SetMinimum(0.00001);
0420   //      Cdepth1hefz225->SetMinimum(0.0001);
0421   Cdepth1hefz225->Draw("COLZ");
0422 
0423   c1->cd(2);
0424 
0425   //         TH1F *adepth1he= (TH1F*)hfile1->Get("h_ADCAmpl_HE");
0426   TH1F *adepth1he = (TH1F *)hfile1->Get("h_ADCAmpl345_HE");
0427   //      TH1F *adepth1he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
0428   gPad->SetLogy();
0429   //      gPad->SetLogx();
0430   adepth1he->SetMarkerStyle(20);
0431   adepth1he->SetMarkerSize(0.8);
0432   adepth1he->GetYaxis()->SetLabelSize(0.04);
0433   adepth1he->SetXTitle("ADCAmpl in each event & cell HE \b");
0434   adepth1he->SetMarkerColor(2);
0435   adepth1he->SetLineColor(2);
0436   adepth1he->Draw("");
0437 
0438   c1->cd(3);
0439 
0440   ///////////////////////////////////////
0441   TH2F *Diffe_Depth1he = (TH2F *)Cdepth1heff->Clone("Diffe_Depth1he");
0442   for (int i = 1; i <= Cdepth1heff->GetXaxis()->GetNbins(); i++) {
0443     for (int j = 1; j <= Cdepth1heff->GetYaxis()->GetNbins(); j++) {
0444       //    TH2F* Cdepth1heff = (TH2F*)twoddepth1he1->Clone("Cdepth1heff");
0445       double ccc1 = Cdepth1heff->GetBinContent(i, j);
0446       //    if(ccc1>0.) cout<<"HE ibin=     "<< i <<" jbin=     "<< j <<" nevents=     "<< ccc1 <<endl;
0447       Diffe_Depth1he->SetBinContent(i, j, 0.);
0448       if (ccc1 < 20000000.)
0449         Diffe_Depth1he->SetBinContent(i, j, ccc1);
0450     }
0451   }
0452   gPad->SetGridy();
0453   gPad->SetGridx();
0454   gPad->SetLogz();
0455   Diffe_Depth1he->SetMarkerStyle(20);
0456   Diffe_Depth1he->SetMarkerSize(0.4);
0457   Diffe_Depth1he->GetZaxis()->SetLabelSize(0.08);
0458   //    Diffe_Depth1he->SetTitle("any Error, HE Depth1he \n");
0459   Diffe_Depth1he->SetXTitle("#eta \b");
0460   Diffe_Depth1he->SetYTitle("#phi \b");
0461   Diffe_Depth1he->SetZTitle("<ADCAmpl> smalle 20000008 - HE Depth1he \b");
0462   Diffe_Depth1he->SetMarkerColor(2);
0463   Diffe_Depth1he->SetLineColor(2);
0464   //      Diffe_Depth1he->SetMaximum(1.000);
0465   //      Diffe_Depth1he->SetMinimum(0.0000001);
0466   Diffe_Depth1he->Draw("COLZ");
0467 
0468   c1->cd(4);
0469   //    TH1F* diffADCAmpl_Depth1he = new TH1F("diffADCAmpl_Depth1he","", 100, 0.,3000.);
0470   TH1F *diffADCAmpl_Depth1he = new TH1F("diffADCAmpl_Depth1he", "", 100, 0., 5000.);
0471   TH2F *Cdepth1hefw = (TH2F *)Cdepth1heff->Clone("diffADCAmpl_Depth1he");
0472   //      TH2F* Cdepth1hefw = (TH2F*)twoddepth1he1->Clone("diffADCAmpl_Depth1he");
0473   for (int i = 1; i <= Cdepth1hefw->GetXaxis()->GetNbins(); i++) {
0474     for (int j = 1; j <= Cdepth1hefw->GetYaxis()->GetNbins(); j++) {
0475       if (Cdepth1hefw->GetBinContent(i, j) != 0) {
0476         double ccc1 = Cdepth1hefw->GetBinContent(i, j);
0477         diffADCAmpl_Depth1he->Fill(ccc1);
0478       }
0479     }
0480   }
0481   gPad->SetLogy();
0482   diffADCAmpl_Depth1he->SetMarkerStyle(20);
0483   diffADCAmpl_Depth1he->SetMarkerSize(0.4);
0484   diffADCAmpl_Depth1he->GetYaxis()->SetLabelSize(0.04);
0485   diffADCAmpl_Depth1he->SetXTitle("<ADCAmpl> in each cell \b");
0486   diffADCAmpl_Depth1he->SetMarkerColor(2);
0487   diffADCAmpl_Depth1he->SetLineColor(2);
0488   diffADCAmpl_Depth1he->Draw("");
0489 
0490   c1->Update();
0491 
0492   delete twoddepth1he0;
0493   delete twoddepth1he1;
0494   delete adepth1he;
0495   delete Cdepth1hefz225;
0496   delete Cdepth1heff;
0497   delete diffADCAmpl_Depth1he;
0498   delete Diffe_Depth1he;
0499 
0500   //=============================================================================================== 4
0501   //======================================================================
0502   //======================================================================
0503   //================
0504   //======================================================================
0505   c1->Clear();
0506   c1->Divide(2, 2);
0507 
0508   TH2F *twoddepth2he1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HE");
0509   TH2F *twoddepth2he0 = (TH2F *)hfile1->Get("h_mapDepth2_HE");
0510   twoddepth2he1->Sumw2();
0511   twoddepth2he0->Sumw2();
0512   //    if(twoddepth2he0->IsA()->InheritsFrom("TH2F")){
0513   TH2F *Cdepth2heff = (TH2F *)twoddepth2he1->Clone("Cdepth2heff");
0514   Cdepth2heff->Divide(twoddepth2he1, twoddepth2he0, 1, 1, "B");
0515   Cdepth2heff->Sumw2();
0516   //    }
0517   c1->cd(1);
0518   TH2F *twoedepth2he1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HE");
0519   TH2F *twoedepth2he0 = (TH2F *)hfile1->Get("h_mapDepth2_HE");
0520   twoedepth2he1->Sumw2();
0521   twoedepth2he0->Sumw2();
0522   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0523   TH2F *Cdepth2hefz225 = (TH2F *)twoedepth2he1->Clone("Cdepth2hefz225");
0524   Cdepth2hefz225->Divide(twoedepth2he1, twoedepth2he0, 1, 1, "B");
0525   Cdepth2hefz225->Sumw2();
0526   //    }
0527   gPad->SetGridy();
0528   gPad->SetGridx();
0529   gPad->SetLogz();
0530   Cdepth2hefz225->SetMarkerStyle(20);
0531   Cdepth2hefz225->SetMarkerSize(0.4);
0532   //    Cdepth2hefz225->GetYaxis()->SetLabelSize(0.04);
0533   Cdepth2hefz225->GetZaxis()->SetLabelSize(0.08);
0534   Cdepth2hefz225->SetXTitle("#eta \b");
0535   Cdepth2hefz225->SetYTitle("#phi \b");
0536   Cdepth2hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth2 \b");
0537   Cdepth2hefz225->SetMarkerColor(2);
0538   Cdepth2hefz225->SetLineColor(2);
0539   Cdepth2hefz225->SetMaximum(1.000);
0540   Cdepth2hefz225->SetMinimum(0.00001);
0541   //      Cdepth2hefz225->SetMinimum(0.0001);
0542   Cdepth2hefz225->Draw("COLZ");
0543 
0544   c1->cd(2);
0545 
0546   TH1F *adepth2he = (TH1F *)hfile1->Get("h_ADCAmpl_HE");
0547   //     TH1F *adepth2he= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HE");
0548   //      TH1F *adepth2he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
0549   gPad->SetLogy();
0550   //      gPad->SetLogx();
0551   adepth2he->SetMarkerStyle(20);
0552   adepth2he->SetMarkerSize(0.8);
0553   adepth2he->GetYaxis()->SetLabelSize(0.04);
0554   adepth2he->SetXTitle("ADCAmpl in each event & cell HE \b");
0555   adepth2he->SetMarkerColor(2);
0556   adepth2he->SetLineColor(2);
0557   adepth2he->Draw("");
0558 
0559   c1->cd(3);
0560 
0561   ///////////////////////////////////////
0562   TH2F *Diffe_Depth2he = (TH2F *)Cdepth2heff->Clone("Diffe_Depth2he");
0563   for (int i = 1; i <= Cdepth2heff->GetXaxis()->GetNbins(); i++) {
0564     for (int j = 1; j <= Cdepth2heff->GetYaxis()->GetNbins(); j++) {
0565       //    TH2F* Cdepth2heff = (TH2F*)twoddepth2he1->Clone("Cdepth2heff");
0566       double ccc1 = Cdepth2heff->GetBinContent(i, j);
0567       Diffe_Depth2he->SetBinContent(i, j, 0.);
0568       if (ccc1 < 20000000.)
0569         Diffe_Depth2he->SetBinContent(i, j, ccc1);
0570     }
0571   }
0572   gPad->SetGridy();
0573   gPad->SetGridx();
0574   gPad->SetLogz();
0575   Diffe_Depth2he->SetMarkerStyle(20);
0576   Diffe_Depth2he->SetMarkerSize(0.4);
0577   Diffe_Depth2he->GetZaxis()->SetLabelSize(0.08);
0578   //    Diffe_Depth2he->SetTitle("any Error, HE Depth2he \n");
0579   Diffe_Depth2he->SetXTitle("#eta \b");
0580   Diffe_Depth2he->SetYTitle("#phi \b");
0581   Diffe_Depth2he->SetZTitle("<ADCAmpl> smalle 20000000 - HE Depth2he \b");
0582   Diffe_Depth2he->SetMarkerColor(2);
0583   Diffe_Depth2he->SetLineColor(2);
0584   //      Diffe_Depth2he->SetMaximum(1.000);
0585   //      Diffe_Depth2he->SetMinimum(0.0000001);
0586   Diffe_Depth2he->Draw("COLZ");
0587 
0588   c1->cd(4);
0589   //    TH1F* diffADCAmpl_Depth2he = new TH1F("diffADCAmpl_Depth2he","", 100, 0.,3000.);
0590   TH1F *diffADCAmpl_Depth2he = new TH1F("diffADCAmpl_Depth2he", "", 100, 0., 5000.);
0591   TH2F *Cdepth2hefw = (TH2F *)Cdepth2heff->Clone("diffADCAmpl_Depth2he");
0592   //      TH2F* Cdepth2hefw = (TH2F*)twoddepth2he1->Clone("diffADCAmpl_Depth2he");
0593   for (int i = 1; i <= Cdepth2hefw->GetXaxis()->GetNbins(); i++) {
0594     for (int j = 1; j <= Cdepth2hefw->GetYaxis()->GetNbins(); j++) {
0595       if (Cdepth2hefw->GetBinContent(i, j) != 0) {
0596         double ccc1 = Cdepth2hefw->GetBinContent(i, j);
0597         diffADCAmpl_Depth2he->Fill(ccc1);
0598       }
0599     }
0600   }
0601   gPad->SetLogy();
0602   diffADCAmpl_Depth2he->SetMarkerStyle(20);
0603   diffADCAmpl_Depth2he->SetMarkerSize(0.4);
0604   diffADCAmpl_Depth2he->GetYaxis()->SetLabelSize(0.04);
0605   diffADCAmpl_Depth2he->SetXTitle("<ADCAmpl> in each cell \b");
0606   diffADCAmpl_Depth2he->SetMarkerColor(2);
0607   diffADCAmpl_Depth2he->SetLineColor(2);
0608   diffADCAmpl_Depth2he->Draw("");
0609 
0610   c1->Update();
0611 
0612   delete twoddepth2he0;
0613   delete twoddepth2he1;
0614   delete adepth2he;
0615   delete Cdepth2hefz225;
0616   delete Cdepth2heff;
0617   delete diffADCAmpl_Depth2he;
0618   delete Diffe_Depth2he;
0619 
0620   //=============================================================================================== 5
0621   //======================================================================
0622   //======================================================================
0623   //================
0624   //======================================================================
0625   c1->Clear();
0626   c1->Divide(2, 2);
0627 
0628   TH2F *twoddepth3he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HE");
0629   TH2F *twoddepth3he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
0630   twoddepth3he1->Sumw2();
0631   twoddepth3he0->Sumw2();
0632   //    if(twoddepth3he0->IsA()->InheritsFrom("TH2F")){
0633   TH2F *Cdepth3heff = (TH2F *)twoddepth3he1->Clone("Cdepth3heff");
0634   Cdepth3heff->Divide(twoddepth3he1, twoddepth3he0, 1, 1, "B");
0635   Cdepth3heff->Sumw2();
0636   //    }
0637   c1->cd(1);
0638   TH2F *twoedepth3he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HE");
0639   TH2F *twoedepth3he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
0640   twoedepth3he1->Sumw2();
0641   twoedepth3he0->Sumw2();
0642   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0643   TH2F *Cdepth3hefz225 = (TH2F *)twoedepth3he1->Clone("Cdepth3hefz225");
0644   Cdepth3hefz225->Divide(twoedepth3he1, twoedepth3he0, 1, 1, "B");
0645   Cdepth3hefz225->Sumw2();
0646   //    }
0647   gPad->SetGridy();
0648   gPad->SetGridx();
0649   gPad->SetLogz();
0650   Cdepth3hefz225->SetMarkerStyle(20);
0651   Cdepth3hefz225->SetMarkerSize(0.4);
0652   //    Cdepth3hefz225->GetYaxis()->SetLabelSize(0.04);
0653   Cdepth3hefz225->GetZaxis()->SetLabelSize(0.08);
0654   Cdepth3hefz225->SetXTitle("#eta \b");
0655   Cdepth3hefz225->SetYTitle("#phi \b");
0656   Cdepth3hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth1 \b");
0657   Cdepth3hefz225->SetMarkerColor(2);
0658   Cdepth3hefz225->SetLineColor(2);
0659   Cdepth3hefz225->SetMaximum(1.000);
0660   Cdepth3hefz225->SetMinimum(0.00001);
0661   //      Cdepth3hefz225->SetMinimum(0.0001);
0662   Cdepth3hefz225->Draw("COLZ");
0663 
0664   c1->cd(2);
0665 
0666   //         TH1F *adepth3he= (TH1F*)hfile1->Get("h_ADCAmpl_HE");
0667   TH1F *adepth3he = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HE");
0668   //      TH1F *adepth3he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
0669   gPad->SetLogy();
0670   //      gPad->SetLogx();
0671   adepth3he->SetMarkerStyle(20);
0672   adepth3he->SetMarkerSize(0.8);
0673   adepth3he->GetYaxis()->SetLabelSize(0.04);
0674   adepth3he->SetXTitle("ADCAmpl in each event & cell HE \b");
0675   adepth3he->SetMarkerColor(2);
0676   adepth3he->SetLineColor(2);
0677   adepth3he->Draw("");
0678 
0679   c1->cd(3);
0680 
0681   ///////////////////////////////////////
0682   TH2F *Diffe_Depth3he = (TH2F *)Cdepth3heff->Clone("Diffe_Depth3he");
0683   for (int i = 1; i <= Cdepth3heff->GetXaxis()->GetNbins(); i++) {
0684     for (int j = 1; j <= Cdepth3heff->GetYaxis()->GetNbins(); j++) {
0685       //    TH2F* Cdepth3heff = (TH2F*)twoddepth3he1->Clone("Cdepth3heff");
0686       double ccc1 = Cdepth3heff->GetBinContent(i, j);
0687       Diffe_Depth3he->SetBinContent(i, j, 0.);
0688       if (ccc1 < 20000000.)
0689         Diffe_Depth3he->SetBinContent(i, j, ccc1);
0690     }
0691   }
0692   gPad->SetGridy();
0693   gPad->SetGridx();
0694   gPad->SetLogz();
0695   Diffe_Depth3he->SetMarkerStyle(20);
0696   Diffe_Depth3he->SetMarkerSize(0.4);
0697   Diffe_Depth3he->GetZaxis()->SetLabelSize(0.08);
0698   //    Diffe_Depth3he->SetTitle("any Error, HE Depth3he \n");
0699   Diffe_Depth3he->SetXTitle("#eta \b");
0700   Diffe_Depth3he->SetYTitle("#phi \b");
0701   Diffe_Depth3he->SetZTitle("<ADCAmpl>smalle 20000000 - HE Depth3he \b");
0702   Diffe_Depth3he->SetMarkerColor(2);
0703   Diffe_Depth3he->SetLineColor(2);
0704   //      Diffe_Depth3he->SetMaximum(1.000);
0705   //      Diffe_Depth3he->SetMinimum(0.0000001);
0706   Diffe_Depth3he->Draw("COLZ");
0707 
0708   c1->cd(4);
0709   //    TH1F* diffADCAmpl_Depth3he = new TH1F("diffADCAmpl_Depth3he","", 100, 0.,3000.);
0710   TH1F *diffADCAmpl_Depth3he = new TH1F("diffADCAmpl_Depth3he", "", 100, 0., 5000.);
0711   TH2F *Cdepth3hefw = (TH2F *)Cdepth3heff->Clone("diffADCAmpl_Depth3he");
0712   //      TH2F* Cdepth3hefw = (TH2F*)twoddepth3he1->Clone("diffADCAmpl_Depth3he");
0713   for (int i = 1; i <= Cdepth3hefw->GetXaxis()->GetNbins(); i++) {
0714     for (int j = 1; j <= Cdepth3hefw->GetYaxis()->GetNbins(); j++) {
0715       if (Cdepth3hefw->GetBinContent(i, j) != 0) {
0716         double ccc1 = Cdepth3hefw->GetBinContent(i, j);
0717         diffADCAmpl_Depth3he->Fill(ccc1);
0718       }
0719     }
0720   }
0721   gPad->SetLogy();
0722   diffADCAmpl_Depth3he->SetMarkerStyle(20);
0723   diffADCAmpl_Depth3he->SetMarkerSize(0.4);
0724   diffADCAmpl_Depth3he->GetYaxis()->SetLabelSize(0.04);
0725   diffADCAmpl_Depth3he->SetXTitle("<ADCAmpl> in each cell \b");
0726   diffADCAmpl_Depth3he->SetMarkerColor(2);
0727   diffADCAmpl_Depth3he->SetLineColor(2);
0728   diffADCAmpl_Depth3he->Draw("");
0729 
0730   c1->Update();
0731 
0732   delete twoddepth3he0;
0733   delete twoddepth3he1;
0734   delete adepth3he;
0735   delete Cdepth3hefz225;
0736   delete Cdepth3heff;
0737   delete diffADCAmpl_Depth3he;
0738   delete Diffe_Depth3he;
0739 
0740   //=============================================================================================== 6
0741   //======================================================================
0742   //======================================================================
0743   //================
0744   //======================================================================
0745   c1->Clear();
0746   c1->Divide(2, 2);
0747 
0748   TH2F *twoddepth1hf1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HF");
0749   TH2F *twoddepth1hf0 = (TH2F *)hfile1->Get("h_mapDepth1_HF");
0750   twoddepth1hf1->Sumw2();
0751   twoddepth1hf0->Sumw2();
0752   //    if(twoddepth1hf0->IsA()->InheritsFrom("TH2F")){
0753   TH2F *Cdepth1hfff = (TH2F *)twoddepth1hf1->Clone("Cdepth1hfff");
0754   Cdepth1hfff->Divide(twoddepth1hf1, twoddepth1hf0, 1, 1, "B");
0755   Cdepth1hfff->Sumw2();
0756   //    }
0757   c1->cd(1);
0758   TH2F *twoedepth1hf1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HF");
0759   TH2F *twoedepth1hf0 = (TH2F *)hfile1->Get("h_mapDepth1_HF");
0760   twoedepth1hf1->Sumw2();
0761   twoedepth1hf0->Sumw2();
0762   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0763   TH2F *Cdepth1hffz225 = (TH2F *)twoedepth1hf1->Clone("Cdepth1hffz225");
0764   Cdepth1hffz225->Divide(twoedepth1hf1, twoedepth1hf0, 1, 1, "B");
0765   Cdepth1hffz225->Sumw2();
0766   //    }
0767   gPad->SetGridy();
0768   gPad->SetGridx();
0769   gPad->SetLogz();
0770   Cdepth1hffz225->SetMarkerStyle(20);
0771   Cdepth1hffz225->SetMarkerSize(0.4);
0772   //    Cdepth1hffz225->GetYaxis()->SetLabelSize(0.04);
0773   Cdepth1hffz225->GetZaxis()->SetLabelSize(0.08);
0774   Cdepth1hffz225->SetXTitle("#eta \b");
0775   Cdepth1hffz225->SetYTitle("#phi \b");
0776   Cdepth1hffz225->SetZTitle("Rate for ADCAmpl in each event & cell - HF Depth1 \b");
0777   Cdepth1hffz225->SetMarkerColor(2);
0778   Cdepth1hffz225->SetLineColor(2);
0779   Cdepth1hffz225->SetMaximum(1.000);
0780   Cdepth1hffz225->SetMinimum(0.00001);
0781   //      Cdepth1hffz225->SetMinimum(0.0001);
0782   Cdepth1hffz225->Draw("COLZ");
0783 
0784   c1->cd(2);
0785 
0786   TH1F *adepth1hf = (TH1F *)hfile1->Get("h_ADCAmpl_HF");
0787   //      TH1F *adepth1hf= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HF");
0788   //      TH1F *adepth1hf= (TH1F*)hfile1->Get("h_ADCAmplZoom_HF");
0789   gPad->SetLogy();
0790   //      gPad->SetLogx();
0791   adepth1hf->SetMarkerStyle(20);
0792   adepth1hf->SetMarkerSize(0.8);
0793   adepth1hf->GetYaxis()->SetLabelSize(0.04);
0794   adepth1hf->SetXTitle("ADCAmpl in each event & cell HF \b");
0795   adepth1hf->SetMarkerColor(2);
0796   adepth1hf->SetLineColor(2);
0797   adepth1hf->Draw("");
0798 
0799   c1->cd(3);
0800 
0801   ///////////////////////////////////////
0802   TH2F *Diffe_Depth1hf = (TH2F *)Cdepth1hfff->Clone("Diffe_Depth1hf");
0803   for (int i = 1; i <= Cdepth1hfff->GetXaxis()->GetNbins(); i++) {
0804     for (int j = 1; j <= Cdepth1hfff->GetYaxis()->GetNbins(); j++) {
0805       //    TH2F* Cdepth1hfff = (TH2F*)twoddepth1hf1->Clone("Cdepth1hfff");
0806       double ccc1 = Cdepth1hfff->GetBinContent(i, j);
0807       Diffe_Depth1hf->SetBinContent(i, j, 0.);
0808       if (ccc1 < 20000000.)
0809         Diffe_Depth1hf->SetBinContent(i, j, ccc1);
0810     }
0811   }
0812   gPad->SetGridy();
0813   gPad->SetGridx();
0814   gPad->SetLogz();
0815   Diffe_Depth1hf->SetMarkerStyle(20);
0816   Diffe_Depth1hf->SetMarkerSize(0.4);
0817   Diffe_Depth1hf->GetZaxis()->SetLabelSize(0.08);
0818   //    Diffe_Depth1hf->SetTitle("any Error, HF Depth1hf \n");
0819   Diffe_Depth1hf->SetXTitle("#eta \b");
0820   Diffe_Depth1hf->SetYTitle("#phi \b");
0821   Diffe_Depth1hf->SetZTitle("<ADCAmpl> smalle 20000000 - HF Depth1hf \b");
0822   Diffe_Depth1hf->SetMarkerColor(2);
0823   Diffe_Depth1hf->SetLineColor(2);
0824   //      Diffe_Depth1hf->SetMaximum(1.000);
0825   //      Diffe_Depth1hf->SetMinimum(0.0000001);
0826   Diffe_Depth1hf->Draw("COLZ");
0827 
0828   c1->cd(4);
0829   TH1F *diffADCAmpl_Depth1hf = new TH1F("diffADCAmpl_Depth1hf", "", 100, 0., 300.);
0830   //      TH1F* diffADCAmpl_Depth1hf = new TH1F("diffADCAmpl_Depth1hf","", 1000, 0.,1000.);
0831   TH2F *Cdepth1hffw = (TH2F *)Cdepth1hfff->Clone("diffADCAmpl_Depth1hf");
0832   //      TH2F* Cdepth1hffw = (TH2F*)twoddepth1hf1->Clone("diffADCAmpl_Depth1hf");
0833   for (int i = 1; i <= Cdepth1hffw->GetXaxis()->GetNbins(); i++) {
0834     for (int j = 1; j <= Cdepth1hffw->GetYaxis()->GetNbins(); j++) {
0835       if (Cdepth1hffw->GetBinContent(i, j) != 0) {
0836         double ccc1 = Cdepth1hffw->GetBinContent(i, j);
0837         diffADCAmpl_Depth1hf->Fill(ccc1);
0838       }
0839     }
0840   }
0841   gPad->SetLogy();
0842   diffADCAmpl_Depth1hf->SetMarkerStyle(20);
0843   diffADCAmpl_Depth1hf->SetMarkerSize(0.4);
0844   diffADCAmpl_Depth1hf->GetYaxis()->SetLabelSize(0.04);
0845   diffADCAmpl_Depth1hf->SetXTitle("<ADCAmpl> in each cell \b");
0846   diffADCAmpl_Depth1hf->SetMarkerColor(2);
0847   diffADCAmpl_Depth1hf->SetLineColor(2);
0848   diffADCAmpl_Depth1hf->Draw("");
0849 
0850   c1->Update();
0851 
0852   delete twoddepth1hf0;
0853   delete twoddepth1hf1;
0854   delete adepth1hf;
0855   delete Cdepth1hffz225;
0856   delete Cdepth1hfff;
0857   delete diffADCAmpl_Depth1hf;
0858   delete Diffe_Depth1hf;
0859 
0860   //=============================================================================================== 7
0861   //======================================================================
0862   //======================================================================
0863   //================
0864   //======================================================================
0865   c1->Clear();
0866   c1->Divide(2, 2);
0867 
0868   TH2F *twoddepth2hf1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HF");
0869   TH2F *twoddepth2hf0 = (TH2F *)hfile1->Get("h_mapDepth2_HF");
0870   twoddepth2hf1->Sumw2();
0871   twoddepth2hf0->Sumw2();
0872   //    if(twoddepth2hf0->IsA()->InheritsFrom("TH2F")){
0873   TH2F *Cdepth2hfff = (TH2F *)twoddepth2hf1->Clone("Cdepth2hfff");
0874   Cdepth2hfff->Divide(twoddepth2hf1, twoddepth2hf0, 1, 1, "B");
0875   Cdepth2hfff->Sumw2();
0876   //    }
0877   c1->cd(1);
0878   TH2F *twoedepth2hf1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HF");
0879   TH2F *twoedepth2hf0 = (TH2F *)hfile1->Get("h_mapDepth2_HF");
0880   twoedepth2hf1->Sumw2();
0881   twoedepth2hf0->Sumw2();
0882   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
0883   TH2F *Cdepth2hffz225 = (TH2F *)twoedepth2hf1->Clone("Cdepth2hffz225");
0884   Cdepth2hffz225->Divide(twoedepth2hf1, twoedepth2hf0, 1, 1, "B");
0885   Cdepth2hffz225->Sumw2();
0886   //    }
0887   gPad->SetGridy();
0888   gPad->SetGridx();
0889   gPad->SetLogz();
0890   Cdepth2hffz225->SetMarkerStyle(20);
0891   Cdepth2hffz225->SetMarkerSize(0.4);
0892   //    Cdepth2hffz225->GetYaxis()->SetLabelSize(0.04);
0893   Cdepth2hffz225->GetZaxis()->SetLabelSize(0.08);
0894   Cdepth2hffz225->SetXTitle("#eta \b");
0895   Cdepth2hffz225->SetYTitle("#phi \b");
0896   Cdepth2hffz225->SetZTitle("Rate for ADCAmpl in each event & cell - HF Depth2 \b");
0897   Cdepth2hffz225->SetMarkerColor(2);
0898   Cdepth2hffz225->SetLineColor(2);
0899   Cdepth2hffz225->SetMaximum(1.000);
0900   Cdepth2hffz225->SetMinimum(0.00001);
0901   //      Cdepth2hffz225->SetMinimum(0.0001);
0902   Cdepth2hffz225->Draw("COLZ");
0903 
0904   c1->cd(2);
0905 
0906   //         TH1F *adepth2hf= (TH1F*)hfile1->Get("h_ADCAmpl_HF");
0907   TH1F *adepth2hf = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HF");
0908   //      TH1F *adepth2hf= (TH1F*)hfile1->Get("h_ADCAmplZoom_HF");
0909   gPad->SetLogy();
0910   //      gPad->SetLogx();
0911   adepth2hf->SetMarkerStyle(20);
0912   adepth2hf->SetMarkerSize(0.8);
0913   adepth2hf->GetYaxis()->SetLabelSize(0.04);
0914   adepth2hf->SetXTitle("ADCAmpl in each event & cell HF \b");
0915   adepth2hf->SetMarkerColor(2);
0916   adepth2hf->SetLineColor(2);
0917   adepth2hf->Draw("");
0918 
0919   c1->cd(3);
0920 
0921   ///////////////////////////////////////
0922   TH2F *Diffe_Depth2hf = (TH2F *)Cdepth2hfff->Clone("Diffe_Depth2hf");
0923   for (int i = 1; i <= Cdepth2hfff->GetXaxis()->GetNbins(); i++) {
0924     for (int j = 1; j <= Cdepth2hfff->GetYaxis()->GetNbins(); j++) {
0925       //    TH2F* Cdepth2hfff = (TH2F*)twoddepth2hf1->Clone("Cdepth2hfff");
0926       double ccc1 = Cdepth2hfff->GetBinContent(i, j);
0927       Diffe_Depth2hf->SetBinContent(i, j, 0.);
0928       if (ccc1 < 20000000.)
0929         Diffe_Depth2hf->SetBinContent(i, j, ccc1);
0930     }
0931   }
0932   gPad->SetGridy();
0933   gPad->SetGridx();
0934   gPad->SetLogz();
0935   Diffe_Depth2hf->SetMarkerStyle(20);
0936   Diffe_Depth2hf->SetMarkerSize(0.4);
0937   Diffe_Depth2hf->GetZaxis()->SetLabelSize(0.08);
0938   //    Diffe_Depth2hf->SetTitle("any Error, HF Depth2hf \n");
0939   Diffe_Depth2hf->SetXTitle("#eta \b");
0940   Diffe_Depth2hf->SetYTitle("#phi \b");
0941   Diffe_Depth2hf->SetZTitle("<ADCAmpl> smalle 20000000 - HF Depth2hf \b");
0942   Diffe_Depth2hf->SetMarkerColor(2);
0943   Diffe_Depth2hf->SetLineColor(2);
0944   //      Diffe_Depth2hf->SetMaximum(1.000);
0945   //      Diffe_Depth2hf->SetMinimum(0.0000001);
0946   Diffe_Depth2hf->Draw("COLZ");
0947 
0948   c1->cd(4);
0949   TH1F *diffADCAmpl_Depth2hf = new TH1F("diffADCAmpl_Depth2hf", "", 100, 0., 300.);
0950   //      TH1F* diffADCAmpl_Depth2hf = new TH1F("diffADCAmpl_Depth2hf","", 1000, 0.,1000.);
0951   TH2F *Cdepth2hffw = (TH2F *)Cdepth2hfff->Clone("diffADCAmpl_Depth2hf");
0952   //      TH2F* Cdepth2hffw = (TH2F*)twoddepth2hf1->Clone("diffADCAmpl_Depth2hf");
0953   for (int i = 1; i <= Cdepth2hffw->GetXaxis()->GetNbins(); i++) {
0954     for (int j = 1; j <= Cdepth2hffw->GetYaxis()->GetNbins(); j++) {
0955       if (Cdepth2hffw->GetBinContent(i, j) != 0) {
0956         double ccc1 = Cdepth2hffw->GetBinContent(i, j);
0957         diffADCAmpl_Depth2hf->Fill(ccc1);
0958       }
0959     }
0960   }
0961   gPad->SetLogy();
0962   diffADCAmpl_Depth2hf->SetMarkerStyle(20);
0963   diffADCAmpl_Depth2hf->SetMarkerSize(0.4);
0964   diffADCAmpl_Depth2hf->GetYaxis()->SetLabelSize(0.04);
0965   diffADCAmpl_Depth2hf->SetXTitle("<ADCAmpl> in each cell \b");
0966   diffADCAmpl_Depth2hf->SetMarkerColor(2);
0967   diffADCAmpl_Depth2hf->SetLineColor(2);
0968   diffADCAmpl_Depth2hf->Draw("");
0969 
0970   c1->Update();
0971 
0972   delete twoddepth2hf0;
0973   delete twoddepth2hf1;
0974   delete adepth2hf;
0975   delete Cdepth2hffz225;
0976   delete Cdepth2hfff;
0977   delete diffADCAmpl_Depth2hf;
0978   delete Diffe_Depth2hf;
0979 
0980   //=============================================================================================== 8
0981   //======================================================================
0982   //======================================================================
0983   //================
0984   //======================================================================
0985   c1->Clear();
0986   c1->Divide(2, 2);
0987 
0988   TH2F *twoddepth4ho1 = (TH2F *)hfile1->Get("h_mapDepth4ADCAmpl_HO");
0989   TH2F *twoddepth4ho0 = (TH2F *)hfile1->Get("h_mapDepth4_HO");
0990   twoddepth4ho1->Sumw2();
0991   twoddepth4ho0->Sumw2();
0992   //    if(twoddepth4ho0->IsA()->InheritsFrom("TH2F")){
0993   TH2F *Cdepth4hoff = (TH2F *)twoddepth4ho1->Clone("Cdepth4hoff");
0994   Cdepth4hoff->Divide(twoddepth4ho1, twoddepth4ho0, 1, 1, "B");
0995   Cdepth4hoff->Sumw2();
0996   //    }
0997   c1->cd(1);
0998   TH2F *twoedepth4ho1 = (TH2F *)hfile1->Get("h_mapDepth4ADCAmpl225_HO");
0999   TH2F *twoedepth4ho0 = (TH2F *)hfile1->Get("h_mapDepth4_HO");
1000   twoedepth4ho1->Sumw2();
1001   twoedepth4ho0->Sumw2();
1002   //    if(twoe0->IsA()->InheritsFrom("TH2F")){
1003   TH2F *Cdepth4hofz225 = (TH2F *)twoedepth4ho1->Clone("Cdepth4hofz225");
1004   Cdepth4hofz225->Divide(twoedepth4ho1, twoedepth4ho0, 1, 1, "B");
1005   Cdepth4hofz225->Sumw2();
1006   //    }
1007   for (int i = 1; i <= Cdepth4hofz225->GetXaxis()->GetNbins(); i++) {
1008     for (int j = 1; j <= Cdepth4hofz225->GetYaxis()->GetNbins(); j++) {
1009       double ccc1 = Cdepth4hofz225->GetBinContent(i, j);
1010       //    if(ccc1> 0.1) cout<<"HO ibin=     "<< i <<" jbin=     "<< j <<" Rate=     "<< ccc1 <<endl;
1011     }
1012   }
1013   gPad->SetGridy();
1014   gPad->SetGridx();
1015   gPad->SetLogz();
1016   Cdepth4hofz225->SetMarkerStyle(20);
1017   Cdepth4hofz225->SetMarkerSize(0.4);
1018   //    Cdepth4hofz225->GetYaxis()->SetLabelSize(0.04);
1019   Cdepth4hofz225->GetZaxis()->SetLabelSize(0.08);
1020   Cdepth4hofz225->SetXTitle("#eta \b");
1021   Cdepth4hofz225->SetYTitle("#phi \b");
1022   Cdepth4hofz225->SetZTitle("Rate for ADCAmpl in each event & cell - HO Depth4 \b");
1023   Cdepth4hofz225->SetMarkerColor(2);
1024   Cdepth4hofz225->SetLineColor(2);
1025   Cdepth4hofz225->SetMaximum(1.000);
1026   Cdepth4hofz225->SetMinimum(0.1);
1027   //      Cdepth4hofz225->SetMinimum(0.0001);
1028   Cdepth4hofz225->Draw("COLZ");
1029 
1030   c1->cd(2);
1031 
1032   TH1F *adepth4ho = (TH1F *)hfile1->Get("h_ADCAmpl_HO");
1033   //     TH1F *adepth4ho= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HO");
1034   //      TH1F *adepth4ho= (TH1F*)hfile1->Get("h_ADCAmplZoom_HO");
1035   gPad->SetLogy();
1036   //      gPad->SetLogx();
1037   adepth4ho->SetMarkerStyle(20);
1038   adepth4ho->SetMarkerSize(0.8);
1039   adepth4ho->GetYaxis()->SetLabelSize(0.04);
1040   adepth4ho->SetXTitle("ADCAmpl in each event & cell HO \b");
1041   adepth4ho->SetMarkerColor(2);
1042   adepth4ho->SetLineColor(2);
1043   adepth4ho->Draw("");
1044 
1045   c1->cd(3);
1046 
1047   ///////////////////////////////////////
1048   TH2F *Diffe_Depth4ho = (TH2F *)Cdepth4hoff->Clone("Diffe_Depth4ho");
1049   for (int i = 1; i <= Cdepth4hoff->GetXaxis()->GetNbins(); i++) {
1050     for (int j = 1; j <= Cdepth4hoff->GetYaxis()->GetNbins(); j++) {
1051       //    TH2F* Cdepth4hoff = (TH2F*)twoddepth4ho1->Clone("Cdepth4hoff");
1052       double ccc1 = Cdepth4hoff->GetBinContent(i, j);
1053       Diffe_Depth4ho->SetBinContent(i, j, 0.);
1054       //    if(ccc1> 1000.|| (ccc1>0.&&ccc1<27.)) cout<<"HO ibin=     "<< i <<" jbin=     "<< j <<" A=     "<< ccc1 <<endl;
1055       //    if(ccc1> 1000.|| (i==46&&j==5)|| (i==56&&j==13)) cout<<"HO ibin=     "<< i <<" jbin=     "<< j <<" A=     "<< ccc1 <<endl;
1056       //    if(ccc1 < 20000000.)  Diffe_Depth4ho->SetBinContent(i,j,ccc1);
1057       if (ccc1 > 160. || ccc1 < 250.)
1058         Diffe_Depth4ho->SetBinContent(i, j, ccc1);
1059       //        if(ccc1 > 250.|| ccc1<400.)  Diffe_Depth4ho->SetBinContent(i,j,ccc1);
1060       //        if(ccc1 > 500.)  Diffe_Depth4ho->SetBinContent(i,j,ccc1);
1061       //    Diffe_Depth4ho->SetBinContent(i,j,ccc1);
1062     }
1063   }
1064   gPad->SetGridy();
1065   gPad->SetGridx();
1066   gPad->SetLogz();
1067   Diffe_Depth4ho->SetMarkerStyle(20);
1068   Diffe_Depth4ho->SetMarkerSize(0.4);
1069   Diffe_Depth4ho->GetZaxis()->SetLabelSize(0.08);
1070   //    Diffe_Depth4ho->SetTitle("any Error, HO Depth4ho \n");
1071   Diffe_Depth4ho->SetXTitle("#eta \b");
1072   Diffe_Depth4ho->SetYTitle("#phi \b");
1073   Diffe_Depth4ho->SetZTitle("<ADCAmpl> smalle 20000000 - HO Depth4ho \b");
1074   Diffe_Depth4ho->SetMarkerColor(2);
1075   Diffe_Depth4ho->SetLineColor(2);
1076   //      Diffe_Depth4ho->SetMaximum(1.000);
1077   //      Diffe_Depth4ho->SetMinimum(0.0000001);
1078   Diffe_Depth4ho->Draw("COLZ");
1079 
1080   c1->cd(4);
1081   //    TH1F* diffADCAmpl_Depth4ho = new TH1F("diffADCAmpl_Depth4ho","", 100, 0.,3000.);
1082   TH1F *diffADCAmpl_Depth4ho = new TH1F("diffADCAmpl_Depth4ho", "", 250, 0., 1000.);
1083 
1084   //      TH2F* Cdepth4hofw = (TH2F*)Cdepth4hoff->Clone("diffADCAmpl_Depth4ho");
1085   //      TH2F* Cdepth4hofw = (TH2F*)twoddepth4ho1->Clone("diffADCAmpl_Depth4ho");
1086 
1087   for (int i = 1; i <= Cdepth4hoff->GetXaxis()->GetNbins(); i++) {
1088     for (int j = 1; j <= Cdepth4hoff->GetYaxis()->GetNbins(); j++) {
1089       if (Cdepth4hoff->GetBinContent(i, j) != 0) {
1090         double ccc1 = Cdepth4hoff->GetBinContent(i, j);
1091         diffADCAmpl_Depth4ho->Fill(ccc1);
1092       }
1093     }
1094   }
1095   gPad->SetLogy();
1096   diffADCAmpl_Depth4ho->SetMarkerStyle(20);
1097   diffADCAmpl_Depth4ho->SetMarkerSize(0.4);
1098   diffADCAmpl_Depth4ho->GetYaxis()->SetLabelSize(0.04);
1099   diffADCAmpl_Depth4ho->SetXTitle("<ADCAmpl> in each cell \b");
1100   diffADCAmpl_Depth4ho->SetMarkerColor(2);
1101   diffADCAmpl_Depth4ho->SetLineColor(2);
1102   diffADCAmpl_Depth4ho->SetMinimum(0.8);
1103   diffADCAmpl_Depth4ho->Draw("");
1104 
1105   c1->Update();
1106 
1107   delete twoddepth4ho0;
1108   delete twoddepth4ho1;
1109   delete adepth4ho;
1110   delete Cdepth4hofz225;
1111   delete Cdepth4hoff;
1112   delete diffADCAmpl_Depth4ho;
1113   delete Diffe_Depth4ho;
1114 
1115   //========================================================================================== 9
1116   //======================================================================
1117   //======================================================================
1118   //================
1119   /*
1120     // fullAmplitude:
1121 ///////////////////////////////////////////////////////////////////////////////////////
1122     h_ADCAmpl345Zoom_HB = new TH1F("h_ADCAmpl345Zoom_HB"," ", 100, 0.,400.);
1123     h_ADCAmpl345Zoom1_HB = new TH1F("h_ADCAmpl345Zoom1_HB"," ", 100, 0.,100.);
1124     h_ADCAmpl345_HB = new TH1F("h_ADCAmpl345_HB"," ", 100, 10.,3000.);
1125 
1126     h_ADCAmplZoom_HB = new TH1F("h_ADCAmplZoom_HB"," ", 100, 0.,400.);
1127     h_ADCAmplZoom1_HB = new TH1F("h_ADCAmplZoom1_HB"," ", 100, -20.,80.);
1128     h_ADCAmpl_HB = new TH1F("h_ADCAmpl_HB"," ", 100, 10.,5000.);
1129 
1130     h_AmplitudeHBrest = new TH1F("h_AmplitudeHBrest"," ", 100, 0.,10000.);
1131     h_AmplitudeHBrest1 = new TH1F("h_AmplitudeHBrest1"," ", 100, 0.,1000000.);
1132     h_AmplitudeHBrest6 = new TH1F("h_AmplitudeHBrest6"," ", 100, 0.,2000000.);
1133 */
1134   //======================================================================
1135   c1->Clear();
1136   c1->Divide(2, 3);
1137 
1138   c1->cd(1);
1139   TH1F *aaaaaa0 = (TH1F *)hfile1->Get("h_ADCAmplZoom_HB");
1140   gPad->SetLogy();
1141   // gPad->SetLogx();
1142   aaaaaa0->SetMarkerStyle(20);
1143   aaaaaa0->SetMarkerSize(0.8);
1144   aaaaaa0->GetYaxis()->SetLabelSize(0.04);
1145   aaaaaa0->SetXTitle("h_ADCAmplZoom_HB \b");
1146   aaaaaa0->SetMarkerColor(2);
1147   aaaaaa0->SetLineColor(2);
1148   aaaaaa0->Draw("");
1149 
1150   c1->cd(2);
1151   TH1F *aaaaaa3 = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HB");
1152   gPad->SetLogy();
1153   //    gPad->SetLogx();
1154   aaaaaa3->SetMarkerStyle(20);
1155   aaaaaa3->SetMarkerSize(0.8);
1156   aaaaaa3->GetYaxis()->SetLabelSize(0.04);
1157   aaaaaa3->SetXTitle("ADCAmpl in each event & cell HB \b");
1158   aaaaaa3->SetMarkerColor(2);
1159   aaaaaa3->SetLineColor(2);
1160   aaaaaa3->Draw("");
1161 
1162   c1->cd(3);
1163   TH1F *aaaaaa5 = (TH1F *)hfile1->Get("h_ADCAmpl_HB");
1164   gPad->SetLogy();
1165   //    gPad->SetLogx();
1166   aaaaaa5->SetMarkerStyle(20);
1167   aaaaaa5->SetMarkerSize(0.8);
1168   aaaaaa5->GetYaxis()->SetLabelSize(0.04);
1169   aaaaaa5->SetXTitle("ADCAmpl in each event & cell HB \b");
1170   aaaaaa5->SetMarkerColor(2);
1171   aaaaaa5->SetLineColor(2);
1172   aaaaaa5->Draw("");
1173 
1174   c1->cd(4);
1175   TH1F *aaaaaa1 = (TH1F *)hfile1->Get("h_AmplitudeHBrest");
1176   gPad->SetLogy();
1177   // gPad->SetLogx();
1178   aaaaaa1->SetMarkerStyle(20);
1179   aaaaaa1->SetMarkerSize(0.8);
1180   aaaaaa1->GetYaxis()->SetLabelSize(0.04);
1181   aaaaaa1->SetXTitle("h_AmplitudeHBrest \b");
1182   aaaaaa1->SetMarkerColor(2);
1183   aaaaaa1->SetLineColor(2);
1184   aaaaaa1->Draw("");
1185 
1186   c1->cd(5);
1187   TH1F *aaaaaa2 = (TH1F *)hfile1->Get("h_AmplitudeHBrest1");
1188   gPad->SetLogy();
1189   // gPad->SetLogx();
1190   aaaaaa2->SetMarkerStyle(20);
1191   aaaaaa2->SetMarkerSize(0.8);
1192   aaaaaa2->GetYaxis()->SetLabelSize(0.04);
1193   aaaaaa2->SetXTitle("h_AmplitudeHBrest1 \b");
1194   aaaaaa2->SetMarkerColor(2);
1195   aaaaaa2->SetLineColor(2);
1196   aaaaaa2->Draw("");
1197 
1198   c1->cd(6);
1199   TH1F *aaaaaa4 = (TH1F *)hfile1->Get("h_AmplitudeHBrest6");
1200   gPad->SetLogy();
1201   // gPad->SetLogx();
1202   aaaaaa4->SetMarkerStyle(20);
1203   aaaaaa4->SetMarkerSize(0.8);
1204   aaaaaa4->GetYaxis()->SetLabelSize(0.04);
1205   aaaaaa4->SetXTitle("h_AmplitudeHBrest6 \b");
1206   aaaaaa4->SetMarkerColor(2);
1207   aaaaaa4->SetLineColor(2);
1208   aaaaaa4->Draw("");
1209 
1210   c1->Update();
1211 
1212   //========================================================================================== 10
1213   //======================================================================
1214   //======================================================================
1215   //================
1216   /*
1217     // fullAmplitude:
1218 ///////////////////////////////////////////////////////////////////////////////////////
1219 
1220     h_ADCAmpl_HE = new TH1F("h_ADCAmpl_HE"," ", 200, 0.,2000000.);
1221     h_ADCAmpl345_HE = new TH1F("h_ADCAmpl345_HE"," ", 70, 0.,700000.);
1222 
1223 // for SiPM:
1224 const int npfit = 220; float anpfit = 220.;
1225     h_ADCAmplZoom1_HE = new TH1F("h_ADCAmplZoom1_HE"," ",npfit, 0.,anpfit);// for amplmaxts 1TS w/ max
1226     h_ADCAmpl345Zoom1_HE = new TH1F("h_ADCAmpl345Zoom1_HE"," ", npfit, 0.,anpfit);// for ampl3ts 3TSs
1227     h_ADCAmpl345Zoom_HE = new TH1F("h_ADCAmpl345Zoom_HE"," ", npfit, 0.,anpfit); // for ampl 4TSs
1228 */
1229   //======================================================================
1230   c1->Clear();
1231   c1->Divide(2, 3);
1232 
1233   c1->cd(1);
1234   TH1F *aaaaab0 = (TH1F *)hfile1->Get("h_ADCAmpl345_HE");
1235   gPad->SetLogy();
1236   // gPad->SetLogx();
1237   aaaaab0->SetMarkerStyle(20);
1238   aaaaab0->SetMarkerSize(0.8);
1239   aaaaab0->GetYaxis()->SetLabelSize(0.04);
1240   aaaaab0->SetXTitle("h_ADCAmpl345_HE \b");
1241   aaaaab0->SetMarkerColor(2);
1242   aaaaab0->SetLineColor(2);
1243   aaaaab0->Draw("");
1244 
1245   c1->cd(2);
1246   TH1F *aaaaab3 = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HE");
1247   gPad->SetLogy();
1248   //    gPad->SetLogx();
1249   aaaaab3->SetMarkerStyle(20);
1250   aaaaab3->SetMarkerSize(0.8);
1251   aaaaab3->GetYaxis()->SetLabelSize(0.04);
1252   aaaaab3->SetXTitle("for amplmaxts 1TS w/ max HE \b");
1253   aaaaab3->SetMarkerColor(2);
1254   aaaaab3->SetLineColor(2);
1255   aaaaab3->Draw("");
1256   TH1F *aaaaab4 = (TH1F *)hfile1->Get("h_ADCAmpl345Zoom1_HE");
1257   gPad->SetLogy();
1258   //    gPad->SetLogx();
1259   aaaaab4->SetMarkerStyle(20);
1260   aaaaab4->SetMarkerSize(0.8);
1261   aaaaab4->GetYaxis()->SetLabelSize(0.04);
1262   aaaaab4->SetMarkerColor(4);
1263   aaaaab4->SetLineColor(4);
1264   aaaaab4->Draw("Same");
1265 
1266   c1->cd(3);
1267   TH1F *aaaaab5 = (TH1F *)hfile1->Get("h_ADCAmpl_HE");
1268   gPad->SetLogy();
1269   //    gPad->SetLogx();
1270   aaaaab5->SetMarkerStyle(20);
1271   aaaaab5->SetMarkerSize(0.8);
1272   aaaaab5->GetYaxis()->SetLabelSize(0.04);
1273   aaaaab5->SetXTitle("ADCAmpl in each event & cell HE \b");
1274   aaaaab5->SetMarkerColor(2);
1275   aaaaab5->SetLineColor(2);
1276   aaaaab5->Draw("");
1277 
1278   c1->cd(4);
1279   TH1F *aaaaab1 = (TH1F *)hfile1->Get("h_ADCAmplrest_HE");
1280   //    TH1F *aaaaab1= (TH1F*)hfile1->Get("h_AmplitudeHEtest");
1281   gPad->SetLogy();
1282   // gPad->SetLogx();
1283   aaaaab1->SetMarkerStyle(20);
1284   aaaaab1->SetMarkerSize(0.8);
1285   aaaaab1->GetYaxis()->SetLabelSize(0.04);
1286   aaaaab1->SetXTitle("h_ADCAmplrest_HE \b");
1287   aaaaab1->SetMarkerColor(2);
1288   aaaaab1->SetLineColor(2);
1289   aaaaab1->Draw("");
1290 
1291   c1->cd(5);
1292   TH1F *aaaaab2 = (TH1F *)hfile1->Get("h_ADCAmplrest1_HE");
1293   //    TH1F *aaaaab2= (TH1F*)hfile1->Get("h_AmplitudeHEtest1");
1294   gPad->SetLogy();
1295   // gPad->SetLogx();
1296   aaaaab2->SetMarkerStyle(20);
1297   aaaaab2->SetMarkerSize(0.8);
1298   aaaaab2->GetYaxis()->SetLabelSize(0.04);
1299   aaaaab2->SetXTitle("h_ADCAmplrest1_HE \b");
1300   aaaaab2->SetMarkerColor(2);
1301   aaaaab2->SetLineColor(2);
1302   aaaaab2->Draw("");
1303 
1304   c1->cd(6);
1305   TH1F *aaaaab6 = (TH1F *)hfile1->Get("h_ADCAmplrest6_HE");
1306   //    TH1F *aaaaab6= (TH1F*)hfile1->Get("h_AmplitudeHEtest6");
1307   gPad->SetLogy();
1308   // gPad->SetLogx();
1309   aaaaab6->SetMarkerStyle(20);
1310   aaaaab6->SetMarkerSize(0.8);
1311   aaaaab6->GetYaxis()->SetLabelSize(0.04);
1312   aaaaab6->SetXTitle("h_ADCAmplrest6_HE \b");
1313   aaaaab6->SetMarkerColor(2);
1314   aaaaab6->SetLineColor(2);
1315   aaaaab6->Draw("");
1316 
1317   c1->Update();
1318 
1319   //========================================================================================== 11
1320   //======================================================================
1321   //======================================================================
1322   //================
1323   /*
1324     // fullAmplitude:
1325 ///////////////////////////////////////////////////////////////////////////////////////
1326     h_ADCAmplZoom1_HF = new TH1F("h_ADCAmplZoom1_HF"," ", 100, 0.,1000000.);
1327     h_ADCAmpl_HF = new TH1F("h_ADCAmpl_HF"," ", 250, 0.,500000.);
1328 
1329 */
1330   //======================================================================
1331   c1->Clear();
1332   c1->Divide(2, 2);
1333 
1334   c1->cd(1);
1335   TH1F *aaaaac0 = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HF");
1336   gPad->SetLogy();
1337   // gPad->SetLogx();
1338   aaaaac0->SetMarkerStyle(20);
1339   aaaaac0->SetMarkerSize(0.8);
1340   aaaaac0->GetYaxis()->SetLabelSize(0.04);
1341   aaaaac0->SetXTitle("h_ADCAmplZoom1_HF \b");
1342   aaaaac0->SetMarkerColor(2);
1343   aaaaac0->SetLineColor(2);
1344   aaaaac0->Draw("");
1345 
1346   c1->cd(2);
1347   TH1F *aaaaac3 = (TH1F *)hfile1->Get("h_ADCAmpl_HF");
1348   gPad->SetLogy();
1349   //    gPad->SetLogx();
1350   aaaaac3->SetMarkerStyle(20);
1351   aaaaac3->SetMarkerSize(0.8);
1352   aaaaac3->GetYaxis()->SetLabelSize(0.04);
1353   aaaaac3->SetXTitle("h_ADCAmpl_HF \b");
1354   aaaaac3->SetMarkerColor(2);
1355   aaaaac3->SetLineColor(2);
1356   aaaaac3->Draw("");
1357 
1358   c1->cd(3);
1359   TH1F *aaaaac5 = (TH1F *)hfile1->Get("h_ADCAmplrest1_HF");
1360   gPad->SetLogy();
1361   //    gPad->SetLogx();
1362   aaaaac5->SetMarkerStyle(20);
1363   aaaaac5->SetMarkerSize(0.8);
1364   aaaaac5->GetYaxis()->SetLabelSize(0.04);
1365   aaaaac5->SetXTitle("h_ADCAmplrest1_HF  \b");
1366   aaaaac5->SetMarkerColor(2);
1367   aaaaac5->SetLineColor(2);
1368   aaaaac5->Draw("");
1369 
1370   c1->cd(4);
1371   TH1F *aaaaac1 = (TH1F *)hfile1->Get("h_ADCAmplrest6_HF");
1372   gPad->SetLogy();
1373   //    gPad->SetLogx();
1374   aaaaac1->SetMarkerStyle(20);
1375   aaaaac1->SetMarkerSize(0.8);
1376   aaaaac1->GetYaxis()->SetLabelSize(0.04);
1377   aaaaac1->SetXTitle("h_ADCAmplrest6_HF \b");
1378   aaaaac1->SetMarkerColor(2);
1379   aaaaac1->SetLineColor(2);
1380   aaaaac1->Draw("");
1381 
1382   c1->Update();
1383 
1384   //========================================================================================== 12
1385   //======================================================================
1386   //======================================================================
1387   //================
1388   /*
1389     // fullAmplitude:
1390 ///////////////////////////////////////////////////////////////////////////////////////
1391     h_ADCAmpl_HO = new TH1F("h_ADCAmpl_HO"," ", 100, 0.,7000.);
1392     h_ADCAmplZoom1_HO = new TH1F("h_ADCAmplZoom1_HO"," ", 100, -20.,280.);
1393     h_ADCAmpl_HO_copy = new TH1F("h_ADCAmpl_HO_copy"," ", 100, 0.,30000.);
1394 */
1395   //======================================================================
1396   c1->Clear();
1397   c1->Divide(2, 3);
1398 
1399   c1->cd(1);
1400   TH1F *aaaaad0 = (TH1F *)hfile1->Get("h_ADCAmplrest1_HO");
1401   gPad->SetLogy();
1402   // gPad->SetLogx();
1403   aaaaad0->SetMarkerStyle(20);
1404   aaaaad0->SetMarkerSize(0.8);
1405   aaaaad0->GetYaxis()->SetLabelSize(0.04);
1406   aaaaad0->SetXTitle("h_ADCAmplrest1_HO \b");
1407   aaaaad0->SetMarkerColor(2);
1408   aaaaad0->SetLineColor(2);
1409   aaaaad0->Draw("");
1410 
1411   c1->cd(2);
1412   TH1F *aaaaad3 = (TH1F *)hfile1->Get("h_ADCAmplrest6_HO");
1413   gPad->SetLogy();
1414   //    gPad->SetLogx();
1415   aaaaad3->SetMarkerStyle(20);
1416   aaaaad3->SetMarkerSize(0.8);
1417   aaaaad3->GetYaxis()->SetLabelSize(0.04);
1418   aaaaad3->SetXTitle("h_ADCAmplrest6_HO \b");
1419   aaaaad3->SetMarkerColor(2);
1420   aaaaad3->SetLineColor(2);
1421   aaaaad3->Draw("");
1422 
1423   c1->cd(3);
1424   TH1F *aaaaad5 = (TH1F *)hfile1->Get("h_ADCAmpl_HO");
1425   gPad->SetLogy();
1426   //    gPad->SetLogx();
1427   aaaaad5->SetMarkerStyle(20);
1428   aaaaad5->SetMarkerSize(0.8);
1429   aaaaad5->GetYaxis()->SetLabelSize(0.04);
1430   aaaaad5->SetXTitle("h_ADCAmpl_HO \b");
1431   aaaaad5->SetMarkerColor(2);
1432   aaaaad5->SetLineColor(2);
1433   aaaaad5->Draw("");
1434 
1435   c1->cd(4);
1436   TH1F *aaaaad1 = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HO");
1437   gPad->SetLogy();
1438   //    gPad->SetLogx();
1439   aaaaad1->SetMarkerStyle(20);
1440   aaaaad1->SetMarkerSize(0.8);
1441   aaaaad1->GetYaxis()->SetLabelSize(0.04);
1442   aaaaad1->SetXTitle("h_ADCAmplZoom1_HO \b");
1443   aaaaad1->SetMarkerColor(2);
1444   aaaaad1->SetLineColor(2);
1445   aaaaad1->Draw("");
1446 
1447   c1->cd(5);
1448   TH1F *aaaaad2 = (TH1F *)hfile1->Get("h_ADCAmpl_HO_copy");
1449   gPad->SetLogy();
1450   //    gPad->SetLogx();
1451   aaaaad2->SetMarkerStyle(20);
1452   aaaaad2->SetMarkerSize(0.8);
1453   aaaaad2->GetYaxis()->SetLabelSize(0.04);
1454   aaaaad2->SetXTitle("h_ADCAmpl_HO_copy \b");
1455   aaaaad2->SetMarkerColor(2);
1456   aaaaad2->SetLineColor(2);
1457   aaaaad2->Draw("");
1458 
1459   c1->Update();
1460 
1461   //========================================================================================== 13
1462   //======================================================================
1463   //======================================================================
1464   c1->Clear();
1465   c1->Divide(2, 2);
1466 
1467   TH1F *time = (TH1F *)hfile1->Get("h_tdc_HE_time");
1468   TH1F *time0 = (TH1F *)hfile1->Get("h_tdc_HE_time0");
1469   TH1F *timer = (TH1F *)time->Clone("timer");
1470   time->Sumw2();
1471   time0->Sumw2();
1472   timer->Divide(time, time0, 1, 1, "B");  //
1473   timer->Sumw2();
1474   TH1F *timet = (TH1F *)hfile1->Get("h_tdc_HE_tdc");
1475   c1->cd(1);
1476   //            gPad->SetLogy();
1477   //    gPad->SetLogx();
1478   time->SetMarkerStyle(20);
1479   time->SetMarkerSize(0.4);
1480   time->GetYaxis()->SetLabelSize(0.04);
1481   time->SetXTitle("time50 weighted by A  HE \b");
1482   time->SetMarkerColor(2);
1483   time->SetLineColor(2);
1484   time->Draw("");
1485   c1->cd(2);
1486   time0->SetMarkerStyle(20);
1487   time0->SetMarkerSize(0.4);
1488   time0->GetYaxis()->SetLabelSize(0.04);
1489   time0->SetXTitle("time50 HE \b");
1490   time0->SetMarkerColor(2);
1491   time0->SetLineColor(2);
1492   time0->Draw("");
1493   c1->cd(3);
1494   timer->SetMarkerStyle(20);
1495   timer->SetMarkerSize(0.4);
1496   timer->GetYaxis()->SetLabelSize(0.04);
1497   timer->SetXTitle("shape HE \b");
1498   timer->SetMarkerColor(2);
1499   timer->SetLineColor(0);
1500   timer->SetMinimum(0.);
1501   timer->SetMaximum(50000.);
1502   timer->Draw("");
1503   c1->cd(4);
1504   gPad->SetLogy();
1505   timet->SetMarkerStyle(20);
1506   timet->SetMarkerSize(0.4);
1507   timet->GetYaxis()->SetLabelSize(0.04);
1508   timet->SetXTitle("initial tdc HE \b");
1509   timet->SetMarkerColor(2);
1510   timet->SetLineColor(2);
1511   timet->Draw("");
1512 
1513   c1->Update();
1514   //======================================================================
1515   //========================================================================================== 14
1516   //======================================================================
1517   //======================================================================
1518   c1->Clear();
1519   c1->Divide(1, 3);
1520 
1521   TH1F *ampldefault = (TH1F *)hfile1->Get("h_tdc_HE_ampldefault");
1522   TH1F *ampldefault50 = (TH1F *)hfile1->Get("h_tdc_HE_ampldefault50");
1523   TH1F *ampldefault63 = (TH1F *)hfile1->Get("h_tdc_HE_ampldefault63");
1524   c1->cd(1);
1525   //            gPad->SetLogy();
1526   //    gPad->SetLogx();
1527   ampldefault->SetMarkerStyle(20);
1528   ampldefault->SetMarkerSize(0.4);
1529   ampldefault->GetYaxis()->SetLabelSize(0.04);
1530   ampldefault->SetXTitle("A_TS HE \b");
1531   ampldefault->SetMarkerColor(2);
1532   ampldefault->SetLineColor(2);
1533   gPad->SetLogy();
1534   ampldefault->Draw("");
1535 
1536   c1->cd(2);
1537   //            gPad->SetLogy();
1538   //    gPad->SetLogx();
1539   ampldefault50->SetMarkerStyle(20);
1540   ampldefault50->SetMarkerSize(0.4);
1541   ampldefault50->GetYaxis()->SetLabelSize(0.04);
1542   ampldefault50->SetXTitle("A_TS_50 HE \b");
1543   ampldefault50->SetMarkerColor(2);
1544   ampldefault50->SetLineColor(2);
1545   gPad->SetLogy();
1546   ampldefault50->Draw("");
1547 
1548   c1->cd(3);
1549   //            gPad->SetLogy();
1550   //    gPad->SetLogx();
1551   ampldefault63->SetMarkerStyle(20);
1552   ampldefault63->SetMarkerSize(0.4);
1553   ampldefault63->GetYaxis()->SetLabelSize(0.04);
1554   ampldefault63->SetXTitle("A_TS_63 HE \b");
1555   ampldefault63->SetMarkerColor(2);
1556   ampldefault63->SetLineColor(2);
1557   gPad->SetLogy();
1558   ampldefault63->Draw("");
1559 
1560   c1->Update();
1561   //======================================================================
1562   //========================================================================================== 15
1563   //======================================================================
1564   //======================================================================
1565   c1->Clear();
1566   c1->Divide(1, 1);
1567 
1568   TH1F *timeVSampldefault = (TH1F *)hfile1->Get("h_tdc_HE_timeVSampldefault");
1569   c1->cd(1);
1570   //            gPad->SetLogy();
1571   //    gPad->SetLogx();
1572   timeVSampldefault->SetMarkerStyle(20);
1573   timeVSampldefault->SetMarkerSize(0.4);
1574   timeVSampldefault->GetYaxis()->SetLabelSize(0.04);
1575   timeVSampldefault->SetXTitle("timeVStampldefault HE \b");
1576   timeVSampldefault->SetMarkerColor(2);
1577   timeVSampldefault->SetLineColor(2);
1578   timeVSampldefault->Draw("box");
1579 
1580   c1->Update();
1581   //========================================================================================== 16
1582   //======================================================================
1583   //======================================================================
1584   c1->Clear();
1585   c1->Divide(2, 3);
1586 
1587   TH1F *shape = (TH1F *)hfile1->Get("h_shape_good_channels_HE");
1588   TH1F *shape0 = (TH1F *)hfile1->Get("h_shape0_good_channels_HE");
1589   TH1F *shaper = (TH1F *)shape->Clone("shaper");
1590   shape->Sumw2();
1591   shape0->Sumw2();
1592   shaper->Divide(shape, shape0, 1, 1, "B");  //
1593   shaper->Sumw2();
1594   c1->cd(1);
1595   //            gPad->SetLogy();
1596   //    gPad->SetLogx();
1597   shape->SetMarkerStyle(20);
1598   shape->SetMarkerSize(0.4);
1599   shape->GetYaxis()->SetLabelSize(0.04);
1600   shape->SetXTitle("TS weighted by A good HE \b");
1601   shape->SetMarkerColor(2);
1602   shape->SetLineColor(2);
1603   shape->Draw("");
1604   c1->cd(3);
1605   shape0->SetMarkerStyle(20);
1606   shape0->SetMarkerSize(0.4);
1607   shape0->GetYaxis()->SetLabelSize(0.04);
1608   shape0->SetXTitle("TS good HE \b");
1609   shape0->SetMarkerColor(2);
1610   shape0->SetLineColor(2);
1611   shape0->Draw("");
1612   c1->cd(5);
1613   shaper->SetMarkerStyle(20);
1614   shaper->SetMarkerSize(0.4);
1615   shaper->GetYaxis()->SetLabelSize(0.04);
1616   shaper->SetXTitle("shape good HE per event, per channel\b");
1617   shaper->SetMarkerColor(2);
1618   shaper->SetLineColor(2);
1619   shaper->Draw("");
1620 
1621   TH1F *badsh = (TH1F *)hfile1->Get("h_shape_bad_channels_HE");
1622   TH1F *badsh0 = (TH1F *)hfile1->Get("h_shape0_bad_channels_HE");
1623   TH1F *badshr = (TH1F *)badsh->Clone("badshr");
1624   badsh->Sumw2();
1625   badsh0->Sumw2();
1626   badshr->Divide(badsh, badsh0, 1, 1, "B");  //
1627   badshr->Sumw2();
1628   c1->cd(2);
1629   //            gPad->SetLogy();
1630   //    gPad->SetLogx();
1631   badsh->SetMarkerStyle(20);
1632   badsh->SetMarkerSize(0.4);
1633   badsh->GetYaxis()->SetLabelSize(0.04);
1634   badsh->SetXTitle("TS weighted by A bad HE \b");
1635   badsh->SetMarkerColor(2);
1636   badsh->SetLineColor(2);
1637   badsh->Draw("");
1638   c1->cd(4);
1639   badsh0->SetMarkerStyle(20);
1640   badsh0->SetMarkerSize(0.4);
1641   badsh0->GetYaxis()->SetLabelSize(0.04);
1642   badsh0->SetXTitle("TS bad HE \b");
1643   badsh0->SetMarkerColor(2);
1644   badsh0->SetLineColor(2);
1645   badsh0->Draw("");
1646   c1->cd(6);
1647   badshr->SetMarkerStyle(20);
1648   badshr->SetMarkerSize(0.4);
1649   badshr->GetYaxis()->SetLabelSize(0.04);
1650   badshr->SetXTitle("shape bad HE per event, per channel\b");
1651   badshr->SetMarkerColor(2);
1652   badshr->SetLineColor(2);
1653   badshr->Draw("");
1654 
1655   c1->Update();
1656   //======================================================================
1657   //========================================================================================== 17
1658   //======================================================================
1659   //======================================================================
1660   c1->Clear();
1661   c1->Divide(2, 2);
1662 
1663   TH1F *timeHFt = (TH1F *)hfile1->Get("h_tdc_HF_tdc");
1664   TH1F *timeHF = (TH1F *)hfile1->Get("h_tdc_HF_time");
1665   TH1F *timeHF0 = (TH1F *)hfile1->Get("h_tdc_HF_time0");
1666   TH1F *timeHFr = (TH1F *)timeHF->Clone("timeHFr");
1667   timeHF->Sumw2();
1668   timeHF0->Sumw2();
1669   timeHFr->Divide(timeHF, timeHF0, 1, 1, "B");  //
1670   timeHFr->Sumw2();
1671   c1->cd(1);
1672   //            gPad->SetLogy();
1673   //    gPad->SetLogx();
1674   timeHF->SetMarkerStyle(20);
1675   timeHF->SetMarkerSize(0.4);
1676   timeHF->GetYaxis()->SetLabelSize(0.04);
1677   timeHF->SetXTitle("time50 weighted by A  HF \b");
1678   timeHF->SetMarkerColor(2);
1679   timeHF->SetLineColor(0);
1680   timeHF->Draw("");
1681   c1->cd(2);
1682   timeHF0->SetMarkerStyle(20);
1683   timeHF0->SetMarkerSize(0.4);
1684   timeHF0->GetYaxis()->SetLabelSize(0.04);
1685   timeHF0->SetXTitle("time50 HF \b");
1686   timeHF0->SetMarkerColor(2);
1687   timeHF0->SetLineColor(0);
1688   timeHF0->Draw("");
1689   c1->cd(3);
1690   timeHFr->SetMarkerStyle(20);
1691   timeHFr->SetMarkerSize(0.4);
1692   timeHFr->GetYaxis()->SetLabelSize(0.04);
1693   timeHFr->SetXTitle("shape HF \b");
1694   timeHFr->SetMarkerColor(2);
1695   timeHFr->SetLineColor(0);
1696   timeHFr->SetMinimum(0.);
1697   timeHFr->SetMaximum(1000.);
1698   timeHFr->Draw("");
1699   c1->cd(4);
1700   gPad->SetLogy();
1701   timeHFt->SetMarkerStyle(20);
1702   timeHFt->SetMarkerSize(0.4);
1703   timeHFt->GetYaxis()->SetLabelSize(0.04);
1704   timeHFt->SetXTitle("initial tdc HF \b");
1705   timeHFt->SetMarkerColor(2);
1706   timeHFt->SetLineColor(2);
1707   timeHFt->Draw("");
1708 
1709   c1->Update();
1710   //======================================================================
1711   //========================================================================================== 18
1712   //======================================================================
1713   //======================================================================
1714   c1->Clear();
1715   c1->Divide(1, 3);
1716 
1717   TH1F *ampldefaultHF = (TH1F *)hfile1->Get("h_tdc_HF_ampldefault");
1718   TH1F *ampldefaultHF50 = (TH1F *)hfile1->Get("h_tdc_HF_ampldefault50");
1719   TH1F *ampldefaultHF63 = (TH1F *)hfile1->Get("h_tdc_HF_ampldefault63");
1720   c1->cd(1);
1721   //            gPad->SetLogy();
1722   //    gPad->SetLogx();
1723   ampldefaultHF->SetMarkerStyle(20);
1724   ampldefaultHF->SetMarkerSize(0.4);
1725   ampldefaultHF->GetYaxis()->SetLabelSize(0.04);
1726   ampldefaultHF->SetXTitle("A_TS HF \b");
1727   ampldefaultHF->SetMarkerColor(2);
1728   ampldefaultHF->SetLineColor(2);
1729   gPad->SetLogy();
1730   ampldefaultHF->Draw("");
1731 
1732   c1->cd(2);
1733   //            gPad->SetLogy();
1734   //    gPad->SetLogx();
1735   ampldefaultHF50->SetMarkerStyle(20);
1736   ampldefaultHF50->SetMarkerSize(0.4);
1737   ampldefaultHF50->GetYaxis()->SetLabelSize(0.04);
1738   ampldefaultHF50->SetXTitle("A_TS_50 HF \b");
1739   ampldefaultHF50->SetMarkerColor(2);
1740   ampldefaultHF50->SetLineColor(2);
1741   gPad->SetLogy();
1742   ampldefaultHF50->Draw("");
1743 
1744   c1->cd(3);
1745   //            gPad->SetLogy();
1746   //    gPad->SetLogx();
1747   ampldefaultHF63->SetMarkerStyle(20);
1748   ampldefaultHF63->SetMarkerSize(0.4);
1749   ampldefaultHF63->GetYaxis()->SetLabelSize(0.04);
1750   ampldefaultHF63->SetXTitle("A_TS_63 HF \b");
1751   ampldefaultHF63->SetMarkerColor(2);
1752   ampldefaultHF63->SetLineColor(2);
1753   gPad->SetLogy();
1754   ampldefaultHF63->Draw("");
1755 
1756   c1->Update();
1757   //======================================================================
1758   //========================================================================================== 19
1759   //======================================================================
1760   //======================================================================
1761   c1->Clear();
1762   c1->Divide(1, 1);
1763 
1764   TH1F *timeHFVSampldefault = (TH1F *)hfile1->Get("h_tdc_HF_timeVSampldefault");
1765   c1->cd(1);
1766   //            gPad->SetLogy();
1767   //    gPad->SetLogx();
1768   timeHFVSampldefault->SetMarkerStyle(20);
1769   timeHFVSampldefault->SetMarkerSize(0.4);
1770   timeHFVSampldefault->GetYaxis()->SetLabelSize(0.04);
1771   timeHFVSampldefault->SetXTitle("timeVStampldefault HF \b");
1772   timeHFVSampldefault->SetMarkerColor(2);
1773   timeHFVSampldefault->SetLineColor(2);
1774   timeHFVSampldefault->Draw("box");
1775 
1776   c1->Update();
1777   //========================================================================================== 20
1778   //======================================================================
1779   //======================================================================
1780   c1->Clear();
1781   c1->Divide(2, 3);
1782 
1783   TH1F *shapeHF = (TH1F *)hfile1->Get("h_shape_good_channels_HF");
1784   TH1F *shapeHF0 = (TH1F *)hfile1->Get("h_shape0_good_channels_HF");
1785   TH1F *shapeHFr = (TH1F *)shapeHF->Clone("shapeHFr");
1786   shapeHF->Sumw2();
1787   shapeHF0->Sumw2();
1788   shapeHFr->Divide(shapeHF, shapeHF0, 1, 1, "B");  //
1789   shapeHFr->Sumw2();
1790   c1->cd(1);
1791   //            gPad->SetLogy();
1792   //    gPad->SetLogx();
1793   shapeHF->SetMarkerStyle(20);
1794   shapeHF->SetMarkerSize(0.4);
1795   shapeHF->GetYaxis()->SetLabelSize(0.04);
1796   shapeHF->SetXTitle("TS weighted by A good HF \b");
1797   shapeHF->SetMarkerColor(2);
1798   shapeHF->SetLineColor(2);
1799   shapeHF->Draw("");
1800   c1->cd(3);
1801   shapeHF0->SetMarkerStyle(20);
1802   shapeHF0->SetMarkerSize(0.4);
1803   shapeHF0->GetYaxis()->SetLabelSize(0.04);
1804   shapeHF0->SetXTitle("TS good HF \b");
1805   shapeHF0->SetMarkerColor(2);
1806   shapeHF0->SetLineColor(2);
1807   shapeHF0->Draw("");
1808   c1->cd(5);
1809   shapeHFr->SetMarkerStyle(20);
1810   shapeHFr->SetMarkerSize(0.4);
1811   shapeHFr->GetYaxis()->SetLabelSize(0.04);
1812   shapeHFr->SetXTitle("shape good HF per event, per channel\b");
1813   shapeHFr->SetMarkerColor(2);
1814   shapeHFr->SetLineColor(2);
1815   shapeHFr->Draw("");
1816 
1817   TH1F *badshHF = (TH1F *)hfile1->Get("h_shape_bad_channels_HF");
1818   TH1F *badshHF0 = (TH1F *)hfile1->Get("h_shape0_bad_channels_HF");
1819   TH1F *badshHFr = (TH1F *)badshHF->Clone("badshHFr");
1820   badshHF->Sumw2();
1821   badshHF0->Sumw2();
1822   badshHFr->Divide(badshHF, badshHF0, 1, 1, "B");  //
1823   badshHFr->Sumw2();
1824   c1->cd(2);
1825   //            gPad->SetLogy();
1826   //    gPad->SetLogx();
1827   badshHF->SetMarkerStyle(20);
1828   badshHF->SetMarkerSize(0.4);
1829   badshHF->GetYaxis()->SetLabelSize(0.04);
1830   badshHF->SetXTitle("TS weighted by A bad HF \b");
1831   badshHF->SetMarkerColor(2);
1832   badshHF->SetLineColor(2);
1833   badshHF->Draw("");
1834   c1->cd(4);
1835   badshHF0->SetMarkerStyle(20);
1836   badshHF0->SetMarkerSize(0.4);
1837   badshHF0->GetYaxis()->SetLabelSize(0.04);
1838   badshHF0->SetXTitle("TS bad HF \b");
1839   badshHF0->SetMarkerColor(2);
1840   badshHF0->SetLineColor(2);
1841   badshHF0->Draw("");
1842   c1->cd(6);
1843   badshHFr->SetMarkerStyle(20);
1844   badshHFr->SetMarkerSize(0.4);
1845   badshHFr->GetYaxis()->SetLabelSize(0.04);
1846   badshHFr->SetXTitle("shape bad HF per event, per channel\b");
1847   badshHFr->SetMarkerColor(2);
1848   badshHFr->SetLineColor(2);
1849   badshHFr->Draw("");
1850 
1851   c1->Update();
1852   //======================================================================
1853   //==================================================================================================== end
1854   //======================================================================
1855   //======================================================================
1856   // close and delete all possible things:
1857 
1858   //   psfile->Close();
1859   psfile.Close();
1860 
1861   hfile1->Close();
1862   //    hfile1->Delete();
1863   hfile1->Close();
1864   //    hfile1->Delete();
1865 
1866   //  Exit Root
1867   gSystem->Exit(0);
1868   //======================================================================
1869 }