void CKM_parameters::to_pdg_convention(Eigen::Matrix<double,3,3>& Vu,
Eigen::Matrix<double,3,3>& Vd,
Eigen::Matrix<double,3,3>& Uu,
Eigen::Matrix<double,3,3>& Ud)
{
Eigen::Matrix<double,3,3> ckm(Vu*Vd.adjoint());
to_pdg_convention(ckm, Vu, Vd, Uu, Ud);
}
CMSSM_soft_parameters MSSMcbs<Two_scale>::calc_beta() const
{
CMSSM_soft_parameters betas(CMSSM<Two_scale>::calc_beta());
if (get_loops() <= 2) return betas;
// 3-loop g3 beta function from http://www.liv.ac.uk/~dij/betas/
Eigen::Matrix<double,3,3> Yt = Yu.transpose();
Eigen::Matrix<double,3,3> Yb = Yd.transpose();
Eigen::Matrix<double,3,3> Yl = Ye.transpose();
Eigen::Matrix<double,3,3> Ytc = Yt.adjoint();
Eigen::Matrix<double,3,3> Ybc = Yb.adjoint();
Eigen::Matrix<double,3,3> Ylc = Yl.adjoint();
double a1 = Sqr(g1);
double a2 = Sqr(g2);
double a3 = Sqr(g3);
double bg33 =
+ Sqr(a1)*a3 * ( - 1702/75.0 )
+ a1*a2*a3 * ( - 3/5.0 )
+ a1*Sqr(a3) * ( 22/15.0 )
+ a1*a3 * ( - 44/15.0*(Ytc*Yt).trace() - 32/15.0*(Ybc*Yb).trace() )
+ Sqr(a2)*a3 * ( - 27 )
+ a2*Sqr(a3) * ( 6 )
+ a2*a3 * ( - 12*(Ytc*Yt).trace() - 12*(Ybc*Yb).trace() )
+ Power(a3,3) * ( 347/3.0 )
+ Sqr(a3) * ( - 104/3.0*(Ytc*Yt).trace() - 104/3.0*(Ybc*Yb).trace() )
+ a3 * ( 18*Sqr((Ytc*Yt).trace()) + 12*(Ytc*Yt*Ytc*Yt).trace()
+ 8*(Ytc*Yt*Ybc*Yb).trace()
+ 18*Sqr((Ybc*Yb).trace()) + 6*(Ybc*Yb).trace()*(Ylc*Yl).trace()
+ 12*(Ybc*Yb*Ybc*Yb).trace() );
betas.set_g3(betas.get_g3() + Power(oneOver16PiSqr,3) * g3 * bg33);
return betas;
}
static T get(const Eigen::Matrix<T, N, 1> &x, const Eigen::Matrix<T, N, 1> &y) {
return x.adjoint() * y;
}
static return_type get(const Eigen::Matrix<T, N, M> &x) {
return x.adjoint();
}