positivity.cpp

Go to the documentation of this file.

/* 
 * Copyright (C) 2015 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#include "positivity.h"
#include "THDM.h"
#include "THDMquantities.h"
 
positivity::positivity(const StandardModel& SM_i)
: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
{
    mylambda1=new lambda1(SM_i);
    mylambda2=new lambda2(SM_i);
};
 
positivity::~positivity()
{
  delete mylambda1;
  delete mylambda2;
};
 
double positivity::computeThValue()
{
    return 0.0;
 
}
 
positivity1::positivity1(const StandardModel& SM_i)
: positivity(SM_i)
{}
 
double positivity1::computeThValue()
{
    double mHl2=myTHDM->getmHl2();
    double mHh2=myTHDM->getmHh2();
    double mHp2=myTHDM->getmHp2();
    double sina=myTHDM->getsina();
    double cosa=myTHDM->getcosa();
    double tanb=myTHDM->gettanb();
    double sinb=myTHDM->getsinb();
    double cosb=myTHDM->getcosb();
    double m12_2=myTHDM->getm12_2();
    double lambda1=mylambda1->computeThValue();
    double lambda2=mylambda2->computeThValue();
 
    if (lambda1>0 && lambda2>0)
    {
        return 0.0001*(2.0*mHp2*cosb*sinb -m12_2 +(mHh2-mHl2)*cosa*sina
                       +sqrt(fabs((mHl2*cosa*cosa +mHh2*sina*sina -m12_2/tanb)
                             *(mHh2*cosa*cosa +mHl2*sina*sina -m12_2*tanb))));
    }
    else
    {
        return -1.;
    }
}
 
positivity2::positivity2(const StandardModel& SM_i)
: positivity(SM_i)
{}
 
double positivity2::computeThValue()
{
    double mHl2=myTHDM->getmHl2();
    double mA2=myTHDM->getmA2();
    double mHh2=myTHDM->getmHh2();
    double sina=myTHDM->getsina();
    double cosa=myTHDM->getcosa();
    double tanb=myTHDM->gettanb();
    double sinb=myTHDM->getsinb();
    double cosb=myTHDM->getcosb();
    double m12_2=myTHDM->getm12_2();
    double lambda1=mylambda1->computeThValue();
    double lambda2=mylambda2->computeThValue();
 
    if (lambda1>0 && lambda2>0)
    {
        return 0.0001*(mA2 -fabs(mA2-m12_2/cosb/sinb)
                       +((mHh2-mHl2)*cosa*sina
                         +sqrt(fabs((mHl2*cosa*cosa +mHh2*sina*sina -m12_2/tanb)
                               *(mHh2*cosa*cosa +mHl2*sina*sina - m12_2*tanb))))
                        /cosb/sinb);
    }
    else
    {
        return -1.;
    }
}