Protection in case of wrong time stemp improved
[u/mrichter/AliRoot.git] / TUHKMgen / AliGenUHKM.h
CommitLineData
7b7936e9 1// This class provides an interface between the HYDJET++ Monte-Carlo model
2// and AliRoot (by inheriting from the AliGenMC class).
3// This class uses the TUHKMgen class (which inherits from TGenerator) to
4// transmit parameters and receive output from the model.
5
b1c2e580 6#ifndef ALIGENUHKM_H
7#define ALIGENUHKM_H
8
9#include "AliGenMC.h"
10#include <TString.h>
b1c2e580 11#include "TUHKMgen.h"
3fa37a65 12#ifndef INITIALSTATEHYDJET_H
b1c2e580 13#include "InitialStateHydjet.h"
14#endif
7b7936e9 15#include "TParticle.h"
b1c2e580 16
17#include <string>
18using namespace std;
19
b1c2e580 20
21class AliGenUHKM : public AliGenMC
22{
7b7936e9 23 protected:
24 Int_t fTrials; // Number of trials
25 TUHKMgen *fUHKMgen; // UHKM
26
27 InitialParamsHydjet_t fHydjetParams; // list of parameters for the initial state
28 // details for the PDG database
29 Char_t fParticleFilename[256]; // particle list filename
30 Char_t fDecayFilename[256]; // decay table filename
31 Int_t fStableFlagPDG[500]; // array of PDG codes flagged
32 Bool_t fStableFlagStatus[500]; // array of decay flag status
33 Int_t fStableFlagged; // number of toggled decay flags
34
35 void SetAllParameters();
36 void CheckPDGTable();
37
38 private:
39 void Copy(TObject &rhs) const;
40 AliGenUHKM(const AliGenUHKM&);
41 AliGenUHKM & operator = (const AliGenUHKM &);
42
b1c2e580 43 public:
44 AliGenUHKM();
45 AliGenUHKM(Int_t npart);
46
47 virtual ~AliGenUHKM();
48 virtual void Generate();
49 virtual void Init();
50 // virtual void AddHeader(AliGenEventHeader* header);
7b7936e9 51
b1c2e580 52 // Setters
53 // set reasonable default parameters suited for central Au+Au collisions at RHIC(200GeV)
54 void SetAllParametersRHIC();
55 // set reasonable default parameters suited for central Pb+Pb collisions at LHC(5.5TeV)
56 void SetAllParametersLHC();
7b7936e9 57
b1c2e580 58 void SetEcms(Double_t value) {fHydjetParams.fSqrtS = value;} // CMS energy per nucleon [GeV] (<2.24 given temperature and ch pot are used)
59 void SetAw(Double_t value) {fHydjetParams.fAw = value;} // nuclei mass number
60 void SetBmin(Double_t value) {fHydjetParams.fBmin = value;} // Minimum impact parameter
61 void SetBmax(Double_t value) {fHydjetParams.fBmax = value;} // Maximum impact parameter
62 void SetChFrzTemperature(Double_t value) {fHydjetParams.fT = value;} // Temperature for the chemical freezeout [GeV]
63 void SetMuB(Double_t value) {fHydjetParams.fMuB = value;} // Baryonic chemical potential [GeV]
64 void SetMuS(Double_t value) {fHydjetParams.fMuS = value;} // Strangeness chemical potential [GeV]
65 void SetMuQ(Double_t value) {fHydjetParams.fMuI3 = value;} // Isospin chemical potential [GeV]
66 void SetThFrzTemperature(Double_t value) {fHydjetParams.fThFO = value;} // Temperature for the thermal freezeout [GeV]
67 void SetMuPionThermal(Double_t value) {fHydjetParams.fMu_th_pip = value;} // Chemical potential for pi+ at thermal freezeout [GeV]
68 void SetSeed(Int_t value) {fHydjetParams.fSeed = value;} //parameter to set the random nuber seed (=0 the current time is used
69 //to set the random generator seed, !=0 the value fSeed is
70 //used to set the random generator seed and then the state of random
71 //number generator in PYTHIA MRPY(1)=fSeed
72 void SetTauB(Double_t value) {fHydjetParams.fTau = value;} // Proper time for the freeze-out hyper surface [fm/c]
73 void SetSigmaTau(Double_t value) {fHydjetParams.fSigmaTau = value;} // Standard deviation for the proper time (emission duration) [fm/c]
74 void SetRmaxB(Double_t value) {fHydjetParams.fR = value;} // Maximal transverse radius [fm]
75 void SetYlMax(Double_t value) {fHydjetParams.fYlmax = value;} // Maximal fireball longitudinal rapidity
76 void SetEtaRMax(Double_t value) {fHydjetParams.fUmax = value;} // Maximal transverse velocity
77 void SetMomAsymmPar(Double_t value) {fHydjetParams.fDelta = value;} // Momentum asymmetry parameter
78 void SetCoordAsymmPar(Double_t value) {fHydjetParams.fEpsilon = value;} // Coordinate asymmetry parameter
79
b1c2e580 80 void SetEtaType(Int_t value) {fHydjetParams.fEtaType = value;} // flag to choose rapidity distribution, if fEtaType<=0,
81 //then uniform rapidity distribution in [-fYlmax,fYlmax] if fEtaType>0,
82 //then Gaussian with dispertion = fYlmax
83 void SetGammaS(Double_t value) {fHydjetParams.fCorrS = value;} // Strangeness suppression parameter (if gamma_s<=0 then it will be calculated)
84
85 //PYQUEN parameters
86 void SetPyquenNhsel(Int_t value) {fHydjetParams.fNhsel = value;} // Flag to choose the type of event to be generated
87 // fNhsel = 0 --> UHKM fireball, no jets
88 // fNhsel = 1 --> UHKM fireball, jets with no quenching
89 // fNhsel = 2 --> UHKM fireball, jets with quenching
90 // fNhsel = 3 --> no UHKM fireball, jets with no quenching
91 // fNhsel = 4 --> no UHKM fireball, jets with quenching
92 void SetPyquenShad(Int_t value) {fHydjetParams.fIshad = value;}//flag to switch on/off impact parameter dependent nuclear
93 // shadowing for gluons and light sea quarks (u,d,s) (0: shadowing off,
94 // 1: shadowing on for fAw=207, 197, 110, 40, default: 1
95 void SetPyquenPtmin(Double_t value) {fHydjetParams.fPtmin = value;} // Pyquen input parameter for minimum Pt of parton-parton scattering (5GeV<pt<500GeV)
96 void SetPyquenT0(Double_t value) {fHydjetParams.fT0 = value;} //proper QGP formation tempereture
97 void SetPyquenTau0(Double_t value) {fHydjetParams.fTau0 = value;} //proper QGP formation time in fm/c (0.01<fTau0<10)
98 void SetPyquenNf(Int_t value) {fHydjetParams.fNf = value;} //number of active quark flavours N_f in QGP fNf=0, 1,2 or 3
99 void SetPyquenIenglu(Int_t value) {fHydjetParams.fIenglu = value;} // flag to fix type of in-medium partonic energy loss
100 //(0: radiative and collisional loss, 1: radiative loss only, 2:
101 //collisional loss only) (default: 0);
102 void SetPyquenIanglu(Int_t value) {fHydjetParams.fIanglu = value;} //flag to fix type of angular distribution of in-medium emitted
103 //gluons (0: small-angular, 1: wide-angular, 2:collinear) (default: 0).
104
105
7b7936e9 106 void SetPDGParticleFile(const Char_t *name) {strcpy(fParticleFilename, name);}//Set the filename containing the particle PDG info
107 void SetPDGDecayFile(const Char_t *name) {strcpy(fDecayFilename, name);} //Set the filename containing the PDG decay channels info
b1c2e580 108 void SetPDGParticleStable(Int_t pdg, Bool_t value) { // Turn on/off the decay flag for a PDG particle
109 fStableFlagPDG[fStableFlagged] = pdg;
110 fStableFlagStatus[fStableFlagged++] = value;
111 }
b1c2e580 112
113 // Getters
7b7936e9 114 Double_t GetEcms() const {return fHydjetParams.fSqrtS;}
115 Double_t GetAw() const {return fHydjetParams.fAw;}
116 Double_t GetBmin() const {return fHydjetParams.fBmin;}
117 Double_t GetBmax() const {return fHydjetParams.fBmax;}
118 Double_t GetChFrzTemperature() const {return fHydjetParams.fT;}
119 Double_t GetMuB() const {return fHydjetParams.fMuB;}
120 Double_t GetMuS() const {return fHydjetParams.fMuS;}
121 Double_t GetMuQ() const {return fHydjetParams.fMuI3;}
122 Double_t GetThFrzTemperature() const {return fHydjetParams.fThFO;}
123 Double_t GetMuPionThermal() const {return fHydjetParams.fMu_th_pip;}
124 Int_t GetSeed() const {return fHydjetParams.fSeed;}
125 Double_t GetTauB() const {return fHydjetParams.fTau;}
126 Double_t GetSigmaTau() const {return fHydjetParams.fSigmaTau;}
127 Double_t GetRmaxB() const {return fHydjetParams.fR;}
128 Double_t GetYlMax() const {return fHydjetParams.fYlmax;}
129 Double_t GetEtaRMax() const {return fHydjetParams.fUmax;}
130 Double_t GetMomAsymmPar() const {return fHydjetParams.fDelta;}
131 Double_t GetCoordAsymmPar() const {return fHydjetParams.fEpsilon;}
132 Int_t GetEtaType() const {return fHydjetParams.fEtaType;}
133 Double_t GetGammaS() const {return fHydjetParams.fCorrS;}
134 Int_t GetPyquenNhsel() const {return fHydjetParams.fNhsel;}
135 Int_t GetPyquenShad() const {return fHydjetParams.fIshad;}
136 Double_t GetPyquenPtmin() const {return fHydjetParams.fPtmin;}
137 Double_t GetPyquenT0() const {return fHydjetParams.fT0;}
138 Double_t GetPyquenTau0() const {return fHydjetParams.fTau0;}
139 Double_t GetPyquenNf() const {return fHydjetParams.fNf;}
140 Double_t GetPyquenIenglu() const {return fHydjetParams.fIenglu;}
141 Double_t GetPyquenIanglu() const {return fHydjetParams.fIanglu;}
142 const Char_t* GetPDGParticleFile() const {return fParticleFilename;}
143 const Char_t* GetPDGDecayFile() const {return fDecayFilename;}
786056a2 144
b1c2e580 145 ClassDef(AliGenUHKM, 6) // AliGenerator interface to UHKM
146};
147#endif
148
149
150
151
152