THDM.h

Go to the documentation of this file.

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

#ifndef THDM_H
#define THDM_H

#include "StandardModel.h"
#include "THDMMatching.h"
#include "THDMcache.h"

//class THDMcache; //forward reference to THDMcache class

class THDM: public StandardModel {
public:

    static const int NTHDMvars = 7;
    static const std::string THDMvars[NTHDMvars];
    
    THDM();
    
    ~THDM();
    
    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);

    
    // Flags

    virtual bool setFlagStr(const std::string name, const std::string value);

    virtual THDMMatching* getMyMatching() const
    {
        return myTHDMMatching;
    }

    
    
    double getv1() const {
        return v() * cosb;
    }

    double getv2() const {
        return v() * sinb;
    }
    

    std::string getModelTypeflag() const {
        return flag_model;
    }

    double getlogtb() const {
        return logtb;
    }

    double gettanb() const {
        return tanb;
    }

    double getsinb() const {
        return sinb;
    }

    double getcosb() const {
        return cosb;
    }

    double getbma() const {
        return bma;
    }

    double getsin_ba() const {
        return sin_ba;
    }

    double getcosa() const{
        return cos(atan(pow(10.,logtb))-bma);
    }

    double getsina() const{
        return sin(atan(pow(10.,logtb))-bma);
    }

    double getmHh2() const {
        return mHh2;
    }

    double getmHh() const {
        return sqrt(mHh2);
    }

    double getmA2() const {
        return mHh2-mHh2mmA2;
    }

    double getmA() const {
        return sqrt(mHh2-mHh2mmA2);
    }

    double getmHp2() const {
        return mHh2-mHh2mmHp2;
    }

    double getmHp() const {
        return sqrt(mHh2-mHh2mmHp2);
    }

    double getm12_2() const {
        return m12_2;
    }

    double getbsgamma_theoryerror() const {
        return bsgamma_theoryerror;
    }

protected: 

    virtual void setParameter(const std::string, const double&);
    THDMcache * mycache;

    bool checkmodelType(const std::string modeltype) const
    {
        if (modeltype.compare("type1") == 0
                || modeltype.compare("type2") == 0
                || modeltype.compare("typeX") == 0
                || modeltype.compare("typeY") == 0)
            return true;
        else
            return false;
    }

private:

    THDMMatching* myTHDMMatching;

    double logtb, tanb, sinb, cosb, bma, sin_ba, mHh2, mHh2mmA2, mHh2mmHp2, m12_2, bsgamma_theoryerror;
    std::string flag_model;
};

#endif  /* THDM_H */