EvolDC1.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EVOLDC1_H
#define EVOLDC1_H
 
class StandardModel;
#include "RGEvolutor.h"
#include <sstream>
 
class EvolDC1 : public RGEvolutor {
    public :
    EvolDC1(unsigned int dim,  schemes scheme, orders order, const StandardModel& model);
    virtual ~EvolDC1();
    gslpp::matrix<double> AnomalousDimension_M(orders order, unsigned int n_u, unsigned int n_d) const;
    gslpp::matrix<double>& DC1Evol(double mu, double M, orders order,  schemes scheme = NDR);
    gslpp::matrix<double> ToRescaledBasis(orders order, unsigned int n_u, unsigned int n_d) const;
    gslpp::matrix<double> ToEffectiveBasis(gslpp::matrix<double> mat)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 DC1Evol(double mu, double M, double nf, schemes scheme);
    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  /* EVOLDC1_H */