Himalaya
H3.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/H3.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  * Constuctor
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 Dmglst1 a double Mgl - Mst1
23  * @param Dmst12 a double Mst1^2 - Mst2^2
24  * @param Dmsqst1 a double Msq^2 - Mst1^2
25  * @param lmMt a double log((renormalization scale / Mt)^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 MuSUSY a double mu parameter
32  * @param s2t a double 2 times the sine of the stop/sbottom quark mixing angle
33  * @param mdrFlag an int 0 for DR and 1 for MDR scheme
34  * @param oneLoopFlag an int flag to consider the one-loop expansion terms
35  * @param twoLoopFlag an int flag to consider the two-loop expansion terms
36  * @param threeLoopFlag an int flag to consider the three-loop expansion terms
37  */
38 H3::H3 (const ExpansionFlags_t& expansionDepth, double Al4p, double beta,
39  double Dmglst1, double Dmst12, double Dmsqst1, double lmMt, double lmMst1,
40  double Mgl, double Mt, double Mst1, double Mst2, double Msq, double MuSUSY,
41  double s2t,
42  int mdrFlag, int oneLoopFlag, int twoLoopFlag, int threeLoopFlag) {
43  // abbrev for tan(beta) and sin(beta)
44  Tbeta = tan(beta);
45  Sbeta = sin(beta);
46  this -> Dmglst1 = Dmglst1;
47  this -> Dmst12 = Dmst12;
48  this -> Dmsqst1 = Dmsqst1;
49  this -> lmMst1 = lmMst1;
50  this -> Mgl = Mgl;
51  this -> Mt = Mt;
52  this -> Mst1 = Mst1;
53  this -> Mst2 = Mst2;
54  this -> Msq = Msq;
55  this -> MuSUSY = MuSUSY;
56  this -> s2t = s2t;
57  this -> lmMt = lmMt;
58  this -> oneLoopFlag = oneLoopFlag;
59  this -> twoLoopFlag = twoLoopFlag;
60  this -> threeLoopFlag = threeLoopFlag;
61  this -> Al4p = Al4p;
62  // mdr flags, indicates if one wants to shift the dr stop mass to the mdr stop mass
63  shiftst1 = mdrFlag;
64  shiftst2 = mdrFlag;
65  shiftst3 = mdrFlag;
66  // expansion flags
67  xDmglst1 = expansionDepth.at(ExpansionDepth::Dmglst1);
68  xDmst12 = expansionDepth.at(ExpansionDepth::Dmglst1);
69  xDmsqst1 = expansionDepth.at(ExpansionDepth::Dmsqst1);
70 }
71 
72 /**
73  * @return The diagonal (1, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3'
74  */
75 double H3::getS1() const {
76  return -(pow2(Mt)*pow2(MuSUSY)*(-4*xDmst12*pow3(Dmst12)*(4050*oneLoopFlag*pow2(
77  s2t) + (8*Al4p*(450*Mgl*Mst1*s2t*twoLoopFlag*(4*Dmglst1*Mgl*((5 - 6*
78  lmMst1)*Mt + 4*(1 - 3*lmMst1)*Mst1*s2t) - 4*(Mt + 6*lmMst1*Mt + 3*(-5 +
79  6*lmMst1)*Mst1*s2t)*pow2(Dmglst1) + (-4*(5 + 6*lmMst1)*Mt + (1 + 6*
80  lmMst1)*Mst1*s2t)*pow2(Mgl)) + (xDmglst1*(-360*Mst1*s2t*(3*(9 + 10*
81  lmMst1)*Mt + 10*(-13 + 12*lmMst1)*Mst1*s2t)*twoLoopFlag*pow2(Msq) +
82  Al4p*threeLoopFlag*(30000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(2*Mst1*Mt*
83  s2t*(37824007 + 770520*lmMst1 - 131400*pow2(lmMst1)) + (59957863 +
84  480000*lmMst1 - 26880*lmMt - 230400*pow2(lmMst1))*pow2(Mt) + 15*(-
85  3044017 - 27472*lmMst1 + 48480*pow2(lmMst1))*pow2(Mst1)*pow2(s2t))))*
86  pow3(Dmglst1))/pow2(Msq)))/(pow2(Mst1)*pow3(Mgl)) + threeLoopFlag*pow2(
87  Al4p)*((Mt*(Mt*(72*pow2(Dmglst1)*(3891491 + 27200*lmMst1 - 960*lmMt -
88  19200*pow2(lmMst1)) + 200*Dmglst1*Mgl*(403559 + 384*(lmMst1 + lmMt) -
89  4608*pow2(lmMst1)) + 15*(-1763661 + 47104*lmMst1 - 5120*lmMt + 24576*
90  pow2(lmMst1))*pow2(Mgl)) + (240*Mst1*s2t*(5*pow2(Mgl)*(840*Dmsqst1 + (-
91  36863 + 80*lmMst1 + 552*pow2(lmMst1))*pow2(Msq)) + 10*pow2(Dmglst1)*(
92  100*Dmsqst1 + (-32829 + 1852*lmMst1 + 660*pow2(lmMst1))*pow2(Msq)) +
93  Dmglst1*Mgl*(1000*Dmsqst1 + (-1282471 + 7264*lmMst1 + 18120*pow2(
94  lmMst1))*pow2(Msq))))/pow2(Msq)))/(pow2(Mgl)*pow2(Mst1)) + 15*(350605 +
95  4320*shiftst1 + 2880*shiftst2 + 8352*shiftst3 - 96*lmMst1*(-115 + 90*
96  shiftst1 + 60*shiftst2 + 54*shiftst3) - 2160*pow2(lmMst1) + (40*
97  Dmglst1*(-84209 - 1264*lmMst1 + 240*pow2(lmMst1)))/Mgl + (8*pow2(
98  Dmglst1)*(-1732531 - 16896*lmMst1 + 24840*pow2(lmMst1)))/pow2(Mgl) + (
99  2400*Dmsqst1*(7 - 24*lmMst1*(-1 + shiftst2) + 36*shiftst2))/pow2(Msq))*
100  pow2(s2t) + (12000*s2t*xDmsqst1*pow2(Dmsqst1)*(20*Dmglst1*Mt*(Dmglst1*
101  Mgl + xDmglst1*pow2(Dmglst1) + pow2(Mgl)) + 3*(28*Mt + Mst1*s2t*(7 + 6*
102  shiftst1 - 24*lmMst1*(-1 + shiftst2) + 30*shiftst2))*pow3(Mgl)))/(Mst1*
103  pow3(Mgl)*pow4(Msq)))) + 8*pow2(Mst2)*((-1800*Al4p*Dmst12*s2t*
104  twoLoopFlag*(pow2(Mgl)*(Dmst12*(4*Mt - 3*Mst1*s2t) + 12*(1 + 2*lmMst1)*
105  Mt*pow2(Mst2)) + pow2(Dmglst1)*(6*Dmst12*(4*Mt + (-5 + 6*lmMst1)*Mst1*
106  s2t) + 4*(-11 + 6*lmMst1)*Mt*pow2(Mst2)) + 4*Dmglst1*Mgl*(2*Dmst12*(-1
107  + 3*lmMst1)*Mst1*s2t + (-5 + 6*lmMst1)*Mt*pow2(Mst2))))/(Mst1*pow2(Mgl)
108  ) + 2025*oneLoopFlag*pow2(Dmst12)*pow2(s2t) - (2*Al4p*xDmglst1*pow3(
109  Dmglst1)*(720*Dmst12*Mst1*s2t*twoLoopFlag*pow2(Msq)*(52*Dmst12*Mt + 5*
110  Dmst12*(-13 + 12*lmMst1)*Mst1*s2t + (-77 + 30*lmMst1)*Mt*pow2(Mst2)) +
111  Al4p*threeLoopFlag*(2*Dmst12*Mt*(30000*Dmsqst1*Mst1*s2t + (4*Mst1*s2t*(
112  31025111 + 290880*lmMst1 - 251100*pow2(lmMst1)) + Mt*(59957863 +
113  480000*lmMst1 - 26880*lmMt - 230400*pow2(lmMst1)))*pow2(Msq))*pow2(
114  Mst2) + pow2(Dmst12)*(-60000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(2*Mst1*
115  Mt*s2t*(-99874229 - 1352280*lmMst1 + 633600*pow2(lmMst1)) + 2*(-
116  59957863 - 480000*lmMst1 + 26880*lmMt + 230400*pow2(lmMst1))*pow2(Mt) +
117  15*(3044017 + 27472*lmMst1 - 48480*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))
118  - 24*(-3877891 - 46400*lmMst1 + 960*lmMt + 19200*pow2(lmMst1))*pow2(
119  Msq)*pow2(Mt)*pow4(Mst2))))/(pow2(Msq)*pow2(Mst1)*pow3(Mgl))) + (
120  threeLoopFlag*pow2(Al4p)*(-480*Mgl*pow2(Mst1)*pow2(s2t)*((Dmst12*pow2(
121  Msq)*pow2(Mst2)*(300*Dmsqst1*pow2(Mgl)*(Dmst12*(-7 - 18*shiftst1 - 18*
122  shiftst2 + 12*lmMst1*(-2 + shiftst1 + shiftst2)) - 12*(-3 + 2*lmMst1)*(
123  shiftst1 - shiftst2)*pow2(Mst2)) + pow2(Msq)*(Dmst12*pow2(Dmglst1)*(
124  1732531 + 16896*lmMst1 - 24840*pow2(lmMst1)) + 5*Dmglst1*Dmst12*Mgl*(
125  84209 + 1264*lmMst1 - 240*pow2(lmMst1)) + 10*pow2(Mgl)*(Dmst12*(1429 -
126  180*shiftst1 - 180*shiftst2 + lmMst1*(-454 + 360*shiftst1 + 360*
127  shiftst2) + 24*pow2(lmMst1)) - 360*(-1 + 2*lmMst1)*(shiftst1 -
128  shiftst2)*pow2(Mst2)))))/2. + pow2(Mgl)*(90*Dmst12*shiftst3*pow2(Mst2)*
129  (Dmst12*(-7 + 6*lmMst1) + 2*(1 - 2*lmMst1)*pow2(Mst2))*pow4(Msq) + 36*
130  z2*(-2*xDmst12*pow3(Dmst12)*(100*Dmsqst1*shiftst2*(Dmsqst1*xDmsqst1 +
131  pow2(Msq)) + (15*shiftst1 + 10*shiftst2 + 9*shiftst3)*pow4(Msq)) +
132  Dmst12*pow2(Mst2)*(5*shiftst3*(3*Dmst12 - 2*pow2(Mst2))*pow4(Msq) - 50*
133  (-(Dmst12*(shiftst1 + shiftst2)) + 2*(shiftst1 - shiftst2)*pow2(Mst2))*
134  (xDmsqst1*pow2(Dmsqst1) + Dmsqst1*pow2(Msq) + pow4(Msq)))))) + pow2(
135  Mst2)*(24000*Dmst12*Mst1*s2t*xDmsqst1*pow2(Dmsqst1)*(40*Dmglst1*Mt*(
136  Dmglst1*Mgl + xDmglst1*pow2(Dmglst1) + pow2(Mgl))*(Dmst12 - pow2(Mst2))
137  - 3*(-(Dmst12*(56*Mt + Mst1*s2t*(7 + 30*shiftst1 + 6*shiftst2 - 12*
138  lmMst1*(-2 + shiftst1 + shiftst2)))) + 8*(7*Mt - 3*(-2 + lmMst1)*Mst1*
139  s2t*(shiftst1 - shiftst2))*pow2(Mst2))*pow3(Mgl)) + 2*Mgl*Mt*pow2(Msq)*
140  (240*Dmst12*Mst1*s2t*(5*pow2(Mgl)*(1680*Dmsqst1*(Dmst12 - pow2(Mst2)) +
141  pow2(Msq)*(3*Dmst12*(-10473 + 40*lmMst1 + 256*pow2(lmMst1)) - 8*(1361 +
142  10*lmMst1 + 54*pow2(lmMst1))*pow2(Mst2))) + Dmglst1*Mgl*(2000*Dmsqst1*(
143  Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(-949861 + 1944*lmMst1 +
144  11520*pow2(lmMst1)) + 20*(-33261 + 532*lmMst1 + 660*pow2(lmMst1))*pow2(
145  Mst2))) + pow2(Dmglst1)*(2000*Dmsqst1*(Dmst12 - pow2(Mst2)) + pow2(Msq)
146  *(Dmst12*(958501 + 24456*lmMst1 - 11520*pow2(lmMst1)) + 2*(-1286791 -
147  5936*lmMst1 + 18120*pow2(lmMst1))*pow2(Mst2)))) + Mt*pow2(Msq)*(16*
148  pow2(Dmglst1)*(-9*Dmst12*(-3891491 - 27200*lmMst1 + 960*lmMt + 19200*
149  pow2(lmMst1))*(Dmst12 - pow2(Mst2)) - 50*(345581 + 4896*lmMst1 + 96*
150  lmMt - 3456*pow2(lmMst1))*pow4(Mst2)) + 400*Dmglst1*Mgl*(Dmst12*(403559
151  + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))*(Dmst12 - pow2(Mst2)) - 24*(
152  9631 + 16*lmMst1 + 48*lmMt - 192*pow2(lmMst1))*pow4(Mst2)) - 15*pow2(
153  Mgl)*(pow2(Dmst12)*(852541 + 9216*lmMst1 - 10240*lmMt + 6144*pow2(
154  lmMst1)) + 160*Dmst12*(11389 - 704*lmMst1 + 192*lmMt - 384*pow2(lmMst1)
155  )*pow2(Mst2) + 1920*(349 - 56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow4(
156  Mst2))))) - 225*z3*(-2*xDmst12*pow3(Dmst12)*((-16*Mgl*Mst1*Mt*s2t*(
157  142987*Dmglst1*Mgl + 37582*pow2(Dmglst1) + 20297*pow2(Mgl)) + 3*Mgl*(
158  197112*Dmglst1*Mgl + 687960*pow2(Dmglst1) - 65963*pow2(Mgl))*pow2(Mt) +
159  8*xDmglst1*(557078*Mst1*Mt*s2t + 442053*pow2(Mt) - 349745*pow2(Mst1)*
160  pow2(s2t))*pow3(Dmglst1))*pow4(Msq) - Mgl*pow2(Mst1)*pow2(s2t)*(-10080*
161  Dmsqst1*pow2(Mgl)*(Dmsqst1*xDmsqst1 + pow2(Msq)) + (403880*Dmglst1*Mgl
162  + 1600920*pow2(Dmglst1) - 37669*pow2(Mgl))*pow4(Msq))) + pow2(Mst2)*(-
163  8*(Mgl*pow2(Dmst12)*pow2(Mst1)*pow2(s2t)*(-1260*Dmsqst1*pow2(Mgl)*(
164  Dmsqst1*xDmsqst1 + pow2(Msq)) + (50485*Dmglst1*Mgl + 200115*pow2(
165  Dmglst1) + 2574*pow2(Mgl))*pow4(Msq)) + xDmglst1*pow3(Dmglst1)*pow4(
166  Msq)*(18*Dmst12*Mt*(49117*Mt + 102024*Mst1*s2t)*pow2(Mst2) + pow2(
167  Dmst12)*(-1475294*Mst1*Mt*s2t - 884106*pow2(Mt) + 349745*pow2(Mst1)*
168  pow2(s2t)) + 687960*pow2(Mt)*pow4(Mst2))) + Mgl*Mt*pow4(Msq)*(16*
169  Dmst12*Mst1*s2t*(5*Dmst12*(-21081*Dmglst1*Mgl + 21081*pow2(Dmglst1) -
170  3457*pow2(Mgl)) - 2*(37582*Dmglst1*Mgl + 142987*pow2(Dmglst1) + 3012*
171  pow2(Mgl))*pow2(Mst2)) + 3*Mt*(48*Dmglst1*Mgl*(8213*pow2(Dmst12) -
172  8213*Dmst12*pow2(Mst2) - 4664*pow4(Mst2)) + 432*pow2(Dmglst1)*(3185*
173  pow2(Dmst12) - 3185*Dmst12*pow2(Mst2) - 1566*pow4(Mst2)) - pow2(Mgl)*(
174  31963*pow2(Dmst12) + 68000*Dmst12*pow2(Mst2) + 24064*pow4(Mst2))))))))/
175  (pow2(Mst1)*pow3(Mgl)*pow4(Msq))))/(777600.*pow6(Mst2));
176 }
177 
178 /**
179  * @return The diagonal (2, 2) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3'
180  */
181 double H3::getS2() const {
182  return (pow2(Mst2)*(2*Al4p*Mgl*Mt*pow2(Msq)*(-70560*Mst1*Mt*twoLoopFlag*pow2(
183  Msq)*(2*MuSUSY*pow2(Sbeta)*(-3*pow2(Dmglst1)*(2*Dmst12*Mt*((477 - 330*
184  lmMst1 + 30*lmMt)*Mt + 50*(13 - 18*lmMst1)*Mst1*s2t)*pow2(Mst2) - pow2(
185  Dmst12)*(70*Mst1*Mt*s2t + 8*(92 - 30*lmMst1 + 5*lmMt)*pow2(Mt) + 25*(11
186  - 6*lmMst1)*pow2(Mst1)*pow2(s2t)) + 400*(4 - 3*lmMst1)*pow2(Mt)*pow4(
187  Mst2)) + 25*pow2(Mgl)*(4*Dmst12*Mt*(2*(5 + 6*lmMst1 + 3*lmMt)*Mt - 9*(1
188  + 2*lmMst1)*Mst1*s2t)*pow2(Mst2) - pow2(Dmst12)*(24*(1 - 3*lmMst1)*
189  Mst1*Mt*s2t + 2*(5 + 6*(lmMst1 + lmMt))*pow2(Mt) + 9*(1 + 2*lmMst1)*
190  pow2(Mst1)*pow2(s2t)) + 72*(1 + lmMst1 + lmMt)*pow2(Mt)*pow4(Mst2)) -
191  Dmglst1*Mgl*(-300*Dmst12*Mt*((-3 + 6*lmMst1)*Mt + 2*(-1 + 6*lmMst1)*
192  Mst1*s2t)*pow2(Mst2) + 3*pow2(Dmst12)*(-100*Mst1*Mt*s2t + 2*(-41 + 90*
193  lmMst1 + 10*lmMt)*pow2(Mt) + 25*(-5 + 6*lmMst1)*pow2(Mst1)*pow2(s2t)) -
194  200*(-7 + 15*lmMst1 + 3*lmMt)*pow2(Mt)*pow4(Mst2))) + Tbeta*(-25*
195  Dmst12*s2t*(4*Dmglst1*Mgl*((-5 + 6*lmMst1)*Mt*pow2(Mst2)*pow2(MuSUSY) -
196  2*Dmst12*Mst1*s2t*((-1 + 3*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*
197  (-1 + 6*lmMst1)*pow2(Mst1)*pow2(Sbeta))) + 2*pow2(Dmglst1)*(2*(-11 + 6*
198  lmMst1)*Mt*pow2(Mst2)*pow2(MuSUSY) + 3*Dmst12*(4*Mt*pow2(MuSUSY) - (-5
199  + 6*lmMst1)*Mst1*s2t*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 2*(13 - 18*
200  lmMst1)*s2t*pow2(Sbeta)*pow3(Mst1))) + pow2(Mgl)*(12*(1 + 2*lmMst1)*Mt*
201  pow2(Mst2)*pow2(MuSUSY) + Dmst12*(4*Mt*pow2(MuSUSY) + 3*Mst1*s2t*pow2(
202  MuSUSY)*(-1 + pow2(Sbeta)) + 36*(1 + 2*lmMst1)*s2t*pow2(Sbeta)*pow3(
203  Mst1)))) + pow2(Sbeta)*(4*Dmst12*Mt*s2t*(-25*pow2(Mgl)*(Dmst12*(4*(1 +
204  3*lmMst1 + 6*lmMt)*pow2(Mst1) - pow2(MuSUSY)) - 3*pow2(Mst2)*(6*(1 + 2*
205  (lmMst1 + lmMt))*pow2(Mst1) + (1 + 2*lmMst1)*pow2(MuSUSY))) + pow2(
206  Dmglst1)*(-6*Dmst12*((-169 + 135*lmMst1 + 15*lmMt)*pow2(Mst1) - 25*
207  pow2(MuSUSY)) + 25*pow2(Mst2)*(24*(-4 + 3*lmMst1)*pow2(Mst1) + (-11 +
208  6*lmMst1)*pow2(MuSUSY))) + 25*Dmglst1*Mgl*(-4*Dmst12*(-4 + 6*lmMst1 +
209  3*lmMt)*pow2(Mst1) + pow2(Mst2)*(2*(-17 + 30*lmMst1 + 6*lmMt)*pow2(
210  Mst1) + (-5 + 6*lmMst1)*pow2(MuSUSY)))) - 2*Mst1*pow2(Mt)*(4*Dmglst1*
211  Mgl*((11 + 60*lmMst1 - 60*lmMt)*pow2(Dmst12) + 25*Dmst12*(1 + 30*lmMst1
212  + 6*lmMt)*pow2(Mst2) + 100*(1 + 12*lmMst1 + 6*lmMt)*pow4(Mst2)) + 2*
213  pow2(Dmglst1)*((-83 + 120*lmMst1 - 120*lmMt)*pow2(Dmst12) + Dmst12*(-
214  1381 + 2490*lmMst1 + 210*lmMt)*pow2(Mst2) + 200*(-11 + 21*lmMst1 + 6*
215  lmMt)*pow4(Mst2)) + 25*pow2(Mgl)*((5 + 42*lmMst1 + 30*lmMt)*pow2(
216  Dmst12) - 8*Dmst12*(4 + 3*lmMst1 + 6*lmMt)*pow2(Mst2) - 72*(1 + lmMst1
217  - lmMt + 2*lmMst1*lmMt + pow2(lmMst1) - 3*pow2(lmMt))*pow4(Mst2)))))) +
218  Al4p*threeLoopFlag*(784*Mst1*MuSUSY*pow2(Sbeta)*(225*Mt*pow2(Dmst12)*(
219  pow2(Dmglst1)*(1000*Dmsqst1 + (1286791 + 5936*lmMst1 - 18120*pow2(
220  lmMst1))*pow2(Msq)) + 10*Dmglst1*Mgl*(100*Dmsqst1 + (33261 - 532*lmMst1
221  - 660*pow2(lmMst1))*pow2(Msq)) + 20*pow2(Mgl)*(210*Dmsqst1 + (1361 +
222  10*lmMst1 + 54*pow2(lmMst1))*pow2(Msq)))*pow2(Mst1)*pow2(s2t) + pow3(
223  Mt)*(2*Dmglst1*(1200*Dmglst1*Dmsqst1*(Dmst12*(557 + 120*lmMst1 - 120*
224  lmMt)*(Dmst12 - pow2(Mst2)) - 16*(4 + 15*lmMst1 - 15*lmMt)*pow4(Mst2))
225  + 1200*Dmsqst1*Mgl*(Dmst12*(557 + 120*lmMst1 - 120*lmMt)*(Dmst12 -
226  pow2(Mst2)) - 50*(26 + 3*lmMst1 - 3*lmMt)*pow4(Mst2)) + Dmglst1*pow2(
227  Msq)*(-3*pow2(Dmst12)*(-129193181 - 401100*lmMt + 100*lmMst1*(-7351 +
228  1164*lmMt) + 1336800*pow2(lmMst1) + 28800*pow2(lmMt)) + 2*Dmst12*(-
229  327941741 - 686700*lmMt + 1080*lmMst1*(-44 + 445*lmMt) + 3531600*pow2(
230  lmMst1) + 64800*pow2(lmMt))*pow2(Mst2) + 80*(-3777727 - 24255*lmMt +
231  486*lmMst1*(118 + 45*lmMt) + 52380*pow2(lmMst1))*pow4(Mst2)) - Mgl*
232  pow2(Msq)*(-3*pow2(Dmst12)*(-39474953 + lmMst1*(52780 - 87600*lmMt) -
233  3300*lmMt + 319200*pow2(lmMst1) + 14400*pow2(lmMt)) + 40*Dmst12*(
234  3746977 + 4005*lmMt - 18*lmMst1*(2711 + 1215*lmMt) - 52380*pow2(lmMst1)
235  )*pow2(Mst2) + 6000*(9961 + 66*lmMt - 10*lmMst1*(67 + 24*lmMt) + 42*
236  pow2(lmMst1) + 72*pow2(lmMt))*pow4(Mst2))) + 15*pow2(Mgl)*(4000*
237  Dmsqst1*(Dmst12*(65 - 6*lmMst1 + 6*lmMt)*(Dmst12 - pow2(Mst2)) + 6*(-20
238  + 3*lmMst1 - 3*lmMt)*pow4(Mst2)) - pow2(Msq)*(-(pow2(Dmst12)*(-2284899
239  + 49840*lmMt - 32*lmMst1*(-1793 + 555*lmMt) + 87360*pow2(lmMst1) +
240  28800*pow2(lmMt))) + 200*Dmst12*(11697 + 448*lmMst1 + 330*lmMt - 372*
241  lmMst1*lmMt + 408*pow2(lmMst1) + 288*pow2(lmMt))*pow2(Mst2) + 1600*(434
242  - 83*lmMst1 + 174*lmMt - 66*lmMst1*lmMt + 183*pow2(lmMst1) + 108*pow2(
243  lmMt))*pow4(Mst2)))) + 450*Mst1*s2t*(Dmst12*(pow2(Dmglst1)*(19160*
244  Dmsqst1*(Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(-586073 - 9268*
245  lmMst1 - 448*lmMt + 19200*pow2(lmMst1)) + 2*(13289 - 916*lmMst1 + 80*
246  lmMt - 60*pow2(lmMst1))*pow2(Mst2))) + 5*pow2(Mgl)*(5*Dmsqst1*(161*
247  Dmst12 + 24*(1 + 12*lmMst1)*pow2(Mst2)) + pow2(Msq)*(Dmst12*(2071 +
248  296*lmMst1 - 96*lmMt - 600*pow2(lmMst1)) + 8*(631 - lmMst1 + 36*lmMt +
249  21*pow2(lmMst1))*pow2(Mst2))) + 4*Dmglst1*Mgl*(4700*Dmsqst1*(Dmst12 -
250  pow2(Mst2)) - pow2(Msq)*(Dmst12*(17539 + 574*lmMst1 + 160*lmMt - 1920*
251  pow2(lmMst1)) + 5*(-4539 + 596*lmMst1 - 48*lmMt + 516*pow2(lmMst1))*
252  pow2(Mst2))))*pow2(Mt) + 45*pow2(Mgl)*(shiftst3*pow2(Msq)*(-8*Dmst12*(-
253  3 + 2*lmMst1)*pow2(Mst2)*pow2(Mt) + pow2(Dmst12)*(16*(-2 + lmMst1)*
254  pow2(Mt) + (1 - 2*lmMst1)*pow2(Mst1)*pow2(s2t)) + 8*(-1 + 2*lmMst1)*
255  pow2(Mt)*pow4(Mst2)) + 10*(Dmsqst1*(-3 + 2*lmMst1) + (-1 + 2*lmMst1)*
256  pow2(Msq))*(-8*Dmst12*shiftst2*pow2(Mst2)*pow2(Mt) + (-shiftst1 +
257  shiftst2)*pow2(Dmst12)*pow2(Mst1)*pow2(s2t) + 8*(shiftst1 - shiftst2)*
258  pow2(Mt)*pow4(Mst2))))) + Mt*Tbeta*(245*(pow2(MuSUSY)*(8*pow2(Dmglst1)*
259  (6*Dmst12*Mt*(10000*Dmsqst1*Mst1*s2t + (3*Mt*(3891491 + 27200*lmMst1 -
260  960*lmMt - 19200*pow2(lmMst1)) + 10*Mst1*s2t*(1286791 + 5936*lmMst1 -
261  18120*pow2(lmMst1)))*pow2(Msq))*pow2(Mst2) - 3*pow2(Dmst12)*(20000*
262  Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(10*Mst1*Mt*s2t*(958501 + 24456*lmMst1
263  - 11520*pow2(lmMst1)) + 6*(3891491 + 27200*lmMst1 - 960*lmMt - 19200*
264  pow2(lmMst1))*pow2(Mt) + 5*(-1732531 - 16896*lmMst1 + 24840*pow2(
265  lmMst1))*pow2(Mst1)*pow2(s2t))) + 100*(345581 + 4896*lmMst1 + 96*lmMt -
266  3456*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow4(Mst2)) + 40*Dmglst1*Mgl*(10*
267  Dmst12*Mt*(1200*Dmsqst1*Mst1*s2t + (Mt*(403559 + 384*(lmMst1 + lmMt) -
268  4608*pow2(lmMst1)) - 12*Mst1*s2t*(-33261 + 532*lmMst1 + 660*pow2(
269  lmMst1)))*pow2(Msq))*pow2(Mst2) - pow2(Dmst12)*(12000*Dmsqst1*Mst1*Mt*
270  s2t + pow2(Msq)*(6*Mst1*Mt*s2t*(-949861 + 1944*lmMst1 + 11520*pow2(
271  lmMst1)) + 10*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))*pow2(
272  Mt) + 15*(-84209 - 1264*lmMst1 + 240*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)
273  )) + 240*(9631 + 16*lmMst1 + 48*lmMt - 192*pow2(lmMst1))*pow2(Msq)*
274  pow2(Mt)*pow4(Mst2)) + 15*pow2(Mgl)*(160*Dmst12*Mt*(840*Dmsqst1*Mst1*
275  s2t + (Mt*(11389 - 704*lmMst1 + 192*lmMt - 384*pow2(lmMst1)) + 4*Mst1*
276  s2t*(1361 + 10*lmMst1 + 54*pow2(lmMst1)))*pow2(Msq))*pow2(Mst2) + pow2(
277  Dmst12)*(-2400*Dmsqst1*Mst1*s2t*(56*Mt + (7 + 24*lmMst1)*Mst1*s2t) +
278  pow2(Msq)*(-240*Mst1*Mt*s2t*(-10473 + 40*lmMst1 + 256*pow2(lmMst1)) + (
279  852541 + 9216*lmMst1 - 10240*lmMt + 6144*pow2(lmMst1))*pow2(Mt) + 80*(
280  1429 - 454*lmMst1 + 24*pow2(lmMst1))*pow2(Mst1)*pow2(s2t))) + 1920*(349
281  - 56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow4(Mst2))
282  ) + 21600*pow2(Mgl)*pow2(Mst1)*(-(shiftst3*pow2(Msq)*(-2*Dmst12*pow2(
283  Mst2)*((-1 + 2*lmMst1)*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 12*(
284  (3 - 2*lmMst1)*pow2(Mt) + (-1 + 2*lmMst1)*pow2(Mst1)*pow2(s2t))*pow2(
285  Sbeta)) + pow2(Dmst12)*((-7 + 6*lmMst1)*pow2(MuSUSY)*pow2(s2t)*(-1 +
286  pow2(Sbeta)) + 6*(-8*(-2 + lmMst1)*pow2(Mt) + (-15 + 14*lmMst1)*pow2(
287  Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(1 - 2*lmMst1)*pow2(Mt)*pow2(Sbeta)*
288  pow4(Mst2))) - 10*(Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*(
289  pow2(Dmst12)*pow2(s2t)*(-((shiftst1 + shiftst2)*pow2(MuSUSY)*(-1 +
290  pow2(Sbeta))) + 6*(-3*shiftst1 + shiftst2)*pow2(Mst1)*pow2(Sbeta)) + 2*
291  Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)*pow2(s2t)*(-1 +
292  pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 - shiftst2)*pow2(Mst1)
293  *pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)*pow2(Mt)*pow2(
294  Sbeta)*pow4(Mst2)))) + pow2(Sbeta)*(784*Dmst12*Mst1*Mt*s2t*(375*pow2(
295  Mgl)*(-80*Dmsqst1*(Dmst12*((-98 + 24*lmMst1 - 24*lmMt)*pow2(Mst1) - 21*
296  pow2(MuSUSY)) + 3*pow2(Mst2)*((74 - 12*lmMst1 + 12*lmMt)*pow2(Mst1) +
297  7*pow2(MuSUSY))) + pow2(Msq)*(-8*pow2(Mst2)*(8*(347 - 50*lmMst1 + 66*
298  lmMt - 66*lmMst1*lmMt + 183*pow2(lmMst1) + 108*pow2(lmMt))*pow2(Mst1) +
299  (1361 + 10*lmMst1 + 54*pow2(lmMst1))*pow2(MuSUSY)) + Dmst12*(16*(-4378
300  - 517*lmMst1 + 243*lmMt - 78*lmMst1*lmMt + 528*pow2(lmMst1) + 288*pow2(
301  lmMt))*pow2(Mst1) + 3*(-10473 + 40*lmMst1 + 256*pow2(lmMst1))*pow2(
302  MuSUSY)))) - pow2(Dmglst1)*(1200*Dmsqst1*(Dmst12*((866 - 240*lmMst1 +
303  240*lmMt)*pow2(Mst1) - 125*pow2(MuSUSY)) + pow2(Mst2)*(16*(23 + 30*
304  lmMst1 - 30*lmMt)*pow2(Mst1) + 125*pow2(MuSUSY))) + pow2(Msq)*(10*pow2(
305  Mst2)*(8*(7531199 + 48510*lmMt - 486*lmMst1*(191 + 90*lmMt) - 104760*
306  pow2(lmMst1))*pow2(Mst1) + 15*(1286791 + 5936*lmMst1 - 18120*pow2(
307  lmMst1))*pow2(MuSUSY)) + Dmst12*(-4*(-176974411 + 218700*lmMt - 540*
308  lmMst1*(3971 + 730*lmMt) + 1436400*pow2(lmMst1) + 64800*pow2(lmMt))*
309  pow2(Mst1) + 75*(-958501 - 24456*lmMst1 + 11520*pow2(lmMst1))*pow2(
310  MuSUSY)))) - 5*Dmglst1*Mgl*(240*Dmsqst1*(25*pow2(Mst2)*(2*(55 + 6*
311  lmMst1 - 6*lmMt)*pow2(Mst1) + 5*pow2(MuSUSY)) - Dmst12*(4*(379 + 15*
312  lmMst1 - 15*lmMt)*pow2(Mst1) + 125*pow2(MuSUSY))) + pow2(Msq)*(300*
313  pow2(Mst2)*(16*(4964 - 3*lmMt - 5*lmMst1*(55 + 24*lmMt) + 21*pow2(
314  lmMst1) + 36*pow2(lmMt))*pow2(Mst1) + (33261 - 532*lmMst1 - 660*pow2(
315  lmMst1))*pow2(MuSUSY)) + Dmst12*(8*(4511549 + 9810*lmMt + 6*lmMst1*(
316  14879 + 4710*lmMt) - 117360*pow2(lmMst1) - 21600*pow2(lmMt))*pow2(Mst1)
317  - 15*(-949861 + 1944*lmMst1 + 11520*pow2(lmMst1))*pow2(MuSUSY))))) -
318  29400*pow2(Dmst12)*pow2(Mst1)*(pow2(Dmglst1)*(114960*Dmsqst1*pow2(Mst1)
319  + pow2(Msq)*(12*(-13249 + 916*lmMst1 - 80*lmMt + 60*pow2(lmMst1))*pow2(
320  Mst1) + (1732531 + 16896*lmMst1 - 24840*pow2(lmMst1))*pow2(MuSUSY))) +
321  5*Dmglst1*Mgl*(22560*Dmsqst1*pow2(Mst1) + pow2(Msq)*(24*(-4515 + 596*
322  lmMst1 - 48*lmMt + 516*pow2(lmMst1))*pow2(Mst1) + (84209 + 1264*lmMst1
323  - 240*pow2(lmMst1))*pow2(MuSUSY))) - 10*pow2(Mgl)*(30*Dmsqst1*(12*(1 +
324  12*lmMst1)*pow2(Mst1) + (7 + 24*lmMst1)*pow2(MuSUSY)) + pow2(Msq)*(24*(
325  613 - lmMst1 + 36*lmMt + 21*pow2(lmMst1))*pow2(Mst1) + (-1429 + 454*
326  lmMst1 - 24*pow2(lmMst1))*pow2(MuSUSY))))*pow2(s2t) - pow2(Mt)*(16*
327  pow2(Mst1)*(pow2(Dmst12)*(980*pow2(Mgl)*(15*Dmsqst1*(4561 + 510*lmMst1
328  - 510*lmMt) + (93973 - 61305*lmMt + 54*lmMst1*(2337 + 1255*lmMt) -
329  48420*pow2(lmMst1) - 58050*pow2(lmMt))*pow2(Msq)) + 98*Dmglst1*Mgl*(
330  10800*Dmsqst1*(423 - 20*lmMst1 + 20*lmMt) + (-22574599 - 21060*lmMt +
331  60*lmMst1*(6677 + 8880*lmMt) + 1598400*pow2(lmMst1) + 172800*pow2(lmMt)
332  )*pow2(Msq)) + 3*pow2(Dmglst1)*(11760*Dmsqst1*(13147 - 90*lmMst1 + 90*
333  lmMt) + (-5259064991 - 32822300*lmMt + 60*lmMst1*(-104761 + 256200*
334  lmMt) + 148327200*pow2(lmMst1) + 5644800*pow2(lmMt))*pow2(Msq))) - 196*
335  Dmst12*(10*Dmglst1*Mgl*(540*Dmsqst1*(423 - 20*lmMst1 + 20*lmMt) + (-
336  5136871 - 29790*lmMt + 12*lmMst1*(8942 + 4305*lmMt) + 86040*pow2(
337  lmMst1) + 21600*pow2(lmMt))*pow2(Msq)) + pow2(Dmglst1)*(180*Dmsqst1*(
338  13147 - 90*lmMst1 + 90*lmMt) + (-125277461 - 1327590*lmMt + 60*lmMst1*(
339  -2303 + 5955*lmMt) + 153000*pow2(lmMst1) + 151200*pow2(lmMt))*pow2(Msq)
340  ) - 750*pow2(Mgl)*(5*Dmsqst1*(79 + 204*lmMst1 + 12*lmMt) + (-4*lmMst1*(
341  131 + 264*lmMt) + 264*pow2(lmMst1) + 21*(783 + 14*lmMt + 30*pow2(lmMt))
342  )*pow2(Msq)))*pow2(Mst2) + 7840*(-15*Dmsqst1*(50*Dmglst1*(146 - 15*
343  lmMst1 + 15*lmMt)*Mgl + (8014 - 435*lmMst1 + 435*lmMt)*pow2(Dmglst1) +
344  225*(14 + 6*lmMt - 6*lmMst1*(3 + lmMt) + 3*(pow2(lmMst1) + pow2(lmMt)))
345  *pow2(Mgl)) + pow2(Msq)*(-75*Dmglst1*Mgl*(-15571 + 852*lmMt + lmMst1*(-
346  508 + 48*lmMt) + 978*pow2(lmMst1) + 612*pow2(lmMt)) - pow2(Dmglst1)*(-
347  3223748 + 2205*lmMt + 6*lmMst1*(-16808 + 4605*lmMt) + 79695*pow2(
348  lmMst1) + 33750*pow2(lmMt)) + 75*pow2(Mgl)*(623 - 555*lmMt - 3*(221 +
349  129*lmMt)*pow2(lmMst1) - 486*pow2(lmMt) + lmMst1*(1066 + 843*lmMt +
350  756*pow2(lmMt)) + 252*pow3(lmMst1) - 621*pow3(lmMt))))*pow4(Mst2)) -
351  245*pow2(Msq)*pow2(MuSUSY)*(16*pow2(Dmglst1)*(-9*Dmst12*(-3891491 -
352  27200*lmMst1 + 960*lmMt + 19200*pow2(lmMst1))*(Dmst12 - pow2(Mst2)) -
353  50*(345581 + 4896*lmMst1 + 96*lmMt - 3456*pow2(lmMst1))*pow4(Mst2)) +
354  400*Dmglst1*Mgl*(Dmst12*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(
355  lmMst1))*(Dmst12 - pow2(Mst2)) - 24*(9631 + 16*lmMst1 + 48*lmMt - 192*
356  pow2(lmMst1))*pow4(Mst2)) - 15*pow2(Mgl)*(pow2(Dmst12)*(852541 + 9216*
357  lmMst1 - 10240*lmMt + 6144*pow2(lmMst1)) + 160*Dmst12*(11389 - 704*
358  lmMst1 + 192*lmMt - 384*pow2(lmMst1))*pow2(Mst2) + 1920*(349 - 56*
359  lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow4(Mst2)))))))) + pow2(Mt)*(
360  3969000*oneLoopFlag*pow2(Mst1)*pow3(Mgl)*pow4(Msq)*(s2t*pow2(Dmst12)*(-
361  12*Mt*MuSUSY*pow2(Sbeta) + s2t*Tbeta*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
362  12*pow2(Mst1)*pow2(Sbeta))) + 12*Mt*pow2(Sbeta)*(2*Dmst12*MuSUSY*s2t*
363  pow2(Mst2) + Mt*Tbeta*(pow2(Dmst12) - 2*Dmst12*pow2(Mst2) + 4*(-lmMst1
364  + lmMt)*pow4(Mst2)))) + 16*Al4p*xDmglst1*pow2(Msq)*pow3(Dmglst1)*(
365  17640*Mst1*twoLoopFlag*pow2(Msq)*(-(pow2(Dmst12)*(8*((476 - 90*lmMst1 +
366  15*lmMt)*MuSUSY - 3*(51 + 10*lmMst1 - 10*lmMt)*Mst1*Tbeta)*pow2(Mt)*
367  pow2(Sbeta) + 5*Mst1*pow2(s2t)*(10*(-13 + 12*lmMst1)*Tbeta*pow2(MuSUSY)
368  *(-1 + pow2(Sbeta)) + 3*Mst1*((77 - 30*lmMst1)*MuSUSY + 4*(-71 + 60*
369  lmMst1)*Mst1*Tbeta)*pow2(Sbeta)) + 4*Mt*s2t*(130*Tbeta*pow2(MuSUSY)*(-1
370  + pow2(Sbeta)) + 3*Mst1*(120*MuSUSY + (277 - 155*lmMst1 + 5*lmMt)*Mst1*
371  Tbeta)*pow2(Sbeta)))) + 2*Dmst12*Mt*pow2(Mst2)*(Mt*((2231 - 990*lmMst1
372  + 90*lmMt)*MuSUSY + 12*(-301 + 290*lmMst1 + 10*lmMt)*Mst1*Tbeta)*pow2(
373  Sbeta) + s2t*(-5*(-77 + 30*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)
374  ) + (60*(71 - 60*lmMst1)*Mst1*MuSUSY + 8*(481 - 240*lmMst1 + 15*lmMt)*
375  Tbeta*pow2(Mst1))*pow2(Sbeta))) + 8*((977 - 480*lmMst1 + 30*lmMt)*
376  MuSUSY + 3*(-509 + 510*lmMst1 + 90*lmMt)*Mst1*Tbeta)*pow2(Mt)*pow2(
377  Sbeta)*pow4(Mst2)) + Al4p*threeLoopFlag*(4*Mst1*MuSUSY*pow2(Sbeta)*(
378  2940*Dmsqst1*(-20*Dmst12*Mt*((557 + 120*lmMst1 - 120*lmMt)*Mt + 3660*
379  Mst1*s2t)*pow2(Mst2) + 5*pow2(Dmst12)*(14640*Mst1*Mt*s2t + 4*(557 +
380  120*lmMst1 - 120*lmMt)*pow2(Mt) + 375*pow2(Mst1)*pow2(s2t)) + 24*(463 -
381  135*lmMst1 + 135*lmMt)*pow2(Mt)*pow4(Mst2)) - pow2(Msq)*(2*Dmst12*Mt*(
382  7350*Mst1*s2t*(548999 + 10980*lmMst1 + 288*lmMt - 19080*pow2(lmMst1)) +
383  Mt*(49723877243 + 60*lmMst1*(4936063 - 389970*lmMt) + 57352680*lmMt -
384  342543600*pow2(lmMst1) - 3175200*pow2(lmMt)))*pow2(Mst2) - 15*pow2(
385  Dmst12)*(490*Mst1*Mt*s2t*(-611423 + 9984*lmMst1 - 768*lmMt + 23040*
386  pow2(lmMst1)) + (5753390765 + 580*lmMst1*(79969 - 1932*lmMt) + 7091364*
387  lmMt - 35700000*pow2(lmMst1) - 282240*pow2(lmMt))*pow2(Mt) + 49*(
388  31025111 + 290880*lmMst1 - 251100*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)) +
389  392*(122282257 + 60*lmMst1*(8318 - 3885*lmMt) + 479550*lmMt - 1351800*
390  pow2(lmMst1) - 21600*pow2(lmMt))*pow2(Mt)*pow4(Mst2))) + Tbeta*(245*
391  pow2(MuSUSY)*(2*Dmst12*Mt*(30000*Dmsqst1*Mst1*s2t + (4*Mst1*s2t*(
392  31025111 + 290880*lmMst1 - 251100*pow2(lmMst1)) + Mt*(59957863 +
393  480000*lmMst1 - 26880*lmMt - 230400*pow2(lmMst1)))*pow2(Msq))*pow2(
394  Mst2) + pow2(Dmst12)*(-60000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(2*Mst1*
395  Mt*s2t*(-99874229 - 1352280*lmMst1 + 633600*pow2(lmMst1)) + 2*(-
396  59957863 - 480000*lmMst1 + 26880*lmMt + 230400*pow2(lmMst1))*pow2(Mt) +
397  15*(3044017 + 27472*lmMst1 - 48480*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))
398  - 24*(-3877891 - 46400*lmMst1 + 960*lmMt + 19200*pow2(lmMst1))*pow2(
399  Msq)*pow2(Mt)*pow4(Mst2)) + pow2(Sbeta)*(2*Dmst12*Mst1*Mt*s2t*(11760*
400  Dmsqst1*(4*Dmst12*(-2729 + 105*lmMst1 - 105*lmMt)*pow2(Mst1) + 6*(791 -
401  270*lmMst1 + 270*lmMt)*pow2(Mst1)*pow2(Mst2) + 625*Dmst12*pow2(MuSUSY)
402  - 625*pow2(Mst2)*pow2(MuSUSY)) - pow2(Msq)*(196*pow2(Mst2)*(8*(61021241
403  + 15*lmMst1*(20521 - 7770*lmMt) + 250575*lmMt - 675900*pow2(lmMst1) -
404  10800*pow2(lmMt))*pow2(Mst1) + 5*(31025111 + 290880*lmMst1 - 251100*
405  pow2(lmMst1))*pow2(MuSUSY)) + Dmst12*(4*(25774874431 - 37697520*lmMt +
406  600*lmMst1*(311999 + 37149*lmMt) - 77590800*pow2(lmMst1) + 1058400*
407  pow2(lmMt))*pow2(Mst1) + 245*(-99874229 - 1352280*lmMst1 + 633600*pow2(
408  lmMst1))*pow2(MuSUSY)))) - 3675*pow2(Dmst12)*pow2(Mst1)*(8*(14640*
409  Dmsqst1 + (548855 + 10980*lmMst1 + 288*lmMt - 19080*pow2(lmMst1))*pow2(
410  Msq))*pow2(Mst1) + (3044017 + 27472*lmMst1 - 48480*pow2(lmMst1))*pow2(
411  Msq)*pow2(MuSUSY))*pow2(s2t) - 2*pow2(Mt)*(47040*Dmsqst1*pow2(Mst1)*(3*
412  Dmst12*(3401 + 105*lmMst1 - 105*lmMt)*(Dmst12 - pow2(Mst2)) + (-19241 +
413  420*lmMst1 - 420*lmMt)*pow4(Mst2)) + pow2(Msq)*(Dmst12*pow2(Mst2)*(16*(
414  2524164367 + 8198205*lmMst1 + 19269705*lmMt - 2463300*lmMst1*lmMt +
415  27997200*pow2(lmMst1) - 1058400*pow2(lmMt))*pow2(Mst1) - 245*(-59957863
416  - 480000*lmMst1 + 26880*lmMt + 230400*pow2(lmMst1))*pow2(MuSUSY)) +
417  pow2(Dmst12)*(4*(-7672052891 - 10289580*lmMt + 900*lmMst1*(15649 +
418  11284*lmMt) + 98506800*pow2(lmMst1) + 4233600*pow2(lmMt))*pow2(Mst1) +
419  245*(-59957863 - 480000*lmMst1 + 26880*lmMt + 230400*pow2(lmMst1))*
420  pow2(MuSUSY)) - 588*(4*(-21126629 - 218510*lmMt + 20*lmMst1*(-28958 +
421  5055*lmMt) + 194100*pow2(lmMst1) + 91800*pow2(lmMt))*pow2(Mst1) + 5*(-
422  3877891 - 46400*lmMst1 + 960*lmMt + 19200*pow2(lmMst1))*pow2(MuSUSY))*
423  pow4(Mst2)))))))) + 23520*Mst1*threeLoopFlag*xDmsqst1*pow2(Al4p)*pow2(
424  Dmsqst1)*(4500*Mst1*Tbeta*pow2(Mt)*pow3(Mgl)*(24*shiftst2*(pow2(Dmst12)
425  + 2*(-2 + lmMst1)*pow2(Mst2)*(Dmst12 + pow2(Mst2)))*pow2(Mt)*pow2(
426  Sbeta) + Dmst12*shiftst2*pow2(s2t)*(-(Dmst12*(-1 + 2*lmMst1)*pow2(
427  MuSUSY)*(-1 + pow2(Sbeta))) + 12*Dmst12*(-2 + lmMst1)*pow2(Mst1)*pow2(
428  Sbeta) - 4*(-2 + lmMst1)*pow2(Mst2)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
429  12*pow2(Mst1)*pow2(Sbeta))) - Dmst12*shiftst1*pow2(s2t)*(Dmst12*(-5 +
430  2*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 36*Dmst12*(-2 + lmMst1)*
431  pow2(Mst1)*pow2(Sbeta) - 4*(-2 + lmMst1)*pow2(Mst2)*(pow2(MuSUSY)*(-1 +
432  pow2(Sbeta)) + 12*pow2(Mst1)*pow2(Sbeta))) + 48*(-2 + lmMst1)*shiftst1*
433  pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) + 8*xDmglst1*pow2(Mt)*pow3(Dmglst1)*(-
434  2*Dmst12*Mt*pow2(Mst2)*(2*Mt*(5*(557 + 120*lmMst1 - 120*lmMt)*MuSUSY -
435  6*(3401 + 105*lmMst1 - 105*lmMt)*Mst1*Tbeta)*pow2(Sbeta) + s2t*(625*
436  Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(36600*MuSUSY + (-7121 +
437  870*lmMst1 - 870*lmMt)*Mst1*Tbeta)*pow2(Sbeta))) + pow2(Dmst12)*((4*(5*
438  (557 + 120*lmMst1 - 120*lmMt)*MuSUSY - 6*(3401 + 105*lmMst1 - 105*lmMt)
439  *Mst1*Tbeta)*pow2(Mt) + 75*(25*MuSUSY - 488*Mst1*Tbeta)*pow2(Mst1)*
440  pow2(s2t))*pow2(Sbeta) + 2*Mt*s2t*(625*Tbeta*pow2(MuSUSY)*(-1 + pow2(
441  Sbeta)) + Mst1*(36600*MuSUSY + (-13291 - 330*lmMst1 + 330*lmMt)*Mst1*
442  Tbeta)*pow2(Sbeta))) + 4*((3778 - 435*lmMst1 + 435*lmMt)*MuSUSY + 25*(
443  1225 - 39*lmMst1 + 39*lmMt)*Mst1*Tbeta)*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)
444  ) + Mgl*(5*Tbeta*pow2(Mt)*(5*Dmst12*s2t*(-30*Dmst12*(7 + 24*lmMst1)*
445  Mst1*s2t*pow2(Mgl)*pow2(MuSUSY) - 80*Dmst12*Mt*(5*Dmglst1*Mgl + 5*pow2(
446  Dmglst1) + 21*pow2(Mgl))*pow2(MuSUSY) + 80*Mt*(5*Dmglst1*Mgl + 5*pow2(
447  Dmglst1) + 21*pow2(Mgl))*pow2(Mst2)*pow2(MuSUSY) - 3*Dmst12*Mst1*s2t*((
448  3760*Dmglst1*Mgl + 3832*pow2(Dmglst1) - 5*(229 + 288*lmMst1)*pow2(Mgl))
449  *pow2(Mst1) - 10*(7 + 24*lmMst1)*pow2(Mgl)*pow2(MuSUSY))*pow2(Sbeta)) +
450  4*Mt*pow2(Sbeta)*(4*Dmst12*s2t*pow2(Mst1)*(50*pow2(Mgl)*(Dmst12*(23 -
451  6*lmMst1 + 6*lmMt) + (-85 + 12*lmMst1 - 12*lmMt)*pow2(Mst2)) + Dmglst1*
452  Mgl*(Dmst12*(1041 - 90*lmMst1 + 90*lmMt) - 25*(91 + 6*lmMst1 - 6*lmMt)*
453  pow2(Mst2)) + pow2(Dmglst1)*(9*Dmst12*(-149 + 10*lmMst1 - 10*lmMt) + (
454  107 - 330*lmMst1 + 330*lmMt)*pow2(Mst2))) + 100*Dmst12*s2t*(5*Dmglst1*
455  Mgl + 5*pow2(Dmglst1) + 21*pow2(Mgl))*(Dmst12 - pow2(Mst2))*pow2(
456  MuSUSY) - (4*Mst1*Mt*pow2(Dmglst1)*(3*Dmst12*(13147 - 90*lmMst1 + 90*
457  lmMt)*(Dmst12 - pow2(Mst2)) + (-64033 + 3360*lmMst1 - 3360*lmMt)*pow4(
458  Mst2)))/5. - 8*Dmglst1*Mgl*Mst1*Mt*(9*Dmst12*(423 - 20*lmMst1 + 20*
459  lmMt)*(Dmst12 - pow2(Mst2)) + 25*(-226 + 21*lmMst1 - 21*lmMt)*pow4(
460  Mst2)) + Mst1*Mt*pow2(Mgl)*(6*(419 - 60*lmMst1 + 60*lmMt)*pow2(Dmst12)
461  - 225*Dmst12*(49 + 46*lmMst1 + 2*lmMt)*pow2(Mst2) + 50*(205 + 174*lmMt
462  - 6*lmMst1*(101 + 18*lmMt) + 54*(pow2(lmMst1) + pow2(lmMt)))*pow4(Mst2)
463  ))) + 50*MuSUSY*pow2(Sbeta)*(60*pow2(Dmst12)*(5*Dmglst1*Mgl + 5*pow2(
464  Dmglst1) + 21*pow2(Mgl))*pow2(Mst1)*pow2(Mt)*pow2(s2t) + Mst1*(3*
465  Dmst12*s2t*(8*Dmglst1*(479*Dmglst1 + 470*Mgl)*(Dmst12 - pow2(Mst2)) -
466  5*pow2(Mgl)*(44*Dmst12 - (229 + 288*lmMst1)*pow2(Mst2)))*pow3(Mt) +
467  540*Mt*s2t*pow2(Mgl)*(-8*Dmst12*(-2 + lmMst1)*shiftst2*pow2(Mst2)*pow2(
468  Mt) + pow2(Dmst12)*(-4*shiftst2*pow2(Mt) - (-2 + lmMst1)*(shiftst1 -
469  shiftst2)*pow2(Mst1)*pow2(s2t)) + 8*(-2 + lmMst1)*(shiftst1 - shiftst2)
470  *pow2(Mt)*pow4(Mst2))) + (16*(25*pow2(Mgl)*(Dmst12*(65 - 6*lmMst1 + 6*
471  lmMt)*(Dmst12 - pow2(Mst2)) + (-91 + 12*lmMst1 - 12*lmMt)*pow4(Mst2)) +
472  Dmglst1*Mgl*(Dmst12*(557 + 120*lmMst1 - 120*lmMt)*(Dmst12 - pow2(Mst2))
473  - 25*(44 + 3*lmMst1 - 3*lmMt)*pow4(Mst2)) + pow2(Dmglst1)*(Dmst12*(557
474  + 120*lmMst1 - 120*lmMt)*(Dmst12 - pow2(Mst2)) + (136 - 165*lmMst1 +
475  165*lmMt)*pow4(Mst2)))*pow4(Mt))/5.)))) + 11025*threeLoopFlag*pow2(
476  Al4p)*(384*z2*pow2(Mst1)*pow3(Mgl)*(pow2(Mst2)*(-50*Mt*xDmsqst1*pow2(
477  Dmsqst1)*(-(pow2(Dmst12)*pow2(s2t)*(-(Mt*(shiftst1 + shiftst2)*Tbeta*
478  pow2(MuSUSY)*(-1 + pow2(Sbeta))) - 3*(MuSUSY*s2t*(shiftst1 - shiftst2)
479  + 2*Mt*(3*shiftst1 - shiftst2)*Tbeta)*pow2(Mst1)*pow2(Sbeta))) + 2*
480  Dmst12*Mt*pow2(Mst2)*(12*Mt*MuSUSY*s2t*shiftst2*pow2(Sbeta) - 12*
481  shiftst2*Tbeta*pow2(Mt)*pow2(Sbeta) - (shiftst1 - shiftst2)*Tbeta*pow2(
482  s2t)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*pow2(Mst1)*pow2(Sbeta))) -
483  24*(MuSUSY*s2t*(shiftst1 - shiftst2) + Mt*(shiftst1 + shiftst2)*Tbeta)*
484  pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) + pow2(Msq)*(50*Tbeta*(Dmsqst1 + pow2(
485  Msq))*pow2(Mt)*(-(pow2(Dmst12)*pow2(s2t)*((shiftst1 + shiftst2)*pow2(
486  MuSUSY)*(-1 + pow2(Sbeta)) + 6*(3*shiftst1 - shiftst2)*pow2(Mst1)*pow2(
487  Sbeta))) + 2*Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(MuSUSY)*
488  pow2(s2t)*(-1 + pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (shiftst1 -
489  shiftst2)*pow2(Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1 + shiftst2)
490  *pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) + 15*MuSUSY*pow2(Sbeta)*(10*(Dmsqst1
491  + pow2(Msq))*(-8*Dmst12*s2t*shiftst2*pow2(Mst2)*pow3(Mt) + Mt*(-
492  shiftst1 + shiftst2)*pow2(Dmst12)*pow2(Mst1)*pow3(s2t) + 8*s2t*(
493  shiftst1 - shiftst2)*pow3(Mt)*pow4(Mst2)) + shiftst3*pow2(Msq)*(-(Mt*
494  pow2(Dmst12)*pow2(Mst1)*pow3(s2t)) + 8*s2t*pow3(Mt)*(pow2(Dmst12) -
495  Dmst12*pow2(Mst2) + pow4(Mst2))))) + 5*shiftst3*Tbeta*pow4(Msq)*(-(
496  Dmst12*pow2(Mt)*pow2(s2t)*((3*Dmst12 - 2*pow2(Mst2))*pow2(MuSUSY)*(-1 +
497  pow2(Sbeta)) + 6*pow2(Mst1)*(7*Dmst12 - 4*pow2(Mst2))*pow2(Sbeta))) +
498  24*pow2(Sbeta)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2))*pow4(Mt)
499  )) - xDmst12*pow3(Dmst12)*(50*Dmsqst1*(Dmsqst1*xDmsqst1 + pow2(Msq))*(-
500  4*shiftst2*pow2(Mt)*(Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) +
501  6*Mt*(MuSUSY*s2t - Mt*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(2*Mt*(MuSUSY*
502  s2t*(-2*shiftst1 + shiftst2) + 2*Mt*(-4*shiftst1 + shiftst2)*Tbeta)*
503  pow2(Mst1)*pow2(s2t) + (shiftst1 - shiftst2)*Tbeta*pow4(Mst1)*pow4(s2t)
504  )) + pow4(Msq)*(2*pow2(Mt)*(-((15*shiftst1 + 10*shiftst2 + 9*shiftst3)*
505  Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta))) + 60*Mt*shiftst3*(
506  MuSUSY*s2t + Mt*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-5*Mt*(MuSUSY*s2t*(
507  40*shiftst1 - 20*shiftst2 + 7*shiftst3) + 2*Mt*(80*shiftst1 - 20*
508  shiftst2 + 29*shiftst3)*Tbeta)*pow2(Mst1)*pow2(s2t) + 5*(10*shiftst1 -
509  10*shiftst2 + shiftst3)*Tbeta*pow4(Mst1)*pow4(s2t))))) - 5*Mt*z3*(Mt*
510  pow2(Mst2)*(8*Dmst12*Mgl*Mst1*s2t*Tbeta*pow2(Msq)*(Dmst12*Mst1*s2t*(12*
511  (420*Dmsqst1*pow2(Mgl) + (2898*Dmglst1*Mgl - 3295*pow2(Dmglst1) + 660*
512  pow2(Mgl))*pow2(Msq))*pow2(Mst1) + (1260*Dmsqst1*pow2(Mgl) - (50485*
513  Dmglst1*Mgl + 200115*pow2(Dmglst1) + 2574*pow2(Mgl))*pow2(Msq))*pow2(
514  MuSUSY)) - 2*Mt*pow2(Msq)*(12*pow2(Mst1)*(Dmst12*(22100*Dmglst1*Mgl +
515  86801*pow2(Dmglst1) + 3290*pow2(Mgl)) + 8*(1963*Dmglst1*Mgl + 9451*
516  pow2(Dmglst1) + 106*pow2(Mgl))*pow2(Mst2)) - (5*Dmst12*(-21081*Dmglst1*
517  Mgl + 21081*pow2(Dmglst1) - 3457*pow2(Mgl)) - 2*(37582*Dmglst1*Mgl +
518  142987*pow2(Dmglst1) + 3012*pow2(Mgl))*pow2(Mst2))*pow2(MuSUSY)))*pow2(
519  Sbeta) - 5040*xDmsqst1*pow2(Dmsqst1)*pow2(Mst1)*pow3(Mgl)*(6*Dmst12*Mt*
520  (-5*MuSUSY*s2t + 6*Mt*Tbeta)*pow2(Mst2)*pow2(Sbeta) + Tbeta*pow2(
521  Dmst12)*pow2(s2t)*(-2*pow2(MuSUSY)*(-1 + pow2(Sbeta)) - 15*pow2(Mst1)*
522  pow2(Sbeta)) + 24*Tbeta*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) - 8*xDmglst1*
523  pow3(Dmglst1)*pow4(Msq)*(6*Dmst12*Mt*pow2(Mst2)*(147351*Mt*Tbeta*pow2(
524  MuSUSY)*(-1 + pow2(Sbeta)) + 48*(3977*MuSUSY*s2t + 7908*Mt*Tbeta)*pow2(
525  Mst1)*pow2(Sbeta) + 4*Mst1*MuSUSY*(76518*MuSUSY*s2t*Tbeta*(-1 + pow2(
526  Sbeta)) + 500581*Mt*pow2(Sbeta)) + 966992*s2t*Tbeta*pow2(Sbeta)*pow3(
527  Mst1)) + pow2(Dmst12)*(-2*Mst1*Mt*MuSUSY*(737647*MuSUSY*s2t*Tbeta*(-1 +
528  pow2(Sbeta)) + 5205992*Mt*pow2(Sbeta)) + pow2(Mst1)*(349745*Tbeta*pow2(
529  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) - 24*Mt*(-22939*MuSUSY*s2t +
530  77043*Mt*Tbeta)*pow2(Sbeta)) - 24*s2t*(114777*MuSUSY*s2t - 258833*Mt*
531  Tbeta)*pow2(Sbeta)*pow3(Mst1) + Tbeta*(-884106*pow2(Mt)*pow2(MuSUSY)*(-
532  1 + pow2(Sbeta)) + 572688*pow2(s2t)*pow2(Sbeta)*pow4(Mst1))) + 24*pow2(
533  Mt)*(28665*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 4*Mst1*(60437*MuSUSY
534  + 28466*Mst1*Tbeta)*pow2(Sbeta))*pow4(Mst2)) + Mgl*pow2(Msq)*(-3*Tbeta*
535  pow2(Mt)*pow2(Sbeta)*(20160*Dmsqst1*Dmst12*pow2(Mgl)*pow2(Mst1)*(Dmst12
536  + 2*pow2(Mst2)) + pow2(Msq)*(pow2(Mgl)*(160*Dmst12*pow2(Mst2)*(2052*
537  pow2(Mst1) + 425*pow2(MuSUSY)) + pow2(Dmst12)*(29472*pow2(Mst1) +
538  31963*pow2(MuSUSY)) + 512*((-214 + 24*lmMst1 - 24*lmMt)*pow2(Mst1) +
539  47*pow2(MuSUSY))*pow4(Mst2)) - 48*Dmglst1*Mgl*(pow2(Dmst12)*(6970*pow2(
540  Mst1) + 8213*pow2(MuSUSY)) - Dmst12*pow2(Mst2)*(30912*pow2(Mst1) +
541  8213*pow2(MuSUSY)) - 8*(3356*pow2(Mst1) + 583*pow2(MuSUSY))*pow4(Mst2))
542  - 48*pow2(Dmglst1)*(7*pow2(Dmst12)*(7481*pow2(Mst1) + 4095*pow2(MuSUSY)
543  ) - 5*Dmst12*pow2(Mst2)*(15248*pow2(Mst1) + 5733*pow2(MuSUSY)) - 2*(
544  37624*pow2(Mst1) + 7047*pow2(MuSUSY))*pow4(Mst2)))) - Tbeta*pow2(
545  MuSUSY)*(8*pow2(Dmst12)*(1260*Dmsqst1*pow2(Mgl) - (50485*Dmglst1*Mgl +
546  200115*pow2(Dmglst1) + 2574*pow2(Mgl))*pow2(Msq))*pow2(Mst1)*pow2(s2t)
547  + Mt*pow2(Msq)*(16*Dmst12*Mst1*s2t*(5*Dmst12*(-21081*Dmglst1*Mgl +
548  21081*pow2(Dmglst1) - 3457*pow2(Mgl)) - 2*(37582*Dmglst1*Mgl + 142987*
549  pow2(Dmglst1) + 3012*pow2(Mgl))*pow2(Mst2)) + 3*Mt*(48*Dmglst1*Mgl*(
550  8213*pow2(Dmst12) - 8213*Dmst12*pow2(Mst2) - 4664*pow4(Mst2)) + 432*
551  pow2(Dmglst1)*(3185*pow2(Dmst12) - 3185*Dmst12*pow2(Mst2) - 1566*pow4(
552  Mst2)) - pow2(Mgl)*(31963*pow2(Dmst12) + 68000*Dmst12*pow2(Mst2) +
553  24064*pow4(Mst2))))) + 16*Mst1*MuSUSY*pow2(Sbeta)*(6*Dmst12*Mst1*Mt*
554  s2t*(105*Dmsqst1*pow2(Mgl)*(7*Dmst12 + 8*pow2(Mst2)) + pow2(Msq)*(3*
555  pow2(Mgl)*(403*Dmst12 + 440*pow2(Mst2)) - 2*pow2(Dmglst1)*(32498*Dmst12
556  + 3295*pow2(Mst2)) + Dmglst1*(-7642*Dmst12*Mgl + 5796*Mgl*pow2(Mst2))))
557  + pow2(Msq)*(3*pow2(Dmst12)*(37582*Dmglst1*Mgl + 142987*pow2(Dmglst1) +
558  3012*pow2(Mgl))*pow2(Mst1)*pow2(s2t) + pow2(Mt)*(4*pow2(Dmglst1)*(
559  286982*pow2(Dmst12) - 487227*Dmst12*pow2(Mst2) - 226824*pow4(Mst2)) -
560  pow2(Mgl)*(51181*pow2(Dmst12) + 49656*Dmst12*pow2(Mst2) + 10176*pow4(
561  Mst2)) - 4*Dmglst1*Mgl*(86833*pow2(Dmst12) + 113412*Dmst12*pow2(Mst2) +
562  47112*pow4(Mst2))))))) + 2*xDmst12*pow3(Dmst12)*(pow3(Mgl)*(-2520*
563  Dmsqst1*pow2(Mst1)*(Dmsqst1*xDmsqst1*(Mt*Tbeta*pow2(s2t)*(4*pow2(
564  MuSUSY)*(-1 + pow2(Sbeta)) + 15*pow2(Mst1)*pow2(Sbeta)) - 36*Tbeta*
565  pow2(Sbeta)*pow3(Mt) + 2*MuSUSY*pow2(Sbeta)*(15*s2t*pow2(Mt) + pow2(
566  Mst1)*pow3(s2t))) + pow2(Msq)*(Mt*Tbeta*pow2(s2t)*(4*pow2(MuSUSY)*(-1 +
567  pow2(Sbeta)) + pow2(Mst1)*pow2(Sbeta)) - 48*Tbeta*pow2(Sbeta)*pow3(Mt)
568  + 2*MuSUSY*pow2(Sbeta)*(22*s2t*pow2(Mt) + pow2(Mst1)*pow3(s2t)))) +
569  pow4(Msq)*(16*Mst1*MuSUSY*pow2(Mt)*(20297*MuSUSY*s2t*Tbeta*(-1 + pow2(
570  Sbeta)) + 76009*Mt*pow2(Sbeta)) - Mt*pow2(Mst1)*(37669*Tbeta*pow2(
571  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) - 4*Mt*(-32783*MuSUSY*s2t + 33933*
572  Mt*Tbeta)*pow2(Sbeta)) + 4*s2t*(Mt*(33783*MuSUSY*s2t - 23402*Mt*Tbeta)
573  + 18*Mst1*s2t*(143*MuSUSY*s2t - 37*Mt*Tbeta))*pow2(Sbeta)*pow3(Mst1) +
574  Tbeta*(197889*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 24096*pow2(
575  Sbeta)*pow3(s2t)*pow5(Mst1)))) + 4*Dmglst1*pow4(Msq)*(pow2(Mgl)*(4*
576  Mst1*MuSUSY*pow2(Mt)*(142987*MuSUSY*s2t*Tbeta*(-1 + pow2(Sbeta)) +
577  574156*Mt*pow2(Sbeta)) + 2*Mt*pow2(Mst1)*(50485*Tbeta*pow2(MuSUSY)*
578  pow2(s2t)*(-1 + pow2(Sbeta)) + 12*Mt*(4744*MuSUSY*s2t + 12729*Mt*Tbeta)
579  *pow2(Sbeta)) + s2t*(Mt*(90723*MuSUSY*s2t - 164264*Mt*Tbeta) + Mst1*
580  s2t*(50485*MuSUSY*s2t - 80628*Mt*Tbeta))*pow2(Sbeta)*pow3(Mst1) + 86*
581  Tbeta*(-1719*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 874*pow2(Sbeta)
582  *pow3(s2t)*pow5(Mst1))) + Dmglst1*Mgl*(-8*Mst1*MuSUSY*pow2(Mt)*(-18791*
583  MuSUSY*s2t*Tbeta*(-1 + pow2(Sbeta)) + 86737*Mt*pow2(Sbeta)) + 6*Mt*
584  pow2(Mst1)*(66705*Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + (
585  273164*Mt*MuSUSY*s2t - 85482*Tbeta*pow2(Mt))*pow2(Sbeta)) + s2t*(-3*Mt*
586  s2t*(391379*MuSUSY + 116812*Mst1*Tbeta) + 4379080*Tbeta*pow2(Mt) +
587  200115*Mst1*MuSUSY*pow2(s2t))*pow2(Sbeta)*pow3(Mst1) + Tbeta*(-515970*
588  pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 285974*pow2(Sbeta)*pow3(s2t)
589  *pow5(Mst1))) + xDmglst1*pow2(Dmglst1)*(-4*Mst1*MuSUSY*pow2(Mt)*(
590  278539*MuSUSY*s2t*Tbeta*(-1 + pow2(Sbeta)) + 2202506*Mt*pow2(Sbeta)) +
591  2*Mt*pow2(Mst1)*(349745*Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta))
592  - 24*Mt*(-46801*MuSUSY*s2t + 29595*Mt*Tbeta)*pow2(Sbeta)) + s2t*(-
593  4967589*Mt*MuSUSY*s2t + 297420*Mst1*Mt*s2t*Tbeta + 16623976*Tbeta*pow2(
594  Mt) + 349745*Mst1*MuSUSY*pow2(s2t))*pow2(Sbeta)*pow3(Mst1) + Tbeta*(-
595  884106*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 918216*pow2(Sbeta)*
596  pow3(s2t)*pow5(Mst1))))))) + 4*xDmst12*pow3(Dmst12)*(5880*Mst1*
597  threeLoopFlag*xDmsqst1*pow2(Al4p)*pow2(Dmsqst1)*(8*Dmglst1*Mt*(5*pow2(
598  Mgl)*(-2*s2t*pow2(Mt)*(125*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
599  Mst1*(7050*MuSUSY + (-193 - 330*lmMst1 + 330*lmMt)*Mst1*Tbeta)*pow2(
600  Sbeta)) + pow2(Sbeta)*(150*Mt*(-5*MuSUSY + 94*Mst1*Tbeta)*pow2(Mst1)*
601  pow2(s2t) - 4*((557 + 120*lmMst1 - 120*lmMt)*MuSUSY + 9*(-423 + 20*
602  lmMst1 - 20*lmMt)*Mst1*Tbeta)*pow3(Mt) + 125*Tbeta*pow3(s2t)*pow4(Mst1)
603  )) + xDmglst1*pow2(Dmglst1)*(-2*s2t*pow2(Mt)*(625*Tbeta*pow2(MuSUSY)*(-
604  1 + pow2(Sbeta)) + (36600*Mst1*MuSUSY - 3*(6487 + 510*lmMst1 - 510*
605  lmMt)*Tbeta*pow2(Mst1))*pow2(Sbeta)) + pow2(Sbeta)*(150*Mt*(-25*MuSUSY
606  + 488*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t) - 4*(5*(557 + 120*lmMst1 - 120*
607  lmMt)*MuSUSY - 6*(3401 + 105*lmMst1 - 105*lmMt)*Mst1*Tbeta)*pow3(Mt) +
608  625*Tbeta*pow3(s2t)*pow4(Mst1))) + Dmglst1*Mgl*(-50*s2t*pow2(Mt)*(25*
609  Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(1437*MuSUSY - 5*(103 + 6*
610  lmMst1 - 6*lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(150*Mt*(-25*
611  MuSUSY + 479*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t) + (-20*(557 + 120*lmMst1
612  - 120*lmMt)*MuSUSY + 6*(13147 - 90*lmMst1 + 90*lmMt)*Mst1*Tbeta)*pow3(
613  Mt) + 625*Tbeta*pow3(s2t)*pow4(Mst1)))) + 5*pow3(Mgl)*(-75*Mst1*pow2(
614  Mt)*pow2(s2t)*(4*(7 + 6*shiftst1 - 24*lmMst1*(-1 + shiftst2) + 30*
615  shiftst2)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*(112*MuSUSY +
616  Mst1*(91 + 240*shiftst1 + 32*lmMst1*(3 - 4*shiftst1 + shiftst2))*Tbeta)
617  *pow2(Sbeta)) - 150*s2t*(56*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
618  Mst1*(3*MuSUSY*(47 + 96*(lmMst1 + shiftst2 - lmMst1*shiftst2)) - 208*
619  Mst1*Tbeta)*pow2(Sbeta))*pow3(Mt) + pow2(Sbeta)*(-150*Mt*(MuSUSY*(7 +
620  108*shiftst1 - 72*shiftst2 + 24*lmMst1*(1 - 2*shiftst1 + shiftst2)) -
621  28*Mst1*Tbeta)*pow3(Mst1)*pow3(s2t) - 4*(200*(65 - 6*lmMst1 + 6*lmMt)*
622  MuSUSY - 3*Mst1*(1999 - 90*lmMt + 90*lmMst1*(41 - 40*shiftst2) + 3600*
623  shiftst2)*Tbeta)*pow4(Mt) - 1800*(-2 + lmMst1)*(shiftst1 - shiftst2)*
624  Tbeta*pow4(s2t)*pow5(Mst1)))) + Al4p*Mgl*pow2(Msq)*(-35280*Mst1*Mt*
625  twoLoopFlag*pow2(Msq)*(pow2(Mgl)*(25*Mst1*Mt*pow2(s2t)*((1 + 6*lmMst1)*
626  Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*Mst1*(MuSUSY + 3*lmMst1*
627  MuSUSY + (1 + 12*lmMst1)*Mst1*Tbeta)*pow2(Sbeta)) + 100*s2t*pow2(Mt)*(-
628  ((5 + 6*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta))) + 2*Mst1*(-12*(-
629  1 + lmMst1)*MuSUSY + (1 + 3*lmMst1 + 9*lmMt)*Mst1*Tbeta)*pow2(Sbeta)) +
630  pow2(Sbeta)*(-4*(50*(5 + 6*lmMst1)*MuSUSY + (137 - 330*lmMst1 - 270*
631  lmMt)*Mst1*Tbeta)*pow3(Mt) + 75*(MuSUSY - 2*(1 + 2*lmMst1)*Mst1*Tbeta)*
632  pow3(Mst1)*pow3(s2t))) + 2*Dmglst1*(Mgl*(-25*Mst1*Mt*pow2(s2t)*(8*(-1 +
633  3*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*((5 - 6*
634  lmMst1)*MuSUSY + 6*(-1 + 4*lmMst1)*Mst1*Tbeta)*pow2(Sbeta)) - 2*s2t*
635  pow2(Mt)*(30*lmMst1*(5*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(
636  60*MuSUSY - 11*Mst1*Tbeta)*pow2(Sbeta)) + Tbeta*(-125*pow2(MuSUSY)*(-1
637  + pow2(Sbeta)) + 2*(131 - 135*lmMt)*pow2(Mst1)*pow2(Sbeta))) + pow2(
638  Sbeta)*(4*(6*(17 - 30*lmMst1 + 5*lmMt)*MuSUSY + (47 + 870*lmMst1 + 30*
639  lmMt)*Mst1*Tbeta)*pow3(Mt) + 25*((4 - 12*lmMst1)*MuSUSY + (5 - 6*
640  lmMst1)*Mst1*Tbeta)*pow3(Mst1)*pow3(s2t))) + Dmglst1*(-15*Mst1*Mt*pow2(
641  s2t)*(10*(-5 + 6*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(
642  (85 - 30*lmMst1)*MuSUSY + (-137 + 180*lmMst1)*Mst1*Tbeta)*pow2(Sbeta))
643  + 10*s2t*pow2(Mt)*(-5*(1 + 6*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(
644  Sbeta)) + 6*Mst1*((58 - 90*lmMst1)*MuSUSY + (2 + 15*lmMst1 + 5*lmMt)*
645  Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-2*(3*(259 + 90*lmMst1 + 10*
646  lmMt)*MuSUSY + (1547 - 2730*lmMst1 + 30*lmMt)*Mst1*Tbeta)*pow3(Mt) +
647  25*(3*(5 - 6*lmMst1)*MuSUSY + (11 - 6*lmMst1)*Mst1*Tbeta)*pow3(Mst1)*
648  pow3(s2t))))) - Al4p*threeLoopFlag*(196*Mst1*Mt*MuSUSY*pow2(Sbeta)*(
649  Dmglst1*Mgl*(2400*Dmsqst1*Mt*(7050*Mst1*Mt*s2t + 2*(557 + 120*lmMst1 -
650  120*lmMt)*pow2(Mt) + 375*pow2(Mst1)*pow2(s2t)) - pow2(Msq)*(3600*Mst1*
651  s2t*(12383 + 4128*lmMst1 + 80*lmMt - 1260*pow2(lmMst1))*pow2(Mt) + 225*
652  Mt*(284641 + 8696*lmMst1 + 1680*pow2(lmMst1))*pow2(Mst1)*pow2(s2t) + 8*
653  (193364399 + 90000*lmMt - 300*lmMst1*(3781 + 582*lmMt) - 2005200*pow2(
654  lmMst1) - 43200*pow2(lmMt))*pow3(Mt) + 375*(84209 + 1264*lmMst1 - 240*
655  pow2(lmMst1))*pow3(Mst1)*pow3(s2t))) + 15*pow2(Mgl)*(500*Dmsqst1*(6*
656  Mst1*s2t*(173 - 144*lmMst1*(-1 + shiftst2) + 216*shiftst2)*pow2(Mt) +
657  504*Mt*pow2(Mst1)*pow2(s2t) + 16*(65 - 6*lmMst1 + 6*lmMt)*pow3(Mt) + 3*
658  (7 + 72*shiftst1 - 36*shiftst2 + 24*lmMst1*(1 - 2*shiftst1 + shiftst2))
659  *pow3(Mst1)*pow3(s2t)) - pow2(Msq)*(15*Mst1*s2t*(-102747 + 640*lmMt +
660  6720*shiftst3 - 32*lmMst1*(331 + 90*shiftst3) + 13888*pow2(lmMst1))*
661  pow2(Mt) - 75*Mt*(-20531 + 200*lmMst1 + 1200*pow2(lmMst1))*pow2(Mst1)*
662  pow2(s2t) - 4*(-3454599 + 16840*lmMt + 48*lmMst1*(262 + 405*lmMt) +
663  46560*pow2(lmMst1))*pow3(Mt) + 50*(1429 - 720*shiftst1 + 360*shiftst2 -
664  234*shiftst3 + 2*lmMst1*(-227 + 720*shiftst1 - 360*shiftst2 + 126*
665  shiftst3) + 24*pow2(lmMst1))*pow3(Mst1)*pow3(s2t))) + pow2(Dmglst1)*(
666  2400*Dmsqst1*Mt*(7185*Mst1*Mt*s2t + 2*(557 + 120*lmMst1 - 120*lmMt)*
667  pow2(Mt) + 375*pow2(Mst1)*pow2(s2t)) + pow2(Msq)*(-7200*Mst1*s2t*(
668  143196 + 2546*lmMst1 + 92*lmMt - 4785*pow2(lmMst1))*pow2(Mt) + 225*Mt*(
669  3532083 + 36328*lmMst1 - 47760*pow2(lmMst1))*pow2(Mst1)*pow2(s2t) + 16*
670  (29818901 + 258300*lmMt + 30*lmMst1*(35963 + 2190*lmMt) - 239400*pow2(
671  lmMst1) - 10800*pow2(lmMt))*pow3(Mt) + 75*(-1732531 - 16896*lmMst1 +
672  24840*pow2(lmMst1))*pow3(Mst1)*pow3(s2t)))) + Tbeta*(-245*(24*pow2(
673  Dmglst1)*(10000*Dmsqst1*Mst1*s2t + (3*Mt*(3891491 + 27200*lmMst1 - 960*
674  lmMt - 19200*pow2(lmMst1)) + 100*Mst1*s2t*(-32829 + 1852*lmMst1 + 660*
675  pow2(lmMst1)))*pow2(Msq)) + 40*Dmglst1*Mgl*(6000*Dmsqst1*Mst1*s2t + (5*
676  Mt*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1)) + 6*Mst1*s2t*(-
677  1282471 + 7264*lmMst1 + 18120*pow2(lmMst1)))*pow2(Msq)) + 15*pow2(Mgl)*
678  (67200*Dmsqst1*Mst1*s2t + (80*Mst1*s2t*(-36863 + 80*lmMst1 + 552*pow2(
679  lmMst1)) + Mt*(-1763661 + 47104*lmMst1 - 5120*lmMt + 24576*pow2(lmMst1)
680  ))*pow2(Msq)))*pow2(MuSUSY)*pow3(Mt) + 3675*pow2(Mst1)*(-(pow2(Mt)*
681  pow2(s2t)*((2400*Dmsqst1*(7 + 24*lmMst1)*pow2(Mgl) + (40*Dmglst1*Mgl*(-
682  84209 - 1264*lmMst1 + 240*pow2(lmMst1)) + 8*pow2(Dmglst1)*(-1732531 -
683  16896*lmMst1 + 24840*pow2(lmMst1)) + 5*(70121 + 2208*lmMst1 - 432*pow2(
684  lmMst1))*pow2(Mgl))*pow2(Msq))*pow2(MuSUSY) - 1440*shiftst1*pow2(Mgl)*(
685  80*Dmsqst1*(3 - 2*lmMst1)*pow2(Mst1) + (1 - 2*lmMst1)*pow2(Msq)*(80*
686  pow2(Mst1) + 3*pow2(MuSUSY)))*pow2(Sbeta))) + 144*pow2(Mgl)*(-(
687  shiftst3*pow2(Msq)*(2*pow2(Mt)*pow2(s2t)*((-29 + 18*lmMst1)*pow2(
688  MuSUSY)*(-1 + pow2(Sbeta)) + 5*(-83 + 58*lmMst1)*pow2(Mst1)*pow2(Sbeta)
689  ) - 80*(-7 + 3*lmMst1)*pow2(Sbeta)*pow4(Mt) + 5*(1 - 2*lmMst1)*pow2(
690  Sbeta)*pow4(Mst1)*pow4(s2t))) - 10*(3*(1 - 2*lmMst1)*shiftst1*pow2(Msq)
691  *pow2(Mt)*pow2(MuSUSY)*pow2(s2t) - (1 - 2*lmMst1)*shiftst2*pow2(Msq)*
692  pow2(s2t)*(2*pow2(Mt)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) - 10*pow2(Mst1)*
693  pow2(Sbeta)) + 5*pow2(s2t)*pow2(Sbeta)*pow4(Mst1)) + 5*shiftst1*(
694  Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*pow2(Sbeta)*pow4(
695  Mst1)*pow4(s2t) + 5*Dmsqst1*(3 - 2*lmMst1)*shiftst2*(4*pow2(Mt)*pow2(
696  s2t)*(-(pow2(MuSUSY)*(-1 + pow2(Sbeta))) + pow2(Mst1)*pow2(Sbeta)) +
697  pow2(Sbeta)*(24*pow4(Mt) - pow4(Mst1)*pow4(s2t)))))) + Mt*pow2(Sbeta)*(
698  -392*Mst1*s2t*pow2(Mt)*(15*pow2(Mgl)*(2000*Dmsqst1*(2*(13 + 6*lmMst1 -
699  6*lmMt)*pow2(Mst1) - 21*pow2(MuSUSY)) - pow2(Msq)*((558619 + 76160*lmMt
700  - 224*lmMst1*(1219 + 60*lmMt) + 123840*pow2(lmMst1) + 86400*pow2(lmMt))
701  *pow2(Mst1) + 50*(-36863 + 80*lmMst1 + 552*pow2(lmMst1))*pow2(MuSUSY)))
702  + 10*pow2(Dmglst1)*(3000*Dmsqst1*(84*pow2(Mst1) - 5*pow2(MuSUSY)) +
703  pow2(Msq)*((148185343 + 1333716*lmMst1 + 170460*lmMt + 87840*lmMst1*
704  lmMt - 1376640*pow2(lmMst1) - 43200*pow2(lmMt))*pow2(Mst1) - 150*(-
705  32829 + 1852*lmMst1 + 660*pow2(lmMst1))*pow2(MuSUSY))) - 2*Dmglst1*Mgl*
706  (600*Dmsqst1*(6*(47 - 30*lmMst1 + 30*lmMt)*pow2(Mst1) + 125*pow2(
707  MuSUSY)) + pow2(Msq)*((28188929 - 143100*lmMt - 3780*lmMst1*(549 + 80*
708  lmMt) + 1389600*pow2(lmMst1) + 388800*pow2(lmMt))*pow2(Mst1) + 75*(-
709  1282471 + 7264*lmMst1 + 18120*pow2(lmMst1))*pow2(MuSUSY)))) - 3675*Mt*
710  pow2(Mst1)*(8*pow2(Dmglst1)*(114960*Dmsqst1*pow2(Mst1) + pow2(Msq)*(3*(
711  -612347 - 7436*lmMst1 - 608*lmMt + 19320*pow2(lmMst1))*pow2(Mst1) + (
712  1732531 + 16896*lmMst1 - 24840*pow2(lmMst1))*pow2(MuSUSY))) - 5*pow2(
713  Mgl)*(120*Dmsqst1*((-137 + 288*lmMst1)*pow2(Mst1) + 4*(7 + 24*lmMst1)*
714  pow2(MuSUSY)) + pow2(Msq)*(24*(2785 - 304*lmMst1 + 384*lmMt + 768*pow2(
715  lmMst1))*pow2(Mst1) + (70121 + 2208*lmMst1 - 432*pow2(lmMst1))*pow2(
716  MuSUSY))) + 8*Dmglst1*Mgl*(112800*Dmsqst1*pow2(Mst1) + pow2(Msq)*(24*(-
717  20017 + 1203*lmMst1 - 200*lmMt + 2250*pow2(lmMst1))*pow2(Mst1) + 5*(
718  84209 + 1264*lmMst1 - 240*pow2(lmMst1))*pow2(MuSUSY))))*pow2(s2t) - (
719  40*pow2(Dmglst1)*(7056*Dmsqst1*(13147 - 90*lmMst1 + 90*lmMt)*pow2(Mst1)
720  + 2*(-700000759 + 6327384*lmMt + 12*lmMst1*(-88589 + 185010*lmMt) +
721  85997520*pow2(lmMst1) + 423360*pow2(lmMt))*pow2(Msq)*pow2(Mst1) + 441*(
722  -3891491 - 27200*lmMst1 + 960*lmMt + 19200*pow2(lmMst1))*pow2(Msq)*
723  pow2(MuSUSY)) + 49*pow2(Mgl)*(9600*Dmsqst1*(3268 + 2805*lmMst1 - 105*
724  lmMt)*pow2(Mst1) + pow2(Msq)*((83430364 - 8607840*lmMt + 480*lmMst1*(
725  36107 + 13380*lmMt) - 9273600*pow2(lmMst1) - 6652800*pow2(lmMt))*pow2(
726  Mst1) + 75*(1763661 - 47104*lmMst1 + 5120*lmMt - 24576*pow2(lmMst1))*
727  pow2(MuSUSY))) + 392*Dmglst1*Mgl*(21600*Dmsqst1*(423 - 20*lmMst1 + 20*
728  lmMt)*pow2(Mst1) + pow2(Msq)*(12*(9598037 + 92280*lmMt + 20*lmMst1*(-
729  11207 + 270*lmMt) + 246000*pow2(lmMst1) - 14400*pow2(lmMt))*pow2(Mst1)
730  - 125*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))*pow2(MuSUSY)))
731  )*pow3(Mt) - 29400*(pow2(Dmglst1)*(1000*Dmsqst1 + (1286791 + 5936*
732  lmMst1 - 18120*pow2(lmMst1))*pow2(Msq)) + 10*Dmglst1*Mgl*(100*Dmsqst1 +
733  (33261 - 532*lmMst1 - 660*pow2(lmMst1))*pow2(Msq)) + 20*pow2(Mgl)*(210*
734  Dmsqst1 + (1361 + 10*lmMst1 + 54*pow2(lmMst1))*pow2(Msq)))*pow3(s2t)*
735  pow5(Mst1))))) + Mt*(992250*oneLoopFlag*pow2(Mst1)*pow3(Mgl)*(-(Mt*
736  Tbeta*pow2(s2t)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 18*pow2(Mst1)*pow2(
737  Sbeta))) - pow2(Sbeta)*(-8*MuSUSY*s2t*pow2(Mt) + 8*Tbeta*pow3(Mt) +
738  MuSUSY*pow2(Mst1)*pow3(s2t)))*pow4(Msq) + 4*Al4p*xDmglst1*pow2(Msq)*
739  pow3(Dmglst1)*(-17640*Mst1*twoLoopFlag*pow2(Msq)*(-5*Mst1*Mt*pow2(s2t)*
740  (20*(-13 + 12*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*((
741  129 - 30*lmMst1)*MuSUSY + 4*(-83 + 60*lmMst1)*Mst1*Tbeta)*pow2(Sbeta))
742  + 2*s2t*pow2(Mt)*(-15*(9 + 10*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(
743  Sbeta)) + 2*Mst1*(30*(47 - 60*lmMst1)*MuSUSY + (106 + 285*lmMst1 + 15*
744  lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-2*((1577 + 270*lmMst1 +
745  30*lmMt)*MuSUSY + 12*(199 - 310*lmMst1 + 10*lmMt)*Mst1*Tbeta)*pow3(Mt)
746  + 5*((130 - 120*lmMst1)*MuSUSY + (77 - 30*lmMst1)*Mst1*Tbeta)*pow3(
747  Mst1)*pow3(s2t))) + Al4p*threeLoopFlag*(11760*Dmsqst1*Mst1*(-2*s2t*
748  pow2(Mt)*(625*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + (36600*Mst1*
749  MuSUSY - 2*(8543 + 390*lmMst1 - 390*lmMt)*Tbeta*pow2(Mst1))*pow2(Sbeta)
750  ) + pow2(Sbeta)*(150*Mt*(-25*MuSUSY + 488*Mst1*Tbeta)*pow2(Mst1)*pow2(
751  s2t) - 4*(5*(557 + 120*lmMst1 - 120*lmMt)*MuSUSY - 6*(3401 + 105*lmMst1
752  - 105*lmMt)*Mst1*Tbeta)*pow3(Mt) + 625*Tbeta*pow3(s2t)*pow4(Mst1))) +
753  pow2(Msq)*(-4*Mst1*MuSUSY*pow2(Mt)*(-245*MuSUSY*s2t*Tbeta*(-37824007 -
754  770520*lmMst1 + 131400*pow2(lmMst1))*(-1 + pow2(Sbeta)) + 8*Mt*(
755  9144246058 + 12254445*lmMt + 90*lmMst1*(1109907 + 18305*lmMt) -
756  48239100*pow2(lmMst1) - 264600*pow2(lmMt))*pow2(Sbeta)) + 6*Mt*pow2(
757  Mst1)*(-1225*Tbeta*(-3044017 - 27472*lmMst1 + 48480*pow2(lmMst1))*pow2(
758  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 8*Mt*(2450*MuSUSY*s2t*(580211 +
759  498*lmMst1 + 528*lmMt - 21060*pow2(lmMst1)) + Mt*Tbeta*(-874574719 +
760  9416610*lmMt + 10*lmMst1*(1016017 + 174300*lmMt) + 51500400*pow2(
761  lmMst1) + 705600*pow2(lmMt)))*pow2(Sbeta)) + s2t*(3675*Mst1*s2t*(
762  MuSUSY*s2t*(3044017 + 27472*lmMst1 - 48480*pow2(lmMst1)) + 4*Mt*Tbeta*(
763  486671 + 31944*lmMst1 - 192*lmMt - 15120*pow2(lmMst1))) + Mt*(735*
764  MuSUSY*s2t*(-223974673 - 2515800*lmMst1 + 1638000*pow2(lmMst1)) + 4*Mt*
765  Tbeta*(137797425107 + 35209020*lmMt + 300*lmMst1*(3595111 + 92568*lmMt)
766  - 690681600*pow2(lmMst1) - 2116800*pow2(lmMt))))*pow2(Sbeta)*pow3(Mst1)
767  + 490*Tbeta*((-59957863 - 480000*lmMst1 + 26880*lmMt + 230400*pow2(
768  lmMst1))*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 2*(31025111 +
769  290880*lmMst1 - 251100*pow2(lmMst1))*pow2(Sbeta)*pow3(s2t)*pow5(Mst1)))
770  )))))/(1.90512e8*Tbeta*pow2(Mst1)*pow2(Sbeta)*pow3(Mgl)*pow4(Msq)*pow6(
771  Mst2));
772 }
773 
774 /**
775  * @return The off-diagonal (1, 2) = (2, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3'
776  */
777 double H3::getS12() const {
778  return -(MuSUSY*(Mt*(1984500*Dmst12*oneLoopFlag*s2t*pow2(Mst1)*pow3(Mgl)*pow4(
779  Msq)*(-2*Dmst12*Mt*(MuSUSY*s2t + 6*Mt*Tbeta)*pow2(Mst2) + pow2(Dmst12)*
780  (2*Mt*MuSUSY*s2t + 8*Tbeta*pow2(Mt) - Tbeta*pow2(Mst1)*pow2(s2t)) + 24*
781  Tbeta*pow2(Mt)*pow4(Mst2)) + 8*Al4p*xDmglst1*pow2(Msq)*pow3(Dmglst1)*(
782  17640*Mst1*twoLoopFlag*pow2(Msq)*(Mt*pow2(Dmst12)*pow2(Mst2)*(80*Mt*
783  s2t*(13*MuSUSY - 18*Mst1*Tbeta) + 8*(-476 + 90*lmMst1 - 15*lmMt)*Tbeta*
784  pow2(Mt) + 5*Mst1*(20*(-13 + 12*lmMst1)*MuSUSY + 3*(-77 + 30*lmMst1)*
785  Mst1*Tbeta)*pow2(s2t)) + pow3(Dmst12)*(-60*s2t*((9 + 10*lmMst1)*MuSUSY
786  + 2*(47 - 60*lmMst1)*Mst1*Tbeta)*pow2(Mt) - 5*Mst1*Mt*(40*(-13 + 12*
787  lmMst1)*MuSUSY + 9*(-43 + 10*lmMst1)*Mst1*Tbeta)*pow2(s2t) + 2*(1577 +
788  270*lmMst1 + 30*lmMt)*Tbeta*pow3(Mt) + 50*(-13 + 12*lmMst1)*Tbeta*pow3(
789  Mst1)*pow3(s2t)) - 2*Dmst12*(10*(77 - 30*lmMst1)*MuSUSY*s2t + (-2231 +
790  990*lmMst1 - 90*lmMt)*Mt*Tbeta + 60*(-71 + 60*lmMst1)*Mst1*s2t*Tbeta)*
791  pow2(Mt)*pow4(Mst2) + 8*(977 - 480*lmMst1 + 30*lmMt)*Tbeta*pow3(Mt)*
792  pow6(Mst2)) + Al4p*threeLoopFlag*(490*Mt*MuSUSY*(-(pow2(Dmst12)*pow2(
793  Mst2)*(60000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(2*Mst1*Mt*s2t*(99874229 +
794  1352280*lmMst1 - 633600*pow2(lmMst1)) + (119915726 + 960000*lmMst1 -
795  53760*lmMt - 460800*pow2(lmMst1))*pow2(Mt) + 15*(-3044017 - 27472*
796  lmMst1 + 48480*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))) + 2*((30000*
797  Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(2*Mst1*Mt*s2t*(37824007 + 770520*
798  lmMst1 - 131400*pow2(lmMst1)) + (59957863 + 480000*lmMst1 - 26880*lmMt
799  - 230400*pow2(lmMst1))*pow2(Mt) + 15*(-3044017 - 27472*lmMst1 + 48480*
800  pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))*pow3(Dmst12) + Dmst12*Mt*(30000*
801  Dmsqst1*Mst1*s2t + (4*Mst1*s2t*(31025111 + 290880*lmMst1 - 251100*pow2(
802  lmMst1)) + Mt*(59957863 + 480000*lmMst1 - 26880*lmMt - 230400*pow2(
803  lmMst1)))*pow2(Msq))*pow4(Mst2)) + 24*(3877891 + 46400*lmMst1 - 960*
804  lmMt - 19200*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow6(Mst2)) - Mst1*Tbeta*
805  (735*Mt*pow2(Dmst12)*pow2(Mst1)*(30000*Dmsqst1*(2*Dmst12 - pow2(Mst2))
806  + pow2(Msq)*(Dmst12*(223974673 + 2515800*lmMst1 - 1638000*pow2(lmMst1))
807  + 4*(-31025111 - 290880*lmMst1 + 251100*pow2(lmMst1))*pow2(Mst2)))*
808  pow2(s2t) + 3675*(-3044017 - 27472*lmMst1 + 48480*pow2(lmMst1))*pow2(
809  Msq)*pow3(Dmst12)*pow3(Mst1)*pow3(s2t) + 29400*Dmst12*Mst1*s2t*pow2(Mt)
810  *(29280*Dmsqst1*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(
811  Msq)*(4*pow2(Dmst12)*(-580211 - 498*lmMst1 - 528*lmMt + 21060*pow2(
812  lmMst1)) + Dmst12*(611423 - 9984*lmMst1 + 768*lmMt - 23040*pow2(lmMst1)
813  )*pow2(Mst2) + 2*(548999 + 10980*lmMst1 + 288*lmMt - 19080*pow2(lmMst1)
814  )*pow4(Mst2))) + 4*pow3(Mt)*(11760*Dmsqst1*(5*Dmst12*(557 + 120*lmMst1
815  - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + 6*(-463 +
816  135*lmMst1 - 135*lmMt)*pow6(Mst2)) + pow2(Msq)*(-15*pow2(Dmst12)*(
817  5753390765 + 580*lmMst1*(79969 - 1932*lmMt) + 7091364*lmMt - 35700000*
818  pow2(lmMst1) - 282240*pow2(lmMt))*pow2(Mst2) + 8*(9144246058 +
819  12254445*lmMt + 90*lmMst1*(1109907 + 18305*lmMt) - 48239100*pow2(
820  lmMst1) - 264600*pow2(lmMt))*pow3(Dmst12) + 2*Dmst12*(49723877243 + 60*
821  lmMst1*(4936063 - 389970*lmMt) + 57352680*lmMt - 342543600*pow2(lmMst1)
822  - 3175200*pow2(lmMt))*pow4(Mst2) + 392*(122282257 + 60*lmMst1*(8318 -
823  3885*lmMt) + 479550*lmMt - 1351800*pow2(lmMst1) - 21600*pow2(lmMt))*
824  pow6(Mst2))))))) + 1470*Al4p*Mgl*pow2(Msq)*(48*Mst1*Mt*twoLoopFlag*
825  pow2(Msq)*(-25*s2t*pow3(Dmst12)*(8*MuSUSY*(Dmglst1*(-5 + 6*lmMst1)*Mgl
826  + (1 + 6*lmMst1)*pow2(Dmglst1) + (5 + 6*lmMst1)*pow2(Mgl))*pow2(Mt) +
827  Tbeta*(8*Dmglst1*(1 - 3*lmMst1)*Mgl + 6*(5 - 6*lmMst1)*pow2(Dmglst1) +
828  3*pow2(Mgl))*pow2(s2t)*pow3(Mst1)) + Mt*(50*Dmst12*MuSUSY*s2t*(Mst1*
829  s2t*pow2(Dmst12)*(16*Dmglst1*(1 - 3*lmMst1)*Mgl + (60 - 72*lmMst1)*
830  pow2(Dmglst1) + (1 + 6*lmMst1)*pow2(Mgl)) + 4*Dmst12*Mt*(6*pow2(
831  Dmglst1) + pow2(Mgl))*pow2(Mst2) - Dmst12*Mst1*s2t*(8*Dmglst1*(1 - 3*
832  lmMst1)*Mgl + 6*(5 - 6*lmMst1)*pow2(Dmglst1) + 3*pow2(Mgl))*pow2(Mst2)
833  - 4*Mt*(Dmglst1*(5 - 6*lmMst1)*Mgl + (11 - 6*lmMst1)*pow2(Dmglst1) - 3*
834  (1 + 2*lmMst1)*pow2(Mgl))*pow4(Mst2)) + 2*Tbeta*(3*pow2(Dmglst1)*(-(
835  pow2(Dmst12)*pow2(Mst2)*(70*Mst1*Mt*s2t + 8*(92 - 30*lmMst1 + 5*lmMt)*
836  pow2(Mt) + 25*(11 - 6*lmMst1)*pow2(Mst1)*pow2(s2t))) + (40*(-29 + 45*
837  lmMst1)*Mst1*Mt*s2t + 2*(259 + 90*lmMst1 + 10*lmMt)*pow2(Mt) + 25*(17 -
838  6*lmMst1)*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + 2*Dmst12*Mt*((477 - 330*
839  lmMst1 + 30*lmMt)*Mt + 50*(13 - 18*lmMst1)*Mst1*s2t)*pow4(Mst2) + 400*(
840  4 - 3*lmMst1)*pow2(Mt)*pow6(Mst2)) - 25*pow2(Mgl)*(-(pow2(Dmst12)*pow2(
841  Mst2)*(24*(1 - 3*lmMst1)*Mst1*Mt*s2t + 2*(5 + 6*lmMst1 + 6*lmMt)*pow2(
842  Mt) + 9*(1 + 2*lmMst1)*pow2(Mst1)*pow2(s2t))) + (-48*(-1 + lmMst1)*
843  Mst1*Mt*s2t - 4*(5 + 6*lmMst1)*pow2(Mt) + 6*(1 + 3*lmMst1)*pow2(Mst1)*
844  pow2(s2t))*pow3(Dmst12) + 4*Dmst12*Mt*(2*(5 + 6*lmMst1 + 3*lmMt)*Mt -
845  9*(1 + 2*lmMst1)*Mst1*s2t)*pow4(Mst2) + 72*(1 + lmMst1 + lmMt)*pow2(Mt)
846  *pow6(Mst2)) + Dmglst1*Mgl*(-3*pow2(Dmst12)*(pow2(Mst2)*(100*Mst1*Mt*
847  s2t - 2*(-41 + 90*lmMst1 + 10*lmMt)*pow2(Mt) + 25*(5 - 6*lmMst1)*pow2(
848  Mst1)*pow2(s2t)) + Dmst12*(-1200*lmMst1*Mst1*Mt*s2t + 8*(17 - 30*lmMst1
849  + 5*lmMt)*pow2(Mt) + 25*(-5 + 6*lmMst1)*pow2(Mst1)*pow2(s2t))) + 300*
850  Dmst12*Mt*((3 - 6*lmMst1)*Mt + 2*(1 - 6*lmMst1)*Mst1*s2t)*pow4(Mst2) -
851  200*(-7 + 15*lmMst1 + 3*lmMt)*pow2(Mt)*pow6(Mst2))))) + Al4p*
852  threeLoopFlag*((MuSUSY*pow2(Mt)*(8*pow2(Dmglst1)*(-3*pow2(Dmst12)*pow2(
853  Mst2)*(20000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(10*Mst1*Mt*s2t*(958501 +
854  24456*lmMst1 - 11520*pow2(lmMst1)) + 6*(3891491 + 27200*lmMst1 - 960*
855  lmMt - 19200*pow2(lmMst1))*pow2(Mt) + 5*(-1732531 - 16896*lmMst1 +
856  24840*pow2(lmMst1))*pow2(Mst1)*pow2(s2t))) + 6*((10000*Dmsqst1*Mst1*Mt*
857  s2t + pow2(Msq)*(100*Mst1*Mt*s2t*(-32829 + 1852*lmMst1 + 660*pow2(
858  lmMst1)) + 3*(3891491 + 27200*lmMst1 - 960*lmMt - 19200*pow2(lmMst1))*
859  pow2(Mt) + 5*(-1732531 - 16896*lmMst1 + 24840*pow2(lmMst1))*pow2(Mst1)*
860  pow2(s2t)))*pow3(Dmst12) + Dmst12*Mt*(10000*Dmsqst1*Mst1*s2t + (3*Mt*(
861  3891491 + 27200*lmMst1 - 960*lmMt - 19200*pow2(lmMst1)) + 10*Mst1*s2t*(
862  1286791 + 5936*lmMst1 - 18120*pow2(lmMst1)))*pow2(Msq))*pow4(Mst2)) +
863  100*(345581 + 4896*lmMst1 + 96*lmMt - 3456*pow2(lmMst1))*pow2(Msq)*
864  pow2(Mt)*pow6(Mst2)) + 40*Dmglst1*Mgl*(-(pow2(Dmst12)*pow2(Mst2)*(
865  12000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(6*Mst1*Mt*s2t*(-949861 + 1944*
866  lmMst1 + 11520*pow2(lmMst1)) + 10*(403559 + 384*(lmMst1 + lmMt) - 4608*
867  pow2(lmMst1))*pow2(Mt) + 15*(-84209 - 1264*lmMst1 + 240*pow2(lmMst1))*
868  pow2(Mst1)*pow2(s2t)))) + 2*(6000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(6*
869  Mst1*Mt*s2t*(-1282471 + 7264*lmMst1 + 18120*pow2(lmMst1)) + 5*(403559 +
870  384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))*pow2(Mt) + 15*(-84209 - 1264*
871  lmMst1 + 240*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))*pow3(Dmst12) + 10*
872  Dmst12*Mt*(1200*Dmsqst1*Mst1*s2t + (Mt*(403559 + 384*(lmMst1 + lmMt) -
873  4608*pow2(lmMst1)) - 12*Mst1*s2t*(-33261 + 532*lmMst1 + 660*pow2(
874  lmMst1)))*pow2(Msq))*pow4(Mst2) + 240*(9631 + 16*lmMst1 + 48*lmMt -
875  192*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow6(Mst2)) + 15*pow2(Mgl)*(-(
876  pow2(Dmst12)*pow2(Mst2)*(2400*Dmsqst1*Mst1*s2t*(56*Mt + Mst1*s2t*(7 -
877  12*lmMst1*(-2 + shiftst1 + shiftst2) + 18*(shiftst1 + shiftst2))) -
878  pow2(Msq)*(-240*Mst1*Mt*s2t*(-10473 + 40*lmMst1 + 256*pow2(lmMst1)) + (
879  852541 + 9216*lmMst1 - 10240*lmMt + 6144*pow2(lmMst1))*pow2(Mt) + 80*(
880  1429 - 454*lmMst1 - 180*(shiftst1 + shiftst2) + 360*lmMst1*(shiftst1 +
881  shiftst2) - 126*shiftst3 + 108*lmMst1*shiftst3 + 24*pow2(lmMst1))*pow2(
882  Mst1)*pow2(s2t)))) + 2*(-2400*Dmsqst1*Mst1*s2t*(-28*Mt + Mst1*s2t*(-7 +
883  24*lmMst1*(-1 + shiftst2) - 36*shiftst2)) + pow2(Msq)*(80*Mst1*Mt*s2t*(
884  -36863 + 80*lmMst1 + 552*pow2(lmMst1)) + (-1763661 + 47104*lmMst1 -
885  5120*lmMt + 24576*pow2(lmMst1))*pow2(Mt) + (350605 + 4320*shiftst1 +
886  2880*shiftst2 + 8352*shiftst3 - 96*lmMst1*(-115 + 90*shiftst1 + 60*
887  shiftst2 + 54*shiftst3) - 2160*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))*
888  pow3(Dmst12) + 160*Dmst12*(60*Dmsqst1*Mst1*s2t*(14*Mt + 3*(3 - 2*
889  lmMst1)*Mst1*s2t*(shiftst1 - shiftst2)) + pow2(Msq)*(4*Mst1*Mt*s2t*(
890  1361 + 10*lmMst1 + 54*pow2(lmMst1)) + (11389 - 704*lmMst1 + 192*lmMt -
891  384*pow2(lmMst1))*pow2(Mt) + 18*(1 - 2*lmMst1)*(10*shiftst1 - 10*
892  shiftst2 + shiftst3)*pow2(Mst1)*pow2(s2t)))*pow4(Mst2) + 1920*(349 -
893  56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow6(Mst2))))
894  /3. - 4*Mst1*Tbeta*(15*Dmst12*Mst1*s2t*pow2(Mt)*(Dmst12*Mst1*s2t*(5*
895  pow2(Mgl)*(1680*Dmsqst1*(2*Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(-
896  20531 + 200*lmMst1 + 1200*pow2(lmMst1)) - 8*(1361 + 10*lmMst1 + 54*
897  pow2(lmMst1))*pow2(Mst2))) + Dmglst1*Mgl*(2000*Dmsqst1*(2*Dmst12 -
898  pow2(Mst2)) - pow2(Msq)*(Dmst12*(284641 + 8696*lmMst1 + 1680*pow2(
899  lmMst1)) - 20*(-33261 + 532*lmMst1 + 660*pow2(lmMst1))*pow2(Mst2))) +
900  pow2(Dmglst1)*(2000*Dmsqst1*(2*Dmst12 - pow2(Mst2)) + pow2(Msq)*(
901  Dmst12*(3532083 + 36328*lmMst1 - 47760*pow2(lmMst1)) + 2*(-1286791 -
902  5936*lmMst1 + 18120*pow2(lmMst1))*pow2(Mst2)))) + Mt*(pow2(Mgl)*(100*
903  Dmsqst1*((346 + 288*lmMst1)*pow2(Dmst12) - 161*Dmst12*pow2(Mst2) - 24*(
904  1 + 12*lmMst1)*pow4(Mst2)) + pow2(Msq)*(pow2(Dmst12)*(102747 + 10592*
905  lmMst1 - 640*lmMt - 13888*pow2(lmMst1)) + 20*Dmst12*(-2071 - 296*lmMst1
906  + 96*lmMt + 600*pow2(lmMst1))*pow2(Mst2) - 160*(631 - lmMst1 + 36*lmMt
907  + 21*pow2(lmMst1))*pow4(Mst2))) + 4*pow2(Dmglst1)*(19160*Dmsqst1*(pow2(
908  Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(Msq)*(8*pow2(Dmst12)*(
909  -143196 - 2546*lmMst1 - 92*lmMt + 4785*pow2(lmMst1)) + Dmst12*(586073 +
910  9268*lmMst1 + 448*lmMt - 19200*pow2(lmMst1))*pow2(Mst2) + 2*(-13289 +
911  916*lmMst1 - 80*lmMt + 60*pow2(lmMst1))*pow4(Mst2))) + 16*Dmglst1*Mgl*(
912  4700*Dmsqst1*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(
913  Msq)*(pow2(Dmst12)*(-12383 - 4128*lmMst1 - 80*lmMt + 1260*pow2(lmMst1))
914  + Dmst12*(17539 + 574*lmMst1 + 160*lmMt - 1920*pow2(lmMst1))*pow2(Mst2)
915  + 5*(-4539 + 596*lmMst1 - 48*lmMt + 516*pow2(lmMst1))*pow4(Mst2))))) +
916  900*Mst1*Mt*s2t*pow2(Mgl)*(-10*shiftst2*pow2(Dmst12)*(Dmsqst1*(3 - 2*
917  lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*pow2(Mst1)*(2*Dmst12 - 3*pow2(Mst2)
918  )*pow2(s2t) + shiftst3*pow2(Dmst12)*pow2(Msq)*pow2(Mst1)*(Dmst12*(13 -
919  14*lmMst1) + 3*(-1 + 2*lmMst1)*pow2(Mst2))*pow2(s2t) - 240*shiftst2*
920  pow2(Mt)*((1 - 2*lmMst1)*pow2(Msq)*(Dmst12 + pow2(Mst2))*pow4(Mst2) +
921  Dmsqst1*(-3 + 2*lmMst1)*(pow3(Dmst12) - Dmst12*pow4(Mst2) - pow6(Mst2))
922  ) - 8*shiftst3*pow2(Msq)*pow2(Mt)*(6*(-2 + lmMst1)*pow2(Dmst12)*pow2(
923  Mst2) + (14 - 6*lmMst1)*pow3(Dmst12) + 3*Dmst12*(3 - 2*lmMst1)*pow4(
924  Mst2) + 3*(-1 + 2*lmMst1)*pow6(Mst2))) + Mt*(50*(30*Dmsqst1*(7 + 24*
925  lmMst1)*pow2(Mgl) - ((pow2(Dmglst1)*(1732531 + 16896*lmMst1 - 24840*
926  pow2(lmMst1)) + 5*Dmglst1*Mgl*(84209 + 1264*lmMst1 - 240*pow2(lmMst1))
927  + 10*(1429 - 454*lmMst1 + 24*pow2(lmMst1))*pow2(Mgl))*pow2(Msq))/10.)*
928  pow3(Dmst12)*pow3(Mst1)*pow3(s2t) + 9000*Mst1*s2t*shiftst1*pow2(Mgl)*(
929  Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*(pow2(Dmst12)*pow2(
930  Mst1)*(4*Dmst12 - 3*pow2(Mst2))*pow2(s2t) + 24*pow2(Mt)*pow6(Mst2))) +
931  (2*pow4(Mt)*(2*Dmglst1*(1200*Dmglst1*Dmsqst1*(Dmst12*(557 + 120*lmMst1
932  - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + 16*(4 +
933  15*lmMst1 - 15*lmMt)*pow6(Mst2)) + 1200*Dmsqst1*Mgl*(Dmst12*(557 + 120*
934  lmMst1 - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) +
935  50*(26 + 3*lmMst1 - 3*lmMt)*pow6(Mst2)) + Dmglst1*pow2(Msq)*(3*pow2(
936  Dmst12)*(-129193181 - 401100*lmMt + 100*lmMst1*(-7351 + 1164*lmMt) +
937  1336800*pow2(lmMst1) + 28800*pow2(lmMt))*pow2(Mst2) + 4*(29818901 +
938  258300*lmMt + 30*lmMst1*(35963 + 2190*lmMt) - 239400*pow2(lmMst1) -
939  10800*pow2(lmMt))*pow3(Dmst12) + 2*Dmst12*(327941741 + 1080*lmMst1*(44
940  - 445*lmMt) + 686700*lmMt - 3531600*pow2(lmMst1) - 64800*pow2(lmMt))*
941  pow4(Mst2) + 80*(3777727 + 24255*lmMt - 486*lmMst1*(118 + 45*lmMt) -
942  52380*pow2(lmMst1))*pow6(Mst2)) - Mgl*pow2(Msq)*(3*pow2(Dmst12)*(-
943  39474953 + lmMst1*(52780 - 87600*lmMt) - 3300*lmMt + 319200*pow2(
944  lmMst1) + 14400*pow2(lmMt))*pow2(Mst2) - 2*(-193364399 - 90000*lmMt +
945  300*lmMst1*(3781 + 582*lmMt) + 2005200*pow2(lmMst1) + 43200*pow2(lmMt))
946  *pow3(Dmst12) + 40*Dmst12*(-3746977 - 4005*lmMt + 18*lmMst1*(2711 +
947  1215*lmMt) + 52380*pow2(lmMst1))*pow4(Mst2) - 6000*(9961 + 66*lmMt -
948  10*lmMst1*(67 + 24*lmMt) + 42*pow2(lmMst1) + 72*pow2(lmMt))*pow6(Mst2))
949  ) + 15*pow2(Mgl)*(4000*Dmsqst1*(Dmst12*(65 - 6*lmMst1 + 6*lmMt)*(pow2(
950  Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + 6*(20 - 3*lmMst1 + 3*lmMt)*
951  pow6(Mst2)) - pow2(Msq)*(pow2(Dmst12)*(-2284899 + 49840*lmMt - 32*
952  lmMst1*(-1793 + 555*lmMt) + 87360*pow2(lmMst1) + 28800*pow2(lmMt))*
953  pow2(Mst2) - 2*(-3454599 + 16840*lmMt + 48*lmMst1*(262 + 405*lmMt) +
954  46560*pow2(lmMst1))*pow3(Dmst12) - 200*Dmst12*(11697 + 448*lmMst1 +
955  330*lmMt - 372*lmMst1*lmMt + 408*pow2(lmMst1) + 288*pow2(lmMt))*pow4(
956  Mst2) - 1600*(434 - 83*lmMst1 + 174*lmMt - 66*lmMst1*lmMt + 183*pow2(
957  lmMst1) + 108*pow2(lmMt))*pow6(Mst2)))))/15.))) + 735*threeLoopFlag*
958  pow2(Al4p)*(-2880*z2*pow2(Mst1)*pow3(Mgl)*(-50*Mt*s2t*xDmsqst1*pow2(
959  Dmsqst1)*(s2t*pow2(Dmst12)*(2*Mt*MuSUSY*(shiftst1 + shiftst2) + 3*s2t*(
960  -shiftst1 + shiftst2)*Tbeta*pow2(Mst1))*pow2(Mst2) + (-8*Mt*MuSUSY*s2t*
961  shiftst2 + 24*shiftst2*Tbeta*pow2(Mt) + 2*(2*shiftst1 - shiftst2)*
962  Tbeta*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) - 4*Dmst12*Mt*(MuSUSY*s2t*(
963  shiftst1 - shiftst2) + 6*Mt*shiftst2*Tbeta)*pow4(Mst2) + 24*(shiftst1 -
964  shiftst2)*Tbeta*pow2(Mt)*pow6(Mst2)) + pow2(Msq)*(2*Dmst12*MuSUSY*pow2(
965  Mt)*pow2(s2t)*(2*pow2(Dmst12)*(100*Dmsqst1*shiftst2 + (15*shiftst1 +
966  10*shiftst2 + 9*shiftst3)*pow2(Msq)) - 5*Dmst12*(10*Dmsqst1*(shiftst1 +
967  shiftst2) + (10*(shiftst1 + shiftst2) + 3*shiftst3)*pow2(Msq))*pow2(
968  Mst2) + 10*(10*Dmsqst1*(shiftst1 - shiftst2) + (10*shiftst1 - 10*
969  shiftst2 + shiftst3)*pow2(Msq))*pow4(Mst2)) + 5*shiftst3*Tbeta*pow2(
970  Msq)*(-(Mt*pow2(Dmst12)*pow2(Mst1)*(7*Dmst12 - 3*pow2(Mst2))*pow3(s2t))
971  + 24*s2t*pow3(Mt)*(-(pow2(Dmst12)*pow2(Mst2)) + pow3(Dmst12) + Dmst12*
972  pow4(Mst2) - pow6(Mst2))) - 50*Tbeta*(shiftst2*(-(Mt*pow2(Dmst12)*(
973  Dmsqst1 + pow2(Msq))*pow2(Mst1)*(2*Dmst12 - 3*pow2(Mst2))*pow3(s2t)) +
974  24*s2t*pow3(Mt)*(Dmsqst1*pow3(Dmst12) - (Dmsqst1 + pow2(Msq))*(Dmst12 +
975  pow2(Mst2))*pow4(Mst2))) + shiftst1*(Dmsqst1 + pow2(Msq))*(Mt*pow2(
976  Dmst12)*pow2(Mst1)*(4*Dmst12 - 3*pow2(Mst2))*pow3(s2t) + 24*s2t*pow3(
977  Mt)*pow6(Mst2))))) + 32*Mst1*xDmsqst1*pow2(Dmsqst1)*(4*xDmglst1*pow2(
978  Mt)*pow3(Dmglst1)*(2500*Dmst12*Mt*MuSUSY*s2t*(pow2(Dmst12) - Dmst12*
979  pow2(Mst2) + pow4(Mst2)) - Tbeta*(1875*pow2(Dmst12)*pow2(Mst1)*(2*
980  Dmst12 - pow2(Mst2))*pow2(s2t) + 73200*Dmst12*Mst1*Mt*s2t*(pow2(Dmst12)
981  - Dmst12*pow2(Mst2) + pow4(Mst2)) + 4*pow2(Mt)*(5*Dmst12*(557 + 120*
982  lmMst1 - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + (-
983  3778 + 435*lmMst1 - 435*lmMt)*pow6(Mst2)))) + Mgl*(250*Dmst12*MuSUSY*
984  s2t*pow2(Mt)*(40*Dmglst1*(Dmglst1 + Mgl)*Mt*(pow2(Dmst12) - Dmst12*
985  pow2(Mst2) + pow4(Mst2)) + 3*pow2(Mgl)*(2*(28*Mt + Mst1*s2t*(7 + 6*
986  shiftst1 - 24*lmMst1*(-1 + shiftst2) + 30*shiftst2))*pow2(Dmst12) -
987  Dmst12*(56*Mt + Mst1*s2t*(7 + 30*shiftst1 + 6*shiftst2 - 12*lmMst1*(-2
988  + shiftst1 + shiftst2)))*pow2(Mst2) + 8*(7*Mt - 3*(-2 + lmMst1)*Mst1*
989  s2t*(shiftst1 - shiftst2))*pow4(Mst2))) - 25*Tbeta*(60*pow2(Dmst12)*(5*
990  Dmglst1*Mgl + 5*pow2(Dmglst1) + 21*pow2(Mgl))*pow2(Mst1)*(2*Dmst12 -
991  pow2(Mst2))*pow2(Mt)*pow2(s2t) + 3*Dmst12*Mst1*s2t*pow3(Mt)*(8*Dmglst1*
992  (479*Dmglst1 + 470*Mgl)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2))
993  + 5*pow2(Mgl)*(3*(47 + 96*lmMst1)*pow2(Dmst12) + 44*Dmst12*pow2(Mst2) -
994  (229 + 288*lmMst1)*pow4(Mst2))) + 15*Mst1*Mt*s2t*pow2(Mgl)*(36*pow2(
995  Dmst12)*pow2(Mst2)*(4*shiftst2*pow2(Mt) + (-2 + lmMst1)*(shiftst1 -
996  shiftst2)*pow2(Mst1)*pow2(s2t)) + (-288*(-1 + lmMst1)*shiftst2*pow2(Mt)
997  + (7 + 108*shiftst1 - 72*shiftst2 + 24*lmMst1*(1 - 2*shiftst1 +
998  shiftst2))*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + 288*Dmst12*(-2 +
999  lmMst1)*shiftst2*pow2(Mt)*pow4(Mst2) - 288*(-2 + lmMst1)*(shiftst1 -
1000  shiftst2)*pow2(Mt)*pow6(Mst2)) + (16*pow4(Mt)*(pow2(Dmglst1)*(Dmst12*(
1001  557 + 120*lmMst1 - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(
1002  Mst2)) + (-136 + 165*lmMst1 - 165*lmMt)*pow6(Mst2)) + Dmglst1*Mgl*(
1003  Dmst12*(557 + 120*lmMst1 - 120*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2)
1004  + pow4(Mst2)) + 25*(44 + 3*lmMst1 - 3*lmMt)*pow6(Mst2)) + 25*pow2(Mgl)*
1005  (Dmst12*(65 - 6*lmMst1 + 6*lmMt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) +
1006  pow4(Mst2)) + (91 - 12*lmMst1 + 12*lmMt)*pow6(Mst2))))/5.))) - 75*Mt*
1007  z3*(-5040*Dmst12*s2t*xDmsqst1*pow2(Dmsqst1)*pow2(Mst1)*pow3(Mgl)*(2*
1008  Dmst12*Mt*MuSUSY*s2t*pow2(Mst2) + pow2(Dmst12)*(-4*Mt*MuSUSY*s2t + 15*
1009  Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*pow2(s2t)) - 15*Tbeta*pow2(Mt)*pow4(
1010  Mst2)) + 4*xDmglst1*pow3(Dmglst1)*pow4(Msq)*(-2*Mt*pow2(Dmst12)*pow2(
1011  Mst2)*(2*Mst1*Mt*s2t*(737647*MuSUSY + 137634*Mst1*Tbeta) + (884106*
1012  MuSUSY - 5205992*Mst1*Tbeta)*pow2(Mt) - (349745*MuSUSY + 1377324*Mst1*
1013  Tbeta)*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(8*Mst1*s2t*(278539*MuSUSY
1014  + 280806*Mst1*Tbeta)*pow2(Mt) - 11*Mt*(127180*MuSUSY + 451599*Mst1*
1015  Tbeta)*pow2(Mst1)*pow2(s2t) + 4*(442053*MuSUSY - 2202506*Mst1*Tbeta)*
1016  pow3(Mt) + 349745*Tbeta*pow3(s2t)*pow4(Mst1)) + 12*Dmst12*(Mt*(147351*
1017  MuSUSY - 1001162*Mst1*Tbeta) + 24*Mst1*s2t*(12753*MuSUSY - 3977*Mst1*
1018  Tbeta))*pow2(Mt)*pow4(Mst2) + 624*(2205*MuSUSY - 9298*Mst1*Tbeta)*pow3(
1019  Mt)*pow6(Mst2)) + Mgl*pow2(Msq)*(Mt*MuSUSY*(2*pow2(Dmst12)*pow2(Mst1)*(
1020  5040*Dmsqst1*pow2(Mgl)*(2*Dmst12 - pow2(Mst2)) - 20*Dmglst1*(40023*
1021  Dmglst1 + 10097*Mgl)*pow2(Msq)*(2*Dmst12 - pow2(Mst2)) + pow2(Mgl)*
1022  pow2(Msq)*(37669*Dmst12 + 10296*pow2(Mst2)))*pow2(s2t) + Mt*pow2(Msq)*(
1023  -16*Dmst12*Mst1*s2t*(2*pow2(Dmst12)*(142987*Dmglst1*Mgl + 37582*pow2(
1024  Dmglst1) + 20297*pow2(Mgl)) + 5*Dmst12*(-21081*Dmglst1*Mgl + 21081*
1025  pow2(Dmglst1) - 3457*pow2(Mgl))*pow2(Mst2) - 2*(37582*Dmglst1*Mgl +
1026  142987*pow2(Dmglst1) + 3012*pow2(Mgl))*pow4(Mst2)) + 3*Mt*(-(pow2(Mgl)*
1027  (-31963*pow2(Dmst12)*pow2(Mst2) + 131926*pow3(Dmst12) - 68000*Dmst12*
1028  pow4(Mst2) - 24064*pow6(Mst2))) + 432*pow2(Dmglst1)*(-3185*pow2(Dmst12)
1029  *pow2(Mst2) + 3185*Dmst12*(pow2(Dmst12) + pow4(Mst2)) + 1566*pow6(Mst2)
1030  ) + 48*Dmglst1*Mgl*(-8213*pow2(Dmst12)*pow2(Mst2) + 8213*Dmst12*(pow2(
1031  Dmst12) + pow4(Mst2)) + 4664*pow6(Mst2))))) + 4*Mst1*Tbeta*(-((1260*
1032  Dmsqst1*pow2(Mgl) - (50485*Dmglst1*Mgl + 200115*pow2(Dmglst1) + 2574*
1033  pow2(Mgl))*pow2(Msq))*pow3(Dmst12)*pow3(Mst1)*pow3(s2t)) - Dmst12*Mst1*
1034  s2t*pow2(Mt)*(1260*Dmsqst1*pow2(Mgl)*(22*pow2(Dmst12) - 7*Dmst12*pow2(
1035  Mst2) - 8*pow4(Mst2)) + pow2(Msq)*(pow2(Mgl)*(32783*pow2(Dmst12) -
1036  14508*Dmst12*pow2(Mst2) - 15840*pow4(Mst2)) - 24*Dmglst1*(Dmglst1*(
1037  68291*pow2(Dmst12) - 32498*Dmst12*pow2(Mst2) - 3295*pow4(Mst2)) + Mgl*(
1038  4744*pow2(Dmst12) - 3821*Dmst12*pow2(Mst2) + 2898*pow4(Mst2))))) +
1039  pow2(Msq)*(-3*Mt*pow2(Dmst12)*pow2(Mst1)*(Dmst12*(-30241*Dmglst1*Mgl +
1040  391379*pow2(Dmglst1) - 11261*pow2(Mgl)) - 2*(37582*Dmglst1*Mgl +
1041  142987*pow2(Dmglst1) + 3012*pow2(Mgl))*pow2(Mst2))*pow2(s2t) + 2*pow3(
1042  Mt)*(4*Dmglst1*Mgl*(-86833*pow2(Dmst12)*pow2(Mst2) + 287078*pow3(
1043  Dmst12) - 113412*Dmst12*pow4(Mst2) - 47112*pow6(Mst2)) + pow2(Mgl)*(-
1044  51181*pow2(Dmst12)*pow2(Mst2) + 152018*pow3(Dmst12) - 49656*Dmst12*
1045  pow4(Mst2) - 10176*pow6(Mst2)) - 4*pow2(Dmglst1)*(-286982*pow2(Dmst12)*
1046  pow2(Mst2) + 86737*pow3(Dmst12) + 487227*Dmst12*pow4(Mst2) + 226824*
1047  pow6(Mst2))))))))))/(1.90512e8*Tbeta*pow2(Mst1)*pow3(Mgl)*pow4(Msq)*
1048  pow6(Mst2));
1049 }
1050 
1051 /**
1052  * @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)
1053  */
1054 double H3::calc_coef_at_as2_no_sm_logs_log0() const {
1055 
1056  const double result =
1057  (49*pow3(Mgl)*(-150*Dmsqst1*pow2(Msq)*(-4*Mt*pow2(Dmst12)*pow2(Mst2)*(
1058  19600*Mst1*Mt*s2t + (-9122 + 14175*z3)*pow2(Mt) + 450*(2 - 21*z3)*pow2(
1059  Mst1)*pow2(s2t)) + pow3(Dmst12)*(-20800*Mst1*s2t*pow2(Mt) - 75*Mt*(274
1060  + 63*z3)*pow2(Mst1)*pow2(s2t) + 16*(-6536 + 14175*z3)*pow3(Mt) - 8400*
1061  pow3(Mst1)*pow3(s2t)) - 200*Dmst12*(-888*Mst1*s2t + Mt*(-158 + 567*z3))
1062  *pow2(Mt)*pow4(Mst2) - 201600*pow3(Mt)*pow6(Mst2)) - 150*pow2(Dmsqst1)*
1063  (-(Mt*pow2(Dmst12)*pow2(Mst2)*(36800*Mst1*Mt*s2t + 20112*pow2(Mt) + 75*
1064  (458 - 945*z3)*pow2(Mst1)*pow2(s2t))) + 3*pow3(Dmst12)*(-20800*Mst1*
1065  s2t*pow2(Mt) - 525*Mt*(-26 + 45*z3)*pow2(Mst1)*pow2(s2t) + 4*(-3998 +
1066  14175*z3)*pow3(Mt) - 2800*pow3(Mst1)*pow3(s2t)) - 100*Dmst12*(-1360*
1067  Mst1*s2t + 63*Mt*(-14 + 27*z3))*pow2(Mt)*pow4(Mst2) - 200*(410 + 567*
1068  z3)*pow3(Mt)*pow6(Mst2)) - pow4(Msq)*(-20*Mt*pow2(Dmst12)*pow2(Mst2)*(
1069  300*Mst1*Mt*s2t*(-17512 + 14805*z3) + (-375892 + 621675*z3)*pow2(Mt) -
1070  900*(-1226 + 495*z3)*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-30*Mst1*
1071  s2t*(-1117238 + 877575*z3)*pow2(Mt) - 2250*Mt*(-5570 + 333*z3)*pow2(
1072  Mst1)*pow2(s2t) + (-41715182 + 38174625*z3)*pow3(Mt) + 3000*(-2722 +
1073  2259*z3)*pow3(Mst1)*pow3(s2t)) - 6000*Dmst12*(81*Mt*(-406 + 285*z3) +
1074  8*Mst1*s2t*(-694 + 477*z3))*pow2(Mt)*pow4(Mst2) + 48000*(623 + 963*z3)*
1075  pow3(Mt)*pow6(Mst2))) + 196*Dmglst1*pow2(Mgl)*(600*pow2(Dmsqst1)*(-6*
1076  Mt*pow2(Dmst12)*pow2(Mst2)*(-347*Mst1*Mt*s2t + 2538*pow2(Mt) + 1175*
1077  pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(386*Mst1*s2t*pow2(Mt) + 14100*Mt*
1078  pow2(Mst1)*pow2(s2t) + 15228*pow3(Mt) + 125*pow3(Mst1)*pow3(s2t)) + 2*
1079  Dmst12*(7614*Mt - 2275*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 22600*pow3(Mt)*
1080  pow6(Mst2)) + 600*Dmsqst1*pow2(Msq)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*(-
1081  1516*Mst1*Mt*s2t + 7614*pow2(Mt) + 3525*pow2(Mst1)*pow2(s2t)) + pow3(
1082  Dmst12)*(-564*Mst1*s2t*pow2(Mt) + 14100*Mt*pow2(Mst1)*pow2(s2t) +
1083  15228*pow3(Mt) + 125*pow3(Mst1)*pow3(s2t)) + 4*Dmst12*(3807*Mt - 1375*
1084  Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 29200*pow3(Mt)*pow6(Mst2)) + pow4(Msq)*
1085  (Mt*pow2(Dmst12)*pow2(Mst2)*(40*Mst1*Mt*s2t*(-4511549 + 3729375*z3) + (
1086  45149198 - 35285625*z3)*pow2(Mt) - 47250*(-430 + 207*z3)*pow2(Mst1)*
1087  pow2(s2t)) + pow3(Dmst12)*(2*Mst1*s2t*(-28188929 + 23099625*z3)*pow2(
1088  Mt) + 225*Mt*(-160136 + 100785*z3)*pow2(Mst1)*pow2(s2t) + (115176444 -
1089  85920750*z3)*pow3(Mt) + 1125*(22174 - 18791*z3)*pow3(Mst1)*pow3(s2t)) +
1090  40*Dmst12*(150*Mst1*s2t*(-19856 + 17667*z3) + Mt*(-5136871 + 3912300*
1091  z3))*pow2(Mt)*pow4(Mst2) + 6000*(-31142 + 22653*z3)*pow3(Mt)*pow6(Mst2)
1092  )) + 2*Mgl*pow2(Dmglst1)*(11760*pow2(Dmsqst1)*(-3*Mt*pow2(Dmst12)*pow2(
1093  Mst2)*(4470*Mst1*Mt*s2t + 26294*pow2(Mt) + 11975*pow2(Mst1)*pow2(s2t))
1094  + pow3(Dmst12)*(25750*Mst1*s2t*pow2(Mt) + 71850*Mt*pow2(Mst1)*pow2(s2t)
1095  + 78882*pow3(Mt) + 625*pow3(Mst1)*pow3(s2t)) + 2*Dmst12*(39441*Mt +
1096  535*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 128066*pow3(Mt)*pow6(Mst2)) +
1097  11760*Dmsqst1*pow2(Msq)*(-(Mt*pow2(Dmst12)*pow2(Mst2)*(8660*Mst1*Mt*s2t
1098  + 78882*pow2(Mt) + 35925*pow2(Mst1)*pow2(s2t))) + pow3(Dmst12)*(21000*
1099  Mst1*s2t*pow2(Mt) + 71850*Mt*pow2(Mst1)*pow2(s2t) + 78882*pow3(Mt) +
1100  625*pow3(Mst1)*pow3(s2t)) + 2*Dmst12*(39441*Mt - 1840*Mst1*s2t)*pow2(
1101  Mt)*pow4(Mst2) + 160280*pow3(Mt)*pow6(Mst2)) + pow4(Msq)*(Mt*pow2(
1102  Dmst12)*pow2(Mst2)*(196*Mst1*Mt*s2t*(-353948822 + 292953375*z3) + (
1103  31554389946 - 25980577875*z3)*pow2(Mt) + 22050*(26498 + 49425*z3)*pow2(
1104  Mst1)*pow2(s2t)) + 5*pow3(Dmst12)*(-196*Mst1*s2t*(-148185343 +
1105  123161625*z3)*pow2(Mt) + 4410*Mt*(-612347 + 438045*z3)*pow2(Mst1)*pow2(
1106  s2t) + (-2800003036 + 2827317150*z3)*pow3(Mt) - 735*(-2573582 +
1107  2144805*z3)*pow3(Mst1)*pow3(s2t)) + 392*Dmst12*(260*Mst1*s2t*(-579323 +
1108  490725*z3) + Mt*(-125277461 + 96491250*z3))*pow2(Mt)*pow4(Mst2) + 3920*
1109  (-12894992 + 9523575*z3)*pow3(Mt)*pow6(Mst2))) + 4*pow3(Dmglst1)*(5880*
1110  pow2(Dmsqst1)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*(13291*Mst1*Mt*s2t +
1111  40812*pow2(Mt) + 18300*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(38922*
1112  Mst1*s2t*pow2(Mt) + 73200*Mt*pow2(Mst1)*pow2(s2t) + 81624*pow3(Mt) +
1113  625*pow3(Mst1)*pow3(s2t)) + 2*Dmst12*(40812*Mt + 7121*Mst1*s2t)*pow2(
1114  Mt)*pow4(Mst2) + 122500*pow3(Mt)*pow6(Mst2)) + 5880*Dmsqst1*pow2(Msq)*(
1115  -8*Mt*pow2(Dmst12)*pow2(Mst2)*(2729*Mst1*Mt*s2t + 10203*pow2(Mt) +
1116  4575*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(34172*Mst1*s2t*pow2(Mt) +
1117  73200*Mt*pow2(Mst1)*pow2(s2t) + 81624*pow3(Mt) + 625*pow3(Mst1)*pow3(
1118  s2t)) + 12*Dmst12*(6802*Mt + 791*Mst1*s2t)*pow2(Mt)*pow4(Mst2) +
1119  153928*pow3(Mt)*pow6(Mst2)) + pow4(Msq)*(2*Mt*pow2(Dmst12)*pow2(Mst2)*(
1120  Mst1*Mt*s2t*(-51549748862 + 42804507375*z3) + (15344105782 -
1121  12740986125*z3)*pow2(Mt) + 36750*(-109771 + 107379*z3)*pow2(Mst1)*pow2(
1122  s2t)) + pow3(Dmst12)*(2*Mst1*s2t*(137797425107 - 114549584625*z3)*pow2(
1123  Mt) - 3675*Mt*(-973342 + 1115325*z3)*pow2(Mst1)*pow2(s2t) + 12*(-
1124  1749149438 + 1631424375*z3)*pow3(Mt) - 6370*(-2386547 + 1986525*z3)*
1125  pow3(Mst1)*pow3(s2t)) + 8*Dmst12*(49*Mst1*s2t*(-244084964 + 203974875*
1126  z3) + Mt*(-5048328734 + 3923356500*z3))*pow2(Mt)*pow4(Mst2) + 2352*(-
1127  21126629 + 16012125*z3)*pow3(Mt)*pow6(Mst2))))/(1.9845e6*pow3(Mgl)*
1128  pow3(Mt)*pow4(Msq)*pow6(Mst2));
1129 
1130  return result;
1131 }
1132 
1133 /**
1134  * @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)
1135  */
1136 double H3::calc_coef_at_as2_no_sm_logs_log1() const {
1137 
1138  const double result =
1139  (4*(2*pow3(Dmglst1)*(2940*pow2(Dmsqst1)*pow2(Mt)*(-((42*Mt + 11*Mst1*s2t)
1140  *pow2(Dmst12)*pow2(Mst2)) + (42*Mt + 51*Mst1*s2t)*pow3(Dmst12) +
1141  Dmst12*(42*Mt - 29*Mst1*s2t)*pow4(Mst2) - 65*Mt*pow6(Mst2)) + 5880*
1142  Dmsqst1*pow2(Msq)*pow2(Mt)*(7*(-3*Mt + Mst1*s2t)*pow2(Dmst12)*pow2(
1143  Mst2) + (21*Mt + 13*Mst1*s2t)*pow3(Dmst12) + 3*Dmst12*(7*Mt - 9*Mst1*
1144  s2t)*pow4(Mst2) - 28*Mt*pow6(Mst2)) + pow4(Msq)*(-10*Mt*pow2(Dmst12)*
1145  pow2(Mst2)*(623998*Mst1*Mt*s2t + 46947*pow2(Mt) + 134505*pow2(Mst1)*
1146  pow2(s2t)) + pow3(Dmst12)*(17975555*Mst1*s2t*pow2(Mt) + 1956570*Mt*
1147  pow2(Mst1)*pow2(s2t) + 2032034*pow3(Mt) + 1187760*pow3(Mst1)*pow3(s2t))
1148  - 2*Dmst12*(546547*Mt + 2011058*Mst1*s2t)*pow2(Mt)*pow4(Mst2) -
1149  11351536*pow3(Mt)*pow6(Mst2))) - 98*Dmglst1*pow2(Mgl)*(600*Dmsqst1*
1150  pow2(Msq)*pow2(Mt)*(-((6*Mt + Mst1*s2t)*pow2(Dmst12)*pow2(Mst2)) + (6*
1151  Mt - 3*Mst1*s2t)*pow3(Dmst12) + Dmst12*(6*Mt + 5*Mst1*s2t)*pow4(Mst2) +
1152  25*Mt*pow6(Mst2)) + 300*pow2(Dmsqst1)*pow2(Mt)*(3*(-4*Mt + Mst1*s2t)*
1153  pow2(Dmst12)*pow2(Mst2) + (12*Mt - 11*Mst1*s2t)*pow3(Dmst12) + Dmst12*(
1154  12*Mt + 5*Mst1*s2t)*pow4(Mst2) + 35*Mt*pow6(Mst2)) + pow4(Msq)*(Mt*
1155  pow2(Dmst12)*pow2(Mst2)*(29758*Mst1*Mt*s2t + 6677*pow2(Mt) + 22350*
1156  pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-34587*Mst1*s2t*pow2(Mt) - 18045*
1157  Mt*pow2(Mst1)*pow2(s2t) + 22414*pow3(Mt) + 3325*pow3(Mst1)*pow3(s2t)) -
1158  8*Dmst12*(4471*Mt + 6875*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 50800*pow3(Mt)
1159  *pow6(Mst2))) - Mgl*pow2(Dmglst1)*(5880*pow2(Dmsqst1)*pow2(Mt)*(-3*(3*
1160  Mt + 5*Mst1*s2t)*pow2(Dmst12)*pow2(Mst2) + (9*Mt - 25*Mst1*s2t)*pow3(
1161  Dmst12) + Dmst12*(9*Mt + 55*Mst1*s2t)*pow4(Mst2) + 112*Mt*pow6(Mst2)) +
1162  5880*Dmsqst1*pow2(Msq)*pow2(Mt)*(-((9*Mt + 40*Mst1*s2t)*pow2(Dmst12)*
1163  pow2(Mst2)) + 9*Mt*pow3(Dmst12) + Dmst12*(9*Mt + 80*Mst1*s2t)*pow4(
1164  Mst2) + 145*Mt*pow6(Mst2)) + pow4(Msq)*(3*Mt*pow2(Dmst12)*pow2(Mst2)*(
1165  2334948*Mst1*Mt*s2t - 104761*pow2(Mt) + 112210*pow2(Mst1)*pow2(s2t)) +
1166  pow3(Dmst12)*(-10892014*Mst1*s2t*pow2(Mt) + 1366365*Mt*pow2(Mst1)*pow2(
1167  s2t) + 177178*pow3(Mt) - 363580*pow3(Mst1)*pow3(s2t)) + 196*Dmst12*(
1168  2303*Mt - 30942*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 13177472*pow3(Mt)*pow6(
1169  Mst2))) - 49*pow3(Mgl)*(-150*Dmsqst1*Mt*pow2(Msq)*(pow2(Dmst12)*pow2(
1170  Mst2)*(-80*Mst1*Mt*s2t - 17*pow2(Mt) + 180*pow2(Mst1)*pow2(s2t)) + 2*(
1171  20*Mst1*Mt*s2t + 187*pow2(Mt) - 90*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) -
1172  20*Dmst12*Mt*(17*Mt - 6*Mst1*s2t)*pow4(Mst2) - 1080*pow2(Mt)*pow6(Mst2)
1173  ) - 150*Mt*pow2(Dmsqst1)*(-4*pow2(Dmst12)*pow2(Mst2)*(10*Mst1*Mt*s2t +
1174  3*pow2(Mt) - 45*pow2(Mst1)*pow2(s2t)) + 9*(41*pow2(Mt) - 20*pow2(Mst1)*
1175  pow2(s2t))*pow3(Dmst12) + 5*Dmst12*Mt*(-69*Mt + 16*Mst1*s2t)*pow4(Mst2)
1176  - 1010*pow2(Mt)*pow6(Mst2)) - pow4(Msq)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*
1177  (25850*Mst1*Mt*s2t + 21033*pow2(Mt) + 75*pow2(Mst1)*pow2(s2t)) + pow3(
1178  Dmst12)*(68264*Mst1*s2t*pow2(Mt) + 5700*Mt*pow2(Mst1)*pow2(s2t) +
1179  36107*pow3(Mt) + 250*pow3(Mst1)*pow3(s2t)) + 200*Dmst12*(131*Mt + 100*
1180  Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 400*(-533 + 54*z3)*pow3(Mt)*pow6(Mst2))
1181  )))/(33075.*pow3(Mgl)*pow3(Mt)*pow4(Msq)*pow6(Mst2));
1182 
1183  return result;
1184 }
1185 
1186 /**
1187  * @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)
1188  */
1189 double H3::calc_coef_at_as2_no_sm_logs_log2() const {
1190 
1191  const double result =
1192  (8*(7*Dmglst1*pow2(Mgl)*pow4(Msq)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*(-1304*
1193  Mst1*Mt*s2t + 888*pow2(Mt) + 645*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(
1194  -1544*Mst1*s2t*pow2(Mt) + 2250*Mt*pow2(Mst1)*pow2(s2t) + 1640*pow3(Mt)
1195  - 275*pow3(Mst1)*pow3(s2t)) + 8*Dmst12*(239*Mt - 35*Mst1*s2t)*pow2(Mt)*
1196  pow4(Mst2) + 6520*pow3(Mt)*pow6(Mst2)) + pow3(Dmglst1)*pow4(Msq)*(-4*
1197  Mt*pow2(Dmst12)*pow2(Mst2)*(-6158*Mst1*Mt*s2t + 7818*pow2(Mt) - 5565*
1198  pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-109632*Mst1*s2t*pow2(Mt) - 8820*
1199  Mt*pow2(Mst1)*pow2(s2t) + 98096*pow3(Mt) - 9765*pow3(Mst1)*pow3(s2t)) +
1200  16*Dmst12*(-2222*Mt + 5257*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 36232*pow3(
1201  Mt)*pow6(Mst2)) + Mgl*pow2(Dmglst1)*pow4(Msq)*(-3*Mt*pow2(Dmst12)*pow2(
1202  Mst2)*(-7448*Mst1*Mt*s2t + 11772*pow2(Mt) + 35*pow2(Mst1)*pow2(s2t)) +
1203  pow3(Dmst12)*(-53536*Mst1*s2t*pow2(Mt) + 16905*Mt*pow2(Mst1)*pow2(s2t)
1204  + 68252*pow3(Mt) - 5285*pow3(Mst1)*pow3(s2t)) + 28*Dmst12*(85*Mt +
1205  1164*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 49588*pow3(Mt)*pow6(Mst2)) + 7*
1206  pow3(Mgl)*(Mt*pow2(Dmst12)*pow2(Mst2)*(1760*Mst1*Mt*s2t + 538*pow2(Mt)
1207  + 105*pow2(Mst1)*pow2(s2t))*pow4(Msq) - pow3(Dmst12)*(1032*Mst1*s2t*
1208  pow2(Mt) + 480*Mt*pow2(Mst1)*pow2(s2t) + 644*pow3(Mt) - 45*pow3(Mst1)*
1209  pow3(s2t))*pow4(Msq) - 40*Dmst12*(11*Mt + 61*Mst1*s2t)*pow2(Mt)*pow4(
1210  Msq)*pow4(Mst2) + 10*pow3(Mt)*(45*pow2(Dmsqst1) + 90*Dmsqst1*pow2(Msq)
1211  + 442*pow4(Msq))*pow6(Mst2))))/(315.*pow3(Mgl)*pow3(Mt)*pow4(Msq)*pow6(
1212  Mst2));
1213 
1214  return result;
1215 }
1216 
1217 /**
1218  * @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)
1219  */
1220 double H3::calc_coef_at_as2_no_sm_logs_log3() const {
1221 
1222  const double result = -298.6666666666667;
1223 
1224  return result;
1225 }
1226 
1227 } // namespace hierarchies
1228 } // namespace himalaya
himalaya::hierarchies::H3::getS2
double getS2() const
Definition: H3.cpp:181
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::ExpansionFlags_t
std::array< int, ExpansionDepth::NUMBER_OF_EXPANSIONS > ExpansionFlags_t
Definition: HierarchyTypes.hpp:17
himalaya::hierarchies::H3::calc_coef_at_as2_no_sm_logs_log1
double calc_coef_at_as2_no_sm_logs_log1() const
Definition: H3.cpp:1136
himalaya
Definition: H3.cpp:14
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::H3::getS1
double getS1() const
Definition: H3.cpp:75
himalaya::hierarchies::H3::lmMt
double lmMt
common variables
Definition: H3.hpp:73
himalaya::hierarchies::H3::calc_coef_at_as2_no_sm_logs_log0
double calc_coef_at_as2_no_sm_logs_log0() const
Definition: H3.cpp:1054
himalaya::hierarchies::H3::calc_coef_at_as2_no_sm_logs_log2
double calc_coef_at_as2_no_sm_logs_log2() const
Definition: H3.cpp:1189
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::H3::getS12
double getS12() const
Definition: H3.cpp:777
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::H3::H3
H3(const ExpansionFlags_t &expansionDepth, double Al4p, double beta, double Dmglst1, double Dmst12, double Dmsqst1, double lmMt, 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: H3.cpp:38
himalaya::hierarchies::H3::xDmsqst1
int xDmsqst1
MDR and truncation flags.
Definition: H3.hpp:74
himalaya::hierarchies::H3::threeLoopFlag
int threeLoopFlag
loop flags
Definition: H3.hpp:75
himalaya::hierarchies::H3::calc_coef_at_as2_no_sm_logs_log3
double calc_coef_at_as2_no_sm_logs_log3() const
Definition: H3.cpp:1220
Generated by   doxygen 1.8.13