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 35 of file GrandCanonical.cc.

{}

Member Function Documentation

double GrandCanonical::BaryonDensity

(

DatabasePDG *

database

)

Definition at line 140 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 43 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 37 of file GrandCanonical.h.

References fBaryonPotential.

{ return fBaryonPotential; }

double GrandCanonical::CharmDensity

(

DatabasePDG *

database

)

Definition at line 185 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 61 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 43 of file GrandCanonical.h.

References fCharmPotential.

{ return fCharmPotential; }

double GrandCanonical::ElectroDensity

(

DatabasePDG *

database

)

Definition at line 169 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 55 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 41 of file GrandCanonical.h.

References fElectroPotential.

{ return fElectroPotential; }

double GrandCanonical::EnergyDensity

(

DatabasePDG *

database

)

Definition at line 125 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 67 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 46 of file GrandCanonical.h.

References fNMax.

{ return fNMax; }

double GrandCanonical::ParticleEnergyDensity

(

ParticlePDG *

particle

)

Definition at line 74 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, gen::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 100 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, gen::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 154 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 49 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 39 of file GrandCanonical.h.

References fStrangePotential.

{ return fStrangePotential; }

void GrandCanonical::Temperature

(

double

value

)

Definition at line 37 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 35 of file GrandCanonical.h.

References fTemperature.

{ return fTemperature; }

Member Data Documentation

double GrandCanonical::fBaryonPotential

private

Definition at line 19 of file GrandCanonical.h.

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

double GrandCanonical::fCharmPotential

private

Definition at line 22 of file GrandCanonical.h.

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

double GrandCanonical::fElectroPotential

private

Definition at line 21 of file GrandCanonical.h.

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

bool GrandCanonical::fInitialized

private

Definition at line 27 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 26 of file GrandCanonical.h.

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

double GrandCanonical::fStrangePotential

private

Definition at line 20 of file GrandCanonical.h.

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

double GrandCanonical::fTemperature

private

Definition at line 18 of file GrandCanonical.h.

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