Added NHard function.
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Oct 2004 15:47:23 +0000 (15:47 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Oct 2004 15:47:23 +0000 (15:47 +0000)
FASTSIM/AliFastGlauber.cxx
FASTSIM/AliFastGlauber.h

index 7b5c3ef..6595b42 100644 (file)
@@ -977,6 +977,19 @@ Double_t AliFastGlauber::FractionOfHardCrossSection(Double_t b1, Double_t b2) co
   return fgWSbinary->Integral(b1, b2)/fgWSbinary->Integral(0., 100.);
 }
 
+Double_t AliFastGlauber::NHard(Double_t b1, Double_t b2) const
+{
+  //
+  //  Number of binary hard collisions 
+  //  as a function of b (nucl/ex/0302016 eq. 19)
+  //
+  const Double_t kshard=HardCrossSection(b1,b2);
+  const Double_t ksgeo=CrossSection(b1,b2); 
+  if(ksgeo>0)
+    return kshard/ksgeo;
+  else return -1; 
+}
+
 Double_t AliFastGlauber::Binaries(Double_t b) const
 {
   //
index 8d817e4..7460fea 100644 (file)
@@ -100,6 +100,7 @@ class AliFastGlauber : public TObject {
     
     Double_t CrossSection(Double_t b1, Double_t b2)               const;
     Double_t HardCrossSection(Double_t b1, Double_t b2)           const;
+    Double_t NHard(Double_t b1, Double_t b2)                      const;
     Double_t FractionOfHardCrossSection(Double_t b1, Double_t b2) const;
     Double_t Binaries(Double_t b)                 const;
     Double_t GetNumberOfBinaries(Double_t b)      const;