SUSYMassInsertionMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef SUSYMASSINSERTIONMATCHING_H
#define SUSYMASSINSERTIONMATCHING_H
 
#include <gslpp.h>
#include <gsl/gsl_sf_dilog.h>
#include <stdexcept>
#include <StandardModelMatching.h>
 
class SUSYMassInsertion;
 
class SUSYMassInsertionMatching: public StandardModelMatching{
public :
    SUSYMassInsertionMatching(const SUSYMassInsertion& SUSYMassInsertion_i);
    
    
    //delta F=2 loop functions
    
    double f6(double x);
    double f6t(double x);
    
    double C0LO(double x);
    
    double C1LO(double x);
    
    double C2LO(double x);
    
    double C3LOA(double x);
    
    double C3LOB(double x);
    
    double C4LOA(double x);
    
    double C4LOB(double x);
    
    double C0NLO(double x, double mumatch2, double Ms2);
    
    double C1NLO(double x, double mumatch2, double Ms2);
    
    double C2NLO(double x, double mumatch2, double Ms2);
    
    double C3NLOA(double x, double mumatch2, double Ms2);
    
    double C3NLOB(double x, double mumatch2, double Ms2);
    
    double C4NLOA(double x, double mumatch2, double Ms2);
    
    double C4NLOB(double x, double mumatch2, double Ms2);
        
    //delta F=1 loop functions
    
    double B1(double x) const;
    
    double B2(double x) const;
    
    double P1(double x) const;
    
    double P2(double x) const;
    
    double M1(double x) const;
    
    double M2(double x) const;
    
    double M3(double x) const;
    
    double M4(double x) const;
    
    
    // Wilson coefficients DeltaF=2
    
    virtual std::vector<WilsonCoefficient>& CMdd2();
    
    virtual std::vector<WilsonCoefficient>& CMdbd2();
    
    virtual std::vector<WilsonCoefficient>& CMdbs2();
    
    virtual std::vector<WilsonCoefficient>& CMdk2();
    
    virtual std::vector<WilsonCoefficient>& CMd1();
    
    gslpp::matrix<double> RtoMisiak() const;
    
    gslpp::matrix<double> EffectiveBase() const;
    
    void LRItoNDR (int i);
    
private:
    const SUSYMassInsertion& SusyMI;
    double MuM2, Ms2;
    unsigned int Nf;
    gslpp::complex DLL, DLR, DRL, DRR;
    gslpp::matrix<double> drNDRLRI;
    
    WilsonCoefficient mcd2, mcd1, mcbd, mcbs, mck2;
    
};
 
#endif  /* SUSYMASSINSERTIONMATCHING_H */