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 "WilsonCoefficient.h"
#include <vector>
#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::matrix<double> RGevolP(int nf, orders order);
 
    gslpp::vector<double> ThresholdCp(orders order);
 
    double C_P(orders order);
 
    double C_P_qed(orders_qed order_qed);
 
    gslpp::matrix<double> RGevolB(int nf, orders order);
    gslpp::vector<double> ThresholdCb(orders order);
 
    gslpp::vector<double> C(orders order,orders_qed order_qed, int contribution);
   
    gslpp::matrix<double> RGevol_J(orders order, int nf, int contribution);
    
    gslpp::matrix<double> RGevol_R(orders_qed order_qed, int nf, int contribution);
 
 
    gslpp::vector<double> CWin_muw(orders order,orders_qed order_qed, int contribution);
 
 
    gslpp::matrix<double> ADM(orders order,orders_qed order_qed , double nf, int contribution);
 
 
    double C_Be(orders order);
 
    double C_Be_qed(orders_qed order_qed);
 
    double C_Bt(orders order);
 
    double C_Bt_qed(orders_qed order_qed);
 
    std::vector<WilsonCoefficient>& EVOCkpnn();
    
    
    double C_TOT(orders order, orders_qed order_qed);
 
private:
    const StandardModel& model;
    WilsonCoefficient evoCkpnn;
    gslpp::vector<double>  dcp,dcb;
    double mc_mc,etab,etacb,etac,kc,xc_mc_qed,L;
    double xi1c,xi2c,xice,xices;
    //box
    gslpp::vector<double> CW0b,CW1b,CW2b,CWeb,CWesb;
    gslpp::matrix<double> U0_4b,U0_5b,J1_4b,J2_4b,J1_5b,J2_5b,R0_4b,R0_5b,R1_4b,R1_5b;
    //penguin
    gslpp::vector<double> CW0p,CW1p,CW2p,CWep,CWesp;
    gslpp::matrix<double> U0_4p,U0_5p,J1_4p,J2_4p,J1_5p,J2_5p,R0_4p,R0_5p,R1_4p,R1_5p;
    
protected:
    std::vector<WilsonCoefficient> vevoCkpnn;
    
    
};
 
#endif /* CHARM_KPNUNU_H */