Added GetIOI1ForPythiaAndXY to get the prod. point into pythia header.
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Jun 2004 08:10:33 +0000 (08:10 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 30 Jun 2004 08:10:33 +0000 (08:10 +0000)
FASTSIM/AliFastGlauber.cxx
FASTSIM/AliFastGlauber.h

index 8eaa227..6088042 100644 (file)
@@ -1596,6 +1596,25 @@ void AliFastGlauber::GetI0I1ForPythia(Int_t n,Double_t* phi,
   return;
 }
 
+void AliFastGlauber::GetI0I1ForPythiaAndXY(Int_t n,Double_t* phi,
+                                     Double_t* integral0,Double_t* integral1,
+                                     Double_t &x,Double_t& y,
+                                     Double_t ellCut,Double_t b)
+{
+  //
+  // Returns I0 and I1 pairs for n partons with azimuthal angles phi[n] 
+  // from random b, x0, y0 and return x0,y0
+  //
+  Double_t x0,y0;
+  if(b<0.) GetRandomBHard(b);
+  GetRandomXY(x0,y0);
+  for(Int_t i=0; i<n; i++) 
+    CalculateI0I1(integral0[i],integral1[i],b,x0,y0,phi[i],ellCut);
+  x=x0;
+  y=y0;
+  return;
+}
+
 void AliFastGlauber::PlotI0I1Distr(Int_t n,Double_t ellCut,
                                   Bool_t save,Char_t *fname)
 {
index cc9986f..ff84978 100644 (file)
@@ -159,6 +159,10 @@ class AliFastGlauber : public TObject {
     void GetI0I1ForPythia(Int_t n,Double_t* phi,
                          Double_t* integral0,Double_t* integral1,
                          Double_t ellCut=20.,Double_t b=-1.);
+    void GetI0I1ForPythiaAndXY(Int_t n,Double_t* phi,
+                              Double_t* integral0,Double_t* integral1,
+                              Double_t&x, Double_t &y,
+                              Double_t ellCut=20.,Double_t b=-1.);
     void PlotI0I1Distr(Int_t n=1000,Double_t ellCut=20.,Bool_t save=kFALSE,
                       Char_t *fname="i0i1.root");
     void PlotI0I1B2BDistr(Int_t n=1000,Double_t ellCut=20.,Bool_t save=kFALSE,