82b9abd0ea48b00300459816600a492e773ad98f
[u/mrichter/AliRoot.git] / 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, kPyZ,
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 Pyevnw();
57     virtual void Pycell(Int_t& nclus);
58     virtual void Pyclus(Int_t& nclus);
59     virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
60     virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
61     virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod);
62     virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
63     virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
64     // return instance of the singleton
65     static  AliPythia* Instance();
66     virtual void Quench();
67  protected:
68     Process_t             fProcess;           // Process type
69     Float_t               fEcms;              // Centre of mass energy
70     StrucFunc_t           fStrucFunc;         // Structure function
71     Int_t                 fDefMDCY[501];      //  ! Default decay switches per particle
72     Int_t                 fDefMDME[2001];     //  ! Default decay switches per mode
73     Double_t              fZQuench[4];        //  ! Quenching fractions for this even
74     Double_t              fXJet;              //  ! Jet production point X
75     Double_t              fYJet;              //  ! Jet production point Y
76     AliFastGlauber*       fGlauber;           //  ! The Glauber model
77     AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
78     static AliPythia*     fgAliPythia;        // Pointer to single instance
79  private: 
80     AliPythia();
81     void ConfigHeavyFlavor();
82     
83     ClassDef(AliPythia,1) //ALICE UI to PYTHIA
84 };
85
86 #endif
87
88
89
90
91