EvolDC1Buras.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EVOLDC1BURAS_H
#define EVOLDC1BURAS_H
 
class StandardModel;
#include "RGEvolutor.h"
#include <sstream>
 
class EvolDC1Buras : public RGEvolutor {
    public:
    EvolDC1Buras(unsigned int dim_i,  schemes scheme, orders order, const StandardModel& model);
    virtual ~EvolDC1Buras();
    gslpp::matrix<double>& DC1EvolBuras(double mu, double M, orders order,  schemes scheme = NDR);
    gslpp::matrix<double> AnomalousDimension_DC1_Buras(orders order, unsigned int n_u,
        unsigned int n_d) const;
    gslpp::matrix<double> StrongThresholds() const;
    
    private:
    int nu,nd;    
    double a[3][10], b[3][10][10][10], c[3][10][10][10], d[3][10][10][10];
    const StandardModel& model;
    void DC1EvolBuras(double mu, double M, double nf, schemes scheme);
    void DC1PenguinThresholds(double M, orders order);
    gslpp::matrix<gslpp::complex> v, vi, js, h, gg, s_s, jssv, jss, jv, vij;
    gslpp::vector<gslpp::complex> e;
    unsigned int dim;
    double alsMZ_cache;
    double Mz_cache;
    
 };
 
 
 
 
#endif  /* EVOLDC1BURAS_H */