CosmicEoS.h

Go to the documentation of this file.

#ifndef COSMICEOS_H
#define COSMICEOS_H
#include <map>
 
#include "HRGBase/ThermalModelBase.h"
#include "HRGBase/Broyden.h"
#include "CosmicEos/EffectiveMassModel.h"
 
namespace thermalfist {

  double IsospinCharge(const ThermalParticle& part);

  struct LeptonFlavor {
    enum Name {
      Electron = 0, 
      Muon = 1,     
      Tau = 2,      
    };

    static const int NumberOfFlavors = 3;

    static double m_e;
    
    static double m_mu;

    static double m_tau;
 
  };

  class CosmicEoS
  {
  public:

    CosmicEoS(ThermalModelBase* THMbase, bool pionsinteract = false);

    virtual ~CosmicEoS(void) { }

    ThermalModelBase* HRGModel() const { return m_modelHRG; }

    const std::vector<double>& ChemicalPotentials() const { return m_ChemCurrent; }

    virtual void SetTemperature(double T);
    
    double Temperature() const { return m_T; }

    virtual void SetBaryonChemicalPotential(double muB) { m_ChemCurrent[0] = muB; m_modelHRG->SetBaryonChemicalPotential(muB); }
    
    double BaryonChemicalPotential() const { return m_ChemCurrent[0]; }

    virtual void SetElectricChemicalPotential(double muQ);// { m_ChemCurrent[1] = muQ; m_modelHRG->SetElectricChemicalPotential(muQ); }
    
    double ElectricChemicalPotential() const { return m_ChemCurrent[1]; }

    virtual void SetLeptonChemicalPotential(LeptonFlavor::Name flavor, double muL) { m_ChemCurrent[2 + static_cast<int>(flavor)] = muL; }
    
    double LeptonChemicalPotential(LeptonFlavor::Name flavor) const { return m_ChemCurrent[2 + static_cast<int>(flavor)]; }

    virtual void SetBaryonAsymmetry(double b) { m_Asymmetries[0] = b; }
    
    double BaryonAsymmetry() const { return m_Asymmetries[0]; }

    virtual void SetChargeAsymmetry(double q) { m_Asymmetries[1] = q; }
    
    double ChargeAsymmetry() const { return m_Asymmetries[1]; }

    virtual void SetLeptonAsymmetry(LeptonFlavor::Name flavor, double l) { m_Asymmetries[2 + static_cast<int>(flavor)] = l; }
    
    double LeptonAsymmetry(LeptonFlavor::Name flavor) const { return m_Asymmetries[2 + static_cast<int>(flavor)]; }

    void SetAsymmetries(const std::vector<double>& asymmetries) { m_Asymmetries = asymmetries; }

    void CalculatePrimordialDensities();

    double EntropyDensity();

    double EntropyDensityHRG();

    double Pressure();

    double PressureHRG();

    double EnergyDensity();

    double EnergyDensityHRG();

    double NetDensityChargedLepton(int iL);

    double PressureChargedLepton(int iL);

    double EnergyDensityChargedLepton(int iL);
    
    double BaryonDensity(bool absolute = false);

    double ElectricChargeDensity(bool absolute = false);

    double ElectricChargeDensityHRG(bool absolute = false);

    double IsospinChargeDensity(bool absolute = false);

    double LeptonFlavorDensity(LeptonFlavor::Name flavor, bool absolute = false);

    std::vector<double> SolveChemicalPotentials(double T, const std::vector<double>& muInit = std::vector<double>());

    static double fpi;

    double GetPionMass() const;

    void SetPionsInteracting(bool pionsinteract = true, double fpiChPT = fpi);
    
    bool InteractingPions() const { return m_InteractingPions; }

    bool InPionCondensedPhase() const;
 
    // std::vector<EffectiveMassModel>& EMMPions() { return m_Pions; }
    // double EMMPionCharge(int ipi) const { return m_PionCharges[ipi]; }

    const ThermalParticle& PhotonParticle() const { return m_Photon; }

    int NumberOfElectroWeakSpecies() const { return 1 + 2 * m_ChargedLeptons.size() + 2 * m_Neutrinos.size(); }

    std::string GetSpeciesName(int id) const;

    double GetDensity(int id) const;
 
  protected:
    ThermalModelBase* m_modelHRG;

    ThermalParticle m_Photon;

    std::vector<ThermalParticle> m_ChargedLeptons;
    
    std::vector<ThermalParticle> m_Neutrinos;
 
    // std::vector<EffectiveMassModel> m_Pions;
    // std::vector<double> m_PionCharges;

    bool m_InteractingPions;

    bool m_IsCalculated;

    double m_T;
    
    std::vector<double> m_ChemCurrent;
    
    std::vector<double> m_Asymmetries;

    void ClearEMMs();

    class BroydenEquationsCosmology : public BroydenEquations
    {
    public:
      BroydenEquationsCosmology(CosmicEoS* model) : BroydenEquations(), m_THM(model) { m_N = 2 + LeptonFlavor::NumberOfFlavors; }
      
      std::vector<double> Equations(const std::vector<double>& x);
    private:
      CosmicEoS* m_THM;
    };
  };
 
} // namespace thermalfist
 
#endif