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 typedef enum
  12 {kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
  13  kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
  14  kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR,
  15  kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR,
  16  kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW,
  17  kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
  18 Process_t;
  19 /*
  20 typedef enum
  21 {
  22     kDOSet1     = 1006,
  23     kGRVLO      = 5005,
  24     kGRVHO      = 5006,
  25     kMRSDminus  = 3031,
  26     kMRSD0      = 3030,
  27     kMRSG       = 3041,
  28     kCTEQ2pM    = 4024,
  29     kCTEQ4L     = 4032,
  30     kCTEQ4M     = 4034,
  31     kMRSTcgLO   = 3072,
  32     kCTEQ5L     = 4046,
  33     kGRVLO98    = 5012
  34 }
  35 StrucFunc_t;
  36 */
  37 class AliFastGlauber;
  38 class AliQuenchingWeights;
  39
  40 class AliPythia : public TPythia6, public AliRndm
  41 {
  42
  43  public:
  44     virtual ~AliPythia(){;}
  45     // convert to compressed code and print result (for debugging only)
  46     virtual Int_t CheckedLuComp(Int_t kf);
  47     // Pythia initialisation for selected processes
  48     virtual void ProcInit
  49         (Process_t process, Float_t energy, StrucFunc_t strucfunc);
  50     // treat protons as inside nuclei
  51     virtual void  SetNuclei(Int_t a1, Int_t a2);
  52     // Print particle properties
  53     virtual void PrintParticles();
  54     virtual void ResetDecayTable();
  55     virtual void SetDecayTable();
  56     virtual void Pycell(Int_t& nclus);
  57     virtual void Pyclus(Int_t& nclus);
  58     virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
  59     virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
  60     virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod);
  61     virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
  62     virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
  63     // return instance of the singleton
  64     static  AliPythia* Instance();
  65     virtual void Quench();
  66  protected:
  67     Process_t             fProcess;           // Process type
  68     Float_t               fEcms;              // Centre of mass energy
  69     StrucFunc_t           fStrucFunc;         // Structure function
  70     Int_t                 fDefMDCY[501];      //  ! Default decay switches per particle
  71     Int_t                 fDefMDME[2001];     //  ! Default decay switches per mode
  72     Double_t              fZQuench[4];        //  ! Quenching fractions for this even
  73     Double_t              fXJet;              //  ! Jet production point X
  74     Double_t              fYJet;              //  ! Jet production point Y
  75     AliFastGlauber*       fGlauber;           //  ! The Glauber model
  76     AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
  77     static AliPythia*     fgAliPythia;        // Pointer to single instance
  78  private:
  79     AliPythia();
  80
  81     ClassDef(AliPythia,1) //ALICE UI to PYTHIA
  82 };
  83
  84 #endif
  85
  86
  87
  88
  89