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