Added GetIOI1ForPythiaAndXY to get the prod. point into pythia header.

diff --git a/FASTSIM/AliFastGlauber.cxx b/FASTSIM/AliFastGlauber.cxx
index 8eaa227694dc5f1ddd04bc9213555a302642d61b..6088042850e755f567e5aa6995fce79f1505ded7 100644 (file)
--- a/FASTSIM/AliFastGlauber.cxx
+++ b/FASTSIM/AliFastGlauber.cxx
@@ -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)
 {

diff --git a/FASTSIM/AliFastGlauber.h b/FASTSIM/AliFastGlauber.h
index cc9986fe38ebf41165e30ae863b9537f573d1513..ff84978a722b6168b249282d29631c97be3b4648 100644 (file)
--- a/FASTSIM/AliFastGlauber.h
+++ b/FASTSIM/AliFastGlauber.h
@@ -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,