Shift soft Higgs masses to meet tree-level EWSB conditions in lattice…

… code to respect the tree-level Slavnov-Taylor identity between the goldstone boson and vector boson masses (in Rxi gauge for xi = 1). See also: commit b4cf167, commit 933400f
FlexibleSUSY · Sep 5, 2014 · e61651b · e61651b
1 parent 4a5217e
commit e61651b
Show file tree

Hide file tree

Showing 3 changed files with 23 additions and 1 deletion.
diff --git a/meta/Lattice.m b/meta/Lattice.m
@@ -150,7 +150,7 @@
 	vevs = Union @ Variables @ vevRules[[All,2]],
 	treeEwsbConstraints = fsEwsbEquations /. sarahOperatorReplacementRules,
 	softHiggsMasses,
-	treeEwsbEquations,
+	treeEwsbEquations, shiftHiggsMasses,
 	ewsbConstraints, ewsbEquations, ewsbDep, ewsbList,
 	fixTsusy, tsusyConstraint = (Exp[t] scl0)^4 - Lattice`Private`M2[Global`Su[{1}]] Lattice`Private`M2[Global`Su[{6}]] /. sarahOperatorReplacementRules
     },
@@ -208,6 +208,7 @@
     treeEwsbEquations = ParametrizeEWSBEquations[
 	treeEwsbConstraints, ConditionPositiveVevs[vevs], softHiggsMasses,
 	parameterRules];
+    shiftHiggsMasses = CShiftHiggsMasses[treeEwsbEquations];
     ewsbConstraints = EWSBConstraintsWithCorrections[treeEwsbConstraints];
     ewsbEquations = ParametrizeEWSBEquations[
 	ewsbConstraints, ConditionPositiveVevs[vevs], softHiggsMasses,
@@ -247,6 +248,7 @@
 	"@nPointDecls@"	    -> IndentText[nPointDecls, 4],
 	"@nPointDefs@"	    -> WrapText[StringJoin@nPointDefs],
 	"@phaseDefs@"	    -> IndentText[phaseDefs, 4],
+	"@shiftHiggsMasses@"-> WrapText@IndentText[shiftHiggsMasses, 2],
 	"@vertexDecls@"	    -> IndentText[vertexDecls, 4],
 	"@vertexDefs@"	    -> WrapText[StringJoin@vertexDefs]
     }]];
@@ -291,6 +293,16 @@
 
 EWSBConditionToC[lhs_ == rhs_] := NConstraintToC[rhs - lhs];
 
+CShiftHiggsMasses[treeEwsbEquations_List] :=
+    StringJoin[CShiftHiggsMass /@ treeEwsbEquations];
+
+CShiftHiggsMass[lhs_ == rhs_] := Module[{
+	clhs = CExpToCFormString @ ToCExp[lhs, x],
+	crhs = CExpToCFormString @ ToCExp[rhs, x]
+    },
+    clhs <> " = " <> crhs <> ";\n"
+];
+
 NConstraintToC[constraint_] :=
     CNConstraint[
 	Dependence -> DependenceList[constraint],

diff --git a/templates/lattice_model.hpp.in b/templates/lattice_model.hpp.in
@@ -68,6 +68,7 @@ public:
 @phaseDefs@
 @dependenceNumDecls@
   private:
+    void eliminate_soft_higgs_masses_at_tree_level(Eigen::VectorXd& x);
   };
 
   // beta functions

diff --git a/templates/lattice_model_interactions.cpp.in b/templates/lattice_model_interactions.cpp.in
@@ -39,10 +39,19 @@ void CLASSNAME::Interactions::set(const Eigen::VectorXd& x_, double scl0_)
 {
   Lattice_model_interactions::set(x_, scl0_);
 
+  Eigen::VectorXd x = Lattice_model_interactions::x;
+  eliminate_soft_higgs_masses_at_tree_level(x);
+
 @matrixStmts@
 @eigenVarStmts@
 }
 
+void CLASSNAME::Interactions::eliminate_soft_higgs_masses_at_tree_level
+(Eigen::VectorXd& x)
+{
+@shiftHiggsMasses@
+}
+
 @nPointDefs@
 @vertexDefs@
 @dependenceNumDefs@