support MZMSbar and MWMSbar symbols for non-SUSY models

FlexibleSUSY · Mar 9, 2017 · e939647 · e939647
1 parent 3f5efc6
commit e939647
Show file tree

Hide file tree

Showing 27 changed files with 48 additions and 42 deletions.
diff --git a/meta/FlexibleSUSY.m b/meta/FlexibleSUSY.m
@@ -106,6 +106,8 @@ FlexibleSUSY model file (FlexibleSUSY.m).
 MT;
 MZDRbar;
 MWDRbar;
+MZMSbar;
+MWMSbar;
 EDRbar;
 ThetaWDRbar;
 SCALE;

diff --git a/model_files/BLSM/FlexibleSUSY.m.in b/model_files/BLSM/FlexibleSUSY.m.in
@@ -40,7 +40,7 @@ LowScaleInput = {
     {Yv, LHInput[Yv]},
     {Yx, LHInput[Yx]},
     {vX, vXinput},
-    {vH, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+    {vH, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
     {Ye, Automatic},
     {Yd, Automatic},
     {Yu, Automatic}

diff --git a/model_files/BLSMlightZp/FlexibleSUSY.m.in b/model_files/BLSMlightZp/FlexibleSUSY.m.in
@@ -40,7 +40,7 @@ LowScaleInput = {
     {Yv, LHInput[Yv]},
     {Yx, LHInput[Yx]},
     {vX, vXinput},
-    {vH, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+    {vH, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
     {Ye, Automatic},
     {Yd, Automatic},
     {Yu, Automatic}

diff --git a/model_files/BetaSM/FlexibleSUSY.m.in b/model_files/BetaSM/FlexibleSUSY.m.in
@@ -35,7 +35,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Ye, Automatic}
 };

diff --git a/model_files/HGTHDMIIMSSMBC/FlexibleSUSY.m.in b/model_files/HGTHDMIIMSSMBC/FlexibleSUSY.m.in
@@ -80,8 +80,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/HGTHDMIIMSSMBCApprox/FlexibleSUSY.m.in b/model_files/HGTHDMIIMSSMBCApprox/FlexibleSUSY.m.in
@@ -140,8 +140,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/HGTHDMIIMSSMBCFull/FlexibleSUSY.m.in b/model_files/HGTHDMIIMSSMBCFull/FlexibleSUSY.m.in
@@ -130,8 +130,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/HSSUSY/FlexibleSUSY.m.in b/model_files/HSSUSY/FlexibleSUSY.m.in
@@ -74,7 +74,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/HTHDMIIMSSMBC/FlexibleSUSY.m.in b/model_files/HTHDMIIMSSMBC/FlexibleSUSY.m.in
@@ -70,8 +70,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/MDM/FlexibleSUSY.m.in b/model_files/MDM/FlexibleSUSY.m.in
@@ -33,7 +33,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},

diff --git a/model_files/NSM/FlexibleSUSY.m.in b/model_files/NSM/FlexibleSUSY.m.in
@@ -38,7 +38,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {vH, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {vH, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SM/FlexibleSUSY.m.in b/model_files/SM/FlexibleSUSY.m.in
@@ -33,7 +33,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SMEWSBAtMZ/FlexibleSUSY.m.in b/model_files/SMEWSBAtMZ/FlexibleSUSY.m.in
@@ -33,7 +33,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},

diff --git a/model_files/SMHighPrecision/FlexibleSUSY.m.in b/model_files/SMHighPrecision/FlexibleSUSY.m.in
@@ -33,7 +33,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SMHighScale/FlexibleSUSY.m.in b/model_files/SMHighScale/FlexibleSUSY.m.in
@@ -35,7 +35,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[gY]^2 gY^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[gY]^2 gY^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SMNoRGEs/FlexibleSUSY.m.in b/model_files/SMNoRGEs/FlexibleSUSY.m.in
@@ -34,7 +34,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SMnoGUT/FlexibleSUSY.m.in b/model_files/SMnoGUT/FlexibleSUSY.m.in
@@ -22,7 +22,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[gY]^2 gY^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[gY]^2 gY^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SMtower/FlexibleSUSY.m.in b/model_files/SMtower/FlexibleSUSY.m.in
@@ -24,7 +24,7 @@ LowScale = LowEnergyConstant[MT];
 LowScaleFirstGuess = LowEnergyConstant[MT];
 
 LowScaleInput = {
-    {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+    {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
     {Yu, Automatic},
     {Yd, Automatic},
     {Ye, Automatic},

diff --git a/model_files/SSM/FlexibleSUSY.m.in b/model_files/SSM/FlexibleSUSY.m.in
@@ -42,7 +42,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/SplitMSSM/FlexibleSUSY.m.in b/model_files/SplitMSSM/FlexibleSUSY.m.in
@@ -63,7 +63,7 @@ LowScaleInput = {
    {MassB , M1Input},
    {MassWB, M2Input},
    {MassG , M3Input},
-   {v, Sqrt[2] MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, Sqrt[2] MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/THDMII/FlexibleSUSY.m.in b/model_files/THDMII/FlexibleSUSY.m.in
@@ -45,8 +45,8 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]},
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/model_files/THDMIIMSSMBC/FlexibleSUSY.m.in b/model_files/THDMIIMSSMBC/FlexibleSUSY.m.in
@@ -69,8 +69,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/THDMIIMSSMBCApprox/FlexibleSUSY.m.in b/model_files/THDMIIMSSMBCApprox/FlexibleSUSY.m.in
@@ -127,8 +127,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/THDMIIMSSMBCFull/FlexibleSUSY.m.in b/model_files/THDMIIMSSMBCFull/FlexibleSUSY.m.in
@@ -117,8 +117,8 @@ LowScaleInput = {
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic},
-   {v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
-   {v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
+   {v1, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
+   {v2, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
 };
 
 InitialGuessAtLowScale = {

diff --git a/model_files/cSM/FlexibleSUSY.m.in b/model_files/cSM/FlexibleSUSY.m.in
@@ -24,7 +24,7 @@ LowScale = LowEnergyConstant[MZ];
 LowScaleFirstGuess = LowEnergyConstant[MZ];
 
 LowScaleInput = {
-   {v, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
+   {v, 2 MZMSbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2]},
    {Yu, Automatic},
    {Yd, Automatic},
    {Ye, Automatic}

diff --git a/templates/two_scale_low_scale_constraint.cpp.in b/templates/two_scale_low_scale_constraint.cpp.in
@@ -53,6 +53,10 @@ namespace flexiblesusy {
 #define SCALE model->get_scale()
 #define MODEL model
 #define MODELCLASSNAME @ModelName@<Two_scale>
+#define MWMSbar mWrun
+#define MWDRbar mWrun
+#define MZMSbar mZrun
+#define MZDRbar mZrun
 #define CKM ckm
 #define PMNS pmns
 #define THETAW theta_w
@@ -135,8 +139,8 @@ void @ModelName@_low_scale_constraint<Two_scale>::clear()
    ckm.setIdentity();
    pmns.setIdentity();
    neutrinoDRbar.setZero();
-   MWDRbar = 0.;
-   MZDRbar = 0.;
+   mWrun = 0.;
+   mZrun = 0.;
    AlphaS = 0.;
    EDRbar = 0.;
    ThetaWDRbar = 0.;
@@ -156,8 +160,8 @@ void @ModelName@_low_scale_constraint<Two_scale>::initialize()
    ckm = qedqcd.get_complex_ckm();
    pmns = qedqcd.get_complex_pmns();
    neutrinoDRbar = Eigen::Matrix<double,3,3>::Zero();
-   MWDRbar = 0.;
-   MZDRbar = 0.;
+   mWrun = 0.;
+   mZrun = 0.;
    AlphaS = 0.;
    EDRbar = 0.;
    ThetaWDRbar = 0.;
@@ -203,12 +207,12 @@ void @ModelName@_low_scale_constraint<Two_scale>::calculate_threshold_correction
    const double e_drbar        = Sqrt(4.0 * Pi * alpha_em_drbar);
 
    // interface variables
-   MZDRbar = mz_pole;
-   MWDRbar = mw_pole;
+   mZrun = mz_pole;
+   mWrun = mw_pole;
 
    if (model->get_thresholds()) {
-      MZDRbar = model->calculate_M@VectorZ@_DRbar(mz_pole);
-      MWDRbar = model->calculate_M@VectorW@_DRbar(mw_pole);
+      mZrun = model->calculate_M@VectorZ@_DRbar(mz_pole);
+      mWrun = model->calculate_M@VectorW@_DRbar(mw_pole);
    }
 
    AlphaS = alpha_s_drbar;

diff --git a/templates/two_scale_low_scale_constraint.hpp.in b/templates/two_scale_low_scale_constraint.hpp.in
@@ -61,8 +61,8 @@ private:
    Eigen::Matrix<std::complex<double>,3,3> ckm{Eigen::Matrix<std::complex<double>,3,3>::Identity()};
    Eigen::Matrix<std::complex<double>,3,3> pmns{Eigen::Matrix<std::complex<double>,3,3>::Identity()};
    Eigen::Matrix<double,3,3> neutrinoDRbar{Eigen::Matrix<double,3,3>::Zero()};
-   double MWDRbar{0.};
-   double MZDRbar{0.};
+   double mWrun{0.};
+   double mZrun{0.};
    double AlphaS{0.};
    double EDRbar{0.};
    double ThetaWDRbar{0.};