ThermalModelCanonicalCharm.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 THERMALMODELCanonicalCharm_H
#define THERMALMODELCanonicalCharm_H

#include <map>
#include <iostream>

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

namespace thermalfist {

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

    virtual ~ThermalModelCanonicalCharm(void);

    void CalculateEnergyDensitiesGCE();

    // Override functions begin

    virtual void SetParameters(const ThermalModelParameters& params);

    virtual void SetCharmChemicalPotential(double muC);

    virtual void ChangeTPS(ThermalParticleSystem *TPS);

    //For charm only Boltzmann supported
    virtual void SetStatistics(bool stats);

    void CalculateDensitiesGCE();

    virtual void FixParameters();

    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

  private:
    std::vector<double> m_densitiesGCE;
    std::vector<double> m_energydensitiesGCE;
    std::vector<int>    m_CharmValues;
    std::map<int, int>  m_CharmMap;
    std::vector<double> m_Zsum;
    std::vector<double> m_partialZ;
  };

} // namespace thermalfist

#endif // THERMALMODELCanonicalCharm_H