GrandCanonical.h

Go to the documentation of this file.

/*                                                                            
                                                                            
Nikolai Amelin, Ludmila Malinina, Timur Pocheptsov (C) JINR/Dubna
amelin@sunhe.jinr.ru, malinina@sunhe.jinr.ru, pocheptsov@sunhe.jinr.ru 
November. 2, 2005                                

*/

#ifndef GRANDCANONICAL_INCLUDED
#define GRANDCANONICAL_INCLUDED
#include "ParticlePDG.h"
#include "DatabasePDG.h"

class GrandCanonical {
private:
  double fTemperature;
  double fBaryonPotential;
  double fStrangePotential;
  double fElectroPotential;
  double fCharmPotential;

  //  Number of terms for summation, if fNMax = 1 then
  //  Maxwell-Boltzmann distribution will be recovered
  int fNMax;
  bool fInitialized;

public:
  GrandCanonical();
  GrandCanonical(int nmax,
                 double temperature,
                 double baryonPotential,
                 double strangePotential,
                 double electroPotential,
                 double charmPotential);
  ~GrandCanonical();

  void Temperature(double value);
  double Temperature() { return fTemperature; }
  void BaryonPotential(double value);
  double BaryonPotential() { return fBaryonPotential; }
  void StrangePotential(double value);
  double StrangePotential() { return fStrangePotential; }
  void ElectroPotential(double value);
  double ElectroPotential() { return fElectroPotential; }
  void CharmPotential(double value);
  double CharmPotential() { return fCharmPotential; }

  void NMax(int value);
  int NMax() { return fNMax; }

  // compute of system baryon number, system strangeness, system charge and
  // system energy
  // calculate system energy density
  double EnergyDensity(DatabasePDG* database);
  // calculate system baryon density
  double BaryonDensity(DatabasePDG* database);
  // calculate system strangeness density
  double StrangeDensity(DatabasePDG* database);
  // calculate system electro density
  double ElectroDensity(DatabasePDG* database);
  // compute of particle number density
  double CharmDensity(DatabasePDG* database);

  // compute of particle number density
  double ParticleNumberDensity(ParticlePDG* particle);
  // compute the particle energy density
  double ParticleEnergyDensity(ParticlePDG* particle);
};

#endif