DIRE/Shower/Alpha_QCD.H

0001 #ifndef DIRE__Shower__Alpha_QCD_H
0002 #define DIRE__Shower__Alpha_QCD_H
0003
0004 #include "DIRE/Shower/Gauge.H"
0005
0006 namespace MODEL {
0007   class Running_AlphaS;
0008 }
0009
0010 namespace DIRE {
0011
0012   class Alpha_QCD: public Gauge {
0013   protected:
0014
0015     MODEL::Running_AlphaS *const &p_cpl;
0016
0017     double m_fac, m_max;
0018     double m_Nc, m_CF, m_CA, m_TR;
0019
0020     int m_lc;
0021
0022     mutable int m_override;
0023
0024     double B0(const double &nf) const;
0025     double B1(const double &nf) const;
0026
0027     double G2(const double &nf) const;
0028     double G3(const double &nf) const;
0029
0030     double CplFac(const double &scale) const;
0031
0032   public:
0033
0034     Alpha_QCD(const Kernel_Key &key);
0035
0036     double TrueScale(const Splitting &s) const;
0037
0038     double Coupling(const Splitting &s) const;
0039     double RenCT(const Splitting &s) const;
0040     double CplMax(const Splitting &s) const;
0041
0042     double K(const Splitting &s) const;
0043     double KMax(const Splitting &s) const;
0044
0045     double Nf(const Splitting &s) const;
0046
0047     void SetLimits();
0048
0049     double Solve(const double &cpl) const;
0050
0051   };// end of class Alpha_QCD
0052
0053 }// end of namespace DIRE
0054
0055 #endif