From f5e004f6c0d67f3a2775da09e3ed1bc7a94f83d5 Mon Sep 17 00:00:00 2001
From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 30 Jan 2002 16:38:36 +0000
Subject: [PATCH] Example for interaction rate calculation added.

---
 LHC/PlotPressure.C | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)

diff --git a/LHC/PlotPressure.C b/LHC/PlotPressure.C
index b2992e1ad6b..7bca21c4e44 100644
--- a/LHC/PlotPressure.C
+++ b/LHC/PlotPressure.C
@@ -122,7 +122,29 @@ void PlotPressure()
     hr5->Draw("L");
     hr6->Draw("L");
     text();
+//
+// One example of interaction rate calculation
+//
+// 3rd year, full intensity
+//
+    Float_t r[20];
+    const Float_t crossSection = 0.094e-28;   //  m^2
+    const Float_t pFlux        = 1.e11/25.e-9; // 1/s
+    
+    for (Int_t i= 0; i < 20; i++) 
+    {
+	r[i] =  g5[i] * crossSection * pFlux; // 1/m/s 
+    }
+
+    TCanvas *c5 = new TCanvas("c5","Interaction Rate Beam 1, 3rd year", 200, 10, 700, 500);
+    gPad->SetLogy();
     
+    TGraph* rr1  = new TGraph(20, z1, r);
+    rr1->SetMaximum(1e6);
+    rr1->SetMinimum(1e1);    
+    rr1->SetLineColor(4);
+    rr1->SetTitle("Ring 1");
+    rr1->Draw("AL");
 }
 
 void text()
-- 
2.31.1