EWSMApproximateFormulae.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef EWSMAPPROXIMATEFORMULAE_H
#define EWSMAPPROXIMATEFORMULAE_H
 
#include <string>
#include <math.h>
#include "EWSMcache.h"
 
class EWSMApproximateFormulae {
public:
 
    EWSMApproximateFormulae(const EWSMcache& cache_i);
 
 
 
    double Mw() const;
 
    double sin2thetaEff(const Particle p) const
    {
        if (p.is("QUARK"))
            return sin2thetaEff_q((QCD::quark) (p.getIndex() - 6));
        else if (p.is("LEPTON"))
            return sin2thetaEff_l((QCD::lepton) p.getIndex());
        else
            throw std::runtime_error("EWSMApproximateFormulae::sin2thetaEff() called with wrong argument");
    }
 
    double DeltaR_TwoLoopEW_rem(const double Mw_i) const;
 
    double DeltaKappa_l_TwoLoopEW_rem(const double Mw_i) const;
 
    double DeltaKappa_b_TwoLoopEW_rem(const double Mw_i) const;
 
    double R0_bottom_OLD() const;
 
    double Gu_over_Gb_OLD() const;
 
    double Gd_over_Gb_OLD() const;
 
    double X(const std::string observable) const;
 
    double X_extended(const std::string observable) const;
    
    
    double X_full_2_loop(const std::string observable) const;
    
    
    double X_full(const std::string observable) const;
    
    
    double sin2thetaEff_b_full() const;
    
    double sin2thetaEff_l_full() const;
 
 
 
private:
    double sin2thetaEff_l(const QCD::lepton l) const;
 
    double sin2thetaEff_q(const QCD::quark q) const;
    
    double sin2thetaEff_b() const;
 
 
 
    const EWSMcache& mycache; 
 
};
 
#endif  /* EWSMAPPROXIMATEFORMULAE_H */