SUSYMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef SUSYMATCHING_H
#define SUSYMATCHING_H
 
#include "gslpp.h"
#include "Polylogarithms.h"
#include "PVfunctions.h"
#include "StandardModelMatching.h"
 
#define LEPS 1.e-5     // tolerance in the limit of S(x,y) to S(x) 
#define SUSYLEPS 5.e-3 // tolerance in the limits of D0[x,y,z,t] and D2[x,y,z,t]
#define SUSYLEPS2 1.e-10 // tolerance in the limits of B0[x,y]
#define SUSYLEPS3 5.e-3 // tolerance in the limits of D0N
 
class SUSY;
 
class SUSYMatching : public StandardModelMatching {
public:
    SUSYMatching(const SUSY & SUSY_i);
 
    virtual std::vector<WilsonCoefficient>& CMDLij(int li_lj);
 
    virtual std::vector<WilsonCoefficient>& CMDLi3j(int li_lj);
 
    virtual std::vector<WilsonCoefficient>& CMmueconv();
 
    virtual std::vector<WilsonCoefficient>& CMgminus2mu();
 
//   /** Calculates gamma penguin amplitudes for m->(3)e (1), t->(3)m (2) and t->(3)e (3)**/
    virtual gslpp::vector<gslpp::complex> AFunctions(int n);
    
    //   /** Calculates Z penguin amplitudes for m->3e (1), t->3m (2) and t->3e (3)**/
    virtual gslpp::vector<gslpp::complex> FFunctions(int n);
 
 
//   /** Calculates box amplitudes for m->3e (1), t->3m (2) and t->3e (3) **/
    virtual gslpp::vector<gslpp::complex> BFunctions(int n);
 
//   /** Calculates Higgs penguin amplitudes for m->3e (1), t->3m (2) and t->3e (3)**/
    virtual gslpp::vector<gslpp::complex> BHFunctions(int n);
 
//   /** Calculates box amplitudes for m->e conversion**/
    virtual gslpp::vector<gslpp::complex> DFunctions();
 
//   /** Calculates the muon g-2**/
    virtual gslpp::vector<gslpp::complex> gminus2mu();
 
//   /** Calculates the muon g-2**/
    virtual double gminus2muNLO();
 
    virtual gslpp::vector<gslpp::complex> C7_Lepton(int n);
 
    virtual gslpp::vector<gslpp::complex> C9_Lepton(int n);
 
    virtual gslpp::vector<gslpp::complex> C10_Lepton(int n);
 
    virtual gslpp::vector<gslpp::complex> CS_Lepton(int n);
 
    virtual gslpp::vector<gslpp::complex> CP_Lepton(int n);
 
    virtual gslpp::vector<gslpp::complex> CT_Lepton(int n);
    
    
    void updateSUSYParameters();
 
private:
    const SUSY & mySUSY;
    const Polylogarithms Polylogs;
    const PVfunctions PV;
 
    WilsonCoefficient mcDLij,mcDLi3j,mcmueconv,mcgminus2mu;
    
    gslpp::vector<double> mym_su_sq,mym_sd_sq,mym_se_sq,mym_sn_sq;
    
    gslpp::matrix<gslpp::complex> myRu,myRd,myRl,myRn;
    
    gslpp::vector<double> MChi;
    gslpp::matrix<gslpp::complex> myV;
    gslpp::matrix<gslpp::complex> myU;
 
    gslpp::vector<double> MChi0,MNeig;
    gslpp::matrix<gslpp::complex> myN,ON;
    
    gslpp::matrix<double> Lepty,Leptz,Leptfa1,Leptfa2,Leptf1,Leptf2,Leptf3,Leptf4;
    std::complex<double> Leptfzn[6][4][4],Leptgzn[6][4][4],Leptfzc[3][2][2],Leptgzc[3][2][2];
 
    gslpp::matrix<gslpp::complex> CRlE,CRlMU,CRlTAU,CRqUP,CRqDOWN,CLlE,CLlMU,CLlTAU,CLqUP,CLqDOWN,NRlE,NRlMU,NRlTAU,NRqUP,NRqDOWN,NLlE,NLlMU,NLlTAU,NLqUP,NLqDOWN;
    
    gslpp::matrix<gslpp::complex> AmpA1LN,AmpA1RN,AmpA1LC,AmpA1RC,AmpTauA1LN,AmpTauA1RN,AmpTauA1LC,AmpTauA1RC,AmpTEA1LN,AmpTEA1RN,AmpTEA1LC,AmpTEA1RC,AmpALN,AmpARN,AmpALC,AmpARC,AmpTauALN,AmpTauARN,AmpTauALC,AmpTauARC,AmpTEALN,AmpTEARN,AmpTEALC,AmpTEARC;
 
 
    gslpp::matrix<gslpp::complex> TUhat;
 
    gslpp::matrix<gslpp::complex> TDhat;
 
    gslpp::matrix<gslpp::complex> TEhat;
 
    double Q_S;
    double mu2R;
    double Als;
    double Mg;
    double MHpm;
    double tanb;
    double sinb;
    double cosb;
    double v;
    double v1;
    double v2;
    double gW,mW,sinthetaW;
    
    gslpp::matrix<gslpp::complex> upQmassM,downQmassM;
    
    void NeutralinoRemixing();
 
    int delta_ab(int a, int b);
 
        // loopfunctions for Barr-Zee corrections:
        double fPS(double x);
        double fS(double x);
        double fft(double x);
        // loopfunction for tanb-enhanced correction
        double It(double a, double b, double c);
        // loopfunctions for the photonic 2loop corrections:
        double F3C(double x);
        double F4C(double x);
        double F3N(double x);
        double F4N(double x);
        double Li2(double x);
        // loopfunctions for the following correction:
        //   Deviation of the smuon-muon-chargino/neutralino couplings
        //   from gauge/Yukawa couplings due to the breaking of SUSY relations
        //   ref. arXiv:1311.1775
        double Fa(double x, double y);
        double Fb(double x, double y);    
 
    gslpp::vector<double> mySUSYMQ;
    
};
 
#endif  /* SUSYMATCHING_H */