ThermalModelCanonicalStrangeness.h

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2014-2019 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#ifndef THERMALMODELCANONICALSTRANGENESS_H
#define THERMALMODELCANONICALSTRANGENESS_H

#include <map>
#include <iostream>

#include "HRGBase/ThermalModelBase.h"
#include "HRGBase/SplineFunction.h"

namespace thermalfist {

  class ThermalModelCanonicalStrangeness : public ThermalModelBase
  {
  public:
    ThermalModelCanonicalStrangeness(ThermalParticleSystem *TPS, const ThermalModelParameters& params = ThermalModelParameters());

    virtual ~ThermalModelCanonicalStrangeness(void);

    virtual void CalculateEnergyDensitiesGCE();

    virtual void CalculatePressuresGCE();

    virtual void CalculateSums(double Vc);

    const std::vector<double>& DensitiesGCE() const { return m_densitiesGCE; }

    // Override functions begin

    virtual void SetParameters(const ThermalModelParameters& params);

    virtual void SetStrangenessChemicalPotential(double muS);

    virtual void ChangeTPS(ThermalParticleSystem *TPS);

    virtual void FixParameters();

    virtual void SetStatistics(bool stats);

    virtual void CalculateDensitiesGCE();

    virtual void CalculatePrimordialDensities();

    void CalculateFluctuations();

    virtual double CalculatePressure();

    virtual double CalculateEnergyDensity();

    virtual double CalculateEntropyDensity();

    // Dummy
    virtual double CalculateBaryonMatterEntropyDensity() { return 0.; }

    virtual double CalculateMesonMatterEntropyDensity() { return 0.; }

    virtual double ParticleScaledVariance(int /*part*/) { return 1.; }

    virtual double ParticleScalarDensity(int /*part*/) { return 0.; }

    virtual bool IsConservedChargeCanonical(ConservedCharge::Name charge) const;

    // Override functions end

  protected:
    std::vector<double> m_densitiesGCE;
    std::vector<double> m_energydensitiesGCE;
    std::vector<double> m_pressuresGCE;
    std::vector<int>    m_StrVals;
    std::map<int, int>  m_StrMap;
    std::vector<double> m_Zsum;
    std::vector<double> m_partialS;
  };

} // namespace thermalfist

#endif // THERMALMODELCANONICALSTRANGENESS_H