ThermalModelPCEAnnihilation.h

Go to the documentation of this file.

#ifndef THERMALMODELPCEANNIHILATION_H
#define THERMALMODELPCEANNIHILATION_H
#include <map>
 
#include "HRGPCE/ThermalModelPCE.h"
#include "HRGBase/Broyden.h"
 
namespace thermalfist {

  class ThermalModelPCEAnnihilation : public ThermalModelPCE
  {
  public:

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

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

    virtual void SetStabilityFlags(const std::vector<int>& StabilityFlags);

    virtual void CalculatePCE(double param, PCEMode mode = AtFixedTemperature);

    void SetPionAnnihilationNumber(double npi) { m_PionAnnihilationNumber = npi; }
 
  protected:

    int StableHadronIndexByGlobalId(int globalid);

    std::vector<double> StableChemsFromBroydenInput(const std::vector<double>& x);
 
  private:
 
    std::vector<int> m_StableNormal;
    std::vector<int> m_StableAnnihilate;
    std::vector<int> m_Pions;
    std::vector<int> m_StableAnnihilateAnti;
 
    double m_PionAnnihilationNumber;
 
    class BroydenEquationsPCEAnnihilation : public BroydenEquations
    {
    public:
      BroydenEquationsPCEAnnihilation(ThermalModelPCEAnnihilation *model, PCEMode mode = PCEMode::AtFixedTemperature) : BroydenEquations(), m_THM(model), m_Mode(mode) { m_N = m_THM->m_StableNormal.size() + m_THM->m_StableAnnihilate.size() + m_THM->m_Pions.size() + 1; }
      std::vector<double> Equations(const std::vector<double> &x);
    private:
      ThermalModelPCEAnnihilation *m_THM;
      PCEMode m_Mode;
    };
 
  };
 
} // namespace thermalfist
 
#endif