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);
17 TFile *
f =
new TFile(filename+
".root",
"read");
19 cout <<
"Opening root file: " << (filename+
".root").Data() << endl;
21 TH1F *NumEventsSurvivedAction = (TH1F*)f->Get(
"NumEventsSurvivedAction");
22 TH2F *NumCombosSurvivedAction = (TH2F*)f->Get(
"NumCombosSurvivedAction");
23 TH1F *BeamEnergy = (TH1F*)f->Get(
"BeamEnergy");
25 TH1F *piplusMomentumMeasured = (TH1F*)f->Get(
"piplusMomentumMeasured");
26 TH1F *piminusMomentumMeasured = (TH1F*)f->Get(
"piminusMomentumMeasured");
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");
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");
42 TH1F *pipDeltap = (TH1F*)f->Get(
"pipDeltap");
43 TH1F *pimDeltap = (TH1F*)f->Get(
"pimDeltap");
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");
52 TH1F *dHist_TaggerAccidentals = (TH1F*)f->Get(
"dHist_TaggerAccidentals");
54 TCanvas *
c0 =
new TCanvas(
"c0",
"c0",200,10,1000,700);
62 Double_t
ymax = 10000;
64 NumEventsSurvivedAction->SetTitle(filename);
67 NumEventsSurvivedAction->GetXaxis()->SetTitleSize(0.05);
68 NumEventsSurvivedAction->GetYaxis()->SetTitleSize(0.05);
70 NumEventsSurvivedAction->SetMarkerColor(4);
71 NumEventsSurvivedAction->Draw();
78 NumCombosSurvivedAction->SetTitle(filename);
81 NumCombosSurvivedAction->GetXaxis()->SetTitleSize(0.05);
82 NumCombosSurvivedAction->GetYaxis()->SetTitleSize(0.05);
85 NumCombosSurvivedAction->Draw(
"colz");
92 BeamEnergy->SetTitle(filename);
95 BeamEnergy->GetXaxis()->SetTitleSize(0.05);
96 BeamEnergy->GetYaxis()->SetTitleSize(0.05);
97 BeamEnergy->GetXaxis()->SetTitle(
"Energy (GeV)");
98 BeamEnergy->SetMarkerColor(4);
106 CosTheta->SetTitle(filename);
109 CosTheta->GetXaxis()->SetTitleSize(0.05);
110 CosTheta->GetYaxis()->SetTitleSize(0.05);
111 CosTheta->GetXaxis()->SetTitle(
"Cos#Theta");
112 CosTheta->SetMarkerColor(4);
120 piplusMomentumMeasured->SetTitle(filename);
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();
134 piminusMomentumMeasured->SetTitle(filename);
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();
144 TCanvas *
c1 =
new TCanvas(
"c1",
"c1",200,10,1000,700);
154 KinFitChiSq->SetTitle(filename);
157 KinFitChiSq->GetXaxis()->SetTitleSize(0.05);
158 KinFitChiSq->GetYaxis()->SetTitleSize(0.05);
159 KinFitChiSq->GetXaxis()->SetTitle(
"#chi^{2}");
160 KinFitChiSq->SetMarkerColor(4);
168 KinFitCL->SetTitle(filename);
171 KinFitCL->GetXaxis()->SetTitleSize(0.05);
172 KinFitCL->GetYaxis()->SetTitleSize(0.05);
173 KinFitCL->GetXaxis()->SetTitle(
"Probability");
174 KinFitCL->SetMarkerColor(4);
182 MissingMassSquared->SetTitle(filename);
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();
196 M2pikin->SetTitle(filename);
199 M2pikin->GetXaxis()->SetTitleSize(0.05);
200 M2pikin->GetYaxis()->SetTitleSize(0.05);
201 M2pikin->GetXaxis()->SetTitle(
"Mass 2pi Kin (GeV)");
202 M2pikin->SetMarkerColor(4);
210 tkin->SetTitle(filename);
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);
225 CosThetadiff->SetTitle(filename);
228 CosThetadiff->GetXaxis()->SetTitleSize(0.05);
229 CosThetadiff->GetYaxis()->SetTitleSize(0.05);
231 CosThetadiff->GetXaxis()->SetTitle(
"Cos#Theta diff Kin-Gen");
232 CosThetadiff->SetMarkerColor(4);
233 CosThetadiff->Draw();
236 TCanvas *
c2 =
new TCanvas(
"c2",
"c2",200,10,1000,700);
245 M2pidiff->SetTitle(filename);
248 M2pidiff->GetXaxis()->SetTitleSize(0.05);
249 M2pidiff->GetYaxis()->SetTitleSize(0.05);
251 M2pidiff->SetMarkerColor(4);
252 M2pidiff->GetXaxis()->SetNdivisions(505);
253 M2pidiff->Fit(
"gaus");
261 pipDeltap_Measured->SetTitle(filename);
264 pipDeltap_Measured->GetXaxis()->SetTitleSize(0.05);
265 pipDeltap_Measured->GetYaxis()->SetTitleSize(0.05);
267 pipDeltap_Measured->SetMarkerColor(4);
268 pipDeltap_Measured->Draw();
275 pimDeltap_Measured->SetTitle(filename);
278 pimDeltap_Measured->GetXaxis()->SetTitleSize(0.05);
279 pimDeltap_Measured->GetYaxis()->SetTitleSize(0.05);
281 pimDeltap_Measured->SetMarkerColor(4);
282 pimDeltap_Measured->Draw();
289 tdiff->SetTitle(filename);
292 tdiff->GetXaxis()->SetTitleSize(0.05);
293 tdiff->GetYaxis()->SetTitleSize(0.05);
295 tdiff->SetMarkerColor(4);
296 tdiff->GetXaxis()->SetNdivisions(505);
305 pipDeltap->SetTitle(filename);
308 pipDeltap->GetXaxis()->SetTitleSize(0.05);
309 pipDeltap->GetYaxis()->SetTitleSize(0.05);
311 pipDeltap->SetMarkerColor(4);
319 pimDeltap->SetTitle(filename);
322 pimDeltap->GetXaxis()->SetTitleSize(0.05);
323 pimDeltap->GetYaxis()->SetTitleSize(0.05);
325 pimDeltap->SetMarkerColor(4);
329 TCanvas *
c3 =
new TCanvas(
"c3",
"c3",200,10,1000,700);
338 Phikin->SetTitle(filename);
341 Phikin->GetXaxis()->SetTitleSize(0.05);
342 Phikin->GetYaxis()->SetTitleSize(0.05);
343 Phikin->GetXaxis()->SetTitle(
"#Phi Kin (Degrees)");
344 Phikin->SetMarkerColor(4);
352 phikin->SetTitle(filename);
355 phikin->GetXaxis()->SetTitleSize(0.05);
356 phikin->GetYaxis()->SetTitleSize(0.05);
357 phikin->GetXaxis()->SetTitle(
"#phi Kin (Degrees)");
358 phikin->SetMarkerColor(4);
366 psikin->SetTitle(filename);
369 psikin->GetXaxis()->SetTitleSize(0.05);
370 psikin->GetYaxis()->SetTitleSize(0.05);
371 psikin->GetXaxis()->SetTitle(
"#psii Kin (Degrees)");
372 psikin->SetMarkerColor(4);
380 Phidiff->SetTitle(filename);
383 Phidiff->GetXaxis()->SetTitleSize(0.05);
384 Phidiff->GetYaxis()->SetTitleSize(0.05);
385 Phidiff->GetXaxis()->SetTitle(
"#Phi Kin-Gen (Degrees)");
386 Phidiff->SetMarkerColor(4);
394 phidiff->SetTitle(filename);
397 phidiff->GetXaxis()->SetTitleSize(0.05);
398 phidiff->GetYaxis()->SetTitleSize(0.05);
399 phidiff->GetXaxis()->SetTitle(
"#phi Kin-Gen (Degrees)");
400 phidiff->SetMarkerColor(4);
408 psidiff->SetTitle(filename);
411 psidiff->GetXaxis()->SetTitleSize(0.05);
412 psidiff->GetYaxis()->SetTitleSize(0.05);
413 psidiff->GetXaxis()->SetTitle(
"#psi Kin-Gen (Degrees)");
414 psidiff->SetMarkerColor(4);
417 TCanvas *
c4 =
new TCanvas(
"c4",
"c4",200,10,1000,700);
426 dHist_TaggerAccidentals->SetTitle(filename);
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();
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)");
void plot_Z2pi_trees(TString filename)