Hall-D Software  alpha
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
BCAL_SiPM_saturation.C
Go to the documentation of this file.
1 void BCAL_SiPM_saturation(void)
2 {
3 // Plot histograms created with the BCAL_SiPM_saturation pluggin
4 //
5 
6 #include <TRandom.h>
7 
8 gROOT->Reset();
9 //TTree *Bfield = (TTree *) gROOT->FindObject("Bfield");
10 gStyle->SetPalette(1,0);
11 // gStyle->SetOptStat(kFALSE);
12 gStyle->SetOptStat(111111);
13 gStyle->SetPadRightMargin(0.15);
14 gStyle->SetPadLeftMargin(0.15);
15 gStyle->SetPadBottomMargin(0.15);
16 //
17 
18  Double_t xmin=0.;
19  Double_t xmax=0.1;
20  Double_t ymin=0.;
21  Double_t ymax=100;
22 
23  // read histograms from file
24 
25  TString histfile_name;
26  TString title;
27 
28  title="gensat_correct0_030746";
29  histfile_name= "hd_root_particle_gun_Photons_"+title+"_000";
30  // histfile_name= "hd_root_particle_gun_Photons_030731";
31  // histfile_name= "BCAL_SiPM_saturation_042433_001";
32  // title = "Data Run 42389_001";
33  // histfile_name= "BCAL_SiPM_data_42389_001";
34  TFile *histfile = new TFile(histfile_name+".root","read");
35 
36 
37  TH1I *Ethrown = (TH1I*)gDirectory->FindObjectAny("Ethrown");
38  if( !Ethrown) { cout << "Can't find Ethrown!" << endl; return; }
39  TH1I *Eshower = (TH1I*)gDirectory->FindObjectAny("Eshower");
40  if( !Eshower) { cout << "Can't find Eshower!" << endl; return; }
41  TH1I *ThrownTheta = (TH1I*)gDirectory->FindObjectAny("Thrown Theta");
42  if( !ThrownTheta) { cout << "Can't find Thrown Theta!" << endl; return; }
43  TH1I *NCell = (TH1I*)gDirectory->FindObjectAny("NCell");
44  if( !NCell) { cout << "Can't find NCell!" << endl; return; }
45  TH1I *layer = (TH1I*)gDirectory->FindObjectAny("layer");
46  if( !layer) { cout << "Can't find layer!" << endl; return; }
47  TH1I *Ept = (TH1I*)gDirectory->FindObjectAny("Ept");
48  if( !Ept) { cout << "Can't find Ept!" << endl; return; }
49  TH1I *EcalcDiff = (TH1I*)gDirectory->FindObjectAny("Ecalc-Ept");
50  if( !EcalcDiff) { cout << "Can't find Ecalc-Ept!" << endl; return; }
51 
52  TH2I *Eshower_vs_Ethrown = (TH2I*)gDirectory->FindObjectAny("Eshower_vs_Ethrown");
53  if( !Eshower_vs_Ethrown) { cout << "Can't find Eshower_vs_Ethrown!" << endl; return; }
54  TH2I *EDiff_vs_Ethrown = (TH2I*)gDirectory->FindObjectAny("EDiff_vs_Ethrown");
55  if( !EDiff_vs_Ethrown) { cout << "Can't find EDiff_vs_Ethrown!" << endl; return; }
56  TH2I *EDiffoverEthrown_vs_Ethrown = (TH2I*)gDirectory->FindObjectAny("EDiff/Ethrown_vs_Ethrown");
57  if( !EDiffoverEthrown_vs_Ethrown) { cout << "Can't find EDiff/Ethrown_vs_Ethrown!" << endl; return; }
58 
59  TH1I *Hit_pulse_peak_layer1 = (TH1I*)gDirectory->FindObjectAny("Hit pulse_peak layer=1");
60  if( !Hit_pulse_peak_layer1) { cout << "Can't find Hit pulse_peak layer=1!" << endl; return; }
61  TH1I *Hit_pulse_peak_layer2 = (TH1I*)gDirectory->FindObjectAny("Hit pulse_peak layer=2");
62  if( !Hit_pulse_peak_layer2) { cout << "Can't find Hit pulse_peak layer=2!" << endl; return; }
63  TH1I *Hit_pulse_peak_layer3 = (TH1I*)gDirectory->FindObjectAny("Hit pulse_peak layer=3");
64  if( !Hit_pulse_peak_layer3) { cout << "Can't find Hit pulse_peak layer=3!" << endl; return; }
65  TH1I *Hit_pulse_peak_layer4 = (TH1I*)gDirectory->FindObjectAny("Hit pulse_peak layer=4");
66  if( !Hit_pulse_peak_layer4) { cout << "Can't find Hit pulse_peak layer=4!" << endl; return; }
67 
68  TH1I *Hit_integral_layer1 = (TH1I*)gDirectory->FindObjectAny("Hit integral layer=1");
69  if( !Hit_integral_layer1) { cout << "Can't find Hit integral layer=1!" << endl; return; }
70  TH1I *Hit_integral_layer2 = (TH1I*)gDirectory->FindObjectAny("Hit integral layer=2");
71  if( !Hit_integral_layer2) { cout << "Can't find Hit integral layer=2!" << endl; return; }
72  TH1I *Hit_integral_layer3 = (TH1I*)gDirectory->FindObjectAny("Hit integral layer=3");
73  if( !Hit_integral_layer3) { cout << "Can't find Hit integral layer=3!" << endl; return; }
74  TH1I *Hit_integral_layer4 = (TH1I*)gDirectory->FindObjectAny("Hit integral layer=4");
75  if( !Hit_integral_layer4) { cout << "Can't find Hit integral layer=4!" << endl; return; }
76 
77  TString histfile_name2= "hd_root_particle_gun_Photons_gensat_correct_030736_000";
78  TFile *histfile_nocorr = new TFile(histfile_name2+".root","read");
79  TH1I *Eshower_nocorr = (TH1I*)gDirectory->FindObjectAny("Eshower");
80 
81  // define derivative histograms
82 
83  int nbins=100;
84  TH1D *Eshower_diff = new TH1D("Eshower_diff","Eshower (Corr - Uncorr)",10*nbins,0,10);
85  Eshower_diff->Sumw2();
86  Eshower_diff->Add(Eshower,Eshower_nocorr,1.,-1.);
87  TH1D *Eshower_ave = new TH1D("Eshower_ave","Eshower (Corr+Uncorr)/2",10*nbins,0,10);
88  Eshower_ave->Add(Eshower,Eshower_nocorr,0.5,0.5);
89  Eshower_ave->Sumw2();
90  TH1D *Eshower_diff_norm = new TH1D("Eshower_diff_norm","Eshower Diff/Ave",10*nbins,0,10);
91  Eshower_diff_norm->Divide(Eshower_diff,Eshower_ave);
92  Eshower_diff_norm->Sumw2();
93 
94  //
95  TCanvas *c1 = new TCanvas("c1","c1 BCAL_SiPM_saturation",200,10,700,700);
96  c1->Divide(2,2);
97 
98  c1->cd(1);
99  Ethrown->SetTitle(title);
100  // Ethrown->GetXaxis()->SetRangeUser(xmin,xmax);
101  // Ethrown->GetYaxis()->SetRangeUser(ymin,ymax);
102  Ethrown->GetXaxis()->SetTitleSize(0.05);
103  Ethrown->GetYaxis()->SetTitleSize(0.05);
104  Ethrown->GetYaxis()->SetTitle("Events");
105  Ethrown->GetXaxis()->SetTitle("Ethrown");
106  Ethrown->SetLineColor(2);
107  Ethrown->Draw("hist");
108 
109  c1->cd(2);
110  Eshower->SetTitle(title);
111  // Eshower->GetXaxis()->SetRangeUser(xmin,xmax);
112  // Eshower->GetYaxis()->SetRangeUser(ymin,ymax);
113  Eshower->GetXaxis()->SetTitleSize(0.05);
114  Eshower->GetYaxis()->SetTitleSize(0.05);
115  Eshower->GetYaxis()->SetTitle("Events");
116  Eshower->GetXaxis()->SetTitle("Eshower");
117  Eshower->SetLineColor(2);
118  Eshower->Draw("hist");
119  // Eshower_nocorr->Draw("histsame");
120  // Eshower_nocorr->SetLineColor(4);
121 
122  c1->cd(3);
123  ThrownTheta->SetTitle(title);
124  // ThrownTheta->GetXaxis()->SetRangeUser(xmin,xmax);
125  // ThrownTheta->GetYaxis()->SetRangeUser(ymin,ymax);
126  ThrownTheta->GetXaxis()->SetTitleSize(0.05);
127  ThrownTheta->GetYaxis()->SetTitleSize(0.05);
128  ThrownTheta->GetYaxis()->SetTitle("Events");
129  ThrownTheta->GetXaxis()->SetTitle("ThrownTheta");
130  ThrownTheta->SetLineColor(2);
131  ThrownTheta->Draw("hist");
132 
133  c1->cd(4);
134  layer->SetTitle(title);
135  // layer->GetXaxis()->SetRangeUser(xmin,xmax);
136  // layer->GetYaxis()->SetRangeUser(ymin,ymax);
137  layer->GetXaxis()->SetTitleSize(0.05);
138  layer->GetYaxis()->SetTitleSize(0.05);
139  layer->GetYaxis()->SetTitle("Events");
140  layer->GetXaxis()->SetTitle("layer");
141  layer->SetLineColor(2);
142  layer->Draw("hist");
143 
144 
145  TCanvas *c2 = new TCanvas("c2","c2 BCAL_SiPM_saturation",200,10,700,700);
146  c2->Divide(2,2);
147 
148  c2->cd(1);
149  NCell->SetTitle("");
150  // NCell->GetXaxis()->SetRangeUser(xmin,xmax);
151  // NCell->GetYaxis()->SetRangeUser(ymin,ymax);
152  NCell->GetXaxis()->SetTitleSize(0.05);
153  NCell->GetYaxis()->SetTitleSize(0.05);
154  NCell->GetYaxis()->SetTitle("Events");
155  NCell->GetXaxis()->SetTitle("NCell");
156  NCell->SetLineColor(2);
157  NCell->Draw("hist");
158 
159  c2->cd(2);
160  gPad->SetLogy();
161  Ept->SetTitle("");
162  // Ept->GetXaxis()->SetRangeUser(xmin,xmax);
163  // Ept->GetYaxis()->SetRangeUser(ymin,ymax);
164  Ept->GetXaxis()->SetTitleSize(0.05);
165  Ept->GetYaxis()->SetTitleSize(0.05);
166  Ept->GetYaxis()->SetTitle("Events");
167  Ept->GetXaxis()->SetTitle("Ept");
168  Ept->SetLineColor(2);
169  Ept->Draw("hist");
170 
171  c2->cd(3);
172  Eshower_nocorr->SetTitle("");
173  // Eshower_nocorr->GetXaxis()->SetRangeUser(xmin,xmax);
174  // Eshower_nocorr->GetYaxis()->SetRangeUser(ymin,ymax);
175  Eshower_nocorr->GetXaxis()->SetTitleSize(0.05);
176  Eshower_nocorr->GetYaxis()->SetTitleSize(0.05);
177  Eshower_nocorr->GetYaxis()->SetTitle("Events");
178  Eshower_nocorr->GetXaxis()->SetTitle("Eshower NO correction ");
179  Eshower_nocorr->SetLineColor(2);
180  Eshower_nocorr->Draw("hist");
181 
182  c2->cd(4);
183  gPad->SetLogy();
184  EcalcDiff->SetTitle("");
185  // EcalcDiff->GetXaxis()->SetRangeUser(xmin,xmax);
186  // EcalcDiff->GetYaxis()->SetRangeUser(ymin,ymax);
187  EcalcDiff->GetXaxis()->SetTitleSize(0.05);
188  EcalcDiff->GetYaxis()->SetTitleSize(0.05);
189  EcalcDiff->GetYaxis()->SetTitle("Events");
190  EcalcDiff->GetXaxis()->SetTitle("EcalcDiff");
191  EcalcDiff->SetLineColor(2);
192  EcalcDiff->SetNdivisions(505);
193  EcalcDiff->Draw("hist");
194 
195 
196  TCanvas *c3 = new TCanvas("c3","c3 BCAL_SiPM_saturation",200,10,700,700);
197  c3->Divide(2,2);
198 
199  c3->cd(1);
200  Eshower_vs_Ethrown->SetTitle(title);
201  // Eshower_vs_Ethrown->GetXaxis()->SetRangeUser(xmin,xmax);
202  // Eshower_vs_Ethrown->GetYaxis()->SetRangeUser(ymin,ymax);
203  Eshower_vs_Ethrown->GetXaxis()->SetTitleSize(0.05);
204  Eshower_vs_Ethrown->GetYaxis()->SetTitleSize(0.05);
205  Eshower_vs_Ethrown->GetYaxis()->SetTitle("Eshower");
206  Eshower_vs_Ethrown->GetXaxis()->SetTitle("Ethrown");
207  Eshower_vs_Ethrown->SetLineColor(2);
208  Eshower_vs_Ethrown->Draw("colz");
209 
210  c3->cd(2);
211  // gPad->SetLogy();
212  EDiff_vs_Ethrown->SetTitle(title);
213  // EDiff_vs_Ethrown->GetXaxis()->SetRangeUser(xmin,xmax);
214  // EDiff_vs_Ethrown->GetYaxis()->SetRangeUser(ymin,ymax);
215  EDiff_vs_Ethrown->GetXaxis()->SetTitleSize(0.05);
216  EDiff_vs_Ethrown->GetYaxis()->SetTitleSize(0.05);
217  EDiff_vs_Ethrown->GetYaxis()->SetTitle("EDiff");
218  EDiff_vs_Ethrown->GetXaxis()->SetTitle("Ethrown");
219  EDiff_vs_Ethrown->SetLineColor(2);
220  EDiff_vs_Ethrown->Draw("colz");
221 
222  c3->cd(3);
223  // gPad->SetLogy();
224  EDiffoverEthrown_vs_Ethrown->SetTitle(title);
225  // EDiff overEThrown_vs_Ethrown->GetXaxis()->SetRangeUser(xmin,xmax);
226  // EDiff overEThrown_vs_Ethrown->GetYaxis()->SetRangeUser(ymin,ymax);
227  EDiffoverEthrown_vs_Ethrown->GetXaxis()->SetTitleSize(0.05);
228  EDiffoverEthrown_vs_Ethrown->GetYaxis()->SetTitleSize(0.05);
229  EDiffoverEthrown_vs_Ethrown->GetXaxis()->SetTitle("Events");
230  EDiffoverEthrown_vs_Ethrown->GetYaxis()->SetTitle("EDiff/EThrown_vs_Ethrown");
231  EDiffoverEthrown_vs_Ethrown->SetLineColor(2);
232  EDiffoverEthrown_vs_Ethrown->Draw("colz");
233 
234 
235  TCanvas *c4 = new TCanvas("c4","c4 BCAL_SiPM_saturation",200,10,700,700);
236  c4->Divide(2,2);
237 
238  c4->cd(1);
239  gPad->SetLogy();
240  Hit_pulse_peak_layer1->SetTitle("");
241  // Hit_pulse_peak_layer1->GetXaxis()->SetRangeUser(xmin,xmax);
242  // Hit_pulse_peak_layer1->GetYaxis()->SetRangeUser(ymin,ymax);
243  Hit_pulse_peak_layer1->GetXaxis()->SetTitleSize(0.05);
244  Hit_pulse_peak_layer1->GetYaxis()->SetTitleSize(0.05);
245  Hit_pulse_peak_layer1->GetYaxis()->SetTitle("Events");
246  Hit_pulse_peak_layer1->GetXaxis()->SetTitle("Hit Peak, Layer 1(FADC counts)");
247  Hit_pulse_peak_layer1->SetLineColor(2);
248  Hit_pulse_peak_layer1->Draw("");
249 
250  c4->cd(2);
251  gPad->SetLogy();
252  Hit_pulse_peak_layer2->SetTitle("");
253  // Hit_pulse_peak_layer2->GetXaxis()->SetRangeUser(xmin,xmax);
254  // Hit_pulse_peak_layer2->GetYaxis()->SetRangeUser(ymin,ymax);
255  Hit_pulse_peak_layer2->GetXaxis()->SetTitleSize(0.05);
256  Hit_pulse_peak_layer2->GetYaxis()->SetTitleSize(0.05);
257  Hit_pulse_peak_layer2->GetYaxis()->SetTitle("Events");
258  Hit_pulse_peak_layer2->GetXaxis()->SetTitle("Hit Peak, Layer 2(FADC counts)");
259  Hit_pulse_peak_layer2->SetLineColor(2);
260  Hit_pulse_peak_layer2->Draw("");
261 
262  c4->cd(3);
263  gPad->SetLogy();
264  Hit_pulse_peak_layer3->SetTitle("");
265  // Hit_pulse_peak_layer3->GetXaxis()->SetRangeUser(xmin,xmax);
266  // Hit_pulse_peak_layer3->GetYaxis()->SetRangeUser(ymin,ymax);
267  Hit_pulse_peak_layer3->GetXaxis()->SetTitleSize(0.05);
268  Hit_pulse_peak_layer3->GetYaxis()->SetTitleSize(0.05);
269  Hit_pulse_peak_layer3->GetYaxis()->SetTitle("Events");
270  Hit_pulse_peak_layer3->GetXaxis()->SetTitle("Hit Peak, Layer 3(FADC counts)");
271  Hit_pulse_peak_layer3->SetLineColor(2);
272  Hit_pulse_peak_layer3->Draw("");
273 
274  c4->cd(4);
275  gPad->SetLogy();
276  Hit_pulse_peak_layer4->SetTitle("");
277  // Hit_pulse_peak_layer4->GetXaxis()->SetRangeUser(xmin,xmax);
278  // Hit_pulse_peak_layer4->GetYaxis()->SetRangeUser(ymin,ymax);
279  Hit_pulse_peak_layer4->GetXaxis()->SetTitleSize(0.05);
280  Hit_pulse_peak_layer4->GetYaxis()->SetTitleSize(0.05);
281  Hit_pulse_peak_layer4->GetYaxis()->SetTitle("Events");
282  Hit_pulse_peak_layer4->GetXaxis()->SetTitle("Hit Peak, Layer 4(FADC counts)");
283  Hit_pulse_peak_layer4->SetLineColor(2);
284  Hit_pulse_peak_layer4->Draw("");
285 
286 
287  TCanvas *c5 = new TCanvas("c5","c5 BCAL_SiPM_saturation",200,10,700,700);
288  c5->Divide(2,2);
289 
290  c5->cd(1);
291  gPad->SetLogy();
292  Hit_integral_layer1->SetTitle("");
293  // Hit_integral_layer1->GetXaxis()->SetRangeUser(xmin,xmax);
294  // Hit_integral_layer1->GetYaxis()->SetRangeUser(ymin,ymax);
295  Hit_integral_layer1->GetXaxis()->SetTitleSize(0.05);
296  Hit_integral_layer1->GetYaxis()->SetTitleSize(0.05);
297  Hit_integral_layer1->GetYaxis()->SetTitle("Events");
298  Hit_integral_layer1->GetXaxis()->SetTitle("Hit Integral, Layer 1(GeV)");
299  Hit_integral_layer1->SetLineColor(2);
300  Hit_integral_layer1->Draw("");
301 
302  c5->cd(2);
303  gPad->SetLogy();
304  Hit_integral_layer2->SetTitle("");
305  // Hit_integral_layer2->GetXaxis()->SetRangeUser(xmin,xmax);
306  // Hit_integral_layer2->GetYaxis()->SetRangeUser(ymin,ymax);
307  Hit_integral_layer2->GetXaxis()->SetTitleSize(0.05);
308  Hit_integral_layer2->GetYaxis()->SetTitleSize(0.05);
309  Hit_integral_layer2->GetYaxis()->SetTitle("Events");
310  Hit_integral_layer2->GetXaxis()->SetTitle("Hit Integral, Layer 2(GeV)");
311  Hit_integral_layer2->SetLineColor(2);
312  Hit_integral_layer2->Draw("");
313 
314  c5->cd(3);
315  gPad->SetLogy();
316  Hit_integral_layer3->SetTitle("");
317  // Hit_integral_layer3->GetXaxis()->SetRangeUser(xmin,xmax);
318  // Hit_integral_layer3->GetYaxis()->SetRangeUser(ymin,ymax);
319  Hit_integral_layer3->GetXaxis()->SetTitleSize(0.05);
320  Hit_integral_layer3->GetYaxis()->SetTitleSize(0.05);
321  Hit_integral_layer3->GetYaxis()->SetTitle("Events");
322  Hit_integral_layer3->GetXaxis()->SetTitle("Hit Integral, Layer 3(GeV)");
323  Hit_integral_layer3->SetLineColor(2);
324  Hit_integral_layer3->Draw("");
325 
326  c5->cd(4);
327  gPad->SetLogy();
328  Hit_integral_layer4->SetTitle("");
329  // Hit_integral_layer4->GetXaxis()->SetRangeUser(xmin,xmax);
330  // Hit_integral_layer4->GetYaxis()->SetRangeUser(ymin,ymax);
331  Hit_integral_layer4->GetXaxis()->SetTitleSize(0.05);
332  Hit_integral_layer4->GetYaxis()->SetTitleSize(0.05);
333  Hit_integral_layer4->GetYaxis()->SetTitle("Events");
334  Hit_integral_layer4->GetXaxis()->SetTitle("Hit Integral, Layer 4(GeV)");
335  Hit_integral_layer4->SetLineColor(2);
336  Hit_integral_layer4->Draw(""); //
337 
338  TCanvas *c6 = new TCanvas("c6","c6 BCAL_SiPM_saturation",200,10,700,700);
339  // c6->Divide(2,2);
340  // c6->cd(1);
341  Eshower->SetTitle("");
342  // Eshower->GetXaxis()->SetRangeUser(xmin,xmax);
343  // Eshower->GetYaxis()->SetRangeUser(ymin,ymax);
344  Eshower->GetXaxis()->SetTitleSize(0.05);
345  Eshower->GetYaxis()->SetTitleSize(0.05);
346  Eshower->GetYaxis()->SetTitle("Events");
347  Eshower->GetXaxis()->SetTitle("Eshower");
348  Eshower->SetLineColor(2);
349  Eshower->Draw("hist");
350  Eshower_nocorr->Draw("histsame");
351  Eshower_nocorr->SetLineColor(4);
352 
353  TCanvas *c7 = new TCanvas("c7","c7 BCAL_SiPM_saturation",200,10,700,700);
354  c7->Divide(2,2);
355 
356  c7->cd(1);
357  gPad->SetLogy();
358  Eshower_ave->SetTitle("");
359  // Eshower_ave->GetXaxis()->SetRangeUser(xmin,xmax);
360  // Eshower_ave->GetYaxis()->SetRangeUser(ymin,ymax);
361  Eshower_ave->GetXaxis()->SetTitleSize(0.05);
362  Eshower_ave->GetYaxis()->SetTitleSize(0.05);
363  Eshower_ave->GetYaxis()->SetTitle("Events");
364  Eshower_ave->GetXaxis()->SetTitle("Eshower Ave (Corr+UnCorr)/2");
365  Eshower_ave->SetLineColor(2);
366  Eshower_ave->RebinX(40);
367  Eshower_ave->Draw("hist");
368 
369  c7->cd(2);
370  Eshower_diff->SetTitle("");
371  // Eshower_diff->GetXaxis()->SetRangeUser(xmin,xmax);
372  // Eshower_diff->GetYaxis()->SetRangeUser(ymin,ymax);
373  Eshower_diff->GetXaxis()->SetTitleSize(0.05);
374  Eshower_diff->GetYaxis()->SetTitleSize(0.05);
375  Eshower_diff->GetYaxis()->SetTitle("Events");
376  Eshower_diff->GetXaxis()->SetTitle("Eshower Diff (Corr-UnCorr)");
377  Eshower_diff->SetLineColor(2);
378  Eshower_diff->RebinX(40);
379  Eshower_diff->Draw("hist");
380 
381  c7->cd(3);
382  Eshower_diff_norm->SetTitle("");
383  // Eshower_diff_norm->GetXaxis()->SetRangeUser(xmin,xmax);
384  // Eshower_diff_norm->GetYaxis()->SetRangeUser(ymin,ymax);
385  Eshower_diff_norm->GetXaxis()->SetTitleSize(0.05);
386  Eshower_diff_norm->GetYaxis()->SetTitleSize(0.05);
387  Eshower_diff_norm->GetYaxis()->SetTitle("Fraction");
388  Eshower_diff_norm->GetXaxis()->SetTitle("Eshower Diff/Ave");
389  Eshower_diff_norm->SetLineColor(2);
390  Eshower_diff_norm->RebinX(40);
391  Eshower_diff_norm->Draw("hist");
392 
393 
394 
395 
396  //
397  TCanvas *c8 = new TCanvas("c8","c8 BCAL_SiPM_saturation",200,10,700,700);
398  c8->Divide(2,2);
399 
400 
401 
402  c8->cd(1);
403  gPad->SetLogy();
404  Eshower->SetTitle(title);
405  // Eshower->GetXaxis()->SetRangeUser(xmin,xmax);
406  // Eshower->GetYaxis()->SetRangeUser(ymin,ymax);
407  Eshower->GetXaxis()->SetTitleSize(0.05);
408  Eshower->GetYaxis()->SetTitleSize(0.05);
409  Eshower->GetYaxis()->SetTitle("Events");
410  Eshower->GetXaxis()->SetTitle("Eshower (GeV)");
411  Eshower->SetLineColor(2);
412  Eshower->Draw("hist");
413  // Eshower_nocorr->Draw("histsame");
414  // Eshower_nocorr->SetLineColor(4);
415 
416 
417  c8->cd(2);
418  NCell->SetTitle("");
419  // NCell->GetXaxis()->SetRangeUser(xmin,xmax);
420  // NCell->GetYaxis()->SetRangeUser(ymin,ymax);
421  NCell->GetXaxis()->SetTitleSize(0.05);
422  NCell->GetYaxis()->SetTitleSize(0.05);
423  NCell->GetYaxis()->SetTitle("Events");
424  NCell->GetXaxis()->SetTitle("NCell");
425  NCell->SetLineColor(2);
426  NCell->Draw("hist");
427 
428  c8->cd(3);
429  gPad->SetLogy();
430  Ept->SetTitle("");
431  // Ept->GetXaxis()->SetRangeUser(xmin,xmax);
432  // Ept->GetYaxis()->SetRangeUser(ymin,ymax);
433  Ept->GetXaxis()->SetTitleSize(0.05);
434  Ept->GetYaxis()->SetTitleSize(0.05);
435  Ept->GetYaxis()->SetTitle("Events");
436  Ept->GetXaxis()->SetTitle("Ept (GeV)");
437  Ept->SetLineColor(2);
438  Ept->Draw("hist");
439 
440 
441  c8->cd(4);
442  layer->SetTitle(title);
443  // layer->GetXaxis()->SetRangeUser(xmin,xmax);
444  // layer->GetYaxis()->SetRangeUser(ymin,ymax);
445  layer->GetXaxis()->SetTitleSize(0.05);
446  layer->GetYaxis()->SetTitleSize(0.05);
447  layer->GetYaxis()->SetTitle("Events");
448  layer->GetXaxis()->SetTitle("layer");
449  layer->SetLineColor(2);
450  layer->Draw("hist");
451 
452 
453 
454 
455 
456  /*TLatex *t1 = new TLatex(0.1,0.92,"");
457  t1->SetNDC();
458  t1->DrawLatex(0.2,0.82,"b)");
459  t1->DrawLatex(0.2,0.75,"Four Configurations");
460  t1->SetTextSize(0.07);*/
461 
462 
463  c1->SaveAs(histfile_name+".pdf(");
464  c2->SaveAs(histfile_name+".pdf");
465  c3->SaveAs(histfile_name+".pdf");
466  c4->SaveAs(histfile_name+".pdf");
467  c5->SaveAs(histfile_name+".pdf");
468  c6->SaveAs(histfile_name+".pdf");
469  c7->SaveAs(histfile_name+".pdf");
470  c8->SaveAs(histfile_name+".pdf)");
471 
472  // in->Close();
473 
474 }
475 
layer
Int_t layer
Definition: Read_bcal_hadronic_eff2.C:50
c1
Double_t c1[2][NMODULES]
Definition: tw_corr.C:68
BCAL_SiPM_saturation
void BCAL_SiPM_saturation(void)
Definition: BCAL_SiPM_saturation.C:1
c2
Double_t c2[2][NMODULES]
Definition: tw_corr.C:69
c6
TCanvas * c6
Definition: ST_Monitoring_Expert_Waveform.C:220
ymin
Double_t ymin
Definition: bcal_hist_eff.C:89
c4
TCanvas * c4
Definition: plot_epics_dump.C:170
c5
TCanvas * c5
Definition: ST_Monitoring_Expert_Waveform.C:185
ymax
Double_t ymax
Definition: bcal_hist_eff.C:91
c3
TCanvas * c3
Definition: plot_epics_dump.C:151