CMS 3D CMS Logo
zadcamplitude.cc
Go to the documentation of this file.
1 #include <iostream>
2 #include <fstream>
3 #include "TH1.h"
4 #include "TH2.h"
5 #include "TCanvas.h"
6 #include "TROOT.h"
7 #include <TMath.h>
8 #include "TStyle.h"
9 #include "TSystem.h"
10 #include "TLegend.h"
11 #include "TText.h"
12 #include "TAxis.h"
13 #include "TFile.h"
14 #include "TLine.h"
15 #include "TGraph.h"
16 #include <TPaveText.h>
17 #include <TPostScript.h>
18 #include <TClass.h>
19 
20 //
21 using namespace std;
22 //
23 //
24 
25 //inline void HERE(const char *msg) { std::cout << msg << std::endl; }
26 
27 int main(int argc, char *argv[]) {
28  //======================================================================
29  printf("reco: gROOT Reset \n");
30  gROOT->Reset();
31  gROOT->SetStyle("Plain");
32  // gStyle->SetOptStat(0); // no statistics _or_
33  // gStyle->SetOptStat(11111111);
34  //gStyle->SetOptStat(1101);// name mean and rms
35  // gStyle->SetOptStat(0101);// name and entries
36  // gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
37  gStyle->SetOptStat(1110000); // und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
38  // gStyle->SetOptStat(101110); // entries, mean, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
39  // gStyle->SetOptStat(100000);// over !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
40  //
41  //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
42  // gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
43  // gStyle->SetOptFit(0010);// constant, mean and sigma only !!
44  // gStyle->SetOptFit(00011);// constant, mean and sigma only !!
45  // gStyle->SetOptFit(1101);
46  // gStyle->SetOptFit(1011);
47  //
48  //gStyle->SetStatX(0.98);
49  //gStyle->SetStatY(0.99);
50  //gStyle->SetStatW(0.30);
51  //gStyle->SetStatH(0.25);
52  gStyle->SetStatX(0.91);
53  gStyle->SetStatY(0.75);
54  gStyle->SetStatW(0.20);
55  gStyle->SetStatH(0.10);
56  // gStyle->SetStatH(0.35);
57  //
58 
59  // Float_t LeftOffset = 0.12;
60  //Float_t TopOffset = 0.12;
61  Float_t LeftOffset = 0.12;
62  Float_t TopOffset = 0.22;
63 
64  gStyle->SetLineWidth(1);
65  gStyle->SetErrorX(0);
66 
67  //---=[ Titles,Labels ]=-----------
68  gStyle->SetOptTitle(0); // title on/off
69  // gStyle->SetTitleColor(0); // title color
70  gStyle->SetTitleColor(1); // title color
71  // gStyle->SetTitleX(0.35); // title x-position
72  gStyle->SetTitleX(0.15); // title x-position
73  gStyle->SetTitleH(0.15); // title height
74  // gStyle->SetTitleW(0.53); // title width
75  gStyle->SetTitleW(0.60); // title width
76  gStyle->SetTitleFont(42); // title font
77  gStyle->SetTitleFontSize(0.07); // title font size
78 
79  gStyle->SetPalette(1);
80  //---=[ Histogram style ]=----------
81  // gStyle->SetHistFillColor(38);
82  // gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
83  // gStyle->SetFrameFillColor(5);// jeltyi
84  // gStyle->SetFrameFillColor(17);// seryi
85  // gStyle->SetFrameFillColor(18);// svetlo seryi
86  // gStyle->SetFrameFillColor(20);// svetlo kori4nev
87  // gStyle->SetFrameFillColor(33);// sine-seryi
88  // gStyle->SetFrameFillColor(40);// fiolet-seryi
89  // gStyle->SetFrameFillColor(23);// sv.kor
90 
91  //---=[ Pad style ]=----------------
92  gStyle->SetPadTopMargin(TopOffset);
93  gStyle->SetPadBottomMargin(LeftOffset);
94  gStyle->SetPadRightMargin(TopOffset);
95  gStyle->SetPadLeftMargin(LeftOffset);
96  //---=[ SetCanvasDef ]=----------------
97  //======================================================================
98  //
99  // Connect the input file and get the 2-d histogram in memory
100  //======================================================================
101  // TBrowser *b = new TBrowser
102 
103  // TFile *hfile1= new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/CMTweb/histos/Global_234556.root", "READ");
104  // TFile *hfile1= new TFile("Global_321177_41.root", "READ");
105  // TFile *hfile1= new TFile("Global_321177_ls1to600.root", "READ");
106  // TFile *hfile1= new TFile("Global_321177_ls1to600.root_no41", "READ");
107  // TFile *hfile1= new TFile("Global_325001_ls1to600.root", "READ");
108  // TFile *hfile1= new TFile("Global_RBX_325001_40.root", "READ");
109  // TFile *hfile1= new TFile("Global_RBX_325001_ls1to600.root", "READ");
111  // TFile *hfile1= new TFile("Global_321177_41_abortgap.root", "READ");
112  // TFile *hfile1= new TFile("Global_321177_ls1to600_abortgap.root", "READ");
113  // TFile *hfile1= new TFile("Global_321177_ls1to600_abortgap_no41.root", "READ");
114  // TFile *hfile1= new TFile("Global_325001_ls1to600_abortgap.root", "READ");
115  // TFile *hfile1= new TFile("Global_321624_1.root", "READ");
116  // TFile *hfile1= new TFile("Global_321625.root", "READ");
117  // TFile *hfile1= new TFile("Global_321313.root", "READ");
118  // TFile *hfile1= new TFile("Global_RBX_325001.root", "READ");
119 
120  // TFile *hfile1= new TFile("Global_RBX_325001test.root", "READ");
121  // TFile *hfile1= new TFile("Global_RBX_321177test.root", "READ");
122  // TFile *hfile1= new TFile("LED_327785test.root", "READ");
123 
124  TFile *hfile1 = new TFile("Global_RBX_325001.root", "READ");
126 
127  // getchar();
128  //
129  TPostScript psfile("zadcamplitude.ps", 111);
130 
131  //
132 
133  TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
134 
135  //=============================================================================================== 1
136  //======================================================================
137  //======================================================================
138  //================
139  //======================================================================
140  c1->Clear();
141  c1->Divide(2, 2);
142 
143  TH2F *twoddepth1hb1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HB");
144  TH2F *twoddepth1hb0 = (TH2F *)hfile1->Get("h_mapDepth1_HB");
145  twoddepth1hb1->Sumw2();
146  twoddepth1hb0->Sumw2();
147  // if(twoddepth1hb0->IsA()->InheritsFrom("TH2F")){
148  TH2F *Cdepth1hbff = (TH2F *)twoddepth1hb1->Clone("Cdepth1hbff");
149  Cdepth1hbff->Divide(twoddepth1hb1, twoddepth1hb0, 1, 1, "B"); // average A
150  Cdepth1hbff->Sumw2();
151  // }
152 
153  c1->cd(1);
154  TH2F *twoedepth1hb1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HB");
155  TH2F *twoedepth1hb0 = (TH2F *)hfile1->Get("h_mapDepth1_HB");
156  twoedepth1hb1->Sumw2();
157  twoedepth1hb0->Sumw2();
158  // if(twoe0->IsA()->InheritsFrom("TH2F")){
159  TH2F *Cdepth1hbfz225 = (TH2F *)twoedepth1hb1->Clone("Cdepth1hbfz225");
160  Cdepth1hbfz225->Divide(twoedepth1hb1, twoedepth1hb0, 1, 1, "B"); // just RATE
161  Cdepth1hbfz225->Sumw2();
162  // }
163  gPad->SetGridy();
164  gPad->SetGridx();
165  gPad->SetLogz();
166  Cdepth1hbfz225->SetMarkerStyle(20);
167  Cdepth1hbfz225->SetMarkerSize(0.4);
168  // Cdepth1hbfz225->GetYaxis()->SetLabelSize(0.04);
169  Cdepth1hbfz225->GetZaxis()->SetLabelSize(0.08);
170  Cdepth1hbfz225->SetXTitle("#eta \b");
171  Cdepth1hbfz225->SetYTitle("#phi \b");
172  Cdepth1hbfz225->SetZTitle("Rate for ADCAmpl in each event & cell - HB Depth1 \b");
173  Cdepth1hbfz225->SetMarkerColor(2);
174  Cdepth1hbfz225->SetLineColor(2);
175  Cdepth1hbfz225->SetMaximum(1.000);
176  Cdepth1hbfz225->SetMinimum(0.00001);
177  // Cdepth1hbfz225->SetMinimum(0.0001);
178  Cdepth1hbfz225->Draw("COLZ");
179 
180  c1->cd(2);
181 
182  // TH1F *adepth1hb= (TH1F*)hfile1->Get("h_ADCAmpl_HB");
183  // TH1F *adepth1hb= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HB");
184  TH1F *adepth1hb = (TH1F *)hfile1->Get("h_ADCAmplZoom_HB");
185  gPad->SetLogy();
186  // gPad->SetLogx();
187  adepth1hb->SetMarkerStyle(20);
188  adepth1hb->SetMarkerSize(0.8);
189  adepth1hb->GetYaxis()->SetLabelSize(0.04);
190  adepth1hb->SetXTitle("ADCAmpl in each event & cell HB \b");
191  adepth1hb->SetMarkerColor(2);
192  adepth1hb->SetLineColor(2);
193  adepth1hb->Draw("");
194 
195  c1->cd(3);
196 
198  TH2F *Diffe_Depth1hb = (TH2F *)Cdepth1hbff->Clone("Diffe_Depth1hb");
199  for (int i = 1; i <= Cdepth1hbff->GetXaxis()->GetNbins(); i++) {
200  for (int j = 1; j <= Cdepth1hbff->GetYaxis()->GetNbins(); j++) {
201  // TH2F* Cdepth1hbff = (TH2F*)twoddepth1hb1->Clone("Cdepth1hbff");
202  double ccc1 = Cdepth1hbff->GetBinContent(i, j);
203  Diffe_Depth1hb->SetBinContent(i, j, 0.);
204  // if(ccc1>0.) cout<<" ibin= "<< i <<" jbin= "<< j <<" nevents= "<< ccc1 <<endl;
205 
206  // if(ccc1 <18|| ccc1>30.) Diffe_Depth1hb->SetBinContent(i,j,ccc1);
207  // if(ccc1 <500.) Diffe_Depth1hb->SetBinContent(i,j,ccc1);
208  Diffe_Depth1hb->SetBinContent(i, j, ccc1);
209  }
210  }
211  gPad->SetGridy();
212  gPad->SetGridx();
213  gPad->SetLogz();
214  Diffe_Depth1hb->SetMarkerStyle(20);
215  Diffe_Depth1hb->SetMarkerSize(0.4);
216  Diffe_Depth1hb->GetZaxis()->SetLabelSize(0.08);
217  // Diffe_Depth1hb->SetTitle("any Error, HB Depth1hb \n");
218  Diffe_Depth1hb->SetXTitle("#eta \b");
219  Diffe_Depth1hb->SetYTitle("#phi \b");
220  Diffe_Depth1hb->SetZTitle("<ADCAmpl> - HB Depth1 \b");
221  Diffe_Depth1hb->SetMarkerColor(2);
222  Diffe_Depth1hb->SetLineColor(2);
223  // Diffe_Depth1hb->SetMaximum(1.000);
224  // Diffe_Depth1hb->SetMinimum(0.0000001);
225  Diffe_Depth1hb->Draw("COLZ");
226 
227  c1->cd(4);
228  // TH1F* diffADCAmpl_Depth1hb = new TH1F("diffADCAmpl_Depth1hb","", 100, 0.,3000.);
229  TH1F *diffADCAmpl_Depth1hb = new TH1F("diffADCAmpl_Depth1hb", "", 250, 0., 1000.);
230  TH2F *Cdepth1hbfw = (TH2F *)Cdepth1hbff->Clone("Cdepth1hbfw");
231  // TH2F* Cdepth1hbfw = (TH2F*)twoddepth1hb1->Clone("diffADCAmpl_Depth1hb");
232  for (int i = 1; i <= Cdepth1hbfw->GetXaxis()->GetNbins(); i++) {
233  for (int j = 1; j <= Cdepth1hbfw->GetYaxis()->GetNbins(); j++) {
234  if (Cdepth1hbfw->GetBinContent(i, j) != 0) {
235  double ccc1 = Cdepth1hbfw->GetBinContent(i, j);
236  diffADCAmpl_Depth1hb->Fill(ccc1);
237  }
238  }
239  }
240  gPad->SetLogy();
241  diffADCAmpl_Depth1hb->SetMarkerStyle(20);
242  diffADCAmpl_Depth1hb->SetMarkerSize(0.4);
243  diffADCAmpl_Depth1hb->GetYaxis()->SetLabelSize(0.04);
244  diffADCAmpl_Depth1hb->SetXTitle("<ADCAmpl> in each cell \b");
245  diffADCAmpl_Depth1hb->SetMarkerColor(2);
246  diffADCAmpl_Depth1hb->SetLineColor(2);
247  diffADCAmpl_Depth1hb->Draw("");
248 
249  c1->Update();
250 
251  delete twoddepth1hb0;
252  delete twoddepth1hb1;
253  delete adepth1hb;
254  delete Cdepth1hbfz225;
255  delete Cdepth1hbff;
256  delete diffADCAmpl_Depth1hb;
257  delete Diffe_Depth1hb;
258 
259  //=============================================================================================== 2
260  //======================================================================
261  //======================================================================
262  //================
263  //======================================================================
264  c1->Clear();
265  c1->Divide(2, 2);
266 
267  TH2F *twoddepth2hb1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HB");
268  TH2F *twoddepth2hb0 = (TH2F *)hfile1->Get("h_mapDepth2_HB");
269  twoddepth2hb1->Sumw2();
270  twoddepth2hb0->Sumw2();
271  // if(twoddepth2hb0->IsA()->InheritsFrom("TH2F")){
272  TH2F *Cdepth2hbff = (TH2F *)twoddepth2hb1->Clone("Cdepth2hbff");
273  Cdepth2hbff->Divide(twoddepth2hb1, twoddepth2hb0, 1, 1, "B");
274  Cdepth2hbff->Sumw2();
275  // }
276  c1->cd(1);
277  TH2F *twoedepth2hb1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HB");
278  TH2F *twoedepth2hb0 = (TH2F *)hfile1->Get("h_mapDepth2_HB");
279  twoedepth2hb1->Sumw2();
280  twoedepth2hb0->Sumw2();
281  // if(twoe0->IsA()->InheritsFrom("TH2F")){
282  TH2F *Cdepth2hbfz225 = (TH2F *)twoedepth2hb1->Clone("Cdepth2hbfz225");
283  Cdepth2hbfz225->Divide(twoedepth2hb1, twoedepth2hb0, 1, 1, "B");
284  Cdepth2hbfz225->Sumw2();
285  // }
286  gPad->SetGridy();
287  gPad->SetGridx();
288  gPad->SetLogz();
289  Cdepth2hbfz225->SetMarkerStyle(20);
290  Cdepth2hbfz225->SetMarkerSize(0.4);
291  // Cdepth2hbfz225->GetYaxis()->SetLabelSize(0.04);
292  Cdepth2hbfz225->GetZaxis()->SetLabelSize(0.08);
293  Cdepth2hbfz225->SetXTitle("#eta \b");
294  Cdepth2hbfz225->SetYTitle("#phi \b");
295  Cdepth2hbfz225->SetZTitle("Rate for ADCAmpl in each event & cell - HB Depth2 \b");
296  Cdepth2hbfz225->SetMarkerColor(2);
297  Cdepth2hbfz225->SetLineColor(2);
298  Cdepth2hbfz225->SetMaximum(1.000);
299  Cdepth2hbfz225->SetMinimum(0.00001);
300  // Cdepth2hbfz225->SetMinimum(0.0001);
301  Cdepth2hbfz225->Draw("COLZ");
302 
303  c1->cd(2);
304 
305  // TH1F *adepth2hb= (TH1F*)hfile1->Get("h_ADCAmpl_HB");
306  TH1F *adepth2hb = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HB");
307  // TH1F *adepth2hb= (TH1F*)hfile1->Get("h_ADCAmplZoom_HB");
308  gPad->SetLogy();
309  // gPad->SetLogx();
310  adepth2hb->SetMarkerStyle(20);
311  adepth2hb->SetMarkerSize(0.8);
312  adepth2hb->GetYaxis()->SetLabelSize(0.04);
313  adepth2hb->SetXTitle("ADCAmpl in each event & cell HB \b");
314  adepth2hb->SetMarkerColor(2);
315  adepth2hb->SetLineColor(2);
316  adepth2hb->Draw("");
317 
318  c1->cd(3);
319 
321  TH2F *Diffe_Depth2hb = (TH2F *)Cdepth2hbff->Clone("Diffe_Depth2hb");
322  for (int i = 1; i <= Cdepth2hbff->GetXaxis()->GetNbins(); i++) {
323  for (int j = 1; j <= Cdepth2hbff->GetYaxis()->GetNbins(); j++) {
324  // TH2F* Cdepth2hbff = (TH2F*)twoddepth2hb1->Clone("Cdepth2hbff");
325  double ccc1 = Cdepth2hbff->GetBinContent(i, j);
326  Diffe_Depth2hb->SetBinContent(i, j, 0.);
327  if (ccc1 < 500.)
328  Diffe_Depth2hb->SetBinContent(i, j, ccc1);
329  }
330  }
331  gPad->SetGridy();
332  gPad->SetGridx();
333  gPad->SetLogz();
334  Diffe_Depth2hb->SetMarkerStyle(20);
335  Diffe_Depth2hb->SetMarkerSize(0.4);
336  Diffe_Depth2hb->GetZaxis()->SetLabelSize(0.08);
337  // Diffe_Depth2hb->SetTitle("any Error, HB Depth2hb \n");
338  Diffe_Depth2hb->SetXTitle("#eta \b");
339  Diffe_Depth2hb->SetYTitle("#phi \b");
340  Diffe_Depth2hb->SetZTitle("<ADCAmpl> smalle 508 - HB Depth2hb \b");
341  Diffe_Depth2hb->SetMarkerColor(2);
342  Diffe_Depth2hb->SetLineColor(2);
343  // Diffe_Depth2hb->SetMaximum(1.000);
344  // Diffe_Depth2hb->SetMinimum(0.0000001);
345  Diffe_Depth2hb->Draw("COLZ");
346 
347  c1->cd(4);
348  // TH1F* diffADCAmpl_Depth2hb = new TH1F("diffADCAmpl_Depth2hb","", 100, 0.,3000.);
349  TH1F *diffADCAmpl_Depth2hb = new TH1F("diffADCAmpl_Depth2hb", "", 250, 0., 1000.);
350  TH2F *Cdepth2hbfw = (TH2F *)Cdepth2hbff->Clone("diffADCAmpl_Depth2hb");
351  // TH2F* Cdepth2hbfw = (TH2F*)twoddepth2hb1->Clone("diffADCAmpl_Depth2hb");
352  for (int i = 1; i <= Cdepth2hbfw->GetXaxis()->GetNbins(); i++) {
353  for (int j = 1; j <= Cdepth2hbfw->GetYaxis()->GetNbins(); j++) {
354  if (Cdepth2hbfw->GetBinContent(i, j) != 0) {
355  double ccc1 = Cdepth2hbfw->GetBinContent(i, j);
356  diffADCAmpl_Depth2hb->Fill(ccc1);
357  }
358  }
359  }
360  gPad->SetLogy();
361  diffADCAmpl_Depth2hb->SetMarkerStyle(20);
362  diffADCAmpl_Depth2hb->SetMarkerSize(0.4);
363  diffADCAmpl_Depth2hb->GetYaxis()->SetLabelSize(0.04);
364  diffADCAmpl_Depth2hb->SetXTitle("<ADCAmpl> in each cell \b");
365  diffADCAmpl_Depth2hb->SetMarkerColor(2);
366  diffADCAmpl_Depth2hb->SetLineColor(2);
367  diffADCAmpl_Depth2hb->Draw("");
368 
369  c1->Update();
370 
371  delete twoddepth2hb0;
372  delete twoddepth2hb1;
373  delete adepth2hb;
374  delete Cdepth2hbfz225;
375  delete Cdepth2hbff;
376  delete diffADCAmpl_Depth2hb;
377  delete Diffe_Depth2hb;
378 
379  //=============================================================================================== 3
380  //======================================================================
381  //======================================================================
382  //================
383  //======================================================================
384  c1->Clear();
385  c1->Divide(2, 2);
386 
387  TH2F *twoddepth1he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HE");
388  TH2F *twoddepth1he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
389  twoddepth1he1->Sumw2();
390  twoddepth1he0->Sumw2();
391  // if(twoddepth1he0->IsA()->InheritsFrom("TH2F")){
392  TH2F *Cdepth1heff = (TH2F *)twoddepth1he1->Clone("Cdepth1heff");
393  Cdepth1heff->Divide(twoddepth1he1, twoddepth1he0, 1, 1, "B");
394  Cdepth1heff->Sumw2();
395  // }
396  c1->cd(1);
397  TH2F *twoedepth1he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HE");
398  TH2F *twoedepth1he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
399  twoedepth1he1->Sumw2();
400  twoedepth1he0->Sumw2();
401  // if(twoe0->IsA()->InheritsFrom("TH2F")){
402  TH2F *Cdepth1hefz225 = (TH2F *)twoedepth1he1->Clone("Cdepth1hefz225");
403  Cdepth1hefz225->Divide(twoedepth1he1, twoedepth1he0, 1, 1, "B");
404  Cdepth1hefz225->Sumw2();
405  // }
406  gPad->SetGridy();
407  gPad->SetGridx();
408  gPad->SetLogz();
409  Cdepth1hefz225->SetMarkerStyle(20);
410  Cdepth1hefz225->SetMarkerSize(0.4);
411  // Cdepth1hefz225->GetYaxis()->SetLabelSize(0.04);
412  Cdepth1hefz225->GetZaxis()->SetLabelSize(0.08);
413  Cdepth1hefz225->SetXTitle("#eta \b");
414  Cdepth1hefz225->SetYTitle("#phi \b");
415  Cdepth1hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth1 \b");
416  Cdepth1hefz225->SetMarkerColor(2);
417  Cdepth1hefz225->SetLineColor(2);
418  Cdepth1hefz225->SetMaximum(1.000);
419  Cdepth1hefz225->SetMinimum(0.00001);
420  // Cdepth1hefz225->SetMinimum(0.0001);
421  Cdepth1hefz225->Draw("COLZ");
422 
423  c1->cd(2);
424 
425  // TH1F *adepth1he= (TH1F*)hfile1->Get("h_ADCAmpl_HE");
426  TH1F *adepth1he = (TH1F *)hfile1->Get("h_ADCAmpl345_HE");
427  // TH1F *adepth1he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
428  gPad->SetLogy();
429  // gPad->SetLogx();
430  adepth1he->SetMarkerStyle(20);
431  adepth1he->SetMarkerSize(0.8);
432  adepth1he->GetYaxis()->SetLabelSize(0.04);
433  adepth1he->SetXTitle("ADCAmpl in each event & cell HE \b");
434  adepth1he->SetMarkerColor(2);
435  adepth1he->SetLineColor(2);
436  adepth1he->Draw("");
437 
438  c1->cd(3);
439 
441  TH2F *Diffe_Depth1he = (TH2F *)Cdepth1heff->Clone("Diffe_Depth1he");
442  for (int i = 1; i <= Cdepth1heff->GetXaxis()->GetNbins(); i++) {
443  for (int j = 1; j <= Cdepth1heff->GetYaxis()->GetNbins(); j++) {
444  // TH2F* Cdepth1heff = (TH2F*)twoddepth1he1->Clone("Cdepth1heff");
445  double ccc1 = Cdepth1heff->GetBinContent(i, j);
446  // if(ccc1>0.) cout<<"HE ibin= "<< i <<" jbin= "<< j <<" nevents= "<< ccc1 <<endl;
447  Diffe_Depth1he->SetBinContent(i, j, 0.);
448  if (ccc1 < 20000000.)
449  Diffe_Depth1he->SetBinContent(i, j, ccc1);
450  }
451  }
452  gPad->SetGridy();
453  gPad->SetGridx();
454  gPad->SetLogz();
455  Diffe_Depth1he->SetMarkerStyle(20);
456  Diffe_Depth1he->SetMarkerSize(0.4);
457  Diffe_Depth1he->GetZaxis()->SetLabelSize(0.08);
458  // Diffe_Depth1he->SetTitle("any Error, HE Depth1he \n");
459  Diffe_Depth1he->SetXTitle("#eta \b");
460  Diffe_Depth1he->SetYTitle("#phi \b");
461  Diffe_Depth1he->SetZTitle("<ADCAmpl> smalle 20000008 - HE Depth1he \b");
462  Diffe_Depth1he->SetMarkerColor(2);
463  Diffe_Depth1he->SetLineColor(2);
464  // Diffe_Depth1he->SetMaximum(1.000);
465  // Diffe_Depth1he->SetMinimum(0.0000001);
466  Diffe_Depth1he->Draw("COLZ");
467 
468  c1->cd(4);
469  // TH1F* diffADCAmpl_Depth1he = new TH1F("diffADCAmpl_Depth1he","", 100, 0.,3000.);
470  TH1F *diffADCAmpl_Depth1he = new TH1F("diffADCAmpl_Depth1he", "", 100, 0., 5000.);
471  TH2F *Cdepth1hefw = (TH2F *)Cdepth1heff->Clone("diffADCAmpl_Depth1he");
472  // TH2F* Cdepth1hefw = (TH2F*)twoddepth1he1->Clone("diffADCAmpl_Depth1he");
473  for (int i = 1; i <= Cdepth1hefw->GetXaxis()->GetNbins(); i++) {
474  for (int j = 1; j <= Cdepth1hefw->GetYaxis()->GetNbins(); j++) {
475  if (Cdepth1hefw->GetBinContent(i, j) != 0) {
476  double ccc1 = Cdepth1hefw->GetBinContent(i, j);
477  diffADCAmpl_Depth1he->Fill(ccc1);
478  }
479  }
480  }
481  gPad->SetLogy();
482  diffADCAmpl_Depth1he->SetMarkerStyle(20);
483  diffADCAmpl_Depth1he->SetMarkerSize(0.4);
484  diffADCAmpl_Depth1he->GetYaxis()->SetLabelSize(0.04);
485  diffADCAmpl_Depth1he->SetXTitle("<ADCAmpl> in each cell \b");
486  diffADCAmpl_Depth1he->SetMarkerColor(2);
487  diffADCAmpl_Depth1he->SetLineColor(2);
488  diffADCAmpl_Depth1he->Draw("");
489 
490  c1->Update();
491 
492  delete twoddepth1he0;
493  delete twoddepth1he1;
494  delete adepth1he;
495  delete Cdepth1hefz225;
496  delete Cdepth1heff;
497  delete diffADCAmpl_Depth1he;
498  delete Diffe_Depth1he;
499 
500  //=============================================================================================== 4
501  //======================================================================
502  //======================================================================
503  //================
504  //======================================================================
505  c1->Clear();
506  c1->Divide(2, 2);
507 
508  TH2F *twoddepth2he1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HE");
509  TH2F *twoddepth2he0 = (TH2F *)hfile1->Get("h_mapDepth2_HE");
510  twoddepth2he1->Sumw2();
511  twoddepth2he0->Sumw2();
512  // if(twoddepth2he0->IsA()->InheritsFrom("TH2F")){
513  TH2F *Cdepth2heff = (TH2F *)twoddepth2he1->Clone("Cdepth2heff");
514  Cdepth2heff->Divide(twoddepth2he1, twoddepth2he0, 1, 1, "B");
515  Cdepth2heff->Sumw2();
516  // }
517  c1->cd(1);
518  TH2F *twoedepth2he1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HE");
519  TH2F *twoedepth2he0 = (TH2F *)hfile1->Get("h_mapDepth2_HE");
520  twoedepth2he1->Sumw2();
521  twoedepth2he0->Sumw2();
522  // if(twoe0->IsA()->InheritsFrom("TH2F")){
523  TH2F *Cdepth2hefz225 = (TH2F *)twoedepth2he1->Clone("Cdepth2hefz225");
524  Cdepth2hefz225->Divide(twoedepth2he1, twoedepth2he0, 1, 1, "B");
525  Cdepth2hefz225->Sumw2();
526  // }
527  gPad->SetGridy();
528  gPad->SetGridx();
529  gPad->SetLogz();
530  Cdepth2hefz225->SetMarkerStyle(20);
531  Cdepth2hefz225->SetMarkerSize(0.4);
532  // Cdepth2hefz225->GetYaxis()->SetLabelSize(0.04);
533  Cdepth2hefz225->GetZaxis()->SetLabelSize(0.08);
534  Cdepth2hefz225->SetXTitle("#eta \b");
535  Cdepth2hefz225->SetYTitle("#phi \b");
536  Cdepth2hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth2 \b");
537  Cdepth2hefz225->SetMarkerColor(2);
538  Cdepth2hefz225->SetLineColor(2);
539  Cdepth2hefz225->SetMaximum(1.000);
540  Cdepth2hefz225->SetMinimum(0.00001);
541  // Cdepth2hefz225->SetMinimum(0.0001);
542  Cdepth2hefz225->Draw("COLZ");
543 
544  c1->cd(2);
545 
546  TH1F *adepth2he = (TH1F *)hfile1->Get("h_ADCAmpl_HE");
547  // TH1F *adepth2he= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HE");
548  // TH1F *adepth2he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
549  gPad->SetLogy();
550  // gPad->SetLogx();
551  adepth2he->SetMarkerStyle(20);
552  adepth2he->SetMarkerSize(0.8);
553  adepth2he->GetYaxis()->SetLabelSize(0.04);
554  adepth2he->SetXTitle("ADCAmpl in each event & cell HE \b");
555  adepth2he->SetMarkerColor(2);
556  adepth2he->SetLineColor(2);
557  adepth2he->Draw("");
558 
559  c1->cd(3);
560 
562  TH2F *Diffe_Depth2he = (TH2F *)Cdepth2heff->Clone("Diffe_Depth2he");
563  for (int i = 1; i <= Cdepth2heff->GetXaxis()->GetNbins(); i++) {
564  for (int j = 1; j <= Cdepth2heff->GetYaxis()->GetNbins(); j++) {
565  // TH2F* Cdepth2heff = (TH2F*)twoddepth2he1->Clone("Cdepth2heff");
566  double ccc1 = Cdepth2heff->GetBinContent(i, j);
567  Diffe_Depth2he->SetBinContent(i, j, 0.);
568  if (ccc1 < 20000000.)
569  Diffe_Depth2he->SetBinContent(i, j, ccc1);
570  }
571  }
572  gPad->SetGridy();
573  gPad->SetGridx();
574  gPad->SetLogz();
575  Diffe_Depth2he->SetMarkerStyle(20);
576  Diffe_Depth2he->SetMarkerSize(0.4);
577  Diffe_Depth2he->GetZaxis()->SetLabelSize(0.08);
578  // Diffe_Depth2he->SetTitle("any Error, HE Depth2he \n");
579  Diffe_Depth2he->SetXTitle("#eta \b");
580  Diffe_Depth2he->SetYTitle("#phi \b");
581  Diffe_Depth2he->SetZTitle("<ADCAmpl> smalle 20000000 - HE Depth2he \b");
582  Diffe_Depth2he->SetMarkerColor(2);
583  Diffe_Depth2he->SetLineColor(2);
584  // Diffe_Depth2he->SetMaximum(1.000);
585  // Diffe_Depth2he->SetMinimum(0.0000001);
586  Diffe_Depth2he->Draw("COLZ");
587 
588  c1->cd(4);
589  // TH1F* diffADCAmpl_Depth2he = new TH1F("diffADCAmpl_Depth2he","", 100, 0.,3000.);
590  TH1F *diffADCAmpl_Depth2he = new TH1F("diffADCAmpl_Depth2he", "", 100, 0., 5000.);
591  TH2F *Cdepth2hefw = (TH2F *)Cdepth2heff->Clone("diffADCAmpl_Depth2he");
592  // TH2F* Cdepth2hefw = (TH2F*)twoddepth2he1->Clone("diffADCAmpl_Depth2he");
593  for (int i = 1; i <= Cdepth2hefw->GetXaxis()->GetNbins(); i++) {
594  for (int j = 1; j <= Cdepth2hefw->GetYaxis()->GetNbins(); j++) {
595  if (Cdepth2hefw->GetBinContent(i, j) != 0) {
596  double ccc1 = Cdepth2hefw->GetBinContent(i, j);
597  diffADCAmpl_Depth2he->Fill(ccc1);
598  }
599  }
600  }
601  gPad->SetLogy();
602  diffADCAmpl_Depth2he->SetMarkerStyle(20);
603  diffADCAmpl_Depth2he->SetMarkerSize(0.4);
604  diffADCAmpl_Depth2he->GetYaxis()->SetLabelSize(0.04);
605  diffADCAmpl_Depth2he->SetXTitle("<ADCAmpl> in each cell \b");
606  diffADCAmpl_Depth2he->SetMarkerColor(2);
607  diffADCAmpl_Depth2he->SetLineColor(2);
608  diffADCAmpl_Depth2he->Draw("");
609 
610  c1->Update();
611 
612  delete twoddepth2he0;
613  delete twoddepth2he1;
614  delete adepth2he;
615  delete Cdepth2hefz225;
616  delete Cdepth2heff;
617  delete diffADCAmpl_Depth2he;
618  delete Diffe_Depth2he;
619 
620  //=============================================================================================== 5
621  //======================================================================
622  //======================================================================
623  //================
624  //======================================================================
625  c1->Clear();
626  c1->Divide(2, 2);
627 
628  TH2F *twoddepth3he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HE");
629  TH2F *twoddepth3he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
630  twoddepth3he1->Sumw2();
631  twoddepth3he0->Sumw2();
632  // if(twoddepth3he0->IsA()->InheritsFrom("TH2F")){
633  TH2F *Cdepth3heff = (TH2F *)twoddepth3he1->Clone("Cdepth3heff");
634  Cdepth3heff->Divide(twoddepth3he1, twoddepth3he0, 1, 1, "B");
635  Cdepth3heff->Sumw2();
636  // }
637  c1->cd(1);
638  TH2F *twoedepth3he1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HE");
639  TH2F *twoedepth3he0 = (TH2F *)hfile1->Get("h_mapDepth1_HE");
640  twoedepth3he1->Sumw2();
641  twoedepth3he0->Sumw2();
642  // if(twoe0->IsA()->InheritsFrom("TH2F")){
643  TH2F *Cdepth3hefz225 = (TH2F *)twoedepth3he1->Clone("Cdepth3hefz225");
644  Cdepth3hefz225->Divide(twoedepth3he1, twoedepth3he0, 1, 1, "B");
645  Cdepth3hefz225->Sumw2();
646  // }
647  gPad->SetGridy();
648  gPad->SetGridx();
649  gPad->SetLogz();
650  Cdepth3hefz225->SetMarkerStyle(20);
651  Cdepth3hefz225->SetMarkerSize(0.4);
652  // Cdepth3hefz225->GetYaxis()->SetLabelSize(0.04);
653  Cdepth3hefz225->GetZaxis()->SetLabelSize(0.08);
654  Cdepth3hefz225->SetXTitle("#eta \b");
655  Cdepth3hefz225->SetYTitle("#phi \b");
656  Cdepth3hefz225->SetZTitle("Rate for ADCAmpl in each event & cell - HE Depth1 \b");
657  Cdepth3hefz225->SetMarkerColor(2);
658  Cdepth3hefz225->SetLineColor(2);
659  Cdepth3hefz225->SetMaximum(1.000);
660  Cdepth3hefz225->SetMinimum(0.00001);
661  // Cdepth3hefz225->SetMinimum(0.0001);
662  Cdepth3hefz225->Draw("COLZ");
663 
664  c1->cd(2);
665 
666  // TH1F *adepth3he= (TH1F*)hfile1->Get("h_ADCAmpl_HE");
667  TH1F *adepth3he = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HE");
668  // TH1F *adepth3he= (TH1F*)hfile1->Get("h_ADCAmplZoom_HE");
669  gPad->SetLogy();
670  // gPad->SetLogx();
671  adepth3he->SetMarkerStyle(20);
672  adepth3he->SetMarkerSize(0.8);
673  adepth3he->GetYaxis()->SetLabelSize(0.04);
674  adepth3he->SetXTitle("ADCAmpl in each event & cell HE \b");
675  adepth3he->SetMarkerColor(2);
676  adepth3he->SetLineColor(2);
677  adepth3he->Draw("");
678 
679  c1->cd(3);
680 
682  TH2F *Diffe_Depth3he = (TH2F *)Cdepth3heff->Clone("Diffe_Depth3he");
683  for (int i = 1; i <= Cdepth3heff->GetXaxis()->GetNbins(); i++) {
684  for (int j = 1; j <= Cdepth3heff->GetYaxis()->GetNbins(); j++) {
685  // TH2F* Cdepth3heff = (TH2F*)twoddepth3he1->Clone("Cdepth3heff");
686  double ccc1 = Cdepth3heff->GetBinContent(i, j);
687  Diffe_Depth3he->SetBinContent(i, j, 0.);
688  if (ccc1 < 20000000.)
689  Diffe_Depth3he->SetBinContent(i, j, ccc1);
690  }
691  }
692  gPad->SetGridy();
693  gPad->SetGridx();
694  gPad->SetLogz();
695  Diffe_Depth3he->SetMarkerStyle(20);
696  Diffe_Depth3he->SetMarkerSize(0.4);
697  Diffe_Depth3he->GetZaxis()->SetLabelSize(0.08);
698  // Diffe_Depth3he->SetTitle("any Error, HE Depth3he \n");
699  Diffe_Depth3he->SetXTitle("#eta \b");
700  Diffe_Depth3he->SetYTitle("#phi \b");
701  Diffe_Depth3he->SetZTitle("<ADCAmpl>smalle 20000000 - HE Depth3he \b");
702  Diffe_Depth3he->SetMarkerColor(2);
703  Diffe_Depth3he->SetLineColor(2);
704  // Diffe_Depth3he->SetMaximum(1.000);
705  // Diffe_Depth3he->SetMinimum(0.0000001);
706  Diffe_Depth3he->Draw("COLZ");
707 
708  c1->cd(4);
709  // TH1F* diffADCAmpl_Depth3he = new TH1F("diffADCAmpl_Depth3he","", 100, 0.,3000.);
710  TH1F *diffADCAmpl_Depth3he = new TH1F("diffADCAmpl_Depth3he", "", 100, 0., 5000.);
711  TH2F *Cdepth3hefw = (TH2F *)Cdepth3heff->Clone("diffADCAmpl_Depth3he");
712  // TH2F* Cdepth3hefw = (TH2F*)twoddepth3he1->Clone("diffADCAmpl_Depth3he");
713  for (int i = 1; i <= Cdepth3hefw->GetXaxis()->GetNbins(); i++) {
714  for (int j = 1; j <= Cdepth3hefw->GetYaxis()->GetNbins(); j++) {
715  if (Cdepth3hefw->GetBinContent(i, j) != 0) {
716  double ccc1 = Cdepth3hefw->GetBinContent(i, j);
717  diffADCAmpl_Depth3he->Fill(ccc1);
718  }
719  }
720  }
721  gPad->SetLogy();
722  diffADCAmpl_Depth3he->SetMarkerStyle(20);
723  diffADCAmpl_Depth3he->SetMarkerSize(0.4);
724  diffADCAmpl_Depth3he->GetYaxis()->SetLabelSize(0.04);
725  diffADCAmpl_Depth3he->SetXTitle("<ADCAmpl> in each cell \b");
726  diffADCAmpl_Depth3he->SetMarkerColor(2);
727  diffADCAmpl_Depth3he->SetLineColor(2);
728  diffADCAmpl_Depth3he->Draw("");
729 
730  c1->Update();
731 
732  delete twoddepth3he0;
733  delete twoddepth3he1;
734  delete adepth3he;
735  delete Cdepth3hefz225;
736  delete Cdepth3heff;
737  delete diffADCAmpl_Depth3he;
738  delete Diffe_Depth3he;
739 
740  //=============================================================================================== 6
741  //======================================================================
742  //======================================================================
743  //================
744  //======================================================================
745  c1->Clear();
746  c1->Divide(2, 2);
747 
748  TH2F *twoddepth1hf1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl_HF");
749  TH2F *twoddepth1hf0 = (TH2F *)hfile1->Get("h_mapDepth1_HF");
750  twoddepth1hf1->Sumw2();
751  twoddepth1hf0->Sumw2();
752  // if(twoddepth1hf0->IsA()->InheritsFrom("TH2F")){
753  TH2F *Cdepth1hfff = (TH2F *)twoddepth1hf1->Clone("Cdepth1hfff");
754  Cdepth1hfff->Divide(twoddepth1hf1, twoddepth1hf0, 1, 1, "B");
755  Cdepth1hfff->Sumw2();
756  // }
757  c1->cd(1);
758  TH2F *twoedepth1hf1 = (TH2F *)hfile1->Get("h_mapDepth1ADCAmpl225_HF");
759  TH2F *twoedepth1hf0 = (TH2F *)hfile1->Get("h_mapDepth1_HF");
760  twoedepth1hf1->Sumw2();
761  twoedepth1hf0->Sumw2();
762  // if(twoe0->IsA()->InheritsFrom("TH2F")){
763  TH2F *Cdepth1hffz225 = (TH2F *)twoedepth1hf1->Clone("Cdepth1hffz225");
764  Cdepth1hffz225->Divide(twoedepth1hf1, twoedepth1hf0, 1, 1, "B");
765  Cdepth1hffz225->Sumw2();
766  // }
767  gPad->SetGridy();
768  gPad->SetGridx();
769  gPad->SetLogz();
770  Cdepth1hffz225->SetMarkerStyle(20);
771  Cdepth1hffz225->SetMarkerSize(0.4);
772  // Cdepth1hffz225->GetYaxis()->SetLabelSize(0.04);
773  Cdepth1hffz225->GetZaxis()->SetLabelSize(0.08);
774  Cdepth1hffz225->SetXTitle("#eta \b");
775  Cdepth1hffz225->SetYTitle("#phi \b");
776  Cdepth1hffz225->SetZTitle("Rate for ADCAmpl in each event & cell - HF Depth1 \b");
777  Cdepth1hffz225->SetMarkerColor(2);
778  Cdepth1hffz225->SetLineColor(2);
779  Cdepth1hffz225->SetMaximum(1.000);
780  Cdepth1hffz225->SetMinimum(0.00001);
781  // Cdepth1hffz225->SetMinimum(0.0001);
782  Cdepth1hffz225->Draw("COLZ");
783 
784  c1->cd(2);
785 
786  TH1F *adepth1hf = (TH1F *)hfile1->Get("h_ADCAmpl_HF");
787  // TH1F *adepth1hf= (TH1F*)hfile1->Get("h_ADCAmplZoom1_HF");
788  // TH1F *adepth1hf= (TH1F*)hfile1->Get("h_ADCAmplZoom_HF");
789  gPad->SetLogy();
790  // gPad->SetLogx();
791  adepth1hf->SetMarkerStyle(20);
792  adepth1hf->SetMarkerSize(0.8);
793  adepth1hf->GetYaxis()->SetLabelSize(0.04);
794  adepth1hf->SetXTitle("ADCAmpl in each event & cell HF \b");
795  adepth1hf->SetMarkerColor(2);
796  adepth1hf->SetLineColor(2);
797  adepth1hf->Draw("");
798 
799  c1->cd(3);
800 
802  TH2F *Diffe_Depth1hf = (TH2F *)Cdepth1hfff->Clone("Diffe_Depth1hf");
803  for (int i = 1; i <= Cdepth1hfff->GetXaxis()->GetNbins(); i++) {
804  for (int j = 1; j <= Cdepth1hfff->GetYaxis()->GetNbins(); j++) {
805  // TH2F* Cdepth1hfff = (TH2F*)twoddepth1hf1->Clone("Cdepth1hfff");
806  double ccc1 = Cdepth1hfff->GetBinContent(i, j);
807  Diffe_Depth1hf->SetBinContent(i, j, 0.);
808  if (ccc1 < 20000000.)
809  Diffe_Depth1hf->SetBinContent(i, j, ccc1);
810  }
811  }
812  gPad->SetGridy();
813  gPad->SetGridx();
814  gPad->SetLogz();
815  Diffe_Depth1hf->SetMarkerStyle(20);
816  Diffe_Depth1hf->SetMarkerSize(0.4);
817  Diffe_Depth1hf->GetZaxis()->SetLabelSize(0.08);
818  // Diffe_Depth1hf->SetTitle("any Error, HF Depth1hf \n");
819  Diffe_Depth1hf->SetXTitle("#eta \b");
820  Diffe_Depth1hf->SetYTitle("#phi \b");
821  Diffe_Depth1hf->SetZTitle("<ADCAmpl> smalle 20000000 - HF Depth1hf \b");
822  Diffe_Depth1hf->SetMarkerColor(2);
823  Diffe_Depth1hf->SetLineColor(2);
824  // Diffe_Depth1hf->SetMaximum(1.000);
825  // Diffe_Depth1hf->SetMinimum(0.0000001);
826  Diffe_Depth1hf->Draw("COLZ");
827 
828  c1->cd(4);
829  TH1F *diffADCAmpl_Depth1hf = new TH1F("diffADCAmpl_Depth1hf", "", 100, 0., 300.);
830  // TH1F* diffADCAmpl_Depth1hf = new TH1F("diffADCAmpl_Depth1hf","", 1000, 0.,1000.);
831  TH2F *Cdepth1hffw = (TH2F *)Cdepth1hfff->Clone("diffADCAmpl_Depth1hf");
832  // TH2F* Cdepth1hffw = (TH2F*)twoddepth1hf1->Clone("diffADCAmpl_Depth1hf");
833  for (int i = 1; i <= Cdepth1hffw->GetXaxis()->GetNbins(); i++) {
834  for (int j = 1; j <= Cdepth1hffw->GetYaxis()->GetNbins(); j++) {
835  if (Cdepth1hffw->GetBinContent(i, j) != 0) {
836  double ccc1 = Cdepth1hffw->GetBinContent(i, j);
837  diffADCAmpl_Depth1hf->Fill(ccc1);
838  }
839  }
840  }
841  gPad->SetLogy();
842  diffADCAmpl_Depth1hf->SetMarkerStyle(20);
843  diffADCAmpl_Depth1hf->SetMarkerSize(0.4);
844  diffADCAmpl_Depth1hf->GetYaxis()->SetLabelSize(0.04);
845  diffADCAmpl_Depth1hf->SetXTitle("<ADCAmpl> in each cell \b");
846  diffADCAmpl_Depth1hf->SetMarkerColor(2);
847  diffADCAmpl_Depth1hf->SetLineColor(2);
848  diffADCAmpl_Depth1hf->Draw("");
849 
850  c1->Update();
851 
852  delete twoddepth1hf0;
853  delete twoddepth1hf1;
854  delete adepth1hf;
855  delete Cdepth1hffz225;
856  delete Cdepth1hfff;
857  delete diffADCAmpl_Depth1hf;
858  delete Diffe_Depth1hf;
859 
860  //=============================================================================================== 7
861  //======================================================================
862  //======================================================================
863  //================
864  //======================================================================
865  c1->Clear();
866  c1->Divide(2, 2);
867 
868  TH2F *twoddepth2hf1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl_HF");
869  TH2F *twoddepth2hf0 = (TH2F *)hfile1->Get("h_mapDepth2_HF");
870  twoddepth2hf1->Sumw2();
871  twoddepth2hf0->Sumw2();
872  // if(twoddepth2hf0->IsA()->InheritsFrom("TH2F")){
873  TH2F *Cdepth2hfff = (TH2F *)twoddepth2hf1->Clone("Cdepth2hfff");
874  Cdepth2hfff->Divide(twoddepth2hf1, twoddepth2hf0, 1, 1, "B");
875  Cdepth2hfff->Sumw2();
876  // }
877  c1->cd(1);
878  TH2F *twoedepth2hf1 = (TH2F *)hfile1->Get("h_mapDepth2ADCAmpl225_HF");
879  TH2F *twoedepth2hf0 = (TH2F *)hfile1->Get("h_mapDepth2_HF");
880  twoedepth2hf1->Sumw2();
881  twoedepth2hf0->Sumw2();
882  // if(twoe0->IsA()->InheritsFrom("TH2F")){
883  TH2F *Cdepth2hffz225 = (TH2F *)twoedepth2hf1->Clone("Cdepth2hffz225");
884  Cdepth2hffz225->Divide(twoedepth2hf1, twoedepth2hf0, 1, 1, "B");
885  Cdepth2hffz225->Sumw2();
886  // }
887  gPad->SetGridy();
888  gPad->SetGridx();
889  gPad->SetLogz();
890  Cdepth2hffz225->SetMarkerStyle(20);
891  Cdepth2hffz225->SetMarkerSize(0.4);
892  // Cdepth2hffz225->GetYaxis()->SetLabelSize(0.04);
893  Cdepth2hffz225->GetZaxis()->SetLabelSize(0.08);
894  Cdepth2hffz225->SetXTitle("#eta \b");
895  Cdepth2hffz225->SetYTitle("#phi \b");
896  Cdepth2hffz225->SetZTitle("Rate for ADCAmpl in each event & cell - HF Depth2 \b");
897  Cdepth2hffz225->SetMarkerColor(2);
898  Cdepth2hffz225->SetLineColor(2);
899  Cdepth2hffz225->SetMaximum(1.000);
900  Cdepth2hffz225->SetMinimum(0.00001);
901  // Cdepth2hffz225->SetMinimum(0.0001);
902  Cdepth2hffz225->Draw("COLZ");
903 
904  c1->cd(2);
905 
906  // TH1F *adepth2hf= (TH1F*)hfile1->Get("h_ADCAmpl_HF");
907  TH1F *adepth2hf = (TH1F *)hfile1->Get("h_ADCAmplZoom1_HF");
908  // TH1F *adepth2hf= (TH1F*)hfile1->Get("h_ADCAmplZoom_HF");
909  gPad->SetLogy();
910  // gPad->SetLogx();
911  adepth2hf->SetMarkerStyle(20);
912  adepth2hf->SetMarkerSize(0.8);
913  adepth2hf->GetYaxis()->SetLabelSize(0.04);
914  adepth2hf->SetXTitle("ADCAmpl in each event & cell HF \b");
915  adepth2hf->SetMarkerColor(2);
916  adepth2hf->SetLineColor(2);
917  adepth2hf->Draw("");
918 
919  c1->cd(3);
920 
922  TH2F *Diffe_Depth2hf = (TH2F *)Cdepth2hfff->Clone("Diffe_Depth2hf");
923  for (int i = 1; i <= Cdepth2hfff->GetXaxis()->GetNbins(); i++) {
924  for (int j = 1; j <= Cdepth2hfff->GetYaxis()->GetNbins(); j++) {
925  // TH2F* Cdepth2hfff = (TH2F*)twoddepth2hf1->Clone("Cdepth2hfff");
926  double ccc1 = Cdepth2hfff->GetBinContent(i, j);
927  Diffe_Depth2hf->SetBinContent(i, j, 0.);
928  if (ccc1 < 20000000.)
929  Diffe_Depth2hf->SetBinContent(i, j, ccc1);
930  }
931  }
932  gPad->SetGridy();
933  gPad->SetGridx();
934  gPad->SetLogz();
935  Diffe_Depth2hf->SetMarkerStyle(20);
936  Diffe_Depth2hf->SetMarkerSize(0.4);
937  Diffe_Depth2hf->GetZaxis()->SetLabelSize(0.08);
938  // Diffe_Depth2hf->SetTitle("any Error, HF Depth2hf \n");
939  Diffe_Depth2hf->SetXTitle("#eta \b");
940  Diffe_Depth2hf->SetYTitle("#phi \b");
941  Diffe_Depth2hf->SetZTitle("<ADCAmpl> smalle 20000000 - HF Depth2hf \b");
942  Diffe_Depth2hf->SetMarkerColor(2);
943  Diffe_Depth2hf->SetLineColor(2);
944  // Diffe_Depth2hf->SetMaximum(1.000);
945  // Diffe_Depth2hf->SetMinimum(0.0000001);
946  Diffe_Depth2hf->Draw("COLZ");
947 
948  c1->cd(4);
949  TH1F *diffADCAmpl_Depth2hf = new TH1F("diffADCAmpl_Depth2hf", "", 100, 0., 300.);
950  // TH1F* diffADCAmpl_Depth2hf = new TH1F("diffADCAmpl_Depth2hf","", 1000, 0.,1000.);
951  TH2F *Cdepth2hffw = (TH2F *)Cdepth2hfff->Clone("diffADCAmpl_Depth2hf");
952  // TH2F* Cdepth2hffw = (TH2F*)twoddepth2hf1->Clone("diffADCAmpl_Depth2hf");
953  for (int i = 1; i <= Cdepth2hffw->GetXaxis()->GetNbins(); i++) {
954  for (int j = 1; j <= Cdepth2hffw->GetYaxis()->GetNbins(); j++) {
955  if (Cdepth2hffw->GetBinContent(i, j) != 0) {
956  double ccc1 = Cdepth2hffw->GetBinContent(i, j);
957  diffADCAmpl_Depth2hf->Fill(ccc1);
958  }
959  }
960  }
961  gPad->SetLogy();
962  diffADCAmpl_Depth2hf->SetMarkerStyle(20);
963  diffADCAmpl_Depth2hf->SetMarkerSize(0.4);
964  diffADCAmpl_Depth2hf->GetYaxis()->SetLabelSize(0.04);
965  diffADCAmpl_Depth2hf->SetXTitle("<ADCAmpl> in each cell \b");
966  diffADCAmpl_Depth2hf->SetMarkerColor(2);
967  diffADCAmpl_Depth2hf->SetLineColor(2);
968  diffADCAmpl_Depth2hf->Draw("");
969 
970  c1->Update();
971 
972  delete twoddepth2hf0;
973  delete twoddepth2hf1;
974  delete adepth2hf;
975  delete Cdepth2hffz225;
976  delete Cdepth2hfff;
977  delete diffADCAmpl_Depth2hf;
978  delete Diffe_Depth2hf;
979 
980  //=============================================================================================== 8
981  //======================================================================
982  //======================================================================
983  //================
984  //======================================================================
985  c1->Clear();
986  c1->Divide(2, 2);
987 
988  TH2F *twoddepth4ho1 = (TH2F *)hfile1->Get("h_mapDepth4ADCAmpl_HO");
989  TH2F *twoddepth4ho0 = (TH2F *)hfile1->Get("h_mapDepth4_HO");
990  twoddepth4ho1->Sumw2();
991  twoddepth4ho0->Sumw2();
992  // if(twoddepth4ho0->IsA()->InheritsFrom("TH2F")){
993  TH2F *Cdepth4hoff = (TH2F *)twoddepth4ho1->Clone("Cdepth4hoff");
994  Cdepth4hoff->Divide(twoddepth4ho1, twoddepth4ho0, 1, 1, "B");
995  Cdepth4hoff->Sumw2();
996  // }
997  c1->cd(1);
998  TH2F *twoedepth4ho1 = (TH2F *)hfile1->Get("h_mapDepth4ADCAmpl225_HO");
999  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 
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:
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:
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:
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:
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 }
mps_fire.i
i
Definition: mps_fire.py:429
cond::time::time0
const boost::posix_time::ptime time0
Definition: TimeConversions.h:37
l1trig_cff.shape
shape
Definition: l1trig_cff.py:152
main
int main(int argc, char *argv[])
Definition: zadcamplitude.cc:27
hcalRecHitTable_cff.time
time
Definition: hcalRecHitTable_cff.py:14
Generated for CMSSW Reference Manual by doxygen 1.8.14