EWSMApproximateFormulae.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * 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((StandardModel::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;



private:
    double sin2thetaEff_l(const StandardModel::lepton l) const;

    double sin2thetaEff_q(const QCD::quark q) const;



    const EWSMcache& mycache; 

};

#endif  /* EWSMAPPROXIMATEFORMULAE_H */