a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Member Functions | Public Attributes | Private Member Functions | Private Attributes | List of all members
GeneralTHDMZ2Runner Class Reference

An RGE running algorithm for the GeneralTHDMZ2 parameters. More...

#include <GeneralTHDMZ2Runner.h>

Detailed Description

An RGE running algorithm for the GeneralTHDMZ2 parameters.

Author
HEPfit Collaboration
Copyright
GNU General Public License

Renormalization group evolution of the relevant SM and GeneralTHDMZ2 parameters

Definition at line 21 of file GeneralTHDMZ2Runner.h.

Public Member Functions

double betalambda1_Z2 (double la1, double la3, double la4, double la5, double Yb1, double Ytau1)
 The beta function of \(\lambda_1\) appearing in unitarity conditions. More...
 
double betalambda2_Z2 (double la2, double la3, double la4, double la5, double Yt, double Yb2, double Ytau2)
 The beta function of \(\lambda_2\) appearing in unitarity conditions. More...
 
double betalambda3_Z2 (double la1, double la2, double la3, double la4, double la5, double Yt, double Yb1, double Yb2, double Ytau1, double Ytau2)
 The beta function of \(\lambda_3\) appearing in unitarity conditions. More...
 
double betalambda4_Z2 (double la1, double la2, double la3, double la4, double la5, double Yt, double Yb1, double Yb2, double Ytau1, double Ytau2)
 The beta function of \(\lambda_4\) appearing in unitarity conditions. More...
 
double betalambda5_Z2 (double la1, double la2, double la3, double la4, double la5, double Yt, double Yb1, double Yb2, double Ytau1, double Ytau2)
 The beta function of \(\lambda_5\) appearing in unitarity conditions. More...
 
 GeneralTHDMZ2Runner (const StandardModel &SM_i)
 GeneralTHDMZ2Runner constructor. More...
 
gslpp::matrix< double > getGTHDMZ2_at_Q ()
 The public function which contains all relevant GTHDMZ2 parameter after running. More...
 
virtual double RGEGeneralTHDMZ2Runner (double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order)
 
virtual ~GeneralTHDMZ2Runner ()
 Runner destructor. More...
 

Public Attributes

const GeneralTHDMZ2 * myGTHDMZ2
 

Private Member Functions

void computeWFR_Z2 ()
 
void runGeneralTHDMZ2parameters ()
 

Private Attributes

double Ale
 
double alpha
 
double Als
 
double beta
 
double bma
 
double cosb
 
double cW2
 
double g1_at_Q
 
double g2_at_Q
 
double g3_at_Q
 
double lambda1_at_Q
 
double lambda2_at_Q
 
double lambda3_at_Q
 
double lambda4_at_Q
 
double lambda5_at_Q
 
double m11_2_at_Q
 
double m12_2
 
double m12_2_at_Q
 
double m22_2_at_Q
 
double mA2
 
double mHh2
 
double mHl2
 
double mHp2
 
gslpp::matrix< double > myZ2_at_Q
 
double MZ
 
double MZ2
 
double NLOuniscale
 
double Q_cutoff
 
double Q_GTHDM
 
double Rpeps
 
double sinb
 
double tanb
 
double vev
 
double WFRcomb1
 
double WFRcomb2
 
double WFRcomb3
 
double WFRcomb4
 
double Ybottom1_at_Q
 
double Ybottom2_at_Q
 
double Ytau1_at_Q
 
double Ytau2_at_Q
 
double Ytop_at_Q
 

Constructor & Destructor Documentation

◆ GeneralTHDMZ2Runner()

GeneralTHDMZ2Runner::GeneralTHDMZ2Runner ( const StandardModel &  SM_i)

GeneralTHDMZ2Runner constructor.

Definition at line 14 of file GeneralTHDMZ2Runner.cpp.

15: myGTHDMZ2(static_cast<const GeneralTHDMZ2*> (&SM_i)), myZ2_at_Q(3, 5, 0.)
16{}
GeneralTHDMZ2Runner::myZ2_at_Q
gslpp::matrix< double > myZ2_at_Q
Definition: GeneralTHDMZ2Runner.h:107
GeneralTHDMZ2Runner::myGTHDMZ2
const GeneralTHDMZ2 * myGTHDMZ2
Definition: GeneralTHDMZ2Runner.h:99

◆ ~GeneralTHDMZ2Runner()

GeneralTHDMZ2Runner::~GeneralTHDMZ2Runner ( )
virtual

Runner destructor.

Definition at line 18 of file GeneralTHDMZ2Runner.cpp.

19{};

Member Function Documentation

◆ betalambda1_Z2()

double GeneralTHDMZ2Runner::betalambda1_Z2 ( double  la1,
double  la3,
double  la4,
double  la5,
double  Yb1,
double  Ytau1 
)
inline

The beta function of \(\lambda_1\) appearing in unitarity conditions.

Returns
\(\beta_{\lambda_1}\) as defined in Cacchio:2016qyh

Definition at line 41 of file GeneralTHDMZ2Runner.h.

41 {
42 return ((12.*la1*la1 + 4.*la3*la3 + 4.*la3*la4 + 2.*la4*la4 + 2.*la5*la5 +
43 12.*la1*Yb1*Yb1 + 4.*la1*Ytau1*Ytau1 - 12.*Yb1*Yb1*Yb1*Yb1 -
44 4.*Ytau1*Ytau1*Ytau1*Ytau1)/16./M_PI/M_PI);
45 }

◆ betalambda2_Z2()

double GeneralTHDMZ2Runner::betalambda2_Z2 ( double  la2,
double  la3,
double  la4,
double  la5,
double  Yt,
double  Yb2,
double  Ytau2 
)
inline

The beta function of \(\lambda_2\) appearing in unitarity conditions.

Returns
\(\beta_{\lambda_2}\) as defined in Cacchio:2016qyh

Definition at line 52 of file GeneralTHDMZ2Runner.h.

52 {
53 return ((12.*la2*la2 + 4.*la3*la3 + 4.*la3*la4 + 2.*la4*la4 + 2.*la5*la5 +
54 12.*la2*Yb2*Yb2 + 4.*la2*Ytau2*Ytau2 + 12.*la2*Yt*Yt - 12.*Yb2*Yb2*Yb2*Yb2 -
55 4.*Ytau2*Ytau2*Ytau2*Ytau2 - 12.*Yt*Yt*Yt*Yt)/16./M_PI/M_PI);
56 }

◆ betalambda3_Z2()

double GeneralTHDMZ2Runner::betalambda3_Z2 ( double  la1,
double  la2,
double  la3,
double  la4,
double  la5,
double  Yt,
double  Yb1,
double  Yb2,
double  Ytau1,
double  Ytau2 
)
inline

The beta function of \(\lambda_3\) appearing in unitarity conditions.

Returns
\(\beta_{\lambda_3}\) as defined in Cacchio:2016qyh

Definition at line 63 of file GeneralTHDMZ2Runner.h.

63 {
64 return ((4.*la3*la3 + 2.*la4*la4 + (la1 + la2)*(6.*la3 + 2.*la4) + 2.*la5*la5 +
65 6.*la3*Yb1*Yb1 + 6.*la3*Yb2*Yb2 + 2.*la3*Ytau1*Ytau1 + 2.*la3*Ytau2*Ytau2 +
66 6.*la3*Yt*Yt - 12.*Yb1*Yb1*(Yb2*Yb2 + Yt*Yt) - 4.*Ytau1*Ytau1*Ytau2*Ytau2)/16./M_PI/M_PI);
67 }

◆ betalambda4_Z2()

double GeneralTHDMZ2Runner::betalambda4_Z2 ( double  la1,
double  la2,
double  la3,
double  la4,
double  la5,
double  Yt,
double  Yb1,
double  Yb2,
double  Ytau1,
double  Ytau2 
)
inline

The beta function of \(\lambda_4\) appearing in unitarity conditions.

Returns
\(\beta_{\lambda_4}\) as defined in Cacchio:2016qyh

Definition at line 74 of file GeneralTHDMZ2Runner.h.

74 {
75 return (((2.*la1 + 2.*la2 + 8.*la3)*la4 + 4.*la4*la4 + 8.*la5*la5 +
76 6.*la4*Yb1*Yb1 + 6.*la4*Yb2*Yb2 + 2.*la4*Ytau1*Ytau1 + 2.*la4*Ytau2*Ytau2 +
77 6.*la4*Yt*Yt - 4.*Ytau1*Ytau1*Ytau2*Ytau2 - 12.*Yb1*Yb1*(Yb2*Yb2 - Yt*Yt))/16./M_PI/M_PI);
78 }

◆ betalambda5_Z2()

double GeneralTHDMZ2Runner::betalambda5_Z2 ( double  la1,
double  la2,
double  la3,
double  la4,
double  la5,
double  Yt,
double  Yb1,
double  Yb2,
double  Ytau1,
double  Ytau2 
)
inline

The beta function of \(\lambda_5\) appearing in unitarity conditions.

Returns
\(\beta_{\lambda_5}\) as defined in Cacchio:2016qyh

Definition at line 85 of file GeneralTHDMZ2Runner.h.

85 {
86 return (((2.*la1 + 2.*la2 + 8.*la3 + 12.*la4)*la5 +
87 6.*la5*Yb1*Yb1 + 2.*la5*Ytau1*Ytau1 + 2.*la5*(3.*Yb2*Yb2 + Ytau2*Ytau2 + 3.*Yt*Yt) -
88 12.*Yb1*Yb1*Yb2*Yb2 - 4.*Ytau1*Ytau1*Ytau2*Ytau2)/16./M_PI/M_PI);
89 }

◆ computeWFR_Z2()

void GeneralTHDMZ2Runner::computeWFR_Z2 ( )
private

Definition at line 737 of file GeneralTHDMZ2Runner.cpp.

738{
739 double WFRcomb1a;
740 double WFRcomb1b;
741 double WFRcomb2a;
742 double WFRcomb3a;
743 double WFRcomb3b;
744 double WFRcomb4a;
745
746 double B000mh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_0_mHh2(MZ2,mHl2).real();
747 double B000mH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_0_mHh2(MZ2,mHh2).real();
748 double B00mHpmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_mHp2_mHl2(MZ2,mHp2,mHl2).real();
749 double B00mHpmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_mHp2_mHh2(MZ2,mHp2,mHh2).real();
750 double B00mAmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_mA2_mHl2(MZ2,mA2,mHl2).real();
751 double B00mAmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_0_mA2_mHh2(MZ2,mA2,mHh2).real();
752 double B0mh00 = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_0_0(MZ2,mHl2).real();
753 double B0mh0mHp = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_0_mHp2(MZ2,mHl2,mHp2).real();
754 double B0mh0mA = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_0_mA2(MZ2,mHl2,mA2).real();
755 double B0mhmhmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_mHl2_mHl2(MZ2,mHl2).real();
756 double B0mhmHmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_mHh2_mHl2(MZ2,mHl2,mHh2).real();
757 double B0mhmHmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_mHh2_mHh2(MZ2,mHl2,mHh2).real();
758 double B0mhmHpmHp = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_mHp2_mHp2(MZ2,mHl2,mHp2).real();
759 double B0mhmAmA = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHl2_mA2_mA2(MZ2,mHl2,mA2).real();
760 double B0mH00 = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_0_0(MZ2,mHh2).real();
761 double B0mH0mHp = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_0_mHp2(MZ2,mHh2,mHp2).real();
762 double B0mH0mA = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_0_mA2(MZ2,mHh2,mA2).real();
763 double B0mHmhmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_mHl2_mHl2(MZ2,mHh2,mHl2).real();
764 double B0mHmHmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_mHh2_mHl2(MZ2,mHh2,mHl2).real();
765 double B0mHmHmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_mHh2_mHh2(MZ2,mHh2).real();
766 double B0mHmHpmHp = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_mHp2_mHp2(MZ2,mHh2,mHp2).real();
767 double B0mHmAmA = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHh2_mA2_mA2(MZ2,mHh2,mA2).real();
768 double B0mHp0mh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHp2_0_mHl2(MZ2,mHp2,mHl2).real();
769 double B0mHp0mH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHp2_0_mHh2(MZ2,mHp2,mHh2).real();
770 double B0mHpmHpmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHp2_mHp2_mHl2(MZ2,mHp2,mHl2).real();
771 double B0mHpmHpmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mHp2_mHp2_mHh2(MZ2,mHp2,mHh2).real();
772 double B0mA0mh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mA2_0_mHl2(MZ2,mA2,mHl2).real();
773 double B0mA0mH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mA2_0_mHh2(MZ2,mA2,mHh2).real();
774 double B0mAmAmh = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mA2_mA2_mHl2(MZ2,mA2,mHl2).real();
775 double B0mAmAmH = myGTHDMZ2->getMyGTHDMCache()->B0_MZ2_mA2_mA2_mHh2(MZ2,mA2,mHh2).real();
776
777 double ddpB000mh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_0_mHl2(MZ2,mHl2).real();
778 double ddpB000mH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_0_mHh2(MZ2,mHh2).real();
779 double ddpB00mHpmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_mHp2_mHl2(MZ2,mHp2,mHl2).real();
780 double ddpB00mHpmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_mHp2_mHh2(MZ2,mHp2,mHh2).real();
781 double ddpB00mHpmA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_mHp2_mA2(MZ2,mHp2,mA2).real();
782 double ddpB00mAmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_mA2_mHl2(MZ2,mA2,mHl2).real();
783 double ddpB00mAmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_0_mA2_mHh2(MZ2,mA2,mHh2).real();
784 double ddpB0mh00 = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_0_0(MZ2,mHl2).real();
785 double ddpB0mh0mHp = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_0_mHp2(MZ2,mHl2,mHp2).real();
786 double ddpB0mh0mA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_0_mA2(MZ2,mHl2,mA2).real();
787 double ddpB0mhmhmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_mHl2_mHl2(MZ2,mHl2).real();
788 double ddpB0mhmHmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_mHh2_mHl2(MZ2,mHl2,mHh2).real();
789 double ddpB0mhmHmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_mHh2_mHh2(MZ2,mHl2,mHh2).real();
790 double ddpB0mhmHpmHp = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_mHp2_mHp2(MZ2,mHl2,mHp2).real();
791 double ddpB0mhmAmA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHl2_mA2_mA2(MZ2,mHl2,mA2).real();
792 double ddpB0mH00 = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_0_0(MZ2,mHh2).real();
793 double ddpB0mH0mHp = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_0_mHp2(MZ2,mHh2,mHp2).real();
794 double ddpB0mH0mA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_0_mA2(MZ2,mHh2,mA2).real();
795 double ddpB0mHmhmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_mHl2_mHl2(MZ2,mHh2,mHl2).real();
796 double ddpB0mHmHmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_mHh2_mHl2(MZ2,mHh2,mHl2).real();
797 double ddpB0mHmHmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_mHh2_mHh2(MZ2,mHh2).real();
798 double ddpB0mHmHpmHp = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_mHp2_mHp2(MZ2,mHh2,mHp2).real();
799 double ddpB0mHmAmA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHh2_mA2_mA2(MZ2,mHh2,mA2).real();
800 double ddpB0mHp0mh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHp2_0_mHl2(MZ2,mHp2,mHl2).real();
801 double ddpB0mHp0mH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHp2_0_mHh2(MZ2,mHp2,mHh2).real();
802 double ddpB0mHp0mA = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHp2_0_mA2(MZ2,mHp2,mA2).real();
803 double ddpB0mHpmHpmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHp2_mHp2_mHl2(MZ2,mHp2,mHl2).real();
804 double ddpB0mHpmHpmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mHp2_mHp2_mHh2(MZ2,mHp2,mHh2).real();
805 double ddpB0mA0mh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mA2_0_mHl2(MZ2,mA2,mHl2).real();
806 double ddpB0mA0mH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mA2_0_mHh2(MZ2,mA2,mHh2).real();
807 double ddpB0mA0mHp = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mA2_0_mHp2(MZ2,mA2,mHp2).real();
808 double ddpB0mAmAmh = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mA2_mA2_mHl2(MZ2,mA2,mHl2).real();
809 double ddpB0mAmAmH = myGTHDMZ2->getMyGTHDMCache()->B0p_MZ2_mA2_mA2_mHh2(MZ2,mA2,mHh2).real();
810
811 WFRcomb1a = 3.0*mHl2*mHl2*cosb*cosb*sin(bma)*sin(bma) * ddpB000mh
812 + 3.0*mHh2*mHh2*cos(bma)*cos(bma)*cosb*cosb * ddpB000mH
813 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*cosb*cosb * ddpB00mHpmh
814 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma)*cosb*cosb * ddpB00mHpmH
815 + 2.0*(mA2-mHp2)*(mA2-mHp2)*cosb*cosb * ddpB00mHpmA
816 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*cosb*cosb * ddpB00mAmh
817 + (mA2-mHh2)*(mA2-mHh2)*cosb*cosb*sin(bma)*sin(bma) * ddpB00mAmH
818 + 1.5*mHl2*mHl2*sin(alpha)*sin(alpha)*sin(bma)*sin(bma) * ddpB0mh00
819 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*sin(alpha)*sin(alpha) * ddpB0mh0mHp
820 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*sin(alpha)*sin(alpha) * ddpB0mh0mA
821 + 9.0*sin(alpha)*sin(alpha)*pow(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
822 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
823 + sin(alpha)*sin(alpha)*pow((cos(alpha)/sinb+sin(alpha)/cosb)*(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-(3.0*m12_2)/(cosb*sinb))),2) * ddpB0mhmHmh
824 + sin(alpha)*sin(alpha)*sin(bma)*sin(bma)*pow(-2.0*m12_2+(2.0*mHh2+mHl2-(3.0*m12_2)/(cosb*sinb))*sin(2.0*alpha),2)/(8.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
825 + (sin(alpha)*sin(alpha)*pow((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
826 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb),2))/(64.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
827 + (sin(alpha)*sin(alpha)*pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb),2))/(128.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
828 + 1.5*mHh2*mHh2*cos(alpha)*cos(alpha)*cos(bma)*cos(bma) * ddpB0mH00
829 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*cos(alpha)*cos(alpha)*sin(bma)*sin(bma) * ddpB0mH0mHp
830 + (mA2-mHh2)*(mA2-mHh2)*cos(alpha)*cos(alpha)*sin(bma)*sin(bma) * ddpB0mH0mA
831 + cos(alpha)*cos(alpha)*cos(bma)*cos(bma)*pow(m12_2+cos(alpha)*sin(alpha)*(mHh2-3.0*m12_2/(cosb*sinb))+mHl2*sin(2.0*alpha),2)/(2.0*cosb*cosb*sinb*sinb) * ddpB0mHmhmh
832 + cos(alpha)*cos(alpha)*sin(bma)*sin(bma)*pow(m12_2*cosb*sinb+0.5*sin(2.0*alpha)*(3.0*m12_2-(2.0*mHh2+mHl2)*cosb*sinb),2)/pow(cosb*sinb,4) * ddpB0mHmHmh
833 + 9.0*cos(alpha)*cos(alpha)*pow(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
834 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
835 + (cos(alpha)*cos(alpha)*pow((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)
836 +2.0*(mHh2+2.0*mHp2)*cos(bma)-(3.0*mHh2+2.0*mHp2)*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta),2))/(256.0*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
837 + (cos(alpha)*cos(alpha)*pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)
838 -2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta),2))/(512.0*pow(cosb*sinb,4)) * ddpB0mHmAmA
839 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*sinb*sinb * ddpB0mHp0mh
840 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma)*sinb*sinb * ddpB0mHp0mH
841 + 2.0*(mA2-mHp2)*(mA2-mHp2)*sinb*sinb * ddpB0mHp0mA
842 + 2.0*pow((m12_2*cos(alpha+beta))/(sinb*sinb*cosb*cosb)-(mHl2*cos(bma)*cos(2.0*beta))/(cosb*sinb)-(mHl2+2.0*mHp2)*sin(bma),2)*sinb*sinb * ddpB0mHpmHpmh
843 + 2.0*pow(sinb*((mHh2+2.0*mHp2)*cos(bma)-mHh2*cos(2.0*beta)*sin(bma)/(cosb*sinb)-m12_2*sin(alpha+beta)/(sinb*sinb*cosb*cosb)),2) * ddpB0mHpmHpmH
844 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*sinb*sinb * ddpB0mA0mh
845 + (mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma)*sinb*sinb * ddpB0mA0mH
846 + 2.0*(mA2-mHp2)*(mA2-mHp2)*sinb*sinb * ddpB0mA0mHp
847 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb),2)/(64.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
848 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
849 + 3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta),2)/(256.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
850
851 WFRcomb1b = (mHl2*(mA2*(2.0*mHl2-3.0*mHp2)+mHl2*mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/(2.0*mA2*mHp2) * B000mh
852 - (mHh2*(mA2*(2.0*mHh2-3.0*mHp2)+mHh2*mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/(2.0*mA2*mHp2) * B000mH
853 + ((mHl2-mHp2)*cos(bma)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
854 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb)))/(2.0*mHp2*sinb*cosb) * B00mHpmh
855 + ((mHp2-mHh2)*sin(bma)*((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
856 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta)))/(4.0*mHp2*sinb*cosb) * B00mHpmH
857 + ((mHl2-mA2)*cos(bma)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
858 +cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb)))/(4.0*mA2*sinb*cosb) * B00mAmh
859 + ((mHh2-mA2)*sin(bma)*((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
860 +3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta)))/(8.0*mA2*sinb*cosb) * B00mAmH
861 + (3.0*mHh2*mHl2*sin(2.0*alpha)*sin(2.0*bma))/(4.0*(mHh2-mHl2)) * B0mh00
862 + ((mHl2-mHp2)*(mHp2-mHh2)*sin(2.0*alpha)*sin(2.0*bma))/(mHh2-mHl2) * B0mh0mHp
863 - ((mA2-mHh2)*(mA2-mHl2)*sin(2.0*alpha)*sin(2.0*bma))/(2.0*(mHh2-mHl2)) * B0mh0mA
864 + 3.0*cos(bma)*sin(2.0*alpha)*(m12_2+cos(alpha)*sin(alpha)*(mHh2-(3.0*m12_2)/(sinb*cosb))+mHl2*sin(2.0*alpha))
865 *(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
866 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta))/(32.0*(mHl2-mHh2)*pow(sinb*cosb,3)) * B0mhmhmh
867 + (sin(2.0*bma)*sin(2.0*alpha)
868 *(4.0*m12_2*m12_2+2.0*m12_2*(mHl2-mHh2)*sin(2.0*alpha)
869 -(2.0*mHh2*mHh2+5.0*mHh2*mHl2+2.0*mHl2*mHl2-(9.0*m12_2*(mHh2+mHl2))/(sinb*cosb)
870 +(9.0*m12_2*m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))/(8.0*(mHh2-mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
871 - 3.0*sin(2.0*alpha)*sin(bma)*(m12_2-cos(alpha)*sin(alpha)*(2.0*mHh2+mHl2-(3.0*m12_2)/(sinb*cosb)))
872 *(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
873 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mhmHmH
874 + sin(2.0*alpha)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb))
875 *((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
876 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(128.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mhmHpmHp
877 + sin(2.0*alpha)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb))
878 *((mHh2-2.0*mA2)*cos(alpha-5.0*beta)+2.0*(2.0*mA2+mHh2)*cos(bma)-2.0*mA2*cos(alpha+3.0*beta)
879 -3.0*mHh2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(256.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mhmAmA
880 - (3.0*mHh2*mHl2*sin(2.0*alpha)*sin(2.0*bma))/(4.0*(mHh2-mHl2)) * B0mH00
881 + ((mHh2-mHp2)*(mHl2-mHp2)*sin(2.0*alpha)*sin(2.0*bma))/(mHh2-mHl2) * B0mH0mHp
882 + ((mA2-mHh2)*(mA2-mHl2)*sin(2.0*alpha)*sin(2.0*bma))/(2.0*(mHh2-mHl2)) * B0mH0mA
883 + 3.0*cos(bma)*sin(2.0*alpha)*(m12_2+cos(alpha)*sin(alpha)*(mHh2-(3.0*m12_2)/(sinb*cosb))+mHl2*sin(2.0*alpha))
884 *(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
885 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mHmhmh
886 + (sin(2.0*bma)*sin(2.0*alpha)*(-4.0*m12_2*m12_2+2.0*m12_2*(mHh2-mHl2)*sin(2.0*alpha)
887 +(2.0*mHh2*mHh2+5.0*mHh2*mHl2+2.0*mHl2*mHl2
888 -(9.0*m12_2*(mHh2+mHl2))/(sinb*cosb)
889 +(9.0*m12_2*m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))
890 /(8.0*(mHh2-mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
891 + 3.0*sin(bma)*sin(2.0*alpha)*(m12_2-cos(alpha)*sin(alpha)*(2.0*mHh2+mHl2-(3.0*m12_2)/(sinb*cosb)))
892 *(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
893 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mHmHmH
894 - sin(2.0*alpha)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb))
895 *((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
896 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(128.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mHmHpmHp
897 - sin(2.0*alpha)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mA2)*2.0*sinb*cosb))
898 *((mHh2-2.0*mA2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mA2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
899 -2.0*mA2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(256.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mHmAmA
900 - (mHl2*(mHl2-mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/mHp2 * B0mHp0mh
901 + (mHh2*(mHh2-mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/mHp2 * B0mHp0mH
902 + ((mHp2-mHl2)*cos(bma)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
903 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb)))/(2.0*mHp2*sinb*cosb) * B0mHpmHpmh
904 + ((mHh2-mHp2)*sin(bma)*((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
905 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta)))/(4.0*mHp2*sinb*cosb) * B0mHpmHpmH
906 + (mHl2*(mA2-mHl2)*sin(bma)*cos(bma)*2.0*sinb*cosb)/mA2 * B0mA0mh
907 + ((mHh2-mA2)*mHh2*sin(bma)*cos(bma)*2.0*sinb*cosb)/mA2 * B0mA0mH
908 + ((mA2-mHl2)*cos(bma)*((-2.0*mA2+mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
909 +cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb)))/(4.0*mA2*sinb*cosb) * B0mAmAmh
910 +((mA2-mHh2)*sin(bma)*((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
911 +3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta)))/(8.0*mA2*sinb*cosb) * B0mAmAmH;
912
913 WFRcomb2a = 1.5*mHl2*mHl2*sin(bma)*sin(bma) * ddpB000mh
914 + 1.5*mHh2*mHh2*cos(bma)*cos(bma) * ddpB000mH
915 + (mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma) * ddpB00mHpmh
916 + (mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma) * ddpB00mHpmH
917 + (mA2-mHp2)*(mA2-mHp2) * ddpB00mHpmA
918 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB00mAmh
919 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB00mAmH
920 + 0.75*mHl2*mHl2*sin(bma)*sin(bma) * ddpB0mh00
921 + (mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma) * ddpB0mh0mHp
922 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB0mh0mA
923 + 9.0*pow(16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta)
924 -mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta)),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
925 + 0.5*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
926 *pow(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-3.0*m12_2/(cosb*sinb)),2) * ddpB0mhmHmh
927 + (pow(-2.0*m12_2+(2.0*mHh2+mHl2-3.0*m12_2/(cosb*sinb))*sin(2.0*alpha),2)
928 *sin(bma)*sin(bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
929 + pow((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*cosb*sinb
930 +cos(alpha+beta)*(-4.0*m12_2+(3.0*mHl2+2.0*mHp2)*cosb*sinb),2)/(32.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
931 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*cosb*sinb
932 +cos(alpha+beta)*(4.0*m12_2-(2.0*mA2+3.0*mHl2)*cosb*sinb),2)/(64.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
933 + 0.75*mHh2*mHh2*cos(bma)*cos(bma) * ddpB0mH00
934 + (mHh2 - mHp2)*(mHh2 - mHp2)*sin(bma)*sin(bma) * ddpB0mH0mHp
935 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB0mH0mA
936 + 0.25*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
937 *pow(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-3.0*m12_2/(cosb*sinb)),2) * ddpB0mHmhmh
938 + (pow(-2.0*m12_2+(2.0*mHh2+mHl2-3.0*m12_2/(cosb*sinb))*sin(2.0*alpha),2)
939 *sin(bma)*sin(bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
940 + 9.0*pow(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
941 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
942 + pow((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-(3.0*mHh2+2.0*mHp2)*cos(alpha+3.0*beta)
943 -16.0*m12_2*sin(alpha+beta),2)/(512*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
944 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)
945 +16.0*m12_2*sin(alpha+beta),2)/(1024*pow(cosb*sinb,4)) * ddpB0mHmAmA
946 + (mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma) * ddpB0mHp0mh
947 + (mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma) * ddpB0mHp0mH
948 + (mA2-mHp2)*(mA2-mHp2) * ddpB0mHp0mA
949 + pow((m12_2*cos(alpha+beta))/(cosb*cosb*sinb*sinb)-(mHl2*cos(bma)*cos(2.0*beta))/(cosb*sinb)
950 -(mHl2+2.0*mHp2)*sin(bma),2) * ddpB0mHpmHpmh
951 + pow((mHh2+2.0*mHp2)*cos(bma)-(mHh2*cos(2.0*beta)*sin(bma))/(cosb*sinb)
952 -(m12_2*sin(alpha+beta))/(cosb*cosb*sinb*sinb),2) * ddpB0mHpmHpmH
953 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB0mA0mh
954 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB0mA0mH
955 + (mA2-mHp2)*(mA2-mHp2) * ddpB0mA0mHp
956 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*cosb*sinb
957 +cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*cosb*sinb),2)/(128.0*pow(cosb*sinb,4)) * ddpB0mAmAmh
958 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)
959 +16.0*m12_2*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mAmAmH;
960
961 WFRcomb3a = 0.5*mHl2*mHl2*(3.0-cos(2.0*beta))*sin(bma)*sin(bma) * ddpB000mh
962 + 0.5*mHh2*mHh2*(3.0-cos(2.0*beta))*cos(bma)*cos(bma) * ddpB000mH
963 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*sinb*sinb * ddpB00mHpmh
964 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma)*sinb*sinb * ddpB00mHpmH
965 + 2.0*(mA2-mHp2)*(mA2-mHp2)*sinb*sinb * ddpB00mHpmA
966 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*cosb*cosb * ddpB00mAmh
967 + (mA2-mHh2)*(mA2-mHh2)*cosb*cosb*sin(bma)*sin(bma) * ddpB00mAmH
968 + 1.5*mHl2*mHl2*sin(alpha)*sin(alpha)*sin(bma)*sin(bma) * ddpB0mh00
969 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*sin(alpha)*sin(alpha) * ddpB0mh0mHp
970 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*sin(alpha)*sin(alpha) * ddpB0mh0mA
971 + 9.0*sin(alpha)*sin(alpha)*pow(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
972 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
973 + sin(alpha)*sin(alpha)*pow((cos(alpha)/sinb+sin(alpha)/cosb)*(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-(3.0*m12_2)/(cosb*sinb))),2) * ddpB0mhmHmh
974 + sin(alpha)*sin(alpha)*sin(bma)*sin(bma)*pow(-2.0*m12_2+(2.0*mHh2+mHl2-(3.0*m12_2)/(cosb*sinb))*sin(2.0*alpha),2)/(8.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
975 + (sin(alpha)*sin(alpha)*pow((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
976 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb),2))/(64.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
977 + (sin(alpha)*sin(alpha)*pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb),2))/(128.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
978 + 1.5*mHh2*mHh2*cos(alpha)*cos(alpha)*cos(bma)*cos(bma) * ddpB0mH00
979 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*cos(alpha)*cos(alpha)*sin(bma)*sin(bma) * ddpB0mH0mHp
980 + (mA2-mHh2)*(mA2-mHh2)*cos(alpha)*cos(alpha)*sin(bma)*sin(bma) * ddpB0mH0mA
981 + cos(alpha)*cos(alpha)*cos(bma)*cos(bma)*pow(m12_2+cos(alpha)*sin(alpha)*(mHh2-3.0*m12_2/(cosb*sinb))+mHl2*sin(2.0*alpha),2)/(2.0*cosb*cosb*sinb*sinb) * ddpB0mHmhmh
982 + cos(alpha)*cos(alpha)*sin(bma)*sin(bma)*pow(m12_2*cosb*sinb+0.5*sin(2.0*alpha)*(3.0*m12_2-(2.0*mHh2+mHl2)*cosb*sinb),2)/pow(cosb*sinb,4) * ddpB0mHmHmh
983 + 9.0*cos(alpha)*cos(alpha)*pow(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
984 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
985 + (cos(alpha)*cos(alpha)*pow((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)
986 +2.0*(mHh2+2.0*mHp2)*cos(bma)-(3.0*mHh2+2.0*mHp2)*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta),2))/(256.0*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
987 + (cos(alpha)*cos(alpha)*pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)
988 -2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta),2))/(512.0*pow(cosb*sinb,4)) * ddpB0mHmAmA
989 + 2.0*(mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma)*cosb*cosb * ddpB0mHp0mh
990 + 2.0*(mHh2-mHp2)*(mHh2-mHp2)*sin(bma)*sin(bma)*cosb*cosb * ddpB0mHp0mH
991 + 2.0*(mA2-mHp2)*(mA2-mHp2)*cosb*cosb * ddpB0mHp0mA
992 + 2.0*cosb*cosb*pow((m12_2*cos(alpha+beta))/(cosb*cosb*sinb*sinb)-(mHl2*cos(bma)*cos(2.0*beta))/(cosb*sinb)-(mHl2+2.0*mHp2)*sin(bma),2) * ddpB0mHpmHpmh
993 + 2.0*cosb*cosb*pow((mHh2+2.0*mHp2)*cos(bma)-(mHh2*cos(2.0*beta)*sin(bma))/(cosb*sinb)-(m12_2*sin(alpha+beta))/(cosb*cosb*sinb*sinb),2) * ddpB0mHpmHpmH
994 + (mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma)*sinb*sinb * ddpB0mA0mh
995 + (mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma)*sinb*sinb * ddpB0mA0mH
996 + 2.0*(mA2-mHp2)*(mA2-mHp2)*sinb*sinb * ddpB0mA0mHp
997 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb),2)/(64.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
998 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
999 + 3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta),2)/(256.0*pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
1000
1001 WFRcomb3b = ((mHl2*mHl2*mHp2+mA2*(-2.0*mHl2*mHl2+mHl2*mHp2))*sin(2.0*bma)*sinb*cosb)/(mA2*mHp2) * B000mh
1002 + ((-mHh2*mHh2*mHp2+mA2*(2.0*mHh2*mHh2-mHh2*mHp2))*sin(2.0*bma)*sinb*cosb)/(mA2*mHp2) * B000mH
1003 + ((mHp2-mHl2)*cos(bma)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
1004 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb)))/(2.0*mHp2*sinb*cosb) * B00mHpmh
1005 + ((mHh2-mHp2)*sin(bma)*((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
1006 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta)))/(4.0*mHp2*sinb*cosb) * B00mHpmH
1007 + ((mHl2-mA2)*cos(bma)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
1008 +cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb)))/(4.0*mA2*sinb*cosb) * B00mAmh
1009 + ((mHh2-mA2)*sin(bma)*((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
1010 +3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta)))/(8.0*mA2*sinb*cosb) * B00mAmH
1011 + (3.0*mHh2*mHl2*sin(2.0*alpha)*sin(2.0*bma))/(4.0*(mHh2-mHl2)) * B0mh00
1012 + ((mHl2-mHp2)*(mHp2-mHh2)*sin(2.0*alpha)*sin(2.0*bma))/(mHh2-mHl2) * B0mh0mHp
1013 - ((mA2-mHh2)*(mA2-mHl2)*sin(2.0*alpha)*sin(2.0*bma))/(2.0*(mHh2-mHl2)) * B0mh0mA
1014 + 3.0*cos(bma)*sin(2.0*alpha)*(m12_2+cos(alpha)*sin(alpha)*(mHh2-(3.0*m12_2)/(sinb*cosb))+mHl2*sin(2.0*alpha))
1015 *(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
1016 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta))/(32.0*(mHl2-mHh2)*pow(sinb*cosb,3)) * B0mhmhmh
1017 + (sin(2.0*bma)*sin(2.0*alpha)
1018 *(4.0*m12_2*m12_2+2.0*m12_2*(mHl2-mHh2)*sin(2.0*alpha)
1019 -(2.0*mHh2*mHh2+5.0*mHh2*mHl2+2.0*mHl2*mHl2-(9.0*m12_2*(mHh2+mHl2))/(sinb*cosb)
1020 +(9.0*m12_2*m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))/(8.0*(mHh2-mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
1021 - 3.0*sin(2.0*alpha)*sin(bma)*(m12_2-cos(alpha)*sin(alpha)*(2.0*mHh2+mHl2-(3.0*m12_2)/(sinb*cosb)))
1022 *(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
1023 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mhmHmH
1024 + sin(2.0*alpha)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb))
1025 *((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
1026 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(128.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mhmHpmHp
1027 + sin(2.0*alpha)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb))
1028 *((mHh2-2.0*mA2)*cos(alpha-5.0*beta)+2.0*(2.0*mA2+mHh2)*cos(bma)-2.0*mA2*cos(alpha+3.0*beta)
1029 -3.0*mHh2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(256.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mhmAmA
1030 - (3.0*mHh2*mHl2*sin(2.0*alpha)*sin(2.0*bma))/(4.0*(mHh2-mHl2)) * B0mH00
1031 + ((mHh2-mHp2)*(mHl2-mHp2)*sin(2.0*alpha)*sin(2.0*bma))/(mHh2-mHl2) * B0mH0mHp
1032 + ((mA2-mHh2)*(mA2-mHl2)*sin(2.0*alpha)*sin(2.0*bma))/(2.0*(mHh2-mHl2)) * B0mH0mA
1033 + 3.0*cos(bma)*sin(2.0*alpha)*(m12_2+cos(alpha)*sin(alpha)*(mHh2-(3.0*m12_2)/(sinb*cosb))+mHl2*sin(2.0*alpha))
1034 *(-mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta))
1035 +16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mHmhmh
1036 + (sin(2.0*bma)*sin(2.0*alpha)*(-4.0*m12_2*m12_2+2.0*m12_2*(mHh2-mHl2)*sin(2.0*alpha)
1037 +(2.0*mHh2*mHh2+5.0*mHh2*mHl2+2.0*mHl2*mHl2
1038 -(9.0*m12_2*(mHh2+mHl2))/(sinb*cosb)
1039 +(9.0*m12_2*m12_2)/(sinb*sinb*cosb*cosb))*sin(2.0*alpha)*sin(2.0*alpha)))
1040 /(8.0*(mHh2-mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
1041 + 3.0*sin(bma)*sin(2.0*alpha)*(m12_2-cos(alpha)*sin(alpha)*(2.0*mHh2+mHl2-(3.0*m12_2)/(sinb*cosb)))
1042 *(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
1043 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta))/(32.0*(mHh2-mHl2)*pow(sinb*cosb,3)) * B0mHmHmH
1044 - sin(2.0*alpha)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb))
1045 *((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
1046 -2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(128.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mHmHpmHp
1047 - sin(2.0*alpha)*((mHl2-2.0*mA2)*cos(alpha-3.0*beta)*2.0*sinb*cosb+cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mA2)*2.0*sinb*cosb))
1048 *((mHh2-2.0*mA2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mA2)*cos(bma)-3.0*mHh2*cos(alpha+3.0*beta)
1049 -2.0*mA2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(256.0*(mHh2-mHl2)*pow(cosb*sinb,4)) * B0mHmAmA
1050 + (mHl2*(mHl2-mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/mHp2 * B0mHp0mh
1051 - (mHh2*(mHh2-mHp2)*sin(2.0*bma)*2.0*sinb*cosb)/mHp2 * B0mHp0mH
1052 + (mHl2-mHp2)*cos(bma)*((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
1053 +cos(alpha+beta)*(-8.0*m12_2+(3.0*mHl2+2.0*mHp2)*2.0*sinb*cosb))/(2.0*mHp2*sinb*cosb) * B0mHpmHpmh
1054 - (mHh2-mHp2)*sin(bma)*((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)
1055 -3.0*mHh2*cos(alpha+3.0*beta)-2.0*mHp2*cos(alpha+3.0*beta)-16.0*m12_2*sin(alpha+beta))/(4.0*mHp2*sinb*cosb) * B0mHpmHpmH
1056 + (mHl2*(mA2-mHl2)*sin(bma)*cos(bma)*2.0*sinb*cosb)/mA2 * B0mA0mh
1057 + ((mHh2-mA2)*mHh2*sin(bma)*cos(bma)*2.0*sinb*cosb)/mA2 * B0mA0mH
1058 + ((mA2-mHl2)*cos(bma)*((-2.0*mA2+mHl2)*cos(alpha-3.0*beta)*2.0*sinb*cosb
1059 +cos(alpha+beta)*(-8.0*m12_2+(2.0*mA2+3.0*mHl2)*2.0*sinb*cosb)))/(4.0*mA2*sinb*cosb) * B0mAmAmh
1060 +((mA2-mHh2)*sin(bma)*((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+2.0*mA2*cos(alpha+3.0*beta)
1061 +3.0*mHh2*cos(alpha+3.0*beta)+16.0*m12_2*sin(alpha+beta)))/(8.0*mA2*sinb*cosb) * B0mAmAmH;
1062
1063 WFRcomb4a = 0.5*mHl2*mHl2*sin(bma)*sin(bma) * ddpB000mh
1064 + 0.5*mHh2*mHh2*cos(bma)*cos(bma) * ddpB000mH
1065 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB00mAmh
1066 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB00mAmH
1067 + 0.75*mHl2*mHl2*sin(bma)*sin(bma) * ddpB0mh00
1068 + (mHl2-mHp2)*(mHl2-mHp2)*cos(bma)*cos(bma) * ddpB0mh0mHp
1069 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB0mh0mA
1070 + 9.0*pow(16.0*m12_2*cos(bma)*cos(bma)*cos(alpha+beta)
1071 -mHl2*(3.0*sin(bma)+sin(3.0*bma)+sin(3.0*alpha+beta)+3.0*sin(alpha+3.0*beta)),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mhmhmh
1072 + 0.5*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
1073 *pow(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-3.0*m12_2/(cosb*sinb)),2) * ddpB0mhmHmh
1074 + (pow(-2.0*m12_2+(2.0*mHh2+mHl2-3.0*m12_2/(cosb*sinb))*sin(2.0*alpha),2)
1075 *sin(bma)*sin(bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
1076 + pow((mHl2-2.0*mHp2)*cos(alpha-3.0*beta)*cosb*sinb
1077 +cos(alpha+beta)*(-4.0*m12_2+(3.0*mHl2+2.0*mHp2)*cosb*sinb),2)/(32.0*pow(cosb*sinb,4)) * ddpB0mhmHpmHp
1078 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*cosb*sinb
1079 +cos(alpha+beta)*(4.0*m12_2-(2.0*mA2+3.0*mHl2)*cosb*sinb),2)/(64.0*pow(cosb*sinb,4)) * ddpB0mhmAmA
1080 + 0.75*mHh2*mHh2*cos(bma)*cos(bma) * ddpB0mH00
1081 + (mHh2 - mHp2)*(mHh2 - mHp2)*sin(bma)*sin(bma) * ddpB0mH0mHp
1082 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB0mH0mA
1083 + 0.25*pow(cos(alpha)/sinb + sin(alpha)/cosb,2)
1084 *pow(m12_2+cos(alpha)*sin(alpha)*(mHh2+2.0*mHl2-3.0*m12_2/(cosb*sinb)),2) * ddpB0mHmhmh
1085 + (pow(-2.0*m12_2+(2.0*mHh2+mHl2-3.0*m12_2/(cosb*sinb))*sin(2.0*alpha),2)
1086 *sin(bma)*sin(bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
1087 + 9.0*pow(mHh2*(-3.0*cos(bma)+cos(3.0*bma)-cos(3.0*alpha+beta)+3.0*cos(alpha+3.0*beta))
1088 +16.0*m12_2*sin(bma)*sin(bma)*sin(alpha+beta),2)/(1024.0*pow(cosb*sinb,4)) * ddpB0mHmHmH
1089 + pow((mHh2-2.0*mHp2)*cos(alpha-5.0*beta)+2.0*(mHh2+2.0*mHp2)*cos(bma)-(3.0*mHh2+2.0*mHp2)*cos(alpha+3.0*beta)
1090 -16.0*m12_2*sin(alpha+beta),2)/(512*pow(cosb*sinb,4)) * ddpB0mHmHpmHp
1091 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)
1092 +16.0*m12_2*sin(alpha+beta),2)/(1024*pow(cosb*sinb,4)) * ddpB0mHmAmA
1093 + 0.5*(mA2-mHl2)*(mA2-mHl2)*cos(bma)*cos(bma) * ddpB0mA0mh
1094 + 0.5*(mA2-mHh2)*(mA2-mHh2)*sin(bma)*sin(bma) * ddpB0mA0mH
1095 + (mA2-mHp2)*(mA2-mHp2) * ddpB0mA0mHp
1096 + pow((2.0*mA2-mHl2)*cos(alpha-3.0*beta)*2.0*cosb*sinb
1097 +cos(alpha+beta)*(8.0*m12_2-(2.0*mA2+3.0*mHl2)*2.0*cosb*sinb),2)/(128.0*pow(cosb*sinb,4)) * ddpB0mAmAmh
1098 + pow((2.0*mA2-mHh2)*cos(alpha-5.0*beta)-2.0*(2.0*mA2+mHh2)*cos(bma)+(2.0*mA2+3.0*mHh2)*cos(alpha+3.0*beta)
1099 +16.0*m12_2*sin(alpha+beta),2)/(512.0*pow(cosb*sinb,4)) * ddpB0mAmAmH;
1100
1101 WFRcomb1 = -(WFRcomb1a+WFRcomb1b)/(vev*vev);
1102 WFRcomb2 = -WFRcomb2a/(vev*vev);
1103 WFRcomb3 = -(WFRcomb3a+WFRcomb3b)/(vev*vev);
1104 WFRcomb4 = -WFRcomb4a/(vev*vev);
1105}

◆ getGTHDMZ2_at_Q()

gslpp::matrix< double > GeneralTHDMZ2Runner::getGTHDMZ2_at_Q ( )

The public function which contains all relevant GTHDMZ2 parameter after running.

Returns
\(\big((\lambda_1,\lambda_2,\lambda_3,\lambda_4,\lambda_5),(Y_t,Y_{b,1},Y_{b,2},Y_{\tau,1},Y_{\tau,2}),(\text{WFR}_1,\text{WFR}_2,\text{WFR}_3,\text{WFR}_4,0)\big)\)

Definition at line 1108 of file GeneralTHDMZ2Runner.cpp.

1109{
1110 vev = myGTHDMZ2->v();
1111 cW2 = myGTHDMZ2->c02();
1112 Ale = myGTHDMZ2->getAle();
1113 Als = myGTHDMZ2->getAlsMz();
1114 MZ = myGTHDMZ2->getMz();
1115 MZ2 = MZ*MZ;
1116
1117 m12_2 = myGTHDMZ2->getm12sq_Z2();
1118 mHl2 = myGTHDMZ2->getmH1sq();
1119 mHh2 = myGTHDMZ2->getmH2sq();
1120 mA2 = myGTHDMZ2->getmH3sq();
1121 mHp2 = myGTHDMZ2->getmHp2();
1122
1123 beta = myGTHDMZ2->getbeta_Z2();
1124 tanb = myGTHDMZ2->gettanb_Z2();
1125 sinb = myGTHDMZ2->getsinb_Z2();
1126 cosb = myGTHDMZ2->getcosb_Z2();
1127 bma = myGTHDMZ2->getbma_Z2();
1128 alpha = beta-bma;
1129
1130 Q_GTHDM = myGTHDMZ2->getQ_GTHDM();
1131 Rpeps = myGTHDMZ2->getRpepsGTHDM();
1132 NLOuniscale = myGTHDMZ2->getNLOuniscaleGTHDM();
1133
1134 runGeneralTHDMZ2parameters();
1135
1136 // Wavefunction renormalization contributions to unitarity, from Grinstein:2015rtl
1137 WFRcomb1 = 0.;
1138 WFRcomb2 = 0.;
1139 WFRcomb3 = 0.;
1140 WFRcomb4 = 0.;
1141
1142 if(myGTHDMZ2->getWFRflag_Z2())
1143 computeWFR_Z2();
1144
1145 // 1st row: (lambda1, lambda2, lambda3, lambda4, lambda5)
1146 // 2nd row: (Ytop, Ybottom1, Ybottom2, Ytau1, Ytau2)
1147 // 3rd row: (WFRcomb1, WFRcomb2, WFRcomb3, WFRcomb4, 0.0)
1148 myZ2_at_Q.assign(0, 0, lambda1_at_Q);
1149 myZ2_at_Q.assign(0, 1, lambda2_at_Q);
1150 myZ2_at_Q.assign(0, 2, lambda3_at_Q);
1151 myZ2_at_Q.assign(0, 3, lambda4_at_Q);
1152 myZ2_at_Q.assign(0, 4, lambda5_at_Q);
1153 myZ2_at_Q.assign(1, 0, Ytop_at_Q);
1154 myZ2_at_Q.assign(1, 1, Ybottom1_at_Q);
1155 myZ2_at_Q.assign(1, 2, Ybottom2_at_Q);
1156 myZ2_at_Q.assign(1, 3, Ytau1_at_Q);
1157 myZ2_at_Q.assign(1, 4, Ytau2_at_Q);
1158 myZ2_at_Q.assign(2, 0, WFRcomb1);
1159 myZ2_at_Q.assign(2, 1, WFRcomb2);
1160 myZ2_at_Q.assign(2, 2, WFRcomb3);
1161 myZ2_at_Q.assign(2, 3, WFRcomb4);
1162
1163 return (myZ2_at_Q);
1164}
GeneralTHDMZ2::getsinb_Z2
double getsinb_Z2() const
A getter for the sine of .
Definition: GeneralTHDMZ2.h:109
GeneralTHDMZ2::getcosb_Z2
double getcosb_Z2() const
A getter for the cosine of .
Definition: GeneralTHDMZ2.h:118
GeneralTHDMZ2::gettanb_Z2
double gettanb_Z2() const
A getter for the tangent of .
Definition: GeneralTHDMZ2.h:100
GeneralTHDMZ2::getWFRflag_Z2
bool getWFRflag_Z2() const
A getter for the chosen option to include WFR in NLO unitarity.
Definition: GeneralTHDMZ2.h:199
GeneralTHDMZ2::getbeta_Z2
double getbeta_Z2() const
A getter for .
Definition: GeneralTHDMZ2.h:91
GeneralTHDMZ2::getbma_Z2
double getbma_Z2() const
A getter for the difference between and .
Definition: GeneralTHDMZ2.h:127
GeneralTHDMZ2::getm12sq_Z2
double getm12sq_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:136

◆ RGEGeneralTHDMZ2Runner()

double GeneralTHDMZ2Runner::RGEGeneralTHDMZ2Runner ( double  InitialValues[],
unsigned long int  NumberOfRGEs,
double  Q1,
double  Q2,
int  order 
)
virtual

Definition at line 532 of file GeneralTHDMZ2Runner.cpp.

533{
534 //Define which stepping function should be used
535 const gsl_odeiv2_step_type * T = gsl_odeiv2_step_rk4;
536
537 //Allocate space for the stepping function
538 gsl_odeiv2_step * s = gsl_odeiv2_step_alloc(T, NumberOfRGEs);
539
540 //Define the absolute (A) and relative (R) error on y at each step.
541 //The real error will be compared to the following error estimate:
542 // A + R * |y_i|
543 gsl_odeiv2_control * c = gsl_odeiv2_control_y_new(1e-6, 0.0);
544
545 //Allocate space for the evolutor
546 gsl_odeiv2_evolve * e = gsl_odeiv2_evolve_alloc(NumberOfRGEs);
547
548 //Definition of the RGE system (the Jacobian is not necessary for the RK4 method; it's an empty function here)
549 gsl_odeiv2_system RGEsystem = {RGEsZ2, JacobianZ2, NumberOfRGEs, &order};
550
551 //Set starting and end point as natural logarithmic scales (conversion from decadic log scale)
552 double t1 = Q1*log(10.0);
553 double t2 = Q2*log(10.0);
554 double tNLOuni = NLOuniscale*log(10.0);
555
556 //Set initial step size
557 double InitialStepSize = 1e-6;
558
559 //Run!
560 while (t1 < t2)
561 {
562 int status = gsl_odeiv2_evolve_apply (e, c, s, &RGEsystem, &t1, t2, &InitialStepSize, InitialValues);
563 if(status != GSL_SUCCESS) break;
564
565 //intermediate checks if appropriate
566 if(RGEcheckZ2(InitialValues,t1,Rpeps,tNLOuni) != 0) break;
567 }
568
569 gsl_odeiv2_evolve_free (e);
570 gsl_odeiv2_control_free (c);
571 gsl_odeiv2_step_free (s);
572
573 //Return the decadic log scale at which the evolution stopped
574 return t1/log(10.0);
575}
JacobianZ2
int JacobianZ2(double t, const double y[], double *dfdy, double dfdt[], void *order)
Definition: GeneralTHDMZ2Runner.cpp:397
RGEsZ2
int RGEsZ2(double t, const double y[], double beta[], void *flags)
Definition: GeneralTHDMZ2Runner.cpp:21
RGEcheckZ2
int RGEcheckZ2(const double InitialValues[], const double t1, const double Rpeps, const double tNLOuni)
Definition: GeneralTHDMZ2Runner.cpp:402
s
Test Observable.

◆ runGeneralTHDMZ2parameters()

void GeneralTHDMZ2Runner::runGeneralTHDMZ2parameters ( )
private

Definition at line 577 of file GeneralTHDMZ2Runner.cpp.

578{
579 std::string modelflag = myGTHDMZ2->getZ2ModelType();
580 std::string RGEorder = myGTHDMZ2->getRGEorderflag();
581 int flag;
582 //flag will be used to transport information about model and RGEorder to the Runner:
583 //flag=3*(0 for type I, 1 for type II, 2 for type X, 3 for type Y, 4 for inert) + (0 for LO, 1 for approxNLO and 2 for NLO)
584 if( RGEorder == "LO" )
585 flag = 0;
586 else if( RGEorder == "approxNLO" )
587 flag = 1;
588 else if( RGEorder == "NLO" )
589 flag = 2;
590 else {
591 throw std::runtime_error("GeneralTHDMZ2Runner: RGEorder can be only any of \"LO\", \"approxNLO\" or \"NLO\"");
592 }
593
594 double g1_at_MZ = sqrt(4.0*M_PI*Ale/cW2);
595 double g2_at_MZ = sqrt(4.0*M_PI*Ale/(1-cW2));
596 double g3_at_MZ = sqrt(4.0*M_PI*Als);
597
598 double Ytop_at_MZ = (sqrt(2.0)*myGTHDMZ2->getQuarks(QCD::TOP).getMass())/(vev*sinb);
599 double Ybottom1_at_MZ = 0.0;
600 double Ybottom2_at_MZ = 0.0;
601 double Ytau1_at_MZ = 0.0;
602 double Ytau2_at_MZ = 0.0;
603
604 /*link these to the SM values*/
605 double Mb_at_MZ = 2.96;//GeV
606 double Mtau_at_MZ = 1.75;//GeV
607
608 if( modelflag == "type1" ) {
609 Ybottom2_at_MZ = (sqrt(2.0)*Mb_at_MZ)/(vev*sinb);
610 Ytau2_at_MZ = (sqrt(2.0)*Mtau_at_MZ)/(vev*sinb);
611 }
612 else if( modelflag == "type2" ) {
613 Ybottom1_at_MZ = (sqrt(2.0)*Mb_at_MZ)/(vev*cosb);
614 Ytau1_at_MZ = (sqrt(2.0)*Mtau_at_MZ)/(vev*cosb);
615
616 flag += 3;
617 }
618 else if( modelflag == "typeX" ) {
619 Ybottom2_at_MZ = (sqrt(2.0)*Mb_at_MZ)/(vev*sinb);
620 Ytau1_at_MZ = (sqrt(2.0)*Mtau_at_MZ)/(vev*cosb);
621
622 flag += 6;
623 }
624 else if( modelflag == "typeY" ) {
625 Ybottom1_at_MZ = (sqrt(2.0)*Mb_at_MZ)/(vev*cosb);
626 Ytau2_at_MZ = (sqrt(2.0)*Mtau_at_MZ)/(vev*sinb);
627 flag += 9;
628 }
629 else if( modelflag == "inert" ) {
630 Ytop_at_MZ = 0.0;
631
632 flag += 12;
633 }
634 else {
635 throw std::runtime_error("modelflag can be only any of \"type1\", \"type2\", \"typeX\", \"typeY\", \"inert\"");
636 }
637
638 double lambda1_at_MZ = myGTHDMZ2->getlambda1_Z2();
639 double lambda2_at_MZ = myGTHDMZ2->getlambda2_Z2();
640 double lambda3_at_MZ = myGTHDMZ2->getlambda3_Z2();
641 double lambda4_at_MZ = myGTHDMZ2->getlambda4_Z2();
642 double lambda5_at_MZ = myGTHDMZ2->getlambda5_Z2();
643
644 double lambda345 = lambda3_at_MZ + lambda4_at_MZ + lambda5_at_MZ;
645
646 double m12_2_at_MZ = m12_2;
647 double m11_2_at_MZ = tanb*m12_2_at_MZ - vev*vev*(cosb*cosb*lambda1_at_MZ - sinb*sinb*lambda345)/2.;
648 double m22_2_at_MZ = m12_2_at_MZ/tanb - vev*vev*(sinb*sinb*lambda2_at_MZ - cosb*cosb*lambda345)/2.;
649
650 if(fabs(Q_GTHDM-log10(MZ))<0.005) //at MZ scale
651 {
652 Q_cutoff = log10(MZ);
653
654 g1_at_Q = g1_at_MZ;
655 g2_at_Q = g2_at_MZ;
656 g3_at_Q = g3_at_MZ;
657 Ytop_at_Q = Ytop_at_MZ;
658 Ybottom1_at_Q = Ybottom1_at_MZ;
659 Ybottom2_at_Q = Ybottom2_at_MZ;
660 Ytau1_at_Q = Ytau1_at_MZ;
661 Ytau2_at_Q = Ytau2_at_MZ;
662 m11_2_at_Q = m11_2_at_MZ;
663 m22_2_at_Q = m22_2_at_MZ;
664 m12_2_at_Q = m12_2_at_MZ;
665 lambda1_at_Q = lambda1_at_MZ;
666 lambda2_at_Q = lambda2_at_MZ;
667 lambda3_at_Q = lambda3_at_MZ;
668 lambda4_at_Q = lambda4_at_MZ;
669 lambda5_at_Q = lambda5_at_MZ;
670 }
671 else //at some other scale
672 {
673 double InitVals[14];
674 InitVals[0] = g1_at_MZ;
675 InitVals[1] = g2_at_MZ;
676 InitVals[2] = g3_at_MZ;
677 InitVals[3] = Ytop_at_MZ;
678 InitVals[4] = Ybottom1_at_MZ+Ybottom2_at_MZ;
679 InitVals[5] = Ytau1_at_MZ+Ytau2_at_MZ;
680 InitVals[6] = m11_2_at_MZ;
681 InitVals[7] = m22_2_at_MZ;
682 InitVals[8] = m12_2_at_MZ;
683 InitVals[9] = lambda1_at_MZ;
684 InitVals[10] = lambda2_at_MZ;
685 InitVals[11] = lambda3_at_MZ;
686 InitVals[12] = lambda4_at_MZ;
687 InitVals[13] = lambda5_at_MZ;
688
689 // Running up to Q_cutoff <= Q_GTHDM
690 Q_cutoff = RGEGeneralTHDMZ2Runner(InitVals, 14, log10(MZ), Q_GTHDM, flag);
691
692 g1_at_Q = InitVals[0];
693 g2_at_Q = InitVals[1];
694 g3_at_Q = InitVals[2];
695 Ytop_at_Q = InitVals[3];
696 Ybottom1_at_Q = 0.0;
697 Ybottom2_at_Q = 0.0;
698 Ytau1_at_Q = 0.0;
699 Ytau2_at_Q = 0.0;
700
701 if( modelflag == "type1" ) {
702 Ybottom2_at_Q = InitVals[4];
703 Ytau2_at_Q = InitVals[5];
704 }
705 else if( modelflag == "type2" ) {
706 Ybottom1_at_Q = InitVals[4];
707 Ytau1_at_Q = InitVals[5];
708 }
709 else if( modelflag == "typeX" ) {
710 Ybottom2_at_Q = InitVals[4];
711 Ytau1_at_Q = InitVals[5];
712 }
713 else if( modelflag == "typeY" ) {
714 Ybottom1_at_Q = InitVals[4];
715 Ytau2_at_Q = InitVals[5];
716 }
717 else if( modelflag == "inert" ) {
718 Ytop_at_Q = 0.0;
719 }
720 else {
721 throw std::runtime_error("modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
722 }
723
724 m11_2_at_Q = InitVals[6];
725 m22_2_at_Q = InitVals[7];
726 m12_2_at_Q = InitVals[8];
727 lambda1_at_Q = InitVals[9];
728 lambda2_at_Q = InitVals[10];
729 lambda3_at_Q = InitVals[11];
730 lambda4_at_Q = InitVals[12];
731 lambda5_at_Q = InitVals[13];
732 }
733}
GeneralTHDMZ2::getlambda2_Z2
double getlambda2_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:155
GeneralTHDMZ2::getlambda4_Z2
double getlambda4_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:174
GeneralTHDMZ2::getlambda5_Z2
double getlambda5_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:183
GeneralTHDMZ2::getlambda3_Z2
double getlambda3_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:165
GeneralTHDMZ2::getlambda1_Z2
double getlambda1_Z2() const
A getter for the parameter of the scalar potential .
Definition: GeneralTHDMZ2.h:145
GeneralTHDMZ2::getZ2ModelType
std::string getZ2ModelType() const
A getter for the chosen Z2-symmetric model type.
Definition: GeneralTHDMZ2.h:191
GeneralTHDMZ2Runner::RGEGeneralTHDMZ2Runner
virtual double RGEGeneralTHDMZ2Runner(double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order)
Definition: GeneralTHDMZ2Runner.cpp:532
QCD::TOP
@ TOP
Definition: QCD.h:328
lambda345
An observable class for the quartic Higgs potential coupling combination .
Definition: THDMquantities.h:497

Member Data Documentation

◆ Ale

double GeneralTHDMZ2Runner::Ale
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ alpha

double GeneralTHDMZ2Runner::alpha
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ Als

double GeneralTHDMZ2Runner::Als
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ beta

double GeneralTHDMZ2Runner::beta
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ bma

double GeneralTHDMZ2Runner::bma
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ cosb

double GeneralTHDMZ2Runner::cosb
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ cW2

double GeneralTHDMZ2Runner::cW2
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ g1_at_Q

double GeneralTHDMZ2Runner::g1_at_Q
private

Definition at line 113 of file GeneralTHDMZ2Runner.h.

◆ g2_at_Q

double GeneralTHDMZ2Runner::g2_at_Q
private

Definition at line 113 of file GeneralTHDMZ2Runner.h.

◆ g3_at_Q

double GeneralTHDMZ2Runner::g3_at_Q
private

Definition at line 113 of file GeneralTHDMZ2Runner.h.

◆ lambda1_at_Q

double GeneralTHDMZ2Runner::lambda1_at_Q
private

Definition at line 116 of file GeneralTHDMZ2Runner.h.

◆ lambda2_at_Q

double GeneralTHDMZ2Runner::lambda2_at_Q
private

Definition at line 116 of file GeneralTHDMZ2Runner.h.

◆ lambda3_at_Q

double GeneralTHDMZ2Runner::lambda3_at_Q
private

Definition at line 116 of file GeneralTHDMZ2Runner.h.

◆ lambda4_at_Q

double GeneralTHDMZ2Runner::lambda4_at_Q
private

Definition at line 116 of file GeneralTHDMZ2Runner.h.

◆ lambda5_at_Q

double GeneralTHDMZ2Runner::lambda5_at_Q
private

Definition at line 116 of file GeneralTHDMZ2Runner.h.

◆ m11_2_at_Q

double GeneralTHDMZ2Runner::m11_2_at_Q
private

Definition at line 115 of file GeneralTHDMZ2Runner.h.

◆ m12_2

double GeneralTHDMZ2Runner::m12_2
private

Definition at line 110 of file GeneralTHDMZ2Runner.h.

◆ m12_2_at_Q

double GeneralTHDMZ2Runner::m12_2_at_Q
private

Definition at line 115 of file GeneralTHDMZ2Runner.h.

◆ m22_2_at_Q

double GeneralTHDMZ2Runner::m22_2_at_Q
private

Definition at line 115 of file GeneralTHDMZ2Runner.h.

◆ mA2

double GeneralTHDMZ2Runner::mA2
private

Definition at line 110 of file GeneralTHDMZ2Runner.h.

◆ mHh2

double GeneralTHDMZ2Runner::mHh2
private

Definition at line 110 of file GeneralTHDMZ2Runner.h.

◆ mHl2

double GeneralTHDMZ2Runner::mHl2
private

Definition at line 110 of file GeneralTHDMZ2Runner.h.

◆ mHp2

double GeneralTHDMZ2Runner::mHp2
private

Definition at line 110 of file GeneralTHDMZ2Runner.h.

◆ myGTHDMZ2

const GeneralTHDMZ2* GeneralTHDMZ2Runner::myGTHDMZ2

Definition at line 99 of file GeneralTHDMZ2Runner.h.

◆ myZ2_at_Q

gslpp::matrix<double> GeneralTHDMZ2Runner::myZ2_at_Q
private

Definition at line 107 of file GeneralTHDMZ2Runner.h.

◆ MZ

double GeneralTHDMZ2Runner::MZ
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ MZ2

double GeneralTHDMZ2Runner::MZ2
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ NLOuniscale

double GeneralTHDMZ2Runner::NLOuniscale
private

Definition at line 112 of file GeneralTHDMZ2Runner.h.

◆ Q_cutoff

double GeneralTHDMZ2Runner::Q_cutoff
private

Definition at line 112 of file GeneralTHDMZ2Runner.h.

◆ Q_GTHDM

double GeneralTHDMZ2Runner::Q_GTHDM
private

Definition at line 112 of file GeneralTHDMZ2Runner.h.

◆ Rpeps

double GeneralTHDMZ2Runner::Rpeps
private

Definition at line 112 of file GeneralTHDMZ2Runner.h.

◆ sinb

double GeneralTHDMZ2Runner::sinb
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ tanb

double GeneralTHDMZ2Runner::tanb
private

Definition at line 111 of file GeneralTHDMZ2Runner.h.

◆ vev

double GeneralTHDMZ2Runner::vev
private

Definition at line 109 of file GeneralTHDMZ2Runner.h.

◆ WFRcomb1

double GeneralTHDMZ2Runner::WFRcomb1
private

Definition at line 117 of file GeneralTHDMZ2Runner.h.

◆ WFRcomb2

double GeneralTHDMZ2Runner::WFRcomb2
private

Definition at line 117 of file GeneralTHDMZ2Runner.h.

◆ WFRcomb3

double GeneralTHDMZ2Runner::WFRcomb3
private

Definition at line 117 of file GeneralTHDMZ2Runner.h.

◆ WFRcomb4

double GeneralTHDMZ2Runner::WFRcomb4
private

Definition at line 117 of file GeneralTHDMZ2Runner.h.

◆ Ybottom1_at_Q

double GeneralTHDMZ2Runner::Ybottom1_at_Q
private

Definition at line 114 of file GeneralTHDMZ2Runner.h.

◆ Ybottom2_at_Q

double GeneralTHDMZ2Runner::Ybottom2_at_Q
private

Definition at line 114 of file GeneralTHDMZ2Runner.h.

◆ Ytau1_at_Q

double GeneralTHDMZ2Runner::Ytau1_at_Q
private

Definition at line 114 of file GeneralTHDMZ2Runner.h.

◆ Ytau2_at_Q

double GeneralTHDMZ2Runner::Ytau2_at_Q
private

Definition at line 114 of file GeneralTHDMZ2Runner.h.

◆ Ytop_at_Q

double GeneralTHDMZ2Runner::Ytop_at_Q
private

Definition at line 114 of file GeneralTHDMZ2Runner.h.

The documentation for this class was generated from the following files: