172 return (-(cosa*cosa*mHh2) - mHl2*sina*sina + 2.0*m12_2*tanb + cosa*(mHl2 - mHh2)*sina*tanb)/2.;
189 return (2.0*m12_2 + sina*cosa*(mHl2 - mHh2) - (cosa*cosa*mHl2 + sina*sina*mHh2)*tanb)/(2.0*tanb);
208 return (mHh2*cosa*cosa+mHl2*sina*sina-m12_2*tanb)/(vev*vev*cosb*cosb);
227 return (mHh2*sina*sina+mHl2*cosa*cosa-m12_2/tanb)/(vev*vev*sinb*sinb);
247 return ((mHh2-mHl2)*cosa*sina+2.0*mHp2*sinb*cosb-m12_2)/(vev*vev*sinb*cosb);
264 return ((mA2-2.0*mHp2)*sinb*cosb+m12_2)/(vev*vev*sinb*cosb);
280 return (m12_2-mA2*sinb*cosb)/(vev*vev*sinb*cosb);
299 return (m12_2+(mHh2-mHl2)*cosa*sina)/(vev*vev*sinb*cosb);
318 return 3.0*( -mHl2*(cosa*cosa*cosa/sinb - sina*sina*sina/cosb)
319 +cos(bma)*m12_2*(cosa*cosa/(sinb*sinb) - sina*sina/(cosb*cosb)) )/vev;
339 return -cos(bma)*(sinb*cosb*m12_2
340 +sina*cosa*(sinb*cosb*(2.0*mHl2 + mHh2) - 3.0*m12_2))/(sinb*sinb*cosb*cosb*vev);
360 return -sin(bma)*(sinb*cosb*m12_2
361 -sina*cosa*(sinb*cosb*(mHl2 + 2.0*mHh2) - 3.0*m12_2))/(sinb*sinb*cosb*cosb*vev);
380 return -3.0*( (cosa*cosa*cosa*sinb + sina*sina*sina*cosb)*sinb*cosb*mHh2
381 -(cosa*cosa*sinb*sinb - cosb*cosb*sina*sina)*m12_2*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
401 return ( (cosa*cosb-sina*sinb)*(m12_2-sinb*cosb*mHl2)
402 +(mHl2-2.0*mA2)*sinb*sinb*cosb*cosb*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
422 return ( (cosa*sinb+cosb*sina)*(m12_2-sinb*cosb*mHh2)
423 +(mHh2-2.0*mA2)*sinb*sinb*cosb*cosb*cos(bma) )/(sinb*sinb*cosb*cosb*vev);
443 return ( (cosa*cosb-sina*sinb)*(m12_2-sinb*cosb*mHl2)
444 +(mHl2-2.0*mHp2)*sinb*sinb*cosb*cosb*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
464 return ( (cosa*sinb+cosb*sina)*(m12_2-sinb*cosb*mHh2)
465 +(mHh2-2.0*mHp2)*sinb*sinb*cosb*cosb*cos(bma) )/(sinb*sinb*cosb*cosb*vev);
479 return 0.25*(-mHl2-mHh2+(mHl2-mHh2)*cos(2.0*bma));
492 double sin2b=2.0*sinb*cosb;
497 double sin2a=2.0*sina*cosa;
501 return 0.125*((mHl2-mHh2)*(3.0*sin2a-sin2b*cos(2.0*bma)
502 -sin(2.0*bma)*(cosb*cosb-sinb*sinb))
503 +2.0*(8.0*m12_2-(mHl2+mHh2)*sin2b))/sin2b;
517 return 0.25*(mHh2-mHl2)*sin(2.0*bma);
532 return (mHl2+mHh2+(mHh2-mHl2)*cos(2.0*bma))/(2.0*vev*vev);
546 double sin2b=2.0*sinb*cosb;
547 double cos2b=cosb*cosb-sinb*sinb;
552 return (-128.0*m12_2*cos2b*cos2b
553 +4.0*(mHl2+mHh2)*(5.0+3.0*(cos2b*cos2b-sin2b*sin2b))*sin2b
554 +2.0*(mHl2-mHh2)*((6.0+10.0*(cos2b*cos2b-sin2b*sin2b))*cos(2.0*bma)
555 +16.0*sin2b*cos2b*sin(2.0*bma))*sin2b)
556 /(16.0*sin2b*sin2b*sin2b*vev*vev);
570 double sin2b=2.0*sinb*cosb;
571 double cos2b=cosb*cosb-sinb*sinb;
577 double sin2a=2.0*sina*cosa;
580 return ((mHh2-mHl2)*(3.0*sin2a-sin(2.0*bma)*cos2b-cos(2.0*bma)*sin2b)
581 +2.0*(-8.0*m12_2+(mHh2+mHl2+4.0*mHp2)*sin2b))/(4.0*sin2b*vev*vev);
598 return (2.0*mA2+mHl2+mHh2-4.0*mHp2+(mHl2-mHh2)*cos(2.0*bma))/(2.0*vev*vev);
614 return (-2.0*mA2+mHl2+mHh2+(mHl2-mHh2)*cos(2.0*bma))/(2.0*vev*vev);
629 return ((mHl2-mHh2)*sin(2.0*bma))/(2.0*vev*vev);
643 double sin2b=2.0*sinb*cosb;
644 double cos2b=cosb*cosb-sinb*sinb;
651 return (4.0*cos2b*(-4.0*m12_2+(mHl2+mHh2)*sin2b)
652 +(mHl2-mHh2)*(3.0*(cosa*cosa-sina*sina)
653 +cos(2.0*bma)*cos2b-sin(2.0*bma)*sin2b)*sin2b)
654 /(4.0*sin2b*sin2b*vev*vev);
667 double sin2b=2.0*sinb*cosb;
671 double sin2a=2.0*sina*cosa;
674 return -((mHh2+mHl2)*sin2b+(mHh2-mHl2)*sin2a-4.0*m12_2)/(4.0*sin2b);
699 double sin2b=2.0*sinb*cosb;
700 double cos2b=cosb*cosb-sinb*sinb;
705 return ((mHl2-mHh2)*sin(2.0*bma)-4.0*m12_2*cos2b)/(4.0*sin2b);
719 double sin2b=2.0*sinb*cosb;
720 double cos2b=cosb*cosb-sinb*sinb;
725 double sin2a=2.0*sina*cosa;
726 double cos2a=cosa*cosa-sina*sina;
729 return ((mHh2+mHl2-(mHh2-mHl2)*(cos2a*cos2b-sin2a*sin2b)
730 -4.0*m12_2/sin2b)/(sin2b*sin2b)
731 +2.0*mHp2)/(4.0*vev*vev);
746 double sin2b=2.0*sinb*cosb;
747 double cos2b=cosb*cosb-sinb*sinb;
751 double cos2a=cosa*cosa-sina*sina;
754 return ((mHh2+mHl2)*(0.5*tanb-0.5/tanb)
755 +(mHh2-mHl2)*cos2a/sin2b
756 +4.0*m12_2*cos2b/(sin2b*sin2b))/(4.0*vev*vev*sin2b);
769 double sin2b=2.0*sinb*cosb;
773 return (m12_2/sin2b-0.5*mHp2)/(vev*vev);
787 return (mA2-mHp2)/(2.0*vev*vev);
801 double sin2b=2.0*sinb*cosb;
802 double cos2b=cosb*cosb-sinb*sinb;
808 return ((mHh2+mHl2-(mHh2-mHl2)*cos(2.0*bma)
809 -4.0*m12_2*cos2b*cos2b/sin2b)/(sin2b*sin2b)
810 -2.0*mHp2)/(4.0*vev*vev);
826 double sin2b=2.0*sinb*cosb;
827 double cos2b=cosb*cosb-sinb*sinb;
831 double sin2a=2.0*sina*cosa;
832 double cos2a=cosa*cosa-sina*sina;
835 return (mHh2+mHl2+4.0*m12_2*cos2b/sin2b
836 +(mHh2-mHl2)*(cos2a-sin2a/tanb))/(4.0*vev*vev*cosb*cosb);
851 double sin2b=2.0*sinb*cosb;
852 double cos2b=cosb*cosb-sinb*sinb;
856 double sin2a=2.0*sina*cosa;
857 double cos2a=cosa*cosa-sina*sina;
860 return (mHh2+mHl2-4.0*m12_2*cos2b/sin2b
861 -(mHh2-mHl2)*(cos2a+sin2a*tanb))/(4.0*vev*vev*sinb*sinb);
875 double sin2b=2.0*sinb*cosb;
879 double sin2a=2.0*sina*cosa;
882 return ((mHh2-mHl2)*sin2a-2.0*m12_2)/(2.0*vev*vev*sin2b);
895 return (2.0*mHp2)/(vev*vev);
908 double sin2b=2.0*sinb*cosb;
911 return 4.0*m12_2/(vev*vev*sin2b);
924 return 2.0*mA2/(vev*vev);
double computeThValue()
The quartic coupling .
E11_THDM(const StandardModel &SM_i)
E11_THDM constructor.
double computeThValue()
The quartic coupling .
E22_THDM(const StandardModel &SM_i)
E22_THDM constructor.
E33_THDM(const StandardModel &SM_i)
E33_THDM constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
HHlambda1(const StandardModel &SM_i)
HHlambda1 constructor.
double computeThValue()
The quartic coupling .
HHlambda2(const StandardModel &SM_i)
HHlambda2 constructor.
double computeThValue()
The quartic coupling .
HHlambda3(const StandardModel &SM_i)
HHlambda3 constructor.
HHlambda4(const StandardModel &SM_i)
HHlambda4 constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
HHlambda5(const StandardModel &SM_i)
HHlambda5 constructor.
HHlambda6(const StandardModel &SM_i)
HHlambda6 constructor.
double computeThValue()
The quartic coupling .
A model class for the Standard Model.
const double v() const
The Higgs vacuum expectation value.
A base class for symmetric Two-Higgs-Doublet models.
double getmHp() const
A method get the charged Higgs mass.
double getmHh() const
A method get the mass of the "non-125 GeV" neutral scalar Higgs.
double gettanb() const
A method get .
double getbma() const
A method get .
double getmHh2() const
A method get the squared mass of the "non-125 GeV" neutral scalar Higgs.
double getsinb() const
A method get .
double getmHl2() const
A method get the squared mass of the lighter neutral scalar Higgs.
double getmHp2() const
A method get the squared charged Higgs mass.
double getmA() const
A method get the mass of the pseudoscalar Higgs A.
double getm12_2() const
A method get the parameter of the Higgs potential .
double getcosb() const
A method get .
double getcosa() const
A method get .
double getsina() const
A method get .
double getmA2() const
A method get the squared mass of the pseudoscalar Higgs A.
A class for a model prediction of an observable.
Y1_THDM(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quadratic coupling .
double computeThValue()
The quadratic coupling .
Y2_THDM(const StandardModel &SM_i)
Constructor.
Y3_THDM(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quadratic coupling .
double computeThValue()
The quartic coupling .
Z1_THDM(const StandardModel &SM_i)
Z1_THDM constructor.
double computeThValue()
The quartic coupling .
Z2_THDM(const StandardModel &SM_i)
Z2_THDM constructor.
double computeThValue()
The quartic coupling .
Z3_THDM(const StandardModel &SM_i)
Z3_THDM constructor.
double computeThValue()
The quartic coupling .
Z4_THDM(const StandardModel &SM_i)
Z4_THDM constructor.
Z5_THDM(const StandardModel &SM_i)
Z5_THDM constructor.
double computeThValue()
The quartic coupling .
Z6_THDM(const StandardModel &SM_i)
Z6_THDM constructor.
double computeThValue()
The quartic coupling .
Z7_THDM(const StandardModel &SM_i)
Z7_THDM constructor.
double computeThValue()
The quartic coupling .
eta00_THDM(const StandardModel &SM_i)
eta00_THDM constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
eta3_THDM(const StandardModel &SM_i)
eta3_THDM constructor.
g_HhAA(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
double computeThValue()
The triple Higgs coupling .
g_HhHhHh(const StandardModel &SM_i)
Constructor.
g_HhHpHm(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
g_hAA(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
double computeThValue()
The triple Higgs coupling .
g_hHhHh(const StandardModel &SM_i)
Constructor.
g_hHpHm(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
g_hhHh(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
g_hhh(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
lambda1(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
lambda2(const StandardModel &SM_i)
Constructor.
lambda345(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling combination .
lambda3(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
lambda4(const StandardModel &SM_i)
Constructor.
lambda5(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
m11_2(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
m22_2(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
mass_mA(const StandardModel &SM_i)
mass_mA constructor.
mass_mHh(const StandardModel &SM_i)
mass_mHh constructor.
mass_mHl(const StandardModel &SM_i)
mass_mHl constructor.
mass_mHp(const StandardModel &SM_i)
mass_mHp constructor.
massdifference_mAmmHh(const StandardModel &SM_i)
massdifference_mAmmHh constructor.
massdifference_mAmmHl(const StandardModel &SM_i)
massdifference_mAmmHl constructor.
massdifference_mAmmHp(const StandardModel &SM_i)
massdifference_mAmmHp constructor.
massdifference_mHhmmA(const StandardModel &SM_i)
massdifference_mHhmmA constructor.
massdifference_mHhmmHp(const StandardModel &SM_i)
massdifference_mHhmmHp constructor.
massdifference_mHlmmA(const StandardModel &SM_i)
massdifference_mHlmmA constructor.
massdifference_mHlmmHp(const StandardModel &SM_i)
massdifference_mHlmmHp constructor.
massdifference_mHpmmA(const StandardModel &SM_i)
massdifference_mHpmmA constructor.
massdifference_mHpmmHh(const StandardModel &SM_i)
massdifference_mHpmmHh constructor.
massdifference_mHpmmHl(const StandardModel &SM_i)
massdifference_mHpmmHl constructor.
tanbeta(const StandardModel &SM_i)
tanbeta constructor.
double computeThValue()
The quadratic coupling .
xi0_THDM(const StandardModel &SM_i)
xi0_THDM constructor.
double computeThValue()
The quadratic coupling .
xi1_THDM(const StandardModel &SM_i)
xi1_THDM constructor.
double computeThValue()
The quadratic coupling .
xi3_THDM(const StandardModel &SM_i)
xi3_THDM constructor.