GrandCanonical Class Reference

#include <GrandCanonical.h>

Public Member Functions
double	BaryonDensity (DatabasePDG *database)
void	BaryonPotential (double value)
double	BaryonPotential ()
double	CharmDensity (DatabasePDG *database)
void	CharmPotential (double value)
double	CharmPotential ()
double	ElectroDensity (DatabasePDG *database)
void	ElectroPotential (double value)
double	ElectroPotential ()
double	EnergyDensity (DatabasePDG *database)
	GrandCanonical ()
	GrandCanonical (int nmax, double temperature, double baryonPotential, double strangePotential, double electroPotential, double charmPotential)
void	NMax (int value)
int	NMax ()
double	ParticleEnergyDensity (ParticlePDG *particle)
double	ParticleNumberDensity (ParticlePDG *particle)
double	StrangeDensity (DatabasePDG *database)
void	StrangePotential (double value)
double	StrangePotential ()
void	Temperature (double value)
double	Temperature ()
	~GrandCanonical ()

Private Attributes
double	fBaryonPotential
double	fCharmPotential
double	fElectroPotential
bool	fInitialized
int	fNMax
double	fStrangePotential
double	fTemperature

Detailed Description

Definition at line 14 of file GrandCanonical.h.

Constructor & Destructor Documentation

GrandCanonical::GrandCanonical

(

)

Definition at line 15 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, and fTemperature.

                               {
  fInitialized = kFALSE;
  fNMax = -1111;
  fTemperature = -1111;
  fBaryonPotential = -1111;
  fStrangePotential = -1111;
  fElectroPotential = -1111;
  fCharmPotential = -1111;
}

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

Definition at line 25 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, and fTemperature.

                                                      {
  fNMax = nmax;
  fTemperature = temperature;
  fBaryonPotential = baryonPotential;
  fStrangePotential = strangePotential;
  fElectroPotential = electroPotential;
  fCharmPotential = charmPotential;
  fInitialized = kTRUE;
}

GrandCanonical::~GrandCanonical

(

)

Definition at line 40 of file GrandCanonical.cc.

{}

Member Function Documentation

double GrandCanonical::BaryonDensity

(

DatabasePDG *

database

)

Definition at line 146 of file GrandCanonical.cc.

References fInitialized, ParticlePDG::GetBaryonNumber(), DatabasePDG::GetNParticles(), DatabasePDG::GetPDGParticleByIndex(), and ParticleNumberDensity().

                                                          {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::BaryonDensity", "GrandCanonical object not fully initialized!!");

  double meanBaryonDensity = 0.;

  for (int currParticle = 0; currParticle < database->GetNParticles(); currParticle++) {
    ParticlePDG* particle = database->GetPDGParticleByIndex(currParticle);
    meanBaryonDensity += ParticleNumberDensity(particle) * particle->GetBaryonNumber();
  }
  return meanBaryonDensity;
}

void GrandCanonical::BaryonPotential

(

double

value

)

Definition at line 49 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                                 {
  fBaryonPotential = value;
  if (fNMax != -1111 && fTemperature != -1111 && fStrangePotential != -1111 && fElectroPotential != -1111 &&
      fCharmPotential != -1111)
    fInitialized = kTRUE;
}

double GrandCanonical::BaryonPotential ( )

inline

Definition at line 40 of file GrandCanonical.h.

References fBaryonPotential.

{ return fBaryonPotential; }

double GrandCanonical::CharmDensity

(

DatabasePDG *

database

)

Definition at line 191 of file GrandCanonical.cc.

References fInitialized, ParticlePDG::GetCharmness(), DatabasePDG::GetNParticles(), DatabasePDG::GetPDGParticleByIndex(), and ParticleNumberDensity().

                                                         {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::CharmDensity", "GrandCanonical object not fully initialized!!");

  double meanCharmDensity = 0.;

  for (int currParticle = 0; currParticle < database->GetNParticles(); currParticle++) {
    ParticlePDG* particle = database->GetPDGParticleByIndex(currParticle);
    meanCharmDensity += ParticleNumberDensity(particle) * particle->GetCharmness();
  }

  return meanCharmDensity;
}

void GrandCanonical::CharmPotential

(

double

value

)

Definition at line 70 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                                {
  fCharmPotential = value;
  if (fNMax != -1111 && fTemperature != -1111 && fBaryonPotential != -1111 && fStrangePotential != -1111 &&
      fElectroPotential != -1111)
    fInitialized = kTRUE;
}

double GrandCanonical::CharmPotential ( )

inline

Definition at line 46 of file GrandCanonical.h.

References fCharmPotential.

{ return fCharmPotential; }

double GrandCanonical::ElectroDensity

(

DatabasePDG *

database

)

Definition at line 175 of file GrandCanonical.cc.

References fInitialized, ParticlePDG::GetElectricCharge(), DatabasePDG::GetNParticles(), DatabasePDG::GetPDGParticleByIndex(), and ParticleNumberDensity().

                                                           {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::ElectroDensity", "GrandCanonical object not fully initialized!!");

  double meanElectroDensity = 0.;

  //hadrons
  for (int currParticle = 0; currParticle < database->GetNParticles(); currParticle++) {
    ParticlePDG* particle = database->GetPDGParticleByIndex(currParticle);
    meanElectroDensity += ParticleNumberDensity(particle) * particle->GetElectricCharge();
  }

  return meanElectroDensity;
}

void GrandCanonical::ElectroPotential

(

double

value

)

Definition at line 63 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                                  {
  fElectroPotential = value;
  if (fNMax != -1111 && fTemperature != -1111 && fBaryonPotential != -1111 && fStrangePotential != -1111 &&
      fCharmPotential != -1111)
    fInitialized = kTRUE;
}

double GrandCanonical::ElectroPotential ( )

inline

Definition at line 44 of file GrandCanonical.h.

References fElectroPotential.

{ return fElectroPotential; }

double GrandCanonical::EnergyDensity

(

DatabasePDG *

database

)

Definition at line 131 of file GrandCanonical.cc.

References fInitialized, DatabasePDG::GetNParticles(), DatabasePDG::GetPDGParticleByIndex(), and ParticleEnergyDensity().

                                                          {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::EnergyDensity", "GrandCanonical object not fully initialized!!");

  double meanEnergyDensity = 0.;

  for (int currParticle = 0; currParticle < database->GetNParticles(); currParticle++) {
    ParticlePDG* particle = database->GetPDGParticleByIndex(currParticle);
    meanEnergyDensity += ParticleEnergyDensity(particle);
  }

  return meanEnergyDensity;
}

void GrandCanonical::NMax

(

int

value

)

Definition at line 77 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                   {
  fNMax = value;
  if (fTemperature != -1111 && fBaryonPotential != -1111 && fStrangePotential != -1111 && fElectroPotential != -1111 &&
      fCharmPotential != -1111)
    fInitialized = kTRUE;
}

int GrandCanonical::NMax ( )

inline

Definition at line 49 of file GrandCanonical.h.

References fNMax.

{ return fNMax; }

double GrandCanonical::ParticleEnergyDensity

(

ParticlePDG *

particle

)

Definition at line 84 of file GrandCanonical.cc.

References ztail::d, fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, ParticlePDG::GetBaryonNumber(), ParticlePDG::GetCharmness(), ParticlePDG::GetElectricCharge(), ParticlePDG::GetMass(), ParticlePDG::GetSpin(), ParticlePDG::GetStrangeness(), HankelK1(), HankelKn(), hbarc, ResonanceBuilder::mass, dqmiodumpmetadata::n, and Pi.

Referenced by EnergyDensity().

                                                                  {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::ParticleEnergyDensity", "GrandCanonical object not fully initialized!!");

  // Compute the particle energy density
  double degFactor = 2. * particle->GetSpin() + 1.;        // degeneracy factor
  double mass = particle->GetMass();                       // PDG table mass
  double d = int(2. * particle->GetSpin()) & 1 ? 1. : -1;  // Bose-Einstein/Fermi-Dirac factor
  double preFactor = (degFactor * mass * mass * fTemperature * fTemperature / hbarc / hbarc / hbarc) /
                     (2. * TMath::Pi() * TMath::Pi());

  double postFactor = 0.;
  //compute chemical potential
  double potential = fBaryonPotential * particle->GetBaryonNumber() + fStrangePotential * particle->GetStrangeness() +
                     fElectroPotential * particle->GetElectricCharge() + fCharmPotential * particle->GetCharmness();

  for (int n = 1; n <= fNMax; ++n) {
    postFactor +=
        TMath::Power(-d, n + 1) / (n * n) * TMath::Exp(n * potential / fTemperature) *
        (3. * HankelKn(2, n * mass / fTemperature) + n * mass / fTemperature * HankelK1(n * mass / fTemperature));
  }
  return preFactor * postFactor;
}

double GrandCanonical::ParticleNumberDensity

(

ParticlePDG *

particle

)

Definition at line 109 of file GrandCanonical.cc.

References ztail::d, fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, ParticlePDG::GetBaryonNumber(), ParticlePDG::GetCharmness(), ParticlePDG::GetElectricCharge(), ParticlePDG::GetMass(), ParticlePDG::GetSpin(), ParticlePDG::GetStrangeness(), HankelKn(), hbarc, ResonanceBuilder::mass, dqmiodumpmetadata::n, and Pi.

Referenced by BaryonDensity(), CharmDensity(), ElectroDensity(), gen::Hydjet2Hadronizer::initializeForInternalPartons(), and StrangeDensity().

                                                                  {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::ParticleNumberDensity", "GrandCanonical object not fully initialized!!");

  double degFactor = 2. * particle->GetSpin() + 1.;
  double mass = particle->GetMass();
  double d = int(2 * particle->GetSpin()) & 1 ? 1. : -1.;
  double preFactor =
      (degFactor * mass * mass * fTemperature / hbarc / hbarc / hbarc) / (2. * TMath::Pi() * TMath::Pi());

  double postFactor = 0.;
  double potential = fBaryonPotential * particle->GetBaryonNumber() + fStrangePotential * particle->GetStrangeness() +
                     fElectroPotential * particle->GetElectricCharge() + fCharmPotential * particle->GetCharmness();

  for (int n = 1; n <= fNMax; ++n) {
    postFactor +=
        TMath::Power(-d, n + 1) / n * TMath::Exp(n * potential / fTemperature) * HankelKn(2, n * mass / fTemperature);
  }
  return preFactor * postFactor;
}

double GrandCanonical::StrangeDensity

(

DatabasePDG *

database

)

Definition at line 160 of file GrandCanonical.cc.

References fInitialized, DatabasePDG::GetNParticles(), DatabasePDG::GetPDGParticleByIndex(), ParticlePDG::GetStrangeness(), and ParticleNumberDensity().

                                                           {
  // Check if all the thermodinamic parameters are set
  if (!fInitialized)
    Fatal("GrandCanonical::StrangeDensity", "GrandCanonical object not fully initialized!!");

  double meanStrangeDensity = 0.;

  for (int currParticle = 0; currParticle < database->GetNParticles(); currParticle++) {
    ParticlePDG* particle = database->GetPDGParticleByIndex(currParticle);
    meanStrangeDensity += ParticleNumberDensity(particle) * particle->GetStrangeness();
  }

  return meanStrangeDensity;
}

void GrandCanonical::StrangePotential

(

double

value

)

Definition at line 56 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                                  {
  fStrangePotential = value;
  if (fNMax != -1111 && fTemperature != -1111 && fBaryonPotential != -1111 && fElectroPotential != -1111 &&
      fCharmPotential != -1111)
    fInitialized = kTRUE;
}

double GrandCanonical::StrangePotential ( )

inline

Definition at line 42 of file GrandCanonical.h.

References fStrangePotential.

{ return fStrangePotential; }

void GrandCanonical::Temperature

(

double

value

)

Definition at line 42 of file GrandCanonical.cc.

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and relativeConstraints::value.

                                             {
  fTemperature = value;
  if (fNMax != -1111 && fBaryonPotential != -1111 && fStrangePotential != -1111 && fElectroPotential != -1111 &&
      fCharmPotential != -1111)
    fInitialized = kTRUE;
}

double GrandCanonical::Temperature ( )

inline

Definition at line 38 of file GrandCanonical.h.

References fTemperature.

{ return fTemperature; }

Member Data Documentation

double GrandCanonical::fBaryonPotential

private

Definition at line 17 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().

double GrandCanonical::fCharmPotential

private

Definition at line 20 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().

double GrandCanonical::fElectroPotential

private

Definition at line 19 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().

bool GrandCanonical::fInitialized

private

Definition at line 25 of file GrandCanonical.h.

Referenced by BaryonDensity(), BaryonPotential(), CharmDensity(), CharmPotential(), ElectroDensity(), ElectroPotential(), EnergyDensity(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangeDensity(), StrangePotential(), and Temperature().

int GrandCanonical::fNMax

private

Definition at line 24 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().

double GrandCanonical::fStrangePotential

private

Definition at line 18 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().

double GrandCanonical::fTemperature

private

Definition at line 16 of file GrandCanonical.h.

Referenced by BaryonPotential(), CharmPotential(), ElectroPotential(), GrandCanonical(), NMax(), ParticleEnergyDensity(), ParticleNumberDensity(), StrangePotential(), and Temperature().