easier sign convention for tadpoles

FlexibleSUSY · Aug 10, 2016 · 37873a7 · 37873a7
1 parent 2200ec7
commit 37873a7
Show file tree

Hide file tree

Showing 3 changed files with 15 additions and 21 deletions.
diff --git a/meta/SelfEnergies.m b/meta/SelfEnergies.m
@@ -608,8 +608,8 @@
       costau, scalesq, amu, tanb, vev2);
 }
 
-tadpole_2l(0) *= -" <> CConversion`ToValidCSymbolString[SARAH`VEVSM1] <> ";
-tadpole_2l(1) *= -" <> CConversion`ToValidCSymbolString[SARAH`VEVSM2] <> ";
+tadpole_2l(0) *= " <> CConversion`ToValidCSymbolString[SARAH`VEVSM1] <> ";
+tadpole_2l(1) *= " <> CConversion`ToValidCSymbolString[SARAH`VEVSM2] <> ";
 
 if (!IsFinite(tadpole_2l))
    tadpole_2l.setZero();
@@ -697,9 +697,9 @@
       costau, scalesq, amu, tanb, vev2);
 }
 
-tadpole_2l(0) *= -" <> CConversion`ToValidCSymbolString[SARAH`VEVSM1] <> ";
-tadpole_2l(1) *= -" <> CConversion`ToValidCSymbolString[SARAH`VEVSM2] <> ";
-tadpole_2l(2) *= -" <>  svevStr <> ";
+tadpole_2l(0) *= " <> CConversion`ToValidCSymbolString[SARAH`VEVSM1] <> ";
+tadpole_2l(1) *= " <> CConversion`ToValidCSymbolString[SARAH`VEVSM2] <> ";
+tadpole_2l(2) *= " <>  svevStr <> ";
 
 if (!IsFinite(tadpole_2l))
    tadpole_2l.setZero();

diff --git a/src/mssm_twoloophiggs.cpp b/src/mssm_twoloophiggs.cpp
@@ -45,7 +45,7 @@ Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_at_as_mssm(
    ewsb2loop_(&rmtsq, &mg, &mst1sq, &mst2sq, &sxt, &cxt, &scalesq,
               &amu, &tanb, &vev2, &gs, &result(0), &result(1));
 
-   return result;
+   return -result;
 }
 
 Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_at_at_mssm(
@@ -62,7 +62,7 @@ Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_at_at_mssm(
            &sxt, &cxt, &sxb, &cxb, &scalesq, &amu, &tanb, &vev2,
            &result(0), &result(1));
 
-   return result;
+   return -result;
 }
 
 Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_ab_as_mssm(
@@ -91,7 +91,7 @@ Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_atau_atau_mssm(
    tausqtad_(&rmtausq, &mAsq, &msnusq, &mstau1sq, &mstau2sq, &sintau,
              &costau, &scalesq, &amu, &tanb, &vev2, &result(0), &result(1));
 
-   return result;
+   return -result;
 }
 
 Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_as_mssm_with_tadpoles(
@@ -243,10 +243,10 @@ Eigen::Matrix<double, 2, 2> rotate_scalar(
 }
 
 Eigen::Matrix<double, 2, 2> subtract_mssm_tadpoles_scalar(
-   double dMA, const Eigen::Matrix<double, 2, 1>& tadpoles,
+   double self_energy, const Eigen::Matrix<double, 2, 1>& tadpoles,
    double tanb)
 {
-   return rotate_scalar(dMA, tanb) - Eigen::Matrix<double, 2, 2>(tadpoles.asDiagonal());
+   return rotate_scalar(self_energy, tanb) + Eigen::Matrix<double, 2, 2>(tadpoles.asDiagonal());
 }
 
 Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_at_as_mssm(
@@ -369,13 +369,7 @@ Eigen::Matrix<double, 2, 2> subtract_mssm_tadpoles_pseudoscalar(
    double self_energy, const Eigen::Matrix<double, 2, 1>& tadpoles,
    double tanb)
 {
-   Eigen::Matrix<double, 2, 2> result =
-      rotate_pseudoscalar(self_energy, tanb);
-
-   result(0,0) -= tadpoles(0);
-   result(1,1) -= tadpoles(1);
-
-   return result;
+   return rotate_pseudoscalar(self_energy, tanb) + Eigen::Matrix<double, 2, 2>(tadpoles.asDiagonal());
 }
 
 Eigen::Matrix<double, 2, 2> self_energy_pseudoscalar_2loop_at_as_mssm(

diff --git a/src/nmssm_twoloophiggs.cpp b/src/nmssm_twoloophiggs.cpp
@@ -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(
@@ -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(
@@ -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(
@@ -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(