EvolDF1nlep.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EVOLDF1NLEP_H
#define EVOLDF1NLEP_H
 
 
class StandardModel;
#include "RGEvolutor.h"
#include <sstream>
 
class EvolDF1nlep : public RGEvolutor {
    public:
    EvolDF1nlep(unsigned int dim,  schemes scheme,orders order, orders_qed order_qed, const StandardModel& model);
    virtual ~EvolDF1nlep();
    gslpp::matrix<double> AnomalousDimension_nlep_S(orders order, unsigned int n_u, unsigned int n_d) const;
    gslpp::matrix<double> AnomalousDimension_nlep_EM(orders order, unsigned int n_u, unsigned int n_d) const;
    gslpp::matrix<double>& Df1Evolnlep(double mu, double M, orders order, orders_qed order_qed, schemes scheme = NDR);
    // evolutor involving only 3 flavors and allowing mu > M regardless thresholds
    gslpp::matrix<double>& Df1Evolnlep3flav(double mu, double M, orders order, orders_qed order_qed, schemes scheme = NDR);
    gslpp::matrix<double> Df1threshold_deltarsT(double nf) const;
    gslpp::matrix<double> Df1threshold_deltareT(double nf) const;    
    
    private:
    //int nu,nd;
    double a[4][10], b[4][10][10][10], c[4][10][10][10], d[4][10][10][10],
           m[4][10][10][10], n[4][10][10][10], mn[4][10][10][10], o[4][10][10][10], 
           p[4][10][10][10], op[4][10][10][10], q[4][10][10][10], qq[4][10][10][10], 
           r[4][10][10][10], rr[4][10][10][10], s[4][10][10][10], t[4][10][10][10];
    const StandardModel& model;
    void Df1Evolnlep(double mu, double M, double nf, schemes scheme);
    void Df1threshold_nlep(double M, double nf);
    gslpp::matrix <gslpp::complex>  V, Vi, gs, Js, ge0, K0, ge11, K11, JsK0V, ViK0Js,
                                    Gamma_s0T, Gamma_s1T, Gamma_eT, Gamma_seT, JsV, ViJs,
                                    K0V, ViK0, ge0sing, K0sing, K0singV, K11V, ViK11, ge11sing, 
                                    K11sing, K11singV, JsK0singV;
    gslpp::vector<gslpp::complex> e;
    unsigned int dim;
    double alsMZ_cache;
    double Mz_cache;
    double Ale_cache;
 
 };
 
 
#endif  /* EVOLDF1NLEP_H */