TUHKMgen/UHKM/GrandCanonical.h

   1 /*
   2
   3         Nikolai Amelin, Ludmila Malinina, Timur Pocheptsov (C) JINR/Dubna
   4       amelin@sunhe.jinr.ru, malinina@sunhe.jinr.ru, pocheptsov@sunhe.jinr.ru
   5                            November. 2, 2005
   6
   7 */
   8
   9 #ifndef GRANDCANONICAL_INCLUDED
  10 #define GRANDCANONICAL_INCLUDED
  11
  12 #ifndef PARTICLE_PDG
  13 #include "ParticlePDG.h"
  14 #endif
  15 #ifndef DATABASEPDG_H
  16 #include "DatabasePDG.h"
  17 #endif
  18
  19 class GrandCanonical {
  20  public:
  21   GrandCanonical();
  22   GrandCanonical(Int_t nmax, Double_t temperature, Double_t baryonPotential, Double_t strangePotential, Double_t electroPotential);
  23   ~GrandCanonical();
  24
  25   void     Temperature(Double_t value);
  26   Double_t Temperature() { return fTemperature; }
  27   void     BaryonPotential(Double_t value);
  28   Double_t BaryonPotential() { return fBaryonPotential; }
  29   void     StrangePotential(Double_t value);
  30   Double_t StrangePotential() { return fStrangePotential; }
  31   void     ElectroPotential(Double_t value);
  32   Double_t ElectroPotential() { return fElectroPotential; }
  33   void     NMax(Int_t value);
  34   Int_t    NMax() { return fNMax; }
  35
  36   // compute of system baryon number, system strangeness, system charge and
  37   // system energy
  38   // calculate system energy density
  39   Double_t EnergyDensity(DatabasePDG* database);
  40   // calculate system baryon density
  41   Double_t BaryonDensity(DatabasePDG* database);
  42   // calculate system strangeness density
  43   Double_t StrangeDensity(DatabasePDG* database);
  44   // calculate system electro density
  45   Double_t ElectroDensity(DatabasePDG* database);
  46   // compute of particle number density
  47   Double_t ParticleNumberDensity(ParticlePDG* particle);
  48   // compute the particle energy density
  49   Double_t ParticleEnergyDensity(ParticlePDG* particle);
  50
  51  private:
  52
  53   Double_t    fTemperature;
  54   Double_t    fBaryonPotential;
  55   Double_t    fStrangePotential;
  56   Double_t    fElectroPotential;
  57
  58   //  Number of terms for summation, if fNMax = 1 then
  59   //  Maxwell-Boltzmann distribution will be recovered
  60   Int_t       fNMax;
  61   Bool_t fInitialized;
  62 };
  63
  64 #endif