New developments of the analysis framework - selectorised version of the manager...
[u/mrichter/AliRoot.git] / LHC / PlotPressure.C
1 void PlotPressure()
2 {
3     FILE* file;
4     file=fopen("gasPressure.dat","r");
5     Float_t z1[20];
6     Float_t g1[20], g2[20], g3[20], g4[20], g5[20];
7     Float_t z2[21];
8     Float_t h1[21], h2[21], h3[21], h4[21], h5[21];
9
10     char c[45];
11     Float_t z;
12     
13     for (Int_t i = 0; i < 20; i++)
14     {
15         fscanf(file, "%f %f %f %f %f %f", &z, 
16                &g1[i], &g2[i], &g3[i], &g4[i], &g5[i]);
17 //      printf("%d %f %f %f %f %f %f \n", i, z1[i], 
18 //                     g1[i][0], g1[i][1], g1[i][2], g1[i][3], g1[i][4]);
19         if (i > 0) {
20             z1[i] = z1[i-1] + z;
21         } else {
22             z1[i] = 20.;
23         }
24     }
25
26     
27     for (Int_t i = 0; i < 21; i++)
28     {
29         fscanf(file, "%f %f %f %f %f %f", &z, 
30                &h1[i], &h2[i], &h3[i], &h4[i], &h5[i]);
31         if (i > 0) {
32             z2[i] = z2[i-1] + z;
33         } else {
34             z2[i] = 20.;
35         }
36     }
37
38 //
39 // 
40     TCanvas *c1 = new TCanvas("c1","Gas Pressure Beam 1", 200, 10, 700, 500);
41     gPad->SetLogy();
42     
43     TGraph* gr1  = new TGraph(20, z1, g1);
44     gr1->SetMaximum(1e17);
45     gr1->SetMinimum(1e11);    
46     gr1->SetLineColor(1);
47     gr1->SetTitle("Ring 1: Beginning of Run");
48     
49     TGraph* gr2  = new TGraph(20, z1, g3);
50     gr2->SetLineColor(2);
51
52     TGraph* gr3  = new TGraph(20, z1, g5);
53     gr3->SetLineColor(4);
54
55     gr1->Draw("AL");
56     gr2->Draw("L");
57 //    gr3->Draw("L");
58     text();
59
60 //
61 // 
62     TCanvas *c2 = new TCanvas("c2","Gas Pressure Beam 1", 200, 10, 700, 500);
63     gPad->SetLogy();
64     
65     TGraph* gr4  = new TGraph(20, z1, g2);
66     gr4->SetMaximum(1e17);
67     gr4->SetMinimum(1e11);    
68     gr4->SetLineColor(1);
69     gr4->SetTitle("Ring 1");
70     
71     TGraph* gr5  = new TGraph(20, z1, g4);
72     gr5->SetLineColor(2);
73
74     TGraph* gr6  = new TGraph(20, z1, g5);
75     gr6->SetLineColor(4);
76
77     gr4->Draw("AL");
78     gr5->Draw("L");
79     gr6->Draw("L");
80     text();
81
82 //
83 // 
84     TCanvas *c3 = new TCanvas("c3","Gas Pressure Beam 2", 200, 10, 700, 500);
85     gPad->SetLogy();
86     
87     TGraph* hr1  = new TGraph(21, z2, h1);
88     hr1->SetMaximum(1e17);
89     hr1->SetMinimum(1e11);    
90     hr1->SetLineColor(1);
91     hr1->SetTitle("Ring 2: Beginning of Run");
92     
93     TGraph* hr2  = new TGraph(21, z2, h3);
94     hr2->SetLineColor(2);
95
96     TGraph* hr3  = new TGraph(21, z2, h5);
97     hr3->SetLineColor(4);
98
99     hr1->Draw("AL");
100     hr2->Draw("L");
101 //    hr3->Draw("L");
102     text();
103
104 //
105 // 
106     TCanvas *c4 = new TCanvas("c4","Gas Pressure Beam 2", 200, 10, 700, 500);
107     gPad->SetLogy();
108     
109     TGraph* hr4  = new TGraph(21, z2, h2);
110     hr4->SetMaximum(1e17);
111     hr4->SetMinimum(1e11);    
112     hr4->SetLineColor(1);
113     hr4->SetTitle("Ring 2");
114     
115     TGraph* hr5  = new TGraph(21, z2, h4);
116     hr5->SetLineColor(2);
117
118     TGraph* hr6  = new TGraph(21, z2, h5);
119     hr6->SetLineColor(4);
120
121     hr4->Draw("AL");
122     hr5->Draw("L");
123     hr6->Draw("L");
124     text();
125 //
126 // One example of interaction rate calculation
127 //
128 // 3rd year, full intensity
129 //
130     Float_t r[20];
131     const Float_t crossSection = 0.094e-28;   //  m^2
132     const Float_t pFlux        = 1.e11/25.e-9; // 1/s
133     
134     for (Int_t i= 0; i < 20; i++) 
135     {
136         r[i] =  g5[i] * crossSection * pFlux; // 1/m/s 
137     }
138
139     TCanvas *c5 = new TCanvas("c5","Interaction Rate Beam 1, 3rd year", 200, 10, 700, 500);
140     gPad->SetLogy();
141     
142     TGraph* rr1  = new TGraph(20, z1, r);
143     rr1->SetMaximum(1e6);
144     rr1->SetMinimum(1e1);    
145     rr1->SetLineColor(4);
146     rr1->SetTitle("Ring 1");
147     rr1->Draw("AL");
148 }
149
150 void text()
151 {
152     
153     TPave *pave = new TPave(194.619,11.0495,293.103,12.6485,4,"br");
154     pave->SetFillColor(18);
155     pave->Draw();
156     TLine *line = new TLine(199.702,11.8109,232.106,11.8109);
157     line->SetLineColor(2);
158     line->Draw();
159     line = new TLine(198.431,12.2868,232.106,12.2868);
160     line->Draw();
161     line = new TLine(199.066,11.3731,232.742,11.3731);
162     line->SetLineColor(4);
163     line->Draw();
164     line = new TLine(215.586,13.3338,215.586,13.3147);
165     line->Draw();
166     tex = new TLatex(239.096,11.6777,"2nd year");
167     tex->SetTextSize(0.05);
168     tex->SetLineWidth(2);
169     tex->Draw();
170     tex = new TLatex(236.554,12.1536,"1st year");
171     tex->SetTextSize(0.0507614);
172     tex->SetLineWidth(2);
173     tex->Draw();
174     tex = new TLatex(239.096,11.2589,"3rd year");
175     tex->SetTextSize(0.05);
176     tex->SetLineWidth(2);
177     tex->Draw();
178     c1->Modified();
179     c1->cd();
180 }