Implemented a new version of cluster (with its store and iterator):

[u/mrichter/AliRoot.git] / PYTHIA6 / AliPythia.h

diff --git a/PYTHIA6/AliPythia.h b/PYTHIA6/AliPythia.h
index c812b730b1cec3d65c825b524fcf96178bc4fecb..3b4645de1e1c57d62c1f50ec26afd1449a51f202 100644 (file)
--- a/PYTHIA6/AliPythia.h
+++ b/PYTHIA6/AliPythia.h
@@ -8,6 +8,7 @@
 #include <TPythia6.h>
 #include <AliRndm.h>
 #include <AliStructFuncType.h>
+
 typedef enum
 {kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
  kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
@@ -15,26 +16,10 @@ typedef enum
  kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyDPlusStrangepPbMNR, kPyBeautypPbMNR,
  kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR, 
  kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyW, kPyZ, kPyMbMSEL1,
- kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
+ kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn,
+ kPyLhwgMb,kPyMbDefault}
 Process_t;
-/*
-typedef enum
-{
-    kDOSet1     = 1006,
-    kGRVLO      = 5005,
-    kGRVHO      = 5006,
-    kMRSDminus  = 3031,
-    kMRSD0      = 3030,
-    kMRSG       = 3041,
-    kCTEQ2pM    = 4024,
-    kCTEQ4L     = 4032,
-    kCTEQ4M     = 4034,
-    kMRSTcgLO   = 3072,
-    kCTEQ5L     = 4046,
-    kGRVLO98    = 5012
-}
-StrucFunc_t;
-*/
+
 class AliFastGlauber;
 class AliQuenchingWeights;
 
@@ -59,12 +44,15 @@ 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 k, Int_t iECMethod);
+    virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod, Float_t zmax = 0.97, Int_t ngmax = 30);
     virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
     virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
     // return instance of the singleton
     static  AliPythia* Instance();
     virtual void Quench();
+    // Assignment Operator
+    AliPythia & operator=(const AliPythia & rhs);
+    void Copy(TObject&) const;
  protected:
     Process_t             fProcess;           // Process type
     Float_t               fEcms;              // Centre of mass energy
@@ -74,11 +62,14 @@ class AliPythia : public TPythia6, public AliRndm
     Double_t              fZQuench[4];        //  ! Quenching fractions for this even
     Double_t              fXJet;              //  ! Jet production point X
     Double_t              fYJet;              //  ! Jet production point Y
+    Int_t                 fNGmax;             //    Maximum number of radiated gluons in quenching
+    Float_t               fZmax;              //    Maximum energy loss in quenching
     AliFastGlauber*       fGlauber;           //  ! The Glauber model
     AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
     static AliPythia*     fgAliPythia;        // Pointer to single instance
  private: 
     AliPythia();
+    AliPythia(const AliPythia& pythia);
     void ConfigHeavyFlavor();
     void AtlasTuning();
     ClassDef(AliPythia,1) //ALICE UI to PYTHIA