EvolDF2.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#ifndef EVOLDF2_H
#define EVOLDF2_H

#include "RGEvolutor.h"
#include "StandardModel.h"
#include <gsl/gsl_sf_dilog.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;    
    
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;
    double Mz_cache;
};

#endif  /* EVOLDF2_H */