Hall-D Software  alpha
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
plot_Z2pi_trees.C
Go to the documentation of this file.
1 void plot_Z2pi_trees(TString filename)
2 {
3 // File: Z2pi_trees.C
4 //
5 
6  gStyle->SetPalette(1,0);
7  gStyle->SetOptStat(111111);
8  gStyle->SetOptFit(111111);
9  gStyle->SetPadRightMargin(0.15);
10  gStyle->SetPadLeftMargin(0.15);
11  gStyle->SetPadBottomMargin(0.15);
12 
13  // TString filename = "DSelector_Z2pi_trees_flat";
14  // TString filename = "DSelector_Z2pi_trees_signal";
15 
16 
17  TFile *f = new TFile(filename+".root","read");
18 
19  cout << "Opening root file: " << (filename+".root").Data() << endl;
20 
21  TH1F *NumEventsSurvivedAction = (TH1F*)f->Get("NumEventsSurvivedAction");
22  TH2F *NumCombosSurvivedAction = (TH2F*)f->Get("NumCombosSurvivedAction");
23  TH1F *BeamEnergy = (TH1F*)f->Get("BeamEnergy");
24  // TH1F *pMomentumMeasured = (TH1F*)f->Get("pMomentumMeasured");
25  TH1F *piplusMomentumMeasured = (TH1F*)f->Get("piplusMomentumMeasured");
26  TH1F *piminusMomentumMeasured = (TH1F*)f->Get("piminusMomentumMeasured");
27 
28  TH1F *KinFitChiSq = (TH1F*)f->Get("KinFitChiSq");
29  TH1F *KinFitCL = (TH1F*)f->Get("KinFitCL");
30  TH1F *MissingMassSquared = (TH1F*)f->Get("MissingMassSquared");
31  TH1F *M2pikin = (TH1F*)f->Get("M2pikin");
32  TH1F *tkin = (TH1F*)f->Get("tkin");
33  TH1F *psikin = (TH1F*)f->Get("psikin");
34  TH2F *CosTheta_Psi = (TH2F*)f->Get("CosTheta_Psi");
35 
36  TH1F *M2pidiff = (TH1F*)f->Get("M2pidiff");
37  TH1F *tdiff = (TH1F*)f->Get("tdiff");
38  TH1F *psidiff = (TH1F*)f->Get("psidiff");
39  TH1F *pipDeltap_Measured = (TH1F*)f->Get("pipDeltap_Measured");
40  TH1F *pimDeltap_Measured = (TH1F*)f->Get("pimDeltap_Measured");
41  // TH1F *pDeltap = (TH1F*)f->Get("pDeltap");
42  TH1F *pipDeltap = (TH1F*)f->Get("pipDeltap");
43  TH1F *pimDeltap = (TH1F*)f->Get("pimDeltap");
44 
45  TH1F *Phikin = (TH1F*)f->Get("Phikin");
46  TH1F *phikin = (TH1F*)f->Get("phikin");
47  TH1F *CosTheta = (TH1F*)f->Get("CosTheta");
48  TH1F *Phidiff = (TH1F*)f->Get("Phidiff");
49  TH1F *phidiff = (TH1F*)f->Get("phidiff");
50  TH1F *CosThetadiff = (TH1F*)f->Get("CosThetadiff");
51 
52  TH1F *dHist_TaggerAccidentals = (TH1F*)f->Get("dHist_TaggerAccidentals");
53 
54  TCanvas *c0 = new TCanvas("c0", "c0",200,10,1000,700);
55 
56  c0->Divide(3,2);
57  c0->cd(1);
58  // gPad->SetLogy();
59  Double_t xmin = 0;
60  Double_t xmax = 2;
61  Double_t ymin = 100;
62  Double_t ymax = 10000;
63 
64  NumEventsSurvivedAction->SetTitle(filename);
65  // NumEventsSurvivedAction->GetXaxis()->SetRangeUser(xmin,xmax);
66  // NumEventsSurvivedAction->GetYaxis()->SetRangeUser(ymin,ymax);
67  NumEventsSurvivedAction->GetXaxis()->SetTitleSize(0.05);
68  NumEventsSurvivedAction->GetYaxis()->SetTitleSize(0.05);
69  // NumEventsSurvivedAction->GetXaxis()->SetTitle("Events");
70  NumEventsSurvivedAction->SetMarkerColor(4);
71  NumEventsSurvivedAction->Draw();
72 
73  c0->cd(2);
74  // gPad->SetLogy();
75  ymin = 0;
76  ymax = 5000;
77 
78  NumCombosSurvivedAction->SetTitle(filename);
79  // NumCombosSurvivedAction->GetXaxis()->SetRangeUser(xmin,xmax);
80  // NumCombosSurvivedAction->GetYaxis()->SetRangeUser(ymin,ymax);
81  NumCombosSurvivedAction->GetXaxis()->SetTitleSize(0.05);
82  NumCombosSurvivedAction->GetYaxis()->SetTitleSize(0.05);
83  // NumCombosSurvivedAction->GetXaxis()->SetTitle("Events");
84  // NumCombosSurvivedAction->SetMarkerColor(4);
85  NumCombosSurvivedAction->Draw("colz");
86 
87  c0->cd(3);
88  // gPad->SetLogy();
89  ymin = 0;
90  ymax = 4000;
91 
92  BeamEnergy->SetTitle(filename);
93  // BeamEnergy->GetXaxis()->SetRangeUser(xmin,xmax);
94  // BeamEnergy->GetYaxis()->SetRangeUser(ymin,ymax);
95  BeamEnergy->GetXaxis()->SetTitleSize(0.05);
96  BeamEnergy->GetYaxis()->SetTitleSize(0.05);
97  BeamEnergy->GetXaxis()->SetTitle("Energy (GeV)");
98  BeamEnergy->SetMarkerColor(4);
99  BeamEnergy->Draw();
100 
101  c0->cd(4);
102  // gPad->SetLogy();
103  ymin = 0;
104  ymax = 4000;
105 
106  CosTheta->SetTitle(filename);
107  // CosTheta->GetXaxis()->SetRangeUser(xmin,xmax);
108  // CosTheta->GetYaxis()->SetRangeUser(ymin,ymax);
109  CosTheta->GetXaxis()->SetTitleSize(0.05);
110  CosTheta->GetYaxis()->SetTitleSize(0.05);
111  CosTheta->GetXaxis()->SetTitle("Cos#Theta");
112  CosTheta->SetMarkerColor(4);
113  CosTheta->Draw();
114 
115  c0->cd(5);
116  // gPad->SetLogy();
117  ymin = 0;
118  ymax = 4000;
119 
120  piplusMomentumMeasured->SetTitle(filename);
121  // piplusMomentumMeasured->GetXaxis()->SetRangeUser(xmin,xmax);
122  // piplusMomentumMeasured->GetYaxis()->SetRangeUser(ymin,ymax);
123  piplusMomentumMeasured->GetXaxis()->SetTitleSize(0.05);
124  piplusMomentumMeasured->GetYaxis()->SetTitleSize(0.05);
125  piplusMomentumMeasured->GetXaxis()->SetTitle("PiPlus Momentum (GeV)");
126  piplusMomentumMeasured->SetMarkerColor(4);
127  piplusMomentumMeasured->Draw();
128 
129  c0->cd(6);
130  // gPad->SetLogy();
131  ymin = 0;
132  ymax = 4000;
133 
134  piminusMomentumMeasured->SetTitle(filename);
135  // piminusMomentumMeasured->GetXaxis()->SetRangeUser(xmin,xmax);
136  // piminusMomentumMeasured->GetYaxis()->SetRangeUser(ymin,ymax);
137  piminusMomentumMeasured->GetXaxis()->SetTitleSize(0.05);
138  piminusMomentumMeasured->GetYaxis()->SetTitleSize(0.05);
139  piminusMomentumMeasured->GetXaxis()->SetTitle("PiMinus Momentum (GeV)");
140  piminusMomentumMeasured->SetMarkerColor(4);
141  piminusMomentumMeasured->Draw();
142 
143 
144  TCanvas *c1 = new TCanvas("c1", "c1",200,10,1000,700);
145 
146  c1->Divide(3,2);
147 
148 
149  c1->cd(1);
150  // gPad->SetLogy();
151  ymin = 0;
152  ymax = 4000;
153 
154  KinFitChiSq->SetTitle(filename);
155  // KinFitChiSq->GetXaxis()->SetRangeUser(xmin,xmax);
156  // KinFitChiSq->GetYaxis()->SetRangeUser(ymin,ymax);
157  KinFitChiSq->GetXaxis()->SetTitleSize(0.05);
158  KinFitChiSq->GetYaxis()->SetTitleSize(0.05);
159  KinFitChiSq->GetXaxis()->SetTitle("#chi^{2}");
160  KinFitChiSq->SetMarkerColor(4);
161  KinFitChiSq->Draw();
162 
163  c1->cd(2);
164  gPad->SetLogy();
165  ymin = 0;
166  ymax = 4000;
167 
168  KinFitCL->SetTitle(filename);
169  // KinFitCL->GetXaxis()->SetRangeUser(xmin,xmax);
170  // KinFitCL->GetYaxis()->SetRangeUser(ymin,ymax);
171  KinFitCL->GetXaxis()->SetTitleSize(0.05);
172  KinFitCL->GetYaxis()->SetTitleSize(0.05);
173  KinFitCL->GetXaxis()->SetTitle("Probability");
174  KinFitCL->SetMarkerColor(4);
175  KinFitCL->Draw();
176 
177  c1->cd(3);
178  // gPad->SetLogy();
179  ymin = 0;
180  ymax = 4000;
181 
182  MissingMassSquared->SetTitle(filename);
183  // MissingMassSquared->GetXaxis()->SetRangeUser(xmin,xmax);
184  // MissingMassSquared->GetYaxis()->SetRangeUser(ymin,ymax);
185  MissingMassSquared->GetXaxis()->SetTitleSize(0.05);
186  MissingMassSquared->GetYaxis()->SetTitleSize(0.05);
187  MissingMassSquared->GetXaxis()->SetTitle("Missing Mass (GeV2)");
188  MissingMassSquared->SetMarkerColor(4);
189  MissingMassSquared->Draw();
190 
191  c1->cd(4);
192  // gPad->SetLogy();
193  ymin = 0;
194  ymax = 4000;
195 
196  M2pikin->SetTitle(filename);
197  // M2pikin->GetXaxis()->SetRangeUser(xmin,xmax);
198  // M2pikin->GetYaxis()->SetRangeUser(ymin,ymax);
199  M2pikin->GetXaxis()->SetTitleSize(0.05);
200  M2pikin->GetYaxis()->SetTitleSize(0.05);
201  M2pikin->GetXaxis()->SetTitle("Mass 2pi Kin (GeV)");
202  M2pikin->SetMarkerColor(4);
203  M2pikin->Draw();
204 
205  c1->cd(5);
206  gPad->SetLogy();
207  ymin = 0;
208  ymax = 4000;
209 
210  tkin->SetTitle(filename);
211  // tkin->GetXaxis()->SetRangeUser(xmin,xmax);
212  // tkin->GetYaxis()->SetRangeUser(ymin,ymax);
213  tkin->GetXaxis()->SetTitleSize(0.05);
214  tkin->GetYaxis()->SetTitleSize(0.05);
215  tkin->GetXaxis()->SetTitle("-t Kin (GeV2)");
216  tkin->GetXaxis()->SetNdivisions(505);
217  tkin->SetMarkerColor(4);
218  tkin->Draw();
219 
220  c1->cd(6);
221  // gPad->SetLogy();
222  ymin = 0;
223  ymax = 4000;
224 
225  CosThetadiff->SetTitle(filename);
226  // CosThetadiff->GetXaxis()->SetRangeUser(xmin,xmax);
227  // CosThetadiff->GetYaxis()->SetRangeUser(ymin,ymax);
228  CosThetadiff->GetXaxis()->SetTitleSize(0.05);
229  CosThetadiff->GetYaxis()->SetTitleSize(0.05);
230  // CosThetadiff->GetYaxis()->SetTitle("CosTheta");
231  CosThetadiff->GetXaxis()->SetTitle("Cos#Theta diff Kin-Gen");
232  CosThetadiff->SetMarkerColor(4);
233  CosThetadiff->Draw();
234 
235 
236  TCanvas *c2 = new TCanvas("c2", "c2",200,10,1000,700);
237 
238  c2->Divide(3,2);
239 
240  c2->cd(1);
241  // gPad->SetLogy();
242  ymin = 0;
243  ymax = 4000;
244 
245  M2pidiff->SetTitle(filename);
246  // M2pidiff->GetXaxis()->SetRangeUser(xmin,xmax);
247  // M2pidiff->GetYaxis()->SetRangeUser(ymin,ymax);
248  M2pidiff->GetXaxis()->SetTitleSize(0.05);
249  M2pidiff->GetYaxis()->SetTitleSize(0.05);
250  // M2pidiff->GetXaxis()->SetTitle("");
251  M2pidiff->SetMarkerColor(4);
252  M2pidiff->GetXaxis()->SetNdivisions(505);
253  M2pidiff->Fit("gaus");
254  M2pidiff->Draw();
255 
256  c2->cd(2);
257  // gPad->SetLogy();
258  ymin = 0;
259  ymax = 4000;
260 
261  pipDeltap_Measured->SetTitle(filename);
262  // pipDeltap_Measured->GetXaxis()->SetRangeUser(xmin,xmax);
263  // pipDeltap_Measured->GetYaxis()->SetRangeUser(ymin,ymax);
264  pipDeltap_Measured->GetXaxis()->SetTitleSize(0.05);
265  pipDeltap_Measured->GetYaxis()->SetTitleSize(0.05);
266  // pipDeltap_Measured->GetXaxis()->SetTitle("#chi^{2}");
267  pipDeltap_Measured->SetMarkerColor(4);
268  pipDeltap_Measured->Draw();
269 
270  c2->cd(3);
271  // gPad->SetLogy();
272  ymin = 0;
273  ymax = 4000;
274 
275  pimDeltap_Measured->SetTitle(filename);
276  // pimDeltap_Measured->GetXaxis()->SetRangeUser(xmin,xmax);
277  // pimDeltap_Measured->GetYaxis()->SetRangeUser(ymin,ymax);
278  pimDeltap_Measured->GetXaxis()->SetTitleSize(0.05);
279  pimDeltap_Measured->GetYaxis()->SetTitleSize(0.05);
280  // pimDeltap_Measured->GetXaxis()->SetTitle("#chi^{2}");
281  pimDeltap_Measured->SetMarkerColor(4);
282  pimDeltap_Measured->Draw();
283 
284  c2->cd(4);
285  // gPad->SetLogy();
286  ymin = 0;
287  ymax = 4000;
288 
289  tdiff->SetTitle(filename);
290  // tdiff->GetXaxis()->SetRangeUser(xmin,xmax);
291  // tdiff->GetYaxis()->SetRangeUser(ymin,ymax);
292  tdiff->GetXaxis()->SetTitleSize(0.05);
293  tdiff->GetYaxis()->SetTitleSize(0.05);
294  // tdiff->GetXaxis()->SetTitle("");
295  tdiff->SetMarkerColor(4);
296  tdiff->GetXaxis()->SetNdivisions(505);
297  tdiff->Fit("gaus");
298  tdiff->Draw();
299 
300  c2->cd(5);
301  // gPad->SetLogy();
302  ymin = 0;
303  ymax = 4000;
304 
305  pipDeltap->SetTitle(filename);
306  // pipDeltap->GetXaxis()->SetRangeUser(xmin,xmax);
307  // pipDeltap->GetYaxis()->SetRangeUser(ymin,ymax);
308  pipDeltap->GetXaxis()->SetTitleSize(0.05);
309  pipDeltap->GetYaxis()->SetTitleSize(0.05);
310  // pipDeltap->GetXaxis()->SetTitle("#chi^{2}");
311  pipDeltap->SetMarkerColor(4);
312  pipDeltap->Draw();
313 
314  c2->cd(6);
315  // gPad->SetLogy();
316  ymin = 0;
317  ymax = 4000;
318 
319  pimDeltap->SetTitle(filename);
320  // pimDeltap->GetXaxis()->SetRangeUser(xmin,xmax);
321  // pimDeltap->GetYaxis()->SetRangeUser(ymin,ymax);
322  pimDeltap->GetXaxis()->SetTitleSize(0.05);
323  pimDeltap->GetYaxis()->SetTitleSize(0.05);
324  // pimDeltap->GetXaxis()->SetTitle("#chi^{2}");
325  pimDeltap->SetMarkerColor(4);
326  pimDeltap->Draw();
327 
328 
329  TCanvas *c3 = new TCanvas("c3", "c3",200,10,1000,700);
330 
331  c3->Divide(3,2);
332 
333  c3->cd(1);
334  // gPad->SetLogy();
335  ymin = 0;
336  ymax = 4000;
337 
338  Phikin->SetTitle(filename);
339  // Phikin->GetXaxis()->SetRangeUser(xmin,xmax);
340  // Phikin->GetYaxis()->SetRangeUser(ymin,ymax);
341  Phikin->GetXaxis()->SetTitleSize(0.05);
342  Phikin->GetYaxis()->SetTitleSize(0.05);
343  Phikin->GetXaxis()->SetTitle("#Phi Kin (Degrees)");
344  Phikin->SetMarkerColor(4);
345  Phikin->Draw();
346 
347  c3->cd(2);
348  // gPad->SetLogy();
349  ymin = 0;
350  ymax = 4000;
351 
352  phikin->SetTitle(filename);
353  // phikin->GetXaxis()->SetRangeUser(xmin,xmax);
354  // phikin->GetYaxis()->SetRangeUser(ymin,ymax);
355  phikin->GetXaxis()->SetTitleSize(0.05);
356  phikin->GetYaxis()->SetTitleSize(0.05);
357  phikin->GetXaxis()->SetTitle("#phi Kin (Degrees)");
358  phikin->SetMarkerColor(4);
359  phikin->Draw();
360 
361  c3->cd(3);
362  // gPad->SetLogy();
363  ymin = 0;
364  ymax = 4000;
365 
366  psikin->SetTitle(filename);
367  // psikin->GetXaxis()->SetRangeUser(xmin,xmax);
368  // psikin->GetYaxis()->SetRangeUser(ymin,ymax);
369  psikin->GetXaxis()->SetTitleSize(0.05);
370  psikin->GetYaxis()->SetTitleSize(0.05);
371  psikin->GetXaxis()->SetTitle("#psii Kin (Degrees)");
372  psikin->SetMarkerColor(4);
373  psikin->Draw();
374 
375  c3->cd(4);
376  // gPad->SetLogy();
377  ymin = 0;
378  ymax = 4000;
379 
380  Phidiff->SetTitle(filename);
381  // Phidiff->GetXaxis()->SetRangeUser(xmin,xmax);
382  // Phidiff->GetYaxis()->SetRangeUser(ymin,ymax);
383  Phidiff->GetXaxis()->SetTitleSize(0.05);
384  Phidiff->GetYaxis()->SetTitleSize(0.05);
385  Phidiff->GetXaxis()->SetTitle("#Phi Kin-Gen (Degrees)");
386  Phidiff->SetMarkerColor(4);
387  Phidiff->Draw();
388 
389  c3->cd(5);
390  // gPad->SetLogy();
391  ymin = 0;
392  ymax = 4000;
393 
394  phidiff->SetTitle(filename);
395  // phidiff->GetXaxis()->SetRangeUser(xmin,xmax);
396  // phidiff->GetYaxis()->SetRangeUser(ymin,ymax);
397  phidiff->GetXaxis()->SetTitleSize(0.05);
398  phidiff->GetYaxis()->SetTitleSize(0.05);
399  phidiff->GetXaxis()->SetTitle("#phi Kin-Gen (Degrees)");
400  phidiff->SetMarkerColor(4);
401  phidiff->Draw();
402 
403  c3->cd(6);
404  // gPad->SetLogy();
405  ymin = 0;
406  ymax = 4000;
407 
408  psidiff->SetTitle(filename);
409  // psidiff->GetXaxis()->SetRangeUser(xmin,xmax);
410  // psidiff->GetYaxis()->SetRangeUser(ymin,ymax);
411  psidiff->GetXaxis()->SetTitleSize(0.05);
412  psidiff->GetYaxis()->SetTitleSize(0.05);
413  psidiff->GetXaxis()->SetTitle("#psi Kin-Gen (Degrees)");
414  psidiff->SetMarkerColor(4);
415  psidiff->Draw();
416 
417  TCanvas *c4 = new TCanvas("c4", "c4",200,10,1000,700);
418 
419  c4->Divide(3,2);
420 
421  c4->cd(1);
422  // gPad->SetLogy();
423  ymin = 0;
424  ymax = 4000;
425 
426  dHist_TaggerAccidentals->SetTitle(filename);
427  // dHist_TaggerAccidentals->GetXaxis()->SetRangeUser(xmin,xmax);
428  // dHist_TaggerAccidentals->GetYaxis()->SetRangeUser(ymin,ymax);
429  dHist_TaggerAccidentals->GetXaxis()->SetTitleSize(0.05);
430  dHist_TaggerAccidentals->GetYaxis()->SetTitleSize(0.05);
431  dHist_TaggerAccidentals->GetXaxis()->SetTitle("Vertex Time - RF (ns)");
432  dHist_TaggerAccidentals->SetMarkerColor(4);
433  dHist_TaggerAccidentals->Draw();
434 
435  c0->SaveAs(filename+".pdf(");
436  c1->SaveAs(filename+".pdf");
437  c2->SaveAs(filename+".pdf");
438  c3->SaveAs(filename+".pdf");
439  c4->SaveAs(filename+".pdf)");
440 }
Double_t c0[2][NMODULES]
Definition: tw_corr.C:67
TString filename
void plot_Z2pi_trees(TString filename)
Double_t c1[2][NMODULES]
Definition: tw_corr.C:68
TF1 * f
Definition: FitGains.C:21
Double_t c2[2][NMODULES]
Definition: tw_corr.C:69
Double_t ymin
Definition: bcal_hist_eff.C:89
TCanvas * c4
Double_t ymax
Definition: bcal_hist_eff.C:91
TCanvas * c3