EvolDF2.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EVOLDF2_H
#define EVOLDF2_H
 
class StandardModel;
 
#include "RGEvolutor.h"
 
class EvolDF2 : public RGEvolutor {
public:
    EvolDF2(unsigned int dim_i, schemes scheme, orders order, const StandardModel& model);
    
    virtual ~EvolDF2();
    
    gslpp::matrix<double> AnomalousDimension(orders order, unsigned int nf, int basis = 0) const;
 
    gslpp::matrix<double>& Df2Evol(double mu, double M, orders order, 
            schemes scheme = NDR);
    
    double etacc(double mu) const;
    
    double etact(double mu) const;
    
    double etatt(double mu) const;    
    
    double etabS0(double mu) const;    
    
private:
    // c and d are the coefficient of als(mu) e als(M)
    // first index number of flavours
    // double b[5][5][5], c[3][5][5][5], d[3][5][5][5];
    
    //double S1tt() const;
    void Df2Evol(double mu, double M, double nf, schemes scheme);
    double a[5];
    double b[5][5][5];
    double c[3][5][5][5];
    double d[3][5][5][5];
    const StandardModel& model;
    unsigned int dim;
    double alsMZ_cache = 0.;
    double Mz_cache = 0.;
};
 
#endif  /* EVOLDF2_H */