#include <GrandCanonical.h>

Public Member Functions
double	BaryonDensity (DatabasePDG *database)

double	BaryonPotential ()

void	BaryonPotential (double value)

double	CharmDensity (DatabasePDG *database)

double	CharmPotential ()

void	CharmPotential (double value)

double	ElectroDensity (DatabasePDG *database)

double	ElectroPotential ()

void	ElectroPotential (double value)

double	EnergyDensity (DatabasePDG *database)

	GrandCanonical ()

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

int	NMax ()

void	NMax (int value)

double	ParticleEnergyDensity (ParticlePDG *particle)

double	ParticleNumberDensity (ParticlePDG *particle)

double	StrangeDensity (DatabasePDG *database)

double	StrangePotential ()

void	StrangePotential (double value)

double	Temperature ()

void	Temperature (double value)

	~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() [1/2]

GrandCanonical::GrandCanonical ( )

Definition at line 15 of file GrandCanonical.cc.

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

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

◆ GrandCanonical() [2/2]

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

Definition at line 25 of file GrandCanonical.cc.

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

References fBaryonPotential, fCharmPotential, fElectroPotential, fInitialized, fNMax, fStrangePotential, fTemperature, and jvcParameters_cfi::temperature.

◆ ~GrandCanonical()

GrandCanonical::~GrandCanonical ( )

Definition at line 40 of file GrandCanonical.cc.

40 {}

Member Function Documentation

◆ BaryonDensity()

double GrandCanonical::BaryonDensity ( DatabasePDG * database )

Definition at line 146 of file GrandCanonical.cc.

                                                           {
   // 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;
 }

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

◆ BaryonPotential() [1/2]

double GrandCanonical::BaryonPotential ( )

inline

Definition at line 40 of file GrandCanonical.h.

40 { return fBaryonPotential; }

References fBaryonPotential.

◆ BaryonPotential() [2/2]

void GrandCanonical::BaryonPotential ( double value )

Definition at line 49 of file GrandCanonical.cc.

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

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

◆ CharmDensity()

double GrandCanonical::CharmDensity ( DatabasePDG * database )

Definition at line 191 of file GrandCanonical.cc.

                                                          {
   // 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;
 }

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

◆ CharmPotential() [1/2]

double GrandCanonical::CharmPotential ( )

inline

Definition at line 46 of file GrandCanonical.h.

46 { return fCharmPotential; }

References fCharmPotential.

◆ CharmPotential() [2/2]

void GrandCanonical::CharmPotential ( double value )

Definition at line 70 of file GrandCanonical.cc.

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

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

◆ ElectroDensity()

double GrandCanonical::ElectroDensity ( DatabasePDG * database )

Definition at line 175 of file GrandCanonical.cc.

                                                            {
   // 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;
 }

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

◆ ElectroPotential() [1/2]

double GrandCanonical::ElectroPotential ( )

inline

Definition at line 44 of file GrandCanonical.h.

44 { return fElectroPotential; }

References fElectroPotential.

◆ ElectroPotential() [2/2]

void GrandCanonical::ElectroPotential ( double value )

Definition at line 63 of file GrandCanonical.cc.

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

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

◆ EnergyDensity()

double GrandCanonical::EnergyDensity ( DatabasePDG * database )

Definition at line 131 of file GrandCanonical.cc.

                                                           {
   // 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;
 }

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

◆ NMax() [1/2]

int GrandCanonical::NMax ( )

inline

Definition at line 49 of file GrandCanonical.h.

49 { return fNMax; }

References fNMax.

◆ NMax() [2/2]

void GrandCanonical::NMax ( int value )

Definition at line 77 of file GrandCanonical.cc.

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

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

◆ ParticleEnergyDensity()

double GrandCanonical::ParticleEnergyDensity ( ParticlePDG * particle )

Definition at line 84 of file GrandCanonical.cc.

                                                                   {
   // 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;
 }

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, createfilelist::int, EgHLTOffHistBins_cfi::mass, dqmiodumpmetadata::n, and Pi.

Referenced by EnergyDensity().

◆ ParticleNumberDensity()

double GrandCanonical::ParticleNumberDensity ( ParticlePDG * particle )

Definition at line 109 of file GrandCanonical.cc.

                                                                   {
   // 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;
 }

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

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

◆ StrangeDensity()

double GrandCanonical::StrangeDensity ( DatabasePDG * database )

Definition at line 160 of file GrandCanonical.cc.

                                                            {
   // 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;
 }

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

◆ StrangePotential() [1/2]

double GrandCanonical::StrangePotential ( )

inline

Definition at line 42 of file GrandCanonical.h.

42 { return fStrangePotential; }

References fStrangePotential.

◆ StrangePotential() [2/2]

void GrandCanonical::StrangePotential ( double value )

Definition at line 56 of file GrandCanonical.cc.

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

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

◆ Temperature() [1/2]

double GrandCanonical::Temperature ( )

inline

Definition at line 38 of file GrandCanonical.h.

38 { return fTemperature; }

References fTemperature.

◆ Temperature() [2/2]

void GrandCanonical::Temperature ( double value )

Definition at line 42 of file GrandCanonical.cc.

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

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

Member Data Documentation

◆ fBaryonPotential

double GrandCanonical::fBaryonPotential

private

Definition at line 17 of file GrandCanonical.h.

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

◆ fCharmPotential

double GrandCanonical::fCharmPotential

private

Definition at line 20 of file GrandCanonical.h.

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

◆ fElectroPotential

double GrandCanonical::fElectroPotential

private

Definition at line 19 of file GrandCanonical.h.

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

◆ fInitialized

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().

◆ fNMax

int GrandCanonical::fNMax

private

Definition at line 24 of file GrandCanonical.h.

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

◆ fStrangePotential

double GrandCanonical::fStrangePotential

private

Definition at line 18 of file GrandCanonical.h.

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

◆ fTemperature

double GrandCanonical::fTemperature

private

Definition at line 16 of file GrandCanonical.h.

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GrandCanonical() [1/2]

◆ GrandCanonical() [2/2]

◆ ~GrandCanonical()

Member Function Documentation

◆ BaryonDensity()

◆ BaryonPotential() [1/2]

◆ BaryonPotential() [2/2]

◆ CharmDensity()

◆ CharmPotential() [1/2]

◆ CharmPotential() [2/2]

◆ ElectroDensity()

◆ ElectroPotential() [1/2]

◆ ElectroPotential() [2/2]

◆ EnergyDensity()

◆ NMax() [1/2]

◆ NMax() [2/2]

◆ ParticleEnergyDensity()

◆ ParticleNumberDensity()

◆ StrangeDensity()

◆ StrangePotential() [1/2]

◆ StrangePotential() [2/2]

◆ Temperature() [1/2]

◆ Temperature() [2/2]

Member Data Documentation

◆ fBaryonPotential

◆ fCharmPotential

◆ fElectroPotential

◆ fInitialized

◆ fNMax

◆ fStrangePotential

◆ fTemperature