Skip to content

Commit

Permalink
make all functions return same sign convention
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Aug 10, 2016
1 parent 1f58f1f commit 2200ec7
Show file tree
Hide file tree
Showing 2 changed files with 24 additions and 24 deletions.
32 changes: 16 additions & 16 deletions src/mssm_twoloophiggs.cpp
View file
Expand Up @@ -109,7 +109,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_as_mssm_with_tadpoles(

result(1,0) = result(0,1);

return result;
return -result;
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_at_mssm_with_tadpoles(
Expand All @@ -128,7 +128,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_at_mssm_with_tadpoles(

result(1,0) = result(0,1);

return result;
return -result;
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm_with_tadpoles(
Expand Down Expand Up @@ -160,7 +160,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm_with_tadpoles

result(1,0) = result(0,1);

return result;
return -result;
}

double self_energy_pseudoscalar_2loop_at_as_mssm_with_tadpoles(
Expand All @@ -175,7 +175,7 @@ double self_energy_pseudoscalar_2loop_at_as_mssm_with_tadpoles(
dszodd_(&rmtsq, &mg, &mst1sq, &mst2sq, &sxt, &cxt, &scalesq, &amu,
&tanb, &vev2, &gs, &result);

return result;
return -result;
}

double self_energy_pseudoscalar_2loop_at_at_mssm_with_tadpoles(
Expand All @@ -191,7 +191,7 @@ double self_energy_pseudoscalar_2loop_at_at_mssm_with_tadpoles(
ddsodd_(&rmtsq, &rmbsq, &fmasq, &mst1sq, &mst2sq, &msb1sq, &msb2sq,
&sxt, &cxt, &sxb, &cxb, &scalesq, &amu, &tanb, &vev2, &result);

return result;
return -result;
}

double self_energy_pseudoscalar_2loop_ab_as_mssm_with_tadpoles(
Expand All @@ -216,7 +216,7 @@ double self_energy_pseudoscalar_2loop_atau_atau_mssm_with_tadpoles(
tausqodd_(&rmtausq, &fmasq, &msnusq, &mstau1sq, &mstau2sq, &sintau,
&costau, &scalesq, &amu, &tanb, &vev2, &result);

return result;
return -result;
}

// self-energies without tadpoles
Expand All @@ -234,10 +234,10 @@ Eigen::Matrix<double, 2, 2> rotate_scalar(

Eigen::Matrix<double, 2, 2> result;

result(0,0) = - self_energy * sqr(sinb);
result(0,1) = + self_energy * sinb * cosb;
result(0,0) = self_energy * sqr(sinb);
result(0,1) = - self_energy * sinb * cosb;
result(1,0) = result(0,1);
result(1,1) = - self_energy * sqr(cosb);
result(1,1) = self_energy * sqr(cosb);

return result;
}
Expand Down Expand Up @@ -270,7 +270,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_as_mssm(
const Eigen::Matrix<double, 2, 2> tM =
subtract_mssm_tadpoles_scalar(dMA, tadpoles, tanb);

return -result + tM;
return result + tM;
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_at_mssm(
Expand All @@ -296,7 +296,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_at_mssm(
const Eigen::Matrix<double, 2, 2> tM =
subtract_mssm_tadpoles_scalar(dMA, tadpoles, tanb);

return -result + tM;
return result + tM;
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm(
Expand All @@ -319,7 +319,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm(
const Eigen::Matrix<double, 2, 2> tM =
subtract_mssm_tadpoles_scalar(dMA, tadpoles, 1./cotbeta);

return -result + tM;
return result + tM;
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm(
Expand All @@ -344,7 +344,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm(
const Eigen::Matrix<double, 2, 2> tM =
subtract_mssm_tadpoles_scalar(dMA, tadpoles, tanb);

return -result + tM;
return result + tM;
}

Eigen::Matrix<double, 2, 2> rotate_pseudoscalar(
Expand All @@ -357,10 +357,10 @@ Eigen::Matrix<double, 2, 2> rotate_pseudoscalar(
Eigen::Matrix<double, 2, 2> result;

// see hep-ph/0105096 Eq. (9)
result(0,0) = - self_energy * sqr(sinb);
result(0,1) = - self_energy * sinb * cosb;
result(0,0) = self_energy * sqr(sinb);
result(0,1) = self_energy * sinb * cosb;
result(1,0) = result(0,1);
result(1,1) = - self_energy * sqr(cosb);
result(1,1) = self_energy * sqr(cosb);

return result;
}
Expand Down
16 changes: 8 additions & 8 deletions src/nmssm_twoloophiggs.cpp
View file
Expand Up @@ -90,7 +90,7 @@ Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_at_as_nmssm(
tadpole_higgs_2loop_at_as_nmssm(
sqr(rmt), mg, mst1sq, mst2sq, sxt, cxt, scalesq, amu, tanb, vev2, gs, svevS).asDiagonal();

return -se - tadpoles;
return se - tadpoles;
}

Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_ab_as_nmssm(
Expand All @@ -109,7 +109,7 @@ Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_ab_as_nmssm(
tadpole_higgs_2loop_ab_as_nmssm(
sqr(rmb), mg, msb1sq, msb2sq, sxb, cxb, scalesq, amu, cotb, vev2, gs, svevS).asDiagonal();

return -se - tadpoles;
return se - tadpoles;
}

Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_at_as_nmssm(
Expand All @@ -129,7 +129,7 @@ Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_at_as_nmssm(
tadpole_higgs_2loop_at_as_nmssm(
sqr(rmt), mg, mst1sq, mst2sq, sxt, cxt, scalesq, amu, tanb, vev2, gs, svevS).asDiagonal();

return -se - tadpoles;
return se - tadpoles;
}

Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_ab_as_nmssm(
Expand All @@ -149,7 +149,7 @@ Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_ab_as_nmssm(
tadpole_higgs_2loop_ab_as_nmssm(
sqr(rmb), mg, msb1sq, msb2sq, sxb, cxb, scalesq, amu, cotb, vev2, gs, svevS).asDiagonal();

return -se - tadpoles;
return se - tadpoles;
}

Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_at_as_nmssm_with_tadpoles(
Expand All @@ -170,7 +170,7 @@ Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_at_as_nmssm_with_tadpoles(
DMS[1][0], DMS[1][1], DMS[1][2],
DMS[2][0], DMS[2][1], DMS[2][2];

return result;
return -result;
}

Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_ab_as_nmssm_with_tadpoles(
Expand All @@ -196,7 +196,7 @@ Eigen::Matrix<double, 3, 3> self_energy_higgs_2loop_ab_as_nmssm_with_tadpoles(
DMS[1][0], DMS[1][1], DMS[1][2],
DMS[2][0], DMS[2][1], DMS[2][2];

return result;
return -result;
}

Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_at_as_nmssm_with_tadpoles(
Expand All @@ -217,7 +217,7 @@ Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_at_as_nmssm_with_tadp
DMP[1][0], DMP[1][1], DMP[1][2],
DMP[2][0], DMP[2][1], DMP[2][2];

return result;
return -result;
}

Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_ab_as_nmssm_with_tadpoles(
Expand All @@ -243,7 +243,7 @@ Eigen::Matrix<double, 3, 3> self_energy_pseudoscalar_2loop_ab_as_nmssm_with_tadp
DMP[1][0], DMP[1][1], DMP[1][2],
DMP[2][0], DMP[2][1], DMP[2][2];

return result;
return -result;
}

} // namespace flexiblesusy

0 comments on commit 2200ec7

Please sign in to comment.