Class implementing the crossterms excluded-volume model. More...

#include <ThermalModelEVCrosstermsLegacy.h>

Inheritance diagram for thermalfist::ThermalModelEVCrosstermsLegacy:

Public Member Functions
	ThermalModelEVCrosstermsLegacy (ThermalParticleSystem *TPS, const ThermalModelParameters &params=ThermalModelParameters())
	Construct a new ThermalModelEVCrossterms object.

virtual	~ThermalModelEVCrosstermsLegacy (void)
	Destroy the ThermalModelEVCrossterms object.

virtual void	FillVirial (const std::vector< double > &ri=std::vector< double >(0))
	Fills the excluded volume coefficients \( \tilde{b}_{ij} \) based on the provided radii parameters for all species.

virtual void	ReadInteractionParameters (const std::string &filename)
	Reads the QvdW interaction parameters from a file.

virtual void	WriteInteractionParameters (const std::string &filename)
	Write the QvdW interaction parameters to a file.

void	SetRadius (double rad)
	Set the same excluded volume radius parameter for all species.

double	VirialCoefficient (int i, int j) const
	Excluded volume coefficient \( \tilde{b}_{ij} = 0 \).

void	SetVirial (int i, int j, double b)
	Set the excluded volume coefficient \( \tilde{b}_{ij} \).

virtual void	ChangeTPS (ThermalParticleSystem *TPS)
	Change the particle list.

virtual void	CalculatePrimordialDensities ()
	Calculates the primordial densities of all species.

void	CalculateTwoParticleCorrelations ()
	Computes the fluctuations and correlations of the primordial particle numbers.

void	CalculateFluctuations ()
	Computes the fluctuation observables.

virtual std::vector< double >	CalculateChargeFluctuations (const std::vector< double > &chgs, int order=4, bool dimensionfull=false)
	Calculates fluctuations (diagonal susceptibilities) of an arbitrary "conserved" charge.

virtual double	CalculatePressure ()

virtual double	CalculateEnergyDensity ()

virtual double	CalculateEntropyDensity ()

virtual double	CalculateBaryonMatterEntropyDensity ()

virtual double	CalculateMesonMatterEntropyDensity ()

virtual double	CalculateEnergyDensityDerivativeT ()

virtual double	CalculateEntropyDensityDerivativeT ()

virtual double	ParticleScaledVariance (int)

virtual double	ParticleSkewness (int part)

virtual double	ParticleKurtosis (int part)

virtual double	ParticleScalarDensity (int)

const std::vector< std::vector< int > > &	EVComponentIndices () const

virtual double	DeltaMu (int i) const

const std::vector< std::vector< double > > &	VirialMatrix () const

Public Member Functions inherited from thermalfist::ThermalModelBase
	ThermalModelBase (ThermalParticleSystem *TPS, const ThermalModelParameters &params=ThermalModelParameters())
	Construct a new ThermalModelBase object.

virtual	~ThermalModelBase (void)

int	ComponentsNumber () const
	Number of different particle species in the list.

bool	UseWidth () const
	Whether finite resonance widths are considered.

void	SetUseWidth (bool useWidth)
	Sets whether finite resonance widths are used. Deprecated.

void	SetUseWidth (ThermalParticle::ResonanceWidthIntegration type)
	Sets the finite resonance widths scheme to use.

void	SetNormBratio (bool normBratio)
	Whether branching ratios are renormalized to 100%.

bool	NormBratio () const

void	SetOMP (bool openMP)
	OpenMP support. Currently not used.

virtual void	SetParameters (const ThermalModelParameters &params)
	The thermal parameters.

const ThermalModelParameters &	Parameters () const

void	UpdateParameters ()
	Calls SetParameters() with current m_Parameters.

virtual void	SetTemperature (double T)
	Set the temperature.

virtual void	SetBaryonChemicalPotential (double muB)
	Set the baryon chemical potential.

virtual void	SetElectricChemicalPotential (double muQ)
	Set the electric chemical potential.

virtual void	SetStrangenessChemicalPotential (double muS)
	Set the strangeness chemical potential.

virtual void	SetCharmChemicalPotential (double muC)
	Set the charm chemical potential.

virtual void	SetGammaq (double gammaq)
	Set the light quark fugacity factor.

virtual void	SetGammaS (double gammaS)
	Set the strange quark fugacity factor.

virtual void	SetGammaC (double gammaC)
	Set the charm quark fugacity factor.

virtual void	SetBaryonCharge (int B)
	Set the total baryon number (for canonical ensemble only)

virtual void	SetElectricCharge (int Q)
	Set the total electric charge (for canonical ensemble only)

virtual void	SetStrangeness (int S)
	Set the total strangeness (for canonical ensemble only)

virtual void	SetCharm (int C)
	Set the total charm (for canonical ensemble only)

virtual void	SetRadius (int, double)
	Set the radius parameter for particle species i.

virtual void	SetRepulsion (int i, int j, double b)
	Same as SetVirial() but with a more clear name on what is actually does.

virtual void	SetAttraction (int, int, double)
	Set the vdW mean field attraction coefficient \( a_{ij} \).

virtual void	DisableMesonMesonVirial ()

void	DisableMesonMesonRepulsion ()

virtual void	DisableMesonMesonAttraction ()

virtual void	DisableMesonBaryonVirial ()

void	DisableMesonBaryonRepulsion ()

virtual void	DisableMesonBaryonAttraction ()

virtual void	DisableBaryonBaryonVirial ()

void	DisableBaryonBaryonRepulsion ()

virtual void	DisableBaryonBaryonAttraction ()

virtual void	DisableBaryonAntiBaryonVirial ()

void	DisableBaryonAntiBaryonRepulsion ()

virtual void	DisableBaryonAntiBaryonAttraction ()

double	RepulsionCoefficient (int i, int j) const

virtual double	AttractionCoefficient (int, int) const
	QvdW mean field attraction coefficient \( a_{ij} \).

bool	QuantumStatistics () const

virtual void	SetStatistics (bool stats)

virtual void	SetCalculationType (IdealGasFunctions::QStatsCalculationType type)
	Sets the CalculationType() method to evaluate quantum statistics. Calls the corresponding method in TPS().

virtual void	SetClusterExpansionOrder (int order)
	Set the number of terms in the cluster expansion method. Calls the corresponding method in TPS().

void	SetResonanceWidthShape (ThermalParticle::ResonanceWidthShape shape)
	Set the ThermalParticle::ResonanceWidthShape for all particles. Calls the corresponding method in TPS().

void	SetResonanceWidthIntegrationType (ThermalParticle::ResonanceWidthIntegration type)
	Set the ThermalParticle::ResonanceWidthIntegration scheme for all particles. Calls the corresponding method in TPS().

virtual void	FillChemicalPotentials ()
	Sets the chemical potentials of all particles.

virtual void	SetChemicalPotentials (const std::vector< double > &chem=std::vector< double >(0))
	Sets the chemical potentials of all particles.

const std::vector< double > &	ChemicalPotentials () const
	A vector of chemical potentials of all particles.

double	ChemicalPotential (int i) const
	Chemical potential of particle species i.

virtual void	SetChemicalPotential (int i, double chem)
	Sets the chemical potential of particle species i.

virtual double	FullIdealChemicalPotential (int i) const
	Chemical potential entering the ideal gas expressions of particle species i.

bool	ConstrainMuB () const

void	ConstrainMuB (bool constrain)

bool	ConstrainMuQ () const

void	ConstrainMuQ (bool constrain)

bool	ConstrainMuS () const

void	ConstrainMuS (bool constrain)

bool	ConstrainMuC () const

void	ConstrainMuC (bool constrain)

void	UsePartialChemicalEquilibrium (bool usePCE)
	Sets whether partial chemical equilibrium with additional chemical potentials is used.

bool	UsePartialChemicalEquilibrium ()
	Whether partial chemical equilibrium with additional chemical potentials is used.

void	SetSoverB (double SB)
	The entropy per baryon ratio to be used to constrain the baryon chemical potential.

double	SoverB () const

void	SetQoverB (double QB)
	The electric-to-baryon charge ratio to be used to constrain the electric chemical potential.

double	QoverB () const

void	SetVolume (double Volume)
	Sets the system volume.

double	Volume () const
	System volume (fm \(^3\))

void	SetVolumeRadius (double R)
	Sets the system radius.

double	CanonicalVolume () const
	The canonical correlation volume V \(_c\) (fm \(^3\))

void	SetCanonicalVolume (double Volume)
	Set the canonical correlation volume V \(_c\).

void	SetCanonicalVolumeRadius (double Rc)
	Set the canonical correlation system radius.

void	ConstrainChemicalPotentials (bool resetInitialValues=true)
	Constrains the chemical potentials \( \mu_B,\,\mu_Q,\,\mu_S,\,\mu_C \) by the conservation laws imposed.

virtual void	FixParameters ()
	Method which actually implements ConstrainChemicalPotentials() (for backward compatibility).

virtual void	FixParametersNoReset ()
	Method which actually implements ConstrainChemicalPotentialsNoReset() (for backward compatibility).

virtual bool	SolveChemicalPotentials (double totB=0., double totQ=0., double totS=0., double totC=0., double muBinit=0., double muQinit=0., double muSinit=0., double muCinit=0., bool ConstrMuB=true, bool ConstrMuQ=true, bool ConstrMuS=true, bool ConstrMuC=true)
	The procedure which calculates the chemical potentials \( \mu_B,\,\mu_Q,\,\mu_S,\,\mu_C \) which reproduce the specified total baryon, electric, strangeness, and charm charges of the system.

virtual bool	FixChemicalPotentialsThroughDensities (double rhoB=0., double rhoQ=0., double rhoS=0., double rhoC=0., double muBinit=0., double muQinit=0., double muSinit=0., double muCinit=0., bool ConstrMuB=true, bool ConstrMuQ=true, bool ConstrMuS=true, bool ConstrMuC=true)
	The procedure which calculates the chemical potentials \( \mu_B,\,\mu_Q,\,\mu_S,\,\mu_C \) which reproduce the specified baryon, electric, strangeness, and charm densities.

virtual void	CalculateDensities ()
	Calculates the primordial and total (after decays) densities of all species.

virtual void	ValidateCalculation ()
	Checks whether issues have occured during the calculation of particle densities in the CalculateDensities() method.

std::string	ValidityCheckLog () const
	All messaged which occured during the validation procedure in the ValidateCalculation() method.

virtual void	CalculateDensitiesGCE ()
	Calculates the particle densities in a grand-canonical ensemble.

virtual void	CalculateFeeddown ()
	Calculates the total densities which include feeddown contributions.

virtual void	CalculateTemperatureDerivatives ()
	Computes the temperature derivatives of densities, shifted chemical potentials, and primordial hadron number susceptibilities.

virtual void	CalculateTwoParticleFluctuationsDecays ()
	Computes particle number correlations and fluctuations for all final-state particles which are marked stable.

virtual double	TwoParticleSusceptibilityPrimordial (int i, int j) const
	Returns the computed primordial particle number (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with ids i and j. CalculateFluctuations() must be called beforehand.

virtual double	TwoParticleSusceptibilityPrimordialByPdg (long long id1, long long id2)
	Returns the computed primordial particle number (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with pdg codes id1 and id2. CalculateFluctuations() must be called beforehand.

virtual double	TwoParticleSusceptibilityTemperatureDerivativePrimordial (int i, int j) const
	Returns the computed temperature derivative of the primordial particle number (cross-)susceptibility \( \frac{\partial}{\partial T} \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with ids i and j. CalculateFluctuations() must be called beforehand.

virtual double	TwoParticleSusceptibilityTemperatureDerivativePrimordialByPdg (long long id1, long long id2)
	Returns the computed temperature derivative of the primordial particle number (cross-)susceptibility \( \frac{\partial}{\partial T} \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with pdg codes id1 and id2. CalculateFluctuations() must be called beforehand.

virtual double	NetParticleSusceptibilityPrimordialByPdg (long long id1, long long id2)
	Returns the computed (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) between primordial net-particle numbers for pdg codes id1 and id2. CalculateFluctuations() must be called beforehand.

virtual double	TwoParticleSusceptibilityFinal (int i, int j) const
	Returns the computed final particle number (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with ids i and j. CalculateFluctuations() must be called beforehand. Both particle species must be those marked stable.

virtual double	TwoParticleSusceptibilityFinalByPdg (long long id1, long long id2)
	Returns the computed final particle number (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) for particles with pdg codes id1 and id2. CalculateFluctuations() must be called beforehand. Both particle species must be those marked stable.

virtual double	NetParticleSusceptibilityFinalByPdg (long long id1, long long id2)
	Returns the computed (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_j \rangle \) between final net-particle numbers for pdg codes id1 and id2. CalculateFluctuations() must be called beforehand.

virtual double	PrimordialParticleChargeSusceptibility (int i, ConservedCharge::Name chg) const
	Returns the computed primordial particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with id i and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	PrimordialParticleChargeSusceptibilityByPdg (long long id1, ConservedCharge::Name chg)
	Returns the computed primordial particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with pdg code id1 and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	PrimordialNetParticleChargeSusceptibilityByPdg (long long id1, ConservedCharge::Name chg)
	Returns the computed primordial net particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with pdg code id1 and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	FinalParticleChargeSusceptibility (int i, ConservedCharge::Name chg) const
	Returns the computed final (after decays) particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with id i and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	FinalParticleChargeSusceptibilityByPdg (long long id1, ConservedCharge::Name chg)
	Returns the computed final (after decays) particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with pdg code id1 and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	FinalNetParticleChargeSusceptibilityByPdg (long long id1, ConservedCharge::Name chg)
	Returns the computed final (after decays) net particle vs conserved charge (cross-)susceptibility \( \frac{1}{VT^3} \, \langle \Delta N_i \Delta N_chg \rangle \) for particle with pdg code id1 and conserved charge chg. CalculateFluctuations() must be called beforehand.

virtual double	SusceptibilityDimensionfull (ConservedCharge::Name i, ConservedCharge::Name j) const
	A 2nd order susceptibility of conserved charges.

virtual void	CalculateSusceptibilityMatrix ()
	Calculates the conserved charges susceptibility matrix.

virtual void	CalculateProxySusceptibilityMatrix ()
	Calculates the susceptibility matrix of conserved charges proxies.

virtual void	CalculateParticleChargeCorrelationMatrix ()
	Calculates the matrix of correlators between primordial (and also final) particle numbers and conserved charges.

Protected Member Functions
virtual void	SolvePressure (bool resetPartials=true)
	Solves the system of transcdental equations for the pressure using the Broyden's method.

void	SolvePressureIter ()

virtual void	CalculatePrimordialDensitiesNoReset ()

virtual void	CalculatePrimordialDensitiesIter ()

virtual void	SolveDiagonal ()
	Solves the transcendental equation of the corresponding diagonal EV model.

double	PartialPressureDiagonal (int i, double P)
	The "partial pressure" of hadron species i for the given total pressure in the diagonal model.

double	PressureDiagonalTotal (double P)
	The total pressure for the given input pressure in the diagonal model.

virtual double	DensityId (int i, const std::vector< double > &pstars=std::vector< double >())
	Calculate the ideal gas density of particle species i for the given values of partial pressures.

virtual double	Pressure (int i, const std::vector< double > &pstars=std::vector< double >())
	Calculate the ideal gas pressure of particle species i for the given values of partial pressures.

double	ScaledVarianceId (int ind)
	Calculate the ideal gas scaled variance of particle species i number fluctuations for the given values of partial pressures.

virtual double	MuShift (int i) const
	The shift in the chemical potential of particle species i due to the excluded volume interactions.

Protected Attributes
std::vector< double >	m_densitiesid

std::vector< double >	m_Ps

std::vector< std::vector< double > >	m_Virial

double	m_Pressure

double	m_TotalEntropyDensity

std::vector< int >	m_MapToEVComponent

std::vector< int >	m_MapFromEVComponent

std::vector< std::vector< int > >	m_EVComponentIndices

Additional Inherited Members
Public Types inherited from thermalfist::ThermalModelBase
enum	ThermalModelEnsemble { GCE = 0 , CE = 1 , SCE = 2 , CCE = 3 }
	The list of statistical ensembles. More...

enum	ThermalModelInteraction { Ideal = 0 , DiagonalEV = 1 , CrosstermsEV = 2 , QvdW = 3 , RealGas = 4 , MeanField = 5 }
	Type of interactions included in the HRG model. More...

Detailed Description

Class implementing the crossterms excluded-volume model.

The model formulation can be found in

M.I. Gorenstein, A.P. Kostyuk, Y.D. Krivenko, J.Phys. G 25, L75 (1999), http://arxiv.org/pdf/nucl-th/9906068.pdf

and in

V. Vovchenko, H. Stoecker, Phys. Rev. C 95, 044904 (2017), http://arxiv.org/pdf/1606.06218.pdf

The system of transcendental equations for "partial pressures" of hadrons is solved using the Broyden's method.

Deprecated: Since version 1.4 ThermalModelEVCrossterms() is implemented as partial case of ThermalModelVDW() The implementation in this file is not a legacy implementation

Definition at line 38 of file ThermalModelEVCrosstermsLegacy.h.

Constructor & Destructor Documentation

◆ ThermalModelEVCrosstermsLegacy()

thermalfist::ThermalModelEVCrosstermsLegacy::ThermalModelEVCrosstermsLegacy	(	ThermalParticleSystem *	TPS,
		const ThermalModelParameters &	params = ThermalModelParameters() )

Construct a new ThermalModelEVCrossterms object.

Parameters

TPS	A pointer to the ThermalParticleSystem object containing the particle list
params	ThermalModelParameters object with current thermal parameters

Definition at line 33 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ ~ThermalModelEVCrosstermsLegacy()

thermalfist::ThermalModelEVCrosstermsLegacy::~ThermalModelEVCrosstermsLegacy ( void )

virtual

Destroy the ThermalModelEVCrossterms object.

Definition at line 46 of file ThermalModelEVCrosstermsLegacy.cpp.

Member Function Documentation

◆ CalculateBaryonMatterEntropyDensity()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateBaryonMatterEntropyDensity ( )

inlinevirtual

Definition at line 86 of file ThermalModelEVCrosstermsLegacy.h.

◆ CalculateChargeFluctuations()

std::vector< double > thermalfist::ThermalModelEVCrosstermsLegacy::CalculateChargeFluctuations	(	const std::vector< double > &	chgs,
		int	order = 4,
		bool	dimensionfull = false )

virtual

Calculates fluctuations (diagonal susceptibilities) of an arbitrary "conserved" charge.

Each particle specie is assumed to carry a conserved charge with a value provided by an input vector.

Restricted to the grand canonical ensemble.

Parameters

chgs	A vector with conserved charge values for all species. 0-based indices of the vector must correspond to the 0-based indices of the particle list TPS()
order	Up to which order the susceptibilities are computed

Returns: std::vector<double> A vector with computed values of diagonal susceptibilities

Parameters

dimensionfull If true, return d^n P / dmu^n (dimensionfull). If false (default), return d^n (P/T^4) / d(mu/T)^n (dimensionless). The dimensionfull option avoids divergences at T=0.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 600 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculateEnergyDensity()

double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateEnergyDensity ( )

virtual

Definition at line 796 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculateEnergyDensityDerivativeT()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateEnergyDensityDerivativeT ( )

inlinevirtual

Definition at line 89 of file ThermalModelEVCrosstermsLegacy.h.

◆ CalculateEntropyDensity()

double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateEntropyDensity ( )

virtual

Definition at line 805 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculateEntropyDensityDerivativeT()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateEntropyDensityDerivativeT ( )

inlinevirtual

Definition at line 91 of file ThermalModelEVCrosstermsLegacy.h.

◆ CalculateFluctuations()

void thermalfist::ThermalModelEVCrosstermsLegacy::CalculateFluctuations ( )

virtual

Computes the fluctuation observables.

Includes the matrix of 2nd order susceptibilities of conserved charges, as well as particle number correlations and fluctuations, both for primordial yields and after decays.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 583 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculateMesonMatterEntropyDensity()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::CalculateMesonMatterEntropyDensity ( )

inlinevirtual

Definition at line 87 of file ThermalModelEVCrosstermsLegacy.h.

◆ CalculatePressure()

double thermalfist::ThermalModelEVCrosstermsLegacy::CalculatePressure ( )

virtual

Definition at line 810 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculatePrimordialDensities()

void thermalfist::ThermalModelEVCrosstermsLegacy::CalculatePrimordialDensities ( )

virtual

Calculates the primordial densities of all species.

Implements thermalfist::ThermalModelBase.

Definition at line 240 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculatePrimordialDensitiesIter()

void thermalfist::ThermalModelEVCrosstermsLegacy::CalculatePrimordialDensitiesIter ( )

protectedvirtual

Definition at line 426 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculatePrimordialDensitiesNoReset()

void thermalfist::ThermalModelEVCrosstermsLegacy::CalculatePrimordialDensitiesNoReset ( )

protectedvirtual

Definition at line 332 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ CalculateTwoParticleCorrelations()

void thermalfist::ThermalModelEVCrosstermsLegacy::CalculateTwoParticleCorrelations ( )

virtual

Computes the fluctuations and correlations of the primordial particle numbers.

More specifically, computes the susceptibility matrix \( \frac{1}{VT^3} \, \langle \Delta N_i^* \Delta N_j^* \rangle \), where \( N_i^* \) is the primordial yield of species i.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 489 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ ChangeTPS()

void thermalfist::ThermalModelEVCrosstermsLegacy::ChangeTPS ( ThermalParticleSystem * TPS )

virtual

Change the particle list.

Parameters

TPS	A pointer to new particle list.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 143 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ DeltaMu()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::DeltaMu ( int i ) const

inlinevirtual

Definition at line 109 of file ThermalModelEVCrosstermsLegacy.h.

◆ DensityId()

double thermalfist::ThermalModelEVCrosstermsLegacy::DensityId	(	int	i,
		const std::vector< double > &	pstars = std::vector<double>() )

protectedvirtual

Calculate the ideal gas density of particle species i for the given values of partial pressures.

Parameters

i	0-based particle specie index
Pressure	Input vector of partial pressures (GeV fm \(^{-3}\))

Returns: Ideal gas density (fm \(^{-3}\))

Definition at line 149 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ EVComponentIndices()

const std::vector< std::vector< int > > & thermalfist::ThermalModelEVCrosstermsLegacy::EVComponentIndices ( ) const

inline

Definition at line 108 of file ThermalModelEVCrosstermsLegacy.h.

◆ FillVirial()

void thermalfist::ThermalModelEVCrosstermsLegacy::FillVirial ( const std::vector< double > & ri = std::vector<double>(0) )

virtual

Fills the excluded volume coefficients \( \tilde{b}_{ij} \) based on the provided radii parameters for all species.

Fills the coefficients in accordance with Eqs. (5) and (7) here https://arxiv.org/pdf/1606.06218.pdf

Parameters

ri	A vector with radii parameters for all species. 0-based indices of the vector must corresponds to the 0-based indices of the particle list TPS()

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 50 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ MuShift()

double thermalfist::ThermalModelEVCrosstermsLegacy::MuShift ( int i ) const

protectedvirtual

The shift in the chemical potential of particle species i due to the excluded volume interactions.

Equal to \(- \sum_j \tilde{b_{ij}} p_j\)

Parameters

i	0-based particle specie index

Returns: The shift in the chemical potential

Definition at line 816 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ PartialPressureDiagonal()

double thermalfist::ThermalModelEVCrosstermsLegacy::PartialPressureDiagonal	(	int	i,
		double	P )

protected

The "partial pressure" of hadron species i for the given total pressure in the diagonal model.

Parameters

i	0-based particle species index.
P	Input pressure (GeV fm \(^{-3}\))

Returns: Computed partial pressure (GeV fm \(^{-3}\))

Definition at line 182 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ ParticleKurtosis()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::ParticleKurtosis ( int part )

inlinevirtual

Definition at line 100 of file ThermalModelEVCrosstermsLegacy.h.

◆ ParticleScalarDensity()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::ParticleScalarDensity ( int )

inlinevirtual

Definition at line 103 of file ThermalModelEVCrosstermsLegacy.h.

◆ ParticleScaledVariance()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::ParticleScaledVariance ( int )

inlinevirtual

Definition at line 94 of file ThermalModelEVCrosstermsLegacy.h.

◆ ParticleSkewness()

virtual double thermalfist::ThermalModelEVCrosstermsLegacy::ParticleSkewness ( int part )

inlinevirtual

Definition at line 97 of file ThermalModelEVCrosstermsLegacy.h.

◆ Pressure()

double thermalfist::ThermalModelEVCrosstermsLegacy::Pressure	(	int	i,
		const std::vector< double > &	pstars = std::vector<double>() )

protectedvirtual

Calculate the ideal gas pressure of particle species i for the given values of partial pressures.

Parameters

i	0-based particle specie index
Pressure	Input vector of partial pressures (GeV fm \(^{-3}\))

Returns: Ideal gas pressure (fm \(^{-3}\))

Definition at line 162 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ PressureDiagonalTotal()

double thermalfist::ThermalModelEVCrosstermsLegacy::PressureDiagonalTotal ( double P )

protected

The total pressure for the given input pressure in the diagonal model.

Parameters

Input pressure (GeV fm \(^{-3}\))

Returns: Computed pressure (GeV fm \(^{-3}\))

Definition at line 190 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ ReadInteractionParameters()

void thermalfist::ThermalModelEVCrosstermsLegacy::ReadInteractionParameters ( const std::string & )

virtual

Reads the QvdW interaction parameters from a file.

Actual implementation is in a derived class.

Parameters

filename File with interaction parameters.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 71 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ ScaledVarianceId()

double thermalfist::ThermalModelEVCrosstermsLegacy::ScaledVarianceId ( int ind )

protected

Calculate the ideal gas scaled variance of particle species i number fluctuations for the given values of partial pressures.

Parameters

i	0-based particle specie index
Pressure	Input vector of partial pressures (GeV fm \(^{-3}\))

Returns: Ideal gas scaled variance

Definition at line 175 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ SetRadius()

void thermalfist::ThermalModelEVCrosstermsLegacy::SetRadius ( double )

virtual

Set the same excluded volume radius parameter for all species.

Parameters

rad	Radius parameter (fm)

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 127 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ SetVirial()

void thermalfist::ThermalModelEVCrosstermsLegacy::SetVirial	(	int	,
		int	,
		double	)

virtual

Set the excluded volume coefficient \( \tilde{b}_{ij} \).

Excluded parameter for repulsive interaction between particle species i and j.

Parameters

i	0-based index of the first particle species
j	0-based index of the second particle species
b	Excluded volume parameter \( \tilde{b}_{ij} \) (fm \(^3\))

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 137 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ SolveDiagonal()

void thermalfist::ThermalModelEVCrosstermsLegacy::SolveDiagonal ( )

protectedvirtual

Solves the transcendental equation of the corresponding diagonal EV model.

The diagonal EV model here has \(v_i = \tilde{b}_{ii}\).

The partial pressures of the diagonal model will be recorded into m_Ps.

Definition at line 197 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ SolvePressure()

void thermalfist::ThermalModelEVCrosstermsLegacy::SolvePressure ( bool resetPartials = true )

protectedvirtual

Solves the system of transcdental equations for the pressure using the Broyden's method.

Solves \( p_i(T,\mu) = p_i^{\rm id} (T, \mu_i - \sum_j \tilde{b_{ij}} p_j). \)

Definition at line 216 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ SolvePressureIter()

void thermalfist::ThermalModelEVCrosstermsLegacy::SolvePressureIter ( )

protected

Definition at line 401 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ VirialCoefficient()

double thermalfist::ThermalModelEVCrosstermsLegacy::VirialCoefficient	(	int	,
		int	) const

virtual

Excluded volume coefficient \( \tilde{b}_{ij} = 0 \).

Excluded parameter for repulsive interaction between particle species i and j.

Parameters

i	0-based index of the first particle species
j	0-based index of the second particle species

Returns: double Coefficient \( \tilde{b}_{ij} = 0 \)

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 131 of file ThermalModelEVCrosstermsLegacy.cpp.

◆ VirialMatrix()

const std::vector< std::vector< double > > & thermalfist::ThermalModelEVCrosstermsLegacy::VirialMatrix ( ) const

inline

Definition at line 110 of file ThermalModelEVCrosstermsLegacy.h.

◆ WriteInteractionParameters()

void thermalfist::ThermalModelEVCrosstermsLegacy::WriteInteractionParameters ( const std::string & )

virtual

Write the QvdW interaction parameters to a file.

Actual implementation is in a derived class.

Parameters

filename Output file.

Reimplemented from thermalfist::ThermalModelBase.

Definition at line 96 of file ThermalModelEVCrosstermsLegacy.cpp.

Member Data Documentation

◆ m_densitiesid

std::vector<double> thermalfist::ThermalModelEVCrosstermsLegacy::m_densitiesid

protected

Vector of ideal gas densities with shifted chemical potentials

Definition at line 205 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_EVComponentIndices

std::vector< std::vector<int> > thermalfist::ThermalModelEVCrosstermsLegacy::m_EVComponentIndices

protected

Definition at line 214 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_MapFromEVComponent

std::vector<int> thermalfist::ThermalModelEVCrosstermsLegacy::m_MapFromEVComponent

protected

Definition at line 213 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_MapToEVComponent

std::vector<int> thermalfist::ThermalModelEVCrosstermsLegacy::m_MapToEVComponent

protected

Definition at line 212 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_Pressure

double thermalfist::ThermalModelEVCrosstermsLegacy::m_Pressure

protected

The (solved) total pressure

Definition at line 208 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_Ps

std::vector<double> thermalfist::ThermalModelEVCrosstermsLegacy::m_Ps

protected

Vector of (solved) partial pressures

Definition at line 206 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_TotalEntropyDensity

double thermalfist::ThermalModelEVCrosstermsLegacy::m_TotalEntropyDensity

protected

The (solved) entropy pressure

Definition at line 209 of file ThermalModelEVCrosstermsLegacy.h.

◆ m_Virial

std::vector< std::vector<double> > thermalfist::ThermalModelEVCrosstermsLegacy::m_Virial

protected

Matrix of virial (excluded-volume) coefficients \( \tilde{b}_{ij} \)

Definition at line 207 of file ThermalModelEVCrosstermsLegacy.h.

The documentation for this class was generated from the following files:

include/HRGEV/ThermalModelEVCrosstermsLegacy.h
src/library/HRGEV/ThermalModelEVCrosstermsLegacy.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ThermalModelEVCrosstermsLegacy()

◆ ~ThermalModelEVCrosstermsLegacy()

Member Function Documentation

◆ CalculateBaryonMatterEntropyDensity()

◆ CalculateChargeFluctuations()

◆ CalculateEnergyDensity()

◆ CalculateEnergyDensityDerivativeT()

◆ CalculateEntropyDensity()

◆ CalculateEntropyDensityDerivativeT()

◆ CalculateFluctuations()

◆ CalculateMesonMatterEntropyDensity()

◆ CalculatePressure()

◆ CalculatePrimordialDensities()

◆ CalculatePrimordialDensitiesIter()

◆ CalculatePrimordialDensitiesNoReset()

◆ CalculateTwoParticleCorrelations()

◆ ChangeTPS()

◆ DeltaMu()

◆ DensityId()

◆ EVComponentIndices()

◆ FillVirial()

◆ MuShift()

◆ PartialPressureDiagonal()

◆ ParticleKurtosis()

◆ ParticleScalarDensity()

◆ ParticleScaledVariance()

◆ ParticleSkewness()

◆ Pressure()

◆ PressureDiagonalTotal()

◆ ReadInteractionParameters()

◆ ScaledVarianceId()

◆ SetRadius()

◆ SetVirial()

◆ SolveDiagonal()

◆ SolvePressure()

◆ SolvePressureIter()

◆ VirialCoefficient()

◆ VirialMatrix()

◆ WriteInteractionParameters()

Member Data Documentation

◆ m_densitiesid

◆ m_EVComponentIndices

◆ m_MapFromEVComponent

◆ m_MapToEVComponent

◆ m_Pressure

◆ m_Ps

◆ m_TotalEntropyDensity

◆ m_Virial