LoopMediators.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2018 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef LOOPMEDIATORS_H
#define LOOPMEDIATORS_H
 
#include "StandardModel.h"
#include "gslpp.h"
#include "LoopMediatorsMatching.h"
#include "ThObservable.h"
 
class LoopMediators: public StandardModel {
public:
 
    static const int NLoopMediatorsvars = 10;
 
    static const std::string LoopMediatorsvars[NLoopMediatorsvars];
    
    LoopMediators();
    
    ~LoopMediators();
    
    virtual bool InitializeModel();
    
    virtual bool Init(const std::map<std::string, double>& DPars);
    
    virtual bool PreUpdate();
    
    virtual bool Update(const std::map<std::string, double>& DPars);
    
    virtual bool PostUpdate();
    
    virtual bool CheckParameters(const std::map<std::string, double>& DPars);
    
    virtual bool setFlag(const std::string name, const bool value);
    
    double F9(double x, double y);
    double F7(double x);
    double F7t(double x);
    double G7(double x);
    double G7t(double x);
    
    virtual LoopMediatorsMatching& getMatching() const
    {
        return LoopMediatorsM.getObj();
    }
    
    gslpp::complex getC1() const {
        return C1;
    }
    
    gslpp::complex getC2() const {
        return C2;
    }
    
    gslpp::complex getC3() const {
        return C3;
    }
    
    gslpp::complex getC4() const {
        return C4;
    }
    
    gslpp::complex getC5() const {
        return C5;
    }
    
    gslpp::complex getC1p() const {
        return C1p;
    }
    
    gslpp::complex getC2p() const {
        return C2p;
    }
    
    gslpp::complex getC3p() const {
        return C3p;
    }
    
    gslpp::complex getC7() const {
        return C7;
    }
    
    gslpp::complex getC8() const {
        return C8;
    }
    
    gslpp::complex getC9() const {
        return C9;
    }
    
    gslpp::complex getC10() const {
        return C10;
    }
    
    gslpp::complex getCS() const {
        return CS;
    }
    
    gslpp::complex getCP() const {
        return CP;
    }
    
    gslpp::complex getC7p() const {
        return C7p;
    }
    
    gslpp::complex getC8p() const {
        return C8p;
    }
    
    gslpp::complex getC9p() const {
        return C9p;
    }
    
    gslpp::complex getC10p() const {
        return C10p;
    }
    
    gslpp::complex getCSp() const {
        return CSp;
    }
    
    gslpp::complex getCPp() const {
        return CPp;
    }
    
    double getDeltaamu() const {
        return Deltaamu;
    }
    
    double getWCscale() const {
        return WCscale;
    }
    
protected: 
    
    virtual void setParameter(const std::string, const double&);
    mutable Matching<LoopMediatorsMatching,LoopMediators> LoopMediatorsM;
 
private:
    
    gslpp::complex C1;
    gslpp::complex C2;
    gslpp::complex C3;
    gslpp::complex C4;
    gslpp::complex C5;
    
    gslpp::complex C1p;
    gslpp::complex C2p;
    gslpp::complex C3p;
 
    gslpp::complex C7;
    gslpp::complex C8;
    gslpp::complex C9;
    gslpp::complex C10;
    gslpp::complex CS;
    gslpp::complex CP;
    
    gslpp::complex C7p;
    gslpp::complex C8p;
    gslpp::complex C9p;
    gslpp::complex C10p;
    gslpp::complex CSp;
    gslpp::complex CPp;
    
    double Deltaamu;
    
    double GammaL;
    double GammaR;
    double GammamuL;
    double GammamuR;
    double lambdaE;
    double mphi;
    double yD;
    double yE;
    double charge;
    
    double WCscale;
    
      
};
 
 
 
class Deltaamu : public ThObservable {
public:
   Deltaamu(const StandardModel& SM_i);
     
   ~Deltaamu();
 
    double computeThValue();
    const LoopMediators * myLM;
};
 
#endif /* LOOPMEDIATORS_H */