New quenching algorithm.
[u/mrichter/AliRoot.git] / PYTHIA6 / AliPythia.h
index c91fa98..9ce2811 100644 (file)
@@ -33,6 +33,9 @@ typedef enum
 }
 StrucFunc_t;
 */
+class AliFastGlauber;
+class AliQuenchingWeights;
+
 class AliPythia : public TPythia6, public AliRndm
 {
 
@@ -53,17 +56,20 @@ class AliPythia : public TPythia6, public AliRndm
     virtual void Pyclus(Int_t& nclus);
     virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
     virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
+    virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t qTransport, Float_t maxLength, Int_t iECMethod);
     
     // return instance of the singleton
     static  AliPythia* Instance();
     virtual void Quench();
  protected:
-    Process_t     fProcess;           // Process type
-    Float_t       fEcms;              // Centre of mass energy
-    StrucFunc_t   fStrucFunc;         // Structure function
-    Int_t         fDefMDCY[501];      //  ! Default decay switches per particle
-    Int_t         fDefMDME[2001];     //  ! Default decay switches per mode
-    static AliPythia*    fgAliPythia; // Pointer to single instance
+    Process_t             fProcess;           // Process type
+    Float_t               fEcms;              // Centre of mass energy
+    StrucFunc_t           fStrucFunc;         // Structure function
+    Int_t                 fDefMDCY[501];      //  ! Default decay switches per particle
+    Int_t                 fDefMDME[2001];     //  ! Default decay switches per mode
+    AliFastGlauber*       fGlauber;           //  ! The Glauber model
+    AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
+    static AliPythia*     fgAliPythia;        // Pointer to single instance
  private: 
     AliPythia();