adding preliminary expressions for matching scale variation in HSSUSY

meta/SM/HSSUSY_scale_variation.m

@@ -0,0 +1,37 @@
+deltaLambdaQmatch2L = With[{
+    k = 1/(4 Pi)^2
+    },
+    k^2 Get["meta/SM/HSSUSY_scale_variation_2L.m"] /. {
+        Derivative[1][TCF[1]] -> TCD1F[1],
+        Derivative[1][TCF[2]] -> TCD1F[2],
+        Derivative[1][TCF[3]] -> TCD1F[3],
+        Derivative[1][TCF[4]] -> TCD1F[4],
+        Derivative[1][TCF[5]] -> TCD1F[5],
+        Derivative[1][TCf[1]] -> TCD1f[1],
+        Derivative[1][TCf[2]] -> TCD1f[2],
+        Derivative[1][TCf[3]] -> TCD1f[3],
+        Derivative[1][TCf[4]] -> TCD1f[4],
+        Derivative[1][TCf0]   -> TCD1f0,
+        Derivative[1][TCg0]   -> TCD1g0,
+        Derivative[1,0][TCf[i_]] :> TCD10f[i],
+        Derivative[0,1][TCf[i_]] :> TCD01f[i],
+        ml2[a__] :> msl2[a],
+        me2[a__] :> mse2[a],
+        mq2[a__] :> msq2[a],
+        mu2[a__] :> msu2[a],
+        md2[a__] :> msd2[a],
+        MS -> MSUSY,
+        Q -> Qmatch,
+        gt -> Yu[3,3],
+        gb -> Yd[3,3],
+        gtau -> Ye[3,3],
+        tb -> TanBeta,
+        mA -> mAInput,
+        M1 -> M1Input,
+        M2 -> M2Input,
+        M3 -> M3Input,
+        \[Mu] -> MuInput
+    }
+];
+
+deltaLambdaQmatch3L = 0;

meta/module.mk

@@ -49,7 +49,10 @@ META_SM_SRC  := \
 		$(DIR)/SM/beta_gt.m \
 		$(DIR)/SM/beta_lambda.m \
 		$(DIR)/SM/beta_m2.m \
-		$(DIR)/SM/HSSUSY_corrections.m
+		$(DIR)/SM/HSSUSY_corrections.m \
+		$(DIR)/SM/HSSUSY_scale_variation.m \
+		$(DIR)/SM/HSSUSY_scale_variation_2L.m \
+		$(DIR)/SM/HSSUSY_scale_variation_3L.m
 
 META_THDM_SRC:= \
 		$(DIR)/THDM/Thresholds_1L_full.m

model_files/HSSUSY/FlexibleSUSY.m.in

@@ -24,7 +24,8 @@ EXTPAR = {
     {103, TwoLoopAtAb},
     {104, TwoLoopAtauAtau},
     {105, TwoLoopAtAt},
-    {200, DeltaEFT}
+    {200, DeltaEFT},
+    {201, Qmatch}
 };
 
 FSAuxiliaryParameterInfo = {
@@ -114,6 +115,13 @@ HighScaleInput = {
                       ), 0]
                     + If[TwoLoopAtauAtau >= 1, lambda2LHSSTau, 0]
                 )
+                + If[Abs[Qmatch] > 1,
+                     Which[LambdaLoopOrder < 2, deltaLambdaQmatch2L,
+                           LambdaLoopOrder < 3, deltaLambdaQmatch3L,
+                           True, 0
+                          ],
+                     0
+                    ]
     }
 };
 
@@ -169,3 +177,4 @@ SARAH`ParameterDefinitions = Append[
 ];
 
 Get[FileNameJoin[{Directory[], "meta", "SM", "HSSUSY_corrections.m"}]];
+Get[FileNameJoin[{Directory[], "meta", "SM", "HSSUSY_scale_variation.m"}]];

model_files/HSSUSY/LesHouches.in.HSSUSY

@@ -69,6 +69,7 @@ Block EXTPAR                 # Input parameters
   104   1                    # TwoLoopAtauAtau
   105   1                    # TwoLoopAtAt
   200   0                    # DeltaEFT
+  201   0                    # Qmatch
 Block MSQ2IN
   1  1     4.00000000E+06   # mq2(1,1)
   2  2     4.00000000E+06   # mq2(2,2)

src/threshold_loop_functions.hpp

@@ -33,6 +33,12 @@
 #define TCW threshold_loop_functions::W
 #define TDelta threshold_loop_functions::delta_xyz
 #define TPhi threshold_loop_functions::phi_xyz
+#define TCD1F(n) threshold_loop_functions::D1F ## n
+#define TCD1f(n) threshold_loop_functions::D1f ## n
+#define TCD10f(n) threshold_loop_functions::D10f ## n
+#define TCD01f(n) threshold_loop_functions::D01f ## n
+#define TCD1f0 threshold_loop_functions::D1f
+#define TCD1g0 threshold_loop_functions::D1g
 
 namespace flexiblesusy {
 namespace threshold_loop_functions {
@@ -61,6 +67,28 @@ double f6(double, double);
 double f7(double, double);
 double f8(double, double);
 
+double D1F1(double);
+double D1F2(double);
+double D1F3(double);
+double D1F4(double);
+double D1F5(double);
+double D1F6(double);
+double D1F7(double);
+double D1f(double);
+double D1g(double);
+double D1f1(double);
+double D1f2(double);
+double D1f3(double);
+double D1f4(double);
+double D10f5(double, double);
+double D01f5(double, double);
+double D10f6(double, double);
+double D01f6(double, double);
+double D10f7(double, double);
+double D01f7(double, double);
+double D10f8(double, double);
+double D01f8(double, double);
+
 /// \f$I_{abc}(a,b,c)\f$ (arguments are interpreted as unsquared)
 double Iabc(double, double, double);