lowe.h

Go to the documentation of this file.

 
#ifndef LOWE_H
#define LOWE_H
 
#include <iostream>
#include <fstream>
#include <sstream>
#include <cstring>
#include "utils.h"
#include "linalg.h"
#include "rge.h"
 
namespace softsusy {
const double MUP = 2.2e-3; 
const double MDOWN = 4.7e-3; 
const double MSTRANGE = 0.096; 
const double MCHARM = 1.27; 
const double MBOTTOM = 4.18; 
const double MTOP = 165.0; 
const double MELECTRON = 5.109989461e-4; 
const double MMUON = 1.0565835745e-1; 
const double MTAU = 1.77686; 
const double ALPHASMZ = 0.1181; 
const double ALPHAMZ = 1.0 / 127.950; 
 
const double PMTOP = 173.21; 
const double PMBOTTOM = 4.9; 
const double THETA12CKM = 0.229206; 
const double THETA13CKM = 0.003960; 
const double THETA23CKM = 0.042223; 
 
typedef enum {mUp=1, mCharm, mTop, mDown, mStrange, mBottom, mElectron,
              mMuon, mTau} mass;
typedef enum {ALPHA=1, ALPHAS} leGauge;
 
DoubleVector gaugeDerivs(double, const DoubleVector &);
 
  class QedQcd: public RGE, public Approx {
  private:
    DoubleVector a;   
    DoubleVector mf;  
    double mtPole, mbPole; 
    double mbMb;      
    double mtauPole;  
    DoubleVector aMz; 
    
public:
  QedQcd(); 
  QedQcd(const QedQcd &); 
  const QedQcd& operator=(const QedQcd & m); 
  virtual ~QedQcd() {};
  
  void setPoleMt(double mt) { mtPole = mt; }; 
  void setPoleMb(double mb) { mbPole = mb; }; 
  void setPoleMtau(double mtau) { mtauPole = mtau; }; 
  void setMbMb(double mb)   { mbMb = mb;   }; 
  void setMass(mass mno, double m) { mf(mno) = m; }; 
    void setAlpha(leGauge ai, double ap) { a(ai) = ap; };
    void setAlphaMz(leGauge ai, double ap) { aMz(ai) = ap; }; 
  void set(const DoubleVector &); 
  
  double displayPoleMt() const { return mtPole; };
  double displayPoleMtau() const { return mtauPole; };
  double displayPoleMb() const { return mbPole; };
  const DoubleVector & displayMass() const { return mf; };
  double displayMass(mass mno) const { return mf.display(mno); };
  double displayAlpha(leGauge ai) const { return a.display(ai); };
    double displayAlphaMz(leGauge ai) const { return aMz.display(ai); };
  const DoubleVector display() const;
  double displayMbMb() const { return mbMb; }
  
  int flavours(double) const;  
  
  double qedBeta() const;   
  double qcdBeta() const;   
  void massBeta(DoubleVector &) const; 
  DoubleVector beta() const; 
  
  void runGauge(double start, double end);
  double extractPoleMb(double asMb);
  double extractRunningMb(double asMb);
  void calcRunningMb();
  void calcPoleMb();
 
  void toMt();
  void toMz();
  // alpha1 is in the GUT normalisation. sinth = sin^2 thetaW(Q) in MSbar
  // scheme
  DoubleVector  getGaugeMu(const double m2, const
                     double sinth) const;
};
 
ostream & operator <<(ostream &, const QedQcd &);
istream & operator >>(istream &left, QedQcd &m);
 
double getRunMt(double poleMt, double asmt);
double getAsmt(double mtop, double alphasMz);
double getRunMtFromMz(double poleMt, double asMZ);
 
inline QedQcd::QedQcd(const QedQcd &m)
  : RGE(), Approx(m.displayApprox()), a(m.a), mf(m.mf), mtPole(m.mtPole), 
    mbPole(m.mbPole), mbMb(m.mbMb), 
    mtauPole(m.mtauPole), aMz(m.aMz) { 
  setPars(11); 
  setMu(m.displayMu());
}
 
void massFermions(const QedQcd & r, DoubleMatrix & mDon, 
                  DoubleMatrix & mUpq, DoubleMatrix & mEle);
 
} // namespace softsusy
 
#endif