CPenguinBox.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef CPENGUINBOX_H
#define CPENGUINBOX_H
 
 
class StandardModel;
#include "OrderScheme.h"
#include "gslpp.h"
#include <sstream>
#include <gsl/gsl_sf_dilog.h>
#include <gsl/gsl_sf_zeta.h>
#include <gsl/gsl_sf_clausen.h>
 
class CPenguinBox {
public:
 
    CPenguinBox(const StandardModel& model_i);
 
    ~CPenguinBox();
 
    gslpp::matrix<double> RGevolP(int nf);
 
    gslpp::matrix<double> RGevolB(int nf);
 
    gslpp::vector<double> WMatchP(orders order);
 
    gslpp::vector<double> WMatchB(orders order);
 
    double Cmatch(orders order);
 
    double CT_tot(orders order, orders_qed order_qed);
 
protected:
 
private:
    const StandardModel& model;
 
    gslpp::matrix<double> up, j5p1, j4p1, j5p2, j4p2;
    gslpp::matrix<double> ub, j5b1, j4b1, j5b2, j4b2;
 
    gslpp::matrix<double> up4, up5, ub4, ub5;
 
    gslpp::vector<double> p, b;
 
    gslpp::vector<double> p0, p1, p2;
    gslpp::vector<double> b0, b1, b2;
 
    gslpp::vector<double> r1, r2;
 
    gslpp::vector<double> CPM0, CPM1, CPM2;
    gslpp::vector<double> CBM0, CBM1, CBM2;
 
    double r1tau, r2tau, r1e, r2e;
    //double tbp, tbb;
    double eta;
 
    double cpm, cbmt, cbme;
 
    double BmatchP(orders order);
 
    double BmatchB(orders order);
};
 
#endif /* CPENGUINBOX_H */