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 double CalculateEnergyDensityDerivativeT() { exit(1); return 0.; } // Not implemented
 
    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