Charm_Kpnunu.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef CHARM_KPNUNU_H
#define CHARM_KPNUNU_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 Charm_Kpnunu {
public:
 
    Charm_Kpnunu(const StandardModel& model_i);
 
    ~Charm_Kpnunu();
 
    gslpp::vector<double> Cp(orders order);
 
    gslpp::matrix<double> RGevolP(int nf, orders order);
 
    gslpp::vector<double> ThresholdCp(orders order);
 
    gslpp::vector<double> C_p(orders order);
 
    double C_P(orders order);
 
    gslpp::vector<double> Cb(orders order);
 
    gslpp::matrix<double> RGevolB(int nf, orders order);
    gslpp::vector<double> ThresholdCb(orders order);
 
    gslpp::vector<double> C_b(orders order);
 
    double C_Be(orders order);
 
    double C_Bt(orders order);
 
    double P_C(orders order);
 
    double C_TOT(orders order, orders_qed order_qed);
 
private:
    const StandardModel& model;
    gslpp::vector<double> cp, dcp, c_p, cpmuW0, cpmuW1, cpmuW2, cb, dcb, c_b, cbmuW0,
    cbmuW1, cbmuW2;
    gslpp::matrix<double> U4p, U5p, J5p1, J4p1, J5p2, J4p2, dc_p,
    U4b, U5b, J5b1, J4b1, J5b2, J4b2, dc_b;
    double etab, etacb, etac, mc, kc, xi1c, xi2c, xc, CP, CBe, CBt;
};
 
#endif /* CHARM_KPNUNU_H */