CMS 3D CMS Logo

ztsmean.cc
Go to the documentation of this file.
1 #include <string>
2 #include <iomanip>
3 //#####
4 #include <iostream>
5 #include <fstream>
6 #include "TH1.h"
7 #include "TH2.h"
8 #include "TCanvas.h"
9 #include "TROOT.h"
10 #include <TMath.h>
11 #include "TStyle.h"
12 #include "TSystem.h"
13 #include "TLegend.h"
14 #include "TText.h"
15 #include "TAxis.h"
16 #include "TFile.h"
17 #include "TLine.h"
18 #include "TGraph.h"
19 #include <TPaveText.h>
20 #include <TPostScript.h>
21 #include <TClass.h>
22 
23 //
24 using namespace std;
25 //
26 //
27 
28 //inline void HERE(const char *msg) { std::cout << msg << std::endl; }
29 
30 int main(int argc, char *argv[]) {
31  std::string dirnm = "Analyzer";
32  //======================================================================
33  printf("reco: gROOT Reset \n");
34  gROOT->Reset();
35  gROOT->SetStyle("Plain");
36  gStyle->SetOptStat(0); // no statistics _or_
37  // gStyle->SetOptStat(11111111);
38  //gStyle->SetOptStat(1101);// mame mean and rms
39  // gStyle->SetOptStat(0101);// name and entries
40  // gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
41  // gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
42  // gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
43  // gStyle->SetOptStat(101110); // entries, meam, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
44  //
45  //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
46  // gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
47  gStyle->SetOptFit(0010); // constant, mean and sigma only !!
48  // gStyle->SetOptFit(00011);// constant, mean and sigma only !!
49  // gStyle->SetOptFit(1101);
50  // gStyle->SetOptFit(1011);
51  //
52  //gStyle->SetStatX(0.98);
53  //gStyle->SetStatY(0.99);
54  //gStyle->SetStatW(0.30);
55  //gStyle->SetStatH(0.25);
56  gStyle->SetStatX(0.91);
57  gStyle->SetStatY(0.75);
58  gStyle->SetStatW(0.20);
59  gStyle->SetStatH(0.10);
60  // gStyle->SetStatH(0.35);
61  //
62 
63  // Float_t LeftOffset = 0.12;
64  //Float_t TopOffset = 0.12;
65  Float_t LeftOffset = 0.12;
66  Float_t TopOffset = 0.22;
67 
68  gStyle->SetLineWidth(1);
69  gStyle->SetErrorX(0);
70 
71  //---=[ Titles,Labels ]=-----------
72  gStyle->SetOptTitle(0); // title on/off
73  // gStyle->SetTitleColor(0); // title color
74  gStyle->SetTitleColor(1); // title color
75  // gStyle->SetTitleX(0.35); // title x-position
76  gStyle->SetTitleX(0.15); // title x-position
77  gStyle->SetTitleH(0.15); // title height
78  // gStyle->SetTitleW(0.53); // title width
79  gStyle->SetTitleW(0.60); // title width
80  gStyle->SetTitleFont(42); // title font
81  gStyle->SetTitleFontSize(0.07); // title font size
82 
83  gStyle->SetPalette(1);
84  //---=[ Histogram style ]=----------
85  // gStyle->SetHistFillColor(38);
86  // gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
87  // gStyle->SetFrameFillColor(5);// jeltyi
88  // gStyle->SetFrameFillColor(17);// seryi
89  // gStyle->SetFrameFillColor(18);// svetlo seryi
90  // gStyle->SetFrameFillColor(20);// svetlo kori4nev
91  // gStyle->SetFrameFillColor(33);// sine-seryi
92  // gStyle->SetFrameFillColor(40);// fiolet-seryi
93  // gStyle->SetFrameFillColor(23);// sv.kor
94 
95  //---=[ Pad style ]=----------------
96  gStyle->SetPadTopMargin(TopOffset);
97  gStyle->SetPadBottomMargin(LeftOffset);
98  gStyle->SetPadRightMargin(TopOffset);
99  gStyle->SetPadLeftMargin(LeftOffset);
100  //---=[ SetCanvasDef ]=----------------
101  //======================================================================
102  //
103  // Connect the input file and get the 2-d histogram in memory
104  //======================================================================
105  // TBrowser *b = new TBrowser
106 
107  // TFile *hfile1= new TFile("/afs/cern.ch/user/z/zhokin/public/html/RMT/lastMarch2014/cmt/Global_229713.root", "READ");
108  // TFile *hfile1= new TFile("test.root", "READ");
109  //
110 
111  TFile *hfile1 = new TFile("Global_362596.root", "READ");
112  //TFile *hfile1 = new TFile("Global_362365.root", "READ");
113 
114  TPostScript psfile("ztsmean.ps", 111);
115 
116  //
117 
118  TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
119 
120  hfile1->ls();
121  TDirectory *dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
122 
123  //=============================================================================================== 1
124  //======================================================================
125  //======================================================================
126  //================
127  //======================================================================
128  c1->Clear();
129  c1->Divide(2, 2);
130 
131  TH2F *twod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HB");
132  TH2F *twod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
133  twod1->Sumw2();
134  twod0->Sumw2();
135  // if(twod0->IsA()->InheritsFrom("TH2F")){
136  TH2F *Ceff = (TH2F *)twod1->Clone("Ceff");
137  Ceff->Divide(twod1, twod0, 1, 1, "B");
138  Ceff->Sumw2();
139  // }
140  c1->cd(1);
141  TH2F *twop1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HB");
142  TH2F *twop0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
143  twop1->Sumw2();
144  twop0->Sumw2();
145  // if(twop0->IsA()->InheritsFrom("TH2F")){
146  TH2F *Cefz225 = (TH2F *)twop1->Clone("Cefz225");
147  Cefz225->Divide(twop1, twop0, 1, 1, "B");
148  Cefz225->Sumw2();
149  // }
150  gPad->SetGridy();
151  gPad->SetGridx();
152  gPad->SetLogz();
153  Cefz225->SetMarkerStyle(20);
154  Cefz225->SetMarkerSize(0.4);
155  // Cefz225->GetYaxis()->SetLabelSize(0.04);
156  Cefz225->GetZaxis()->SetLabelSize(0.08);
157  Cefz225->SetXTitle("#eta \b");
158  Cefz225->SetYTitle("#phi \b");
159  Cefz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth1 \b");
160  Cefz225->SetMarkerColor(2);
161  Cefz225->SetLineColor(2);
162  Cefz225->SetMaximum(1.000);
163  Cefz225->SetMinimum(0.0001);
164  Cefz225->Draw("COLZ");
165 
166  c1->cd(2);
167  TH1F *aaaaaa1 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HB");
168  gPad->SetLogy();
169  aaaaaa1->SetMarkerStyle(20);
170  aaaaaa1->SetMarkerSize(0.8);
171  aaaaaa1->GetYaxis()->SetLabelSize(0.04);
172  aaaaaa1->SetXTitle("TSmeanA in each event & cell HB \b");
173  aaaaaa1->SetMarkerColor(2);
174  aaaaaa1->SetLineColor(2);
175  aaaaaa1->Draw("");
176 
177  c1->cd(3);
179  TH2F *Diffe_Depth1_HB = (TH2F *)Ceff->Clone("Diffe_Depth1_HB");
180  for (int i = 1; i <= Ceff->GetXaxis()->GetNbins(); i++) {
181  for (int j = 1; j <= Ceff->GetYaxis()->GetNbins(); j++) {
182  double ccc1 = Ceff->GetBinContent(i, j);
183  Diffe_Depth1_HB->SetBinContent(i, j, 0.);
184  // Diffe_Depth1_HB->SetBinContent(i,j,ccc1);
185  if (ccc1 < 2.5 || ccc1 > 4.5)
186  Diffe_Depth1_HB->SetBinContent(i, j, ccc1);
187  }
188  }
189  gPad->SetGridy();
190  gPad->SetGridx();
191  gPad->SetLogz();
192  Diffe_Depth1_HB->SetMarkerStyle(20);
193  Diffe_Depth1_HB->SetMarkerSize(0.4);
194  Diffe_Depth1_HB->GetZaxis()->SetLabelSize(0.08);
195  // Diffe_Depth1_HB->SetTitle("any Error, HB Depth1 \n");
196  Diffe_Depth1_HB->SetXTitle("#eta \b");
197  Diffe_Depth1_HB->SetYTitle("#phi \b");
198  Diffe_Depth1_HB->SetZTitle("<TSmeanA> out 2.5-4.5- HB Depth1 \b");
199  Diffe_Depth1_HB->SetMarkerColor(2);
200  Diffe_Depth1_HB->SetLineColor(2);
201  Diffe_Depth1_HB->Draw("COLZ");
202 
203  c1->cd(4);
204  TH1F *diffTSmeanA_Depth1_HB = new TH1F("diffTSmeanA_Depth1_HB", "", 100, 1.0, 6.0);
205  for (int i = 1; i <= Ceff->GetXaxis()->GetNbins(); i++) {
206  for (int j = 1; j <= Ceff->GetYaxis()->GetNbins(); j++) {
207  if (Ceff->GetBinContent(i, j) != 0) {
208  double ccc1 = Ceff->GetBinContent(i, j);
209  diffTSmeanA_Depth1_HB->Fill(ccc1);
210  }
211  }
212  }
213  gPad->SetLogy();
214  diffTSmeanA_Depth1_HB->SetMarkerStyle(20);
215  diffTSmeanA_Depth1_HB->SetMarkerSize(0.4);
216  diffTSmeanA_Depth1_HB->GetYaxis()->SetLabelSize(0.04);
217  diffTSmeanA_Depth1_HB->SetXTitle("<TSmeanA> in each cell \b");
218  diffTSmeanA_Depth1_HB->SetMarkerColor(2);
219  diffTSmeanA_Depth1_HB->SetLineColor(2);
220  diffTSmeanA_Depth1_HB->Draw("");
221 
222  c1->Update();
223 
224  //=============================================================================================== 2
225  //======================================================================
226  //======================================================================
227  //================
228  //======================================================================
229  c1->Clear();
230  c1->Divide(2, 2);
231 
232  TH2F *awod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HB");
233  TH2F *awod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
234  awod1->Sumw2();
235  awod0->Sumw2();
236  // if(awod0->IsA()->InheritsFrom("TH2F")){
237  TH2F *C2ff = (TH2F *)awod1->Clone("C2ff");
238  C2ff->Divide(awod1, awod0, 1, 1, "B");
239  C2ff->Sumw2();
240  // }
241  c1->cd(1);
242  TH2F *bwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HB");
243  TH2F *bwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
244  bwod1->Sumw2();
245  bwod0->Sumw2();
246  // if(bwod0->IsA()->InheritsFrom("TH2F")){
247  TH2F *C2fz225 = (TH2F *)bwod1->Clone("C2fz225");
248  C2fz225->Divide(bwod1, bwod0, 1, 1, "B");
249  C2fz225->Sumw2();
250  // }
251  gPad->SetGridy();
252  gPad->SetGridx();
253  gPad->SetLogz();
254  C2fz225->SetMarkerStyle(20);
255  C2fz225->SetMarkerSize(0.4);
256  C2fz225->GetZaxis()->SetLabelSize(0.08);
257  C2fz225->SetXTitle("#eta \b");
258  C2fz225->SetYTitle("#phi \b");
259  C2fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HB Depth2 \b");
260  C2fz225->SetMarkerColor(2);
261  C2fz225->SetLineColor(2);
262  C2fz225->SetMaximum(1.000);
263  C2fz225->SetMinimum(0.0001);
264  C2fz225->Draw("COLZ");
265 
266  c1->cd(2);
267  TH1F *aaaaaa2 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HB");
268  gPad->SetLogy();
269  aaaaaa2->SetMarkerStyle(20);
270  aaaaaa2->SetMarkerSize(0.8);
271  aaaaaa2->GetYaxis()->SetLabelSize(0.04);
272  aaaaaa2->SetXTitle("TSmeanA in each event & cell HB \b");
273  aaaaaa2->SetMarkerColor(2);
274  aaaaaa2->SetLineColor(2);
275  aaaaaa2->Draw("");
276 
277  c1->cd(3);
279  TH2F *Diffe_Depth2_HB = (TH2F *)C2ff->Clone("Diffe_Depth2_HB");
280  for (int i = 1; i <= C2ff->GetXaxis()->GetNbins(); i++) {
281  for (int j = 1; j <= C2ff->GetYaxis()->GetNbins(); j++) {
282  double ccc1 = C2ff->GetBinContent(i, j);
283  Diffe_Depth2_HB->SetBinContent(i, j, 0.);
284  if (ccc1 < 2.4 || ccc1 > 4.4)
285  Diffe_Depth2_HB->SetBinContent(i, j, ccc1);
286  }
287  }
288  gPad->SetGridy();
289  gPad->SetGridx();
290  gPad->SetLogz();
291  Diffe_Depth2_HB->SetMarkerStyle(20);
292  Diffe_Depth2_HB->SetMarkerSize(0.4);
293  Diffe_Depth2_HB->GetZaxis()->SetLabelSize(0.08);
294  // Diffe_Depth2_HB->SetTitle("any Error, HB Depth2 \n");
295  Diffe_Depth2_HB->SetXTitle("#eta \b");
296  Diffe_Depth2_HB->SetYTitle("#phi \b");
297  Diffe_Depth2_HB->SetZTitle("<TSmeanA> out 2.4-4.4 - HB Depth2 \b");
298  Diffe_Depth2_HB->SetMarkerColor(2);
299  Diffe_Depth2_HB->SetLineColor(2);
300  Diffe_Depth2_HB->Draw("COLZ");
301 
302  c1->cd(4);
303  TH1F *diffTSmeanA_Depth2_HB = new TH1F("diffTSmeanA_Depth2_HB", "", 100, 1.0, 6.0);
304  for (int i = 1; i <= C2ff->GetXaxis()->GetNbins(); i++) {
305  for (int j = 1; j <= C2ff->GetYaxis()->GetNbins(); j++) {
306  if (C2ff->GetBinContent(i, j) != 0) {
307  double ccc1 = C2ff->GetBinContent(i, j);
308  diffTSmeanA_Depth2_HB->Fill(ccc1);
309  }
310  }
311  }
312  gPad->SetLogy();
313  diffTSmeanA_Depth2_HB->SetMarkerStyle(20);
314  diffTSmeanA_Depth2_HB->SetMarkerSize(0.4);
315  diffTSmeanA_Depth2_HB->GetYaxis()->SetLabelSize(0.04);
316  diffTSmeanA_Depth2_HB->SetXTitle("<TSmeanA> in each cell \b");
317  diffTSmeanA_Depth2_HB->SetMarkerColor(2);
318  diffTSmeanA_Depth2_HB->SetLineColor(2);
319  diffTSmeanA_Depth2_HB->Draw("");
320 
321  c1->Update();
322 
323  //=============================================================================================== 3
324  //======================================================================
325  //======================================================================
326  //================
327  //======================================================================
328  c1->Clear();
329  c1->Divide(2, 2);
330 
331  TH2F *cwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HE");
332  TH2F *cwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
333  cwod1->Sumw2();
334  cwod0->Sumw2();
335  // if(cwod0->IsA()->InheritsFrom("TH2F")){
336  TH2F *C3ff = (TH2F *)cwod1->Clone("C3ff");
337  C3ff->Divide(cwod1, cwod0, 1, 1, "B");
338  C3ff->Sumw2();
339  // }
340  c1->cd(1);
341  TH2F *dwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HE");
342  TH2F *dwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
343  dwod1->Sumw2();
344  dwod0->Sumw2();
345  // if(dwod0->IsA()->InheritsFrom("TH2F")){
346  TH2F *C3fz225 = (TH2F *)dwod1->Clone("C3fz225");
347  C3fz225->Divide(dwod1, dwod0, 1, 1, "B");
348  C3fz225->Sumw2();
349  // }
350  gPad->SetGridy();
351  gPad->SetGridx();
352  gPad->SetLogz();
353  C3fz225->SetMarkerStyle(20);
354  C3fz225->SetMarkerSize(0.4);
355  C3fz225->GetZaxis()->SetLabelSize(0.08);
356  C3fz225->SetXTitle("#eta \b");
357  C3fz225->SetYTitle("#phi \b");
358  C3fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth1 \b");
359  C3fz225->SetMarkerColor(2);
360  C3fz225->SetLineColor(2);
361  C3fz225->SetMaximum(1.000);
362  C3fz225->SetMinimum(0.0001);
363  C3fz225->Draw("COLZ");
364 
365  c1->cd(2);
366  TH1F *aaaaaa3 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
367  gPad->SetLogy();
368  aaaaaa3->SetMarkerStyle(20);
369  aaaaaa3->SetMarkerSize(0.8);
370  aaaaaa3->GetYaxis()->SetLabelSize(0.04);
371  aaaaaa3->SetXTitle("TSmeanA in each event & cell HE \b");
372  aaaaaa3->SetMarkerColor(2);
373  aaaaaa3->SetLineColor(2);
374  aaaaaa3->Draw("");
375 
376  c1->cd(3);
378  TH2F *Diffe_Depth1_HE = (TH2F *)C3ff->Clone("Diffe_Depth1_HE");
379  for (int i = 1; i <= C3ff->GetXaxis()->GetNbins(); i++) {
380  for (int j = 1; j <= C3ff->GetYaxis()->GetNbins(); j++) {
381  double ccc1 = C3ff->GetBinContent(i, j);
382  Diffe_Depth1_HE->SetBinContent(i, j, 0.);
383  if (ccc1 < 2.4 || ccc1 > 4.6)
384  Diffe_Depth1_HE->SetBinContent(i, j, ccc1);
385  }
386  }
387  gPad->SetGridy();
388  gPad->SetGridx();
389  gPad->SetLogz();
390  Diffe_Depth1_HE->SetMarkerStyle(20);
391  Diffe_Depth1_HE->SetMarkerSize(0.4);
392  Diffe_Depth1_HE->GetZaxis()->SetLabelSize(0.08);
393  // Diffe_Depth1_HE->SetTitle("any Error, HE Depth1 \n");
394  Diffe_Depth1_HE->SetXTitle("#eta \b");
395  Diffe_Depth1_HE->SetYTitle("#phi \b");
396  Diffe_Depth1_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth1 \b");
397  Diffe_Depth1_HE->SetMarkerColor(2);
398  Diffe_Depth1_HE->SetLineColor(2);
399  Diffe_Depth1_HE->Draw("COLZ");
400 
401  c1->cd(4);
402  TH1F *diffTSmeanA_Depth1_HE = new TH1F("diffTSmeanA_Depth1_HE", "", 100, 0.0, 7.0);
403  for (int i = 1; i <= C3ff->GetXaxis()->GetNbins(); i++) {
404  for (int j = 1; j <= C3ff->GetYaxis()->GetNbins(); j++) {
405  if (C3ff->GetBinContent(i, j) != 0) {
406  double ccc1 = C3ff->GetBinContent(i, j);
407  diffTSmeanA_Depth1_HE->Fill(ccc1);
408  }
409  }
410  }
411  gPad->SetLogy();
412  diffTSmeanA_Depth1_HE->SetMarkerStyle(20);
413  diffTSmeanA_Depth1_HE->SetMarkerSize(0.4);
414  diffTSmeanA_Depth1_HE->GetYaxis()->SetLabelSize(0.04);
415  diffTSmeanA_Depth1_HE->SetXTitle("<TSmeanA> in each cell \b");
416  diffTSmeanA_Depth1_HE->SetMarkerColor(2);
417  diffTSmeanA_Depth1_HE->SetLineColor(2);
418  diffTSmeanA_Depth1_HE->Draw("");
419 
420  c1->Update();
421 
422  //=============================================================================================== 4
423  //======================================================================
424  //======================================================================
425  //================
426  //======================================================================
427  c1->Clear();
428  c1->Divide(2, 2);
429 
430  TH2F *ewod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HE");
431  TH2F *ewod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
432  ewod1->Sumw2();
433  ewod0->Sumw2();
434  // if(ewod0->IsA()->InheritsFrom("TH2F")){
435  TH2F *C4ff = (TH2F *)ewod1->Clone("C4ff");
436  C4ff->Divide(ewod1, ewod0, 1, 1, "B");
437  C4ff->Sumw2();
438  // }
439  c1->cd(1);
440  TH2F *fwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HE");
441  TH2F *fwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
442  fwod1->Sumw2();
443  fwod0->Sumw2();
444  // if(fwod0->IsA()->InheritsFrom("TH2F")){
445  TH2F *C4fz225 = (TH2F *)fwod1->Clone("C4fz225");
446  C4fz225->Divide(fwod1, fwod0, 1, 1, "B");
447  C4fz225->Sumw2();
448  // }
449  gPad->SetGridy();
450  gPad->SetGridx();
451  gPad->SetLogz();
452  C4fz225->SetMarkerStyle(20);
453  C4fz225->SetMarkerSize(0.4);
454  C4fz225->GetZaxis()->SetLabelSize(0.08);
455  C4fz225->SetXTitle("#eta \b");
456  C4fz225->SetYTitle("#phi \b");
457  C4fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth2 \b");
458  C4fz225->SetMarkerColor(2);
459  C4fz225->SetLineColor(2);
460  C4fz225->SetMaximum(1.000);
461  C4fz225->SetMinimum(0.0001);
462  C4fz225->Draw("COLZ");
463 
464  c1->cd(2);
465  TH1F *aaaaaa4 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
466  gPad->SetLogy();
467  aaaaaa4->SetMarkerStyle(20);
468  aaaaaa4->SetMarkerSize(0.8);
469  aaaaaa4->GetYaxis()->SetLabelSize(0.04);
470  aaaaaa4->SetXTitle("TSmeanA in each event & cell HE \b");
471  aaaaaa4->SetMarkerColor(2);
472  aaaaaa4->SetLineColor(2);
473  aaaaaa4->Draw("");
474 
475  c1->cd(3);
477  TH2F *Diffe_Depth2_HE = (TH2F *)C4ff->Clone("Diffe_Depth2_HE");
478  for (int i = 1; i <= C4ff->GetXaxis()->GetNbins(); i++) {
479  for (int j = 1; j <= C4ff->GetYaxis()->GetNbins(); j++) {
480  double ccc1 = C4ff->GetBinContent(i, j);
481  Diffe_Depth2_HE->SetBinContent(i, j, 0.);
482  if (ccc1 < 2.4 || ccc1 > 4.6)
483  Diffe_Depth2_HE->SetBinContent(i, j, ccc1);
484  }
485  }
486  gPad->SetGridy();
487  gPad->SetGridx();
488  gPad->SetLogz();
489  Diffe_Depth2_HE->SetMarkerStyle(20);
490  Diffe_Depth2_HE->SetMarkerSize(0.4);
491  Diffe_Depth2_HE->GetZaxis()->SetLabelSize(0.08);
492  // Diffe_Depth2_HE->SetTitle("any Error, HE Depth2 \n");
493  Diffe_Depth2_HE->SetXTitle("#eta \b");
494  Diffe_Depth2_HE->SetYTitle("#phi \b");
495  Diffe_Depth2_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth2 \b");
496  Diffe_Depth2_HE->SetMarkerColor(2);
497  Diffe_Depth2_HE->SetLineColor(2);
498  Diffe_Depth2_HE->Draw("COLZ");
499 
500  c1->cd(4);
501  TH1F *diffTSmeanA_Depth2_HE = new TH1F("diffTSmeanA_Depth2_HE", "", 100, 0.0, 7.0);
502  for (int i = 1; i <= C4ff->GetXaxis()->GetNbins(); i++) {
503  for (int j = 1; j <= C4ff->GetYaxis()->GetNbins(); j++) {
504  if (C4ff->GetBinContent(i, j) != 0) {
505  double ccc1 = C4ff->GetBinContent(i, j);
506  diffTSmeanA_Depth2_HE->Fill(ccc1);
507  }
508  }
509  }
510  gPad->SetLogy();
511  diffTSmeanA_Depth2_HE->SetMarkerStyle(20);
512  diffTSmeanA_Depth2_HE->SetMarkerSize(0.4);
513  diffTSmeanA_Depth2_HE->GetYaxis()->SetLabelSize(0.04);
514  diffTSmeanA_Depth2_HE->SetXTitle("<TSmeanA> in each cell \b");
515  diffTSmeanA_Depth2_HE->SetMarkerColor(2);
516  diffTSmeanA_Depth2_HE->SetLineColor(2);
517  diffTSmeanA_Depth2_HE->Draw("");
518 
519  c1->Update();
520 
521  //=============================================================================================== 5
522  //======================================================================
523  //======================================================================
524  //================
525  //======================================================================
526  c1->Clear();
527  c1->Divide(2, 2);
528 
529  TH2F *gwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HE");
530  TH2F *gwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
531  gwod1->Sumw2();
532  gwod0->Sumw2();
533  // if(gwod0->IsA()->InheritsFrom("TH2F")){
534  TH2F *C5ff = (TH2F *)gwod1->Clone("C5ff");
535  C5ff->Divide(gwod1, gwod0, 1, 1, "B");
536  C5ff->Sumw2();
537  // }
538  c1->cd(1);
539  TH2F *jwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA225_HE");
540  TH2F *jwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
541  jwod1->Sumw2();
542  jwod0->Sumw2();
543  // if(jwod0->IsA()->InheritsFrom("TH2F")){
544  TH2F *C5fz225 = (TH2F *)jwod1->Clone("C5fz225");
545  C5fz225->Divide(jwod1, jwod0, 1, 1, "B");
546  C5fz225->Sumw2();
547  // }
548  gPad->SetGridy();
549  gPad->SetGridx();
550  gPad->SetLogz();
551  C5fz225->SetMarkerStyle(20);
552  C5fz225->SetMarkerSize(0.4);
553  C5fz225->GetZaxis()->SetLabelSize(0.08);
554  C5fz225->SetXTitle("#eta \b");
555  C5fz225->SetYTitle("#phi \b");
556  C5fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HE Depth3 \b");
557  C5fz225->SetMarkerColor(2);
558  C5fz225->SetLineColor(2);
559  C5fz225->SetMaximum(1.000);
560  C5fz225->SetMinimum(0.0001);
561  C5fz225->Draw("COLZ");
562 
563  c1->cd(2);
564  TH1F *aaaaaa5 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
565  gPad->SetLogy();
566  aaaaaa5->SetMarkerStyle(20);
567  aaaaaa5->SetMarkerSize(0.8);
568  aaaaaa5->GetYaxis()->SetLabelSize(0.04);
569  aaaaaa5->SetXTitle("TSmeanA in each event & cell HE \b");
570  aaaaaa5->SetMarkerColor(2);
571  aaaaaa5->SetLineColor(2);
572  aaaaaa5->Draw("");
573 
574  c1->cd(3);
576  TH2F *Diffe_Depth3_HE = (TH2F *)C5ff->Clone("Diffe_Depth3_HE");
577  for (int i = 1; i <= C5ff->GetXaxis()->GetNbins(); i++) {
578  for (int j = 1; j <= C5ff->GetYaxis()->GetNbins(); j++) {
579  double ccc1 = C5ff->GetBinContent(i, j);
580  Diffe_Depth3_HE->SetBinContent(i, j, 0.);
581  if (ccc1 < 2.4 || ccc1 > 4.6)
582  Diffe_Depth3_HE->SetBinContent(i, j, ccc1);
583  }
584  }
585  gPad->SetGridy();
586  gPad->SetGridx();
587  gPad->SetLogz();
588  Diffe_Depth3_HE->SetMarkerStyle(20);
589  Diffe_Depth3_HE->SetMarkerSize(0.4);
590  Diffe_Depth3_HE->GetZaxis()->SetLabelSize(0.08);
591  // Diffe_Depth3_HE->SetTitle("any Error, HE Depth3 \n");
592  Diffe_Depth3_HE->SetXTitle("#eta \b");
593  Diffe_Depth3_HE->SetYTitle("#phi \b");
594  Diffe_Depth3_HE->SetZTitle("<TSmeanA> out 2.4-4.6 - HE Depth3 \b");
595  Diffe_Depth3_HE->SetMarkerColor(2);
596  Diffe_Depth3_HE->SetLineColor(2);
597  Diffe_Depth3_HE->Draw("COLZ");
598 
599  c1->cd(4);
600  TH1F *diffTSmeanA_Depth3_HE = new TH1F("diffTSmeanA_Depth3_HE", "", 100, 0.0, 7.0);
601  for (int i = 1; i <= C5ff->GetXaxis()->GetNbins(); i++) {
602  for (int j = 1; j <= C5ff->GetYaxis()->GetNbins(); j++) {
603  if (C5ff->GetBinContent(i, j) != 0) {
604  double ccc1 = C5ff->GetBinContent(i, j);
605  diffTSmeanA_Depth3_HE->Fill(ccc1);
606  }
607  }
608  }
609  gPad->SetLogy();
610  diffTSmeanA_Depth3_HE->SetMarkerStyle(20);
611  diffTSmeanA_Depth3_HE->SetMarkerSize(0.4);
612  diffTSmeanA_Depth3_HE->GetYaxis()->SetLabelSize(0.04);
613  diffTSmeanA_Depth3_HE->SetXTitle("<TSmeanA> in each cell \b");
614  diffTSmeanA_Depth3_HE->SetMarkerColor(2);
615  diffTSmeanA_Depth3_HE->SetLineColor(2);
616  diffTSmeanA_Depth3_HE->Draw("");
617 
618  c1->Update();
619 
620  //=============================================================================================== 6
621  //======================================================================
622  //======================================================================
623  //================
624  //======================================================================
625  c1->Clear();
626  c1->Divide(2, 2);
627 
628  TH2F *iwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HF");
629  TH2F *iwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
630  iwod1->Sumw2();
631  iwod0->Sumw2();
632  // if(iwod0->IsA()->InheritsFrom("TH2F")){
633  TH2F *C6ff = (TH2F *)iwod1->Clone("C6ff");
634  C6ff->Divide(iwod1, iwod0, 1, 1, "B");
635  C6ff->Sumw2();
636  // }
637  c1->cd(1);
638  TH2F *kwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HF");
639  TH2F *kwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
640  kwod1->Sumw2();
641  kwod0->Sumw2();
642  // if(kwod0->IsA()->InheritsFrom("TH2F")){
643  TH2F *C6fz225 = (TH2F *)kwod1->Clone("C6fz225");
644  C6fz225->Divide(kwod1, kwod0, 1, 1, "B");
645  C6fz225->Sumw2();
646  // }
647  gPad->SetGridy();
648  gPad->SetGridx();
649  gPad->SetLogz();
650  C6fz225->SetMarkerStyle(20);
651  C6fz225->SetMarkerSize(0.4);
652  // C6fz225->GetYaxis()->SetLabelSize(0.04);
653  C6fz225->GetZaxis()->SetLabelSize(0.08);
654  C6fz225->SetXTitle("#eta \b");
655  C6fz225->SetYTitle("#phi \b");
656  C6fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth1 \b");
657  C6fz225->SetMarkerColor(2);
658  C6fz225->SetLineColor(2);
659  C6fz225->SetMaximum(1.000);
660  C6fz225->SetMinimum(0.0001);
661  C6fz225->Draw("COLZ");
662 
663  c1->cd(2);
664  TH1F *aaaaaa6 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HF");
665  gPad->SetLogy();
666  aaaaaa6->SetMarkerStyle(20);
667  aaaaaa6->SetMarkerSize(0.8);
668  aaaaaa6->GetYaxis()->SetLabelSize(0.04);
669  aaaaaa6->SetXTitle("TSmeanA in each event & cell HF \b");
670  aaaaaa6->SetMarkerColor(2);
671  aaaaaa6->SetLineColor(2);
672  aaaaaa6->Draw("");
673 
674  c1->cd(3);
676  TH2F *Diffe_Depth1_HF = (TH2F *)C6ff->Clone("Diffe_Depth1_HF");
677  for (int i = 1; i <= C6ff->GetXaxis()->GetNbins(); i++) {
678  for (int j = 1; j <= C6ff->GetYaxis()->GetNbins(); j++) {
679  double ccc1 = C6ff->GetBinContent(i, j);
680  Diffe_Depth1_HF->SetBinContent(i, j, 0.);
681  if (ccc1 < 0.5 || ccc1 > 1.5)
682  Diffe_Depth1_HF->SetBinContent(i, j, ccc1);
683  }
684  }
685  gPad->SetGridy();
686  gPad->SetGridx();
687  gPad->SetLogz();
688  Diffe_Depth1_HF->SetMarkerStyle(20);
689  Diffe_Depth1_HF->SetMarkerSize(0.4);
690  Diffe_Depth1_HF->GetZaxis()->SetLabelSize(0.08);
691  // Diffe_Depth1_HF->SetTitle("any Error, HF Depth1 \n");
692  Diffe_Depth1_HF->SetXTitle("#eta \b");
693  Diffe_Depth1_HF->SetYTitle("#phi \b");
694  Diffe_Depth1_HF->SetZTitle("<TSmeanA> out 0.5-1.5 - HF Depth1 \b");
695  Diffe_Depth1_HF->SetMarkerColor(2);
696  Diffe_Depth1_HF->SetLineColor(2);
697  Diffe_Depth1_HF->Draw("COLZ");
698 
699  c1->cd(4);
700  TH1F *diffTSmeanA_Depth1_HF = new TH1F("diffTSmeanA_Depth1_HF", "", 100, 0.0, 2.0);
701  for (int i = 1; i <= C6ff->GetXaxis()->GetNbins(); i++) {
702  for (int j = 1; j <= C6ff->GetYaxis()->GetNbins(); j++) {
703  if (C6ff->GetBinContent(i, j) != 0) {
704  double ccc1 = C6ff->GetBinContent(i, j);
705  diffTSmeanA_Depth1_HF->Fill(ccc1);
706  }
707  }
708  }
709  gPad->SetLogy();
710  diffTSmeanA_Depth1_HF->SetMarkerStyle(20);
711  diffTSmeanA_Depth1_HF->SetMarkerSize(0.4);
712  diffTSmeanA_Depth1_HF->GetYaxis()->SetLabelSize(0.04);
713  diffTSmeanA_Depth1_HF->SetXTitle("<TSmeanA> in each cell \b");
714  diffTSmeanA_Depth1_HF->SetMarkerColor(2);
715  diffTSmeanA_Depth1_HF->SetLineColor(2);
716  diffTSmeanA_Depth1_HF->Draw("");
717 
718  c1->Update();
719 
720  //=============================================================================================== 7
721  //======================================================================
722  //======================================================================
723  //================
724  //======================================================================
725  c1->Clear();
726  c1->Divide(2, 2);
727 
728  TH2F *lwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HF");
729  TH2F *lwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
730  lwod1->Sumw2();
731  lwod0->Sumw2();
732  // if(lwod0->IsA()->InheritsFrom("TH2F")){
733  TH2F *C7ff = (TH2F *)lwod1->Clone("C7ff");
734  C7ff->Divide(lwod1, lwod0, 1, 1, "B");
735  C7ff->Sumw2();
736  // }
737  c1->cd(1);
738  TH2F *mwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HF");
739  TH2F *mwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
740  mwod1->Sumw2();
741  mwod0->Sumw2();
742  // if(mwod0->IsA()->InheritsFrom("TH2F")){
743  TH2F *C7fz225 = (TH2F *)mwod1->Clone("C7fz225");
744  C7fz225->Divide(mwod1, mwod0, 1, 1, "B");
745  C7fz225->Sumw2();
746  // }
747  gPad->SetGridy();
748  gPad->SetGridx();
749  gPad->SetLogz();
750  C7fz225->SetMarkerStyle(20);
751  C7fz225->SetMarkerSize(0.4);
752  C7fz225->GetZaxis()->SetLabelSize(0.08);
753  C7fz225->SetXTitle("#eta \b");
754  C7fz225->SetYTitle("#phi \b");
755  C7fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5 - HF Depth2 \b");
756  C7fz225->SetMarkerColor(2);
757  C7fz225->SetLineColor(2);
758  C7fz225->SetMaximum(1.000);
759  C7fz225->SetMinimum(0.0001);
760  C7fz225->Draw("COLZ");
761 
762  c1->cd(2);
763  TH1F *aaaaaa7 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HF");
764  gPad->SetLogy();
765  aaaaaa7->SetMarkerStyle(20);
766  aaaaaa7->SetMarkerSize(0.8);
767  aaaaaa7->GetYaxis()->SetLabelSize(0.04);
768  aaaaaa7->SetXTitle("TSmeanA in each event & cell HF \b");
769  aaaaaa7->SetMarkerColor(2);
770  aaaaaa7->SetLineColor(2);
771  aaaaaa7->Draw("");
772 
773  c1->cd(3);
775  TH2F *Diffe_Depth2_HF = (TH2F *)C7ff->Clone("Diffe_Depth2_HF");
776  for (int i = 1; i <= C7ff->GetXaxis()->GetNbins(); i++) {
777  for (int j = 1; j <= C7ff->GetYaxis()->GetNbins(); j++) {
778  double ccc1 = C7ff->GetBinContent(i, j);
779  Diffe_Depth2_HF->SetBinContent(i, j, 0.);
780  if (ccc1 < 0.5 || ccc1 > 1.5)
781  Diffe_Depth2_HF->SetBinContent(i, j, ccc1);
782  }
783  }
784  gPad->SetGridy();
785  gPad->SetGridx();
786  gPad->SetLogz();
787  Diffe_Depth2_HF->SetMarkerStyle(20);
788  Diffe_Depth2_HF->SetMarkerSize(0.4);
789  Diffe_Depth2_HF->GetZaxis()->SetLabelSize(0.08);
790  // Diffe_Depth2_HF->SetTitle("any Error, HF Depth2 \n");
791  Diffe_Depth2_HF->SetXTitle("#eta \b");
792  Diffe_Depth2_HF->SetYTitle("#phi \b");
793  Diffe_Depth2_HF->SetZTitle("<TSmeanA> out 0.5-1.5 -HF Depth2 \b");
794  Diffe_Depth2_HF->SetMarkerColor(2);
795  Diffe_Depth2_HF->SetLineColor(2);
796  Diffe_Depth2_HF->Draw("COLZ");
797 
798  c1->cd(4);
799  TH1F *diffTSmeanA_Depth2_HF = new TH1F("diffTSmeanA_Depth2_HF", "", 100, 0.0, 2.0);
800  for (int i = 1; i <= C7ff->GetXaxis()->GetNbins(); i++) {
801  for (int j = 1; j <= C7ff->GetYaxis()->GetNbins(); j++) {
802  if (C7ff->GetBinContent(i, j) != 0) {
803  double ccc1 = C7ff->GetBinContent(i, j);
804  diffTSmeanA_Depth2_HF->Fill(ccc1);
805  }
806  }
807  }
808  gPad->SetLogy();
809  diffTSmeanA_Depth2_HF->SetMarkerStyle(20);
810  diffTSmeanA_Depth2_HF->SetMarkerSize(0.4);
811  diffTSmeanA_Depth2_HF->GetYaxis()->SetLabelSize(0.04);
812  diffTSmeanA_Depth2_HF->SetXTitle("<TSmeanA> in each cell \b");
813  diffTSmeanA_Depth2_HF->SetMarkerColor(2);
814  diffTSmeanA_Depth2_HF->SetLineColor(2);
815  diffTSmeanA_Depth2_HF->Draw("");
816 
817  c1->Update();
818 
819  //=============================================================================================== 8
820  //======================================================================
821  //======================================================================
822  //================
823  //======================================================================
824  c1->Clear();
825  c1->Divide(2, 2);
826 
827  TH2F *nwod1 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HO");
828  TH2F *nwod0 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
829  nwod1->Sumw2();
830  nwod0->Sumw2();
831  // if(nwod0->IsA()->InheritsFrom("TH2F")){
832  TH2F *C8ff = (TH2F *)nwod1->Clone("C8ff");
833  C8ff->Divide(nwod1, nwod0, 1, 1, "B");
834  C8ff->Sumw2();
835  // }
836  c1->cd(1);
837  TH2F *owod1 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HO");
838  TH2F *owod0 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
839  owod1->Sumw2();
840  owod0->Sumw2();
841  // if(owod0->IsA()->InheritsFrom("TH2F")){
842  TH2F *C8fz225 = (TH2F *)owod1->Clone("C8fz225");
843  C8fz225->Divide(owod1, owod0, 1, 1, "B");
844  C8fz225->Sumw2();
845  // }
846  gPad->SetGridy();
847  gPad->SetGridx();
848  gPad->SetLogz();
849  C8fz225->SetMarkerStyle(20);
850  C8fz225->SetMarkerSize(0.4);
851  C8fz225->GetZaxis()->SetLabelSize(0.08);
852  C8fz225->SetXTitle("#eta \b");
853  C8fz225->SetYTitle("#phi \b");
854  C8fz225->SetZTitle("Rate for TSmeanA in each event & cell out 2.3-5.5 - HO Depth4 \b");
855  C8fz225->SetMarkerColor(2);
856  C8fz225->SetLineColor(2);
857  C8fz225->SetMaximum(1.000);
858  C8fz225->SetMinimum(0.000005);
859  C8fz225->Draw("COLZ");
860 
861  c1->cd(2);
862  TH1F *aaaaaa8 = (TH1F *)dir->FindObjectAny("h_TSmeanA_HO");
863  gPad->SetLogy();
864  aaaaaa8->SetMarkerStyle(20);
865  aaaaaa8->SetMarkerSize(0.8);
866  aaaaaa8->GetYaxis()->SetLabelSize(0.04);
867  aaaaaa8->SetXTitle("TSmeanA in each event & cell HO \b");
868  aaaaaa8->SetMarkerColor(2);
869  aaaaaa8->SetLineColor(2);
870  aaaaaa8->Draw("");
871 
872  c1->cd(3);
874  TH2F *Diffe_Depth4_HO = (TH2F *)C8ff->Clone("Diffe_Depth4_HO");
875  for (int i = 1; i <= C8ff->GetXaxis()->GetNbins(); i++) {
876  for (int j = 1; j <= C8ff->GetYaxis()->GetNbins(); j++) {
877  double ccc1 = C8ff->GetBinContent(i, j);
878  Diffe_Depth4_HO->SetBinContent(i, j, 0.);
879  if (ccc1 < 4.0 || ccc1 > 5.0)
880  Diffe_Depth4_HO->SetBinContent(i, j, ccc1);
881  }
882  }
883  gPad->SetGridy();
884  gPad->SetGridx();
885  gPad->SetLogz();
886  Diffe_Depth4_HO->SetMarkerStyle(20);
887  Diffe_Depth4_HO->SetMarkerSize(0.4);
888  Diffe_Depth4_HO->GetZaxis()->SetLabelSize(0.08);
889  // Diffe_Depth4_HO->SetTitle("any Error, HO Depth4 \n");
890  Diffe_Depth4_HO->SetXTitle("#eta \b");
891  Diffe_Depth4_HO->SetYTitle("#phi \b");
892  Diffe_Depth4_HO->SetZTitle("<TSmeanA> out 4.0-5.0 - HO Depth4 \b");
893  Diffe_Depth4_HO->SetMarkerColor(2);
894  Diffe_Depth4_HO->SetLineColor(2);
895  Diffe_Depth4_HO->Draw("COLZ");
896 
897  c1->cd(4);
898  TH1F *diffTSmeanA_Depth4_HO = new TH1F("diffTSmeanA_Depth4_HO", "", 100, 4.0, 5.0);
899  for (int i = 1; i <= C8ff->GetXaxis()->GetNbins(); i++) {
900  for (int j = 1; j <= C8ff->GetYaxis()->GetNbins(); j++) {
901  if (C8ff->GetBinContent(i, j) != 0) {
902  double ccc1 = C8ff->GetBinContent(i, j);
903  diffTSmeanA_Depth4_HO->Fill(ccc1);
904  }
905  }
906  }
907  gPad->SetLogy();
908  diffTSmeanA_Depth4_HO->SetMarkerStyle(20);
909  diffTSmeanA_Depth4_HO->SetMarkerSize(0.4);
910  diffTSmeanA_Depth4_HO->GetYaxis()->SetLabelSize(0.04);
911  diffTSmeanA_Depth4_HO->SetXTitle("<TSmeanA> in each cell \b");
912  diffTSmeanA_Depth4_HO->SetMarkerColor(2);
913  diffTSmeanA_Depth4_HO->SetLineColor(2);
914  diffTSmeanA_Depth4_HO->Draw("");
915 
916  c1->Update();
917 
918  //======================================================================
919  //======================================================================================== end
920  //======================================================================
921  //======================================================================
922  // close and delete all possible things:
923 
924  // psfile->Close();
925  psfile.Close();
926 
927  hfile1->Close();
928  // hfile1->Delete();
929  hfile1->Close();
930  // hfile1->Delete();
931 
932  // Exit Root
933  gSystem->Exit(0);
934  //======================================================================
935 }
int main(int argc, char *argv[])
Definition: ztsmean.cc:30