Himalaya
H9q2.cpp
Go to the documentation of this file.
1 // ====================================================================
2 // This file is part of Himalaya.
3 //
4 // Himalaya is licenced under the GNU General Public License (GNU GPL)
5 // version 3.
6 // ====================================================================
7 
8 #include "himalaya/hierarchies/H9q2.hpp"
9 #include "himalaya/hierarchies/HierarchyFlags.hpp"
10 #include "himalaya/misc/Constants.hpp"
11 #include "himalaya/misc/Powers.hpp"
12 #include <cmath>
13 
14 namespace himalaya {
15 namespace hierarchies {
16 
17 /**
18  * Constructor
19  * @param expansionDepth the flagMap for the truncation of expansion variables
20  * @param Al4p a double alpha_s/4/Pi
21  * @param beta a double which is the mixing angle beta
22  * @param Dmst12 a double Mst1^2 - Mst2^2
23  * @param Dmsqst1 a double Msq1^2 - Mst1^2
24  * @param lmMt a double log((<renormalization scale> / Mt)^2)
25  * @param lmMgl a double log((<renormalization scale> / Mgl)^2)
26  * @param lmMst1 a double log((<renormalization scale> / Mst1)^2)
27  * @param Mgl a double gluino mass
28  * @param Mt a double top/bottom quark mass
29  * @param Mst1 a double stop 1 mass
30  * @param Mst2 a double stop 2 mass
31  * @param Msq a double average squark mass w/o the stop quark
32  * @param MuSUSY a double mu parameter
33  * @param s2t a double 2 times the sine of the stop/sbottom quark mixing angle
34  * @param mdrFlag an int 0 for DR and 1 for MDR scheme
35  * @param oneLoopFlag an int flag to consider the one-loop expansion terms
36  * @param twoLoopFlag an int flag to consider the two-loop expansion terms
37  * @param threeLoopFlag an int flag to consider the three-loop expansion terms
38  */
39 H9q2::H9q2(const ExpansionFlags_t& expansionDepth, double Al4p, double beta, double Dmst12, double Dmsqst1,
40  double lmMt, double lmMgl, double lmMst1,
41  double Mgl, double Mt, double Mst1, double Mst2, double Msq, double MuSUSY,
42  double s2t,
43  int mdrFlag, int oneLoopFlag, int twoLoopFlag, int threeLoopFlag){
44  // abbrev for tan(beta) and sin(beta)
45  Tbeta = tan(beta);
46  Sbeta = sin(beta);
47  this -> Dmst12 = Dmst12;
48  this -> Dmsqst1 = Dmsqst1;
49  this -> lmMt = lmMt;
50  this -> lmMgl = lmMgl;
51  this -> lmMst1 = lmMst1;
52  this -> Mgl = Mgl;
53  this -> Mt = Mt;
54  this -> Mst1 = Mst1;
55  this -> Mst2 = Mst2;
56  this -> Msq = Msq;
57  this -> MuSUSY = MuSUSY;
58  this -> s2t = s2t;
59  this -> oneLoopFlag = oneLoopFlag;
60  this -> twoLoopFlag = twoLoopFlag;
61  this -> threeLoopFlag = threeLoopFlag;
62  this -> Al4p = Al4p;
63  // mdr flags, indicates if one wants to shift the dr stop mass to the mdr stop mass
64  shiftst1 = mdrFlag;
65  shiftst2 = mdrFlag;
66  shiftst3 = mdrFlag;
67  // expansion flags
68  x = expansionDepth.at(ExpansionDepth::threeLoop);
69  xDmst12 = expansionDepth.at(ExpansionDepth::Dmglst1);
70  xDmsqst1 = expansionDepth.at(ExpansionDepth::Dmsqst1);
71  xMgl = expansionDepth.at(ExpansionDepth::Mgl);
72 }
73 
74 /**
75  * @return The diagonal (1, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H9q2'
76  */
77 double H9q2::getS1() const {
78  return -(pow2(MuSUSY)*((Al4p*(10*xMgl*pow4(Mgl)*(2*xDmst12*pow2(Mt)*(216*(1 - 3*
79  lmMgl + 3*lmMst1)*twoLoopFlag*pow2(Mst1)*pow2(s2t) + Al4p*
80  threeLoopFlag*(288*(262 + 335*lmMst1 - 96*lmMt - 72*lmMst1*lmMt +
81  lmMgl*(-199 - 120*lmMst1 + 72*lmMt) + 120*pow2(lmMst1))*pow2(Mt) + (
82  41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl)
83  + 540*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))*pow3(Dmst12) + pow2(Mst2)*(
84  216*(-1 + 3*lmMgl - 3*lmMst1)*twoLoopFlag*pow2(Dmst12)*pow2(Mst1)*pow2(
85  Mt)*pow2(s2t) + Al4p*threeLoopFlag*(-(pow2(Dmst12)*pow2(Mt)*(576*(262 +
86  335*lmMst1 - 96*lmMt - 72*lmMst1*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*
87  lmMt) + 120*pow2(lmMst1))*pow2(Mt) + (41653 + 11784*lmMst1 + 48*lmMgl*(
88  -313 + 180*lmMst1) - 12636*pow2(lmMgl) + 540*pow2(lmMst1))*pow2(Mst1)*
89  pow2(s2t))) + 576*pow2(Mst2)*(Dmst12*(262 + 335*lmMst1 - 96*lmMt - 72*
90  lmMst1*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt) + 120*pow2(lmMst1)) +
91  2*(59 + 85*lmMst1 - 24*lmMt - 24*lmMst1*lmMt + lmMgl*(-53 - 40*lmMst1 +
92  24*lmMt) + 40*pow2(lmMst1))*pow2(Mst2))*pow4(Mt)))) + (9*Al4p*
93  threeLoopFlag*(2*pow2(Mst1)*pow2(Mst2)*pow2(Mt)*(5*Dmst12*s2t*(-4*s2t*
94  pow2(Msq)*pow2(Mst1)*(15*Dmsqst1*(48*Dmst12*pow2(Mgl) + pow2(Mst1)*(
95  Dmst12*(3 - 6*shiftst1 - 6*shiftst2 + 4*lmMst1*(-2 + shiftst1 +
96  shiftst2)) + 4*(3 - 2*lmMst1)*(shiftst1 - shiftst2)*pow2(Mst2))) +
97  pow2(Msq)*(2*Dmst12*(200 + 18*lmMgl + 265*lmMst1 + 91*pow2(lmMst1))*
98  pow2(Mgl) + pow2(Mst1)*(-(Dmst12*(32 + 30*(shiftst1 + shiftst2) + 2*
99  lmMst1*(115 - 30*(shiftst1 + shiftst2) - 9*shiftst3) + 21*shiftst3 +
100  82*pow2(lmMst1))) + 6*(1 - 2*lmMst1)*(10*shiftst1 - 10*shiftst2 +
101  shiftst3)*pow2(Mst2)))) - 15*xDmsqst1*pow2(Dmsqst1)*(32*pow2(Mst2)*(20*
102  Mgl*Mt*pow2(Mst1) + 8*(3 - lmMgl + lmMst1)*Mt*pow3(Mgl) - (-2 + lmMst1)
103  *s2t*(shiftst1 - shiftst2)*pow4(Mst1)) + Dmst12*(192*Mgl*Mt*pow2(Mst1)
104  + 192*s2t*pow2(Mgl)*pow2(Mst1) + 256*(-3 + lmMgl - lmMst1)*Mt*pow3(Mgl)
105  + s2t*(33 + 16*lmMst1*(-2 + shiftst1 + shiftst2) - 32*(shiftst1 +
106  shiftst2))*pow4(Mst1)))) + 32*Mgl*Mt*pow2(Msq)*(5*Dmst12*s2t*(-30*
107  Dmsqst1*(Dmst12*(4*(-3 + lmMgl - lmMst1)*pow2(Mgl) + 7*pow2(Mst1)) + 2*
108  ((6 - 2*lmMgl + 2*lmMst1)*pow2(Mgl) + 3*pow2(Mst1))*pow2(Mst2)) + pow2(
109  Msq)*(Dmst12*(2*(-249 + 164*lmMgl - 176*lmMst1 - 20*lmMgl*lmMst1 + 27*
110  pow2(lmMgl) + pow2(lmMst1))*pow2(Mgl) - (149 + 41*lmMst1)*pow2(Mst1)) +
111  (4*(-11 + 154*lmMgl - 148*lmMst1 - 20*lmMgl*lmMst1 + 27*pow2(lmMgl) +
112  pow2(lmMst1))*pow2(Mgl) - (215 + 118*lmMst1 + 41*pow2(lmMst1))*pow2(
113  Mst1))*pow2(Mst2))) + Mgl*Mt*pow2(Msq)*(pow2(Dmst12)*(163 + 15*lmMgl -
114  27*lmMst1 + 20*lmMt - 8*pow2(lmMst1)) + 20*Dmst12*(16 - 9*lmMgl + 41*
115  lmMst1 - 12*lmMt + 4*pow2(lmMst1))*pow2(Mst2) - 10*(5 + 18*lmMgl - 74*
116  lmMst1 + 24*lmMt - 8*pow2(lmMst1))*pow4(Mst2)))) - 16*z2*pow2(Mt)*(
117  xDmst12*pow2(Mst1)*pow3(Dmst12)*(pow2(Msq)*(4*Mgl*Mt*(498*Mgl*Mt*pow2(
118  Msq) - 40*s2t*pow2(Mgl)*(15*Dmsqst1 + (-94 + 45*lmMgl - 45*lmMst1)*
119  pow2(Msq)) + 3*s2t*(500*Dmsqst1 - 79*pow2(Msq))*pow2(Mst1)) + pow2(
120  Mst1)*(30*Dmsqst1*(30*pow2(Mgl) - (3*shiftst1 + 2*shiftst2)*pow2(Mst1))
121  - 2*pow2(Msq)*(455*pow2(Mgl) + (10 + 45*shiftst1 + 30*shiftst2 + 27*
122  shiftst3)*pow2(Mst1)))*pow2(s2t)) - 30*s2t*xDmsqst1*pow2(Dmsqst1)*(-
123  140*Mgl*Mt*pow2(Mst1) - 30*s2t*pow2(Mgl)*pow2(Mst1) + 80*Mt*pow3(Mgl) +
124  s2t*(-15 + 3*shiftst1 + 2*shiftst2)*pow4(Mst1))) + pow2(Mst2)*(4*Mt*
125  pow2(Msq)*pow2(Mst1)*(Mt*pow2(Mgl)*pow2(Msq)*(11*pow2(Dmst12) - 520*
126  pow2(Mst2)*(Dmst12 + pow2(Mst2))) + 5*Dmst12*s2t*(Dmst12*(-30*Mgl*(4*
127  Dmsqst1 + pow2(Msq))*pow2(Mst1) + (120*Dmsqst1 + (-263 + 90*lmMgl - 90*
128  lmMst1)*pow2(Msq))*pow3(Mgl)) - 2*pow2(Mst2)*(Mgl*(30*Dmsqst1 - 17*
129  pow2(Msq))*pow2(Mst1) + (60*Dmsqst1 + (113 - 90*lmMgl + 90*lmMst1)*
130  pow2(Msq))*pow3(Mgl)))) + Dmst12*(-5*pow2(Msq)*(30*Dmsqst1*(3*Dmst12*
131  pow2(Mgl) + pow2(Mst1)*(-(Dmst12*(-3 + shiftst1 + shiftst2)) + 2*(
132  shiftst1 - shiftst2)*pow2(Mst2))) - pow2(Msq)*(274*Dmst12*pow2(Mgl) +
133  pow2(Mst1)*(-41*Dmst12 + 30*Dmst12*(shiftst1 + shiftst2) + 9*Dmst12*
134  shiftst3 - 6*(10*shiftst1 - 10*shiftst2 + shiftst3)*pow2(Mst2))))*pow2(
135  s2t)*pow4(Mst1) - 75*s2t*xDmsqst1*pow2(Dmsqst1)*pow2(Mst1)*(Dmst12*(8*
136  Mgl*Mt*pow2(Mst1) + 6*s2t*pow2(Mgl)*pow2(Mst1) - 32*Mt*pow3(Mgl) + s2t*
137  (9 - 2*shiftst1 - 2*shiftst2)*pow4(Mst1)) + 4*pow2(Mst2)*(10*Mgl*Mt*
138  pow2(Mst1) + 8*Mt*pow3(Mgl) + s2t*(shiftst1 - shiftst2)*pow4(Mst1)))))
139  + 20*xMgl*pow4(Mgl)*pow4(Msq)*(-6*xDmst12*(12*(17 - 4*lmMgl + 4*lmMst1)
140  *pow2(Mt) + (49 - 17*lmMgl + 17*lmMst1)*pow2(Mst1)*pow2(s2t))*pow3(
141  Dmst12) + pow2(Mst2)*(3*(49 - 17*lmMgl + 17*lmMst1)*pow2(Dmst12)*pow2(
142  Mst1)*pow2(s2t) + 8*pow2(Mt)*(9*Dmst12*(17 - 4*lmMgl + 4*lmMst1)*(
143  Dmst12 - pow2(Mst2)) + 2*(-47 + 12*lmMgl - 12*lmMst1)*pow4(Mst2))))) +
144  3*z3*(3*pow2(Dmst12)*pow2(Mst1)*pow2(Mt)*(-32*Dmst12*Mgl*Mt*s2t*
145  xDmst12*pow2(Mst1)*pow4(Msq) + 40*pow2(Mgl)*(-14*Dmst12*xDmst12*pow2(
146  Mt) + 15*pow2(Mst1)*(-(Dmst12*xDmst12) + 2*pow2(Mst2))*pow2(s2t))*pow4(
147  Msq) + 34560*Dmst12*Mt*s2t*xDmst12*pow3(Mgl)*pow4(Msq) - 5*pow2(s2t)*(
148  35*Dmsqst1*pow2(Msq)*(Dmst12*xDmst12 + 4*pow2(Mst2)) + 35*xDmsqst1*
149  pow2(Dmsqst1)*(-5*Dmst12*xDmst12 + 7*pow2(Mst2)) + 2*(11*Dmst12*xDmst12
150  + 24*pow2(Mst2))*pow4(Msq))*pow4(Mst1)) + 40*Mgl*pow4(Msq)*(pow2(Mst1)*
151  pow2(Mst2)*(56*Dmst12*s2t*pow2(Mst1)*pow2(Mst2) - Mgl*(648*Dmst12*Mgl*
152  s2t*(Dmst12 + 2*pow2(Mst2)) + 7*Mt*(pow2(Dmst12) - 8*pow2(Mst2)*(Dmst12
153  + pow2(Mst2)))))*pow3(Mt) + xMgl*pow3(Mgl)*(-6*xDmst12*pow3(Dmst12)*(
154  259*pow2(Mst1)*pow2(Mt)*pow2(s2t) + 484*pow4(Mt)) + pow2(Mst2)*(777*
155  pow2(Dmst12)*pow2(Mst1)*pow2(Mt)*pow2(s2t) + 968*(3*pow2(Dmst12) - 3*
156  Dmst12*pow2(Mst2) - 2*pow4(Mst2))*pow4(Mt)))))))/pow4(Msq)))/pow6(Mst1)
157  + pow2(Mt)*(-3*xDmst12*pow3(Dmst12)*(135*oneLoopFlag*pow2(s2t) + (2*
158  threeLoopFlag*pow2(Al4p)*(16*Mgl*Mt*(12*Mgl*(20*Mgl*s2t*(-130 + lmMgl*(
159  159 - 20*lmMst1) - 162*lmMst1 + 27*pow2(lmMgl) + pow2(lmMst1)) + Mt*(
160  323 - 75*lmMgl + 383*lmMst1 - 100*lmMt + 32*pow2(lmMst1))) + (3600*
161  Dmsqst1*(3 - lmMgl + lmMst1)*s2t*pow2(Mgl))/pow2(Msq) + s2t*(365 -
162  1049*lmMst1 - 123*pow2(lmMst1) - (18000*Dmsqst1)/pow2(Msq))*pow2(Mst1))
163  - (3*pow2(Mst1)*(15*Dmsqst1*(1920*pow2(Mgl) + (5 - 72*shiftst1 - 48*
164  shiftst2 + 16*lmMst1*(-5 + 3*shiftst1 + 2*shiftst2))*pow2(Mst1)) + 2*
165  pow2(Msq)*(10*(-425 + 18*lmMgl + 83*lmMst1 + 91*pow2(lmMst1))*pow2(Mgl)
166  + (605 - 180*shiftst1 - 120*shiftst2 - 348*shiftst3 + 8*lmMst1*(-110 +
167  45*shiftst1 + 30*shiftst2 + 27*shiftst3) - 820*pow2(lmMst1))*pow2(Mst1)
168  ))*pow2(s2t))/pow2(Msq) - (45*s2t*xDmsqst1*pow2(Dmsqst1)*(5120*Mgl*Mt*
169  pow2(Mst1) + 1920*s2t*pow2(Mgl)*pow2(Mst1) + 1280*(-3 + lmMgl - lmMst1)
170  *Mt*pow3(Mgl) + s2t*(215 - 72*shiftst1 - 208*shiftst2 + 16*lmMst1*(-5 +
171  3*shiftst1 + 2*shiftst2))*pow4(Mst1)))/pow4(Msq)))/pow4(Mst1) - (72*
172  Al4p*s2t*twoLoopFlag*(160*(1 - lmMgl + lmMst1)*Mt*pow2(Mst1)*pow3(Mgl)
173  + 2*(9 + 2*lmMst1)*Mgl*Mt*pow4(Mst1) + 5*(1 + 2*lmMst1)*s2t*pow2(Mgl)*
174  pow4(Mst1) + 40*(-1 + 3*lmMgl - 3*lmMst1)*Mt*x*pow5(Mgl) - 20*lmMst1*
175  s2t*pow6(Mst1)))/pow6(Mst1)) + 135*Dmst12*s2t*pow2(Mst2)*(3*Dmst12*
176  oneLoopFlag*s2t - (16*Al4p*twoLoopFlag*(4*Mgl*Mt*pow2(Mst2)*((1 - 2*
177  lmMgl + 2*lmMst1)*pow2(Mgl)*pow2(Mst1) + (1 - 3*lmMgl + 3*lmMst1)*x*
178  pow4(Mgl) + (2 + lmMst1)*pow4(Mst1)) + Dmst12*(2*(3 - 2*lmMgl + 2*
179  lmMst1)*Mt*pow2(Mst1)*pow3(Mgl) + 2*Mgl*Mt*pow4(Mst1) + (3 + 2*lmMst1)*
180  s2t*pow2(Mgl)*pow4(Mst1) + 4*(-1 + 3*lmMgl - 3*lmMst1)*Mt*x*pow5(Mgl) -
181  (1 + 2*lmMst1)*s2t*pow6(Mst1))))/pow6(Mst1)))))/(19440.*pow6(Mst2));
182 }
183 
184 /**
185  * @return The diagonal (2, 2) matrix element of the Higgs mass matrix as a double for the hierarchy 'H9q2'
186  */
187 double H9q2::getS2() const {
188  return (576*Al4p*z2*(1440*twoLoopFlag*pow2(Sbeta)*pow3(Mgl)*pow3(Mt)*pow4(Msq)*(
189  xDmst12*pow2(Mst1)*(3*x*pow2(Mgl)*(Mt*(MuSUSY - Mgl*Tbeta) + s2t*Tbeta*
190  pow2(Mgl)) - (Mt*MuSUSY + 4*s2t*Tbeta*x*pow2(Mgl))*pow2(Mst1))*pow3(
191  Dmst12) + (Dmst12*s2t*Tbeta*pow2(Mst1) + Mt*MuSUSY*(Dmst12 + pow2(Mst2)
192  ))*pow4(Mst1)*pow4(Mst2) + x*pow2(Mgl)*pow2(Mst2)*(Dmst12*pow2(Mst1)*(
193  Dmst12*(-3*Mt*(MuSUSY - Mgl*Tbeta) + s2t*Tbeta*(-3*pow2(Mgl) + pow2(
194  Mst1))) + (3*Mt*(MuSUSY - Mgl*Tbeta) + s2t*Tbeta*(3*pow2(Mgl) + 2*pow2(
195  Mst1)))*pow2(Mst2)) + Mt*(MuSUSY - Mgl*Tbeta)*(3*pow2(Mgl) + 2*pow2(
196  Mst1))*pow4(Mst2))) - 5*pow2(Mst1)*(Mt*xMgl*pow4(Mgl)*(pow2(Mst1)*pow2(
197  Sbeta)*(-9*Mt*xDmst12*pow3(Dmst12)*(400*Al4p*threeLoopFlag*xDmsqst1*
198  pow2(Dmsqst1)*(Mt*MuSUSY*s2t + 6*(4 + lmMst1 - lmMt)*Tbeta*pow2(Mt) -
199  Tbeta*pow2(Mst1)*pow2(s2t)) + 64*Tbeta*twoLoopFlag*pow2(Mt)*pow4(Msq))
200  + 144*Tbeta*twoLoopFlag*pow2(Mst2)*pow3(Mt)*pow4(Msq)*(pow2(Dmst12) +
201  2*Dmst12*pow2(Mst2) + 2*pow4(Mst2))) - Al4p*threeLoopFlag*(xDmst12*
202  pow2(Msq)*pow3(Dmst12)*(3*Mt*MuSUSY*s2t*pow2(Mst1)*(8*(-49 + 17*lmMgl -
203  17*lmMst1)*MuSUSY*s2t*Tbeta*pow2(Msq) + 16*Mt*(75*Dmsqst1 + 2*(170 -
204  27*lmMgl + 27*lmMst1)*pow2(Msq))*pow2(Sbeta)) + Tbeta*pow2(Sbeta)*(-3*
205  Mt*pow2(Mst1)*(4*(300*Dmsqst1 + (-242 + 27*lmMgl - 27*lmMst1)*pow2(Msq)
206  )*pow2(Mst1) + 8*(-49 + 17*lmMgl - 17*lmMst1)*pow2(Msq)*pow2(MuSUSY))*
207  pow2(s2t) + (21600*Dmsqst1*(4 + lmMst1 - lmMt)*pow2(Mst1) - 16*pow2(
208  Msq)*((2363 + 324*lmMgl + 564*lmMst1 - 2004*lmMt)*pow2(Mst1) + 18*(-17
209  + 4*lmMgl - 4*lmMst1)*pow2(MuSUSY)))*pow3(Mt)) + 12*MuSUSY*pow2(Msq)*(
210  24*(-17 + 4*lmMgl - 4*lmMst1)*MuSUSY*Tbeta*pow3(Mt) + (49 - 17*lmMgl +
211  17*lmMst1)*pow2(Sbeta)*pow3(s2t)*pow4(Mst1))) + 1800*Mt*xDmsqst1*pow2(
212  Dmsqst1)*pow2(Mst1)*pow2(Mst2)*pow2(Sbeta)*(2*Dmst12*Mt*(MuSUSY*s2t +
213  6*(4 + lmMst1 - lmMt)*Mt*Tbeta)*pow2(Mst2) + pow2(Dmst12)*(-2*Mt*
214  MuSUSY*s2t - 12*(4 + lmMst1 - lmMt)*Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*
215  pow2(s2t)) + 2*(19 + 5*lmMst1 - 5*lmMt)*Tbeta*pow2(Mt)*pow4(Mst2)))) +
216  Al4p*threeLoopFlag*(-(xMgl*pow2(Msq)*pow2(Mst2)*pow2(Mt)*pow4(Mgl)*(3*
217  Dmst12*s2t*pow2(Mst1)*(-600*Dmsqst1*(2*Dmst12*Mt*MuSUSY - Dmst12*s2t*
218  Tbeta*pow2(Mst1) - 2*Mt*MuSUSY*pow2(Mst2))*pow2(Sbeta) + 4*pow2(Msq)*(
219  2*Dmst12*(-146 + 27*lmMgl - 27*lmMst1)*Mt*MuSUSY*pow2(Sbeta) + 4*(-194
220  + 27*lmMgl - 27*lmMst1)*Mt*MuSUSY*pow2(Mst2)*pow2(Sbeta) + Dmst12*s2t*
221  Tbeta*((-49 + 17*lmMgl - 17*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
222  2*(-194 + 27*lmMgl - 27*lmMst1)*pow2(Mst1)*pow2(Sbeta)))) + Tbeta*pow2(
223  Mt)*(32*pow2(Msq)*pow2(MuSUSY)*(9*Dmst12*(17 - 4*lmMgl + 4*lmMst1)*(
224  Dmst12 - pow2(Mst2)) + 2*(-47 + 12*lmMgl - 12*lmMst1)*pow4(Mst2)) -
225  pow2(Sbeta)*(32*pow2(Msq)*pow2(MuSUSY)*(9*Dmst12*(17 - 4*lmMgl + 4*
226  lmMst1)*(Dmst12 - pow2(Mst2)) + 2*(-47 + 12*lmMgl - 12*lmMst1)*pow4(
227  Mst2)) + pow2(Mst1)*(7200*Dmsqst1*(4 + lmMst1 - lmMt)*(3*pow2(Dmst12) -
228  3*Dmst12*pow2(Mst2) - pow4(Mst2)) - 4*pow2(Msq)*(3*(2997 + 108*lmMgl +
229  798*lmMst1 - 1214*lmMt)*pow2(Dmst12) + 2*Dmst12*(-4265 + 324*lmMgl -
230  1266*lmMst1 - 366*lmMt)*pow2(Mst2) + 2*(-5207 + 324*lmMgl - 1716*lmMst1
231  + 84*lmMt)*pow4(Mst2))))))) + 3*xDmsqst1*xDmst12*pow2(Dmsqst1)*pow2(
232  Mst1)*pow3(Dmst12)*(Mt*(10*Mgl*pow2(Mst1)*(4*s2t*pow2(Mt)*(7*Tbeta*
233  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 63*Mgl*MuSUSY*pow2(Sbeta) - 12*(13 +
234  4*lmMst1 - 4*lmMt)*Tbeta*pow2(Mst1)*pow2(Sbeta)) - 6*Mt*pow2(s2t)*(-(
235  Mgl*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta))) + 4*(MuSUSY - 3*Mgl*Tbeta)*
236  pow2(Mst1)*pow2(Sbeta)) - 12*((78 + 20*lmMst1 - 20*lmMt)*MuSUSY - 3*(81
237  + 22*lmMst1 - 22*lmMt)*Mgl*Tbeta)*pow2(Sbeta)*pow3(Mt) + pow2(Mst1)*(3*
238  Mgl*MuSUSY + 10*Tbeta*pow2(Mst1))*pow2(Sbeta)*pow3(s2t)) + 80*pow3(Mgl)
239  *(MuSUSY*pow2(Sbeta)*(-6*Mt*pow2(Mst1)*pow2(s2t) + 72*(4 + lmMst1 -
240  lmMt)*pow3(Mt)) + Tbeta*(-2*s2t*pow2(Mt)*(pow2(MuSUSY)*(-1 + pow2(
241  Sbeta)) + 3*(65 + 14*lmMst1 - 14*lmMt)*pow2(Mst1)*pow2(Sbeta)) + pow2(
242  Sbeta)*pow3(s2t)*pow4(Mst1)))) + pow4(Mst1)*(2*Tbeta*pow2(Mt)*pow2(s2t)
243  *(-((-15 + 3*shiftst1 + 2*shiftst2)*pow2(MuSUSY)*(-1 + pow2(Sbeta))) +
244  10*(9 - 8*shiftst1 + 2*shiftst2)*pow2(Mst1)*pow2(Sbeta)) + 5*pow2(
245  Sbeta)*(96*MuSUSY*s2t*pow3(Mt) + Mt*MuSUSY*(9 - 8*shiftst1 + 4*
246  shiftst2)*pow2(Mst1)*pow3(s2t) + 16*(11 + 6*lmMst1 - 6*lmMt)*Tbeta*
247  pow4(Mt) + 2*(shiftst1 - shiftst2)*Tbeta*pow4(Mst1)*pow4(s2t)))))) -
248  Al4p*threeLoopFlag*pow4(Mst1)*(75*xDmsqst1*pow2(Dmsqst1)*pow2(Mst2)*(
249  Tbeta*pow2(Mt)*(Dmst12*s2t*(8*Mgl*Mt*(4*pow2(Mgl) + 5*pow2(Mst1))*pow2(
250  Mst2)*pow2(MuSUSY) - Dmst12*(-8*Mgl*Mt*pow2(Mst1)*pow2(MuSUSY) + 6*s2t*
251  pow2(Mgl)*pow2(Mst1)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 30*pow2(Mst1)*
252  pow2(Sbeta)) + 32*Mt*pow2(MuSUSY)*pow3(Mgl) + 9*s2t*(pow2(MuSUSY)*(-1 +
253  pow2(Sbeta)) + 8*pow2(Mst1)*pow2(Sbeta))*pow4(Mst1))) - 2*pow4(Mst1)*(
254  pow2(Dmst12)*pow2(s2t)*(-((shiftst1 + shiftst2)*pow2(MuSUSY)*(-1 +
255  pow2(Sbeta))) + 6*(-3*shiftst1 + shiftst2)*pow2(Mst1)*pow2(Sbeta)) + 2*
256  Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)*pow2(s2t)*(-1 +
257  pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 - shiftst2)*pow2(Mst1)
258  *pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)*pow2(Mt)*pow2(
259  Sbeta)*pow4(Mst2)) + 4*Mt*pow2(Sbeta)*(2*Dmst12*s2t*(Mgl*(Dmst12*(4*
260  pow2(Mgl) - pow2(Mst1)) - (4*pow2(Mgl) + 5*pow2(Mst1))*pow2(Mst2))*
261  pow2(MuSUSY) + 12*pow2(Mst1)*(Mgl*((43 + 10*lmMst1 - 10*lmMt)*pow2(Mgl)
262  + (15 + 4*lmMst1 - 4*lmMt)*pow2(Mst1))*pow2(Mst2) + Dmst12*((11 + 3*
263  lmMst1 - 3*lmMt)*Mgl*pow2(Mst1) + (11 + 2*lmMst1 - 2*lmMt)*pow3(Mgl))))
264  - 9*Mt*pow2(Mst1)*(pow2(Dmst12)*(4*(9 + 2*lmMst1 - 2*lmMt)*pow2(Mgl) +
265  (5 + 2*lmMst1 - 2*lmMt)*pow2(Mst1)) + 2*Dmst12*((45 + 14*lmMst1 - 14*
266  lmMt)*pow2(Mgl) + (3 + 2*lmMst1 - 2*lmMt)*pow2(Mst1))*pow2(Mst2) + 4*((
267  13 + 4*lmMst1 - 4*lmMt)*pow2(Mgl) + (1 + lmMst1 - lmMt)*pow2(Mst1))*
268  pow4(Mst2)))) + 6*MuSUSY*pow2(Sbeta)*(pow3(Mt)*(4*Dmst12*pow2(Mst2)*(2*
269  (45 + 14*lmMst1 - 14*lmMt)*Mgl*Mt*pow2(Mst1) - 15*s2t*pow2(Mgl)*pow2(
270  Mst1) + 48*(4 + lmMst1 - lmMt)*Mt*pow3(Mgl) + 2*s2t*(-3 + shiftst2)*
271  pow4(Mst1)) + 8*(2*(13 + 4*lmMst1 - 4*lmMt)*Mgl*Mt*pow2(Mst1) + 4*(19 +
272  5*lmMst1 - 5*lmMt)*Mt*pow3(Mgl) + s2t*(-shiftst1 + shiftst2)*pow4(Mst1)
273  )*pow4(Mst2)) - Mt*pow2(Dmst12)*(4*Mgl*((6 + 4*lmMst1 - 4*lmMt)*Mt + 3*
274  Mgl*s2t)*pow2(Mst1)*pow2(Mt) + 8*pow3(Mgl)*(-(Mt*pow2(Mst1)*pow2(s2t))
275  + 24*(4 + lmMst1 - lmMt)*pow3(Mt)) + 2*Mt*s2t*(6*Mt - 5*Mgl*s2t)*pow4(
276  Mst1) - (shiftst1 - shiftst2)*pow3(s2t)*pow6(Mst1)))) - xDmst12*pow2(
277  Msq)*pow3(Dmst12)*(Tbeta*(4*Mgl*(498*Mgl*Mt*pow2(Msq) - 40*s2t*pow2(
278  Mgl)*(15*Dmsqst1 + (-94 + 45*lmMgl - 45*lmMst1)*pow2(Msq)) + 3*s2t*(
279  500*Dmsqst1 - 79*pow2(Msq))*pow2(Mst1))*pow2(MuSUSY)*pow3(Mt) + pow2(
280  s2t)*(5*pow2(Mst1)*pow2(Mt)*(2*(90*Dmsqst1*pow2(Mgl) - pow2(Msq)*(91*
281  pow2(Mgl) + 2*pow2(Mst1)))*pow2(MuSUSY) + 6*shiftst1*(Dmsqst1 + pow2(
282  Msq))*pow2(Mst1)*(80*pow2(Mst1) + 3*pow2(MuSUSY))*pow2(Sbeta)) - 30*(
283  Dmsqst1 + pow2(Msq))*pow4(Mst1)*(pow2(Mt)*(20*shiftst2*pow2(Mst1)*pow2(
284  Sbeta) + pow2(MuSUSY)*(3*shiftst1 + 2*shiftst2 - 2*shiftst2*pow2(Sbeta)
285  )) + 5*(shiftst1 - shiftst2)*pow2(s2t)*pow2(Sbeta)*pow4(Mst1))) + Mt*
286  pow2(Sbeta)*(-20*(Mgl*(30*Dmsqst1 - 17*pow2(Msq))*pow2(Mst1) + (60*
287  Dmsqst1 + (113 - 90*lmMgl + 90*lmMst1)*pow2(Msq))*pow3(Mgl))*pow3(s2t)*
288  pow4(Mst1) + 10*Mt*pow2(Mst1)*pow2(s2t)*(pow2(Msq)*(2*pow2(Mst1)*pow2(
289  MuSUSY) + pow2(Mgl)*(1332*pow2(Mst1) + 91*pow2(MuSUSY)) - 90*pow4(Mst1)
290  ) - 90*Dmsqst1*(pow2(Mgl)*(-6*pow2(Mst1) + pow2(MuSUSY)) + 2*pow4(Mst1)
291  )) - 4*Mgl*s2t*pow2(Mt)*(300*Dmsqst1*(5*pow2(Mst1)*pow2(MuSUSY) - 2*
292  pow2(Mgl)*(30*(9 + 2*lmMst1 - 2*lmMt)*pow2(Mst1) + pow2(MuSUSY)) + 6*
293  pow4(Mst1)) + pow2(Msq)*(-237*pow2(Mst1)*pow2(MuSUSY) + 20*pow2(Mgl)*(
294  21*(56 + 11*lmMst1 - 11*lmMt)*pow2(Mst1) + 2*(94 - 45*lmMgl + 45*
295  lmMst1)*pow2(MuSUSY)) + 60*(32 + 5*lmMst1 - 5*lmMt)*pow4(Mst1))) -
296  pow3(Mt)*(1200*Dmsqst1*(63*(7 + 2*lmMst1 - 2*lmMt)*pow2(Mgl)*pow2(Mst1)
297  - pow4(Mst1)) - 8*pow2(Msq)*(pow2(Mgl)*(5*(929 + 210*lmMst1 - 210*lmMt)
298  *pow2(Mst1) - 249*pow2(MuSUSY)) + 5*(17 - 24*lmMst1 + 24*lmMt)*pow4(
299  Mst1)))) - 3*shiftst3*pow2(Msq)*pow4(Mst1)*(-2*pow2(Mt)*pow2(s2t)*(9*
300  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 145*pow2(Mst1)*pow2(Sbeta)) + 5*pow2(
301  Sbeta)*(24*pow4(Mt) + pow4(Mst1)*pow4(s2t)))) + Mt*MuSUSY*pow2(Sbeta)*(
302  pow2(Msq)*(-80*Mgl*Mt*pow2(Mst1)*((16 - 9*lmMst1 + 9*lmMt)*pow2(Mt) +
303  24*pow2(Mst1)*pow2(s2t)) + 10*Mt*(32*(380 + 147*lmMst1 - 249*lmMt)*
304  pow2(Mt) - 3*(37 + 90*lmMgl - 90*lmMst1)*pow2(Mst1)*pow2(s2t))*pow3(
305  Mgl) + 5*s2t*(-8*(-16 + 9*shiftst3)*pow2(Mt) + (-41 + 120*shiftst1 -
306  60*shiftst2 + 21*shiftst3)*pow2(Mst1)*pow2(s2t))*pow4(Mst1) + 2*pow2(
307  Mgl)*(4556*s2t*pow2(Mst1)*pow2(Mt) + 685*pow3(s2t)*pow4(Mst1))) - 150*
308  Dmsqst1*(12*Mgl*Mt*pow2(Mst1)*(-16*(7 + 2*lmMst1 - 2*lmMt)*pow2(Mt) +
309  pow2(Mst1)*pow2(s2t)) + 48*pow3(Mgl)*(-(Mt*pow2(Mst1)*pow2(s2t)) + 12*(
310  4 + lmMst1 - lmMt)*pow3(Mt)) + 3*pow2(Mgl)*(120*s2t*pow2(Mst1)*pow2(Mt)
311  + pow3(s2t)*pow4(Mst1)) + (3 - 4*shiftst1 + 2*shiftst2)*pow3(s2t)*pow6(
312  Mst1)))))) + threeLoopFlag*pow2(Al4p)*(-576*pow4(Mst1)*(900*xDmsqst1*
313  xDmst12*xMgl*pow2(Dmsqst1)*pow2(Mt)*(2*(29 - 10*lmMgl + 10*lmMst1)*Mt*
314  MuSUSY*s2t + Tbeta*(307 + 4*lmMst1*(67 - 5*lmMt) - 10*lmMgl*(3 + 2*
315  lmMst1 - 2*lmMt) - 238*lmMt + 20*pow2(lmMst1))*pow2(Mt) + 2*(-29 + 10*
316  lmMgl - 10*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t))*pow2(Sbeta)*pow3(Dmst12)
317  *pow4(Mgl) - z2*pow2(Msq)*pow2(Mst2)*(5*MuSUSY*pow2(Sbeta)*(12*Dmst12*
318  s2t*pow2(Mst1)*pow2(Mt)*(6*Mt*(10*Dmsqst1*(Dmst12*(6*pow2(Mgl) + pow2(
319  Mst1)) + (3*pow2(Mgl) + 2*pow2(Mst1))*pow2(Mst2)) + pow2(Msq)*(5*
320  Dmst12*(3*pow2(Mgl) + pow2(Mst1)) + (-59*pow2(Mgl) + 10*pow2(Mst1))*
321  pow2(Mst2))) - Dmst12*s2t*(Mgl*(30*Dmsqst1 - 17*pow2(Msq))*pow2(Mst1) +
322  (60*Dmsqst1 + (113 - 90*lmMgl + 90*lmMst1)*pow2(Msq))*pow3(Mgl))) +
323  pow4(Mst1)*(90*(Dmsqst1 + pow2(Msq))*(-8*Dmst12*s2t*shiftst2*pow2(Mst2)
324  *pow3(Mt) + Mt*(-shiftst1 + shiftst2)*pow2(Dmst12)*pow2(Mst1)*pow3(s2t)
325  + 8*s2t*(shiftst1 - shiftst2)*pow3(Mt)*pow4(Mst2)) + 9*shiftst3*pow2(
326  Msq)*(-(Mt*pow2(Dmst12)*pow2(Mst1)*pow3(s2t)) + 8*s2t*pow3(Mt)*(pow2(
327  Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)))) + (-16*Dmst12*Mgl*(180*
328  Dmsqst1*(6*(4 + lmMst1 - lmMt)*pow2(Mgl) + (7 + 2*lmMst1 - 2*lmMt)*
329  pow2(Mst1)) - pow2(Msq)*((1627 + 234*lmMst1 + 174*lmMt)*pow2(Mgl) + (-
330  457 - 87*lmMst1 + 87*lmMt)*pow2(Mst1)))*pow2(Mst2) + pow2(Dmst12)*(-8*
331  Mgl*(180*Dmsqst1*(7 + 2*lmMst1 - 2*lmMt) + (691 + 216*lmMst1 - 216*
332  lmMt)*pow2(Msq))*pow2(Mst1) + 24*(720*Dmsqst1*(4 + lmMst1 - lmMt) + (-
333  1049 - 274*lmMst1 + 274*lmMt)*pow2(Msq))*pow3(Mgl)) - 32*(90*Dmsqst1*(4
334  + lmMst1 - lmMt)*Mgl*(2*pow2(Mgl) + pow2(Mst1)) + Mgl*pow2(Msq)*((-1070
335  - 207*lmMst1 + 3*lmMt)*pow2(Mgl) + (100 - 3*lmMst1 + 3*lmMt)*pow2(Mst1)
336  ))*pow4(Mst2))*pow4(Mt)) + Tbeta*(pow2(Mt)*pow2(MuSUSY)*(4*pow2(Mgl)*
337  pow2(Msq)*(11*pow2(Dmst12) - 520*pow2(Mst2)*(Dmst12 + pow2(Mst2)))*
338  pow2(Mt) - 5*pow2(Dmst12)*pow2(Mst1)*(90*Dmsqst1*(pow2(Mgl) + pow2(
339  Mst1)) + pow2(Msq)*(-274*pow2(Mgl) + 41*pow2(Mst1)))*pow2(s2t) + 20*
340  Dmst12*Mt*s2t*(Dmst12*(-30*Mgl*(4*Dmsqst1 + pow2(Msq))*pow2(Mst1) + (
341  120*Dmsqst1 + (-263 + 90*lmMgl - 90*lmMst1)*pow2(Msq))*pow3(Mgl)) - 2*
342  pow2(Mst2)*(Mgl*(30*Dmsqst1 - 17*pow2(Msq))*pow2(Mst1) + (60*Dmsqst1 +
343  (113 - 90*lmMgl + 90*lmMst1)*pow2(Msq))*pow3(Mgl)))) + pow2(Mt)*pow2(
344  Sbeta)*(-5*Dmst12*Mt*s2t*(4*pow2(MuSUSY)*(Dmst12*(-30*Mgl*(4*Dmsqst1 +
345  pow2(Msq))*pow2(Mst1) + (120*Dmsqst1 + (-263 + 90*lmMgl - 90*lmMst1)*
346  pow2(Msq))*pow3(Mgl)) - 2*pow2(Mst2)*(Mgl*(30*Dmsqst1 - 17*pow2(Msq))*
347  pow2(Mst1) + (60*Dmsqst1 + (113 - 90*lmMgl + 90*lmMst1)*pow2(Msq))*
348  pow3(Mgl))) + pow2(Mst1)*(16*Mgl*(90*Dmsqst1*(9 + 2*lmMst1 - 2*lmMt)*(
349  2*pow2(Mgl) + pow2(Mst1)) + pow2(Msq)*((-2143 - 414*lmMst1 + 6*lmMt)*
350  pow2(Mgl) + (197 - 6*lmMst1 + 6*lmMt)*pow2(Mst1)))*pow2(Mst2) + Dmst12*
351  (8*Mgl*(180*Dmsqst1*(7 + 2*lmMst1 - 2*lmMt) + (607 + 186*lmMst1 - 186*
352  lmMt)*pow2(Msq))*pow2(Mst1) + 144*(40*Dmsqst1*(9 + 2*lmMst1 - 2*lmMt) +
353  (67 + 20*lmMst1 - 20*lmMt)*pow2(Msq))*pow3(Mgl)))) + 5*pow2(Dmst12)*
354  pow2(Mst1)*pow2(s2t)*(36*pow2(Mgl)*(30*Dmsqst1 - 59*pow2(Msq))*pow2(
355  Mst1) + (90*Dmsqst1*(pow2(Mgl) + pow2(Mst1)) + pow2(Msq)*(-274*pow2(
356  Mgl) + 41*pow2(Mst1)))*pow2(MuSUSY) + 360*(2*Dmsqst1 + pow2(Msq))*pow4(
357  Mst1)) + pow2(Mt)*(-4*pow2(Mgl)*pow2(Msq)*(11*pow2(Dmst12) - 520*pow2(
358  Mst2)*(Dmst12 + pow2(Mst2)))*pow2(MuSUSY) + (7200*Dmsqst1*(1 + lmMst1 -
359  lmMt) + 40*(223 + 180*lmMst1 - 180*lmMt)*pow2(Msq))*pow4(Mst1)*pow4(
360  Mst2) + 40*pow2(Mst1)*(pow2(Dmst12)*(pow2(Msq)*(4*(613 + 165*lmMst1 -
361  165*lmMt)*pow2(Mgl) + 23*(5 + 3*lmMst1 - 3*lmMt)*pow2(Mst1)) + 45*
362  Dmsqst1*(15*(7 + 2*lmMst1 - 2*lmMt)*pow2(Mgl) + (5 + 2*lmMst1 - 2*lmMt)
363  *pow2(Mst1))) + 2*Dmst12*(45*Dmsqst1*(6*(7 + 2*lmMst1 - 2*lmMt)*pow2(
364  Mgl) + (3 + 2*lmMst1 - 2*lmMt)*pow2(Mst1)) + pow2(Msq)*((679 + 150*
365  lmMst1 - 150*lmMt)*pow2(Mgl) + 45*(2 + lmMst1 - lmMt)*pow2(Mst1)))*
366  pow2(Mst2) + 4*pow2(Mgl)*(135*Dmsqst1*(4 + lmMst1 - lmMt) + (124 + 15*
367  lmMst1 - 15*lmMt)*pow2(Msq))*pow4(Mst2)))) + pow4(Mst1)*(150*(Dmsqst1 +
368  pow2(Msq))*pow2(Mt)*(-(pow2(Dmst12)*pow2(s2t)*((shiftst1 + shiftst2)*
369  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 6*(3*shiftst1 - shiftst2)*pow2(Mst1)*
370  pow2(Sbeta))) + 2*Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)
371  *pow2(s2t)*(-1 + pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 -
372  shiftst2)*pow2(Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)
373  *pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) + 15*shiftst3*pow2(Msq)*(-(Dmst12*
374  pow2(Mt)*pow2(s2t)*((3*Dmst12 - 2*pow2(Mst2))*pow2(MuSUSY)*(-1 + pow2(
375  Sbeta)) + 6*pow2(Mst1)*(7*Dmst12 - 4*pow2(Mst2))*pow2(Sbeta))) + 24*
376  pow2(Sbeta)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2))*pow4(Mt))))
377  )) + 27*Mt*z3*pow2(Mst1)*(75*Dmsqst1*pow2(Mst1)*(Dmsqst1*xDmsqst1*(-4*
378  Mt*pow2(Dmst12)*pow2(Mst2)*((3*Mt*(70*MuSUSY*s2t + 449*Mt*Tbeta)*pow2(
379  Sbeta) + 7*Tbeta*pow2(s2t)*(7*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 30*
380  pow2(Mst1)*pow2(Sbeta)))*pow4(Mst1) - 2048*Mgl*Mt*pow2(Sbeta)*(-(Mt*
381  MuSUSY*pow2(Mst1)) - 3*Mgl*Mt*Tbeta*pow2(Mst1) + 2*pow2(Mgl)*(-6*Mt*
382  MuSUSY + s2t*Tbeta*pow2(Mst1)) + 15*Mt*Tbeta*xMgl*pow3(Mgl) + 3*s2t*
383  Tbeta*pow4(Mst1))) + xDmst12*pow3(Dmst12)*((140*Mt*Tbeta*pow2(s2t)*(
384  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*pow2(Mst1)*pow2(Sbeta)) + 7184*
385  Tbeta*pow2(Sbeta)*pow3(Mt) + 7*MuSUSY*pow2(Sbeta)*(159*s2t*pow2(Mt) +
386  28*pow2(Mst1)*pow3(s2t)))*pow4(Mst1) - 4096*Mgl*pow2(Mt)*pow2(Sbeta)*(
387  10*Mt*MuSUSY*pow2(Mst1) - 33*Mgl*Mt*Tbeta*pow2(Mst1) - 4*pow2(Mgl)*(6*
388  Mt*MuSUSY - 7*s2t*Tbeta*pow2(Mst1)) + 30*Mt*Tbeta*xMgl*pow3(Mgl) + 8*
389  s2t*Tbeta*pow4(Mst1))) + 16*pow2(Mt)*pow2(Sbeta)*(1024*Dmst12*(6*Mt*
390  MuSUSY + 5*s2t*Tbeta*pow2(Mst1))*pow3(Mgl) - 768*Dmst12*Mt*Tbeta*(7*
391  pow2(Mgl)*pow2(Mst1) + 10*xMgl*pow4(Mgl)) - 3*Dmst12*(35*MuSUSY*s2t +
392  214*Mt*Tbeta)*pow4(Mst1) - 4*Mt*pow2(Mst2)*(-512*Mgl*MuSUSY*pow2(Mst1)
393  + 768*Tbeta*pow2(Mgl)*pow2(Mst1) - 1280*MuSUSY*pow3(Mgl) + 1600*Tbeta*
394  xMgl*pow4(Mgl) + 171*Tbeta*pow4(Mst1)) + 512*Dmst12*Mgl*(7*Mt*MuSUSY*
395  pow2(Mst1) + 4*s2t*Tbeta*pow4(Mst1)))*pow4(Mst2)) + pow2(Msq)*(-16*Mt*
396  pow2(Dmst12)*pow2(Mst2)*(-256*Mgl*Mt*pow2(Sbeta)*(-15*Mgl*Mt*Tbeta*
397  pow2(Mst1) - 8*pow2(Mgl)*(3*Mt*MuSUSY - 2*s2t*Tbeta*pow2(Mst1)) + 4*
398  pow2(Mst1)*(Mt*MuSUSY + s2t*Tbeta*pow2(Mst1)) + 30*Mt*Tbeta*xMgl*pow3(
399  Mgl)) + (Mt*(49*MuSUSY*s2t + 214*Mt*Tbeta)*pow2(Sbeta) + 7*Tbeta*pow2(
400  s2t)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 4*pow2(Mst1)*pow2(Sbeta)))*
401  pow4(Mst1)) - xDmst12*pow3(Dmst12)*(8192*Mgl*pow2(Mt)*pow2(Sbeta)*(8*
402  Mt*MuSUSY*pow2(Mst1) - 21*Mgl*Mt*Tbeta*pow2(Mst1) - 4*pow2(Mgl)*(3*Mt*
403  MuSUSY - 5*s2t*Tbeta*pow2(Mst1)) + 15*Mt*Tbeta*xMgl*pow3(Mgl)) - 7*(4*
404  Mt*Tbeta*pow2(s2t)*(pow2(MuSUSY) + (2*pow2(Mst1) - pow2(MuSUSY))*pow2(
405  Sbeta)) - 24*Tbeta*pow2(Sbeta)*pow3(Mt) + MuSUSY*pow2(Sbeta)*(31*s2t*
406  pow2(Mt) + 16*pow2(Mst1)*pow3(s2t)))*pow4(Mst1)) + 64*pow2(Mt)*pow2(
407  Sbeta)*(-128*Mt*pow2(Mst2)*(-2*Mgl*MuSUSY*pow2(Mst1) + 3*Tbeta*pow2(
408  Mgl)*pow2(Mst1) - 4*MuSUSY*pow3(Mgl) + 5*Tbeta*xMgl*pow4(Mgl) + Tbeta*
409  pow4(Mst1)) - Dmst12*(-512*(3*Mt*MuSUSY + s2t*Tbeta*pow2(Mst1))*pow3(
410  Mgl) + 384*Mt*Tbeta*(2*pow2(Mgl)*pow2(Mst1) + 5*xMgl*pow4(Mgl)) + (14*
411  MuSUSY*s2t + 107*Mt*Tbeta)*pow4(Mst1) - 256*Mgl*(2*Mt*MuSUSY*pow2(Mst1)
412  + s2t*Tbeta*pow4(Mst1))))*pow4(Mst2))) - 8*pow4(Msq)*(-80*Dmst12*pow2(
413  Mt)*pow4(Mst2)*(2*pow2(Mgl)*pow2(Mst1)*(7*Mt*Tbeta*pow2(MuSUSY)*(-1 +
414  pow2(Sbeta)) + 200*Mgl*Mt*MuSUSY*pow2(Sbeta) + 4*(3*MuSUSY*s2t - 199*
415  Mt*Tbeta)*pow2(Mst1)*pow2(Sbeta) - 2*Mgl*s2t*Tbeta*(81*pow2(MuSUSY)*(-1
416  + pow2(Sbeta)) + 380*pow2(Mst1)*pow2(Sbeta))) + 2*(xMgl*(-363*Mt*Tbeta*
417  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 2*(279*MuSUSY*s2t + 58*Mt*Tbeta)*
418  pow2(Mst1)*pow2(Sbeta))*pow4(Mgl) + Mgl*(496*Mt*MuSUSY*pow2(Sbeta) +
419  s2t*Tbeta*(7*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 32*pow2(Mst1)*pow2(
420  Sbeta)))*pow4(Mst1)) + 3*(31*MuSUSY*s2t - 152*Mt*Tbeta)*pow2(Sbeta)*
421  pow6(Mst1)) + 20*Mt*pow2(Dmst12)*pow2(Mst2)*(pow2(Mgl)*pow2(Mst1)*(-8*
422  Mgl*(MuSUSY*(2192*pow2(Mt) + 243*pow2(Mst1)*pow2(s2t))*pow2(Sbeta) + 3*
423  Mt*s2t*Tbeta*(-27*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 320*pow2(Mst1)*
424  pow2(Sbeta))) + Tbeta*(-6*pow2(Mst1)*pow2(s2t)*(15*pow2(MuSUSY)*(-1 +
425  pow2(Sbeta)) + 8*pow2(Mst1)*pow2(Sbeta)) + pow2(Mt)*(7*pow2(MuSUSY)*(-1
426  + pow2(Sbeta)) + 14080*pow2(Mst1)*pow2(Sbeta)))) - 4*Mgl*(928*Mt*s2t*
427  Tbeta*pow2(Mst1) + 1136*MuSUSY*pow2(Mt) - 21*MuSUSY*pow2(Mst1)*pow2(
428  s2t))*pow2(Sbeta)*pow4(Mst1) - xMgl*pow4(Mgl)*(pow2(Mst1)*(777*Tbeta*
429  pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) - 8*Mt*(-279*MuSUSY*s2t +
430  4442*Mt*Tbeta)*pow2(Sbeta)) + Tbeta*(2904*pow2(Mt)*pow2(MuSUSY)*(-1 +
431  pow2(Sbeta)) + 2232*pow2(s2t)*pow2(Sbeta)*pow4(Mst1))) + (Mt*(396*
432  MuSUSY*s2t + 1141*Mt*Tbeta)*pow2(Sbeta) - 6*Tbeta*pow2(s2t)*(-3*pow2(
433  MuSUSY)*(-1 + pow2(Sbeta)) + 31*pow2(Mst1)*pow2(Sbeta)))*pow6(Mst1)) +
434  xDmst12*pow3(Dmst12)*(-20*xMgl*pow4(Mgl)*(-2*Mt*pow2(Mst1)*(777*Tbeta*
435  pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 16*Mt*(279*MuSUSY*s2t -
436  2192*Mt*Tbeta)*pow2(Sbeta)) - 2904*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)
437  )*pow3(Mt) - 3*(259*MuSUSY*s2t + 372*Mt*Tbeta)*pow2(s2t)*pow2(Sbeta)*
438  pow4(Mst1)) + 16*Mgl*pow4(Mst1)*(-(s2t*Tbeta*pow2(Mt)*(-3*pow2(MuSUSY)*
439  (-1 + pow2(Sbeta)) + 1040*pow2(Mst1)*pow2(Sbeta))) + 5*MuSUSY*pow2(
440  Sbeta)*(-21*Mt*pow2(Mst1)*pow2(s2t) + 88*pow3(Mt)) + 35*Tbeta*pow2(
441  Sbeta)*pow3(s2t)*pow4(Mst1)) + 20*pow2(Mgl)*pow2(Mst1)*(-(Mt*pow2(Mst1)
442  *(-45*Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 14*Mt*(MuSUSY*
443  s2t - 320*Mt*Tbeta)*pow2(Sbeta))) + 42*Tbeta*pow2(MuSUSY)*(-1 + pow2(
444  Sbeta))*pow3(Mt) + 6*(15*MuSUSY*s2t + 8*Mt*Tbeta)*pow2(s2t)*pow2(Sbeta)
445  *pow4(Mst1) + 4*Mgl*(-8*s2t*Tbeta*pow2(Mt)*(81*pow2(MuSUSY)*(-1 + pow2(
446  Sbeta)) + 308*pow2(Mst1)*pow2(Sbeta)) + 3*MuSUSY*pow2(Sbeta)*(81*Mt*
447  pow2(Mst1)*pow2(s2t) + 3056*pow3(Mt)) - 162*Tbeta*pow2(Sbeta)*pow3(s2t)
448  *pow4(Mst1))) - 15*(-(Mt*Tbeta*pow2(s2t)*(11*pow2(MuSUSY)*(-1 + pow2(
449  Sbeta)) + 512*pow2(Mst1)*pow2(Sbeta))) + 500*Tbeta*pow2(Sbeta)*pow3(Mt)
450  + MuSUSY*pow2(Sbeta)*(85*s2t*pow2(Mt) + 24*pow2(Mst1)*pow3(s2t)))*pow6(
451  Mst1)) + 160*pow3(Mt)*(8*MuSUSY*pow2(Sbeta)*(95*pow2(Mst1)*pow3(Mgl) -
452  4*Mgl*pow4(Mst1)) + Tbeta*(2*xMgl*(121*pow2(MuSUSY)*(-1 + pow2(Sbeta))
453  - 658*pow2(Mst1)*pow2(Sbeta))*pow4(Mgl) + pow2(Mgl)*(-7*pow2(Mst1)*
454  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 152*pow2(Sbeta)*pow4(Mst1)) + 2*(395
455  - 12*lmMst1 + 12*lmMt)*pow2(Sbeta)*pow6(Mst1)))*pow6(Mst2)))) + pow2(
456  Mst2)*(72*threeLoopFlag*pow2(Al4p)*pow4(Mst1)*(-4*pow2(Msq)*(Tbeta*
457  pow2(Mt)*(-(pow2(MuSUSY)*(-40*Dmst12*Mgl*Mt*(-4*Mgl*Mt*(-16 + 9*lmMgl -
458  41*lmMst1 + 12*lmMt - 4*pow2(lmMst1))*pow2(Msq) + 4*s2t*pow2(Mgl)*(30*
459  Dmsqst1*(-3 + lmMgl - lmMst1) + (-11 + 154*lmMgl - 148*lmMst1 - 20*
460  lmMgl*lmMst1 + 27*pow2(lmMgl) + pow2(lmMst1))*pow2(Msq)) - s2t*(180*
461  Dmsqst1 + (215 + 118*lmMst1 + 41*pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*
462  pow2(Mst2) - pow2(Dmst12)*(-2*pow2(Mgl)*(1800*Dmsqst1*pow2(Mst1)*pow2(
463  s2t) + pow2(Msq)*(-4*(163 + 15*lmMgl - 27*lmMst1 + 20*lmMt - 8*pow2(
464  lmMst1))*pow2(Mt) + 5*(200 + 18*lmMgl + 265*lmMst1 + 91*pow2(lmMst1))*
465  pow2(Mst1)*pow2(s2t))) + Mt*s2t*(-40*Mgl*(210*Dmsqst1 + (149 + 41*
466  lmMst1)*pow2(Msq))*pow2(Mst1) + 80*(60*Dmsqst1*(3 - lmMgl + lmMst1) + (
467  -249 - 176*lmMst1 - 4*lmMgl*(-41 + 5*lmMst1) + 27*pow2(lmMgl) + pow2(
468  lmMst1))*pow2(Msq))*pow3(Mgl)) - 5*(15*Dmsqst1*(3 - 8*lmMst1) - 2*(16 +
469  115*lmMst1 + 41*pow2(lmMst1))*pow2(Msq))*pow2(s2t)*pow4(Mst1)) + 80*(5
470  + 18*lmMgl - 74*lmMst1 + 24*lmMt - 8*pow2(lmMst1))*pow2(Mgl)*pow2(Msq)*
471  pow2(Mt)*pow4(Mst2))) + pow4(Mst1)*(15*shiftst3*pow2(Msq)*(-2*Dmst12*
472  pow2(Mst2)*((-1 + 2*lmMst1)*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) +
473  12*((3 - 2*lmMst1)*pow2(Mt) + (-1 + 2*lmMst1)*pow2(Mst1)*pow2(s2t))*
474  pow2(Sbeta)) + pow2(Dmst12)*((-7 + 6*lmMst1)*pow2(MuSUSY)*pow2(s2t)*(-1
475  + pow2(Sbeta)) + 6*(-8*(-2 + lmMst1)*pow2(Mt) + (-15 + 14*lmMst1)*pow2(
476  Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(1 - 2*lmMst1)*pow2(Mt)*pow2(Sbeta)*
477  pow4(Mst2)) + 150*(Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*(
478  pow2(Dmst12)*pow2(s2t)*(-((shiftst1 + shiftst2)*pow2(MuSUSY)*(-1 +
479  pow2(Sbeta))) + 6*(-3*shiftst1 + shiftst2)*pow2(Mst1)*pow2(Sbeta)) + 2*
480  Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)*pow2(s2t)*(-1 +
481  pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 - shiftst2)*pow2(Mst1)
482  *pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)*pow2(Mt)*pow2(
483  Sbeta)*pow4(Mst2))) + pow2(Sbeta)*(5*pow2(Dmst12)*pow2(Mst1)*pow2(s2t)*
484  (15*Dmsqst1*((3 - 8*lmMst1)*pow2(Mst1)*pow2(MuSUSY) + 48*pow2(Mgl)*(9*
485  pow2(Mst1) + pow2(MuSUSY)) + (60 - 48*lmMst1)*pow4(Mst1)) + 2*pow2(Msq)
486  *(6*(525 - 18*lmMgl + 504*lmMst1 - 48*lmMt + 107*pow2(lmMst1))*pow2(
487  Mgl)*pow2(Mst1) + (200 + 18*lmMgl + 265*lmMst1 + 91*pow2(lmMst1))*pow2(
488  Mgl)*pow2(MuSUSY) - (16 + 115*lmMst1 + 41*pow2(lmMst1))*pow2(Mst1)*
489  pow2(MuSUSY) - 12*(-6 + 79*lmMst1 + 41*pow2(lmMst1))*pow4(Mst1))) - 40*
490  Dmst12*Mgl*Mt*s2t*(30*Dmsqst1*(Dmst12*(-7*pow2(Mst1)*pow2(MuSUSY) + 4*
491  pow2(Mgl)*(-3*(-24 + lmMgl - 15*lmMst1 + 14*lmMt)*pow2(Mst1) + (3 -
492  lmMgl + lmMst1)*pow2(MuSUSY)) + 6*(5 + 6*lmMst1 - 6*lmMt)*pow4(Mst1)) +
493  2*pow2(Mst2)*(-3*pow2(Mst1)*pow2(MuSUSY) - 2*pow2(Mgl)*((-57 - 39*
494  lmMst1 + lmMgl*(9 + 6*lmMst1 - 6*lmMt) + 30*lmMt + 6*lmMst1*lmMt - 6*
495  pow2(lmMst1))*pow2(Mst1) + (3 - lmMgl + lmMst1)*pow2(MuSUSY)) + 24*(1 +
496  lmMst1 - lmMt)*pow4(Mst1))) + pow2(Msq)*(pow2(Mst2)*(-((215 + 118*
497  lmMst1 + 41*pow2(lmMst1))*pow2(Mst1)*pow2(MuSUSY)) + 2*pow2(Mgl)*(2*(-
498  11 + lmMgl*(154 - 20*lmMst1) - 148*lmMst1 + 27*pow2(lmMgl) + pow2(
499  lmMst1))*pow2(MuSUSY) + pow2(Mst1)*(2086 - 996*lmMt - 81*(9 + 4*lmMst1
500  - 6*lmMt)*pow2(lmMgl) + (663 + 366*lmMt)*pow2(lmMst1) - 144*pow2(lmMt)
501  - 6*lmMst1*(-521 + 74*lmMt + 24*pow2(lmMt)) + 6*lmMgl*(-238 + lmMst1*(
502  73 - 142*lmMt) + 52*lmMt + 91*pow2(lmMst1) + 24*pow2(lmMt)) - 54*pow3(
503  lmMgl) - 168*pow3(lmMst1))) - 4*(-407 + 276*lmMt + 3*lmMst1*(-65 + 2*
504  lmMt) + 45*pow2(lmMst1) + 72*pow2(lmMt))*pow4(Mst1)) + Dmst12*(-((149 +
505  41*lmMst1)*pow2(Mst1)*pow2(MuSUSY)) + 2*pow2(Mgl)*(3*(1091 + 611*lmMst1
506  + lmMgl*(-17 + 122*lmMst1 - 82*lmMt) + (-374 + 62*lmMst1)*lmMt - 54*
507  pow2(lmMgl) - 24*(pow2(lmMst1) + pow2(lmMt)))*pow2(Mst1) + (-249 - 176*
508  lmMst1 - 4*lmMgl*(-41 + 5*lmMst1) + 27*pow2(lmMgl) + pow2(lmMst1))*
509  pow2(MuSUSY)) + 2*(148 - 90*lmMt + lmMst1*(75 + 6*lmMt) + 45*pow2(
510  lmMst1) + 72*pow2(lmMt))*pow4(Mst1)))) + 2*pow2(Mt)*(5*Dmst12*pow2(
511  Mst1)*(Dmst12*(45*Dmsqst1*(12*(19 + 33*lmMst1 - 33*lmMt)*pow2(Mgl) + (-
512  3 + 28*lmMst1 - 28*lmMt)*pow2(Mst1)) + pow2(Msq)*((5308 - 7972*lmMt +
513  4*lmMst1*(1977 + 64*lmMt) - 128*pow2(lmMst1))*pow2(Mgl) - (131 + 1500*
514  lmMt - 60*lmMst1*(29 + 2*lmMt) + 306*(pow2(lmMst1) + pow2(lmMt)))*pow2(
515  Mst1))) + (90*Dmsqst1*(24*(3 + 7*lmMst1 - 7*lmMt)*pow2(Mgl) + (-19 +
516  32*lmMst1 - 24*lmMt)*pow2(Mst1)) + 4*pow2(Msq)*((3049 - 1908*lmMt - 54*
517  lmMgl*(1 + 12*lmMt) + 6*lmMst1*(215 + 49*lmMt) + 324*pow2(lmMgl) - 195*
518  pow2(lmMst1) - 144*pow2(lmMt))*pow2(Mgl) + (-196 + lmMst1*(288 - 60*
519  lmMt) + 54*lmMt + 153*pow2(lmMst1) + 153*pow2(lmMt))*pow2(Mst1)))*pow2(
520  Mst2)) + (40*pow2(Mgl)*(540*Dmsqst1*(1 + 2*lmMst1 - 2*lmMt) + (1909 +
521  600*lmMst1 + 6*(-187 + 49*lmMst1)*lmMt - 54*lmMgl*(1 + 12*lmMt) + 324*
522  pow2(lmMgl) - 195*pow2(lmMst1) - 144*pow2(lmMt))*pow2(Msq))*pow2(Mst1)
523  + (-200*Dmsqst1*(92 + 66*lmMt - 6*lmMst1*(17 + 3*lmMt) + 9*(pow2(
524  lmMst1) + pow2(lmMt))) + 30*pow2(Msq)*(-1375 + 64*lmMt - 24*(2 + 3*
525  lmMt)*pow2(lmMgl) - 4*(37 + 25*lmMt)*pow2(lmMst1) - 120*pow2(lmMt) + 8*
526  lmMgl*(26 + 12*lmMt + 9*pow2(lmMt)) + 8*lmMst1*(44 + 64*lmMt + 33*pow2(
527  lmMt)) + 24*pow3(lmMgl) + 88*pow3(lmMst1) - 276*pow3(lmMt)))*pow4(Mst1)
528  )*pow4(Mst2) - 4*pow2(Mgl)*pow2(Msq)*pow2(MuSUSY)*(pow2(Dmst12)*(163 +
529  15*lmMgl - 27*lmMst1 + 20*lmMt - 8*pow2(lmMst1)) + 20*Dmst12*(16 - 9*
530  lmMgl + 41*lmMst1 - 12*lmMt + 4*pow2(lmMst1))*pow2(Mst2) - 10*(5 + 18*
531  lmMgl - 74*lmMst1 + 24*lmMt - 8*pow2(lmMst1))*pow4(Mst2))))) + 5*
532  MuSUSY*pow2(Sbeta)*(3*Dmst12*s2t*pow2(Mst1)*(-4*Dmst12*Mgl*s2t*(60*
533  Dmsqst1*(2*(3 - lmMgl + lmMst1)*pow2(Mgl) + 3*pow2(Mst1)) + pow2(Msq)*(
534  -4*(-11 + lmMgl*(154 - 20*lmMst1) - 148*lmMst1 + 27*pow2(lmMgl) + pow2(
535  lmMst1))*pow2(Mgl) + (215 + 118*lmMst1 + 41*pow2(lmMst1))*pow2(Mst1)))
536  + Mt*(15*Dmsqst1*(Dmst12*(336*pow2(Mgl) + 31*pow2(Mst1)) + 8*(36*pow2(
537  Mgl) + (5 - 4*lmMst1)*pow2(Mst1))*pow2(Mst2)) + 4*pow2(Msq)*(Dmst12*(2*
538  (522 + 107*lmMst1)*pow2(Mgl) + (-95 - 66*lmMst1 + 82*pow2(lmMst1))*
539  pow2(Mst1)) + 2*((501 - 18*lmMgl + 504*lmMst1 - 48*lmMt + 107*pow2(
540  lmMst1))*pow2(Mgl) - 2*(-6 + 79*lmMst1 + 41*pow2(lmMst1))*pow2(Mst1))*
541  pow2(Mst2))))*pow2(Mt) + pow4(Mst1)*(-9*Mt*s2t*shiftst3*pow2(Msq)*(-8*
542  Dmst12*(-3 + 2*lmMst1)*pow2(Mst2)*pow2(Mt) + pow2(Dmst12)*(16*(-2 +
543  lmMst1)*pow2(Mt) + (1 - 2*lmMst1)*pow2(Mst1)*pow2(s2t)) + 8*(-1 + 2*
544  lmMst1)*pow2(Mt)*pow4(Mst2)) + 90*Mt*s2t*(Dmsqst1*(3 - 2*lmMst1) + (1 -
545  2*lmMst1)*pow2(Msq))*(-8*Dmst12*shiftst2*pow2(Mst2)*pow2(Mt) + (-
546  shiftst1 + shiftst2)*pow2(Dmst12)*pow2(Mst1)*pow2(s2t) + 8*(shiftst1 -
547  shiftst2)*pow2(Mt)*pow4(Mst2))) + 16*Mgl*(30*Dmsqst1*(pow2(Dmst12)*(6*(
548  31 + 27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*
549  lmMt) + 2*pow2(lmMst1))*pow2(Mgl) + (-13 - 19*lmMst1 + 19*lmMt)*pow2(
550  Mst1)) - 6*Dmst12*((31 + 27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-
551  3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))*pow2(Mgl) + (3 + 7*lmMst1 - 7*
552  lmMt)*pow2(Mst1))*pow2(Mst2) - 12*((7 + 6*lmMst1 - 5*lmMt - lmMst1*lmMt
553  + lmMgl*(-1 - lmMst1 + lmMt) + pow2(lmMst1))*pow2(Mgl) + (1 + 2*lmMst1
554  - 2*lmMt)*pow2(Mst1))*pow4(Mst2)) - pow2(Msq)*(pow2(Dmst12)*(-((-3749 -
555  2374*lmMst1 + 1770*lmMt - 117*lmMst1*lmMt + lmMgl*(96 - 237*lmMst1 +
556  157*lmMt) + 108*pow2(lmMgl) + 41*pow2(lmMst1) + 48*pow2(lmMt))*pow2(
557  Mgl)) + (lmMst1*(297 + 2*lmMt) + 15*pow2(lmMst1) + 3*(47 - 79*lmMt + 8*
558  pow2(lmMt)))*pow2(Mst1)) + Dmst12*pow2(Mst2)*(-((-605 + 6*lmMst1*(-77 +
559  lmMt) + 714*lmMt + 45*pow2(lmMst1) + 72*pow2(lmMt))*pow2(Mst1)) + 2*
560  pow2(Mgl)*(2096 - 1167*lmMt - 81*(4 + 2*lmMst1 - 3*lmMt)*pow2(lmMgl) +
561  3*(129 + 61*lmMt)*pow2(lmMst1) - 108*pow2(lmMt) - 3*lmMst1*(-793 + 67*
562  lmMt + 24*pow2(lmMt)) + 3*lmMgl*(-229 + lmMst1*(63 - 142*lmMt) + 35*
563  lmMt + 91*pow2(lmMst1) + 24*pow2(lmMt)) - 27*pow3(lmMgl) - 84*pow3(
564  lmMst1))) + 2*(-((-233 + 6*lmMst1*(-32 + lmMt) + 348*lmMt + 45*pow2(
565  lmMst1) + 72*pow2(lmMt))*pow2(Mst1)) + pow2(Mgl)*(758 - 570*lmMt - 81*(
566  3 + 2*lmMst1 - 3*lmMt)*pow2(lmMgl) + 3*(141 + 61*lmMt)*pow2(lmMst1) -
567  72*pow2(lmMt) - 6*lmMst1*(-236 + 49*lmMt + 12*pow2(lmMt)) + 3*lmMgl*(-
568  200 + 2*lmMst1 + 76*lmMt - 142*lmMst1*lmMt + 91*pow2(lmMst1) + 24*pow2(
569  lmMt)) - 27*pow3(lmMgl) - 84*pow3(lmMst1)))*pow4(Mst2)))*pow4(Mt))) -
570  25*xDmsqst1*pow2(Dmsqst1)*(Tbeta*pow2(Mt)*(3*Dmst12*s2t*(128*Mgl*Mt*(2*
571  (-3 + lmMgl - lmMst1)*pow2(Mgl) - 5*pow2(Mst1))*pow2(Mst2)*pow2(MuSUSY)
572  + Dmst12*(-192*Mgl*Mt*pow2(Mst1)*pow2(MuSUSY) + 192*s2t*pow2(Mgl)*pow2(
573  Mst1)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 15*pow2(Mst1)*pow2(Sbeta)) -
574  256*(-3 + lmMgl - lmMst1)*Mt*pow2(MuSUSY)*pow3(Mgl) + s2t*(-((-33 + 32*
575  lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta))) + 6*(87 - 32*lmMst1)*pow2(
576  Mst1)*pow2(Sbeta))*pow4(Mst1))) + 48*(2 - lmMst1)*pow4(Mst1)*(pow2(
577  Dmst12)*pow2(s2t)*(-((shiftst1 + shiftst2)*pow2(MuSUSY)*(-1 + pow2(
578  Sbeta))) + 6*(-3*shiftst1 + shiftst2)*pow2(Mst1)*pow2(Sbeta)) + 2*
579  Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)*pow2(s2t)*(-1 +
580  pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 - shiftst2)*pow2(Mst1)
581  *pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)*pow2(Mt)*pow2(
582  Sbeta)*pow4(Mst2)) + Mt*pow2(Sbeta)*(-96*Dmst12*Mgl*s2t*(Dmst12*(-6*
583  pow2(Mst1)*pow2(MuSUSY) + 4*pow2(Mgl)*((96 + 51*lmMst1 - 3*lmMgl*(5 +
584  2*lmMst1 - 2*lmMt) - 36*lmMt - 6*lmMst1*lmMt + 6*pow2(lmMst1))*pow2(
585  Mst1) + 2*(3 - lmMgl + lmMst1)*pow2(MuSUSY)) + 3*(33 + 38*lmMst1 - 38*
586  lmMt)*pow4(Mst1)) + 4*pow2(Mst2)*(-5*pow2(Mst1)*pow2(MuSUSY) + pow2(
587  Mgl)*(3*(18 - lmMgl + lmMst1)*(3 + 2*lmMst1 - 2*lmMt)*pow2(Mst1) + 2*(-
588  3 + lmMgl - lmMst1)*pow2(MuSUSY)) + 3*(9 + 14*lmMst1 - 14*lmMt)*pow4(
589  Mst1))) + Mt*pow2(Mst1)*(pow2(Dmst12)*(1728*(10 + 7*lmMst1 - 7*lmMt)*
590  pow2(Mgl) + (-425 + 2928*lmMst1 - 2928*lmMt)*pow2(Mst1)) + 72*Dmst12*(
591  12*(11 + 47*lmMst1 - 47*lmMt)*pow2(Mgl) + (-73 + 92*lmMst1 - 76*lmMt)*
592  pow2(Mst1))*pow2(Mst2) + 16*(216*(1 + 7*lmMst1 - 7*lmMt)*pow2(Mgl) - (
593  559 + 384*lmMt - 12*lmMst1*(44 + 3*lmMt) + 18*(pow2(lmMst1) + pow2(
594  lmMt)))*pow2(Mst1))*pow4(Mst2)))) + 6*MuSUSY*pow2(Sbeta)*(pow3(Mt)*(-6*
595  Dmst12*pow2(Mst2)*(16*(11 + 47*lmMst1 - 47*lmMt)*Mgl*Mt*pow2(Mst1) -
596  480*s2t*pow2(Mgl)*pow2(Mst1) + 64*Mt*(31 + 27*lmMst1 - 24*lmMt - 2*
597  lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))*pow3(
598  Mgl) + s2t*(-87 - 32*lmMst1*(-1 + shiftst2) + 64*shiftst2)*pow4(Mst1))
599  - 192*(2*Mt*((1 + 7*lmMst1 - 7*lmMt)*Mgl*pow2(Mst1) + (18 - lmMgl +
600  lmMst1)*(1 + lmMst1 - lmMt)*pow3(Mgl)) + (-2 + lmMst1)*s2t*(shiftst1 -
601  shiftst2)*pow4(Mst1))*pow4(Mst2)) + Mt*pow2(Dmst12)*(864*s2t*pow2(Mgl)*
602  pow2(Mst1)*pow2(Mt) - 192*pow3(Mgl)*((3 - lmMgl + lmMst1)*Mt*pow2(Mst1)
603  *pow2(s2t) + (-62 - 54*lmMst1 + lmMgl*(6 + 4*lmMst1 - 4*lmMt) + 48*lmMt
604  + 4*lmMst1*lmMt - 4*pow2(lmMst1))*pow3(Mt)) + 321*s2t*pow2(Mt)*pow4(
605  Mst1) - 32*Mgl*((29 - 19*lmMst1 + 19*lmMt)*pow2(Mst1)*pow3(Mt) + 15*Mt*
606  pow2(s2t)*pow4(Mst1)) + 24*(-2 + lmMst1)*(shiftst1 - shiftst2)*pow3(
607  s2t)*pow6(Mst1))))) + 540*pow2(Mt)*pow4(Msq)*(16*Al4p*twoLoopFlag*(
608  MuSUSY*pow4(Mst1)*(Dmst12*MuSUSY*s2t*Tbeta*(4*Mt*pow2(Mst2)*((2 +
609  lmMst1)*Mgl*pow2(Mst1) + (1 - 2*lmMgl + 2*lmMst1)*pow3(Mgl)) + Dmst12*
610  s2t*((3 + 2*lmMst1)*pow2(Mgl)*pow2(Mst1) - (1 + 2*lmMst1)*pow4(Mst1)))
611  + 2*pow2(Sbeta)*(12*Dmst12*Mt*pow2(Mst2)*(-((1 + lmMt)*Mgl*Mt*pow2(
612  Mst1)) - (2 + lmMst1)*s2t*pow2(Mgl)*pow2(Mst1) + Mt*(2 + 6*lmMst1 - 3*
613  lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))
614  *pow3(Mgl) + (1 + 2*lmMst1)*s2t*pow4(Mst1)) - pow2(Dmst12)*(6*Mt*s2t*
615  pow2(Mgl)*pow2(Mst1) + (8*(-3 + lmMgl - 2*lmMst1 + lmMt)*pow2(Mt) + 3*(
616  -1 + 2*lmMgl - 2*lmMst1)*pow2(Mst1)*pow2(s2t))*pow3(Mgl) + 6*(-1 + 2*
617  lmMst1)*Mt*s2t*pow4(Mst1) - Mgl*(4*(1 + lmMt)*pow2(Mst1)*pow2(Mt) + 3*(
618  2 + lmMst1)*pow2(s2t)*pow4(Mst1))) + 24*Mgl*((-(lmMst1*(-2 + lmMt)) -
619  lmMt + lmMgl*(-1 - lmMst1 + lmMt) + pow2(lmMst1))*pow2(Mgl) - (1 +
620  lmMt)*pow2(Mst1))*pow2(Mt)*pow4(Mst2))) + Tbeta*pow4(Mst1)*(s2t*pow2(
621  Dmst12)*(2*Mt*pow2(MuSUSY)*(Mgl*pow2(Mst1) + (3 - 2*lmMgl + 2*lmMst1)*
622  pow3(Mgl)) + s2t*pow2(Sbeta)*((1 + 2*lmMst1)*(12*pow2(Mst1) + pow2(
623  MuSUSY))*pow4(Mst1) - pow2(Mgl)*((3 + 2*lmMst1)*pow2(Mst1)*pow2(MuSUSY)
624  + 12*(2 + lmMst1)*pow4(Mst1)))) + 2*Mt*pow2(Sbeta)*(6*(lmMst1 + lmMt)*
625  Mt*pow2(Dmst12)*pow4(Mst1) + Mgl*s2t*pow2(Dmst12)*(-(pow2(Mst1)*pow2(
626  MuSUSY)) + pow2(Mgl)*(6*(5 + 4*lmMst1 - 2*(lmMgl + lmMt))*pow2(Mst1) +
627  (-3 + 2*lmMgl - 2*lmMst1)*pow2(MuSUSY)) + 6*(1 + 2*lmMt)*pow4(Mst1)) -
628  2*Dmst12*Mgl*s2t*pow2(Mst2)*((2 + lmMst1)*pow2(Mst1)*pow2(MuSUSY) +
629  pow2(Mgl)*(6*(-1 - 5*lmMst1 + lmMgl*(3 + 2*lmMst1 - 2*lmMt) + 2*lmMt +
630  2*lmMst1*lmMt - 2*pow2(lmMst1))*pow2(Mst1) + (1 - 2*lmMgl + 2*lmMst1)*
631  pow2(MuSUSY)) + 6*(1 + 2*lmMt)*pow4(Mst1)) + 12*Mt*pow2(Mst1)*(Dmst12*(
632  (1 + lmMt)*pow2(Mgl) - (1 + lmMst1 + lmMt)*pow2(Mst1))*pow2(Mst2) + (2*
633  (1 + lmMt)*pow2(Mgl) - (2 + 2*lmMst1*(1 + lmMt) + pow2(lmMst1) - 3*
634  pow2(lmMt))*pow2(Mst1))*pow4(Mst2)))) - x*(3*MuSUSY*pow2(Sbeta)*(2*
635  Dmst12*(-1 + 3*lmMgl - 3*lmMst1)*s2t*pow2(Mst1)*(4*Dmst12*Mgl*Mt +
636  Dmst12*s2t*pow2(Mst1) - 4*Mgl*Mt*pow2(Mst2)) + 8*pow2(Mt)*(Dmst12*(
637  lmMst1*(22 - 6*lmMt) - 5*lmMt + lmMgl*(-17 - 6*lmMst1 + 6*lmMt) + 6*
638  pow2(lmMst1))*pow2(Mst1)*(Dmst12 - pow2(Mst2)) + ((9 + lmMgl*(20 + 6*
639  lmMst1 - 6*lmMt) + 2*lmMt + lmMst1*(-22 + 6*lmMt) - 6*pow2(lmMst1))*
640  pow2(Mgl) + 2*(2 + lmMgl*(5 + 2*lmMst1 - 2*lmMt) + 2*lmMst1*(-3 + lmMt)
641  + lmMt - 2*pow2(lmMst1))*pow2(Mst1))*pow4(Mst2))) + Tbeta*(4*Dmst12*(1
642  - 3*lmMgl + 3*lmMst1)*s2t*pow2(Mst1)*(Dmst12*Mt*pow2(MuSUSY) - Mt*pow2(
643  Mst2)*pow2(MuSUSY) + 3*Dmst12*Mgl*s2t*pow2(Mst1)*pow2(Sbeta)) + 4*Mt*
644  pow2(Sbeta)*(Dmst12*s2t*pow2(Mst1)*(pow2(Mst2)*(6*(8 - 25*lmMst1 +
645  lmMgl*(23 + 6*lmMst1 - 6*lmMt) + 2*lmMt + 6*lmMst1*lmMt - 6*pow2(
646  lmMst1))*pow2(Mgl) + 6*(3 - 14*lmMst1 + 4*lmMgl*(3 + lmMst1 - lmMt) +
647  2*lmMt + 4*lmMst1*lmMt - 4*pow2(lmMst1))*pow2(Mst1) + (1 - 3*lmMgl + 3*
648  lmMst1)*pow2(MuSUSY)) - Dmst12*(6*(8 - 25*lmMst1 + lmMgl*(23 + 6*lmMst1
649  - 6*lmMt) + 2*lmMt + 6*lmMst1*lmMt - 6*pow2(lmMst1))*pow2(Mgl) + 3*(11
650  + lmMst1*(22 - 4*lmMt) - 4*lmMgl*(4 + lmMst1 - lmMt) - 6*lmMt + 4*pow2(
651  lmMst1))*pow2(Mst1) + (1 - 3*lmMgl + 3*lmMst1)*pow2(MuSUSY))) - 6*Mgl*
652  Mt*(Dmst12*(lmMst1*(22 - 6*lmMt) - 5*lmMt + lmMgl*(-17 - 6*lmMst1 + 6*
653  lmMt) + 6*pow2(lmMst1))*pow2(Mst1)*(Dmst12 - pow2(Mst2)) + ((9 + lmMgl*
654  (20 + 6*lmMst1 - 6*lmMt) + 2*lmMt + lmMst1*(-22 + 6*lmMt) - 6*pow2(
655  lmMst1))*pow2(Mgl) + 2*(2 + lmMgl*(5 + 2*lmMst1 - 2*lmMt) + 2*lmMst1*(-
656  3 + lmMt) + lmMt - 2*pow2(lmMst1))*pow2(Mst1))*pow4(Mst2)))))*pow5(Mgl)
657  ) + 3*oneLoopFlag*(s2t*pow2(Dmst12)*(-12*Mt*MuSUSY*pow2(Sbeta) + s2t*
658  Tbeta*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*pow2(Mst1)*pow2(Sbeta))) +
659  12*Mt*pow2(Sbeta)*(2*Dmst12*MuSUSY*s2t*pow2(Mst2) + Mt*Tbeta*(pow2(
660  Dmst12) - 2*Dmst12*pow2(Mst2) + 4*(-lmMst1 + lmMt)*pow4(Mst2))))*pow8(
661  Mst1))) + pow2(Mst1)*(40*Al4p*Mt*xMgl*pow4(Mgl)*(216*twoLoopFlag*
662  xDmst12*pow2(Mst1)*pow3(Dmst12)*(Tbeta*(-2*Mt*pow2(s2t)*((-1 + 3*lmMgl
663  - 3*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*(1 + 2*lmMgl - 2*
664  lmMst1)*pow2(Mst1)*pow2(Sbeta)) + 48*(5 - 2*lmMgl*(2 + lmMst1 - lmMt) +
665  2*lmMst1*(4 + lmMst1 - lmMt) - 4*lmMt)*pow2(Sbeta)*pow3(Mt)) + MuSUSY*
666  pow2(Sbeta)*(96*(1 - lmMgl + lmMst1)*s2t*pow2(Mt) + (1 - 3*lmMgl + 3*
667  lmMst1)*pow2(Mst1)*pow3(s2t)))*pow4(Msq) + Al4p*threeLoopFlag*(xDmst12*
668  pow2(Msq)*pow3(Dmst12)*(-2*Mt*MuSUSY*s2t*pow2(Mst1)*(MuSUSY*s2t*Tbeta*(
669  41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl)
670  + 540*pow2(lmMst1))*pow2(Msq) + 12960*Dmsqst1*(29 - 10*lmMgl + 10*
671  lmMst1)*Mt*pow2(Sbeta) + 864*Mt*(-1797 - 1711*lmMst1 + lmMgl*(1315 -
672  198*lmMst1 - 96*lmMt) + 144*lmMt + 96*lmMst1*lmMt + 333*pow2(lmMgl) -
673  103*pow2(lmMst1))*pow2(Msq)*pow2(Sbeta)) + Tbeta*pow2(Sbeta)*(2*Mt*
674  pow2(Mst1)*(12960*Dmsqst1*(29 - 10*lmMgl + 10*lmMst1)*pow2(Mst1) - 108*
675  (1695 + 8*lmMst1*(-125 + 12*lmMt) - 2*lmMgl*(-605 + 99*lmMst1 + 48*
676  lmMt) + 333*pow2(lmMgl) - 103*pow2(lmMst1))*pow2(Msq)*pow2(Mst1) + (
677  41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl)
678  + 540*pow2(lmMst1))*pow2(Msq)*pow2(MuSUSY))*pow2(s2t) - 288*(45*
679  Dmsqst1*(307 + 4*lmMst1*(67 - 5*lmMt) - 10*lmMgl*(3 + 2*lmMst1 - 2*
680  lmMt) - 238*lmMt + 20*pow2(lmMst1))*pow2(Mst1) + pow2(Msq)*(-2*(262 +
681  lmMst1*(335 - 72*lmMt) - 96*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt)
682  + 120*pow2(lmMst1))*pow2(MuSUSY) + pow2(Mst1)*(72*(49 + 18*lmMst1 - 39*
683  lmMt)*pow2(lmMgl) - 16*(181 + 108*lmMt)*pow2(lmMst1) + 8*lmMgl*(868 +
684  287*lmMt + 9*lmMst1*(-59 + 63*lmMt) - 306*pow2(lmMst1) - 72*pow2(lmMt))
685  + 8*lmMst1*(-4817 - 3*lmMt + 72*pow2(lmMt)) + 3*(-17921 + 7796*lmMt +
686  384*pow2(lmMt)) + 504*pow3(lmMgl) + 648*pow3(lmMst1))))*pow3(Mt)) +
687  MuSUSY*pow2(Msq)*(576*MuSUSY*Tbeta*(-262 - 335*lmMst1 + lmMgl*(199 +
688  120*lmMst1 - 72*lmMt) + 96*lmMt + 72*lmMst1*lmMt - 120*pow2(lmMst1))*
689  pow3(Mt) + (41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) -
690  12636*pow2(lmMgl) + 540*pow2(lmMst1))*pow2(Sbeta)*pow3(s2t)*pow4(Mst1))
691  ) - 12960*Mt*xDmsqst1*pow2(Dmsqst1)*pow2(Mst1)*pow2(Mst2)*pow2(Sbeta)*(
692  Dmst12*Mt*(2*(29 - 10*lmMgl + 10*lmMst1)*MuSUSY*s2t + Mt*Tbeta*(307 +
693  268*lmMst1 - 10*lmMgl*(3 + 2*lmMst1 - 2*lmMt) - 238*lmMt - 20*lmMst1*
694  lmMt + 20*pow2(lmMst1)))*pow2(Mst2) - pow2(Dmst12)*(2*(29 - 10*lmMgl +
695  10*lmMst1)*Mt*MuSUSY*s2t + Tbeta*(307 + 4*lmMst1*(67 - 5*lmMt) - 10*
696  lmMgl*(3 + 2*lmMst1 - 2*lmMt) - 238*lmMt + 20*pow2(lmMst1))*pow2(Mt) +
697  (-29 + 10*lmMgl - 10*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t)) + (179 - 10*
698  lmMgl + 10*lmMst1)*(1 + lmMst1 - lmMt)*Tbeta*pow2(Mt)*pow4(Mst2))) +
699  Mt*pow2(Mst2)*(-216*twoLoopFlag*pow2(Mst1)*(24*Mt*pow2(Mst2)*(Dmst12*((
700  1 - 2*lmMgl + 2*lmMst1)*MuSUSY*s2t + Mt*Tbeta*(2 + 6*lmMst1 - 3*lmMt -
701  2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))) - 2*
702  Mt*Tbeta*(lmMgl*(1 + lmMst1 - lmMt) + lmMst1*(-2 + lmMt) + lmMt - pow2(
703  lmMst1))*pow2(Mst2))*pow2(Sbeta) + pow2(Dmst12)*(12*Mt*((3 - 2*lmMgl +
704  2*lmMst1)*MuSUSY*s2t + Mt*Tbeta*(8 - 2*lmMst1*(-5 + lmMt) - 5*lmMt +
705  lmMgl*(-5 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1)))*pow2(Sbeta) + Tbeta*
706  pow2(s2t)*(-((-1 + 3*lmMgl - 3*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)))
707  + 12*(1 - 2*lmMgl + 2*lmMst1)*pow2(Mst1)*pow2(Sbeta))))*pow4(Msq) +
708  Al4p*threeLoopFlag*pow2(Msq)*(Dmst12*s2t*pow2(Mst1)*(-(Dmst12*s2t*
709  Tbeta*(12960*Dmsqst1*(29 - 10*lmMgl + 10*lmMst1)*pow2(Mst1) + pow2(Msq)
710  *(432*(681 + 1506*lmMst1 - 96*(1 + lmMst1)*lmMt + 2*lmMgl*(-656 + 99*
711  lmMst1 + 48*lmMt) - 333*pow2(lmMgl) + 103*pow2(lmMst1))*pow2(Mst1) + (
712  41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl)
713  + 540*pow2(lmMst1))*pow2(MuSUSY)))*pow2(Sbeta)) + MuSUSY*(Dmst12*
714  MuSUSY*s2t*Tbeta*(41653 + 11784*lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) -
715  12636*pow2(lmMgl) + 540*pow2(lmMst1))*pow2(Msq) + 432*Mt*(-60*Dmsqst1*(
716  -29 + 10*lmMgl - 10*lmMst1)*(Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(
717  -2961 - 1964*lmMst1 + 192*lmMt + 96*lmMst1*lmMt - 2*lmMgl*(-683 + 99*
718  lmMst1 + 48*lmMt) + 333*pow2(lmMgl) - 103*pow2(lmMst1)) + 2*(-633 -
719  1458*lmMst1 + 96*lmMt + 96*lmMst1*lmMt - 2*lmMgl*(-632 + 99*lmMst1 +
720  48*lmMt) + 333*pow2(lmMgl) - 103*pow2(lmMst1))*pow2(Mst2)))*pow2(Sbeta)
721  )) + 144*Tbeta*pow2(Mt)*(4*pow2(Msq)*pow2(MuSUSY)*(Dmst12*(262 +
722  lmMst1*(335 - 72*lmMt) - 96*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt)
723  + 120*pow2(lmMst1))*(Dmst12 - pow2(Mst2)) - 2*(59 + 85*lmMst1 - lmMgl*(
724  53 + 40*lmMst1 - 24*lmMt) - 24*(1 + lmMst1)*lmMt + 40*pow2(lmMst1))*
725  pow4(Mst2)) + pow2(Sbeta)*(90*Dmsqst1*pow2(Mst1)*(Dmst12*(307 + 4*
726  lmMst1*(67 - 5*lmMt) - 10*lmMgl*(3 + 2*lmMst1 - 2*lmMt) - 238*lmMt +
727  20*pow2(lmMst1))*(Dmst12 - pow2(Mst2)) - 2*(69 + lmMst1*(59 - 10*lmMt)
728  - 10*lmMgl*(1 + lmMst1 - lmMt) - 49*lmMt + 10*pow2(lmMst1))*pow4(Mst2))
729  - pow2(Msq)*(pow2(Dmst12)*(4*(262 + lmMst1*(335 - 72*lmMt) - 96*lmMt +
730  lmMgl*(-199 - 120*lmMst1 + 72*lmMt) + 120*pow2(lmMst1))*pow2(MuSUSY) +
731  pow2(Mst1)*(40279 - 17847*lmMt - 6*(303 + 108*lmMst1 - 203*lmMt)*pow2(
732  lmMgl) + (952 + 870*lmMt)*pow2(lmMst1) + lmMst1*(26599 + 804*lmMt -
733  288*pow2(lmMt)) - 720*pow2(lmMt) + 2*lmMgl*(-1474 - 977*lmMt - 12*
734  lmMst1*(-118 + 87*lmMt) + 615*pow2(lmMst1) + 144*pow2(lmMt)) - 190*
735  pow3(lmMgl) - 392*pow3(lmMst1))) + 2*(Dmst12*pow2(Mst2)*(-2*(262 +
736  lmMst1*(335 - 72*lmMt) - 96*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt)
737  + 120*pow2(lmMst1))*pow2(MuSUSY) + pow2(Mst1)*(13484 - 5541*lmMt + 6*(
738  324 + 143*lmMt)*pow2(lmMst1) - 6*lmMgl*(666 + lmMgl*(285 + 108*lmMst1 -
739  265*lmMt) + 57*lmMt + 4*lmMst1*(-59 + 102*lmMt) - 203*pow2(lmMst1) -
740  48*pow2(lmMt)) - 432*pow2(lmMt) - 3*lmMst1*(-3979 + 260*lmMt + 96*pow2(
741  lmMt)) - 314*pow3(lmMgl) - 256*pow3(lmMst1))) + (-4*(59 + lmMst1*(85 -
742  24*lmMt) - 24*lmMt + lmMgl*(-53 - 40*lmMst1 + 24*lmMt) + 40*pow2(
743  lmMst1))*pow2(MuSUSY) + pow2(Mst1)*(6677 - 2742*lmMt + 6*(313 + 143*
744  lmMt)*pow2(lmMst1) - 6*lmMgl*(645 + lmMgl*(231 + 108*lmMst1 - 265*lmMt)
745  - 72*lmMt + 12*lmMst1*(-9 + 34*lmMt) - 203*pow2(lmMst1) - 48*pow2(lmMt)
746  ) - 288*pow2(lmMt) - 12*lmMst1*(-617 + 99*lmMt + 24*pow2(lmMt)) - 314*
747  pow3(lmMgl) - 256*pow3(lmMst1)))*pow4(Mst2)))))))) + 6*xDmst12*pow3(
748  Dmst12)*(Mt*pow4(Msq)*(270*oneLoopFlag*(-(Mt*Tbeta*pow2(s2t)*(pow2(
749  MuSUSY)*(-1 + pow2(Sbeta)) + 18*pow2(Mst1)*pow2(Sbeta))) - pow2(Sbeta)*
750  (-8*MuSUSY*s2t*pow2(Mt) + 8*Tbeta*pow3(Mt) + MuSUSY*pow2(Mst1)*pow3(
751  s2t)))*pow6(Mst1) - 144*Al4p*twoLoopFlag*(-20*x*(Tbeta*((-12*Mgl*(-((
752  lmMgl*(17 + 6*lmMst1 - 6*lmMt) + 5*lmMt + lmMst1*(-22 + 6*lmMt))*Mt) +
753  (8 - 25*lmMst1 + lmMgl*(23 + 6*lmMst1 - 6*lmMt) + 2*lmMt + 6*lmMst1*
754  lmMt)*Mgl*s2t + 6*(Mt - Mgl*s2t)*pow2(lmMst1))*pow2(Mt) + (1 - 3*lmMgl
755  + 3*lmMst1)*pow2(Mst1)*(12*Mgl*Mt + s2t*pow2(Mst1))*pow2(s2t))*pow2(
756  Sbeta) + s2t*pow2(Mt)*(2*(-1 + 3*lmMgl - 3*lmMst1)*pow2(MuSUSY)*(-1 +
757  pow2(Sbeta)) + 48*(-2 - 9*lmMst1 + lmMgl*(7 + 2*lmMst1 - 2*lmMt) + 2*
758  lmMt + 2*lmMst1*lmMt - 2*pow2(lmMst1))*pow2(Mst1)*pow2(Sbeta))) +
759  MuSUSY*pow2(Sbeta)*(6*(-1 + 3*lmMgl - 3*lmMst1)*Mt*s2t*(2*Mgl*Mt + s2t*
760  pow2(Mst1)) + 12*(lmMst1*(22 - 6*lmMt) - 5*lmMt + lmMgl*(-17 - 6*lmMst1
761  + 6*lmMt) + 6*pow2(lmMst1))*pow3(Mt)))*pow5(Mgl) + Tbeta*pow2(Sbeta)*(
762  2*Mgl*s2t*pow2(Mst1)*pow2(Mt)*(-((9 + 2*lmMst1)*pow2(Mst1)*pow2(MuSUSY)
763  ) + 20*pow2(Mgl)*((1 + 4*lmMst1 - 2*(lmMgl + lmMt))*pow2(Mst1) + 4*(-1
764  + lmMgl - lmMst1)*pow2(MuSUSY)) + 40*(1 + 2*lmMt)*pow4(Mst1)) - 5*Mt*
765  pow2(s2t)*pow4(Mst1)*(pow2(Mgl)*(12*(3 + 2*lmMst1)*pow2(Mst1) + (1 + 2*
766  lmMst1)*pow2(MuSUSY)) - 4*(lmMst1*pow2(Mst1)*pow2(MuSUSY) + 3*(1 + 6*
767  lmMst1)*pow4(Mst1))) + (40*(-1 + 2*lmMst1 + 2*lmMt)*pow3(Mt) - 20*((2 +
768  lmMst1)*Mgl*pow2(Mst1) + (1 - 2*lmMgl + 2*lmMst1)*pow3(Mgl))*pow3(s2t))
769  *pow6(Mst1)) + MuSUSY*pow2(Mst1)*(Mt*MuSUSY*s2t*Tbeta*(2*(9 + 2*lmMst1)
770  *Mgl*Mt*pow2(Mst1) + 5*(1 + 2*lmMst1)*s2t*pow2(Mgl)*pow2(Mst1) + 160*(1
771  - lmMgl + lmMst1)*Mt*pow3(Mgl) - 20*lmMst1*s2t*pow4(Mst1)) + 10*pow2(
772  Sbeta)*(pow3(Mgl)*(-3*(1 + 2*lmMgl - 2*lmMst1)*Mt*pow2(Mst1)*pow2(s2t)
773  + 8*(18 + 26*lmMst1 - 13*lmMt - 6*lmMst1*lmMt + lmMgl*(-13 - 6*lmMst1 +
774  6*lmMt) + 6*pow2(lmMst1))*pow3(Mt)) + 2*(9 - 8*lmMst1)*s2t*pow2(Mt)*
775  pow4(Mst1) + Mgl*(4*(1 + lmMt)*pow2(Mst1)*pow3(Mt) + 3*(3 + 2*lmMst1)*
776  Mt*pow2(s2t)*pow4(Mst1)) - pow2(Mgl)*(4*s2t*pow2(Mst1)*pow2(Mt) + (3 +
777  2*lmMst1)*pow3(s2t)*pow4(Mst1)) + (1 + 2*lmMst1)*pow3(s2t)*pow6(Mst1)))
778  )) + threeLoopFlag*pow2(Al4p)*pow2(Mst1)*(5*xDmsqst1*pow2(Dmsqst1)*(
779  5760*Mgl*Mt*(pow2(Mst1)*(4*s2t*pow2(Mt)*(8*Tbeta*pow2(MuSUSY)*(-1 +
780  pow2(Sbeta)) + 72*Mgl*MuSUSY*pow2(Sbeta) + (-15 - 38*lmMst1 + 38*lmMt)*
781  Tbeta*pow2(Mst1)*pow2(Sbeta)) - 3*Mt*pow2(s2t)*(-4*Mgl*Tbeta*pow2(
782  MuSUSY)*(-1 + pow2(Sbeta)) + 7*(MuSUSY - 3*Mgl*Tbeta)*pow2(Mst1)*pow2(
783  Sbeta)) + (-2*(91 + 103*lmMst1 - 103*lmMt)*MuSUSY + 27*(17 + 25*lmMst1
784  - 25*lmMt)*Mgl*Tbeta)*pow2(Sbeta)*pow3(Mt) + 2*pow2(Mst1)*(3*Mgl*MuSUSY
785  + 5*Tbeta*pow2(Mst1))*pow2(Sbeta)*pow3(s2t)) + 4*pow2(Mgl)*(MuSUSY*
786  pow2(Sbeta)*(6*(-3 + lmMgl - lmMst1)*Mt*pow2(Mst1)*pow2(s2t) + 6*(31 +
787  27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) +
788  2*pow2(lmMst1))*pow3(Mt)) + Tbeta*(-(s2t*pow2(Mt)*(-2*(-3 + lmMgl -
789  lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*(118 + lmMst1*(73 - 6*lmMt)
790  - 60*lmMt + lmMgl*(-13 - 6*lmMst1 + 6*lmMt) + 6*pow2(lmMst1))*pow2(
791  Mst1)*pow2(Sbeta))) + (3 - lmMgl + lmMst1)*pow2(Sbeta)*pow3(s2t)*pow4(
792  Mst1)))) + pow4(Mst1)*(36*Tbeta*pow2(Mt)*pow2(s2t)*((215 - 72*shiftst1
793  - 208*shiftst2 + 16*lmMst1*(-5 + 3*shiftst1 + 2*shiftst2))*pow2(MuSUSY)
794  *(-1 + pow2(Sbeta)) + 5*(201 - 512*shiftst1 + 128*shiftst2 - 64*lmMst1*
795  (3 - 4*shiftst1 + shiftst2))*pow2(Mst1)*pow2(Sbeta)) + 5*pow2(Sbeta)*(
796  27*MuSUSY*s2t*(1389 - 256*shiftst2)*pow3(Mt) + 36*Mt*MuSUSY*(33 - 128*
797  shiftst1 + 64*shiftst2 - 32*lmMst1*(1 - 2*shiftst1 + shiftst2))*pow2(
798  Mst1)*pow3(s2t) + 16*(-2657 + 2928*lmMst1 - 2928*lmMt + 432*shiftst2)*
799  Tbeta*pow4(Mt) - 576*(-2 + lmMst1)*(shiftst1 - shiftst2)*Tbeta*pow4(
800  Mst1)*pow4(s2t)))) + pow2(Msq)*(3*MuSUSY*pow2(Sbeta)*(s2t*pow2(Mst1)*(
801  Mt*(225*Dmsqst1*(15360*pow2(Mgl) - 83*pow2(Mst1)) - 8*pow2(Msq)*(8*(
802  4113 + 90*lmMgl - 1792*lmMst1 + 120*lmMt - 48*pow2(lmMst1))*pow2(Mgl) +
803  5*(2615 + 2600*lmMst1 - 1312*pow2(lmMst1))*pow2(Mst1))) - 960*Mgl*s2t*(
804  -30*Dmsqst1*(8*(-3 + lmMgl - lmMst1)*pow2(Mgl) + pow2(Mst1)) + pow2(
805  Msq)*(-2*(227 + 144*lmMgl - 120*lmMst1 - 20*lmMgl*lmMst1 + 27*pow2(
806  lmMgl) + pow2(lmMst1))*pow2(Mgl) + (66 + 77*lmMst1 + 41*pow2(lmMst1))*
807  pow2(Mst1))))*pow2(Mt) + 80*(Mt*(15*Dmsqst1*(48*pow2(Mgl) + (3 - 24*
808  shiftst1 + 12*shiftst2 - 8*lmMst1*(1 - 2*shiftst1 + shiftst2))*pow2(
809  Mst1)) + 2*pow2(Msq)*((200 + 18*lmMgl + 265*lmMst1 + 91*pow2(lmMst1))*
810  pow2(Mgl) - (16 + 60*shiftst1 - 30*shiftst2 + lmMst1*(115 - 120*
811  shiftst1 + 60*shiftst2) + 41*pow2(lmMst1))*pow2(Mst1)))*pow3(s2t) + 3*
812  shiftst3*pow2(Msq)*(16*(7 - 3*lmMst1)*s2t*pow3(Mt) + (-13 + 14*lmMst1)*
813  Mt*pow2(Mst1)*pow3(s2t)))*pow4(Mst1) + 640*Mgl*(360*Dmsqst1*(31 +
814  lmMst1*(27 - 2*lmMt) - 24*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*
815  pow2(lmMst1))*pow2(Mgl) - 240*Dmsqst1*(11 + 20*lmMst1 - 20*lmMt)*pow2(
816  Mst1) - (-233 + 47*lmMt + lmMst1*(-3 + 2*lmMt) + 15*pow2(lmMst1) + 24*
817  pow2(lmMt))*pow2(Msq)*pow2(Mst1) - 4*pow2(Mgl)*pow2(Msq)*(5845 - 2937*
818  lmMt - 27*(16 + lmMgl + 6*lmMst1 - 9*lmMt)*pow2(lmMgl) + (346 + 183*
819  lmMt)*pow2(lmMst1) + lmMst1*(4753 - 84*lmMt - 72*pow2(lmMt)) - 156*
820  pow2(lmMt) + lmMgl*(-783 - 426*lmMst1*(-1 + lmMt) - 52*lmMt + 273*pow2(
821  lmMst1) + 72*pow2(lmMt)) - 84*pow3(lmMst1)))*pow4(Mt)) + Tbeta*(-60*
822  pow2(Mst1)*pow2(Mt)*pow2(s2t)*((15*Dmsqst1*(384*pow2(Mgl) + (1 - 16*
823  lmMst1)*pow2(Mst1)) + 2*pow2(Msq)*(2*(-425 + 18*lmMgl + 83*lmMst1 + 91*
824  pow2(lmMst1))*pow2(Mgl) + (121 - 176*lmMst1 - 164*pow2(lmMst1))*pow2(
825  Mst1)))*pow2(MuSUSY) + 24*shiftst1*(Dmsqst1*(3 - 2*lmMst1) + (1 - 2*
826  lmMst1)*pow2(Msq))*pow2(Mst1)*(80*pow2(Mst1) + 3*pow2(MuSUSY))*pow2(
827  Sbeta)) + 64*Mgl*(12*Mgl*Mt*(323 - 75*lmMgl + 383*lmMst1 - 100*lmMt +
828  32*pow2(lmMst1))*pow2(Msq) + 240*s2t*pow2(Mgl)*(15*Dmsqst1*(3 - lmMgl +
829  lmMst1) + (-130 + lmMgl*(159 - 20*lmMst1) - 162*lmMst1 + 27*pow2(lmMgl)
830  + pow2(lmMst1))*pow2(Msq)) - s2t*(18000*Dmsqst1 + (-365 + 1049*lmMst1 +
831  123*pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*pow2(MuSUSY)*pow3(Mt) + 144*
832  pow4(Mst1)*(shiftst3*pow2(Msq)*(2*pow2(Mt)*pow2(s2t)*((-29 + 18*lmMst1)
833  *pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 5*(-83 + 58*lmMst1)*pow2(Mst1)*pow2(
834  Sbeta)) - 80*(-7 + 3*lmMst1)*pow2(Sbeta)*pow4(Mt) + 5*(1 - 2*lmMst1)*
835  pow2(Sbeta)*pow4(Mst1)*pow4(s2t)) + 10*(Dmsqst1*(3 - 2*lmMst1) + (1 -
836  2*lmMst1)*pow2(Msq))*(pow2(Mt)*pow2(s2t)*(20*shiftst2*pow2(Mst1)*pow2(
837  Sbeta) + pow2(MuSUSY)*(3*shiftst1 + 2*shiftst2 - 2*shiftst2*pow2(Sbeta)
838  )) - 5*shiftst2*pow2(Sbeta)*pow4(Mst1)*pow4(s2t))) + pow2(Sbeta)*(60*
839  pow2(Mst1)*pow2(Mt)*pow2(s2t)*(15*Dmsqst1*(-96*pow2(Mgl)*(3*pow2(Mst1)
840  - 4*pow2(MuSUSY)) + (1 - 16*lmMst1)*pow2(Mst1)*pow2(MuSUSY) + 6*(9 -
841  32*lmMst1)*pow4(Mst1)) + 2*pow2(Msq)*((121 - 176*lmMst1 - 164*pow2(
842  lmMst1))*pow2(Mst1)*pow2(MuSUSY) + 2*pow2(Mgl)*(12*(3 - 18*lmMgl + 397*
843  lmMst1 - 48*lmMt + 107*pow2(lmMst1))*pow2(Mst1) + (-425 + 18*lmMgl +
844  83*lmMst1 + 91*pow2(lmMst1))*pow2(MuSUSY)) - 12*(-119 + 250*lmMst1 +
845  246*pow2(lmMst1))*pow4(Mst1))) + Mgl*(960*Mt*(60*Dmsqst1*(2*(3 - lmMgl
846  + lmMst1)*pow2(Mgl) + 3*pow2(Mst1)) + pow2(Msq)*(-4*(-11 + lmMgl*(154 -
847  20*lmMst1) - 148*lmMst1 + 27*pow2(lmMgl) + pow2(lmMst1))*pow2(Mgl) + (
848  215 + 118*lmMst1 + 41*pow2(lmMst1))*pow2(Mst1)))*pow3(s2t)*pow4(Mst1) +
849  64*s2t*pow3(Mt)*(900*Dmsqst1*(20*pow2(Mst1)*pow2(MuSUSY) + 4*pow2(Mgl)*
850  (3*(-67 - 43*lmMst1 + lmMgl*(5 + 2*lmMst1 - 2*lmMt) + 38*lmMt + 2*
851  lmMst1*lmMt - 2*pow2(lmMst1))*pow2(Mst1) + (-3 + lmMgl - lmMst1)*pow2(
852  MuSUSY)) + (15 + 2*lmMst1 - 2*lmMt)*pow4(Mst1)) + pow2(Msq)*((-365 +
853  1049*lmMst1 + 123*pow2(lmMst1))*pow2(Mst1)*pow2(MuSUSY) + 30*pow2(Mgl)*
854  ((2722 + 965*lmMst1 - 6*(200 + 23*lmMst1)*lmMt + lmMgl*(67 - 258*lmMst1
855  + 178*lmMt) + 108*pow2(lmMgl) + 62*pow2(lmMst1) + 48*pow2(lmMt))*pow2(
856  Mst1) - 8*(-130 + lmMgl*(159 - 20*lmMst1) - 162*lmMst1 + 27*pow2(lmMgl)
857  + pow2(lmMst1))*pow2(MuSUSY)) - 30*(-223 + 342*lmMt + lmMst1*(-442 + 8*
858  lmMt) + 60*pow2(lmMst1) + 96*pow2(lmMt))*pow4(Mst1)))) + 8*(25*Dmsqst1*
859  (2592*(25 + 47*lmMst1 - 47*lmMt)*pow2(Mgl)*pow2(Mst1) + (-907 + 96*
860  lmMst1 - 96*lmMt)*pow4(Mst1)) - 4*pow2(Msq)*(12*pow2(Mgl)*(5*(1140 +
861  lmMst1*(747 - 32*lmMt) - 747*lmMt + 16*pow2(lmMst1))*pow2(Mst1) + 2*(
862  323 - 75*lmMgl + 383*lmMst1 - 100*lmMt + 32*pow2(lmMst1))*pow2(MuSUSY))
863  + 5*(2243 + 2052*lmMt - 24*lmMst1*(129 + 10*lmMt) + 612*(pow2(lmMst1) +
864  pow2(lmMt)))*pow4(Mst1)))*pow4(Mt) + 7200*shiftst1*(Dmsqst1*(3 - 2*
865  lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*pow4(s2t)*pow8(Mst1))))))))/(77760.
866  *Tbeta*pow2(Sbeta)*pow4(Msq)*pow6(Mst2)*pow8(Mst1));
867 }
868 
869 /**
870  * @return The off-diagonal (1, 2) = (2, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H9q2'
871  */
872 double H9q2::getS12() const {
873  return (MuSUSY*(-6*xDmst12*pow3(Dmst12)*(threeLoopFlag*pow2(Al4p)*(-3*s2t*(
874  104600 + 104000*lmMst1 - 52480*pow2(lmMst1) + (18675*Dmsqst1)/pow2(Msq)
875  + (64*pow2(Mgl)*(4113 + 90*lmMgl - 1792*lmMst1 + 120*lmMt - 48*pow2(
876  lmMst1) - (54000*Dmsqst1)/pow2(Msq)))/pow2(Mst1))*pow3(Mt) + 720*
877  shiftst3*(16*(7 - 3*lmMst1)*s2t*pow3(Mt) + (-13 + 14*lmMst1)*Mt*pow2(
878  Mst1)*pow3(s2t)) + (240*Mt*((15*Dmsqst1*(48*pow2(Mgl) + (3 - 24*
879  shiftst1 + 12*shiftst2 - 8*lmMst1*(1 - 2*shiftst1 + shiftst2))*pow2(
880  Mst1)) + 2*pow2(Msq)*((200 + 18*lmMgl + 265*lmMst1 + 91*pow2(lmMst1))*
881  pow2(Mgl) - (16 + 60*shiftst1 - 30*shiftst2 + lmMst1*(115 - 120*
882  shiftst1 + 60*shiftst2) + 41*pow2(lmMst1))*pow2(Mst1)))*pow3(s2t) + (4*
883  Mgl*Mt*(-3*pow2(Mst1)*(-30*Dmsqst1*(8*(-3 + lmMgl - lmMst1)*pow2(Mgl) +
884  pow2(Mst1)) + pow2(Msq)*(-2*(227 + 144*lmMgl - 120*lmMst1 - 20*lmMgl*
885  lmMst1 + 27*pow2(lmMgl) + pow2(lmMst1))*pow2(Mgl) + (66 + 77*lmMst1 +
886  41*pow2(lmMst1))*pow2(Mst1)))*pow2(s2t) + 2*pow2(Mt)*(360*Dmsqst1*(31 +
887  lmMst1*(27 - 2*lmMt) - 24*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*
888  pow2(lmMst1))*pow2(Mgl) - 240*Dmsqst1*(11 + 20*lmMst1 - 20*lmMt)*pow2(
889  Mst1) - (-233 + 47*lmMt + lmMst1*(-3 + 2*lmMt) + 15*pow2(lmMst1) + 24*
890  pow2(lmMt))*pow2(Msq)*pow2(Mst1) - 4*pow2(Mgl)*pow2(Msq)*(5845 - 2937*
891  lmMt - 27*(16 + lmMgl + 6*lmMst1 - 9*lmMt)*pow2(lmMgl) + (346 + 183*
892  lmMt)*pow2(lmMst1) + lmMst1*(4753 - 84*lmMt - 72*pow2(lmMt)) - 156*
893  pow2(lmMt) + lmMgl*(-783 - 426*lmMst1*(-1 + lmMt) - 52*lmMt + 273*pow2(
894  lmMst1) + 72*pow2(lmMt)) - 84*pow3(lmMst1)))))/pow4(Mst1)))/pow2(Msq) +
895  (Mt*(8*Mt*MuSUSY*pow2(Msq)*(-12*pow2(Mgl)*(7200*Dmsqst1*pow2(Mst1)*
896  pow2(s2t) + pow2(Msq)*(16*(-323 + 75*lmMgl - 383*lmMst1 + 100*lmMt -
897  32*pow2(lmMst1))*pow2(Mt) + 5*(-425 + 18*lmMgl + 83*lmMst1 + 91*pow2(
898  lmMst1))*pow2(Mst1)*pow2(s2t))) + Mt*s2t*(-16*Mgl*(18000*Dmsqst1 + (-
899  365 + 1049*lmMst1 + 123*pow2(lmMst1))*pow2(Msq))*pow2(Mst1) + 3840*(15*
900  Dmsqst1*(3 - lmMgl + lmMst1) + (-130 + lmMgl*(159 - 20*lmMst1) - 162*
901  lmMst1 + 27*pow2(lmMgl) + pow2(lmMst1))*pow2(Msq))*pow3(Mgl)) - 3*(15*
902  Dmsqst1*(5 - 72*shiftst1 - 48*shiftst2 + 16*lmMst1*(-5 + 3*shiftst1 +
903  2*shiftst2)) - 2*(-605 + 180*shiftst1 + 120*shiftst2 + 348*shiftst3 -
904  8*lmMst1*(-110 + 45*shiftst1 + 30*shiftst2 + 27*shiftst3) + 820*pow2(
905  lmMst1))*pow2(Msq))*pow2(s2t)*pow4(Mst1)) - 45*xDmsqst1*pow2(Dmsqst1)*(
906  -5120*Mt*(2*(3 - lmMgl + lmMst1)*Mt*MuSUSY*s2t + 3*Tbeta*(31 + 27*
907  lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*
908  pow2(lmMst1))*pow2(Mt) + 3*(-3 + lmMgl - lmMst1)*Tbeta*pow2(Mst1)*pow2(
909  s2t))*pow3(Mgl) + 640*Mgl*pow2(Mst1)*(32*s2t*(2*MuSUSY - 9*Mgl*Tbeta)*
910  pow2(Mt) + 3*Mt*(8*Mgl*MuSUSY + 7*Tbeta*pow2(Mst1))*pow2(s2t) + 2*(91 +
911  103*lmMst1 - 103*lmMt)*Tbeta*pow3(Mt) - 6*Mgl*Tbeta*pow2(Mst1)*pow3(
912  s2t)) + s2t*(8*Mt*MuSUSY*s2t*(215 - 72*shiftst1 - 208*shiftst2 + 16*
913  lmMst1*(-5 + 3*shiftst1 + 2*shiftst2)) + Tbeta*(15*(-1389 + 256*
914  shiftst2)*pow2(Mt) - 20*(33 - 128*shiftst1 + 64*shiftst2 - 32*lmMst1*(1
915  - 2*shiftst1 + shiftst2))*pow2(Mst1)*pow2(s2t)))*pow4(Mst1))))/(Tbeta*
916  pow4(Msq)*pow4(Mst1))) + (Mt*(-270*oneLoopFlag*s2t*(-2*Mt*MuSUSY*s2t -
917  8*Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*pow2(s2t)) - (288*Al4p*twoLoopFlag*
918  (5*Mt*pow2(Mst1)*(32*(1 - lmMgl + lmMst1)*Mt*MuSUSY*s2t + 8*Tbeta*(18 +
919  lmMst1*(26 - 6*lmMt) - 13*lmMt + lmMgl*(-13 - 6*lmMst1 + 6*lmMt) + 6*
920  pow2(lmMst1))*pow2(Mt) - 3*(1 + 2*lmMgl - 2*lmMst1)*Tbeta*pow2(Mst1)*
921  pow2(s2t))*pow3(Mgl) + 5*s2t*pow2(Mgl)*((1 + 2*lmMst1)*Mt*MuSUSY*s2t -
922  4*Tbeta*pow2(Mt) - (3 + 2*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t))*pow4(
923  Mst1) + Mgl*Mt*(2*(9 + 2*lmMst1)*Mt*MuSUSY*s2t + 20*(1 + lmMt)*Tbeta*
924  pow2(Mt) + 15*(3 + 2*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t))*pow4(Mst1) -
925  20*Mt*x*(2*(1 - 3*lmMgl + 3*lmMst1)*Mt*MuSUSY*s2t - 6*Tbeta*(-22*lmMst1
926  + lmMgl*(17 + 6*lmMst1 - 6*lmMt) + 5*lmMt + 6*lmMst1*lmMt - 6*pow2(
927  lmMst1))*pow2(Mt) + 3*(-1 + 3*lmMgl - 3*lmMst1)*Tbeta*pow2(Mst1)*pow2(
928  s2t))*pow5(Mgl) + 120*(1 - 3*lmMgl + 3*lmMst1)*s2t*Tbeta*x*pow2(Mt)*
929  pow6(Mgl) - 5*s2t*(4*lmMst1*Mt*MuSUSY*s2t + 2*(-9 + 8*lmMst1)*Tbeta*
930  pow2(Mt) - (1 + 2*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t))*pow6(Mst1)))/
931  pow6(Mst1)))/Tbeta) + ((-9*threeLoopFlag*pow2(Al4p)*(Mt*(1200*xDmsqst1*
932  pow2(Dmsqst1)*pow2(Mst1)*pow2(Mst2)*(-2*Dmst12*Mt*pow2(Mst2)*(Mt*(16*
933  Mgl*(-20*MuSUSY*s2t - 3*(11 + 47*lmMst1 - 47*lmMt)*Mt*Tbeta + 90*Mgl*
934  s2t*Tbeta)*pow2(Mst1) - 64*(2*(3 - lmMgl + lmMst1)*MuSUSY*s2t + 3*Mt*
935  Tbeta*(31 + 27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1
936  + 2*lmMt) + 2*pow2(lmMst1)))*pow3(Mgl)) - s2t*(-16*(-2 + lmMst1)*
937  MuSUSY*s2t*(shiftst1 - shiftst2) + 3*Mt*(-87 - 32*lmMst1*(-1 +
938  shiftst2) + 64*shiftst2)*Tbeta)*pow4(Mst1)) + pow2(Dmst12)*(Mt*(32*Mgl*
939  pow2(Mst1)*(3*Mt*s2t*(2*MuSUSY - 9*Mgl*Tbeta) + (29 - 19*lmMst1 + 19*
940  lmMt)*Tbeta*pow2(Mt) + 3*(2*Mgl*MuSUSY + 5*Tbeta*pow2(Mst1))*pow2(s2t))
941  - 64*(4*(3 - lmMgl + lmMst1)*Mt*MuSUSY*s2t + 6*Tbeta*(31 + 27*lmMst1 -
942  24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(
943  lmMst1))*pow2(Mt) + 3*(-3 + lmMgl - lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t))
944  *pow3(Mgl)) + s2t*(Mt*MuSUSY*s2t*(33 + 16*lmMst1*(-2 + shiftst1 +
945  shiftst2) - 32*(shiftst1 + shiftst2)) + Tbeta*(-321*pow2(Mt) - 24*(-2 +
946  lmMst1)*(shiftst1 - shiftst2)*pow2(Mst1)*pow2(s2t)))*pow4(Mst1)) + 192*
947  Tbeta*pow2(Mt)*(2*Mt*((1 + 7*lmMst1 - 7*lmMt)*Mgl*pow2(Mst1) + (18 -
948  lmMgl + lmMst1)*(1 + lmMst1 - lmMt)*pow3(Mgl)) + (-2 + lmMst1)*s2t*(
949  shiftst1 - shiftst2)*pow4(Mst1))*pow4(Mst2)) - 3*z3*(80*Dmst12*pow2(Mt)
950  *pow4(Mst2)*(-96*xMgl*(121*Mt*MuSUSY + 93*s2t*Tbeta*pow2(Mst1))*pow4(
951  Mgl)*pow4(Msq) + 32*pow2(Mgl)*pow2(Mst1)*((7*Mt*MuSUSY - 6*s2t*Tbeta*
952  pow2(Mst1))*pow4(Msq) - 2*Mgl*(1440*Dmsqst1*Mt*Tbeta*(Dmsqst1*xDmsqst1
953  + pow2(Msq)) + (81*MuSUSY*s2t + 50*Mt*Tbeta)*pow4(Msq))) - 32*Mgl*(
954  1680*Mt*Tbeta*xDmsqst1*pow2(Dmsqst1) + 960*Dmsqst1*Mt*Tbeta*pow2(Msq) +
955  (-7*MuSUSY*s2t + 248*Mt*Tbeta)*pow4(Msq))*pow4(Mst1) + 3*s2t*Tbeta*(
956  525*xDmsqst1*pow2(Dmsqst1) + 280*Dmsqst1*pow2(Msq) - 248*pow4(Msq))*
957  pow6(Mst1)) - 40*Mt*pow2(Dmst12)*pow2(Mst2)*(-24*xMgl*(-372*Mt*s2t*
958  Tbeta*pow2(Mst1) + 968*MuSUSY*pow2(Mt) + 259*MuSUSY*pow2(Mst1)*pow2(
959  s2t))*pow4(Mgl)*pow4(Msq) + 8*pow2(Mgl)*pow2(Mst1)*(-23040*Dmsqst1*Mgl*
960  Tbeta*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mt) + MuSUSY*(7*pow2(Mt) -
961  90*pow2(Mst1)*pow2(s2t))*pow4(Msq) + 4*Mgl*(162*Mt*MuSUSY*s2t + 2192*
962  Tbeta*pow2(Mt) + 243*Tbeta*pow2(Mst1)*pow2(s2t))*pow4(Msq)) + 16*Mgl*
963  Tbeta*(-960*Dmsqst1*(Dmsqst1*xDmsqst1 - 2*pow2(Msq))*pow2(Mt) + (1136*
964  pow2(Mt) - 21*pow2(Mst1)*pow2(s2t))*pow4(Msq))*pow4(Mst1) + 3*s2t*(35*(
965  7*MuSUSY*s2t - 15*Mt*Tbeta)*xDmsqst1*pow2(Dmsqst1) + 70*Dmsqst1*(2*
966  MuSUSY*s2t - 7*Mt*Tbeta)*pow2(Msq) + 48*(MuSUSY*s2t - 11*Mt*Tbeta)*
967  pow4(Msq))*pow6(Mst1)) - xDmst12*pow3(Dmst12)*(Mt*(1920*pow2(Mst1)*
968  pow3(Mgl)*(3840*Dmsqst1*Tbeta*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mt) -
969  (432*Mt*MuSUSY*s2t + 3056*Tbeta*pow2(Mt) + 81*Tbeta*pow2(Mst1)*pow2(
970  s2t))*pow4(Msq)) - 128*Mgl*(4800*Dmsqst1*Tbeta*(5*Dmsqst1*xDmsqst1 + 8*
971  pow2(Msq))*pow2(Mt) + (-6*Mt*MuSUSY*s2t + 440*Tbeta*pow2(Mt) - 105*
972  Tbeta*pow2(Mst1)*pow2(s2t))*pow4(Msq))*pow4(Mst1)) + 15*s2t*(35*
973  Dmsqst1*pow2(Msq)*(8*Mt*MuSUSY*s2t + 31*Tbeta*pow2(Mt) + 16*Tbeta*pow2(
974  Mst1)*pow2(s2t)) + 35*xDmsqst1*pow2(Dmsqst1)*(-40*Mt*MuSUSY*s2t + 159*
975  Tbeta*pow2(Mt) + 28*Tbeta*pow2(Mst1)*pow2(s2t)) + 8*(22*Mt*MuSUSY*s2t +
976  85*Tbeta*pow2(Mt) + 24*Tbeta*pow2(Mst1)*pow2(s2t))*pow4(Msq))*pow6(
977  Mst1) + pow4(Msq)*(480*xMgl*pow4(Mgl)*(-2976*s2t*Tbeta*pow2(Mst1)*pow2(
978  Mt) + 1036*Mt*MuSUSY*pow2(Mst1)*pow2(s2t) + 1936*MuSUSY*pow3(Mt) - 259*
979  Tbeta*pow3(s2t)*pow4(Mst1)) - 320*pow2(Mgl)*(-42*MuSUSY*pow2(Mst1)*
980  pow3(Mt) - Mt*s2t*(45*MuSUSY*s2t + 7*Mt*Tbeta)*pow4(Mst1) + 45*Tbeta*
981  pow3(s2t)*pow6(Mst1)))) + 2560*Mgl*pow3(Mt)*(-20*Tbeta*pow2(Mgl)*pow2(
982  Mst1)*(120*xDmsqst1*pow2(Dmsqst1) + 48*Dmsqst1*pow2(Msq) - 19*pow4(Msq)
983  ) + MuSUSY*(7*Mgl*pow2(Mst1) - 242*xMgl*pow3(Mgl))*pow4(Msq) - 16*
984  Tbeta*(60*xDmsqst1*pow2(Dmsqst1) + 30*Dmsqst1*pow2(Msq) + pow4(Msq))*
985  pow4(Mst1))*pow6(Mst2)) - 320*z2*(3*xDmsqst1*xDmst12*pow2(Dmsqst1)*
986  pow2(Mst1)*pow3(Dmst12)*(30*s2t*pow2(Mgl)*pow2(Mst1)*(-4*Mt*MuSUSY*s2t
987  + 84*Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*pow2(s2t)) - 80*Mgl*Mt*pow2(
988  Mst1)*(7*Mt*MuSUSY*s2t + 3*(39 + 10*lmMst1 - 10*lmMt)*Tbeta*pow2(Mt) +
989  3*Tbeta*pow2(Mst1)*pow2(s2t)) + 160*Mt*(2*Mt*MuSUSY*s2t + 36*(4 +
990  lmMst1 - lmMt)*Tbeta*pow2(Mt) - 3*Tbeta*pow2(Mst1)*pow2(s2t))*pow3(Mgl)
991  + s2t*(4*Mt*MuSUSY*s2t*(-15 + 3*shiftst1 + 2*shiftst2) + 480*Tbeta*
992  pow2(Mt) + 5*(9 - 8*shiftst1 + 4*shiftst2)*Tbeta*pow2(Mst1)*pow2(s2t))*
993  pow4(Mst1)) + 4*xMgl*pow4(Mgl)*(6*Mt*pow2(Dmst12)*pow2(Mst2)*(150*
994  Dmsqst1*Mt*s2t*Tbeta*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mst1) + (146 -
995  27*lmMgl + 27*lmMst1)*Mt*s2t*Tbeta*pow2(Mst1)*pow4(Msq) - MuSUSY*(24*(
996  17 - 4*lmMgl + 4*lmMst1)*pow2(Mt) + (49 - 17*lmMgl + 17*lmMst1)*pow2(
997  Mst1)*pow2(s2t))*pow4(Msq)) - 3*xDmst12*pow3(Dmst12)*(300*Dmsqst1*s2t*
998  Tbeta*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mst1)*pow2(Mt) - pow4(Msq)*(
999  8*(-170 + 27*lmMgl - 27*lmMst1)*s2t*Tbeta*pow2(Mst1)*pow2(Mt) + 4*(49 -
1000  17*lmMgl + 17*lmMst1)*Mt*MuSUSY*pow2(Mst1)*pow2(s2t) + 48*(17 - 4*lmMgl
1001  + 4*lmMst1)*MuSUSY*pow3(Mt) + (-49 + 17*lmMgl - 17*lmMst1)*Tbeta*pow3(
1002  s2t)*pow4(Mst1))) - 12*Dmst12*pow2(Mt)*(75*Dmsqst1*s2t*Tbeta*(Dmsqst1*
1003  xDmsqst1 + pow2(Msq))*pow2(Mst1) - (12*(17 - 4*lmMgl + 4*lmMst1)*Mt*
1004  MuSUSY + (194 - 27*lmMgl + 27*lmMst1)*s2t*Tbeta*pow2(Mst1))*pow4(Msq))*
1005  pow4(Mst2) + 32*(47 - 12*lmMgl + 12*lmMst1)*MuSUSY*pow3(Mt)*pow4(Msq)*
1006  pow6(Mst2)))) + pow2(Mst1)*(64*MuSUSY*pow2(Msq)*pow2(Mst2)*pow2(Mt)*(-
1007  10*Dmst12*pow2(Mst2)*(16*(16 - 9*lmMgl + 41*lmMst1 - 12*lmMt + 26*z2 +
1008  4*pow2(lmMst1))*pow2(Mgl)*pow2(Msq)*pow2(Mt) + Mt*s2t*(4*Mgl*pow2(Msq)*
1009  (2*(-22 - 296*lmMst1 + lmMgl*(308 - 40*lmMst1 - 90*z2) + 113*z2 + 90*
1010  lmMst1*z2 + 54*pow2(lmMgl) + 2*pow2(lmMst1))*pow2(Mgl) - (215 + 118*
1011  lmMst1 + 34*z2 + 41*pow2(lmMst1))*pow2(Mst1)) + 240*Dmsqst1*(Mgl*(-3 +
1012  z2)*pow2(Mst1) + 2*(-3 + lmMgl - lmMst1 + z2)*pow3(Mgl))) - 3*(10*
1013  Dmsqst1*(shiftst1 - shiftst2)*(3 - 2*(lmMst1 + z2)) + (10*shiftst1 -
1014  10*shiftst2 + shiftst3)*(1 - 2*(lmMst1 + z2))*pow2(Msq))*pow2(s2t)*
1015  pow4(Mst1)) + pow2(Dmst12)*(2*pow2(Mgl)*(-450*Dmsqst1*(-4 + z2)*pow2(
1016  Mst1)*pow2(s2t) + pow2(Msq)*(-4*(163 + 15*lmMgl - 27*lmMst1 + 20*lmMt -
1017  11*z2 - 8*pow2(lmMst1))*pow2(Mt) + 5*(200 + 18*lmMgl + 265*lmMst1 +
1018  274*z2 + 91*pow2(lmMst1))*pow2(Mst1)*pow2(s2t))) + Mt*s2t*(40*Mgl*(30*
1019  Dmsqst1*(7 - 4*z2) + (149 + 41*lmMst1 - 30*z2)*pow2(Msq))*pow2(Mst1) -
1020  40*(-120*Dmsqst1*(-3 + lmMgl - lmMst1 + z2) + (-498 - 352*lmMst1 +
1021  lmMgl*(328 - 40*lmMst1 - 90*z2) + 263*z2 + 90*lmMst1*z2 + 54*pow2(
1022  lmMgl) + 2*pow2(lmMst1))*pow2(Msq))*pow3(Mgl)) + 5*(15*Dmsqst1*(3 + 4*
1023  lmMst1*(-2 + shiftst1 + shiftst2) - 6*(shiftst1 + shiftst2) + 4*(-3 +
1024  shiftst1 + shiftst2)*z2) - (32 + 30*(shiftst1 + shiftst2) + 2*lmMst1*(
1025  115 - 30*(shiftst1 + shiftst2) - 9*shiftst3) + 21*shiftst3 + (82 - 60*(
1026  shiftst1 + shiftst2) - 18*shiftst3)*z2 + 82*pow2(lmMst1))*pow2(Msq))*
1027  pow2(s2t)*pow4(Mst1)) + 80*(5 + 18*lmMgl - 74*lmMst1 + 24*lmMt - 52*z2
1028  - 8*pow2(lmMst1))*pow2(Mgl)*pow2(Msq)*pow2(Mt)*pow4(Mst2)) - 160*Tbeta*
1029  (Mt*pow2(Msq)*pow2(Mst2)*(3*Dmst12*Mt*s2t*pow2(Mst1)*(-4*Dmst12*Mgl*
1030  s2t*(60*Dmsqst1*(2*(3 - lmMgl + lmMst1)*pow2(Mgl) + 3*pow2(Mst1)) +
1031  pow2(Msq)*(-4*(-11 + lmMgl*(154 - 20*lmMst1) - 148*lmMst1 + 27*pow2(
1032  lmMgl) + pow2(lmMst1))*pow2(Mgl) + (215 + 118*lmMst1 + 41*pow2(lmMst1))
1033  *pow2(Mst1))) + Mt*(15*Dmsqst1*(Dmst12*(336*pow2(Mgl) + 31*pow2(Mst1))
1034  + 8*(36*pow2(Mgl) + (5 - 4*lmMst1)*pow2(Mst1))*pow2(Mst2)) + 4*pow2(
1035  Msq)*(Dmst12*(2*(522 + 107*lmMst1)*pow2(Mgl) + (-95 - 66*lmMst1 + 82*
1036  pow2(lmMst1))*pow2(Mst1)) + 2*((501 - 18*lmMgl + 504*lmMst1 - 48*lmMt +
1037  107*pow2(lmMst1))*pow2(Mgl) - 2*(-6 + 79*lmMst1 + 41*pow2(lmMst1))*
1038  pow2(Mst1))*pow2(Mst2)))) - 9*s2t*shiftst3*pow2(Msq)*pow4(Mst1)*(-8*
1039  Dmst12*(-3 + 2*lmMst1)*pow2(Mst2)*pow2(Mt) + pow2(Dmst12)*(16*(-2 +
1040  lmMst1)*pow2(Mt) + (1 - 2*lmMst1)*pow2(Mst1)*pow2(s2t)) + 8*(-1 + 2*
1041  lmMst1)*pow2(Mt)*pow4(Mst2)) + 90*s2t*(Dmsqst1*(3 - 2*lmMst1) + (1 - 2*
1042  lmMst1)*pow2(Msq))*pow4(Mst1)*(-8*Dmst12*shiftst2*pow2(Mst2)*pow2(Mt) +
1043  (-shiftst1 + shiftst2)*pow2(Dmst12)*pow2(Mst1)*pow2(s2t) + 8*(shiftst1
1044  - shiftst2)*pow2(Mt)*pow4(Mst2)) + 16*Mgl*pow3(Mt)*(30*Dmsqst1*(pow2(
1045  Dmst12)*(6*(31 + 27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt + lmMgl*(-3 - 2*
1046  lmMst1 + 2*lmMt) + 2*pow2(lmMst1))*pow2(Mgl) + (-13 - 19*lmMst1 + 19*
1047  lmMt)*pow2(Mst1)) - 6*Dmst12*((31 + 27*lmMst1 - 24*lmMt - 2*lmMst1*lmMt
1048  + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))*pow2(Mgl) + (3 + 7*
1049  lmMst1 - 7*lmMt)*pow2(Mst1))*pow2(Mst2) - 12*((7 + 6*lmMst1 - 5*lmMt -
1050  lmMst1*lmMt + lmMgl*(-1 - lmMst1 + lmMt) + pow2(lmMst1))*pow2(Mgl) + (1
1051  + 2*lmMst1 - 2*lmMt)*pow2(Mst1))*pow4(Mst2)) - pow2(Msq)*(pow2(Dmst12)*
1052  (-((-3749 - 2374*lmMst1 + 1770*lmMt - 117*lmMst1*lmMt + lmMgl*(96 -
1053  237*lmMst1 + 157*lmMt) + 108*pow2(lmMgl) + 41*pow2(lmMst1) + 48*pow2(
1054  lmMt))*pow2(Mgl)) + (lmMst1*(297 + 2*lmMt) + 15*pow2(lmMst1) + 3*(47 -
1055  79*lmMt + 8*pow2(lmMt)))*pow2(Mst1)) + Dmst12*pow2(Mst2)*(-((-605 + 6*
1056  lmMst1*(-77 + lmMt) + 714*lmMt + 45*pow2(lmMst1) + 72*pow2(lmMt))*pow2(
1057  Mst1)) + 2*pow2(Mgl)*(2096 - 1167*lmMt - 81*(4 + 2*lmMst1 - 3*lmMt)*
1058  pow2(lmMgl) + 3*(129 + 61*lmMt)*pow2(lmMst1) - 108*pow2(lmMt) - 3*
1059  lmMst1*(-793 + 67*lmMt + 24*pow2(lmMt)) + 3*lmMgl*(-229 + lmMst1*(63 -
1060  142*lmMt) + 35*lmMt + 91*pow2(lmMst1) + 24*pow2(lmMt)) - 27*pow3(lmMgl)
1061  - 84*pow3(lmMst1))) + 2*(-((-233 + 6*lmMst1*(-32 + lmMt) + 348*lmMt +
1062  45*pow2(lmMst1) + 72*pow2(lmMt))*pow2(Mst1)) + pow2(Mgl)*(758 - 570*
1063  lmMt - 81*(3 + 2*lmMst1 - 3*lmMt)*pow2(lmMgl) + 3*(141 + 61*lmMt)*pow2(
1064  lmMst1) - 72*pow2(lmMt) - 6*lmMst1*(-236 + 49*lmMt + 12*pow2(lmMt)) +
1065  3*lmMgl*(-200 + 2*lmMst1 + 76*lmMt - 142*lmMst1*lmMt + 91*pow2(lmMst1)
1066  + 24*pow2(lmMt)) - 27*pow3(lmMgl) - 84*pow3(lmMst1)))*pow4(Mst2))) + 2*
1067  z2*(90*Dmsqst1*(8*pow2(Mt)*(Dmst12*pow2(Mst2)*(4*(7 + 2*lmMst1 - 2*
1068  lmMt)*Mgl*Mt*pow2(Mst1) - 3*s2t*pow2(Mgl)*pow2(Mst1) + 24*(4 + lmMst1 -
1069  lmMt)*Mt*pow3(Mgl) + s2t*(-2 + shiftst2)*pow4(Mst1)) + pow2(Mst1)*(4*(4
1070  + lmMst1 - lmMt)*Mgl*Mt + s2t*(-shiftst1 + shiftst2)*pow2(Mst1))*pow4(
1071  Mst2)) - pow2(Dmst12)*(8*pow3(Mgl)*(-(Mt*pow2(Mst1)*pow2(s2t)) + 24*(4
1072  + lmMst1 - lmMt)*pow3(Mt)) + 8*s2t*pow2(Mt)*(6*pow2(Mgl)*pow2(Mst1) +
1073  pow4(Mst1)) - 4*Mgl*(4*(7 + 2*lmMst1 - 2*lmMt)*pow2(Mst1)*pow3(Mt) +
1074  Mt*pow2(s2t)*pow4(Mst1)) - (shiftst1 - shiftst2)*pow3(s2t)*pow6(Mst1)))
1075  + pow2(Msq)*(pow2(Mt)*(-8*Dmst12*pow2(Mst2)*(-2*(457 + 87*lmMst1 - 87*
1076  lmMt)*Mgl*Mt*pow2(Mst1) - 531*s2t*pow2(Mgl)*pow2(Mst1) + 2*(1627 + 234*
1077  lmMst1 + 174*lmMt)*Mt*pow3(Mgl) - 9*s2t*(-10 + 10*shiftst2 + shiftst3)*
1078  pow4(Mst1)) + 8*pow2(Mst1)*(4*(100 - 3*lmMst1 + 3*lmMt)*Mgl*Mt - 9*s2t*
1079  (10*shiftst1 - 10*shiftst2 + shiftst3)*pow2(Mst1))*pow4(Mst2)) - pow2(
1080  Dmst12)*(4*Mgl*Mt*pow2(Mst1)*(270*Mgl*Mt*s2t - 2*(691 + 216*lmMst1 -
1081  216*lmMt)*pow2(Mt) + 51*pow2(Mst1)*pow2(s2t)) - 12*Mt*((2098 + 548*
1082  lmMst1 - 548*lmMt)*pow2(Mt) + (113 - 90*lmMgl + 90*lmMst1)*pow2(Mst1)*
1083  pow2(s2t))*pow3(Mgl) + 72*s2t*(5 + shiftst3)*pow2(Mt)*pow4(Mst1) - 9*(
1084  10*shiftst1 - 10*shiftst2 + shiftst3)*pow3(s2t)*pow6(Mst1))))) + 16*
1085  pow3(Mgl)*pow3(Mt)*(45*(29 - 10*lmMgl + 10*lmMst1)*Mgl*s2t*xDmsqst1*
1086  xDmst12*xMgl*pow2(Dmsqst1)*pow3(Dmst12) + Mt*z2*(720*Dmsqst1*(4 +
1087  lmMst1 - lmMt)*pow2(Msq) - 4*(1070 + 207*lmMst1 - 3*lmMt)*pow4(Msq))*
1088  pow6(Mst2))))))/(pow4(Msq)*pow6(Mst1)) + 1080*pow2(Mst2)*pow2(Mt)*(3*
1089  Dmst12*oneLoopFlag*s2t*(Dmst12*MuSUSY*s2t + 6*Dmst12*Mt*Tbeta - 12*Mt*
1090  Tbeta*pow2(Mst2)) - (16*Al4p*twoLoopFlag*(pow4(Mst1)*(Dmst12*s2t*(4*
1091  Mgl*Mt*MuSUSY*((1 - 2*lmMgl + 2*lmMst1)*pow2(Mgl) + (2 + lmMst1)*pow2(
1092  Mst1))*pow2(Mst2) + Dmst12*((3 + 2*lmMst1)*MuSUSY*s2t*pow2(Mgl)*pow2(
1093  Mst1) + (2*(3 - 2*lmMgl + 2*lmMst1)*Mt*MuSUSY + 3*(1 - 2*lmMgl + 2*
1094  lmMst1)*s2t*Tbeta*pow2(Mst1))*pow3(Mgl) - (1 + 2*lmMst1)*MuSUSY*s2t*
1095  pow4(Mst1) + Mgl*(2*Mt*MuSUSY*pow2(Mst1) + 3*(2 + lmMst1)*s2t*Tbeta*
1096  pow4(Mst1)))) + Mt*Tbeta*(-3*Dmst12*s2t*pow2(Mst1)*(2*Dmst12*(pow2(Mgl)
1097  + (-1 + 2*lmMst1)*pow2(Mst1)) + 4*((2 + lmMst1)*pow2(Mgl) - (1 + 2*
1098  lmMst1)*pow2(Mst1))*pow2(Mst2)) + Mt*(pow2(Dmst12)*(4*(1 + lmMt)*Mgl*
1099  pow2(Mst1) - 8*(-3 + lmMgl - 2*lmMst1 + lmMt)*pow3(Mgl)) + 3*Dmst12*
1100  pow2(Mst2)*(-4*(1 + lmMt)*Mgl*pow2(Mst1) + 4*(2 + 6*lmMst1 - 3*lmMt -
1101  2*lmMst1*lmMt + lmMgl*(-3 - 2*lmMst1 + 2*lmMt) + 2*pow2(lmMst1))*pow3(
1102  Mgl)) + (-24*(1 + lmMt)*Mgl*pow2(Mst1) - 24*(lmMgl*(1 + lmMst1 - lmMt)
1103  + lmMst1*(-2 + lmMt) + lmMt - pow2(lmMst1))*pow3(Mgl))*pow4(Mst2)))) -
1104  x*(Dmst12*(-1 + 3*lmMgl - 3*lmMst1)*s2t*pow2(Mst1)*(-4*Dmst12*Mt*MuSUSY
1105  + 3*Dmst12*s2t*Tbeta*pow2(Mst1) + 4*Mt*MuSUSY*pow2(Mst2)) + 3*Mt*Tbeta*
1106  (4*Mt*(9 + lmMgl*(20 + 6*lmMst1 - 6*lmMt) + 2*lmMt + lmMst1*(-22 + 6*
1107  lmMt) - 6*pow2(lmMst1))*pow2(Mgl)*pow4(Mst2) + pow2(Mst1)*(4*Dmst12*(-1
1108  + 3*lmMgl - 3*lmMst1)*Mgl*s2t*(Dmst12 - pow2(Mst2)) + 4*Dmst12*Mt*(
1109  lmMst1*(22 - 6*lmMt) - 5*lmMt + lmMgl*(-17 - 6*lmMst1 + 6*lmMt) + 6*
1110  pow2(lmMst1))*(Dmst12 - pow2(Mst2)) + 8*Mt*(2 + lmMgl*(5 + 2*lmMst1 -
1111  2*lmMt) + 2*lmMst1*(-3 + lmMt) + lmMt - 2*pow2(lmMst1))*pow4(Mst2))))*
1112  pow5(Mgl)))/pow8(Mst1)))/Tbeta + (8*Al4p*((5*xMgl*pow2(Mst1)*pow4(Mgl)*
1113  (216*Dmst12*Mt*s2t*twoLoopFlag*pow2(Mst1)*(-2*Dmst12*Mt*((1 - 3*lmMgl +
1114  3*lmMst1)*MuSUSY*s2t + 6*(-3 + 2*lmMgl - 2*lmMst1)*Mt*Tbeta)*pow2(Mst2)
1115  + xDmst12*pow2(Dmst12)*(4*(1 - 3*lmMgl + 3*lmMst1)*Mt*MuSUSY*s2t + 96*(
1116  -1 + lmMgl - lmMst1)*Tbeta*pow2(Mt) + (-1 + 3*lmMgl - 3*lmMst1)*Tbeta*
1117  pow2(Mst1)*pow2(s2t)) + 24*(1 - 2*lmMgl + 2*lmMst1)*Tbeta*pow2(Mt)*
1118  pow4(Mst2)) + (Al4p*threeLoopFlag*(xDmst12*pow2(Msq)*pow3(Dmst12)*(
1119  25920*Dmsqst1*(29 - 10*lmMgl + 10*lmMst1)*s2t*Tbeta*pow2(Mst1)*pow3(Mt)
1120  + Mt*pow2(Msq)*(1728*s2t*Tbeta*(-1797 - 1711*lmMst1 + lmMgl*(1315 -
1121  198*lmMst1 - 96*lmMt) + 144*lmMt + 96*lmMst1*lmMt + 333*pow2(lmMgl) -
1122  103*pow2(lmMst1))*pow2(Mst1)*pow2(Mt) + 4*Mt*MuSUSY*(41653 + 11784*
1123  lmMst1 + 48*lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl) + 540*pow2(
1124  lmMst1))*pow2(Mst1)*pow2(s2t) + 1152*MuSUSY*(262 + 335*lmMst1 - 96*lmMt
1125  - 72*lmMst1*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt) + 120*pow2(
1126  lmMst1))*pow3(Mt) + Tbeta*(-41653 - 11784*lmMst1 - 48*lmMgl*(-313 +
1127  180*lmMst1) + 12636*pow2(lmMgl) - 540*pow2(lmMst1))*pow3(s2t)*pow4(
1128  Mst1))) + pow2(Mst2)*(-432*Dmst12*s2t*Tbeta*pow2(Mst1)*pow3(Mt)*(60*
1129  Dmsqst1*(29 - 10*lmMgl + 10*lmMst1)*(Dmsqst1*xDmsqst1 + pow2(Msq))*(
1130  Dmst12 - pow2(Mst2)) - (Dmst12*(2961 + lmMst1*(1964 - 96*lmMt) - 192*
1131  lmMt + 2*lmMgl*(-683 + 99*lmMst1 + 48*lmMt) - 333*pow2(lmMgl) + 103*
1132  pow2(lmMst1)) + 2*(633 + 1458*lmMst1 - 96*(1 + lmMst1)*lmMt + 2*lmMgl*(
1133  -632 + 99*lmMst1 + 48*lmMt) - 333*pow2(lmMgl) + 103*pow2(lmMst1))*pow2(
1134  Mst2))*pow4(Msq)) + 2*MuSUSY*pow4(Msq)*(-(pow2(Dmst12)*pow2(Mt)*(576*(
1135  262 + 335*lmMst1 - 96*lmMt - 72*lmMst1*lmMt + lmMgl*(-199 - 120*lmMst1
1136  + 72*lmMt) + 120*pow2(lmMst1))*pow2(Mt) + (41653 + 11784*lmMst1 + 48*
1137  lmMgl*(-313 + 180*lmMst1) - 12636*pow2(lmMgl) + 540*pow2(lmMst1))*pow2(
1138  Mst1)*pow2(s2t))) + 576*pow2(Mst2)*(Dmst12*(262 + 335*lmMst1 - 96*lmMt
1139  - 72*lmMst1*lmMt + lmMgl*(-199 - 120*lmMst1 + 72*lmMt) + 120*pow2(
1140  lmMst1)) + 2*(59 + 85*lmMst1 - 24*lmMt - 24*lmMst1*lmMt + lmMgl*(-53 -
1141  40*lmMst1 + 24*lmMt) + 40*pow2(lmMst1))*pow2(Mst2))*pow4(Mt)))))/pow4(
1142  Msq)))/Tbeta + 72*z2*((Al4p*threeLoopFlag*pow4(Mst1)*(-(Mt*xDmst12*
1143  pow2(Msq)*pow3(Dmst12)*(Mt*(8*Mgl*pow2(Mst1)*(75*Dmsqst1*(20*Mt*MuSUSY*
1144  s2t + 48*(7 + 2*lmMst1 - 2*lmMt)*Tbeta*pow2(Mt) - 3*Tbeta*pow2(Mst1)*
1145  pow2(s2t)) - pow2(Msq)*(237*Mt*MuSUSY*s2t + 10*(16 - 9*lmMst1 + 9*lmMt)
1146  *Tbeta*pow2(Mt) + 240*Tbeta*pow2(Mst1)*pow2(s2t))) - 10*(240*Dmsqst1*(
1147  2*Mt*MuSUSY*s2t + 36*(4 + lmMst1 - lmMt)*Tbeta*pow2(Mt) - 3*Tbeta*pow2(
1148  Mst1)*pow2(s2t)) - pow2(Msq)*(32*Mt*((94 - 45*lmMgl + 45*lmMst1)*
1149  MuSUSY*s2t + (380 + 147*lmMst1 - 249*lmMt)*Mt*Tbeta) - 3*(37 + 90*lmMgl
1150  - 90*lmMst1)*Tbeta*pow2(Mst1)*pow2(s2t)))*pow3(Mgl)) - s2t*(30*Dmsqst1*
1151  s2t*(Mt*MuSUSY*(6*shiftst1 + 4*shiftst2) + 5*s2t*(3 - 4*shiftst1 + 2*
1152  shiftst2)*Tbeta*pow2(Mst1)) + pow2(Msq)*(4*Mt*MuSUSY*s2t*(10 + 45*
1153  shiftst1 + 30*shiftst2 + 27*shiftst3) + 40*(-16 + 9*shiftst3)*Tbeta*
1154  pow2(Mt) + 5*(41 - 120*shiftst1 + 60*shiftst2 - 21*shiftst3)*Tbeta*
1155  pow2(Mst1)*pow2(s2t)))*pow4(Mst1) + 2*pow2(Mgl)*(-225*Dmsqst1*s2t*pow2(
1156  Mst1)*(-4*Mt*MuSUSY*s2t + 120*Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*pow2(
1157  s2t)) + pow2(Msq)*(4556*s2t*Tbeta*pow2(Mst1)*pow2(Mt) - 910*Mt*MuSUSY*
1158  pow2(Mst1)*pow2(s2t) + 1992*MuSUSY*pow3(Mt) + 685*Tbeta*pow3(s2t)*pow4(
1159  Mst1))))) + 30*xDmsqst1*pow2(Dmsqst1)*pow2(Mst2)*(5*Dmst12*MuSUSY*s2t*
1160  pow2(Mt)*(Dmst12*(8*Mgl*Mt*pow2(Mst1) + 6*s2t*pow2(Mgl)*pow2(Mst1) -
1161  32*Mt*pow3(Mgl) + s2t*(9 - 2*shiftst1 - 2*shiftst2)*pow4(Mst1)) + 4*
1162  pow2(Mst2)*(10*Mgl*Mt*pow2(Mst1) + 8*Mt*pow3(Mgl) + s2t*(shiftst1 -
1163  shiftst2)*pow4(Mst1))) + 15*Tbeta*(pow3(Mt)*(4*Dmst12*pow2(Mst2)*(2*(45
1164  + 14*lmMst1 - 14*lmMt)*Mgl*Mt*pow2(Mst1) - 15*s2t*pow2(Mgl)*pow2(Mst1)
1165  + 48*(4 + lmMst1 - lmMt)*Mt*pow3(Mgl) + 2*s2t*(-3 + shiftst2)*pow4(
1166  Mst1)) + 8*(2*(13 + 4*lmMst1 - 4*lmMt)*Mgl*Mt*pow2(Mst1) + 4*(19 + 5*
1167  lmMst1 - 5*lmMt)*Mt*pow3(Mgl) + s2t*(-shiftst1 + shiftst2)*pow4(Mst1))*
1168  pow4(Mst2)) - Mt*pow2(Dmst12)*(4*Mgl*((6 + 4*lmMst1 - 4*lmMt)*Mt + 3*
1169  Mgl*s2t)*pow2(Mst1)*pow2(Mt) + 8*pow3(Mgl)*(-(Mt*pow2(Mst1)*pow2(s2t))
1170  + 24*(4 + lmMst1 - lmMt)*pow3(Mt)) + 2*Mt*s2t*(6*Mt - 5*Mgl*s2t)*pow4(
1171  Mst1) - (shiftst1 - shiftst2)*pow3(s2t)*pow6(Mst1))))))/(Tbeta*pow4(
1172  Msq)) - 1440*twoLoopFlag*pow4(Mt)*((Dmst12 + pow2(Mst2))*pow3(Mgl)*
1173  pow4(Mst1)*pow4(Mst2) + x*pow2(Mst1)*pow2(Mst2)*(-3*pow2(Dmst12) + 3*
1174  Dmst12*pow2(Mst2) + 2*pow4(Mst2))*pow5(Mgl) - xDmst12*pow2(Mst1)*pow3(
1175  Dmst12)*(pow2(Mst1)*pow3(Mgl) - 3*x*pow5(Mgl)) + 3*x*pow6(Mst2)*pow7(
1176  Mgl)))))/pow8(Mst1)))/(155520.*pow6(Mst2));
1177 }
1178 
1179 /**
1180  * @return returns the susy log^0 term of Mh^2 @ O(at*as^2) without any log(mu^2) terms normalized to DO (H3m*12/Mt^4/Sbeta^2)
1181  */
1182 double H9q2::calc_coef_at_as2_no_sm_logs_log0() const {
1183 
1184  const double result =
1185  (-432*pow2(log(pow2(Mst1)/pow2(Mgl)))*pow4(Msq)*(-2*Mt*pow2(Dmst12)*pow2(
1186  Mst2)*(-72*Mt*s2t*pow2(Mst1) + 202*Mgl*pow2(Mt) + 111*Mgl*pow2(Mst1)*
1187  pow2(s2t))*pow3(Mgl) + pow3(Dmst12)*pow3(Mgl)*(-48*s2t*pow2(Mst1)*pow2(
1188  Mt) + 111*Mgl*Mt*pow2(Mst1)*pow2(s2t) + 1568*Mgl*pow3(Mt) + 24*pow3(
1189  s2t)*pow4(Mst1)) - 8*Dmst12*pow2(Mgl)*(95*Mt*pow2(Mgl) - 18*Mt*pow2(
1190  Mst1) - 81*Mgl*s2t*pow2(Mst1))*pow2(Mt)*pow4(Mst2) - 8*pow3(Mt)*(-36*
1191  pow2(Mgl)*pow2(Mst1) + 77*pow4(Mgl) + 4*pow4(Mst1))*pow6(Mst2)) - 192*
1192  pow2(Mt)*pow3(log(pow2(Mst1)/pow2(Mgl)))*pow4(Msq)*(-95*Mt*pow2(Dmst12)
1193  *pow2(Mst2)*pow4(Mgl) + 504*Mt*pow3(Dmst12)*pow4(Mgl) + 2*Dmst12*(-157*
1194  Mgl*Mt + 54*s2t*pow2(Mst1))*pow3(Mgl)*pow4(Mst2) + 2*Mt*(-157*pow4(Mgl)
1195  + 18*pow4(Mst1))*pow6(Mst2)) + 5*pow2(Dmsqst1)*(216*Dmst12*pow2(Mt)*(
1196  24*Mt*(-11 + 90*z2 - 112*z3)*pow2(Mgl)*pow2(Mst1) - 64*s2t*(-27 + 43*z2
1197  - 40*z3)*pow2(Mst1)*pow3(Mgl) + 8*Mt*(-307 + 480*z2 - 480*z3)*pow4(Mgl)
1198  + Mt*(146 + 144*z2 - 321*z3)*pow4(Mst1) - 32*Mgl*s2t*(-9 + 30*z2 - 32*
1199  z3)*pow4(Mst1))*pow4(Mst2) - 3*Mt*pow2(Dmst12)*pow2(Mst2)*(4608*Mt*s2t*
1200  (11*z2 - 8*(2 + z3))*pow2(Mst1)*pow3(Mgl) + 576*((-307 + 480*z2 - 480*
1201  z3)*pow2(Mt) + (29 - 10*z2)*pow2(Mst1)*pow2(s2t))*pow4(Mgl) + 576*Mgl*
1202  Mt*s2t*(-33 + 88*z2 - 96*z3)*pow4(Mst1) + (-850 - 8640*z2 + 12123*z3)*
1203  pow2(Mt)*pow4(Mst1) - 1728*pow2(Mgl)*(4*(-5 + 9*z2 - 8*z3)*pow2(Mst1)*
1204  pow2(Mt) + 5*(-2 + z2)*pow2(s2t)*pow4(Mst1)) - 54*(-58 + 64*z2 - 35*z3)
1205  *pow2(s2t)*pow6(Mst1)) + pow3(Dmst12)*(1728*(2*Mt*(29 - 10*z2)*pow2(
1206  Mst1)*pow2(s2t) + (-307 + 480*z2 - 480*z3)*pow3(Mt))*pow4(Mgl) - 4*(
1207  5314 + 6336*z2 - 12123*z3)*pow3(Mt)*pow4(Mst1) - 5184*pow2(Mgl)*((-51 +
1208  162*z2 - 176*z3)*pow2(Mst1)*pow3(Mt) + Mt*(-7 + 4*z2)*pow2(s2t)*pow4(
1209  Mst1)) - 2304*pow3(Mgl)*(-6*s2t*(-59 + 65*z2 - 56*z3)*pow2(Mst1)*pow2(
1210  Mt) + (-3 + z2)*pow3(s2t)*pow4(Mst1)) + 27*Mt*(134 - 192*z2 + 105*z3)*
1211  pow2(s2t)*pow6(Mst1) - 576*Mgl*(12*s2t*(5 - 26*z2 + 32*z3)*pow2(Mt)*
1212  pow4(Mst1) + 5*(-2 + z2)*pow3(s2t)*pow6(Mst1))) + 48*pow3(Mt)*(432*(-1
1213  + 13*z2 - 16*z3)*pow2(Mgl)*pow2(Mst1) + 36*(-179 + 380*z2 - 400*z3)*
1214  pow4(Mgl) + (1118 + 432*z2 - 1539*z3)*pow4(Mst1))*pow6(Mst2)) + 10*
1215  Dmsqst1*pow2(Msq)*(216*Dmst12*pow2(Mt)*(48*Mt*(-3 + 14*z2 - 16*z3)*
1216  pow2(Mgl)*pow2(Mst1) + 32*s2t*(19 - 18*z2 + 16*z3)*pow2(Mst1)*pow3(Mgl)
1217  + 4*Mt*(-307 + 480*z2 - 480*z3)*pow4(Mgl) + Mt*(38 + 48*z2 - 107*z3)*
1218  pow4(Mst1) - 32*Mgl*s2t*(-4 + 9*z2 - 8*z3)*pow4(Mst1))*pow4(Mst2) -
1219  108*Mt*pow2(Dmst12)*pow2(Mst2)*(256*Mt*s2t*(-6 + 9*z2 - 8*z3)*pow2(
1220  Mst1)*pow3(Mgl) + 8*((-307 + 480*z2 - 480*z3)*pow2(Mt) + (29 - 10*z2)*
1221  pow2(Mst1)*pow2(s2t))*pow4(Mgl) + 32*Mgl*Mt*s2t*(-5 + 14*z2 - 16*z3)*
1222  pow4(Mst1) + (-6 - 80*z2 + 107*z3)*pow2(Mt)*pow4(Mst1) - 24*pow2(Mgl)*(
1223  (-19 + 70*z2 - 80*z3)*pow2(Mst1)*pow2(Mt) + 2*(-3 + z2)*pow2(s2t)*pow4(
1224  Mst1)) + 2*(10 - 16*z2 + 7*z3)*pow2(s2t)*pow6(Mst1)) + pow3(Dmst12)*(
1225  864*(2*Mt*(29 - 10*z2)*pow2(Mst1)*pow2(s2t) + (-307 + 480*z2 - 480*z3)*
1226  pow3(Mt))*pow4(Mgl) + (-1814 + 1152*z2 - 567*z3)*pow3(Mt)*pow4(Mst1) +
1227  2592*pow2(Mgl)*(2*(25 - 98*z2 + 112*z3)*pow2(Mst1)*pow3(Mt) + Mt*(-1 +
1228  2*z2)*pow2(s2t)*pow4(Mst1)) - 1152*pow3(Mgl)*(-6*s2t*(-67 + 90*z2 - 80*
1229  z3)*pow2(Mst1)*pow2(Mt) + (-3 + z2)*pow3(s2t)*pow4(Mst1)) + 27*Mt*(18 -
1230  64*z2 + 7*z3)*pow2(s2t)*pow6(Mst1) - 576*Mgl*(3*s2t*(-5 + 4*z2)*pow2(
1231  Mt)*pow4(Mst1) + (-3 + z2)*pow3(s2t)*pow6(Mst1))) + 192*pow3(Mt)*(108*(
1232  -1 + 4*z2 - 4*z3)*pow2(Mgl)*pow2(Mst1) + 9*(-69 + 80*z2 - 80*z3)*pow4(
1233  Mgl) + 4*(23 + 9*z2 - 36*z3)*pow4(Mst1))*pow6(Mst2)) + 8*pow4(Msq)*(-
1234  24*Dmst12*pow2(Mt)*(Mt*(3049 - 2716*z2 + 2388*z3)*pow2(Mgl)*pow2(Mst1)
1235  - 4*s2t*(1043 + 2143*z2 - 570*z3)*pow2(Mst1)*pow3(Mgl) + 2*Mt*(6742 +
1236  4265*z2 - 174*z3)*pow4(Mgl) - 4*Mt*(49 + 90*z2 - 171*z3)*pow4(Mst1) +
1237  4*Mgl*s2t*(-407 + 197*z2 - 24*z3)*pow4(Mst1))*pow4(Mst2) + 3*Mt*pow2(
1238  Dmst12)*pow2(Mst2)*(-48*Mt*s2t*(-1091 + 402*z2 - 480*z3)*pow2(Mst1)*
1239  pow3(Mgl) + 4*((-40279 + 17982*z2 - 26652*z3)*pow2(Mt) - 3*(681 + 776*
1240  z2 - 558*z3)*pow2(Mst1)*pow2(s2t))*pow4(Mgl) - 16*Mgl*Mt*s2t*(-148 +
1241  607*z2 - 696*z3)*pow4(Mst1) + (262 + 1840*z2 - 3423*z3)*pow2(Mt)*pow4(
1242  Mst1) - 4*pow2(Mgl)*(2*(1327 - 4904*z2 + 5280*z3)*pow2(Mst1)*pow2(Mt) +
1243  9*(175 + 118*z2 - 4*z3)*pow2(s2t)*pow4(Mst1)) + 18*(-8 + 40*z2 + 31*z3)
1244  *pow2(s2t)*pow6(Mst1)) + pow3(Dmst12)*(-6*(3*Mt*(1695 - 968*z2 + 558*
1245  z3)*pow2(Mst1)*pow2(s2t) + 4*(-53763 + 9452*z2 - 26304*z3)*pow3(Mt))*
1246  pow4(Mgl) + (-4486 + 816*z2 + 3375*z3)*pow3(Mt)*pow4(Mst1) + 12*pow2(
1247  Mgl)*(4*(-570 + 929*z2 - 840*z3)*pow2(Mst1)*pow3(Mt) + 9*Mt*(1 + 148*z2
1248  - 4*z3)*pow2(s2t)*pow4(Mst1)) - 24*pow3(Mgl)*(s2t*(-2722 + 4704*z2 -
1249  3696*z3)*pow2(Mst1)*pow2(Mt) + (-22 + 113*z2 - 243*z3)*pow3(s2t)*pow4(
1250  Mst1)) - 18*Mt*(-119 + 60*z2 + 192*z3)*pow2(s2t)*pow6(Mst1) + 12*Mgl*(
1251  s2t*(446 - 768*z2 + 624*z3)*pow2(Mt)*pow4(Mst1) + (215 + 34*z2 - 21*z3)
1252  *pow3(s2t)*pow6(Mst1))) - 12*pow3(Mt)*(4*(1909 - 496*z2 + 228*z3)*pow2(
1253  Mgl)*pow2(Mst1) + 2*(6677 + 10414*z2 - 3948*z3)*pow4(Mgl) + (-4125 -
1254  892*z2 + 4740*z3)*pow4(Mst1))*pow6(Mst2)) + 96*log(pow2(Mst1)/pow2(Mgl)
1255  )*(60*pow2(Dmsqst1)*pow3(Mgl)*(-15*Mt*pow2(Dmst12)*pow2(Mst2)*(2*Mt*
1256  s2t*pow2(Mst1) + 3*Mgl*pow2(Mt) - Mgl*pow2(Mst1)*pow2(s2t)) + pow3(
1257  Dmst12)*(78*s2t*pow2(Mst1)*pow2(Mt) - 30*Mgl*Mt*pow2(Mst1)*pow2(s2t) +
1258  45*Mgl*pow3(Mt) - 2*pow3(s2t)*pow4(Mst1)) + 9*Dmst12*(5*Mgl*Mt - 2*s2t*
1259  pow2(Mst1))*pow2(Mt)*pow4(Mst2) + 15*Mgl*pow3(Mt)*pow6(Mst2)) + 60*
1260  Dmsqst1*pow2(Msq)*pow3(Mgl)*(-3*Mt*pow2(Dmst12)*pow2(Mst2)*(4*Mt*s2t*
1261  pow2(Mst1) + 15*Mgl*pow2(Mt) - 5*Mgl*pow2(Mst1)*pow2(s2t)) + pow3(
1262  Dmst12)*(60*s2t*pow2(Mst1)*pow2(Mt) - 30*Mgl*Mt*pow2(Mst1)*pow2(s2t) +
1263  45*Mgl*pow3(Mt) - 2*pow3(s2t)*pow4(Mst1)) + 9*Dmst12*(5*Mgl*Mt - 4*s2t*
1264  pow2(Mst1))*pow2(Mt)*pow4(Mst2) + 30*Mgl*pow3(Mt)*pow6(Mst2)) + pow4(
1265  Msq)*(2*Mt*pow2(Dmst12)*pow2(Mgl)*pow2(Mst2)*(-102*Mgl*Mt*s2t*pow2(
1266  Mst1) + 2*pow2(Mgl)*((737 + 162*z2)*pow2(Mt) + 3*(328 + 27*z2)*pow2(
1267  Mst1)*pow2(s2t)) + 27*pow2(s2t)*pow4(Mst1)) - pow2(Mgl)*pow3(Dmst12)*(
1268  Mt*pow2(Mgl)*(32*(434 + 81*z2)*pow2(Mt) + 3*(605 + 54*z2)*pow2(Mst1)*
1269  pow2(s2t)) + 54*Mt*pow2(s2t)*pow4(Mst1) - 2*Mgl*(67*s2t*pow2(Mst1)*
1270  pow2(Mt) + 2*(-154 + 45*z2)*pow3(s2t)*pow4(Mst1))) + 12*Dmst12*pow2(
1271  Mgl)*(18*Mt*(37 + 6*z2)*pow2(Mgl) + 9*Mt*pow2(Mst1) - 476*Mgl*s2t*pow2(
1272  Mst1))*pow2(Mt)*pow4(Mst2) + 12*pow3(Mt)*(18*pow2(Mgl)*pow2(Mst1) + 3*(
1273  215 + 36*z2)*pow4(Mgl) - 52*pow4(Mst1))*pow6(Mst2))))/(108.*pow3(Mt)*
1274  pow4(Msq)*pow4(Mst1)*pow6(Mst2));
1275 
1276  return result;
1277 }
1278 
1279 /**
1280  * @return returns the susy log^1 term of Mh^2 @ O(at*as^2) without any log(mu^2) terms normalized to DO (H3m*12/Mt^4/Sbeta^2)
1281  */
1282 double H9q2::calc_coef_at_as2_no_sm_logs_log1() const {
1283 
1284  const double result =
1285  (8*(-6*pow2(Mt)*pow2(log(pow2(Mst1)/pow2(Mgl)))*pow4(Msq)*(-203*Mt*pow2(
1286  Dmst12)*pow2(Mst2)*pow4(Mgl) + 936*Mt*pow3(Dmst12)*pow4(Mgl) + 2*
1287  Dmst12*(-265*Mgl*Mt + 162*s2t*pow2(Mst1))*pow3(Mgl)*pow4(Mst2) + 2*Mt*(
1288  -265*pow4(Mgl) + 18*pow4(Mst1))*pow6(Mst2)) + 10*Dmsqst1*Mt*pow2(Msq)*(
1289  18*Dmst12*Mt*(6*Mt*(-7 + 4*z2)*pow2(Mgl)*pow2(Mst1) - 8*s2t*(-5 + 2*z2)
1290  *pow2(Mst1)*pow3(Mgl) + Mt*(-119 + 60*z2)*pow4(Mgl) + 4*Mt*(-2 + z2)*
1291  pow4(Mst1) - 8*Mgl*s2t*(-2 + z2)*pow4(Mst1))*pow4(Mst2) + 2*pow3(
1292  Dmst12)*(-27*(-47 + 28*z2)*pow2(Mgl)*pow2(Mst1)*pow2(Mt) + 72*Mt*s2t*(-
1293  19 + 10*z2)*pow2(Mst1)*pow3(Mgl) + 9*(-119 + 60*z2)*pow2(Mt)*pow4(Mgl)
1294  + 6*Mgl*Mt*s2t*pow4(Mst1) + pow2(Mt)*pow4(Mst1) - 9*pow2(s2t)*pow6(
1295  Mst1)) - 9*pow2(Dmst12)*pow2(Mst2)*(3*(33 - 20*z2)*pow2(Mgl)*pow2(Mst1)
1296  *pow2(Mt) + 16*Mt*s2t*(-7 + 4*z2)*pow2(Mst1)*pow3(Mgl) + 2*(-119 + 60*
1297  z2)*pow2(Mt)*pow4(Mgl) + 8*Mgl*Mt*s2t*(-3 + 2*z2)*pow4(Mst1) + (7 - 4*
1298  z2)*pow2(Mt)*pow4(Mst1) - 2*pow2(s2t)*pow6(Mst1)) + 6*pow2(Mt)*(36*(-2
1299  + z2)*pow2(Mgl)*pow2(Mst1) + 3*(-49 + 20*z2)*pow4(Mgl) + 2*(-17 + 6*z2)
1300  *pow4(Mst1))*pow6(Mst2)) + 5*Mt*pow2(Dmsqst1)*(2*pow3(Dmst12)*(-27*(-75
1301  + 44*z2)*pow2(Mgl)*pow2(Mst1)*pow2(Mt) + 144*Mt*s2t*(-15 + 7*z2)*pow2(
1302  Mst1)*pow3(Mgl) + 18*(-119 + 60*z2)*pow2(Mt)*pow4(Mgl) + 24*Mgl*Mt*s2t*
1303  (-19 + 12*z2)*pow4(Mst1) - 2*((-61 + 36*z2)*pow2(Mt) + 9*pow2(Mst1)*
1304  pow2(s2t))*pow4(Mst1)) + 18*Dmst12*Mt*(3*Mt*(-47 + 28*z2)*pow2(Mgl)*
1305  pow2(Mst1) - 16*s2t*(-9 + 5*z2)*pow2(Mst1)*pow3(Mgl) + 2*Mt*(-119 + 60*
1306  z2)*pow4(Mgl) - 8*Mgl*s2t*(-7 + 4*z2)*pow4(Mst1) + Mt*(-23 + 12*z2)*
1307  pow4(Mst1))*pow4(Mst2) - 3*pow2(Dmst12)*pow2(Mst2)*(-36*(-7 + 4*z2)*
1308  pow2(Mgl)*pow2(Mst1)*pow2(Mt) + 96*Mt*s2t*(-3 + z2)*pow2(Mst1)*pow3(
1309  Mgl) + 12*(-119 + 60*z2)*pow2(Mt)*pow4(Mgl) + 12*Mgl*Mt*s2t*(-19 + 12*
1310  z2)*pow4(Mst1) + (61 - 36*z2)*pow2(Mt)*pow4(Mst1) - 12*pow2(s2t)*pow6(
1311  Mst1)) + 6*pow2(Mt)*(36*(-7 + 4*z2)*pow2(Mgl)*pow2(Mst1) + (-537 + 300*
1312  z2)*pow4(Mgl) + 4*(-22 + 9*z2)*pow4(Mst1))*pow6(Mst2)) + pow4(Msq)*(Mt*
1313  pow2(Dmst12)*pow2(Mst2)*(-6*pow2(Mgl)*pow2(Mst1)*((659 - 440*z2)*pow2(
1314  Mt) + 243*pow2(Mst1)*pow2(s2t)) - 72*Mt*s2t*(-99 + 20*z2)*pow2(Mst1)*
1315  pow3(Mgl) + ((-23651 + 5436*z2)*pow2(Mt) - 582*pow2(Mst1)*pow2(s2t))*
1316  pow4(Mgl) + 12*Mgl*Mt*s2t*(25 - 62*z2)*pow4(Mst1) + 6*((-145 + 46*z2)*
1317  pow2(Mt) + 79*pow2(Mst1)*pow2(s2t))*pow4(Mst1)) + pow3(Dmst12)*(3*Mt*
1318  pow2(Mgl)*pow2(Mst1)*((-498 + 280*z2)*pow2(Mt) + 379*pow2(Mst1)*pow2(
1319  s2t)) - 24*s2t*pow2(Mst1)*((-86 + 77*z2)*pow2(Mt) + pow2(Mst1)*pow2(
1320  s2t))*pow3(Mgl) + (-315*Mt*pow2(Mst1)*pow2(s2t) + (63184 - 7104*z2)*
1321  pow3(Mt))*pow4(Mgl) + 2*Mgl*s2t*((442 - 60*z2)*pow2(Mt) + 59*pow2(Mst1)
1322  *pow2(s2t))*pow4(Mst1) - 3*Mt*(4*(-43 + 8*z2)*pow2(Mt) + 125*pow2(Mst1)
1323  *pow2(s2t))*pow4(Mst1)) - 6*Dmst12*pow2(Mt)*(4*Mt*(103 - 50*z2)*pow2(
1324  Mgl)*pow2(Mst1) - 4*s2t*(283 + 138*z2)*pow2(Mst1)*pow3(Mgl) + Mt*(2647
1325  + 628*z2)*pow4(Mgl) + 12*Mt*(8 - 5*z2)*pow4(Mst1) - 4*Mgl*s2t*(65 + 2*
1326  z2)*pow4(Mst1))*pow4(Mst2) - 12*pow3(Mt)*(2*(91 - 10*z2)*pow2(Mgl)*
1327  pow2(Mst1) + (589 + 464*z2)*pow4(Mgl) - 4*(-35 + 15*z2 + 3*z3)*pow4(
1328  Mst1))*pow6(Mst2)) + 2*log(pow2(Mst1)/pow2(Mgl))*(180*Dmsqst1*pow2(Msq)
1329  *pow2(Mt)*pow3(Mgl)*(-5*Mgl*Mt*pow2(Dmst12)*pow2(Mst2) + (5*Mgl*Mt + 4*
1330  s2t*pow2(Mst1))*pow3(Dmst12) + Dmst12*(5*Mgl*Mt - 4*s2t*pow2(Mst1))*
1331  pow4(Mst2) + 5*Mgl*Mt*pow6(Mst2)) + 90*pow2(Dmsqst1)*pow2(Mt)*pow3(Mgl)
1332  *(-2*pow2(Dmst12)*(5*Mgl*Mt + 2*s2t*pow2(Mst1))*pow2(Mst2) + 2*(5*Mgl*
1333  Mt + 6*s2t*pow2(Mst1))*pow3(Dmst12) + 2*Dmst12*(5*Mgl*Mt - 2*s2t*pow2(
1334  Mst1))*pow4(Mst2) + 5*Mgl*Mt*pow6(Mst2)) + pow4(Msq)*(6*Mt*pow2(Dmst12)
1335  *pow2(Mst2)*(14*Mt*s2t*pow2(Mst1) + 67*Mgl*pow2(Mt) + 117*Mgl*pow2(
1336  Mst1)*pow2(s2t))*pow3(Mgl) - pow3(Dmst12)*pow3(Mgl)*(42*s2t*pow2(Mst1)*
1337  pow2(Mt) + 351*Mgl*(Mt*pow2(Mst1)*pow2(s2t) + 8*pow3(Mt)) + 68*pow3(
1338  s2t)*pow4(Mst1)) + 12*Dmst12*pow2(Mgl)*(167*Mt*pow2(Mgl) - 54*Mt*pow2(
1339  Mst1) - 170*Mgl*s2t*pow2(Mst1))*pow2(Mt)*pow4(Mst2) + 36*pow3(Mt)*(-36*
1340  pow2(Mgl)*pow2(Mst1) + 59*pow4(Mgl) + 4*pow4(Mst1))*pow6(Mst2)))))/(9.*
1341  pow3(Mt)*pow4(Msq)*pow4(Mst1)*pow6(Mst2));
1342 
1343  return result;
1344 }
1345 
1346 /**
1347  * @return returns the susy log^2 term of Mh^2 @ O(at*as^2) without any log(mu^2) terms normalized to DO (H3m*12/Mt^4/Sbeta^2)
1348  */
1349 double H9q2::calc_coef_at_as2_no_sm_logs_log2() const {
1350 
1351  const double result =
1352  (8*(Mt*pow2(Dmst12)*pow2(Mst2)*pow4(Msq)*(528*Mt*s2t*pow2(Mst1)*pow3(Mgl)
1353  + (-1966*pow2(Mt) + 96*pow2(Mst1)*pow2(s2t))*pow4(Mgl) + 180*Mgl*Mt*
1354  s2t*pow4(Mst1) + 3*(51*pow2(Mt) + 82*pow2(Mst1)*pow2(s2t))*pow4(Mst1) +
1355  pow2(Mgl)*(64*pow2(Mst1)*pow2(Mt) - 321*pow2(s2t)*pow4(Mst1))) - 6*
1356  Dmst12*pow2(Mt)*pow4(Msq)*(43*Mt*pow2(Mgl)*pow2(Mst1) - 248*s2t*pow2(
1357  Mst1)*pow3(Mgl) + 550*Mt*pow4(Mgl) + 51*Mt*pow4(Mst1) + 60*Mgl*s2t*
1358  pow4(Mst1))*pow4(Mst2) - pow3(Dmst12)*pow4(Msq)*((48*Mt*pow2(Mst1)*
1359  pow2(s2t) - 7232*pow3(Mt))*pow4(Mgl) + 102*pow3(Mt)*pow4(Mst1) + pow2(
1360  Mgl)*(32*pow2(Mst1)*pow3(Mt) - 321*Mt*pow2(s2t)*pow4(Mst1)) + 16*pow3(
1361  Mgl)*(11*s2t*pow2(Mst1)*pow2(Mt) + 2*pow3(s2t)*pow4(Mst1)) + 369*Mt*
1362  pow2(s2t)*pow6(Mst1) + Mgl*(120*s2t*pow2(Mt)*pow4(Mst1) - 41*pow3(s2t)*
1363  pow6(Mst1))) + 6*pow3(Mt)*(-86*pow2(Mgl)*pow2(Mst1)*pow4(Msq) - 380*
1364  pow4(Mgl)*pow4(Msq) + (15*pow2(Dmsqst1) + 30*Dmsqst1*pow2(Msq) + 98*
1365  pow4(Msq))*pow4(Mst1))*pow6(Mst2) - 12*log(pow2(Mst1)/pow2(Mgl))*pow2(
1366  Mt)*pow4(Msq)*(-53*Mt*pow2(Dmst12)*pow2(Mst2)*pow4(Mgl) + 276*Mt*pow3(
1367  Dmst12)*pow4(Mgl) + 2*Dmst12*(-85*Mgl*Mt + 44*s2t*pow2(Mst1))*pow3(Mgl)
1368  *pow4(Mst2) + 2*Mt*(-85*pow4(Mgl) + 9*pow4(Mst1))*pow6(Mst2))))/(9.*
1369  pow3(Mt)*pow4(Msq)*pow4(Mst1)*pow6(Mst2));
1370 
1371  return result;
1372 }
1373 
1374 /**
1375  * @return returns the susy log^3 term of Mh^2 @ O(at*as^2) without any log(mu^2) terms normalized to DO (H3m*12/Mt^4/Sbeta^2)
1376  */
1377 double H9q2::calc_coef_at_as2_no_sm_logs_log3() const {
1378 
1379  const double result =
1380  -298.6666666666667;
1381 
1382  return result;
1383 }
1384 
1385 } // namespace hierarchies
1386 } // namespace himalaya
himalaya::hierarchies::ExpansionDepth::Dmsqst1
truncate the expansion depth in the difference of the stop/sbottom 1 mass and the average squark mass...
Definition: HierarchyFlags.hpp:44
himalaya::hierarchies::H9q2::calc_coef_at_as2_no_sm_logs_log3
double calc_coef_at_as2_no_sm_logs_log3() const
Definition: H9q2.cpp:1377
himalaya::hierarchies::H9q2::getS2
double getS2() const
Definition: H9q2.cpp:187
himalaya::hierarchies::ExpansionFlags_t
std::array< int, ExpansionDepth::NUMBER_OF_EXPANSIONS > ExpansionFlags_t
Definition: HierarchyTypes.hpp:17
himalaya
Definition: H3.cpp:14
himalaya::hierarchies::H9q2::threeLoopFlag
int threeLoopFlag
loop flags
Definition: H9q2.hpp:75
himalaya::hierarchies::ExpansionDepth::Dmst12
truncate the expansion depth in the difference of the stop/sbottom masses by one order ...
Definition: HierarchyFlags.hpp:45
himalaya::hierarchies::H9q2::getS12
double getS12() const
Definition: H9q2.cpp:872
himalaya::hierarchies::ExpansionDepth::Mgl
truncate the expansion depth in the gluino mass by one order
Definition: HierarchyFlags.hpp:52
himalaya::hierarchies::ExpansionDepth::Msq
truncate the expansion depth in the average squark mass by one order
Definition: HierarchyFlags.hpp:48
himalaya::hierarchies::H9q2::calc_coef_at_as2_no_sm_logs_log2
double calc_coef_at_as2_no_sm_logs_log2() const
Definition: H9q2.cpp:1349
himalaya::hierarchies::H9q2::calc_coef_at_as2_no_sm_logs_log0
double calc_coef_at_as2_no_sm_logs_log0() const
Definition: H9q2.cpp:1182
himalaya::hierarchies::ExpansionDepth::threeLoop
truncate the two loop expansion at the three loop expansion depth
Definition: HierarchyFlags.hpp:41
himalaya::hierarchies::H9q2::calc_coef_at_as2_no_sm_logs_log1
double calc_coef_at_as2_no_sm_logs_log1() const
Definition: H9q2.cpp:1282
himalaya::hierarchies::ExpansionDepth::Dmglst1
truncate the expansion depth in the difference of stop/sbottom 1 mass and the gluino mass by one orde...
Definition: HierarchyFlags.hpp:43
himalaya::hierarchies::H9q2::getS1
double getS1() const
Definition: H9q2.cpp:77
himalaya::log
Complex< T > log(const Complex< T > &z_) noexcept
Definition: complex.hpp:45
himalaya::hierarchies::H9q2::H9q2
H9q2(const ExpansionFlags_t &expansionDepth, double Al4p, double beta, double Dmst12, double Dmsqst1, double lmMt, double lmMgl, double lmMst1, double Mgl, double Mt, double Mst1, double Mst2, double Msq, double MuSUSY, double s2t, int mdrFlag, int oneLoopFlag, int twoLoopFlag, int threeLoopFlag)
Definition: H9q2.cpp:39
himalaya::hierarchies::H9q2::Al4p
double Al4p
common variables
Definition: H9q2.hpp:73
himalaya::hierarchies::H9q2::xDmsqst1
int xDmsqst1
MDR and truncation flags.
Definition: H9q2.hpp:74
Generated by   doxygen 1.8.13