FreezeoutModels.h

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2019-2019 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#ifndef FREEZEOUTMODELS_H
#define FREEZEOUTMODELS_H
 
//#include "HRGEventGenerator/RandomGenerators.h"
#include <cmath>
 
namespace thermalfist {

  class BoostInvariantFreezeoutParametrization {
  public:
 
    BoostInvariantFreezeoutParametrization() : m_ProbabilityMaximumComputed(false), m_ProbabilityMaximum(1.) {}
    virtual ~BoostInvariantFreezeoutParametrization() {}

    virtual double Rfunc(double zeta) const { return zeta; }

    virtual double dRdZeta(double zeta) const { return (Rfunc(zeta + dzeta) - Rfunc(zeta)) / dzeta; }

    virtual double taufunc(double zeta) const { return 1.; }

    virtual double dtaudZeta(double zeta) const { return (taufunc(zeta + dzeta) - taufunc(zeta)) / dzeta; }

    virtual double etaperp(double zeta) const { return 0.; }
 
    virtual double sinhetaperp(double zeta) const { return sinh(etaperp(zeta)); }
    virtual double coshetaperp(double zeta) const { return cosh(etaperp(zeta)); }
    virtual double tanhetaperp(double zeta) const { return sinhetaperp(zeta) / coshetaperp(zeta); }

    virtual double ZetaProbability(double zeta) const;// { return zeta; }
 
    virtual double ProbabilityMaximum();

    //virtual double GetRandomZeta(MTRand& rangen = RandomGenerators::randgenMT);

    virtual bool InverseZetaDistributionIsExplicit() const { return false; }

    virtual double InverseZetaDistribution(double xi) const { return 0.; }
 
 
  protected:
    virtual double ComputeProbabilitydMaximum();
 
 
  private:
    virtual double TernarySearchForIntegrandMaximum(double zetaMin = 0., double zetaMax = 1.) const;
 
    static const double dzeta;
 
    bool m_ProbabilityMaximumComputed;
    double m_ProbabilityMaximum;
  };

  class CylindricalBlastWaveParametrization : public BoostInvariantFreezeoutParametrization {
  public:

    CylindricalBlastWaveParametrization(double betaSurface = 0.5, double nPower = 1., double tau = 10., double Rmax = 6.);
 
    virtual ~CylindricalBlastWaveParametrization() {}
 
    // Override functions begin
 
    virtual double Rfunc(double zeta) const { return zeta * m_R; }
 
    virtual double dRdZeta(double zeta) const { return m_R; }
 
    virtual double taufunc(double zeta) const { return m_tau; }
 
    virtual double dtaudZeta(double zeta) const { return 0.; }
 
    virtual double etaperp(double zeta) const { return atanh(m_BetaS * pow(zeta, m_n)); }
 
 
    virtual double sinhetaperp(double zeta) const { return m_BetaS * pow(zeta, m_n) / sqrt(1. - m_BetaS * m_BetaS * pow(zeta, 2. * m_n)); }
    virtual double coshetaperp(double zeta) const { return 1. / sqrt(1. - m_BetaS * m_BetaS * pow(zeta, 2. * m_n)); }
    virtual double tanhetaperp(double zeta) const { return m_BetaS * pow(zeta, m_n); }
 
    virtual double ZetaProbability(double zeta) const;
 
  protected:
    virtual double ComputeProbabilitydMaximum() { return ZetaProbability(1.); }
 
    // Override functions end
 
  private:
    double m_BetaS, m_n;
    double m_tau, m_R;
  };

  class CracowFreezeoutParametrization : public BoostInvariantFreezeoutParametrization {
  public:

    CracowFreezeoutParametrization(double RoverTauH = 1., double tauH = 10.);
 
    virtual ~CracowFreezeoutParametrization() {}
 
    // Override functions begin
 
    virtual double Rfunc(double zeta) const { return zeta * Rmax(); }
 
    virtual double dRdZeta(double zeta) const { return Rmax(); }
 
    virtual double taufunc(double zeta) const { return m_tauH * coshetaperp(zeta); }
 
    virtual double dtaudZeta(double zeta) const { return Rmax() * sinhetaperp(zeta) / coshetaperp(zeta); }
 
    virtual double etaperp(double zeta) const { return asinh(zeta * m_RoverTauH); }
 
 
    virtual double sinhetaperp(double zeta) const { return zeta * m_RoverTauH; }
    virtual double coshetaperp(double zeta) const { return sqrt(1. + zeta * zeta * m_RoverTauH * m_RoverTauH); }
 
    virtual double ZetaProbability(double zeta) const;
 
    //virtual double GetRandomZeta(MTRand& rangen = RandomGenerators::randgenMT);
 
  protected:
    virtual double ComputeProbabilitydMaximum() { return ZetaProbability(1.); }
 
    virtual bool InverseZetaDistributionIsExplicit() const { return true; }
 
    virtual double InverseZetaDistribution(double xi) const { return sqrt(xi); }
 
    // Override functions end
 
  private:
    double Rmax() const { return  m_RoverTauH * m_tauH; }
 
    double m_RoverTauH, m_tauH;
  };
 
} // namespace thermalfist
 
#endif