PYTHIA6/AliPythia.h

   1 #ifndef ALIPYTHIA_H
   2 #define ALIPYTHIA_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 #include <TPythia6.h>
   9 #include <AliRndm.h>
  10 #include <AliStructFuncType.h>
  11
  12 typedef enum
  13 {kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
  14  kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
  15  kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyDPlusStrangePbPbMNR, kPyBeautyPbPbMNR,
  16  kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyDPlusStrangepPbMNR, kPyBeautypPbMNR,
  17  kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR,
  18  kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyW, kPyZ, kPyMbMSEL1,
  19  kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn,
  20  kPyLhwgMb,kPyMbDefault}
  21 Process_t;
  22
  23 class AliFastGlauber;
  24 class AliQuenchingWeights;
  25
  26 class AliPythia : public TPythia6, public AliRndm
  27 {
  28
  29  public:
  30     virtual ~AliPythia(){;}
  31     // convert to compressed code and print result (for debugging only)
  32     virtual Int_t CheckedLuComp(Int_t kf);
  33     // Pythia initialisation for selected processes
  34     virtual void ProcInit
  35         (Process_t process, Float_t energy, StrucFunc_t strucfunc);
  36     // treat protons as inside nuclei
  37     virtual void  SetNuclei(Int_t a1, Int_t a2);
  38     // Print particle properties
  39     virtual void PrintParticles();
  40     virtual void ResetDecayTable();
  41     virtual void SetDecayTable();
  42     virtual void Pyevnw();
  43     virtual void Pycell(Int_t& nclus);
  44     virtual void Pyclus(Int_t& nclus);
  45     virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
  46     virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
  47     virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod, Float_t zmax = 0.97, Int_t ngmax = 30);
  48     virtual void SetPyquenParameters(Double_t t0, Double_t tau0, Int_t nf, Int_t iengl, Int_t iangl);
  49     virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
  50     virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
  51     // return instance of the singleton
  52     static  AliPythia* Instance();
  53     virtual void Quench();
  54     // Assignment Operator
  55     AliPythia & operator=(const AliPythia & rhs);
  56     void Copy(TObject&) const;
  57  protected:
  58     Process_t             fProcess;           // Process type
  59     Float_t               fEcms;              // Centre of mass energy
  60     StrucFunc_t           fStrucFunc;         // Structure function
  61     Int_t                 fDefMDCY[501];      //  ! Default decay switches per particle
  62     Int_t                 fDefMDME[2001];     //  ! Default decay switches per mode
  63     Double_t              fZQuench[4];        //  ! Quenching fractions for this even
  64     Double_t              fXJet;              //  ! Jet production point X
  65     Double_t              fYJet;              //  ! Jet production point Y
  66     Int_t                 fNGmax;             //    Maximum number of radiated gluons in quenching
  67     Float_t               fZmax;              //    Maximum energy loss in quenching
  68     AliFastGlauber*       fGlauber;           //  ! The Glauber model
  69     AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
  70     static AliPythia*     fgAliPythia;        // Pointer to single instance
  71  private:
  72     AliPythia();
  73     AliPythia(const AliPythia& pythia);
  74     void ConfigHeavyFlavor();
  75     void AtlasTuning();
  76     ClassDef(AliPythia,1) //ALICE UI to PYTHIA
  77 };
  78
  79 #endif
  80
  81
  82
  83
  84