CMS 3D CMS Logo

Functions
zamplpedcorr.cc File Reference
#include <TChain.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TProfile.h>
#include <TFitResult.h>
#include <TFitResultPtr.h>
#include <TPaveStats.h>
#include <vector>
#include <string>
#include <iomanip>
#include <iostream>
#include <fstream>
#include "TH1.h"
#include "TH2.h"
#include "TCanvas.h"
#include "TROOT.h"
#include <TMath.h>
#include "TStyle.h"
#include "TSystem.h"
#include "TLegend.h"
#include "TText.h"
#include "TAxis.h"
#include "TFile.h"
#include "TLine.h"
#include "TGraph.h"
#include <TPaveText.h>
#include <TPostScript.h>
#include <TClass.h>

Go to the source code of this file.

Functions

int main (int argc, char *argv[])
 

Function Documentation

◆ main()

int main ( int  argc,
char *  argv[] 
)

===============================================================================================================================================================================================


variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk




















































































eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:







Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Sticking ADC counts in neughbouring TS for HB:

Sticking ADC counts in neughbouring TS for HE:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 38 of file zamplpedcorr.cc.

References alignmentValidation::c1, DeadROC_duringRun::dir, and AlCaHLTBitMon_QueryRunRegistry::string.

38  {
39  std::string dirnm = "Analyzer";
40  //======================================================================
41  printf("reco: gROOT Reset \n");
42  gROOT->Reset();
43  gROOT->SetStyle("Plain");
44  // gStyle->SetOptStat(0); // no statistics _or_
45  // gStyle->SetOptStat(11111111);
46  //gStyle->SetOptStat(1101);// name mean and rms
47  // gStyle->SetOptStat(0101);// name and entries
48  // gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
49  // gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
50  gStyle->SetOptStat(101110); // entries, mean, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
51  // gStyle->SetOptStat(100000);// over !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
52  //
53  //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
54  // gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
55  gStyle->SetOptFit(0010); // constant, mean and sigma only !!
56  // gStyle->SetOptFit(00011);// constant, mean and sigma only !!
57  // gStyle->SetOptFit(1101);
58  // gStyle->SetOptFit(1011);
59  //
60  //gStyle->SetStatX(0.98);
61  //gStyle->SetStatY(0.99);
62  //gStyle->SetStatW(0.30);
63  //gStyle->SetStatH(0.25);
64  gStyle->SetStatX(0.91);
65  gStyle->SetStatY(0.75);
66  gStyle->SetStatW(0.20);
67  gStyle->SetStatH(0.10);
68  // gStyle->SetStatH(0.35);
69  //
70 
71  // Float_t LeftOffset = 0.12;
72  //Float_t TopOffset = 0.12;
73  Float_t LeftOffset = 0.12;
74  Float_t TopOffset = 0.22;
75 
76  gStyle->SetLineWidth(1);
77  gStyle->SetErrorX(0);
78 
79  //---=[ Titles,Labels ]=-----------
80  gStyle->SetOptTitle(0); // title on/off
81  // gStyle->SetTitleColor(0); // title color
82  gStyle->SetTitleColor(1); // title color
83  // gStyle->SetTitleX(0.35); // title x-position
84  gStyle->SetTitleX(0.15); // title x-position
85  gStyle->SetTitleH(0.15); // title height
86  // gStyle->SetTitleW(0.53); // title width
87  gStyle->SetTitleW(0.60); // title width
88  gStyle->SetTitleFont(42); // title font
89  gStyle->SetTitleFontSize(0.07); // title font size
90 
91  gStyle->SetPalette(1);
92  //---=[ Histogram style ]=----------
93  // gStyle->SetHistFillColor(38);
94  // gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
95  // gStyle->SetFrameFillColor(5);// jeltyi
96  // gStyle->SetFrameFillColor(17);// seryi
97  // gStyle->SetFrameFillColor(18);// svetlo seryi
98  // gStyle->SetFrameFillColor(20);// svetlo kori4nev
99  // gStyle->SetFrameFillColor(33);// sine-seryi
100  // gStyle->SetFrameFillColor(40);// fiolet-seryi
101  // gStyle->SetFrameFillColor(23);// sv.kor
102 
103  //---=[ Pad style ]=----------------
104  gStyle->SetPadTopMargin(TopOffset);
105  gStyle->SetPadBottomMargin(LeftOffset);
106  gStyle->SetPadRightMargin(TopOffset);
107  gStyle->SetPadLeftMargin(LeftOffset);
108  //---=[ SetCanvasDef ]=----------------
109  //======================================================================
110  //
111  // Connect the input file and get the 2-d histogram in memory
112  //======================================================================
113  // TBrowser *b = new TBrowser
114 
115  // TFile *hfile1= new TFile("test.root", "READ");
116  //
117 
118  // TFile *hfile1= new TFile("testNZS.root", "READ");
119  // TFile *hfile1= new TFile("test.root", "READ");
120 
121  // TFile *hfile1= new TFile("newruns/Global_234034.root", "READ");
122  //
123 
124  // TFile *hfile1= new TFile("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/CMTweb/histos/Global_234556.root", "READ");
125  // TFile *hfile1= new TFile("Global_234034.root", "READ");
126  // TFile *hfile1= new TFile("test211006.root", "READ");
127  // TFile *hfile1= new TFile("test234457.root", "READ");
128 
129  // TFile *hfile1= new TFile("Global_51.root", "READ");
130  // TFile *hfile1= new TFile("Global_235698.root", "READ");
131 
132  // TFile *hfile1= new TFile("Global_39.root", "READ");
133  // TFile *hfile1= new TFile("test237165.root", "READ");
134  //
135  // TFile *hfile1= new TFile("Laser_238187.root", "READ");
136  // TFile *hfile1= new TFile("Laser_238183.root", "READ");
137 
138  // TFile *hfile1= new TFile("Global_255031.root", "READ");
139 
140  // TFile *hfile1= new TFile("Global_256001.root", "READ");
141  // TFile *hfile1= new TFile("Global_256167.root", "READ");
142  // TFile *hfile1= new TFile("Global_256348.root", "READ");
143  // TFile *hfile1= new TFile("Global_256630.root", "READ");
144 
145  // TFile *hfile1= new TFile("../PYTHON_runlist_test/Global_283884_1.root", "READ");
146  // TFile *hfile1= new TFile("Global_test.root", "READ");
147 
148  // TFile *hfile1= new TFile("LED_280702.root", "READ");
149  // TFile *hfile2= new TFile("LED_287824.root", "READ");
150 
151  // TFile *hfile1= new TFile("LED_284352.root", "READ");
152  // TFile *hfile1= new TFile("LEDtest.root", "READ");
153  // TFile *hfile1 = new TFile("Global_346445.root", "READ");
154 
155  TFile *hfile1 = new TFile("Global_362365.root", "READ");
156  // TH1D *hist1(nullptr);
157  // hist1 = (TH1D *)dir->FindObjectAny("h_mapDepth1_HE");
158 
159  // TFile *hfile2= new TFile("LED_284902.root", "READ");
160  // TFile *hfile2= new TFile("LED_284499.root", "READ");
161  // TFile *hfile2= new TFile("LED_284352.root", "READ");
162 
163  // TFile *hfile2= new TFile("LED_286590.root", "READ");
164 
165  // getchar();
166  //
167  TPostScript psfile("zamplpedcorr.ps", 111);
168 
169  //
170 
171  TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
172 
173  hfile1->ls();
174  TDirectory *dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
175 
176  //=============================================================================================== 1
177  //======================================================================
178  //======================================================================
179  //================
180  //======================================================================
181  c1->Clear();
182  c1->Divide(2, 2);
183 
184  c1->cd(1);
185  TH1F *twod1 = (TH1F *)dir->FindObjectAny("h_pedvsampl_HB");
186  TH1F *twod0 = (TH1F *)dir->FindObjectAny("h_pedvsampl0_HB");
187  twod1->Sumw2();
188  twod0->Sumw2();
189  // if(twod0->IsA()->InheritsFrom("TH1F")){
190  TH1F *Cefz225 = (TH1F *)twod1->Clone("Cefz225");
191  Cefz225->Divide(twod1, twod0, 1, 1, "B");
192  Cefz225->Sumw2();
193  // }
194  gPad->SetGridy();
195  gPad->SetGridx();
196  Cefz225->SetMarkerStyle(20);
197  Cefz225->SetMarkerSize(0.4);
198  Cefz225->SetXTitle("Pedestals \b");
199  Cefz225->SetYTitle("<A> \b");
200  Cefz225->SetMarkerColor(2);
201  Cefz225->SetLineColor(2);
202  // Cefz225->SetMaximum(1.000);
203  // Cefz225->SetMinimum(0.0001);
204  Cefz225->Draw("COLZ");
205 
206  c1->cd(2);
207  TH1F *twod61 = (TH1F *)dir->FindObjectAny("h_pedwvsampl_HB");
208  TH1F *twod60 = (TH1F *)dir->FindObjectAny("h_pedwvsampl0_HB");
209  twod61->Sumw2();
210  twod60->Sumw2();
211  // if(twod60->IsA()->InheritsFrom("TH1F")){
212  TH1F *Cefz226 = (TH1F *)twod61->Clone("Cefz226");
213  Cefz226->Divide(twod61, twod60, 1, 1, "B");
214  Cefz226->Sumw2();
215  // }
216  gPad->SetGridy();
217  gPad->SetGridx();
218  Cefz226->SetMarkerStyle(20);
219  Cefz226->SetMarkerSize(0.4);
220  Cefz226->SetXTitle("Width_Pedestals \b");
221  Cefz226->SetYTitle("<A> \b");
222  Cefz226->SetMarkerColor(2);
223  Cefz226->SetLineColor(2);
224  // Cefz226->SetMaximum(1.000);
225  // Cefz226->SetMinimum(0.0001);
226  Cefz226->Draw("COLZ");
227 
228  c1->cd(3);
229  TH1F *twod71 = (TH1F *)dir->FindObjectAny("h_amplvsped_HB");
230  TH1F *twod70 = (TH1F *)dir->FindObjectAny("h_amplvsped0_HB");
231  twod71->Sumw2();
232  twod70->Sumw2();
233  // if(twod70->IsA()->InheritsFrom("TH1F")){
234  TH1F *Cefz227 = (TH1F *)twod71->Clone("Cefz227");
235  Cefz227->Divide(twod71, twod70, 1, 1, "B");
236  Cefz227->Sumw2();
237  // }
238  gPad->SetGridy();
239  gPad->SetGridx();
240  Cefz227->SetMarkerStyle(20);
241  Cefz227->SetMarkerSize(0.4);
242  Cefz227->SetXTitle("Amplitude \b");
243  Cefz227->SetYTitle("<Pedestals> \b");
244  Cefz227->SetMarkerColor(2);
245  Cefz227->SetLineColor(2);
246  // Cefz227->SetMaximum(1.000);
247  // Cefz227->SetMinimum(0.0001);
248  Cefz227->Draw("COLZ");
249 
250  c1->cd(4);
251  TH1F *twod81 = (TH1F *)dir->FindObjectAny("h_amplvspedw_HB");
252  TH1F *twod80 = (TH1F *)dir->FindObjectAny("h_amplvsped0_HB");
253  twod81->Sumw2();
254  twod80->Sumw2();
255  // if(twod80->IsA()->InheritsFrom("TH1F")){
256  TH1F *Cefz228 = (TH1F *)twod81->Clone("Cefz228");
257  Cefz228->Divide(twod81, twod80, 1, 1, "B");
258  Cefz228->Sumw2();
259  // }
260  gPad->SetGridy();
261  gPad->SetGridx();
262  Cefz228->SetMarkerStyle(20);
263  Cefz228->SetMarkerSize(0.4);
264  Cefz228->SetXTitle("Amplitude \b");
265  Cefz228->SetYTitle("<Width_Pedestals> \b");
266  Cefz228->SetMarkerColor(2);
267  Cefz228->SetLineColor(2);
268  // Cefz228->SetMaximum(1.000);
269  // Cefz228->SetMinimum(0.0001);
270  Cefz228->Draw("COLZ");
271 
272  c1->Update();
273 
274  //========================================================================================= 2
275  //======================================================================
276  //======================================================================
277  //================
278  //======================================================================
279  c1->Clear();
280  c1->Divide(2, 2);
281 
282  c1->cd(1);
283  TH2F *two11 = (TH2F *)dir->FindObjectAny("h2_pedvsampl_HB");
284  gPad->SetGridy();
285  gPad->SetGridx();
286  two11->SetMarkerStyle(20);
287  two11->SetMarkerSize(0.4);
288  two11->SetXTitle("Pedestals HB \b");
289  two11->SetYTitle("Amplitude HB\b");
290  two11->SetMarkerColor(2);
291  two11->SetLineColor(2);
292  // gPad->SetLogy();
293  two11->Draw("COLZ");
294 
295  c1->cd(2);
296  TH2F *two12 = (TH2F *)dir->FindObjectAny("h2_pedwvsampl_HB");
297  gPad->SetGridy();
298  gPad->SetGridx();
299  two12->SetMarkerStyle(20);
300  two12->SetMarkerSize(0.4);
301  two12->SetXTitle("Width_Pedestals HB \b");
302  two12->SetYTitle("Amplitude HB\b");
303  two12->SetMarkerColor(2);
304  two12->SetLineColor(2);
305  // gPad->SetLogy();
306  two12->Draw("COLZ");
307 
308  c1->cd(3);
309  TH2F *two22 = (TH2F *)dir->FindObjectAny("h2_amplvsped_HB");
310  gPad->SetGridy();
311  gPad->SetGridx();
312  two22->SetMarkerStyle(20);
313  two22->SetMarkerSize(0.4);
314  two22->SetYTitle("Pedestals HB \b");
315  two22->SetXTitle("Amplitude HB\b");
316  two22->SetMarkerColor(2);
317  two22->SetLineColor(2);
318  two22->Draw("COLZ");
319 
320  c1->cd(4);
321  TH2F *two23 = (TH2F *)dir->FindObjectAny("h2_amplvspedw_HB");
322  gPad->SetGridy();
323  gPad->SetGridx();
324  two23->SetMarkerStyle(20);
325  two23->SetMarkerSize(0.4);
326  two23->SetYTitle("Width_Pedestals HB \b");
327  two23->SetXTitle("Amplitude HB\b");
328  two23->SetMarkerColor(2);
329  two23->SetLineColor(2);
330  two23->Draw("COLZ");
331 
332  c1->Update();
333 
334  //========================================================================================= 3
335  //======================================================================
336  //======================================================================
337  //================
338  //======================================================================
339  c1->Clear();
340  c1->Divide(2, 2);
341 
342  c1->cd(1);
343  TH1F *aaaaaa1 = (TH1F *)dir->FindObjectAny("h_sumADCAmplLS1copy2");
344  gPad->SetLogy();
345  aaaaaa1->SetMarkerStyle(20);
346  aaaaaa1->SetMarkerSize(0.8);
347  aaaaaa1->GetYaxis()->SetLabelSize(0.04);
348  aaaaaa1->SetXTitle("<A>(ev.in LS) in LSs & channels - HB depth1\b");
349  aaaaaa1->SetMarkerColor(4);
350  aaaaaa1->SetLineColor(0);
351  aaaaaa1->Draw("Error");
352 
353  c1->cd(2);
354  TH1F *aaaaaa2 = (TH1F *)dir->FindObjectAny("h_sumADCAmplLS1copy3");
355  gPad->SetLogy();
356  aaaaaa2->SetMarkerStyle(20);
357  aaaaaa2->SetMarkerSize(0.8);
358  aaaaaa2->GetYaxis()->SetLabelSize(0.04);
359  aaaaaa2->SetXTitle("<A>(ev.in LS) in LSs & channels - HB depth1\b");
360  aaaaaa2->SetMarkerColor(4);
361  aaaaaa2->SetLineColor(0);
362  aaaaaa2->Draw("Error");
363 
364  c1->cd(3);
365  TH1F *aaaaaa3 = (TH1F *)dir->FindObjectAny("h_sumADCAmplLS1copy4");
366  gPad->SetLogy();
367  aaaaaa3->SetMarkerStyle(20);
368  aaaaaa3->SetMarkerSize(0.8);
369  aaaaaa3->GetYaxis()->SetLabelSize(0.04);
370  aaaaaa3->SetXTitle("<A>(ev.in LS) in LSs & channels - HB depth1\b");
371  aaaaaa3->SetMarkerColor(4);
372  aaaaaa3->SetLineColor(0);
373  aaaaaa3->Draw("Error");
374  c1->cd(4);
375  TH1F *aaaaaa4 = (TH1F *)dir->FindObjectAny("h_sumADCAmplLS1copy5");
376  gPad->SetLogy();
377  aaaaaa4->SetMarkerStyle(20);
378  aaaaaa4->SetMarkerSize(0.8);
379  aaaaaa4->GetYaxis()->SetLabelSize(0.04);
380  aaaaaa4->SetXTitle("<A>(ev.in LS) in LSs & channels - HB depth1\b");
381  aaaaaa4->SetMarkerColor(4);
382  aaaaaa4->SetLineColor(0);
383  aaaaaa4->Draw("Error");
384 
385  c1->Update();
386 
387  //=============================================================================================== 4
388  //======================================================================
389  //======================================================================
390  //================
391  //======================================================================
392  c1->Clear();
393  c1->Divide(2, 2);
394  /*
395  c1->cd(1);
396  TH1F *twodhe1= (TH1F*)dir->FindObjectAny("h_pedvsampl_HE");
397  TH1F *twodhe0= (TH1F*)dir->FindObjectAny("h_pedvsampl0_HE");
398  twodhe1->Sumw2();
399  twodhe0->Sumw2();
400 // if(twodhe0->IsA()->InheritsFrom("TH1F")){
401  TH1F* Cefzhe225= (TH1F*)twodhe1->Clone("Cefzhe225");
402  Cefzhe225->Divide(twodhe1,twodhe0, 1, 1, "B");
403  Cefzhe225->Sumw2();
404 // }
405  gPad->SetGridy();
406  gPad->SetGridx();
407  Cefzhe225->SetMarkerStyle(20);
408  Cefzhe225->SetMarkerSize(0.4);
409  Cefzhe225->SetXTitle("Pedestals \b");
410  Cefzhe225->SetYTitle("<A> \b");
411  Cefzhe225->SetMarkerColor(2);
412  Cefzhe225->SetLineColor(2);
413  // Cefzhe225->SetMaximum(1.000);
414  // Cefzhe225->SetMinimum(0.0001);
415  Cefzhe225->Draw("COLZ");
416 
417  c1->cd(2);
418  TH1F *twodhe61= (TH1F*)dir->FindObjectAny("h_pedwvsampl_HE");
419  TH1F *twodhe60= (TH1F*)dir->FindObjectAny("h_pedwvsampl0_HE");
420  twodhe61->Sumw2();
421  twodhe60->Sumw2();
422 // if(twodhe60->IsA()->InheritsFrom("TH1F")){
423  TH1F* Cefzhe226= (TH1F*)twodhe61->Clone("Cefzhe226");
424  Cefzhe226->Divide(twodhe61,twodhe60, 1, 1, "B");
425  Cefzhe226->Sumw2();
426 // }
427  gPad->SetGridy();
428  gPad->SetGridx();
429  Cefzhe226->SetMarkerStyle(20);
430  Cefzhe226->SetMarkerSize(0.4);
431  Cefzhe226->SetXTitle("Width_Pedestals \b");
432  Cefzhe226->SetYTitle("<A> \b");
433  Cefzhe226->SetMarkerColor(2);
434  Cefzhe226->SetLineColor(2);
435  // Cefzhe226->SetMaximum(1.000);
436  // Cefzhe226->SetMinimum(0.0001);
437  Cefzhe226->Draw("COLZ");
438 
439  c1->cd(3);
440  TH1F *twodhe71= (TH1F*)dir->FindObjectAny("h_amplvsped_HE");
441  TH1F *twodhe70= (TH1F*)dir->FindObjectAny("h_amplvsped0_HE");
442  twodhe71->Sumw2();
443  twodhe70->Sumw2();
444 // if(twodhe70->IsA()->InheritsFrom("TH1F")){
445  TH1F* Cefzhe227= (TH1F*)twodhe71->Clone("Cefzhe227");
446  Cefzhe227->Divide(twodhe71,twodhe70, 1, 1, "B");
447  Cefzhe227->Sumw2();
448 // }
449  gPad->SetGridy();
450  gPad->SetGridx();
451  Cefzhe227->SetMarkerStyle(20);
452  Cefzhe227->SetMarkerSize(0.4);
453  Cefzhe227->SetXTitle("Amplitude \b");
454  Cefzhe227->SetYTitle("<Pedestals> \b");
455  Cefzhe227->SetMarkerColor(2);
456  Cefzhe227->SetLineColor(2);
457  // Cefzhe227->SetMaximum(1.000);
458  // Cefzhe227->SetMinimum(0.0001);
459  Cefzhe227->Draw("COLZ");
460 
461  c1->cd(4);
462  TH1F *twodhe81= (TH1F*)dir->FindObjectAny("h_amplvspedw_HE");
463  TH1F *twodhe80= (TH1F*)dir->FindObjectAny("h_amplvsped0_HE");
464  twodhe81->Sumw2();
465  twodhe80->Sumw2();
466 // if(twodhe80->IsA()->InheritsFrom("TH1F")){
467  TH1F* Cefzhe228= (TH1F*)twodhe81->Clone("Cefzhe228");
468  Cefzhe228->Divide(twodhe81,twodhe80, 1, 1, "B");
469  Cefzhe228->Sumw2();
470 // }
471  gPad->SetGridy();
472  gPad->SetGridx();
473  Cefzhe228->SetMarkerStyle(20);
474  Cefzhe228->SetMarkerSize(0.4);
475  Cefzhe228->SetXTitle("Amplitude \b");
476  Cefzhe228->SetYTitle("<Width_Pedestals> \b");
477  Cefzhe228->SetMarkerColor(2);
478  Cefzhe228->SetLineColor(2);
479  // Cefzhe228->SetMaximum(1.000);
480  // Cefzhe228->SetMinimum(0.0001);
481  Cefzhe228->Draw("COLZ");
482 
483  */
484  c1->Update();
485 
486  //=============================================================================================== 5
487  //======================================================================
488  //======================================================================
489  //================
490  //======================================================================
491  c1->Clear();
492  c1->Divide(2, 2);
493  /*
494  c1->cd(1);
495  TH1F *twodhf1= (TH1F*)dir->FindObjectAny("h_pedvsampl_HF");
496  TH1F *twodhf0= (TH1F*)dir->FindObjectAny("h_pedvsampl0_HF");
497  twodhf1->Sumw2();
498  twodhf0->Sumw2();
499 // if(twodhf0->IsA()->InheritsFrom("TH1F")){
500  TH1F* Cefzhf225= (TH1F*)twodhf1->Clone("Cefzhf225");
501  Cefzhf225->Divide(twodhf1,twodhf0, 1, 1, "B");
502  Cefzhf225->Sumw2();
503 // }
504  gPad->SetGridy();
505  gPad->SetGridx();
506  Cefzhf225->SetMarkerStyle(20);
507  Cefzhf225->SetMarkerSize(0.4);
508  Cefzhf225->SetXTitle("Pedestals \b");
509  Cefzhf225->SetYTitle("<A> \b");
510  Cefzhf225->SetMarkerColor(2);
511  Cefzhf225->SetLineColor(2);
512  // Cefzhf225->SetMaximum(1.000);
513  // Cefzhf225->SetMinimum(0.0001);
514  Cefzhf225->Draw("COLZ");
515 
516  c1->cd(2);
517  TH1F *twodhf61= (TH1F*)dir->FindObjectAny("h_pedwvsampl_HF");
518  TH1F *twodhf60= (TH1F*)dir->FindObjectAny("h_pedwvsampl0_HF");
519  twodhf61->Sumw2();
520  twodhf60->Sumw2();
521 // if(twodhf60->IsA()->InheritsFrom("TH1F")){
522  TH1F* Cefzhf226= (TH1F*)twodhf61->Clone("Cefzhf226");
523  Cefzhf226->Divide(twodhf61,twodhf60, 1, 1, "B");
524  Cefzhf226->Sumw2();
525 // }
526  gPad->SetGridy();
527  gPad->SetGridx();
528  Cefzhf226->SetMarkerStyle(20);
529  Cefzhf226->SetMarkerSize(0.4);
530  Cefzhf226->SetXTitle("Width_Pedestals \b");
531  Cefzhf226->SetYTitle("<A> \b");
532  Cefzhf226->SetMarkerColor(2);
533  Cefzhf226->SetLineColor(2);
534  // Cefzhf226->SetMaximum(1.000);
535  // Cefzhf226->SetMinimum(0.0001);
536  Cefzhf226->Draw("COLZ");
537 
538  c1->cd(3);
539  TH1F *twodhf71= (TH1F*)dir->FindObjectAny("h_amplvsped_HF");
540  TH1F *twodhf70= (TH1F*)dir->FindObjectAny("h_amplvsped0_HF");
541  twodhf71->Sumw2();
542  twodhf70->Sumw2();
543 // if(twodhf70->IsA()->InheritsFrom("TH1F")){
544  TH1F* Cefzhf227= (TH1F*)twodhf71->Clone("Cefzhf227");
545  Cefzhf227->Divide(twodhf71,twodhf70, 1, 1, "B");
546  Cefzhf227->Sumw2();
547 // }
548  gPad->SetGridy();
549  gPad->SetGridx();
550  Cefzhf227->SetMarkerStyle(20);
551  Cefzhf227->SetMarkerSize(0.4);
552  Cefzhf227->SetXTitle("Amplitude \b");
553  Cefzhf227->SetYTitle("<Pedestals> \b");
554  Cefzhf227->SetMarkerColor(2);
555  Cefzhf227->SetLineColor(2);
556  // Cefzhf227->SetMaximum(1.000);
557  // Cefzhf227->SetMinimum(0.0001);
558  Cefzhf227->Draw("COLZ");
559 
560  c1->cd(4);
561  TH1F *twodhf81= (TH1F*)dir->FindObjectAny("h_amplvspedw_HF");
562  TH1F *twodhf80= (TH1F*)dir->FindObjectAny("h_amplvsped0_HF");
563  twodhf81->Sumw2();
564  twodhf80->Sumw2();
565 // if(twodhf80->IsA()->InheritsFrom("TH1F")){
566  TH1F* Cefzhf228= (TH1F*)twodhf81->Clone("Cefzhf228");
567  Cefzhf228->Divide(twodhf81,twodhf80, 1, 1, "B");
568  Cefzhf228->Sumw2();
569 // }
570  gPad->SetGridy();
571  gPad->SetGridx();
572  Cefzhf228->SetMarkerStyle(20);
573  Cefzhf228->SetMarkerSize(0.4);
574  Cefzhf228->SetXTitle("Amplitude \b");
575  Cefzhf228->SetYTitle("<Width_Pedestals> \b");
576  Cefzhf228->SetMarkerColor(2);
577  Cefzhf228->SetLineColor(2);
578  // Cefzhf228->SetMaximum(1.000);
579  // Cefzhf228->SetMinimum(0.0001);
580  Cefzhf228->Draw("COLZ");
581 
582  */
583 
584  c1->Update();
585 
586  //=============================================================================================== 6
587  //======================================================================
588  //======================================================================
589  //================
590  //======================================================================
591  c1->Clear();
592  c1->Divide(2, 2);
593  /*
594  c1->cd(1);
595  TH1F *twodho1= (TH1F*)dir->FindObjectAny("h_pedvsampl_HO");
596  TH1F *twodho0= (TH1F*)dir->FindObjectAny("h_pedvsampl0_HO");
597  twodho1->Sumw2();
598  twodho0->Sumw2();
599 // if(twodho0->IsA()->InheritsFrom("TH1F")){
600  TH1F* Cefzho225= (TH1F*)twodho1->Clone("Cefzho225");
601  Cefzho225->Divide(twodho1,twodho0, 1, 1, "B");
602  Cefzho225->Sumw2();
603 // }
604  gPad->SetGridy();
605  gPad->SetGridx();
606  Cefzho225->SetMarkerStyle(20);
607  Cefzho225->SetMarkerSize(0.4);
608  Cefzho225->SetXTitle("Pedestals \b");
609  Cefzho225->SetYTitle("<A> \b");
610  Cefzho225->SetMarkerColor(2);
611  Cefzho225->SetLineColor(2);
612  // Cefzho225->SetMaximum(1.000);
613  // Cefzho225->SetMinimum(0.0001);
614  Cefzho225->Draw("COLZ");
615 
616  c1->cd(2);
617  TH1F *twodho61= (TH1F*)dir->FindObjectAny("h_pedwvsampl_HO");
618  TH1F *twodho60= (TH1F*)dir->FindObjectAny("h_pedwvsampl0_HO");
619  twodho61->Sumw2();
620  twodho60->Sumw2();
621 // if(twodho60->IsA()->InheritsFrom("TH1F")){
622  TH1F* Cefzho226= (TH1F*)twodho61->Clone("Cefzho226");
623  Cefzho226->Divide(twodho61,twodho60, 1, 1, "B");
624  Cefzho226->Sumw2();
625 // }
626  gPad->SetGridy();
627  gPad->SetGridx();
628  Cefzho226->SetMarkerStyle(20);
629  Cefzho226->SetMarkerSize(0.4);
630  Cefzho226->SetXTitle("Width_Pedestals \b");
631  Cefzho226->SetYTitle("<A> \b");
632  Cefzho226->SetMarkerColor(2);
633  Cefzho226->SetLineColor(2);
634  // Cefzho226->SetMaximum(1.000);
635  // Cefzho226->SetMinimum(0.0001);
636  Cefzho226->Draw("COLZ");
637 
638  c1->cd(3);
639  TH1F *twodho71= (TH1F*)dir->FindObjectAny("h_amplvsped_HO");
640  TH1F *twodho70= (TH1F*)dir->FindObjectAny("h_amplvsped0_HO");
641  twodho71->Sumw2();
642  twodho70->Sumw2();
643 // if(twodho70->IsA()->InheritsFrom("TH1F")){
644  TH1F* Cefzho227= (TH1F*)twodho71->Clone("Cefzho227");
645  Cefzho227->Divide(twodho71,twodho70, 1, 1, "B");
646  Cefzho227->Sumw2();
647 // }
648  gPad->SetGridy();
649  gPad->SetGridx();
650  Cefzho227->SetMarkerStyle(20);
651  Cefzho227->SetMarkerSize(0.4);
652  Cefzho227->SetXTitle("Amplitude \b");
653  Cefzho227->SetYTitle("<Pedestals> \b");
654  Cefzho227->SetMarkerColor(2);
655  Cefzho227->SetLineColor(2);
656  // Cefzho227->SetMaximum(1.000);
657  // Cefzho227->SetMinimum(0.0001);
658  Cefzho227->Draw("COLZ");
659 
660  c1->cd(4);
661  TH1F *twodho81= (TH1F*)dir->FindObjectAny("h_amplvspedw_HO");
662  TH1F *twodho80= (TH1F*)dir->FindObjectAny("h_amplvsped0_HO");
663  twodho81->Sumw2();
664  twodho80->Sumw2();
665 // if(twodho80->IsA()->InheritsFrom("TH1F")){
666  TH1F* Cefzho228= (TH1F*)twodho81->Clone("Cefzho228");
667  Cefzho228->Divide(twodho81,twodho80, 1, 1, "B");
668  Cefzho228->Sumw2();
669 // }
670  gPad->SetGridy();
671  gPad->SetGridx();
672  Cefzho228->SetMarkerStyle(20);
673  Cefzho228->SetMarkerSize(0.4);
674  Cefzho228->SetXTitle("Amplitude \b");
675  Cefzho228->SetYTitle("<Width_Pedestals> \b");
676  Cefzho228->SetMarkerColor(2);
677  Cefzho228->SetLineColor(2);
678  // Cefzho228->SetMaximum(1.000);
679  // Cefzho228->SetMinimum(0.0001);
680  Cefzho228->Draw("COLZ");
681 
682  */
683 
684  c1->Update();
685 
686  //======================================================================
687  //======================================================================================== end
688  //======================================================================
689  //======================================================================
690  // close and delete all possible things:
691 
692  // psfile->Close();
693  psfile.Close();
694 
695  hfile1->Close();
696  // hfile1->Delete();
697  hfile1->Close();
698  // hfile1->Delete();
699 
700  // Exit Root
701  gSystem->Exit(0);
702  //======================================================================
703 }