From c0382a5718a3d891339f4f9c113ba22846c3947a Mon Sep 17 00:00:00 2001 From: morsch Date: Wed, 30 Jan 2002 16:13:11 +0000 Subject: [PATCH] LHC gas pressures for LSS2 Optics Version 6.3 (LHC Project Note 274) --- LHC/PlotPressure.C | 158 ++++++++++++++++++++++++++++++++++++++++++++ LHC/gasPressure.dat | 41 ++++++++++++ 2 files changed, 199 insertions(+) create mode 100644 LHC/PlotPressure.C create mode 100644 LHC/gasPressure.dat diff --git a/LHC/PlotPressure.C b/LHC/PlotPressure.C new file mode 100644 index 00000000000..b2992e1ad6b --- /dev/null +++ b/LHC/PlotPressure.C @@ -0,0 +1,158 @@ +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(); +} diff --git a/LHC/gasPressure.dat b/LHC/gasPressure.dat new file mode 100644 index 00000000000..be388385b73 --- /dev/null +++ b/LHC/gasPressure.dat @@ -0,0 +1,41 @@ +0.83 1.5e14 2e13 2e15 2e14 1e14 +7.70 2e13 3e12 5e13 8e12 6e12 +1.40 2e13 3e12 3e15 3e14 1.5e14 +12.58 2e13 3e12 3e13 5e12 3e12 +1.90 2e13 3e12 3e15 3e14 1.5e14 +8.40 2e13 3e12 5e13 8e12 6e12 +3.23 2e13 3e12 3e15 3e14 1.5e14 +11.36 2e13 3e12 3e15 3e14 1.5e14 +5.0 1e16 1e16 1e16 1e16 1e16 +47.07 1e12 1e12 1e12 1e12 1e12 +1.50 3e15 1e14 1.5e14 5e13 1e12 +4.0 1e12 1e12 7e16 6e14 5e12 +11.67 1e12 1e12 6e15 3e14 1.5e14 +12.50 3e12 9e11 6e14 4e13 2e13 +16.75 1e15 1e12 6e16 6e14 5e12 +12.99 3e12 9e11 6e14 4e13 2e13 +60.96 3e15 1e13 6e16 6e14 5e12 +10.36 3e13 2e12 4e14 4e13 2e13 +14.7 1e14 7e13 6e14 3e14 1.5e14 +9.00 1.5e14 1.2e13 1.8e13 1.6e13 2e13 +0.83 1.5e14 2e13 2e15 2e14 1e14 +7.70 2e13 3e12 5e13 8e12 6e12 +1.40 2e13 3e12 3e15 3e14 1.5e14 +12.58 2e13 3e12 3e13 5e12 3e12 +1.90 2e13 3e12 3e15 3e14 1.5e14 +8.40 2e13 3e12 5e13 8e12 6e12 +3.23 2e13 3e12 3e15 3e14 1.5e14 +11.36 2e13 3e12 3e15 3e14 1.5e14 +5.0 1e16 1e16 1e16 1e16 1e16 +47.07 1e12 1e12 1e12 1e12 1e12 +1.50 3e15 1e14 1.5e14 5e13 1e12 +4.0 1e12 1e12 7e16 6e14 5e12 +11.67 1e12 1e12 6e15 3e14 1.5e14 +12.50 3e12 9e11 6e14 4e13 2e13 +16.75 2e15 1e12 6e16 6e14 5e12 +12.99 3e12 1.2e12 5e14 4e13 2e13 +60.96 2e14 1e13 6e16 6e14 5e12 +10.36 2e12 1e12 5e14 4e13 2e13 +14.7 1e14 1e12 6e14 3e14 1.5e14 +9.00 1e14 1e12 1.8e13 1.6e13 2e13 +170.4 2e14 1e13 7e14 3e13 4e12 -- 2.31.1