adding complete 1L bottom/tau and 2L O((at+ab)^2 + atau^2) corrections

from [arXiv:1703.08166].

Many thanks to Emanuele Bagnaschi and Pietro Slavich for providing the
expressions!

.gitattributes

@@ -14,8 +14,6 @@ examples/lattice_numerical_fmssm.cpp        export-ignore
 examples/lattice_numerical_fmssm_fmssmn.cpp export-ignore
 examples/switch_MSSM.cpp                    export-ignore
 
-meta/SM/HSSUSY_corrections.m    export-ignore
-
 model_files/BetaSM              export-ignore
 model_files/cCMSSM              export-ignore
 model_files/complexMSSM         export-ignore

ChangeLog

@@ -252,6 +252,39 @@ FlexibleSUSY 2.0.0 [not released yet]
           higgs3loopCorrectionAbAsAs -> 1,   (* FlexibleSUSY[27] *)
       }
 
+ * Feature: Adding complete 1-loop O(alpha_tau + alpha_b) and complete
+   2-loop O((alpha_t + alpha_b)^2 + atau^2) threshold corrections for
+   lambda(MSUSY) to the HSSUSY model file from [arXiv:1703.08166].
+   Many thanks to Pietro Slavich and Emanuele Bagnaschi for providing
+   the expressions.
+
+   Note: 5 new flags are introduced to enable/disable the individual
+   2-loop corrections.  In the SLHA input the flags which control the
+   inclusion of 2-loop corrections are:
+
+      Block EXTPAR                 # Input parameters
+        100   2                    # LambdaLoopOrder
+        101   1                    # TwoLoopAtAs
+        102   1                    # TwoLoopAbAs
+        103   1                    # TwoLoopAtAb
+        104   1                    # TwoLoopAtauAtau
+        105   1                    # TwoLoopAtAt
+
+   In the Mathematica interface the flags which control the inclusion
+   of 2-loop corrections are:
+
+      handle = FSHSSUSYOpenHandle[
+         fsModelParameters -> {
+            ...
+            LambdaLoopOrder -> 2,
+            TwoLoopAtAs -> 1,
+            TwoLoopAbAs -> 1,
+            TwoLoopAtAb -> 1,
+            TwoLoopAtauAtau -> 1,
+            TwoLoopAtAt -> 1
+         }
+      ]
+
  * Change: The interface for adding constraints and matching
    conditions to the Two_scale_solver class has been simplified.
    Constrains and matching conditions are now added using the add()

doc/hssusy.dox

@@ -9,8 +9,8 @@ Model, matched to the MSSM at the SUSY scale, \f$M_\text{SUSY}\f$.  At
 the SUSY scale, the quartic Higgs coupling,
 \f$\lambda(M_\text{SUSY})\f$ is predicted from the matching to the
 MSSM using the threshold corrections of [arXiv:1407.4081,
-arxiv:1504.05200].  The 3-loop renormalization group equations of
-[arxiv:1504.05200, arxiv:1303.4364] are used to run
+arxiv:1504.05200, arXiv:1703.08166].  The 3-loop renormalization group
+equations of [arxiv:1504.05200, arxiv:1303.4364] are used to run
 \f$\lambda(M_\text{SUSY})\f$ down to the electroweak scale.  At the
 electroweak scale, the gauge and Yukawa couplings as well as the SM
 vacuum expectation value (VEV) are calculated at the 1-loop level from
@@ -56,7 +56,8 @@ _High scale_
 The high scale (`HighScale`) is fixed \f$M_{\text{SUSY}}\f$.  At
 \f$M_{\text{SUSY}}\f$, the quartic Higgs coupling
 \f$\lambda(M_\text{SUSY})\f$ is fixed using the 1- and leading 2-loop
-threshold corrections of [arXiv:1407.4081, arxiv:1504.05200].
+threshold corrections of [arXiv:1407.4081, arxiv:1504.05200,
+arXiv:1703.08166].
 
 _EWSB scale_
 

doc/librarylink.dox

@@ -69,7 +69,6 @@ CalcMh[TB_, Xtt_, MS_] := Module[{handle, spec},
         },
         fsModelParameters -> {
             TanBeta -> TB,
-            LambdaLoopOrder -> 2,
             MEWSB -> 173.34,
             MSUSY -> MS,
             M1Input -> MS,
@@ -82,7 +81,13 @@ CalcMh[TB_, Xtt_, MS_] := Module[{handle, spec},
             msu2 -> MS^2 IdentityMatrix[3],
             msd2 -> MS^2 IdentityMatrix[3],
             msl2 -> MS^2 IdentityMatrix[3],
-            mse2 -> MS^2 IdentityMatrix[3]
+            mse2 -> MS^2 IdentityMatrix[3],
+            LambdaLoopOrder -> 2,
+            TwoLoopAtAs -> 1,
+            TwoLoopAbAs -> 1,
+            TwoLoopAtAb -> 1,
+            TwoLoopAtauAtau -> 1,
+            TwoLoopAtAt -> 1
         }
     ];
     spec = HSSUSY /. FSHSSUSYCalculateSpectrum[handle];

meta/SM/HSSUSY_corrections.m

@@ -110,7 +110,7 @@
    We express the correction in terms of the bottom yukawa of the MSSM
    to resum tanb enhanced corrections
 *)
-lambda1LbottomMSSM = With[{
+lambda1Lbottom = With[{
     Nc = 3,
     mQ3  = Sqrt[msq2[3,3]],
     mQ32 = msq2[3,3],