Prepare for the new TPC layout
[u/mrichter/AliRoot.git] / EVE / macros / show_scan_results.C
1 // Display some histograms from scanning.
2 //
3 // BIT(1) stores the original selection.
4 // BIT(0) stores the user selection (set to same value as b1 at init).
5 //
6 // This allows to check all possible combinations.
7
8
9 void show_scan_results()
10 {
11   TFile *f = TFile::Open("scan_results.root");
12
13   TTree* t = (TTree*) gDirectory->Get("SR");
14
15   if (t == 0)
16     Error("show_scan_results", "Tree 'SR' with scan results not found.");
17
18   TCanvas *c = 0;
19
20
21   //----------------------------------------------------------------------------
22   // Tracks
23   //----------------------------------------------------------------------------
24
25   c = new TCanvas("Tracks", "Track Scanning Results", 800, 600);
26   c->Divide(2, 3);
27
28   c->cd(1);
29   t->Draw("Sum$(T.fLabel & 1)");
30
31   c->cd(2);
32   t->Draw("T.GetSign()", "T.fLabel & 1");
33
34   c->cd(3);
35   t->Draw("T.Pt()", "T.fLabel & 1");
36
37   c->cd(4);
38   t->Draw("T.Eta()", "T.fLabel & 1");
39
40   c->cd(5);
41   t->Draw("T.Phi()", "T.fLabel & 1");
42
43   c->Modified();
44   c->Update();
45
46
47   //----------------------------------------------------------------------------
48   // Trackelts
49   //----------------------------------------------------------------------------
50
51   c = new TCanvas("Tracklets", "Tracklet Scanning Results", 800, 600);
52   c->Divide(2, 3);
53
54   c->cd(1);
55   t->Draw("Sum$(M.fLabels & 1)");
56
57   c->cd(2);
58   t->Draw("M.fNsingle");
59
60   c->cd(3);
61   t->Draw("M.fFiredChips[1]:Sum$(M.fLabels & 1)");
62
63   c->cd(4);
64   t->Draw("M.fDeltTh", "M.fLabels & 1");
65
66   c->cd(5);
67   t->Draw("M.fDeltPhi", "M.fLabels & 1");
68
69   c->cd(6);
70   t->Draw("M.fPhi", "M.fLabels & 1");
71
72   c->Modified();
73   c->Update();
74
75
76   //----------------------------------------------------------------------------
77   // Vertices
78   //----------------------------------------------------------------------------
79
80   c = new TCanvas("Vertices", "Vertex Scanning Results", 800, 600);
81   c->Divide(3, 3);
82
83   c->cd(1);
84   t->Draw("VT.GetX()", "VT.GetNContributors()>0");
85
86   c->cd(2);
87   t->Draw("VT.GetY()", "VT.GetNContributors()>0");
88
89   c->cd(3);
90   t->Draw("VT.GetZ()", "VT.GetNContributors()>0");
91
92   c->cd(4);
93   t->Draw("VSPD.GetX()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
94
95   c->cd(5);
96   t->Draw("VSPD.GetY()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
97
98   c->cd(6);
99   t->Draw("VSPD.GetZ()", "VSPD.GetNContributors()>0 && VSPD.fTitle.Contains(\"3D\")");
100
101   c->cd(7);
102   t->Draw("VTPC.GetX()", "VTPC.GetNContributors()>0");
103
104   c->cd(8);
105   t->Draw("VTPC.GetY()", "VTPC.GetNContributors()>0");
106
107   c->cd(9);
108   t->Draw("VTPC.GetZ()", "VTPC.GetNContributors()>0");
109
110   //----------------------------------------------------------------------------
111   // End
112   //----------------------------------------------------------------------------
113
114   f->Close();
115   delete f;
116 }