SUSYMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#ifndef SUSYMATCHING_H
#define SUSYMATCHING_H

#include <gslpp.h>
#include <complex>
#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();

    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 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, 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> TEhat;

    double tanb;
    double sinb;
    double cosb;
    double v;
    double v1;
    double v2;
    double gW;

    void NeutralinoRemixing();

    int delta_ab(int a, int b);
};

#endif  /* SUSYMATCHING_H */