20 return 0.0000001*m12_2*(m12_2 +(mHl*mHl-mHh2)*cosa*sina
21 +
sqrt((mHl*mHl*cosa*cosa +mHh2*sina*sina -m12_2/tanb)
22 *(mHh2*cosa*cosa +mHl*mHl*sina*sina -m12_2*tanb)));
double computeThValue()
A condition to assure that our vacuum is global.
double getMHl() const
A get method to retrieve the Higgs mass .
complex sqrt(const complex &z)