Hall-D Software  alpha
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
trig_fcalbcal2.C
Go to the documentation of this file.
1 
2 // File: trig_fcalbcal.C
3 // Created: 05/01/15
4 // Creator: Elton Smith
5 // Purpose: Display basic energy plots in fcal and bcal for trigger
6 // Uses output of TRIG_online/JEventProcessor_TRIG_online.cc
7 
8 // hnamepath: /trig/h1trig_fcal
9 // hnamepath: /trig/h1trig_fcalN
10 // hnamepath: /trig/h1trig_bcal
11 // hnamepath: /trig/h1trig_bcalN
12 // hnamepath: /trig/h2trig_fcalVSbcal
13 
14 {
15  gStyle->SetPalette(1,0);
16  gStyle->SetOptStat(kTRUE);
17  gStyle->SetPadRightMargin(0.15);
18  gStyle->SetPadLeftMargin(0.15);
19  gStyle->SetPadBottomMargin(0.15);
20 
21  char string[256];
22 
23  TString filename = "hd_rawdata_010270_000";
24  TFile* f = new TFile(filename+".root");
25 
26  // get histograms from trig subdirectory
27 
28  TDirectory *dir = (TDirectory*)gDirectory->FindObjectAny("trig");
29  if(!dir){
30  cout << "Can't find trig TDirectory!" << endl;
31  return;
32  }
33 
34  TH1F *h1trig_fcal = (TH1F*)dir->Get("h1trig_fcal" );
35  TH1F *h1trig_fcalN = (TH1F*)dir->Get("h1trig_fcalN" );
36  TH1F *h1trig_bcal = (TH1F*)dir->Get("h1trig_bcal" );
37  TH1F *h1trig_bcalN = (TH1F*)dir->Get("h1trig_bcalN" );
38  TH1F *h2trig_fcalVSbcal = (TH1F*)dir->Get("h2trig_fcalVSbcal" );
39 
40  if( !h1trig_fcal ) { cout << "Can't find h1trig_fcal!" << endl; return; }
41  if( !h1trig_fcalN ) { cout << "Can't find h1trig_fcalN!" << endl; return; }
42  if( !h1trig_bcal ) { cout << "Can't find h1trig_bcal!" << endl; return; }
43  if( !h1trig_bcalN ) { cout << "Can't find h1trig_bcalN!" << endl; return; }
44  if( !h2trig_fcalVSbcal ) { cout << "Can't find h2trig_fcalVSbcal!" << endl; return; }
45 
46  TH1F *h1trig_trgbits = (TH1F*)dir->Get("h1trig_trgbits" );
47 
48  TH1F *h1trig1_fcal = (TH1F*)dir->Get("h1trig1_fcal" );
49  TH1F *h1trig1_fcalN = (TH1F*)dir->Get("h1trig1_fcalN" );
50  TH1F *h1trig1_bcal = (TH1F*)dir->Get("h1trig1_bcal" );
51  TH1F *h1trig1_bcalN = (TH1F*)dir->Get("h1trig1_bcalN" );
52  TH1F *h2trig1_fcalVSbcal = (TH1F*)dir->Get("h2trig1_fcalVSbcal" );
53 
54  TH1F *h1trig3_fcal = (TH1F*)dir->Get("h1trig3_fcal" );
55  TH1F *h1trig3_fcalN = (TH1F*)dir->Get("h1trig3_fcalN" );
56  TH1F *h1trig3_bcal = (TH1F*)dir->Get("h1trig3_bcal" );
57  TH1F *h1trig3_bcalN = (TH1F*)dir->Get("h1trig3_bcalN" );
58  TH1F *h2trig3_fcalVSbcal = (TH1F*)dir->Get("h2trig3_fcalVSbcal" );
59 
60  TH1F *h1trig5_fcal = (TH1F*)dir->Get("h1trig5_fcal" );
61  TH1F *h1trig5_fcalN = (TH1F*)dir->Get("h1trig5_fcalN" );
62  TH1F *h1trig5_bcal = (TH1F*)dir->Get("h1trig5_bcal" );
63  TH1F *h1trig5_bcalN = (TH1F*)dir->Get("h1trig5_bcalN" );
64  TH1F *h2trig5_fcalVSbcal = (TH1F*)dir->Get("h2trig5_fcalVSbcal" );
65 
66  TH1F *h1trig7_fcal = (TH1F*)dir->Get("h1trig7_fcal" );
67  TH1F *h1trig7_fcalN = (TH1F*)dir->Get("h1trig7_fcalN" );
68  TH1F *h1trig7_bcal = (TH1F*)dir->Get("h1trig7_bcal" );
69  TH1F *h1trig7_bcalN = (TH1F*)dir->Get("h1trig7_bcalN" );
70  TH1F *h2trig7_fcalVSbcal = (TH1F*)dir->Get("h2trig7_fcalVSbcal" );
71 
72  double fcal_max = 2;
73  double bcal_max = 1;
74  double thf = 0.7;
75  double thb = 0.2;
76 
77  TCanvas *c0 = new TCanvas("c0", "c0",200,10,700,700);
78 
79  c0->Divide(2,2);
80  c0->cd(1);
81 
82  h1trig_trgbits->SetTitle(filename);
83  // h1trig_trgbits->GetXaxis()->SetRangeUser(xmin,xmax);
84  // h1trig_trgbits->GetYaxis()->SetRangeUser(ymin,ymax);
85  h1trig_trgbits->GetXaxis()->SetTitleSize(0.05);
86  h1trig_trgbits->GetYaxis()->SetTitleSize(0.05);
87  h1trig_trgbits->GetXaxis()->SetTitle("trig_mask || (10+fp_trig_mask)");
88  h1trig_trgbits->SetLineColor(2);
89  h1trig_trgbits->Draw("");
90 
91  c0->cd(2);
92  gPad->SetLogy();
93 
94  // h1trig_trgbits->SetTitle("");
95  // h1trig_trgbits->GetXaxis()->SetRangeUser(xmin,xmax);
96  // h1trig_trgbits->GetYaxis()->SetRangeUser(ymin,ymax);
97  h1trig_trgbits->GetXaxis()->SetTitleSize(0.05);
98  h1trig_trgbits->GetYaxis()->SetTitleSize(0.05);
99  h1trig_trgbits->GetXaxis()->SetTitle("trig_mask || (10+fp_trig_mask)");
100  h1trig_trgbits->SetLineColor(2);
101  h1trig_trgbits->Draw("");
102 
103  c0->cd(3);
104 
105 
106  TH2F *temp = new TH2F("temp","Template",10,0,1,10,0,2);
107  temp->GetXaxis()->SetTitleSize(0.05);
108  temp->GetYaxis()->SetTitleSize(0.05);
109  temp->GetXaxis()->SetTitle("Bcal Energy (GeV)");
110  temp->GetYaxis()->SetTitle("Fcal Energy (GeV)");
111  temp->Draw();
112 
113  TLine *thresh_fb = new TLine(0,0.5,0.14,0);
114  thresh_fb->SetLineWidth(2);
115  thresh_fb->Draw();
116 
117  TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
118  thresh_f->SetLineWidth(2);
119  thresh_f->SetLineStyle(2);
120  thresh_f->Draw();
121 
122  TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
123  thresh_b->SetLineWidth(2);
124  thresh_b->SetLineStyle(2);
125  thresh_b->Draw();
126 
127 
128  sprintf (string,"trg 5 = 101\n");
129  printf("string=%s",string);
130  TLatex *t1 = new TLatex(0.5,0.4,string);
131  // t1->SetNDC();
132  t1->SetTextSize(0.04);
133  t1->Draw();
134 
135  sprintf (string,"trg 7 = 111\n");
136  t1->DrawLatex(0.5,1.4,string);
137 
138  sprintf (string,"trg 3=\n");
139  t1->DrawLatex(0.02,1.4,string);
140 
141  sprintf (string,"011\n");
142  t1->DrawLatex(0.01,1.2,string);
143 
144  sprintf (string,"trg 1=\n");
145  t1->DrawLatex(0.01,0.6,string);
146 
147  sprintf (string,"001\n");
148  t1->DrawLatex(0.1,0.4,string);
149 
150  c0->cd(4);
151 
152 
153  sprintf (string,"Trigger configuration\n");
154  t1->DrawLatex(0.1,0.9,string);
155  sprintf (string,"001 : 10Fcal +3Bcal > 18000 counts\n");
156  t1->DrawLatex(0.1,0.8,string);
157  sprintf (string,"\t\t: Fcal + 3.6 Bcal > 0.49 GeV \n");
158  t1->DrawLatex(0.1,0.75,string);
159  sprintf (string,"010 : Fcal > 1800 counts : Fcal > 0.49 GeV \n");
160  t1->DrawLatex(0.1,0.65,string);
161  sprintf (string,"100 : Bcal > 6000 counts: Bcal > 0.14 GeV\n");
162  t1->DrawLatex(0.1,0.6,string);
163 
164 
165  TCanvas *c1 = new TCanvas("c1", "c1",200,10,700,700);
166 
167  c1->Divide(2,2);
168  c1->cd(1);
169 
170  // h1trig_fcal->SetTitle("");
171  // h1trig_fcal->GetXaxis()->SetRangeUser(xmin,xmax);
172  // h1trig_fcal->GetYaxis()->SetRangeUser(ymin,ymax);
173  h1trig_fcal->GetXaxis()->SetTitleSize(0.05);
174  h1trig_fcal->GetYaxis()->SetTitleSize(0.05);
175  h1trig_fcal->GetXaxis()->SetTitle("Fcal Energy (GeV)");
176  h1trig_fcal->SetLineColor(2);
177  h1trig_fcal->Draw("");
178 
179  double ymax;
180  ymax = h1trig_fcal->GetMaximum();
181  TLine *thrf = new TLine(thf,0,thf,ymax);
182  thrf->SetLineWidth(2);
183  thrf->SetLineStyle(2);
184  thrf->DrawLine(thf,0,thf,ymax);
185 
186 
187  c1->cd(2);
188 
189  temp->SetTitle("");
190  // h2trig_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
191  // h2trig_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
192  h2trig_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
193  h2trig_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
194  h2trig_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
195  h2trig_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
196  h2trig_fcalVSbcal->SetLineColor(2);
197  h2trig_fcalVSbcal->Draw("colz");
198 
199  thresh_fb->Draw();
200  thresh_f->Draw();
201  thresh_b->Draw();
202 
203  c1->cd(3);
204 
205  h1trig_bcal->SetTitle("");
206  // h1trig_bcal->GetXaxis()->SetRangeUser(xmin,xmax);
207  // h1trig_bcal->GetYaxis()->SetRangeUser(ymin,ymax);
208  h1trig_bcal->GetXaxis()->SetTitleSize(0.05);
209  h1trig_bcal->GetYaxis()->SetTitleSize(0.05);
210  h1trig_bcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
211  h1trig_bcal->SetLineColor(2);
212  h1trig_bcal->Draw("");
213 
214  ymax = h1trig_bcal->GetMaximum();
215  TLine *thrb = new TLine(thb,0,thb,ymax);
216  thrb->SetLineWidth(2);
217  thrb->SetLineStyle(2);
218  thrb->DrawLine(thb,0,thb,ymax);
219 
220  c1->cd(4);
221  gPad->SetLogz();
222 
223  h2trig_fcalVSbcal->SetTitle("");
224  // h2trig_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
225  // h2trig_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
226  h2trig_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
227  h2trig_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
228  h2trig_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
229  h2trig_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
230  h2trig_fcalVSbcal->SetLineColor(2);
231  h2trig_fcalVSbcal->Draw("colz");
232 
233  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
234  thresh_fb->SetLineWidth(2);
235  thresh_fb->Draw();
236 
237  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
238  thresh_f->SetLineWidth(2);
239  thresh_f->SetLineStyle(2);
240  thresh_f->Draw();
241 
242  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
243  thresh_b->SetLineWidth(2);
244  thresh_b->SetLineStyle(2);
245  thresh_b->Draw();
246 
247  TCanvas *c2 = new TCanvas("c2", "c2",200,10,700,700);
248 
249  c2->Divide(2,2);
250  c2->cd(1);
251 
252  // h1trig1_fcal->SetTitle("");
253  // h1trig1_fcal->GetXaxis()->SetRangeUser(xmin,xmax);
254  // h1trig1_fcal->GetYaxis()->SetRangeUser(ymin,ymax);
255  h1trig1_fcal->GetXaxis()->SetTitleSize(0.05);
256  h1trig1_fcal->GetYaxis()->SetTitleSize(0.05);
257  h1trig1_fcal->GetXaxis()->SetTitle("Fcal Energy (GeV)");
258  h1trig1_fcal->SetLineColor(2);
259  h1trig1_fcal->Draw("");
260 
261  ymax = h1trig1_fcal->GetMaximum();
262  thrf->DrawLine(thf,0,thf,ymax);
263 
264  c2->cd(2);
265 
266  h2trig1_fcalVSbcal->SetTitle("");
267  // h2trig1_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
268  // h2trig1_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
269  h2trig1_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
270  h2trig1_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
271  h2trig1_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
272  h2trig1_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
273  h2trig1_fcalVSbcal->SetLineColor(2);
274  h2trig1_fcalVSbcal->Draw("colz");
275 
276  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
277  thresh_fb->SetLineWidth(2);
278  thresh_fb->Draw();
279 
280  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
281  thresh_f->SetLineWidth(2);
282  thresh_f->SetLineStyle(2);
283  thresh_f->Draw();
284 
285  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
286  thresh_b->SetLineWidth(2);
287  thresh_b->SetLineStyle(2);
288  thresh_b->Draw();
289 
290  c2->cd(3);
291 
292  h1trig1_bcal->SetTitle("");
293  // h1trig1_bcal->GetXaxis()->SetRangeUser(xmin,xmax);
294  // h1trig1_bcal->GetYaxis()->SetRangeUser(ymin,ymax);
295  h1trig1_bcal->GetXaxis()->SetTitleSize(0.05);
296  h1trig1_bcal->GetYaxis()->SetTitleSize(0.05);
297  h1trig1_bcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
298  h1trig1_bcal->SetLineColor(2);
299  h1trig1_bcal->Draw("");
300 
301  ymax = h1trig1_bcal->GetMaximum();
302  thrb->DrawLine(thb,0,thb,ymax);
303 
304  c2->cd(4);
305  gPad->SetLogz();
306 
307  h2trig1_fcalVSbcal->SetTitle("");
308  // h2trig1_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
309  // h2trig1_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
310  h2trig1_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
311  h2trig1_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
312  h2trig1_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
313  h2trig1_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
314  h2trig1_fcalVSbcal->SetLineColor(2);
315  h2trig1_fcalVSbcal->Draw("colz");
316 
317  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
318  thresh_fb->SetLineWidth(2);
319  thresh_fb->Draw();
320 
321  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
322  thresh_f->SetLineWidth(2);
323  thresh_f->SetLineStyle(2);
324  thresh_f->Draw();
325 
326  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
327  thresh_b->SetLineWidth(2);
328  thresh_b->SetLineStyle(2);
329  thresh_b->Draw();
330 
331  TCanvas *c3 = new TCanvas("c3", "c3",200,10,700,700);
332 
333  c3->Divide(2,2);
334  c3->cd(1);
335 
336  // h1trig3_fcal->SetTitle("");
337  // h1trig3_fcal->GetXaxis()->SetRangeUser(xmin,xmax);
338  // h1trig3_fcal->GetYaxis()->SetRangeUser(ymin,ymax);
339  h1trig3_fcal->GetXaxis()->SetTitleSize(0.05);
340  h1trig3_fcal->GetYaxis()->SetTitleSize(0.05);
341  h1trig3_fcal->GetXaxis()->SetTitle("Fcal Energy (GeV)");
342  h1trig3_fcal->SetLineColor(2);
343  h1trig3_fcal->Draw("");
344 
345  ymax = h1trig3_fcal->GetMaximum();
346  thrf->DrawLine(thf,0,thf,ymax);
347 
348  c3->cd(2);
349 
350  h2trig3_fcalVSbcal->SetTitle("");
351  // h2trig3_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
352  // h2trig3_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
353  h2trig3_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
354  h2trig3_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
355  h2trig3_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
356  h2trig3_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
357  h2trig3_fcalVSbcal->SetLineColor(2);
358  h2trig3_fcalVSbcal->Draw("colz");
359 
360  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
361  thresh_fb->SetLineWidth(2);
362  thresh_fb->Draw();
363 
364  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
365  thresh_f->SetLineWidth(2);
366  thresh_f->SetLineStyle(2);
367  thresh_f->Draw();
368 
369  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
370  thresh_b->SetLineWidth(2);
371  thresh_b->SetLineStyle(2);
372  thresh_b->Draw();
373 
374  c3->cd(3);
375 
376  h1trig3_bcal->SetTitle("");
377  // h1trig3_bcal->GetXaxis()->SetRangeUser(xmin,xmax);
378  // h1trig3_bcal->GetYaxis()->SetRangeUser(ymin,ymax);
379  h1trig3_bcal->GetXaxis()->SetTitleSize(0.05);
380  h1trig3_bcal->GetYaxis()->SetTitleSize(0.05);
381  h1trig3_bcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
382  h1trig3_bcal->SetLineColor(2);
383  h1trig3_bcal->Draw("");
384 
385  ymax = h1trig3_bcal->GetMaximum();
386  thrb->DrawLine(thb,0,thb,ymax);
387 
388 
389  c3->cd(4);
390  gPad->SetLogz();
391 
392  h2trig3_fcalVSbcal->SetTitle("");
393  // h2trig3_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
394  // h2trig3_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
395  h2trig3_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
396  h2trig3_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
397  h2trig3_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
398  h2trig3_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
399  h2trig3_fcalVSbcal->SetLineColor(2);
400  h2trig3_fcalVSbcal->Draw("colz");
401 
402  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
403  thresh_fb->SetLineWidth(2);
404  thresh_fb->Draw();
405 
406  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
407  thresh_f->SetLineWidth(2);
408  thresh_f->SetLineStyle(2);
409  thresh_f->Draw();
410 
411  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
412  thresh_b->SetLineWidth(2);
413  thresh_b->SetLineStyle(2);
414  thresh_b->Draw();
415 
416  TCanvas *c4 = new TCanvas("c4", "c4",200,10,700,700);
417 
418  c4->Divide(2,2);
419  c4->cd(1);
420 
421  // h1trig5_fcal->SetTitle("");
422  // h1trig5_fcal->GetXaxis()->SetRangeUser(xmin,xmax);
423  // h1trig5_fcal->GetYaxis()->SetRangeUser(ymin,ymax);
424  h1trig5_fcal->GetXaxis()->SetTitleSize(0.05);
425  h1trig5_fcal->GetYaxis()->SetTitleSize(0.05);
426  h1trig5_fcal->GetXaxis()->SetTitle("Fcal Energy (GeV)");
427  h1trig5_fcal->SetLineColor(2);
428  h1trig5_fcal->Draw("");
429 
430  ymax = h1trig5_fcal->GetMaximum();
431  thrf->DrawLine(thf,0,thf,ymax);
432 
433  c4->cd(2);
434 
435  h2trig5_fcalVSbcal->SetTitle("");
436  // h2trig5_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
437  // h2trig5_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
438  h2trig5_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
439  h2trig5_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
440  h2trig5_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
441  h2trig5_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
442  h2trig5_fcalVSbcal->SetLineColor(2);
443  h2trig5_fcalVSbcal->Draw("colz");
444 
445  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
446  thresh_fb->SetLineWidth(2);
447  thresh_fb->Draw();
448 
449  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
450  thresh_f->SetLineWidth(2);
451  thresh_f->SetLineStyle(2);
452  thresh_f->Draw();
453 
454  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
455  thresh_b->SetLineWidth(2);
456  thresh_b->SetLineStyle(2);
457  thresh_b->Draw();
458 
459  c4->cd(3);
460 
461  h1trig5_bcal->SetTitle("");
462  // h1trig5_bcal->GetXaxis()->SetRangeUser(xmin,xmax);
463  // h1trig5_bcal->GetYaxis()->SetRangeUser(ymin,ymax);
464  h1trig5_bcal->GetXaxis()->SetTitleSize(0.05);
465  h1trig5_bcal->GetYaxis()->SetTitleSize(0.05);
466  h1trig5_bcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
467  h1trig5_bcal->SetLineColor(2);
468  h1trig5_bcal->Draw("");
469 
470  ymax = h1trig5_bcal->GetMaximum();
471  thrb->DrawLine(thb,0,thb,ymax);
472 
473 
474  c4->cd(4);
475  gPad->SetLogz();
476 
477  h2trig5_fcalVSbcal->SetTitle("");
478  // h2trig5_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
479  // h2trig5_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
480  h2trig5_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
481  h2trig5_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
482  h2trig5_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
483  h2trig5_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
484  h2trig5_fcalVSbcal->SetLineColor(2);
485  h2trig5_fcalVSbcal->Draw("colz");
486 
487  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
488  thresh_fb->SetLineWidth(2);
489  thresh_fb->Draw();
490 
491  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
492  thresh_f->SetLineWidth(2);
493  thresh_f->SetLineStyle(2);
494  thresh_f->Draw();
495 
496  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
497  thresh_b->SetLineWidth(2);
498  thresh_b->SetLineStyle(2);
499  thresh_b->Draw();
500 
501 
502 
503  TCanvas *c5 = new TCanvas("c5", "c5",200,10,700,700);
504 
505  c5->Divide(2,2);
506  c5->cd(1);
507 
508  // h1trig7_fcal->SetTitle("");
509  // h1trig7_fcal->GetXaxis()->SetRangeUser(xmin,xmax);
510  // h1trig7_fcal->GetYaxis()->SetRangeUser(ymin,ymax);
511  h1trig7_fcal->GetXaxis()->SetTitleSize(0.05);
512  h1trig7_fcal->GetYaxis()->SetTitleSize(0.05);
513  h1trig7_fcal->GetXaxis()->SetTitle("Fcal Energy (GeV)");
514  h1trig7_fcal->SetLineColor(2);
515  h1trig7_fcal->Draw("");
516 
517 
518  ymax = h1trig7_fcal->GetMaximum();
519  thrf->DrawLine(thf,0,thf,ymax);
520 
521  c5->cd(2);
522 
523  h2trig7_fcalVSbcal->SetTitle("");
524  // h2trig7_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
525  // h2trig7_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
526  h2trig7_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
527  h2trig7_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
528  h2trig7_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
529  h2trig7_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
530  h2trig7_fcalVSbcal->SetLineColor(2);
531  h2trig7_fcalVSbcal->Draw("colz");
532 
533  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
534  thresh_fb->SetLineWidth(2);
535  thresh_fb->Draw();
536 
537  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
538  thresh_f->SetLineWidth(2);
539  thresh_f->SetLineStyle(2);
540  thresh_f->Draw();
541 
542  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
543  thresh_b->SetLineWidth(2);
544  thresh_b->SetLineStyle(2);
545  thresh_b->Draw();
546 
547  c5->cd(3);
548 
549  h1trig7_bcal->SetTitle("");
550  // h1trig7_bcal->GetXaxis()->SetRangeUser(xmin,xmax);
551  // h1trig7_bcal->GetYaxis()->SetRangeUser(ymin,ymax);
552  h1trig7_bcal->GetXaxis()->SetTitleSize(0.05);
553  h1trig7_bcal->GetYaxis()->SetTitleSize(0.05);
554  h1trig7_bcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
555  h1trig7_bcal->SetLineColor(2);
556  h1trig7_bcal->Draw("");
557 
558  ymax = h1trig7_bcal->GetMaximum();
559  thrb->DrawLine(thb,0,thb,ymax);
560 
561 
562  c5->cd(4);
563  gPad->SetLogz();
564 
565  h2trig7_fcalVSbcal->SetTitle("");
566  // h2trig7_fcalVSbcal->GetXaxis()->SetRangeUser(xmin,xmax);
567  // h2trig7_fcalVSbcal->GetYaxis()->SetRangeUser(ymin,ymax);
568  h2trig7_fcalVSbcal->GetXaxis()->SetTitleSize(0.05);
569  h2trig7_fcalVSbcal->GetYaxis()->SetTitleSize(0.05);
570  h2trig7_fcalVSbcal->GetXaxis()->SetTitle("Bcal Energy (GeV)");
571  h2trig7_fcalVSbcal->GetYaxis()->SetTitle("Fcal Energy (GeV)");
572  h2trig7_fcalVSbcal->SetLineColor(2);
573  h2trig7_fcalVSbcal->Draw("colz");
574 
575  // TLine *thresh_fb = new TLine(0,0.5,0.14,0);
576  thresh_fb->SetLineWidth(2);
577  thresh_fb->Draw();
578 
579  // TLine *thresh_f = new TLine(0,thf,bcal_max,thf);
580  thresh_f->SetLineWidth(2);
581  thresh_f->SetLineStyle(2);
582  thresh_f->Draw();
583 
584  // TLine *thresh_b = new TLine(thb,0,thb,fcal_max);
585  thresh_b->SetLineWidth(2);
586  thresh_b->SetLineStyle(2);
587  thresh_b->Draw();
588 
589  sprintf (string,"trig_fcalbcal2.pdf(");
590  c0->SaveAs(string);
591  sprintf (string,"trig_fcalbcal2.pdf");
592  c1->SaveAs(string);
593  sprintf (string,"trig_fcalbcal2.pdf");
594  c2->SaveAs(string);
595  sprintf (string,"trig_fcalbcal2.pdf");
596  c3->SaveAs(string);
597  sprintf (string,"trig_fcalbcal2.pdf");
598  c4->SaveAs(string);
599  sprintf (string,"trig_fcalbcal2.pdf)");
600  c5->SaveAs(string);
601 
602 
603 
604 }
605 
606 
TH1F * h1trig5_fcal
Double_t c0[2][NMODULES]
Definition: tw_corr.C:67
TLine * thresh_b
TH1F * h2trig3_fcalVSbcal
TH1F * h1trig7_fcalN
TH1F * h1trig1_bcalN
TH1F * h1trig1_fcalN
sprintf(text,"Post KinFit Cut")
TString filename
TH1F * h2trig1_fcalVSbcal
TH1F * h1trig_bcalN
Definition: trig_fcalbcal.C:32
TH1F * h1trig_trgbits
double thb
TH1F * h1trig5_fcalN
TLine * thresh_f
TH1F * h1trig3_fcal
Double_t c1[2][NMODULES]
Definition: tw_corr.C:68
double bcal_max
TF1 * f
Definition: FitGains.C:21
TH1F * h2trig_fcalVSbcal
Definition: trig_fcalbcal.C:33
Double_t c2[2][NMODULES]
Definition: tw_corr.C:69
TH1F * h2trig7_fcalVSbcal
TH1F * h1trig1_bcal
TH1F * h2trig5_fcalVSbcal
TH1F * h1trig7_bcalN
TH1F * h1trig1_fcal
TLatex * t1
TH1F * h1trig3_fcalN
TH1F * h1trig7_bcal
double thf
TH1F * h1trig_fcal
Definition: trig_fcalbcal.C:29
TH1F * h1trig3_bcal
TLine * thresh_fb
TH1F * h1trig3_bcalN
TH1F * h1trig7_fcal
TCanvas * c4
TH1F * h1trig5_bcal
TH1F * h1trig_fcalN
Definition: trig_fcalbcal.C:30
Double_t ymax
Definition: bcal_hist_eff.C:91
double fcal_max
TDirectory * dir
Definition: bcal_hist_eff.C:31
TH1F * h1trig_bcal
Definition: trig_fcalbcal.C:31
TCanvas * c3
printf("string=%s", string)
TH2F * temp
TLine * thrf
TLine * thrb
TH1F * h1trig5_bcalN