Class implementing HRG in partial chemical equilibrium with baryon annihilation. More...

#include <ThermalModelPCEAnnihilation.h>

Inheritance diagram for thermalfist::ThermalModelPCEAnnihilation:

Public Member Functions
	ThermalModelPCEAnnihilation (ThermalModelBase *THMbase, bool FreezeLonglived=false, double LonglivedResoWidthCut=0.015)
	Construct a new ThermalModelPCEAnnihilation object.

virtual	~ThermalModelPCEAnnihilation (void)
	Destroy the ThermalModelPCEAnnihilation object.

std::vector< int >	RecalculateStabilityFlags (const std::vector< long long > &annihilationpdgs={2212, 2112})

virtual void	SetStabilityFlags (const std::vector< int > &StabilityFlags)
	Manually set the PCE stability flags for all species.

virtual void	CalculatePCE (double param, PCEMode mode=AtFixedTemperature)
	Solves the equations of partial chemical equilibrium at a fixed temperature or a fixed volume.

void	SetPionAnnihilationNumber (double npi)
	Set the average number of pion produced in baryon-antibaryon annihilations.

Public Member Functions inherited from thermalfist::ThermalModelPCE
	ThermalModelPCE (ThermalModelBase *THMbase, bool FreezeLonglived=false, double LonglivedResoWidthCut=0.015)
	Construct a new ThermalModelPCE object.

virtual	~ThermalModelPCE (void)
	Destroy the ThermalModelPCE object.

void	UseSahaForNuclei (bool flag)
	Whether the nuclear abundances are evaluated through the Saha equation.

bool	UseSahaForNuclei () const

void	FreezeLonglivedResonances (bool flag)
	Whether long-lived resonances yields should be frozen.

bool	FreezeLonglivedResonances () const

void	SetLonglivedResonanceWidthCut (double width_cut)
	The threshold resonance width value to consider the resonance long-lived and its abundance frozen in the hadronic phase.

double	LonglivedResonanceWidthCut () const

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

void	SetChemicalFreezeout (const ThermalModelParameters &params, const std::vector< double > &ChemInit=std::vector< double >(0))
	Sets the chemical freeze-out conditions to be used as an initial condition for PCE calculations.

void	SetEntropyDensityChem (double sinit)
	Sets the entropy density at the chemical freeze-out.

double	EntropyDensityChem () const

ThermalModelBase *	ThermalModel () const
	Pointer to the HRG model used in calculations.

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

std::vector< double > &	ChemicalPotentials ()

double	Volume () const

void	ApplyFixForBoseCondensation ()
	Modifies the decay threshold masses of bosonic resonances such that the Bose-Condesation does not occur due to large fugacities.

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

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

Protected Member Functions
int	StableHadronIndexByGlobalId (int globalid)
	Returns the PCE-based index of the stable hadron based on its global (particle list) index.

std::vector< double >	StableChemsFromBroydenInput (const std::vector< double > &x)
	Chemical potentials of all PCE-based hadrons from the solution to PCE equations.

Additional Inherited Members
Public Types inherited from thermalfist::ThermalModelPCE
enum	PCEMode { AtFixedTemperature = 0 , AtFixedVolume = 1 }
	Whether partial chemical equilibrium should be calculated at a fixed value of the temperature or a fixed value of the volume. More...

Static Public Member Functions inherited from thermalfist::ThermalModelPCE
static void	PrepareNucleiForPCE (ThermalParticleSystem *TPS)
	Fills the "decay" products of light nuclei in accordance with their baryon content.

static std::vector< int >	ComputePCEStabilityFlags (const ThermalParticleSystem *TPS, bool SahaEquationForNuclei=true, bool FreezeLongLived=false, double WidthCut=0.015)
	Computes the PCE stability flags based on the provided particle list and a number of parameters.

Protected Attributes inherited from thermalfist::ThermalModelPCE
ThermalModelBase *	m_model

bool	m_UseSahaForNuclei
	Whether nuclear abundances are calculated via the Saha equation.

bool	m_FreezeLonglivedResonances
	Whether long-lived resonances are frozen at Tch.

double	m_ResoWidthCut
	Resonance width cut for freezeing the resonance abundances.

bool	m_ChemicalFreezeoutSet
	Whether the chemical freeze-out "initial" condition has been set.

bool	m_IsCalculated
	Whether PCE has been calculated.

bool	m_StabilityFlagsSet
	PCE configuration, list of stable species etc.

std::vector< int >	m_StabilityFlags

int	m_StableComponentsNumber

std::vector< std::vector< double > >	m_EffectiveCharges

std::vector< int >	m_StableMapTo

ThermalModelParameters	m_ParametersInit
	Parameters at the chemical freeze-out.

std::vector< double >	m_ChemInit

std::vector< double >	m_DensitiesInit

std::vector< double >	m_StableDensitiesInit

double	m_EntropyDensityInit

double	m_ParticleDensityInit

ThermalModelParameters	m_ParametersCurrent
	The current PCE thermal paratmeres and chemical potentials.

std::vector< double >	m_ChemCurrent

Detailed Description

Class implementing HRG in partial chemical equilibrium with baryon annihilation.

Implements equilibrium of N\bar{N} <-> <N_pi} pi reactions N is protons and neutrons by default, but can also include other stable baryons like Lambda's

The implementation was used (and described) in a paper V. Vovchenko, V. Koch, Phys. Lett. B 835, 137577 (2022), arXiv:2210.15641

A typical usage pattern is similar to ThermalModelPCE() except one can specify the annihilating (anti)baryons by a value 2 through SetStabilityFlags()

Definition at line 25 of file ThermalModelPCEAnnihilation.h.

Constructor & Destructor Documentation

◆ ThermalModelPCEAnnihilation()

thermalfist::ThermalModelPCEAnnihilation::ThermalModelPCEAnnihilation	(	ThermalModelBase *	THMbase,
		bool	FreezeLonglived = false,
		double	LonglivedResoWidthCut = 0.015 )

Construct a new ThermalModelPCEAnnihilation object.

Parameters

THMbase	A pointer to the ThermalModelBase object containing the HRG model implementation
FreezeLonglived	Whether long-lived resonance abundances should be frozen at Tch along with the stable hadrons
LonglivedResoWidthCut	Threshold width for the long-lived resonances to be considered stable in the PCE

Definition at line 10 of file ThermalModelPCEAnnihilation.cpp.

◆ ~ThermalModelPCEAnnihilation()

virtual thermalfist::ThermalModelPCEAnnihilation::~ThermalModelPCEAnnihilation ( void )

inlinevirtual

Destroy the ThermalModelPCEAnnihilation object.

Definition at line 42 of file ThermalModelPCEAnnihilation.h.

Member Function Documentation

◆ CalculatePCE()

void thermalfist::ThermalModelPCEAnnihilation::CalculatePCE	(	double	param,
		PCEMode	mode = AtFixedTemperature )

virtual

Solves the equations of partial chemical equilibrium at a fixed temperature or a fixed volume.

Parameters

param	Temperature (in GeV) if `mode` correspond to PCE at fixed temperature, Volume (in fm^3) if `mode` correspond to PCE at fixed volume.
mode	Determines whether the PCE calculation is perfored at a fixed temperature (default) or at a fixed volume.

Reimplemented from thermalfist::ThermalModelPCE.

Definition at line 66 of file ThermalModelPCEAnnihilation.cpp.

◆ RecalculateStabilityFlags()

std::vector< int > thermalfist::ThermalModelPCEAnnihilation::RecalculateStabilityFlags ( const std::vector< long long > & annihilationpdgs = {2212, 2112} )

Definition at line 213 of file ThermalModelPCEAnnihilation.cpp.

◆ SetPionAnnihilationNumber()

void thermalfist::ThermalModelPCEAnnihilation::SetPionAnnihilationNumber ( double npi )

inline

Set the average number of pion produced in baryon-antibaryon annihilations.

Parameters

npi	The average number of pion produced in baryon-antibaryon annihilations.

Definition at line 72 of file ThermalModelPCEAnnihilation.h.

◆ SetStabilityFlags()

void thermalfist::ThermalModelPCEAnnihilation::SetStabilityFlags ( const std::vector< int > & StabilityFlags )

virtual

Manually set the PCE stability flags for all species.

Parameters

StabilityFlags Vector of stability flags. Each element corresponds to a particle specie with the same 0-based index in the particle list. Element i equal to zero means particle yield is not frozen in the hadronic phase, if it is equal to one it means its yield is stable, and if it is equal to two it means the corresponding baryon annihilates.

Reimplemented from thermalfist::ThermalModelPCE.

Definition at line 21 of file ThermalModelPCEAnnihilation.cpp.