RGEvolutorNew.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef RGEVOLUTORNEW_H
#define RGEVOLUTORNEW_H
 
#include <gslpp.h>
#include "OrderScheme.h"
#include "WilsonTemplateNew.h"
 
class RGEvolutorNew : public WilsonTemplateNew<gslpp::matrix<double> > {
public:
    
    RGEvolutorNew(unsigned int dim, schemes scheme, qcd_orders order_qcd_i, qed_orders order_qed_i = QED0);
    
    virtual ~RGEvolutorNew() {};
    
    void setEvol(unsigned int i, unsigned int j, double x, qcd_orders order_i, qed_orders order_qed = QED0);
    
    void setEvol(const gslpp::matrix<double>& m, qcd_orders order_qcd_i, qed_orders order_qed_i = QED0);
    
    const Expanded<gslpp::matrix<double> >& getEvol() const;
 
    double getM() const;
 
    void setScales(double mu, double M);
 
    void setM(double M);
    
    void setMu(double mu);
 
    const gslpp::matrix<double>& Evol(qcd_orders order_qcd, qed_orders order_qed = QED0) const;
    
protected:
    double M;
};
 
#endif  /* RGEVOLUTOR_H */