EWSUSY.h

Go to the documentation of this file.

/*
 * Copyright (C) 2013 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EWSUSY_H
#define EWSUSY_H
 
#include <gslpp.h>
#include <PVfunctions.h>
#include "StandardModel.h"
 
class SUSY;
 
class EWSUSY {
public:
 
    EWSUSY(const SUSY& SUSY_in);
 
    void SetRosiekParameters();
 
    gslpp::complex FA(const double mu, const double p2, const double mi, const double mj,
               const gslpp::complex cV_aij, const gslpp::complex cV_bji,
               const gslpp::complex cA_aij, const gslpp::complex cA_bji) const;
 
    gslpp::complex dFA(const double mu, const double p2, const double mi, const double mj,
                const gslpp::complex cV_aij, const gslpp::complex cV_bji,
                const gslpp::complex cA_aij, const gslpp::complex cA_bji) const;
    
    gslpp::complex PiT_Z(const double mu, const double p2, const double Mw_i) const;
 
    gslpp::complex PiT_AZ(const double mu, const double p2, const double Mw_i) const;
 
    gslpp::complex PiT_W(const double mu, const double p2, const double Mw_i) const;
 
    gslpp::complex PiTp_A(const double mu, const double p2, const double Mw_i) const;
 
    double PiThat_W_0(const double Mw_i) const;
 
    double DeltaR_rem_SM(const double Mw_i) const;
 
    double DeltaR_boxLL_SUSY(const double Mw_i) const;
 
    double DeltaR_boxLR_SUSY(const double Mw_i) const;
    
    gslpp::complex v(const double mu, const QCD::lepton M,
              const QCD::lepton J, const double Mw_i) const;
 
    gslpp::complex delta_v(const double mu, const QCD::lepton M,
                    const QCD::lepton J, const double Mw_i) const;
 
    double DeltaR_vertex_SUSY(const double Mw_i) const;
 
    gslpp::complex Sigma_nu_0(const double mu, const QCD::lepton I,
                       const QCD::lepton J, const double Mw_i) const;
 
    double DeltaR_neutrino_SUSY(const double Mw_i) const;
 
    double DeltaAlphaL5q_SM_EW1() const;
 
    double DeltaR_TOTAL_EW1(const double Mw_i) const;
 
    double DeltaR_SUSY_EW1(const double Mw_i) const;
 
    double Mw_MSSM_TMP(const double Mw_i) const;
 
    double Mw_MSSM() const;
 
    double getRenormalizationScaleFactor() const
    {
        return RenormalizationScaleFactor;
    }
 
private:
    static const double Mw_unphysical;
    static const double RenormalizationScaleFactor;
 
    const PVfunctions PV; 
    const SUSY& mySUSY; 
 
    gslpp::matrix<gslpp::complex> Yu; 
    gslpp::matrix<gslpp::complex> Yd; 
    gslpp::matrix<gslpp::complex> Yl; 
 
    gslpp::matrix<gslpp::complex> Au; 
    gslpp::matrix<gslpp::complex> Ad; 
    gslpp::matrix<gslpp::complex> Al; 
 
    gslpp::matrix<gslpp::complex> Zm; 
    gslpp::matrix<gslpp::complex> Zp; 
    gslpp::matrix<gslpp::complex> ZN; 
    gslpp::matrix<gslpp::complex> ZU; 
    gslpp::matrix<gslpp::complex> ZD; 
    gslpp::matrix<gslpp::complex> ZL; 
    gslpp::matrix<gslpp::complex> Zne; 
    gslpp::matrix<double> ZR; 
    gslpp::matrix<double> ZH; 
 
    double m_u[3], m_d[3], m_l[3], mH02[4], Msu2[6], Msd2[6], Mse2[6], Msn2[3], mC[2], mN[4];
 
    gslpp::complex L_esnC(const int N, const int K, const int i, const double Mw_i) const;
    gslpp::complex R_esnC(const int N, const int K, const int i) const;
    gslpp::complex L_nLC(const int I, const int k, const int i, const double Mw_i) const;
    gslpp::complex L_nsnN(const int J, const int K, const int j, const double Mw_i) const;
    gslpp::complex L_eLN(const int N, const int k, const int j, const double Mw_i) const;
    gslpp::complex R_eLN(const int N, const int k, const int j, const double Mw_i) const;
 
    gslpp::complex F(const double m1, const double m2, const double m3, const double m4) const;
 
    gslpp::complex H(const double m1, const double m2, const double m3, const double m4) const;
 
    gslpp::complex f(const double m1, const double m2, const double m3) const;
 
};
 
#endif  /* EWSUSY_H */