HeffDF2.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef HEFFDF2_H
#define HEFFDF2_H
 
class StandardModel;
class EvolDF2;
#include "WilsonCoefficient.h"
#include "gslpp.h"
#include <memory>
 
class HeffDF2 {
public:
    HeffDF2(const StandardModel& SM);
    
    virtual ~HeffDF2();
    
    void ChangeScheme(schemes schout, WilsonCoefficient& c_in, orders order);
    
    gslpp::vector<gslpp::complex>** ComputeCoeffBd(double mu, schemes scheme = NDR);
 
    gslpp::vector<gslpp::complex>** ComputeCoeffBs(double mu, schemes scheme = NDR, bool SM = false);
 
    gslpp::vector<gslpp::complex>** ComputeCoeffdd(double mu, schemes scheme = NDR);
    
    gslpp::vector<gslpp::complex>** ComputeCoeffK(double mu, schemes scheme = NDR);
    
    gslpp::vector<gslpp::complex>** ComputeCoeffmK(double mu, schemes scheme = NDR);
    
    gslpp::matrix<double> AnomalousDimension(orders order, unsigned int nf = 0) const;
 
    WilsonCoefficient getCoeffBd() const {
        return coeffbd;
    }
 
    WilsonCoefficient getCoeffBs() const {
        return coeffbs;
    }
    
    WilsonCoefficient getCoeffDD() const {
        return coeffDd;
    }
    
    WilsonCoefficient getCoeffK() const {
        return coeffk;
    }
    
     WilsonCoefficient getCoeffmK() const {
        return coeffmk;
    }
    
    EvolDF2& getUDF2() const {
        return *evolDF2;
    }
 
 
private:
    gslpp::complex S0tt(double mu) const;
    const StandardModel& model;
    gslpp::matrix<double> drNDRLRI;
    WilsonCoefficient coeffbd;
    WilsonCoefficient coeffbs;
    WilsonCoefficient coeffDd;
    WilsonCoefficient coeffk;
    WilsonCoefficient coeffmk;
    
    std::unique_ptr<EvolDF2> evolDF2;
};
 
#endif  /* HEFFDF2_H */