+void PlotPressure()
+{
+ FILE* file;
+ file=fopen("gasPressure.dat","r");
+ Float_t z1[20];
+ Float_t g1[20], g2[20], g3[20], g4[20], g5[20];
+ Float_t z2[21];
+ Float_t h1[21], h2[21], h3[21], h4[21], h5[21];
+
+ char c[45];
+ Float_t z;
+
+ for (Int_t i = 0; i < 20; i++)
+ {
+ fscanf(file, "%f %f %f %f %f %f", &z,
+ &g1[i], &g2[i], &g3[i], &g4[i], &g5[i]);
+// printf("%d %f %f %f %f %f %f \n", i, z1[i],
+// g1[i][0], g1[i][1], g1[i][2], g1[i][3], g1[i][4]);
+ if (i > 0) {
+ z1[i] = z1[i-1] + z;
+ } else {
+ z1[i] = 20.;
+ }
+ }
+
+
+ for (Int_t i = 0; i < 21; i++)
+ {
+ fscanf(file, "%f %f %f %f %f %f", &z,
+ &h1[i], &h2[i], &h3[i], &h4[i], &h5[i]);
+ if (i > 0) {
+ z2[i] = z2[i-1] + z;
+ } else {
+ z2[i] = 20.;
+ }
+ }
+
+//
+//
+ TCanvas *c1 = new TCanvas("c1","Gas Pressure Beam 1", 200, 10, 700, 500);
+ gPad->SetLogy();
+
+ TGraph* gr1 = new TGraph(20, z1, g1);
+ gr1->SetMaximum(1e17);
+ gr1->SetMinimum(1e11);
+ gr1->SetLineColor(1);
+ gr1->SetTitle("Ring 1: Beginning of Run");
+
+ TGraph* gr2 = new TGraph(20, z1, g3);
+ gr2->SetLineColor(2);
+
+ TGraph* gr3 = new TGraph(20, z1, g5);
+ gr3->SetLineColor(4);
+
+ gr1->Draw("AL");
+ gr2->Draw("L");
+// gr3->Draw("L");
+ text();
+
+//
+//
+ TCanvas *c2 = new TCanvas("c2","Gas Pressure Beam 1", 200, 10, 700, 500);
+ gPad->SetLogy();
+
+ TGraph* gr4 = new TGraph(20, z1, g2);
+ gr4->SetMaximum(1e17);
+ gr4->SetMinimum(1e11);
+ gr4->SetLineColor(1);
+ gr4->SetTitle("Ring 1");
+
+ TGraph* gr5 = new TGraph(20, z1, g4);
+ gr5->SetLineColor(2);
+
+ TGraph* gr6 = new TGraph(20, z1, g5);
+ gr6->SetLineColor(4);
+
+ gr4->Draw("AL");
+ gr5->Draw("L");
+ gr6->Draw("L");
+ text();
+
+//
+//
+ TCanvas *c3 = new TCanvas("c3","Gas Pressure Beam 2", 200, 10, 700, 500);
+ gPad->SetLogy();
+
+ TGraph* hr1 = new TGraph(21, z2, h1);
+ hr1->SetMaximum(1e17);
+ hr1->SetMinimum(1e11);
+ hr1->SetLineColor(1);
+ hr1->SetTitle("Ring 2: Beginning of Run");
+
+ TGraph* hr2 = new TGraph(21, z2, h3);
+ hr2->SetLineColor(2);
+
+ TGraph* hr3 = new TGraph(21, z2, h5);
+ hr3->SetLineColor(4);
+
+ hr1->Draw("AL");
+ hr2->Draw("L");
+// hr3->Draw("L");
+ text();
+
+//
+//
+ TCanvas *c4 = new TCanvas("c4","Gas Pressure Beam 2", 200, 10, 700, 500);
+ gPad->SetLogy();
+
+ TGraph* hr4 = new TGraph(21, z2, h2);
+ hr4->SetMaximum(1e17);
+ hr4->SetMinimum(1e11);
+ hr4->SetLineColor(1);
+ hr4->SetTitle("Ring 2");
+
+ TGraph* hr5 = new TGraph(21, z2, h4);
+ hr5->SetLineColor(2);
+
+ TGraph* hr6 = new TGraph(21, z2, h5);
+ hr6->SetLineColor(4);
+
+ hr4->Draw("AL");
+ hr5->Draw("L");
+ hr6->Draw("L");
+ text();
+
+}
+
+void text()
+{
+
+ TPave *pave = new TPave(194.619,11.0495,293.103,12.6485,4,"br");
+ pave->SetFillColor(18);
+ pave->Draw();
+ TLine *line = new TLine(199.702,11.8109,232.106,11.8109);
+ line->SetLineColor(2);
+ line->Draw();
+ line = new TLine(198.431,12.2868,232.106,12.2868);
+ line->Draw();
+ line = new TLine(199.066,11.3731,232.742,11.3731);
+ line->SetLineColor(4);
+ line->Draw();
+ line = new TLine(215.586,13.3338,215.586,13.3147);
+ line->Draw();
+ tex = new TLatex(239.096,11.6777,"2nd year");
+ tex->SetTextSize(0.05);
+ tex->SetLineWidth(2);
+ tex->Draw();
+ tex = new TLatex(236.554,12.1536,"1st year");
+ tex->SetTextSize(0.0507614);
+ tex->SetLineWidth(2);
+ tex->Draw();
+ tex = new TLatex(239.096,11.2589,"3rd year");
+ tex->SetTextSize(0.05);
+ tex->SetLineWidth(2);
+ tex->Draw();
+ c1->Modified();
+ c1->cd();
+}
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff