AliFastGlauber is a singleton.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 Nov 2008 15:34:55 +0000 (15:34 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 Nov 2008 15:34:55 +0000 (15:34 +0000)
EVGEN/AliGenMUONCocktail.cxx
FASTSIM/AliFastGlauber.cxx
FASTSIM/AliFastGlauber.h
PYTHIA6/AliPythia.cxx
PYTHIA6/AliPythia6.cxx

index 1b8257e..dd263e5 100644 (file)
@@ -82,7 +82,7 @@ void AliGenMUONCocktail::Init()
   Double_t sigmaReaction =   0.072;   // MinBias NN cross section for PbPb LHC energies  http://arxiv.org/pdf/nucl-ex/0302016
 
   // Initialising Fast Glauber object
-  fFastGlauber = new AliFastGlauber();
+  fFastGlauber = AliFastGlauber::Instance();
   fFastGlauber->SetPbPbLHC();
   fFastGlauber->SetNNCrossSection(sigmaReaction*1000.); //Expected NN cross-section in mb at LHC with diffractive part http://arxiv.org/pdf/nucl-ex/0302016 )
   fFastGlauber->Init(1);
index 020b456..eed4aeb 100644 (file)
@@ -85,7 +85,8 @@ TF1*    AliFastGlauber::fgWParticipants  = NULL;
 TF2*    AliFastGlauber::fgWAlmondCurrent = NULL;    
 TF2*    AliFastGlauber::fgWAlmondFixedB[40]; 
 const Int_t AliFastGlauber::fgkMCInts = 100000;
-Int_t AliFastGlauber::fgCounter = 0;       
+AliFastGlauber* AliFastGlauber::fgGlauber = NULL;
+
 
 AliFastGlauber::AliFastGlauber(): 
     fWSr0(0.),
@@ -101,10 +102,6 @@ AliFastGlauber::AliFastGlauber():
     fName()     
 {
   //  Default Constructor 
-  fgCounter++;
-  if(fgCounter>1)
-    Error("AliFastGlauber","More than one instance (%d) is not supported, check your code!",fgCounter);
-
   //  Defaults for Pb
   SetMaxImpact();
   SetLengthDefinition();
@@ -129,10 +126,20 @@ AliFastGlauber::AliFastGlauber(const AliFastGlauber & gl)
     gl.Copy(*this);
 }
 
+AliFastGlauber* AliFastGlauber::Instance()
+{ 
+// Set random number generator 
+    if (fgGlauber) {
+       return fgGlauber;
+    } else {
+       fgGlauber = new AliFastGlauber();
+       return fgGlauber;
+    }
+}
+
 AliFastGlauber::~AliFastGlauber()
 {
 // Destructor
-  fgCounter--;
   for(Int_t k=0; k<40; k++) delete fgWAlmondFixedB[k];
 }
 
index 0036b27..479a6af 100644 (file)
@@ -18,8 +18,7 @@ class TF2;
 
 class AliFastGlauber : public TObject {
  public:
-    AliFastGlauber();
-    AliFastGlauber(const AliFastGlauber& glauber);
+    static AliFastGlauber* Instance();
     virtual ~AliFastGlauber();
     void Init(Int_t mode = 0);
 
@@ -178,11 +177,15 @@ class AliFastGlauber : public TObject {
     static Double_t WEnergyDensity (Double_t *xx, Double_t *par);
 
     void Reset() const;
+ private:
+    AliFastGlauber();
+    AliFastGlauber(const AliFastGlauber& glauber);
 
-    static Float_t fgBMax;           // Maximum Impact Parameter
-    static Int_t fgCounter;          // Counter to protect double instantiation
-    static const Int_t fgkMCInts;    // Number of MC integrations
-
+    static Float_t         fgBMax;       // Maximum Impact Parameter
+    static const Int_t     fgkMCInts;    // Number of MC integrations
+    static AliFastGlauber* fgGlauber;    // Singleton instance
+     
+    
     static TF1*    fgWSb;            // Wood-Saxon Function (b)
     static TF2*    fgWSbz;           // Wood-Saxon Function (b, z)
     static TF1*    fgWSz;            // Wood-Saxon Function (b = b0, z)
index 34c3bf2..2a532d3 100644 (file)
@@ -692,7 +692,7 @@ void AliPythia::InitQuenching(Float_t cMin, Float_t cMax, Float_t k, Int_t iECMe
 // (2) The nuclear geometry using the Glauber Model
 //     
     
-    fGlauber = new AliFastGlauber();
+    fGlauber = AliFastGlauber::Instance();
     fGlauber->Init(2);
     fGlauber->SetCentralityClass(cMin, cMax); 
 
index b0b126d..f8c761f 100644 (file)
@@ -704,7 +704,7 @@ void AliPythia6::InitQuenching(Float_t cMin, Float_t cMax, Float_t k, Int_t iECM
 // (2) The nuclear geometry using the Glauber Model
 //     
     
-    fGlauber = new AliFastGlauber();
+    fGlauber = AliFastGlauber::Instance();
     fGlauber->Init(2);
     fGlauber->SetCentralityClass(cMin, cMax);