fix for Savannah bug #80565
[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 #include "PythiaProcesses.h"
12 #include "AliOmegaDalitz.h"
13 class AliFastGlauber;
14 class AliQuenchingWeights;
15
16 class AliPythia : public TPythia6, public AliRndm
17 {
18
19  public:
20     virtual ~AliPythia(){;}
21     // convert to compressed code and print result (for debugging only)
22     virtual Int_t CheckedLuComp(Int_t kf);
23     // Pythia initialisation for selected processes
24     virtual void ProcInit
25         (Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune = -1);
26     // treat protons as inside nuclei
27     virtual void  SetNuclei(Int_t a1, Int_t a2, Int_t pdf);
28     // Print particle properties
29     virtual void PrintParticles();
30     virtual void ResetDecayTable();
31     virtual void SetDecayTable();
32     virtual void Pyevnw();
33     virtual void Pycell(Int_t& nclus);
34     virtual void Pyclus(Int_t& nclus);
35     virtual void Pyshow(Int_t ip1, Int_t ip2, Double_t qmax);
36     virtual void Pyshowq(Int_t ip1, Int_t ip2, Double_t qmax);
37     virtual void Pyrobo(Int_t imi, Int_t ima, Double_t the, Double_t phi, Double_t bex, Double_t bey, Double_t bez);
38     virtual void Pytune(Int_t itune);
39     virtual void Py2ent(Int_t idx, Int_t pdg1, Int_t pdg2, Double_t p);
40     virtual void InitQuenching(Float_t bmin, Float_t bmax, Float_t k, Int_t iECMethod, Float_t zmax = 0.97, Int_t ngmax = 30);
41     virtual void SetPyquenParameters(Double_t t0, Double_t tau0, Int_t nf, Int_t iengl, Int_t iangl);
42     virtual void Pyquen(Double_t a, Int_t ibf, Double_t b);
43     virtual void Qpygin0();
44     virtual void GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]);
45     // return instance of the singleton
46     static  AliPythia* Instance();
47     virtual void Quench();
48     void DalitzDecays();
49     // Assignment Operator
50     AliPythia & operator=(const AliPythia & rhs);
51     void Copy(TObject&) const;
52  protected:
53     Process_t             fProcess;           // Process type
54     Float_t               fEcms;              // Centre of mass energy
55     StrucFunc_t           fStrucFunc;         // Structure function
56     Int_t                 fDefMDCY[501];      //  ! Default decay switches per particle
57     Int_t                 fDefMDME[2001];     //  ! Default decay switches per mode
58     Double_t              fZQuench[4];        //  ! Quenching fractions for this even
59     Double_t              fXJet;              //  ! Jet production point X
60     Double_t              fYJet;              //  ! Jet production point Y
61     Int_t                 fNGmax;             //    Maximum number of radiated gluons in quenching
62     Float_t               fZmax;              //    Maximum energy loss in quenching
63     AliFastGlauber*       fGlauber;           //  ! The Glauber model
64     AliQuenchingWeights*  fQuenchingWeights;  //  ! The Quenching Weights model
65     Int_t                 fItune;             //  ! Pythia tune 
66     AliOmegaDalitz        fOmegaDalitz;       //  ! omega dalitz decayer
67     static AliPythia*     fgAliPythia;        // Pointer to single instance
68  private: 
69     AliPythia();
70     AliPythia(const AliPythia& pythia);
71     void ConfigHeavyFlavor();
72     void AtlasTuning();
73     void AtlasTuning_MC09();
74     ClassDef(AliPythia,1) //ALICE UI to PYTHIA
75 };
76
77 #endif
78
79
80
81
82