EvolDF1.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EVOLDF1_H
#define EVOLDF1_H
 
#include <string>
#include "RGEvolutorNew.h"
#include "StandardModel.h"
#include "gslpp_special_functions.h"
#include <map>
#include "boost/multi_array.hpp"
 
#define NF 4
 
typedef unsigned int uint;
typedef unsigned int indices;
 
class EvolDF1 : public RGEvolutorNew {
public:
    EvolDF1(std::string reqblocks, schemes scheme, const StandardModel& model_i, qcd_orders order_qcd, qed_orders order_qed);
    virtual ~EvolDF1();
    gslpp::matrix<double> AnomalousDimension(indices nm, uint n_u, uint n_d) const;
 
    const Expanded<gslpp::matrix<double> >& DF1Evol(double mu, double M, schemes scheme = NDR);
 
    //gslpp::matrix<double> ToRescaleBasis(orders order, uint n_u, uint n_d) const;
    //gslpp::matrix<double> ToEffectiveBasis(gslpp::matrix<double> mat)const;
 
    //    std::map<std::string,uint> blocks_nops;
 
    //    {{"C",2},{"CP",6},{"CPM",8},{"L",2},{"CPML",10},{"CPQB",11},{"CPMQB",13},{"CPMLQB",15}};
    //    std::map<std::string,orders> blocks_ord;// = {{"C",NNLO},{"CP",NNLO},{"CPM",NNLO},{"L",NNLO},{"CPML",NNLO},{"CPQB",NLO},{"CPMQB",NLO},{"CPMLQB",NLO}};
 
 
private:
 
    //double a[4][13],
    //double b[4][13][13][13], c[4][13][13][13], d[4][13][13][13];
 
    //    typedef boost::multi_array<double, 4> array_type4;
    //    typedef boost::multi_array<double, 2> array_type2;
    //    array_type2 mn_a;
    //    array_type4 mn_b,mn_c,mn_d;
 
 
    void CheckNf(indices nm, uint nf) const;
 
    //    void DF1Evol(double mu, double M, double nf, schemes scheme);
    void DF1Ev(double mu, double M, int nf, schemes scheme);
 
    friend double gslpp_special_functions::zeta(int i);
 
    gslpp::matrix<double> GammaCC(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaCP(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaCM(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaCL(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaCQ(indices nm, uint n_u, uint n_d) const;
 
    gslpp::matrix<double> GammaPP(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaPM(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaPL(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaPQ(indices nm, uint n_u, uint n_d) const;
 
    gslpp::matrix<double> GammaMM(indices nm, uint n_u, uint n_d) const;
 
    gslpp::matrix<double> GammaLL(indices nm, uint n_u, uint n_d) const;
 
    gslpp::matrix<double> GammaQP(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaQM(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaQL(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaQQ(indices nm, uint n_u, uint n_d) const;
 
    gslpp::matrix<double> GammaBP(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaBL(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaBQ(indices nm, uint n_u, uint n_d) const;
    gslpp::matrix<double> GammaBB(indices nm, uint n_u, uint n_d) const;
 
    double f_f(uint nf, uint i, uint j, int k, double eta);
 
    double f_r(uint nf, uint i, uint j, int k, double eta);
 
    double f_g(uint nf, uint i, uint p, uint j, int k, int l, double eta);
 
    double f_h(uint nf, uint i, uint p, uint q, uint j, int k, int l, int m, double eta);
 
    std::map< uint, double > ai[NF];
    std::map< std::vector<uint>, double > vM0vi[NF], vM1vi[NF], vM2vi[NF], vM3vi[NF], vM4vi[NF], vM5vi[NF],
    vM6vi[NF], vM11vi[NF], vM33vi[NF], vM31vi[NF], vM13vi[NF], vM34vi[NF], vM43vi[NF], vM23vi[NF], vM32vi[NF],
    vM14vi[NF], vM41vi[NF], vM113vi[NF], vM131vi[NF], vM311vi[NF], vM133vi[NF], vM313vi[NF], vM331vi[NF];
 
    const StandardModel& model;
 
    // operators to include         {C, P, M, L, Q, b}            Huber et al., hep-ph/0512066
    uint nops, nfmin, nfmax;
    std::string blocks;
 
    gslpp::matrix<double> evec, evec_i, js, h, gg, s_s, jssv, jss, jv, vij;
    gslpp::vector<double> eval;
    gslpp::matrix<gslpp::complex> evecc;
    gslpp::vector<gslpp::complex> evalc;
    double alsM_cache, MAls_cache;
 
    //caching
#define F_iCacheSize 5
    int f_f_c[4][F_iCacheSize];
    double f_f_d[2][F_iCacheSize];
};
 
#endif /* EVOLDF1_H */