help FS finding the soft Higgs masses

I introduces the new user variable FSSoftHiggsMasses where the user
can specify the list of soft Higgs masses.

For the calculation of the Higgs masses, the tree-level EWSB
eqs. substitution must be done in the Higgs mass matrix.  In FS this
substitution is done like this:

 0) save the soft Higgs masses in temporary varibles

 1) solve the tree-level EWSB eqs. for the soft Higgs masses
    (thereby the soft Higgs masse are overwritten)

 2) calculate the Higgs mass matrix

 3) reset the soft Higgs masses to the saved values from 0)

For these steps, FS needs to know which model parameters are the soft
Higgs masses.  In some models (for example in the munuSSM), this is
difficult to find out.  If FS fails to find the soft Higgs masses, the
user can assist by setting them in the variable FSSoftHiggsMasses.

meta/EWSB.m

@@ -383,7 +383,7 @@
            If[Length[equations] =!= Length[parameters],
               Print["Error: SolveTreeLevelEwsbVia: trying to solve ",
                     Length[equations], " equations for ", Length[parameters],
-                    " parameters ", parameters];
+                    " parameters ", InputForm[parameters]];
               Quit[1];
              ];
            simplifiedEqs = (# == 0)& /@ equations;

meta/FlexibleSUSY.m

@@ -70,6 +70,9 @@
 FSSimplifyBetaFunctionsTimeConstraint = 120;
 FSSolveWeinbergAngleTimeConstraint = 120;
 
+(* list of soft breaking Higgs masses for solving EWSB eqs. *)
+FSSoftHiggsMasses = {};
+
 (* EWSB solvers *)
 GSLHybrid;   (* hybrid method *)
 GSLHybridS;  (* hybrid method with dynamic step size *)
@@ -221,6 +224,11 @@
            If[Head[FlexibleSUSY`ExtraSLHAOutputBlocks] =!= List,
               FlexibleSUSY`ExtraSLHAOutputBlocks = {};
              ];
+           If[Head[FlexibleSUSY`EWSBOutputParameters] =!= List,
+              Print["Error: EWSBOutputParameters has to be set to a list",
+                    " of model parameters chosen to be output of the EWSB eqs."];
+              Quit[1];
+             ];
           ];
 
 ReplaceIndicesInUserInput[] :=
@@ -582,8 +590,14 @@
               softScalarMasses          = DeleteDuplicates[SARAH`ListSoftBreakingScalarMasses];,
               softScalarMasses          = {};
              ];
-           softHiggsMasses              = Select[softScalarMasses, (!FreeQ[ewsbEquations, #])&];
-           softHiggsMasses              = Parameters`ExpandExpressions[Parameters`AppendGenerationIndices[softHiggsMasses]];
+           (* find soft Higgs masses that appear in tree-level EWSB eqs. *)
+           If[Head[FlexibleSUSY`FSSoftHiggsMasses] =!= List ||
+              FlexibleSUSY`FSSoftHiggsMasses === {},
+              softHiggsMasses = Select[softScalarMasses, (!FreeQ[ewsbEquations, #])&];
+              ,
+              softHiggsMasses = FlexibleSUSY`FSSoftHiggsMasses;
+             ];
+           softHiggsMasses              = Parameters`DecreaseIndexLiterals[Parameters`ExpandExpressions[Parameters`AppendGenerationIndices[softHiggsMasses]]];
            saveSoftHiggsMasses          = Parameters`SaveParameterLocally[softHiggsMasses, "old_", ""];
            restoreSoftHiggsMasses       = Parameters`RestoreParameter[softHiggsMasses, "old_", ""];
            If[Head[softHiggsMasses] === List && Length[softHiggsMasses] > 0,
@@ -1213,12 +1227,13 @@
              ];
 
            ewsbEquations = Parameters`ExpandExpressions[ewsbEquations];
+           FlexibleSUSY`EWSBOutputParameters = Parameters`DecreaseIndexLiterals[FlexibleSUSY`EWSBOutputParameters];
 
            If[FlexibleSUSY`TreeLevelEWSBSolution === {},
               (* trying to find an analytic solution for the EWSB eqs. *)
               treeLevelEwsbOutputFile = FileNameJoin[{Global`$flexiblesusyOutputDir,
                                                       FlexibleSUSY`FSModelName <> "_tree_level_EWSB_solution.m"}];
-              Print["Solving EWSB equations ..."];
+              Print["Solving EWSB equations for ", FlexibleSUSY`EWSBOutputParameters," ..."];
               {ewsbSolution, freePhases} = EWSB`FindSolutionAndFreePhases[ewsbEquations,
                                                                           FlexibleSUSY`EWSBOutputParameters,
                                                                           treeLevelEwsbOutputFile];

model_files/munuSSM/FlexibleSUSY.m.in

@@ -23,6 +23,8 @@ EXTPAR = {
 
 EWSBOutputParameters = {mlHd2[1], mlHd2[2], mlHd2[3], mHd2, mHu2, mv2};
 
+FSSoftHiggsMasses = EWSBOutputParameters;
+
 SUSYScale = Sqrt[Product[M[Su[i]]^(Abs[ZU[i,3]]^2 + Abs[ZU[i,6]]^2), {i,6}]];
 
 SUSYScaleFirstGuess = Sqrt[m0^2 + 4 m12^2];