Himalaya
H3q22g.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/H3q22g.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 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 Mt a double top/bottom quark mass
28  * @param Mst1 a double stop 1 mass
29  * @param Mst2 a double stop 2 mass
30  * @param Msq a double average squark mass w/o the stop quarks
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 H3q22g::H3q22g(const ExpansionFlags_t& expansionDepth, double Al4p, double beta,
39  double Dmglst1, double Dmst12, double Dmsqst1, double lmMt, double lmMst1,
40  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 -> Mt = Mt;
51  this -> Mst1 = Mst1;
52  this -> Mst2 = Mst2;
53  this -> Msq = Msq;
54  this -> MuSUSY = MuSUSY;
55  this -> s2t = s2t;
56  this -> lmMt = lmMt;
57  this -> oneLoopFlag = oneLoopFlag;
58  this -> twoLoopFlag = twoLoopFlag;
59  this -> threeLoopFlag = threeLoopFlag;
60  this -> Al4p = Al4p;
61  // mdr flags, indicates if one wants to shift the dr stop mass to the mdr stop mass
62  shiftst1 = mdrFlag;
63  shiftst2 = mdrFlag;
64  shiftst3 = mdrFlag;
65  // expansion flags
66  xDmglst1 = expansionDepth.at(ExpansionDepth::Dmglst1);
67  xDmst12 = expansionDepth.at(ExpansionDepth::Dmglst1);
68  xDmsqst1 = expansionDepth.at(ExpansionDepth::Dmsqst1);
69 }
70 
71 /**
72  * @return The diagonal (1, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3q22g'
73  */
74 double H3q22g::getS1() const {
75  return -(pow2(Mt)*pow2(MuSUSY)*(-4*Mst1*xDmst12*pow3(Dmst12)*(12000*Dmsqst1*s2t*
76  threeLoopFlag*pow2(Al4p)*(20*Dmglst1*Mt*(Dmsqst1*xDmsqst1 + pow2(Msq))
77  + 3*Mst1*(Dmsqst1*(28*Mt + Mst1*s2t*(7 + 6*shiftst1 - 24*lmMst1*(-1 +
78  shiftst2) + 30*shiftst2))*xDmsqst1 + (28*Mt + Mst1*s2t*(7 - 24*lmMst1*(
79  -1 + shiftst2) + 36*shiftst2))*pow2(Msq)))*pow2(Mst1) + pow4(Msq)*(
80  Al4p*(8*Dmglst1*Mt*(4*Dmglst1*(720*Dmglst1*s2t*twoLoopFlag*xDmglst1 +
81  Al4p*threeLoopFlag*(Dmglst1*s2t*xDmglst1*(14217821 + 161940*lmMst1 -
82  28800*pow2(lmMst1)) + Mt*(6233611 + 58800*lmMst1 - 4560*lmMt - 14400*
83  pow2(lmMst1)))) + 5*Mst1*(-2160*Dmglst1*s2t*twoLoopFlag + Al4p*
84  threeLoopFlag*(6*Dmglst1*s2t*(954181 + 11256*lmMst1 - 11520*pow2(
85  lmMst1)) + 5*Mt*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))))) -
86  15*pow2(Mst1)*(960*Dmglst1*s2t*((-5 + 6*lmMst1)*Mt + Dmglst1*(-11 + 6*
87  lmMst1)*s2t)*twoLoopFlag + Al4p*threeLoopFlag*(16*Dmglst1*s2t*(Mt*(
88  1282471 - 7264*lmMst1 - 18120*pow2(lmMst1)) + Dmglst1*s2t*(655743 +
89  5288*lmMst1 - 11820*pow2(lmMst1))) + (1763661 - 47104*lmMst1 + 5120*
90  lmMt - 24576*pow2(lmMst1))*pow2(Mt))) - 600*s2t*(24*((5 + 6*lmMst1)*Mt
91  + 4*Dmglst1*(-1 + 3*lmMst1)*s2t)*twoLoopFlag + Al4p*threeLoopFlag*(
92  Dmglst1*s2t*(84209 + 1264*lmMst1 - 240*pow2(lmMst1)) - 2*Mt*(-36863 +
93  80*lmMst1 + 552*pow2(lmMst1))))*pow3(Mst1)) + 15*(270*oneLoopFlag +
94  Al4p*(240*(twoLoopFlag + 6*lmMst1*twoLoopFlag) + Al4p*threeLoopFlag*(
95  350605 + 4320*shiftst1 + 2880*shiftst2 + 8352*shiftst3 - 96*lmMst1*(-
96  115 + 90*shiftst1 + 60*shiftst2 + 54*shiftst3) - 2160*pow2(lmMst1))))*
97  pow2(s2t)*pow4(Mst1))) - threeLoopFlag*pow2(Al4p)*(480*pow2(s2t)*pow3(
98  Mst1)*((Dmst12*pow2(Msq)*pow2(Mst2)*(2*Dmst12*pow2(Dmglst1)*(655743 +
99  5288*lmMst1 - 11820*pow2(lmMst1))*pow2(Msq) + 5*Dmglst1*Dmst12*Mst1*(
100  84209 + 1264*lmMst1 - 240*pow2(lmMst1))*pow2(Msq) + 10*pow2(Mst1)*(30*
101  Dmsqst1*(Dmst12*(-7 - 18*shiftst1 - 18*shiftst2 + 12*lmMst1*(-2 +
102  shiftst1 + shiftst2)) - 12*(-3 + 2*lmMst1)*(shiftst1 - shiftst2)*pow2(
103  Mst2)) + pow2(Msq)*(Dmst12*(1429 - 180*shiftst1 - 180*shiftst2 +
104  lmMst1*(-454 + 360*shiftst1 + 360*shiftst2) + 24*pow2(lmMst1)) - 360*(-
105  1 + 2*lmMst1)*(shiftst1 - shiftst2)*pow2(Mst2)))))/2. + pow2(Mst1)*(90*
106  Dmst12*shiftst3*pow2(Mst2)*(Dmst12*(-7 + 6*lmMst1) + 2*(1 - 2*lmMst1)*
107  pow2(Mst2))*pow4(Msq) + 36*z2*(-2*xDmst12*pow3(Dmst12)*(100*Dmsqst1*
108  shiftst2*(Dmsqst1*xDmsqst1 + pow2(Msq)) + (15*shiftst1 + 10*shiftst2 +
109  9*shiftst3)*pow4(Msq)) + Dmst12*pow2(Mst2)*(5*shiftst3*(3*Dmst12 - 2*
110  pow2(Mst2))*pow4(Msq) - 50*(-(Dmst12*(shiftst1 + shiftst2)) + 2*(
111  shiftst1 - shiftst2)*pow2(Mst2))*(xDmsqst1*pow2(Dmsqst1) + Dmsqst1*
112  pow2(Msq) + pow4(Msq)))))) + 2*Mst1*pow2(Mst2)*(240*Dmst12*Mst1*s2t*(-
113  50*Mst1*xDmsqst1*pow2(Dmsqst1)*(40*Dmglst1*Mt*(Dmst12 - pow2(Mst2)) +
114  3*Mst1*(Dmst12*(56*Mt + Mst1*s2t*(7 + 30*shiftst1 + 6*shiftst2 - 12*
115  lmMst1*(-2 + shiftst1 + shiftst2))) + 8*(-7*Mt + 3*(-2 + lmMst1)*Mst1*
116  s2t*(shiftst1 - shiftst2))*pow2(Mst2))) - Mt*pow2(Msq)*(2*pow2(Dmglst1)
117  *(954181 + 11256*lmMst1 - 11520*pow2(lmMst1))*pow2(Msq)*(Dmst12 - pow2(
118  Mst2)) + 5*pow2(Mst1)*(1680*Dmsqst1*(Dmst12 - pow2(Mst2)) + pow2(Msq)*(
119  3*Dmst12*(-10473 + 40*lmMst1 + 256*pow2(lmMst1)) - 8*(1361 + 10*lmMst1
120  + 54*pow2(lmMst1))*pow2(Mst2))) + Dmglst1*Mst1*(2000*Dmsqst1*(Dmst12 -
121  pow2(Mst2)) + pow2(Msq)*(Dmst12*(-949861 + 1944*lmMst1 + 11520*pow2(
122  lmMst1)) + 20*(-33261 + 532*lmMst1 + 660*pow2(lmMst1))*pow2(Mst2))))) +
123  pow2(Mt)*pow4(Msq)*(-400*Dmglst1*Mst1*(Dmst12*(403559 + 384*(lmMst1 +
124  lmMt) - 4608*pow2(lmMst1))*(Dmst12 - pow2(Mst2)) - 24*(9631 + 16*lmMst1
125  + 48*lmMt - 192*pow2(lmMst1))*pow4(Mst2)) + 15*pow2(Mst1)*(pow2(Dmst12)
126  *(852541 + 9216*lmMst1 - 10240*lmMt + 6144*pow2(lmMst1)) + 160*Dmst12*(
127  11389 - 704*lmMst1 + 192*lmMt - 384*pow2(lmMst1))*pow2(Mst2) + 1920*(
128  349 - 56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow4(Mst2)) - 32*pow2(
129  Dmglst1)*(-2*Dmst12*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(
130  lmMst1))*(Dmst12 - pow2(Mst2)) + 25*(-230009 - 4704*lmMst1 + 480*lmMt +
131  1152*pow2(lmMst1))*pow4(Mst2)))) + 225*z3*(-2*Mst1*xDmst12*pow3(Dmst12)
132  *(Mt*(-8*Dmglst1*Mst1*(-73917*Mt + 285974*Mst1*s2t) + 48*(30678*Mt +
133  35135*Mst1*s2t)*pow2(Dmglst1) - (197889*Mt + 324752*Mst1*s2t)*pow2(
134  Mst1) + 3371456*s2t*xDmglst1*pow3(Dmglst1))*pow4(Msq) - pow2(Mst1)*
135  pow2(s2t)*(-10080*Dmsqst1*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mst1) + (
136  403880*Dmglst1*Mst1 + 1197040*pow2(Dmglst1) - 37669*pow2(Mst1))*pow4(
137  Msq))) + pow2(Mst2)*(-8*pow2(Dmst12)*pow2(s2t)*pow3(Mst1)*(-1260*
138  Dmsqst1*(Dmsqst1*xDmsqst1 + pow2(Msq))*pow2(Mst1) + (50485*Dmglst1*Mst1
139  + 149630*pow2(Dmglst1) + 2574*pow2(Mst1))*pow4(Msq)) + Mt*pow4(Msq)*(
140  16*Dmst12*s2t*pow2(Mst1)*(210810*pow2(Dmglst1)*(Dmst12 - pow2(Mst2)) -
141  pow2(Mst1)*(17285*Dmst12 + 6024*pow2(Mst2)) - Dmglst1*Mst1*(105405*
142  Dmst12 + 75164*pow2(Mst2))) - 1408*xDmglst1*pow3(Dmglst1)*(-4789*
143  Dmst12*Mst1*s2t*(Dmst12 - pow2(Mst2)) + 1503*Mt*pow4(Mst2)) + 3*Mst1*
144  Mt*(96*pow2(Dmglst1)*(10226*pow2(Dmst12) - 10226*Dmst12*pow2(Mst2) -
145  4715*pow4(Mst2)) + 48*Dmglst1*Mst1*(8213*pow2(Dmst12) - 8213*Dmst12*
146  pow2(Mst2) - 4664*pow4(Mst2)) - pow2(Mst1)*(31963*pow2(Dmst12) + 68000*
147  Dmst12*pow2(Mst2) + 24064*pow4(Mst2))))))) + 8*pow2(Mst2)*pow4(Msq)*(
148  1800*Al4p*Dmst12*s2t*twoLoopFlag*pow2(Mst1)*(Dmst12*Mst1*s2t*(8*
149  Dmglst1*(1 - 3*lmMst1)*Mst1 + (22 - 12*lmMst1)*pow2(Dmglst1) + 3*pow2(
150  Mst1)) + 4*Mt*(-6*pow2(Dmglst1)*(Dmst12 - pow2(Mst2)) + Dmglst1*(5 - 6*
151  lmMst1)*Mst1*pow2(Mst2) - pow2(Mst1)*(Dmst12 + 3*(1 + 2*lmMst1)*pow2(
152  Mst2)))) - 16*Al4p*Mt*xDmglst1*pow3(Dmglst1)*(-720*Dmst12*Mst1*s2t*
153  twoLoopFlag*(Dmst12 - pow2(Mst2)) + Al4p*threeLoopFlag*(Dmst12*Mst1*
154  s2t*(-14217821 - 161940*lmMst1 + 28800*pow2(lmMst1))*(Dmst12 - pow2(
155  Mst2)) + 2*(2219399 + 9600*lmMst1 + 960*lmMt)*Mt*pow4(Mst2))) + 2025*
156  oneLoopFlag*pow2(Dmst12)*pow2(s2t)*pow5(Mst1))))/(777600.*pow4(Msq)*
157  pow5(Mst1)*pow6(Mst2));
158 }
159 
160 /**
161  * @return The diagonal (2, 2) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3q22g'
162  */
163 double H3q22g::getS2() const {
164  return (pow2(Mst2)*(Mst1*threeLoopFlag*pow2(Al4p)*(23520*Mt*xDmsqst1*pow2(
165  Dmsqst1)*(64*xDmglst1*pow2(Mt)*pow2(Sbeta)*pow3(Dmglst1)*(Dmst12*(1331
166  - 420*lmMst1 + 420*lmMt)*s2t*Tbeta*pow2(Mst1)*(Dmst12 - pow2(Mst2)) -
167  Mt*((1541 - 420*lmMst1 + 420*lmMt)*MuSUSY + (2579 + 120*lmMst1 - 120*
168  lmMt)*Mst1*Tbeta)*pow4(Mst2)) - 5*pow3(Mst1)*(50*Dmst12*Mt*pow2(Mst2)*(
169  2*(8*(65 - 6*lmMst1 + 6*lmMt)*MuSUSY + 9*Mst1*(49 + 2*lmMt + lmMst1*(46
170  - 48*shiftst2) + 96*shiftst2)*Tbeta)*pow2(Mt)*pow2(Sbeta) - Mt*s2t*(-
171  168*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(MuSUSY*(687 - 864*
172  lmMst1*(-1 + shiftst2) + 1728*shiftst2) + 16*(-85 + 12*lmMst1 - 12*
173  lmMt)*Mst1*Tbeta)*pow2(Sbeta)) - 72*(-2 + lmMst1)*Mst1*(shiftst1 -
174  shiftst2)*Tbeta*pow2(s2t)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*pow2(
175  Mst1)*pow2(Sbeta))) - pow2(Dmst12)*(-200*s2t*pow2(Mt)*(-42*Tbeta*pow2(
176  MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(3*MuSUSY*(11 + 36*shiftst2) + 4*(-23
177  + 6*lmMst1 - 6*lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + 75*Mst1*Mt*pow2(s2t)*(-
178  2*(7 + 30*shiftst1 + 6*shiftst2 - 12*lmMst1*(-2 + shiftst1 + shiftst2))
179  *Tbeta*pow2(MuSUSY)*(1 - pow2(Sbeta)) + Mst1*(168*MuSUSY + Mst1*(229 +
180  864*shiftst1 - 288*shiftst2 + 144*lmMst1*(2 - 3*shiftst1 + shiftst2))*
181  Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(8*(100*(65 - 6*lmMst1 + 6*lmMt)*
182  MuSUSY + 3*Mst1*(419 - 60*lmMst1 + 60*lmMt + 900*shiftst2)*Tbeta)*pow3(
183  Mt) - 5400*(-2 + lmMst1)*MuSUSY*(shiftst1 - shiftst2)*pow3(Mst1)*pow3(
184  s2t))) + 200*(MuSUSY*((364 - 48*lmMst1 + 48*lmMt)*Mt - 216*(-2 +
185  lmMst1)*Mst1*s2t*(shiftst1 - shiftst2)) - Mst1*Mt*Tbeta*(205 + 174*lmMt
186  - 432*(shiftst1 + shiftst2) - 6*lmMst1*(101 + 18*lmMt - 36*(shiftst1 +
187  shiftst2)) + 54*(pow2(lmMst1) + pow2(lmMt))))*pow2(Mt)*pow2(Sbeta)*
188  pow4(Mst2)) + 8*Dmglst1*Mst1*Mt*(3*Dmglst1*pow2(Sbeta)*(Dmst12*Mst1*(2*
189  Mt*(-225*MuSUSY*s2t + (457 + 510*lmMst1 - 510*lmMt)*Mt*Tbeta + 10*(397
190  - 30*lmMst1 + 30*lmMt)*Mst1*s2t*Tbeta)*pow2(Mst2) - Dmst12*(-10*Mt*s2t*
191  (45*MuSUSY + 2*(-397 + 30*lmMst1 - 30*lmMt)*Mst1*Tbeta) + 2*(457 + 510*
192  lmMst1 - 510*lmMt)*Tbeta*pow2(Mt) + 225*Tbeta*pow2(Mst1)*pow2(s2t))) +
193  2*(20*(206 - 15*lmMst1 + 15*lmMt)*MuSUSY + 9*(279 + 70*lmMst1 - 70*
194  lmMt)*Mst1*Tbeta)*pow2(Mt)*pow4(Mst2)) - 5*Mst1*(2*Dmst12*Mt*pow2(Mst2)
195  *(2*Mt*((557 + 120*lmMst1 - 120*lmMt)*MuSUSY + 9*(-423 + 20*lmMst1 -
196  20*lmMt)*Mst1*Tbeta)*pow2(Sbeta) + 25*s2t*(5*Tbeta*pow2(MuSUSY)*(-1 +
197  pow2(Sbeta)) + Mst1*(282*MuSUSY + (91 + 6*lmMst1 - 6*lmMt)*Mst1*Tbeta)*
198  pow2(Sbeta))) - pow2(Dmst12)*((4*((557 + 120*lmMst1 - 120*lmMt)*MuSUSY
199  + 9*(-423 + 20*lmMst1 - 20*lmMt)*Mst1*Tbeta)*pow2(Mt) + 75*(5*MuSUSY -
200  94*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t))*pow2(Sbeta) + 2*Mt*s2t*(125*Tbeta*
201  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*(2350*MuSUSY + (347 - 30*
202  lmMst1 + 30*lmMt)*Mst1*Tbeta)*pow2(Sbeta))) + 100*((44 + 3*lmMst1 - 3*
203  lmMt)*MuSUSY + (-226 + 21*lmMst1 - 21*lmMt)*Mst1*Tbeta)*pow2(Mt)*pow2(
204  Sbeta)*pow4(Mst2)))) + 2*Mt*pow2(Msq)*(784*Mst1*MuSUSY*pow2(Sbeta)*(
205  450*Dmst12*Mst1*s2t*(5*pow2(Mst1)*(5*Dmsqst1*(161*Dmst12 + 24*(1 + 12*
206  lmMst1)*pow2(Mst2)) + pow2(Msq)*(Dmst12*(2071 + 296*lmMst1 - 96*lmMt -
207  600*pow2(lmMst1)) + 8*(631 - lmMst1 + 36*lmMt + 21*pow2(lmMst1))*pow2(
208  Mst2))) + 4*Dmglst1*Mst1*(4700*Dmsqst1*(Dmst12 - pow2(Mst2)) - pow2(
209  Msq)*(Dmst12*(17539 + 574*lmMst1 + 160*lmMt - 1920*pow2(lmMst1)) + 5*(-
210  4539 + 596*lmMst1 - 48*lmMt + 516*pow2(lmMst1))*pow2(Mst2))) + pow2(
211  Dmglst1)*(360*Dmsqst1*(Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(-
212  515917 - 6972*lmMst1 + 192*lmMt + 11520*pow2(lmMst1)) + 2*(-32101 +
213  5044*lmMst1 - 400*lmMt + 5100*pow2(lmMst1))*pow2(Mst2))))*pow2(Mt) +
214  225*Mt*pow2(Dmst12)*pow2(Mst1)*(pow2(Dmglst1)*(954181 + 11256*lmMst1 -
215  11520*pow2(lmMst1))*pow2(Msq) + 10*Dmglst1*Mst1*(100*Dmsqst1 + (33261 -
216  532*lmMst1 - 660*pow2(lmMst1))*pow2(Msq)) + 20*(210*Dmsqst1 + (1361 +
217  10*lmMst1 + 54*pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*pow2(s2t) + 20250*
218  s2t*pow3(Mst1)*(shiftst3*pow2(Msq)*(-8*Dmst12*(-3 + 2*lmMst1)*pow2(
219  Mst2)*pow2(Mt) + pow2(Dmst12)*(16*(-2 + lmMst1)*pow2(Mt) + (1 - 2*
220  lmMst1)*pow2(Mst1)*pow2(s2t)) + 8*(-1 + 2*lmMst1)*pow2(Mt)*pow4(Mst2))
221  + 10*(Dmsqst1*(-3 + 2*lmMst1) + (-1 + 2*lmMst1)*pow2(Msq))*(-8*Dmst12*
222  shiftst2*pow2(Mst2)*pow2(Mt) + (-shiftst1 + shiftst2)*pow2(Dmst12)*
223  pow2(Mst1)*pow2(s2t) + 8*(shiftst1 - shiftst2)*pow2(Mt)*pow4(Mst2))) +
224  pow3(Mt)*(4*pow2(Dmglst1)*(3*Dmst12*(84334067 + 120*lmMst1*(2843 - 120*
225  lmMt) + 202200*lmMt - 828000*pow2(lmMst1) - 21600*pow2(lmMt))*pow2(Msq)
226  *(Dmst12 - pow2(Mst2)) + 40*(90*Dmsqst1*(206 - 15*lmMst1 + 15*lmMt) + (
227  -3030652 - 19305*lmMt + 6*lmMst1*(1183 + 645*lmMt) + 55530*pow2(lmMst1)
228  + 5400*pow2(lmMt))*pow2(Msq))*pow4(Mst2)) + 2*Dmglst1*Mst1*(1200*
229  Dmsqst1*(Dmst12*(557 + 120*lmMst1 - 120*lmMt)*(Dmst12 - pow2(Mst2)) -
230  50*(26 + 3*lmMst1 - 3*lmMt)*pow4(Mst2)) - pow2(Msq)*(-3*pow2(Dmst12)*(-
231  39474953 + lmMst1*(52780 - 87600*lmMt) - 3300*lmMt + 319200*pow2(
232  lmMst1) + 14400*pow2(lmMt)) + 40*Dmst12*(3746977 + 4005*lmMt - 18*
233  lmMst1*(2711 + 1215*lmMt) - 52380*pow2(lmMst1))*pow2(Mst2) + 6000*(9961
234  + 66*lmMt - 10*lmMst1*(67 + 24*lmMt) + 42*pow2(lmMst1) + 72*pow2(lmMt))
235  *pow4(Mst2))) + 15*pow2(Mst1)*(4000*Dmsqst1*(Dmst12*(65 - 6*lmMst1 + 6*
236  lmMt)*(Dmst12 - pow2(Mst2)) + 6*(-20 + 3*lmMst1 - 3*lmMt)*pow4(Mst2)) -
237  pow2(Msq)*(-(pow2(Dmst12)*(-2284899 + 49840*lmMt - 32*lmMst1*(-1793 +
238  555*lmMt) + 87360*pow2(lmMst1) + 28800*pow2(lmMt))) + 200*Dmst12*(11697
239  + 448*lmMst1 + 330*lmMt - 372*lmMst1*lmMt + 408*pow2(lmMst1) + 288*
240  pow2(lmMt))*pow2(Mst2) + 1600*(434 - 83*lmMst1 + 174*lmMt - 66*lmMst1*
241  lmMt + 183*pow2(lmMst1) + 108*pow2(lmMt))*pow4(Mst2))))) + Mt*Tbeta*(
242  245*(pow2(MuSUSY)*(16*pow2(Dmglst1)*pow2(Msq)*(2*Dmst12*Mt*(15*Mst1*
243  s2t*(954181 + 11256*lmMst1 - 11520*pow2(lmMst1)) - 2*Mt*(-6233611 -
244  58800*lmMst1 + 4560*lmMt + 14400*pow2(lmMst1)))*pow2(Mst2) + pow2(
245  Dmst12)*(30*Mst1*Mt*s2t*(-954181 - 11256*lmMst1 + 11520*pow2(lmMst1)) +
246  4*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(lmMst1))*pow2(Mt) +
247  15*(655743 + 5288*lmMst1 - 11820*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)) -
248  50*(-230009 - 4704*lmMst1 + 480*lmMt + 1152*pow2(lmMst1))*pow2(Mt)*
249  pow4(Mst2)) + 40*Dmglst1*Mst1*(10*Dmst12*Mt*(1200*Dmsqst1*Mst1*s2t + (
250  Mt*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1)) - 12*Mst1*s2t*(-
251  33261 + 532*lmMst1 + 660*pow2(lmMst1)))*pow2(Msq))*pow2(Mst2) - pow2(
252  Dmst12)*(12000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(6*Mst1*Mt*s2t*(-949861
253  + 1944*lmMst1 + 11520*pow2(lmMst1)) + 10*(403559 + 384*(lmMst1 + lmMt)
254  - 4608*pow2(lmMst1))*pow2(Mt) + 15*(-84209 - 1264*lmMst1 + 240*pow2(
255  lmMst1))*pow2(Mst1)*pow2(s2t))) + 240*(9631 + 16*lmMst1 + 48*lmMt -
256  192*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow4(Mst2)) + 15*pow2(Mst1)*(160*
257  Dmst12*Mt*(840*Dmsqst1*Mst1*s2t + (Mt*(11389 - 704*lmMst1 + 192*lmMt -
258  384*pow2(lmMst1)) + 4*Mst1*s2t*(1361 + 10*lmMst1 + 54*pow2(lmMst1)))*
259  pow2(Msq))*pow2(Mst2) + pow2(Dmst12)*(-2400*Dmsqst1*Mst1*s2t*(56*Mt + (
260  7 + 24*lmMst1)*Mst1*s2t) + pow2(Msq)*(-240*Mst1*Mt*s2t*(-10473 + 40*
261  lmMst1 + 256*pow2(lmMst1)) + (852541 + 9216*lmMst1 - 10240*lmMt + 6144*
262  pow2(lmMst1))*pow2(Mt) + 80*(1429 - 454*lmMst1 + 24*pow2(lmMst1))*pow2(
263  Mst1)*pow2(s2t))) + 1920*(349 - 56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*
264  pow2(Msq)*pow2(Mt)*pow4(Mst2))) + 21600*pow4(Mst1)*(-(shiftst3*pow2(
265  Msq)*(-2*Dmst12*pow2(Mst2)*((-1 + 2*lmMst1)*pow2(MuSUSY)*pow2(s2t)*(-1
266  + pow2(Sbeta)) + 12*((3 - 2*lmMst1)*pow2(Mt) + (-1 + 2*lmMst1)*pow2(
267  Mst1)*pow2(s2t))*pow2(Sbeta)) + pow2(Dmst12)*((-7 + 6*lmMst1)*pow2(
268  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 6*(-8*(-2 + lmMst1)*pow2(Mt) + (
269  -15 + 14*lmMst1)*pow2(Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(1 - 2*lmMst1)
270  *pow2(Mt)*pow2(Sbeta)*pow4(Mst2))) - 10*(Dmsqst1*(3 - 2*lmMst1) + (1 -
271  2*lmMst1)*pow2(Msq))*(pow2(Dmst12)*pow2(s2t)*(-((shiftst1 + shiftst2)*
272  pow2(MuSUSY)*(-1 + pow2(Sbeta))) + 6*(-3*shiftst1 + shiftst2)*pow2(
273  Mst1)*pow2(Sbeta)) + 2*Dmst12*pow2(Mst2)*((shiftst1 - shiftst2)*pow2(
274  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 12*(shiftst2*pow2(Mt) + (
275  shiftst1 - shiftst2)*pow2(Mst1)*pow2(s2t))*pow2(Sbeta)) + 24*(shiftst1
276  + shiftst2)*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)))) + pow2(Sbeta)*(784*
277  Dmst12*Mst1*Mt*s2t*(375*pow2(Mst1)*(-80*Dmsqst1*(Dmst12*((-98 + 24*
278  lmMst1 - 24*lmMt)*pow2(Mst1) - 21*pow2(MuSUSY)) + 3*pow2(Mst2)*((74 -
279  12*lmMst1 + 12*lmMt)*pow2(Mst1) + 7*pow2(MuSUSY))) + pow2(Msq)*(-8*
280  pow2(Mst2)*(8*(347 - 50*lmMst1 + 66*lmMt - 66*lmMst1*lmMt + 183*pow2(
281  lmMst1) + 108*pow2(lmMt))*pow2(Mst1) + (1361 + 10*lmMst1 + 54*pow2(
282  lmMst1))*pow2(MuSUSY)) + Dmst12*(16*(-4378 - 517*lmMst1 + 243*lmMt -
283  78*lmMst1*lmMt + 528*pow2(lmMst1) + 288*pow2(lmMt))*pow2(Mst1) + 3*(-
284  10473 + 40*lmMst1 + 256*pow2(lmMst1))*pow2(MuSUSY)))) - 2*pow2(Dmglst1)
285  *(3600*Dmsqst1*(397 - 30*lmMst1 + 30*lmMt)*pow2(Mst1)*(Dmst12 - pow2(
286  Mst2)) + pow2(Msq)*(5*pow2(Mst2)*(-8*(-6041999 - 49410*lmMt + 6*lmMst1*
287  (1721 + 1290*lmMt) + 111060*pow2(lmMst1) + 10800*pow2(lmMt))*pow2(Mst1)
288  + 15*(954181 + 11256*lmMst1 - 11520*pow2(lmMst1))*pow2(MuSUSY)) +
289  Dmst12*((263717842 - 633600*lmMt + 2400*lmMst1*(1043 + 93*lmMt) -
290  525600*pow2(lmMst1) + 302400*pow2(lmMt))*pow2(Mst1) + 75*(-954181 -
291  11256*lmMst1 + 11520*pow2(lmMst1))*pow2(MuSUSY)))) - 5*Dmglst1*Mst1*(
292  240*Dmsqst1*(25*pow2(Mst2)*(2*(55 + 6*lmMst1 - 6*lmMt)*pow2(Mst1) + 5*
293  pow2(MuSUSY)) - Dmst12*(4*(379 + 15*lmMst1 - 15*lmMt)*pow2(Mst1) + 125*
294  pow2(MuSUSY))) + pow2(Msq)*(300*pow2(Mst2)*(16*(4964 - 3*lmMt - 5*
295  lmMst1*(55 + 24*lmMt) + 21*pow2(lmMst1) + 36*pow2(lmMt))*pow2(Mst1) + (
296  33261 - 532*lmMst1 - 660*pow2(lmMst1))*pow2(MuSUSY)) + Dmst12*(8*(
297  4511549 + 9810*lmMt + 6*lmMst1*(14879 + 4710*lmMt) - 117360*pow2(
298  lmMst1) - 21600*pow2(lmMt))*pow2(Mst1) - 15*(-949861 + 1944*lmMst1 +
299  11520*pow2(lmMst1))*pow2(MuSUSY))))) - 29400*pow2(Dmst12)*pow2(Mst1)*
300  pow2(s2t)*(2*pow2(Dmglst1)*(1080*Dmsqst1*pow2(Mst1) + pow2(Msq)*(6*(
301  31901 - 5044*lmMst1 + 400*lmMt - 5100*pow2(lmMst1))*pow2(Mst1) + (
302  655743 + 5288*lmMst1 - 11820*pow2(lmMst1))*pow2(MuSUSY))) - 10*pow2(
303  Mst1)*(30*Dmsqst1*(12*(1 + 12*lmMst1)*pow2(Mst1) + (7 + 24*lmMst1)*
304  pow2(MuSUSY)) + pow2(Msq)*(24*(613 - lmMst1 + 36*lmMt + 21*pow2(lmMst1)
305  )*pow2(Mst1) + (-1429 + 454*lmMst1 - 24*pow2(lmMst1))*pow2(MuSUSY))) +
306  5*Dmglst1*(Mst1*pow2(Msq)*(24*(-4515 + 596*lmMst1 - 48*lmMt + 516*pow2(
307  lmMst1))*pow2(Mst1) + (84209 + 1264*lmMst1 - 240*pow2(lmMst1))*pow2(
308  MuSUSY)) + 22560*Dmsqst1*pow3(Mst1))) - pow2(Mt)*(9408000*(10*Dmglst1*
309  Dmsqst1*(-146 + 15*lmMst1 - 15*lmMt) - 45*Dmsqst1*Mst1*(14 + 6*lmMt -
310  6*lmMst1*(3 + lmMt) + 3*pow2(lmMst1) + 3*pow2(lmMt)) - Dmglst1*(-15571
311  + 852*lmMt + lmMst1*(-508 + 48*lmMt) + 978*pow2(lmMst1) + 612*pow2(
312  lmMt))*pow2(Msq) - Mst1*pow2(Msq)*(-623 + 555*lmMt + (663 + 387*lmMt)*
313  pow2(lmMst1) + 486*pow2(lmMt) - lmMst1*(1066 + 843*lmMt + 756*pow2(
314  lmMt)) - 252*pow3(lmMst1) + 621*pow3(lmMt)))*pow3(Mst1)*pow4(Mst2) +
315  16*pow2(Mst1)*(pow2(Dmst12)*(pow2(Dmglst1)*(35280*Dmsqst1*(457 + 510*
316  lmMst1 - 510*lmMt) + (-13564884271 - 96403020*lmMt - 60*lmMst1*(968629
317  + 101640*lmMt) + 288338400*pow2(lmMst1))*pow2(Msq)) + 98*Dmglst1*Mst1*(
318  10800*Dmsqst1*(423 - 20*lmMst1 + 20*lmMt) + (-22574599 - 21060*lmMt +
319  60*lmMst1*(6677 + 8880*lmMt) + 1598400*pow2(lmMst1) + 172800*pow2(lmMt)
320  )*pow2(Msq)) + 980*(15*Dmsqst1*(4561 + 510*lmMst1 - 510*lmMt) + (93973
321  - 61305*lmMt + 54*lmMst1*(2337 + 1255*lmMt) - 48420*pow2(lmMst1) -
322  58050*pow2(lmMt))*pow2(Msq))*pow2(Mst1)) - 196*Dmst12*(10*Dmglst1*Mst1*
323  (540*Dmsqst1*(423 - 20*lmMst1 + 20*lmMt) + (-5136871 - 29790*lmMt + 12*
324  lmMst1*(8942 + 4305*lmMt) + 86040*pow2(lmMst1) + 21600*pow2(lmMt))*
325  pow2(Msq)) + pow2(Dmglst1)*(180*Dmsqst1*(457 + 510*lmMst1 - 510*lmMt) -
326  (73908751 + 1029690*lmMt + 540*lmMst1*(2243 + 295*lmMt) + 707400*pow2(
327  lmMst1) + 64800*pow2(lmMt))*pow2(Msq)) - 750*(5*Dmsqst1*(79 + 204*
328  lmMst1 + 12*lmMt) + (-4*lmMst1*(131 + 264*lmMt) + 264*pow2(lmMst1) +
329  21*(783 + 14*lmMt + 30*pow2(lmMt)))*pow2(Msq))*pow2(Mst1))*pow2(Mst2) -
330  7840*pow2(Dmglst1)*(315*Dmsqst1*(34 + 15*lmMst1 - 15*lmMt) + (-2055923
331  - 61695*lmMt + 54*lmMst1*(-1162 + 445*lmMt) + 6345*pow2(lmMst1) -
332  12150*pow2(lmMt))*pow2(Msq))*pow4(Mst2)) + 245*pow2(Msq)*pow2(MuSUSY)*(
333  -400*Dmglst1*Mst1*(Dmst12*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(
334  lmMst1))*(Dmst12 - pow2(Mst2)) - 24*(9631 + 16*lmMst1 + 48*lmMt - 192*
335  pow2(lmMst1))*pow4(Mst2)) + 15*pow2(Mst1)*(pow2(Dmst12)*(852541 + 9216*
336  lmMst1 - 10240*lmMt + 6144*pow2(lmMst1)) + 160*Dmst12*(11389 - 704*
337  lmMst1 + 192*lmMt - 384*pow2(lmMst1))*pow2(Mst2) + 1920*(349 - 56*
338  lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow4(Mst2)) - 32*pow2(Dmglst1)*(-2*
339  Dmst12*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(lmMst1))*(
340  Dmst12 - pow2(Mst2)) + 25*(-230009 - 4704*lmMst1 + 480*lmMt + 1152*
341  pow2(lmMst1))*pow4(Mst2)))))))) + pow2(Mt)*(Al4p*(-64*xDmglst1*pow2(
342  Msq)*pow3(Dmglst1)*(8820*Mst1*twoLoopFlag*pow2(Msq)*(Dmst12*Mt*pow2(
343  Mst2)*(40*s2t*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mt*((181 - 90*
344  lmMst1 + 90*lmMt)*MuSUSY + 800*Mst1*Tbeta)*pow2(Sbeta) - 6*Mst1*s2t*(
345  110*MuSUSY + (-17 + 30*lmMst1 - 30*lmMt)*Mst1*Tbeta)*pow2(Sbeta)) -
346  pow2(Dmst12)*((((181 - 90*lmMst1 + 90*lmMt)*MuSUSY + 800*Mst1*Tbeta)*
347  pow2(Mt) + 30*(2*MuSUSY + 11*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t))*pow2(
348  Sbeta) - 2*Mt*s2t*(-20*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(
349  330*MuSUSY + (17 + 45*lmMst1 - 45*lmMt)*Mst1*Tbeta)*pow2(Sbeta))) + 4*(
350  (48 - 45*lmMst1 + 45*lmMt)*MuSUSY + (377 - 30*lmMst1 + 30*lmMt)*Mst1*
351  Tbeta)*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) + Al4p*threeLoopFlag*(2*Dmst12*
352  Mst1*Mt*pow2(Mst2)*(11760*Dmsqst1*(1331 - 420*lmMst1 + 420*lmMt)*s2t*
353  Tbeta*pow2(Mst1)*pow2(Sbeta) + pow2(Msq)*(5*Mt*(4*Mst1*Tbeta*(16826654
354  - 1921395*lmMt - 3*lmMst1*(555463 + 2520*lmMt) + 4241160*pow2(lmMst1))
355  + 3*MuSUSY*(1417174939 - 401268*lmMt + 4*lmMst1*(4061413 + 37800*lmMt)
356  - 3185280*pow2(lmMst1) + 211680*pow2(lmMt)))*pow2(Sbeta) + 49*s2t*(-5*
357  Tbeta*(-14217821 - 161940*lmMst1 + 28800*pow2(lmMst1))*pow2(MuSUSY)*(-1
358  + pow2(Sbeta)) + 2*Mst1*(75*MuSUSY*(520781 + 17172*lmMst1 - 192*lmMt -
359  11520*pow2(lmMst1)) + 4*Mst1*Tbeta*(27088246 + 5775*lmMt + 45*lmMst1*(
360  12571 + 270*lmMt) - 136350*pow2(lmMst1) + 16200*pow2(lmMt)))*pow2(
361  Sbeta)))) - pow2(Dmst12)*(Mst1*pow2(Msq)*((10*(4*Mst1*Tbeta*(16826654 -
362  1921395*lmMt - 3*lmMst1*(555463 + 2520*lmMt) + 4241160*pow2(lmMst1)) +
363  3*MuSUSY*(1417174939 - 401268*lmMt + 4*lmMst1*(4061413 + 37800*lmMt) -
364  3185280*pow2(lmMst1) + 211680*pow2(lmMt)))*pow2(Mt) + 735*(MuSUSY*(
365  14217821 + 161940*lmMst1 - 28800*pow2(lmMst1)) + 10*Mst1*Tbeta*(-520877
366  - 17172*lmMst1 + 192*lmMt + 11520*pow2(lmMst1)))*pow2(Mst1)*pow2(s2t))*
367  pow2(Sbeta) + 2*Mt*s2t*(-245*Tbeta*(-14217821 - 161940*lmMst1 + 28800*
368  pow2(lmMst1))*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(-7350*MuSUSY*(-
369  520781 - 17172*lmMst1 + 192*lmMt + 11520*pow2(lmMst1)) + Mst1*Tbeta*(-
370  10642041163 + 11458020*lmMt + 60*lmMst1*(-346639 + 41580*lmMt) -
371  5670000*pow2(lmMst1) + 3175200*pow2(lmMt)))*pow2(Sbeta))) + 23520*
372  Dmsqst1*(1331 - 420*lmMst1 + 420*lmMt)*Mt*s2t*Tbeta*pow2(Sbeta)*pow3(
373  Mst1)) + 196*pow2(Mt)*(120*Dmsqst1*Mst1*((1541 - 420*lmMst1 + 420*lmMt)
374  *MuSUSY + (2579 + 120*lmMst1 - 120*lmMt)*Mst1*Tbeta)*pow2(Sbeta) +
375  pow2(Msq)*(5*(2219399 + 9600*lmMst1 + 960*lmMt)*Tbeta*pow2(MuSUSY)*(-1
376  + pow2(Sbeta)) + 2*Mst1*(Mst1*Tbeta*(10583177 + 60630*lmMt + 540*
377  lmMst1*(188 + 395*lmMt) + 278100*pow2(lmMst1) - 59400*pow2(lmMt)) +
378  MuSUSY*(54198467 + 43950*lmMt + 180*lmMst1*(6353 + 135*lmMt) - 272700*
379  pow2(lmMst1) + 32400*pow2(lmMt)))*pow2(Sbeta)))*pow4(Mst2))) + 141120*
380  twoLoopFlag*pow2(Mst1)*pow4(Msq)*(2*MuSUSY*pow2(Sbeta)*(-25*pow2(Mst1)*
381  (4*Dmst12*Mt*(2*(5 + 6*lmMst1 + 3*lmMt)*Mt - 9*(1 + 2*lmMst1)*Mst1*s2t)
382  *pow2(Mst2) - pow2(Dmst12)*(24*(1 - 3*lmMst1)*Mst1*Mt*s2t + 2*(5 + 6*(
383  lmMst1 + lmMt))*pow2(Mt) + 9*(1 + 2*lmMst1)*pow2(Mst1)*pow2(s2t)) + 72*
384  (1 + lmMst1 + lmMt)*pow2(Mt)*pow4(Mst2)) + 2*pow2(Dmglst1)*(3*Dmst12*
385  Mt*((327 - 30*lmMst1 + 30*lmMt)*Mt + 50*(11 - 6*lmMst1)*Mst1*s2t)*pow2(
386  Mst2) + 9*pow2(Dmst12)*(5*Mst1*Mt*s2t + (-109 + 10*lmMst1 - 10*lmMt)*
387  pow2(Mt) - 25*pow2(Mst1)*pow2(s2t)) + 100*(17 - 3*lmMst1 + 3*lmMt)*
388  pow2(Mt)*pow4(Mst2)) + Dmglst1*Mst1*(-300*Dmst12*Mt*((-3 + 6*lmMst1)*Mt
389  + 2*(-1 + 6*lmMst1)*Mst1*s2t)*pow2(Mst2) + 3*pow2(Dmst12)*(-100*Mst1*
390  Mt*s2t + 2*(-41 + 90*lmMst1 + 10*lmMt)*pow2(Mt) + 25*(-5 + 6*lmMst1)*
391  pow2(Mst1)*pow2(s2t)) - 200*(-7 + 15*lmMst1 + 3*lmMt)*pow2(Mt)*pow4(
392  Mst2))) + Tbeta*(25*Dmst12*s2t*(2*pow2(Dmglst1)*(12*Dmst12*Mt*pow2(
393  MuSUSY) - 12*Mt*pow2(Mst2)*pow2(MuSUSY) - Dmst12*(-11 + 6*lmMst1)*Mst1*
394  s2t*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 6*pow2(Mst1)*pow2(Sbeta))) - 4*
395  Dmglst1*Mst1*((5 - 6*lmMst1)*Mt*pow2(Mst2)*pow2(MuSUSY) + 2*Dmst12*
396  Mst1*s2t*((-1 + 3*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*(-1 + 6*
397  lmMst1)*pow2(Mst1)*pow2(Sbeta))) + pow2(Mst1)*(12*(1 + 2*lmMst1)*Mt*
398  pow2(Mst2)*pow2(MuSUSY) + Dmst12*(4*Mt*pow2(MuSUSY) + 3*Mst1*s2t*pow2(
399  MuSUSY)*(-1 + pow2(Sbeta)) + 36*(1 + 2*lmMst1)*s2t*pow2(Sbeta)*pow3(
400  Mst1)))) + pow2(Sbeta)*(-4*Dmst12*Mt*s2t*(2*pow2(Dmglst1)*(-25*pow2(
401  Mst2)*((31 - 6*lmMst1 + 6*lmMt)*pow2(Mst1) + 3*pow2(MuSUSY)) + Dmst12*(
402  (307 - 105*lmMst1 + 105*lmMt)*pow2(Mst1) + 75*pow2(MuSUSY))) - 25*pow2(
403  Mst1)*(Dmst12*(4*(1 + 3*lmMst1 + 6*lmMt)*pow2(Mst1) - pow2(MuSUSY)) -
404  3*pow2(Mst2)*(6*(1 + 2*(lmMst1 + lmMt))*pow2(Mst1) + (1 + 2*lmMst1)*
405  pow2(MuSUSY))) + 25*Dmglst1*(Mst1*pow2(Mst2)*(2*(-17 + 30*lmMst1 + 6*
406  lmMt)*pow2(Mst1) + (-5 + 6*lmMst1)*pow2(MuSUSY)) - 4*Dmst12*(-4 + 6*
407  lmMst1 + 3*lmMt)*pow3(Mst1))) - 2*Mst1*pow2(Mt)*(6*pow2(Dmglst1)*(35*
408  pow2(Dmst12) + 3*Dmst12*(159 - 110*lmMst1 + 10*lmMt)*pow2(Mst2) + 200*(
409  4 - 3*lmMst1)*pow4(Mst2)) - 4*Dmglst1*Mst1*((11 + 60*lmMst1 - 60*lmMt)*
410  pow2(Dmst12) + 25*Dmst12*(1 + 30*lmMst1 + 6*lmMt)*pow2(Mst2) + 100*(1 +
411  12*lmMst1 + 6*lmMt)*pow4(Mst2)) - 25*pow2(Mst1)*((5 + 42*lmMst1 + 30*
412  lmMt)*pow2(Dmst12) - 8*Dmst12*(4 + 3*lmMst1 + 6*lmMt)*pow2(Mst2) - 72*(
413  1 + lmMst1 - lmMt + 2*lmMst1*lmMt + pow2(lmMst1) - 3*pow2(lmMt))*pow4(
414  Mst2))))))) + 3969000*oneLoopFlag*pow4(Msq)*(s2t*pow2(Dmst12)*(-12*Mt*
415  MuSUSY*pow2(Sbeta) + s2t*Tbeta*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*
416  pow2(Mst1)*pow2(Sbeta))) + 12*Mt*pow2(Sbeta)*(2*Dmst12*MuSUSY*s2t*pow2(
417  Mst2) + Mt*Tbeta*(pow2(Dmst12) - 2*Dmst12*pow2(Mst2) + 4*(-lmMst1 +
418  lmMt)*pow4(Mst2))))*pow5(Mst1))) + 4*Mst1*xDmst12*pow3(Dmst12)*(Mt*(
419  992250*oneLoopFlag*(-(Mt*Tbeta*pow2(s2t)*(pow2(MuSUSY)*(-1 + pow2(
420  Sbeta)) + 18*pow2(Mst1)*pow2(Sbeta))) - pow2(Sbeta)*(-8*MuSUSY*s2t*
421  pow2(Mt) + 8*Tbeta*pow3(Mt) + MuSUSY*pow2(Mst1)*pow3(s2t)))*pow4(Msq)*
422  pow4(Mst1) + Al4p*(-35280*Mst1*twoLoopFlag*(pow2(Mst1)*(25*Mst1*Mt*
423  pow2(s2t)*((1 + 6*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 12*
424  Mst1*(MuSUSY + 3*lmMst1*MuSUSY + (1 + 12*lmMst1)*Mst1*Tbeta)*pow2(
425  Sbeta)) + 100*s2t*pow2(Mt)*(-((5 + 6*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 +
426  pow2(Sbeta))) + 2*Mst1*(-12*(-1 + lmMst1)*MuSUSY + (1 + 3*lmMst1 + 9*
427  lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-4*(50*(5 + 6*lmMst1)*
428  MuSUSY + (137 - 330*lmMst1 - 270*lmMt)*Mst1*Tbeta)*pow3(Mt) + 75*(
429  MuSUSY - 2*(1 + 2*lmMst1)*Mst1*Tbeta)*pow3(Mst1)*pow3(s2t))) + 2*
430  Dmglst1*(Dmglst1*(-5*Mst1*Mt*pow2(s2t)*(10*(-11 + 6*lmMst1)*Tbeta*pow2(
431  MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*(60*MuSUSY + (-107 + 60*lmMst1)*
432  Mst1*Tbeta)*pow2(Sbeta)) + 4*s2t*pow2(Mt)*(-75*Tbeta*pow2(MuSUSY)*(-1 +
433  pow2(Sbeta)) + Mst1*((870 - 450*lmMst1)*MuSUSY + (161 + 60*lmMst1 - 60*
434  lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-6*((327 - 30*lmMst1 +
435  30*lmMt)*MuSUSY + (547 - 330*lmMst1 + 30*lmMt)*Mst1*Tbeta)*pow3(Mt) +
436  25*((11 - 6*lmMst1)*MuSUSY + 6*Mst1*Tbeta)*pow3(Mst1)*pow3(s2t))) +
437  Mst1*(-25*Mst1*Mt*pow2(s2t)*(8*(-1 + 3*lmMst1)*Tbeta*pow2(MuSUSY)*(-1 +
438  pow2(Sbeta)) + 3*Mst1*((5 - 6*lmMst1)*MuSUSY + 6*(-1 + 4*lmMst1)*Mst1*
439  Tbeta)*pow2(Sbeta)) - 2*s2t*pow2(Mt)*(30*lmMst1*(5*Tbeta*pow2(MuSUSY)*(
440  -1 + pow2(Sbeta)) + Mst1*(60*MuSUSY - 11*Mst1*Tbeta)*pow2(Sbeta)) +
441  Tbeta*(-125*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 2*(131 - 135*lmMt)*pow2(
442  Mst1)*pow2(Sbeta))) + pow2(Sbeta)*(4*(6*(17 - 30*lmMst1 + 5*lmMt)*
443  MuSUSY + (47 + 870*lmMst1 + 30*lmMt)*Mst1*Tbeta)*pow3(Mt) + 25*((4 -
444  12*lmMst1)*MuSUSY + (5 - 6*lmMst1)*Mst1*Tbeta)*pow3(Mst1)*pow3(s2t)))))
445  *pow4(Msq) - 16*xDmglst1*pow2(Msq)*pow3(Dmglst1)*(-8820*twoLoopFlag*
446  pow2(Msq)*(10*s2t*pow2(Mt)*(-4*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
447  Mst1*(66*MuSUSY + 17*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(-60*Mt*(2*
448  MuSUSY + 11*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t) + ((-181 + 90*lmMst1 - 90*
449  lmMt)*MuSUSY - 800*Mst1*Tbeta)*pow3(Mt) + 20*Tbeta*pow3(s2t)*pow4(Mst1)
450  )) + Al4p*threeLoopFlag*(23520*Dmsqst1*(1331 - 420*lmMst1 + 420*lmMt)*
451  s2t*Tbeta*pow2(Mst1)*pow2(Mt)*pow2(Sbeta) + 2*s2t*pow2(Msq)*pow2(Mt)*(-
452  245*Tbeta*(-14217821 - 161940*lmMst1 + 28800*pow2(lmMst1))*pow2(MuSUSY)
453  *(-1 + pow2(Sbeta)) + 2*Mst1*(3675*MuSUSY*(520781 + 17172*lmMst1 - 192*
454  lmMt - 11520*pow2(lmMst1)) + Mst1*Tbeta*(-15951337379 + 10326120*lmMt +
455  120*lmMst1*(-1097288 + 945*lmMt) + 21054600*pow2(lmMst1)))*pow2(Sbeta))
456  - 5*pow2(Msq)*pow2(Sbeta)*(-294*Mt*(MuSUSY*(14217821 + 161940*lmMst1 -
457  28800*pow2(lmMst1)) + 10*Mst1*Tbeta*(-520877 - 17172*lmMst1 + 192*lmMt
458  + 11520*pow2(lmMst1)))*pow2(Mst1)*pow2(s2t) + 2*(4*Mst1*Tbeta*(-
459  16826654 + 1921395*lmMt + 3*lmMst1*(555463 + 2520*lmMt) - 4241160*pow2(
460  lmMst1)) + 3*MuSUSY*(-1417174939 + 401268*lmMt - 4*lmMst1*(4061413 +
461  37800*lmMt) + 3185280*pow2(lmMst1) - 211680*pow2(lmMt)))*pow3(Mt) + 49*
462  Tbeta*(14217821 + 161940*lmMst1 - 28800*pow2(lmMst1))*pow3(s2t)*pow4(
463  Mst1)))))) + threeLoopFlag*pow2(Al4p)*(-(pow2(Msq)*(196*Mst1*Mt*MuSUSY*
464  pow2(Sbeta)*(Dmglst1*Mst1*(2400*Dmsqst1*Mt*(7050*Mst1*Mt*s2t + 2*(557 +
465  120*lmMst1 - 120*lmMt)*pow2(Mt) + 375*pow2(Mst1)*pow2(s2t)) - pow2(Msq)
466  *(3600*Mst1*s2t*(12383 + 4128*lmMst1 + 80*lmMt - 1260*pow2(lmMst1))*
467  pow2(Mt) + 225*Mt*(284641 + 8696*lmMst1 + 1680*pow2(lmMst1))*pow2(Mst1)
468  *pow2(s2t) + 8*(193364399 + 90000*lmMt - 300*lmMst1*(3781 + 582*lmMt) -
469  2005200*pow2(lmMst1) - 43200*pow2(lmMt))*pow3(Mt) + 375*(84209 + 1264*
470  lmMst1 - 240*pow2(lmMst1))*pow3(Mst1)*pow3(s2t))) + 15*pow2(Mst1)*(500*
471  Dmsqst1*(6*Mst1*s2t*(173 - 144*lmMst1*(-1 + shiftst2) + 216*shiftst2)*
472  pow2(Mt) + 504*Mt*pow2(Mst1)*pow2(s2t) + 16*(65 - 6*lmMst1 + 6*lmMt)*
473  pow3(Mt) + 3*(7 + 72*shiftst1 - 36*shiftst2 + 24*lmMst1*(1 - 2*shiftst1
474  + shiftst2))*pow3(Mst1)*pow3(s2t)) - pow2(Msq)*(15*Mst1*s2t*(-102747 +
475  640*lmMt + 6720*shiftst3 - 32*lmMst1*(331 + 90*shiftst3) + 13888*pow2(
476  lmMst1))*pow2(Mt) - 75*Mt*(-20531 + 200*lmMst1 + 1200*pow2(lmMst1))*
477  pow2(Mst1)*pow2(s2t) - 4*(-3454599 + 16840*lmMt + 48*lmMst1*(262 + 405*
478  lmMt) + 46560*pow2(lmMst1))*pow3(Mt) + 50*(1429 - 720*shiftst1 + 360*
479  shiftst2 - 234*shiftst3 + 2*lmMst1*(-227 + 720*shiftst1 - 360*shiftst2
480  + 126*shiftst3) + 24*pow2(lmMst1))*pow3(Mst1)*pow3(s2t))) + 6*pow2(
481  Dmglst1)*(54000*Dmsqst1*Mst1*s2t*pow2(Mt) + pow2(Msq)*(-600*Mst1*s2t*(
482  274009 + 964*lmMst1 + 104*lmMt - 8310*pow2(lmMst1))*pow2(Mt) + 150*Mt*(
483  954181 + 11256*lmMst1 - 11520*pow2(lmMst1))*pow2(Mst1)*pow2(s2t) + 4*(
484  84334067 + 120*lmMst1*(2843 - 120*lmMt) + 202200*lmMt - 828000*pow2(
485  lmMst1) - 21600*pow2(lmMt))*pow3(Mt) + 25*(-655743 - 5288*lmMst1 +
486  11820*pow2(lmMst1))*pow3(Mst1)*pow3(s2t)))) + Tbeta*(3675*pow2(Mst1)*
487  pow2(Mt)*pow2(s2t)*((8*Dmglst1*(2*Dmglst1*(655743 + 5288*lmMst1 -
488  11820*pow2(lmMst1)) + 5*Mst1*(84209 + 1264*lmMst1 - 240*pow2(lmMst1)))*
489  pow2(Msq) - 5*(480*Dmsqst1*(7 + 24*lmMst1) + (70121 + 2208*lmMst1 -
490  432*pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*pow2(MuSUSY) + 1440*shiftst1*
491  pow2(Mst1)*(80*Dmsqst1*(3 - 2*lmMst1)*pow2(Mst1) + (1 - 2*lmMst1)*pow2(
492  Msq)*(80*pow2(Mst1) + 3*pow2(MuSUSY)))*pow2(Sbeta)) - 245*(16*pow2(
493  Dmglst1)*(2*Mt*(6233611 + 58800*lmMst1 - 4560*lmMt - 14400*pow2(lmMst1)
494  ) + 15*Mst1*s2t*(954181 + 11256*lmMst1 - 11520*pow2(lmMst1)))*pow2(Msq)
495  + 40*Dmglst1*Mst1*(6000*Dmsqst1*Mst1*s2t + (5*Mt*(403559 + 384*(lmMst1
496  + lmMt) - 4608*pow2(lmMst1)) + 6*Mst1*s2t*(-1282471 + 7264*lmMst1 +
497  18120*pow2(lmMst1)))*pow2(Msq)) + 15*(67200*Dmsqst1*Mst1*s2t + (80*
498  Mst1*s2t*(-36863 + 80*lmMst1 + 552*pow2(lmMst1)) + Mt*(-1763661 +
499  47104*lmMst1 - 5120*lmMt + 24576*pow2(lmMst1)))*pow2(Msq))*pow2(Mst1))*
500  pow2(MuSUSY)*pow3(Mt) + 529200*pow4(Mst1)*(-(shiftst3*pow2(Msq)*(2*
501  pow2(Mt)*pow2(s2t)*((-29 + 18*lmMst1)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
502  5*(-83 + 58*lmMst1)*pow2(Mst1)*pow2(Sbeta)) - 80*(-7 + 3*lmMst1)*pow2(
503  Sbeta)*pow4(Mt) + 5*(1 - 2*lmMst1)*pow2(Sbeta)*pow4(Mst1)*pow4(s2t))) -
504  10*(3*(1 - 2*lmMst1)*shiftst1*pow2(Msq)*pow2(Mt)*pow2(MuSUSY)*pow2(s2t)
505  - (1 - 2*lmMst1)*shiftst2*pow2(Msq)*pow2(s2t)*(2*pow2(Mt)*(pow2(MuSUSY)
506  *(-1 + pow2(Sbeta)) - 10*pow2(Mst1)*pow2(Sbeta)) + 5*pow2(s2t)*pow2(
507  Sbeta)*pow4(Mst1)) + 5*shiftst1*(Dmsqst1*(3 - 2*lmMst1) + (1 - 2*
508  lmMst1)*pow2(Msq))*pow2(Sbeta)*pow4(Mst1)*pow4(s2t) + 5*Dmsqst1*(3 - 2*
509  lmMst1)*shiftst2*(4*pow2(Mt)*pow2(s2t)*(-(pow2(MuSUSY)*(-1 + pow2(
510  Sbeta))) + pow2(Mst1)*pow2(Sbeta)) + pow2(Sbeta)*(24*pow4(Mt) - pow4(
511  Mst1)*pow4(s2t))))) + Mt*pow2(Sbeta)*(-392*Mst1*s2t*pow2(Mt)*(15*pow2(
512  Mst1)*(2000*Dmsqst1*(2*(13 + 6*lmMst1 - 6*lmMt)*pow2(Mst1) - 21*pow2(
513  MuSUSY)) - pow2(Msq)*((558619 + 76160*lmMt - 224*lmMst1*(1219 + 60*
514  lmMt) + 123840*pow2(lmMst1) + 86400*pow2(lmMt))*pow2(Mst1) + 50*(-36863
515  + 80*lmMst1 + 552*pow2(lmMst1))*pow2(MuSUSY))) + 2*pow2(Dmglst1)*(3600*
516  Dmsqst1*(397 - 30*lmMst1 + 30*lmMt)*pow2(Mst1) + pow2(Msq)*(4*(
517  192278911 + 177300*lmMt + 60*lmMst1*(19139 + 570*lmMt) - 1373400*pow2(
518  lmMst1) + 43200*pow2(lmMt))*pow2(Mst1) + 75*(-954181 - 11256*lmMst1 +
519  11520*pow2(lmMst1))*pow2(MuSUSY))) - 2*Dmglst1*Mst1*(600*Dmsqst1*(6*(47
520  - 30*lmMst1 + 30*lmMt)*pow2(Mst1) + 125*pow2(MuSUSY)) + pow2(Msq)*((
521  28188929 - 143100*lmMt - 3780*lmMst1*(549 + 80*lmMt) + 1389600*pow2(
522  lmMst1) + 388800*pow2(lmMt))*pow2(Mst1) + 75*(-1282471 + 7264*lmMst1 +
523  18120*pow2(lmMst1))*pow2(MuSUSY)))) - 3675*Mt*pow2(Mst1)*pow2(s2t)*(8*
524  pow2(Dmglst1)*(2160*Dmsqst1*pow2(Mst1) + pow2(Msq)*(3*(-452211 - 17060*
525  lmMst1 + 992*lmMt + 1320*pow2(lmMst1))*pow2(Mst1) + 2*(655743 + 5288*
526  lmMst1 - 11820*pow2(lmMst1))*pow2(MuSUSY))) - 5*pow2(Mst1)*(120*
527  Dmsqst1*((-137 + 288*lmMst1)*pow2(Mst1) + 4*(7 + 24*lmMst1)*pow2(
528  MuSUSY)) + pow2(Msq)*(24*(2785 - 304*lmMst1 + 384*lmMt + 768*pow2(
529  lmMst1))*pow2(Mst1) + (70121 + 2208*lmMst1 - 432*pow2(lmMst1))*pow2(
530  MuSUSY))) + 8*Dmglst1*(Mst1*pow2(Msq)*(24*(-20017 + 1203*lmMst1 - 200*
531  lmMt + 2250*pow2(lmMst1))*pow2(Mst1) + 5*(84209 + 1264*lmMst1 - 240*
532  pow2(lmMst1))*pow2(MuSUSY)) + 112800*Dmsqst1*pow3(Mst1))) - pow3(Mt)*(
533  16*pow2(Dmglst1)*(17640*Dmsqst1*(457 + 510*lmMst1 - 510*lmMt)*pow2(
534  Mst1) + pow2(Msq)*((-6321826673 + 60581820*lmMst1 + 4506600*lmMt +
535  9513000*lmMst1*lmMt + 357663600*pow2(lmMst1) + 6350400*pow2(lmMt))*
536  pow2(Mst1) + 490*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(
537  lmMst1))*pow2(MuSUSY))) + 392*Dmglst1*(Mst1*pow2(Msq)*(12*(9598037 +
538  92280*lmMt + 20*lmMst1*(-11207 + 270*lmMt) + 246000*pow2(lmMst1) -
539  14400*pow2(lmMt))*pow2(Mst1) - 125*(403559 + 384*(lmMst1 + lmMt) -
540  4608*pow2(lmMst1))*pow2(MuSUSY)) + 21600*Dmsqst1*(423 - 20*lmMst1 + 20*
541  lmMt)*pow3(Mst1)) + 49*(pow2(Msq)*pow2(Mst1)*((83430364 - 8607840*lmMt
542  + 480*lmMst1*(36107 + 13380*lmMt) - 9273600*pow2(lmMst1) - 6652800*
543  pow2(lmMt))*pow2(Mst1) + 75*(1763661 - 47104*lmMst1 + 5120*lmMt -
544  24576*pow2(lmMst1))*pow2(MuSUSY)) + 9600*Dmsqst1*(3268 + 2805*lmMst1 -
545  105*lmMt)*pow4(Mst1))) - 29400*(pow2(Dmglst1)*(954181 + 11256*lmMst1 -
546  11520*pow2(lmMst1))*pow2(Msq) + 10*Dmglst1*Mst1*(100*Dmsqst1 + (33261 -
547  532*lmMst1 - 660*pow2(lmMst1))*pow2(Msq)) + 20*(210*Dmsqst1 + (1361 +
548  10*lmMst1 + 54*pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*pow3(s2t)*pow5(
549  Mst1))))) + 5880*xDmsqst1*pow2(Dmsqst1)*pow2(Mst1)*(16*pow2(Dmglst1)*
550  pow2(Mt)*(Mt*s2t*(-675*MuSUSY - 30*(-397 + 30*lmMst1 - 30*lmMt)*Mst1*
551  Tbeta + 4*Dmglst1*(-1331 + 420*lmMst1 - 420*lmMt)*Tbeta*xDmglst1) + 3*(
552  457 + 510*lmMst1 - 510*lmMt)*Tbeta*pow2(Mt) + 675*Tbeta*pow2(Mst1)*
553  pow2(s2t))*pow2(Sbeta) + 40*Dmglst1*Mt*(-2*s2t*pow2(Mt)*(125*Tbeta*
554  pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(7050*MuSUSY + (-193 - 330*
555  lmMst1 + 330*lmMt)*Mst1*Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(150*Mt*(-5*
556  MuSUSY + 94*Mst1*Tbeta)*pow2(Mst1)*pow2(s2t) - 4*((557 + 120*lmMst1 -
557  120*lmMt)*MuSUSY + 9*(-423 + 20*lmMst1 - 20*lmMt)*Mst1*Tbeta)*pow3(Mt)
558  + 125*Tbeta*pow3(s2t)*pow4(Mst1))) + 5*Mst1*(-75*Mst1*pow2(Mt)*pow2(
559  s2t)*(4*(7 + 6*shiftst1 - 24*lmMst1*(-1 + shiftst2) + 30*shiftst2)*
560  Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 3*Mst1*(112*MuSUSY + Mst1*(91 +
561  240*shiftst1 + 32*lmMst1*(3 - 4*shiftst1 + shiftst2))*Tbeta)*pow2(
562  Sbeta)) - 150*s2t*(56*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + Mst1*(3*
563  MuSUSY*(47 + 96*(lmMst1 + shiftst2 - lmMst1*shiftst2)) - 208*Mst1*
564  Tbeta)*pow2(Sbeta))*pow3(Mt) + pow2(Sbeta)*(-150*Mt*(MuSUSY*(7 + 108*
565  shiftst1 - 72*shiftst2 + 24*lmMst1*(1 - 2*shiftst1 + shiftst2)) - 28*
566  Mst1*Tbeta)*pow3(Mst1)*pow3(s2t) - 4*(200*(65 - 6*lmMst1 + 6*lmMt)*
567  MuSUSY - 3*Mst1*(1999 - 90*lmMt + 90*lmMst1*(41 - 40*shiftst2) + 3600*
568  shiftst2)*Tbeta)*pow4(Mt) - 1800*(-2 + lmMst1)*(shiftst1 - shiftst2)*
569  Tbeta*pow4(s2t)*pow5(Mst1)))))) + 11025*threeLoopFlag*pow2(Al4p)*(384*
570  z2*(pow2(Mst2)*(-50*Mt*xDmsqst1*pow2(Dmsqst1)*(-(pow2(Dmst12)*pow2(s2t)
571  *(-(Mt*(shiftst1 + shiftst2)*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta))) -
572  3*(MuSUSY*s2t*(shiftst1 - shiftst2) + 2*Mt*(3*shiftst1 - shiftst2)*
573  Tbeta)*pow2(Mst1)*pow2(Sbeta))) + 2*Dmst12*Mt*pow2(Mst2)*(12*Mt*MuSUSY*
574  s2t*shiftst2*pow2(Sbeta) - 12*shiftst2*Tbeta*pow2(Mt)*pow2(Sbeta) - (
575  shiftst1 - shiftst2)*Tbeta*pow2(s2t)*(pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
576  12*pow2(Mst1)*pow2(Sbeta))) - 24*(MuSUSY*s2t*(shiftst1 - shiftst2) +
577  Mt*(shiftst1 + shiftst2)*Tbeta)*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) +
578  pow2(Msq)*(50*Tbeta*(Dmsqst1 + pow2(Msq))*pow2(Mt)*(-(pow2(Dmst12)*
579  pow2(s2t)*((shiftst1 + shiftst2)*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 6*(
580  3*shiftst1 - shiftst2)*pow2(Mst1)*pow2(Sbeta))) + 2*Dmst12*pow2(Mst2)*(
581  (shiftst1 - shiftst2)*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 12*(
582  shiftst2*pow2(Mt) + (shiftst1 - shiftst2)*pow2(Mst1)*pow2(s2t))*pow2(
583  Sbeta)) + 24*(shiftst1 + shiftst2)*pow2(Mt)*pow2(Sbeta)*pow4(Mst2)) +
584  15*MuSUSY*pow2(Sbeta)*(10*(Dmsqst1 + pow2(Msq))*(-8*Dmst12*s2t*
585  shiftst2*pow2(Mst2)*pow3(Mt) + Mt*(-shiftst1 + shiftst2)*pow2(Dmst12)*
586  pow2(Mst1)*pow3(s2t) + 8*s2t*(shiftst1 - shiftst2)*pow3(Mt)*pow4(Mst2))
587  + shiftst3*pow2(Msq)*(-(Mt*pow2(Dmst12)*pow2(Mst1)*pow3(s2t)) + 8*s2t*
588  pow3(Mt)*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2))))) + 5*
589  shiftst3*Tbeta*pow4(Msq)*(-(Dmst12*pow2(Mt)*pow2(s2t)*((3*Dmst12 - 2*
590  pow2(Mst2))*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 6*pow2(Mst1)*(7*Dmst12 -
591  4*pow2(Mst2))*pow2(Sbeta))) + 24*pow2(Sbeta)*(pow2(Dmst12) - Dmst12*
592  pow2(Mst2) + pow4(Mst2))*pow4(Mt))) - xDmst12*pow3(Dmst12)*(50*Dmsqst1*
593  (Dmsqst1*xDmsqst1 + pow2(Msq))*(-4*shiftst2*pow2(Mt)*(Tbeta*pow2(
594  MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) + 6*Mt*(MuSUSY*s2t - Mt*Tbeta)*
595  pow2(Sbeta)) + pow2(Sbeta)*(2*Mt*(MuSUSY*s2t*(-2*shiftst1 + shiftst2) +
596  2*Mt*(-4*shiftst1 + shiftst2)*Tbeta)*pow2(Mst1)*pow2(s2t) + (shiftst1 -
597  shiftst2)*Tbeta*pow4(Mst1)*pow4(s2t))) + pow4(Msq)*(2*pow2(Mt)*(-((15*
598  shiftst1 + 10*shiftst2 + 9*shiftst3)*Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 +
599  pow2(Sbeta))) + 60*Mt*shiftst3*(MuSUSY*s2t + Mt*Tbeta)*pow2(Sbeta)) +
600  pow2(Sbeta)*(-5*Mt*(MuSUSY*s2t*(40*shiftst1 - 20*shiftst2 + 7*shiftst3)
601  + 2*Mt*(80*shiftst1 - 20*shiftst2 + 29*shiftst3)*Tbeta)*pow2(Mst1)*
602  pow2(s2t) + 5*(10*shiftst1 - 10*shiftst2 + shiftst3)*Tbeta*pow4(Mst1)*
603  pow4(s2t)))))*pow5(Mst1) - 5*Mt*z3*(Mt*pow2(Mst2)*(8*Dmst12*s2t*Tbeta*
604  pow2(Msq)*pow2(Mst1)*pow2(Sbeta)*(-2*Mt*pow2(Msq)*(pow2(Mst1)*(12*
605  Dmst12*(22100*Dmglst1*Mst1 + 64701*pow2(Dmglst1) + 3290*pow2(Mst1)) +
606  96*(1963*Dmglst1*Mst1 + 7488*pow2(Dmglst1) + 106*pow2(Mst1))*pow2(Mst2)
607  ) + (-210810*pow2(Dmglst1)*(Dmst12 - pow2(Mst2)) + pow2(Mst1)*(17285*
608  Dmst12 + 6024*pow2(Mst2)) + Dmglst1*Mst1*(105405*Dmst12 + 75164*pow2(
609  Mst2)))*pow2(MuSUSY)) + Dmst12*Mst1*s2t*(18*pow2(Mst1)*(11*pow2(Msq)*(
610  40*pow2(Mst1) - 13*pow2(MuSUSY)) + 70*Dmsqst1*(4*pow2(Mst1) + pow2(
611  MuSUSY))) - Dmglst1*pow2(Msq)*(74316*Dmglst1*pow2(Mst1) + 149630*
612  Dmglst1*pow2(MuSUSY) + 50485*Mst1*pow2(MuSUSY) - 34776*pow3(Mst1)))) -
613  64*xDmglst1*pow3(Dmglst1)*pow4(Msq)*(Dmst12*Mst1*(-(Dmst12*((9*(71189*
614  MuSUSY + 1240*Mst1*Tbeta)*pow2(Mt) + 33*(4789*MuSUSY - 1738*Mst1*Tbeta)
615  *pow2(Mst1)*pow2(s2t))*pow2(Sbeta) + Mt*s2t*(105358*Tbeta*pow2(MuSUSY)*
616  (-1 + pow2(Sbeta)) + 3*Mst1*(38236*MuSUSY - 107331*Mst1*Tbeta)*pow2(
617  Sbeta)))) + Mt*pow2(Mst2)*(9*Mt*(71189*MuSUSY + 1240*Mst1*Tbeta)*pow2(
618  Sbeta) + 2*s2t*(52679*Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 6*Mst1*(
619  9559*MuSUSY + 26559*Mst1*Tbeta)*pow2(Sbeta)))) + 6*pow2(Mt)*(5511*
620  Tbeta*pow2(MuSUSY)*(-1 + pow2(Sbeta)) + 2*Mst1*(26559*MuSUSY + 5282*
621  Mst1*Tbeta)*pow2(Sbeta))*pow4(Mst2)) + Mst1*pow2(Msq)*(-3*Tbeta*pow2(
622  Mt)*pow2(Sbeta)*(20160*Dmsqst1*Dmst12*(Dmst12 + 2*pow2(Mst2))*pow4(
623  Mst1) + pow2(Msq)*pow2(Mst1)*(160*Dmst12*pow2(Mst2)*(2052*pow2(Mst1) +
624  425*pow2(MuSUSY)) + pow2(Dmst12)*(29472*pow2(Mst1) + 31963*pow2(MuSUSY)
625  ) + 512*((-214 + 24*lmMst1 - 24*lmMt)*pow2(Mst1) + 47*pow2(MuSUSY))*
626  pow4(Mst2)) + 48*Dmglst1*pow2(Msq)*(Mst1*(-(pow2(Dmst12)*(6970*pow2(
627  Mst1) + 8213*pow2(MuSUSY))) + Dmst12*pow2(Mst2)*(30912*pow2(Mst1) +
628  8213*pow2(MuSUSY)) + 8*(3356*pow2(Mst1) + 583*pow2(MuSUSY))*pow4(Mst2))
629  + Dmglst1*(4*Dmst12*pow2(Mst2)*(11332*pow2(Mst1) + 5113*pow2(MuSUSY)) -
630  pow2(Dmst12)*(45397*pow2(Mst1) + 20452*pow2(MuSUSY)) + 10*(4840*pow2(
631  Mst1) + 943*pow2(MuSUSY))*pow4(Mst2)))) - Tbeta*pow2(MuSUSY)*(-8*pow2(
632  Dmst12)*pow2(Mst1)*(5*Dmglst1*(29926*Dmglst1 + 10097*Mst1)*pow2(Msq) -
633  18*(70*Dmsqst1 - 143*pow2(Msq))*pow2(Mst1))*pow2(s2t) + Mt*pow2(Msq)*(
634  16*Dmst12*Mst1*s2t*(210810*pow2(Dmglst1)*(Dmst12 - pow2(Mst2)) - pow2(
635  Mst1)*(17285*Dmst12 + 6024*pow2(Mst2)) - Dmglst1*Mst1*(105405*Dmst12 +
636  75164*pow2(Mst2))) + 3*Mt*(96*pow2(Dmglst1)*(10226*pow2(Dmst12) -
637  10226*Dmst12*pow2(Mst2) - 4715*pow4(Mst2)) + 48*Dmglst1*Mst1*(8213*
638  pow2(Dmst12) - 8213*Dmst12*pow2(Mst2) - 4664*pow4(Mst2)) - pow2(Mst1)*(
639  31963*pow2(Dmst12) + 68000*Dmst12*pow2(Mst2) + 24064*pow4(Mst2))))) +
640  16*Mst1*MuSUSY*pow2(Sbeta)*(-6*Dmst12*Mst1*Mt*s2t*(2*Dmglst1*pow2(Msq)*
641  (28677*Dmglst1*Dmst12 + 3821*Dmst12*Mst1 + 6193*Dmglst1*pow2(Mst2) -
642  2898*Mst1*pow2(Mst2)) - 3*pow2(Mst1)*(35*Dmsqst1*(7*Dmst12 + 8*pow2(
643  Mst2)) + pow2(Msq)*(403*Dmst12 + 440*pow2(Mst2)))) + pow2(Msq)*(3*pow2(
644  Dmst12)*pow2(Mst1)*(37582*Dmglst1*Mst1 + 105405*pow2(Dmglst1) + 3012*
645  pow2(Mst1))*pow2(s2t) + pow2(Mt)*(1404*pow2(Dmglst1)*(1065*pow2(Dmst12)
646  - 1065*Dmst12*pow2(Mst2) - 512*pow4(Mst2)) - pow2(Mst1)*(51181*pow2(
647  Dmst12) + 49656*Dmst12*pow2(Mst2) + 10176*pow4(Mst2)) - 4*Dmglst1*Mst1*
648  (86833*pow2(Dmst12) + 113412*Dmst12*pow2(Mst2) + 47112*pow4(Mst2))))))
649  - 5040*xDmsqst1*pow2(Dmsqst1)*(6*Dmst12*Mt*(-5*MuSUSY*s2t + 6*Mt*Tbeta)
650  *pow2(Mst2)*pow2(Sbeta) + Tbeta*pow2(Dmst12)*pow2(s2t)*(-2*pow2(MuSUSY)
651  *(-1 + pow2(Sbeta)) - 15*pow2(Mst1)*pow2(Sbeta)) + 24*Tbeta*pow2(Mt)*
652  pow2(Sbeta)*pow4(Mst2))*pow5(Mst1)) - 2*Mst1*xDmst12*pow3(Dmst12)*(
653  2520*Dmsqst1*(Dmsqst1*xDmsqst1*(Mt*Tbeta*pow2(s2t)*(4*pow2(MuSUSY)*(-1
654  + pow2(Sbeta)) + 15*pow2(Mst1)*pow2(Sbeta)) - 36*Tbeta*pow2(Sbeta)*
655  pow3(Mt) + 2*MuSUSY*pow2(Sbeta)*(15*s2t*pow2(Mt) + pow2(Mst1)*pow3(s2t)
656  )) + pow2(Msq)*(Mt*Tbeta*pow2(s2t)*(4*pow2(MuSUSY)*(-1 + pow2(Sbeta)) +
657  pow2(Mst1)*pow2(Sbeta)) - 48*Tbeta*pow2(Sbeta)*pow3(Mt) + 2*MuSUSY*
658  pow2(Sbeta)*(22*s2t*pow2(Mt) + pow2(Mst1)*pow3(s2t))))*pow4(Mst1) +
659  pow4(Msq)*(-32*xDmglst1*pow3(Dmglst1)*(-2*s2t*pow2(Mt)*(52679*Tbeta*
660  pow2(MuSUSY)*(-1 + pow2(Sbeta)) - 237*Mst1*(-242*MuSUSY + 2031*Mst1*
661  Tbeta)*pow2(Sbeta)) + pow2(Sbeta)*(66*Mt*(-4789*MuSUSY + 1738*Mst1*
662  Tbeta)*pow2(Mst1)*pow2(s2t) - 9*(71189*MuSUSY + 1240*Mst1*Tbeta)*pow3(
663  Mt) + 52679*Tbeta*pow3(s2t)*pow4(Mst1))) - 4*Dmglst1*Mst1*(4*Mst1*
664  MuSUSY*pow2(Mt)*(142987*MuSUSY*s2t*Tbeta*(-1 + pow2(Sbeta)) + 574156*
665  Mt*pow2(Sbeta)) + 2*Mt*pow2(Mst1)*(50485*Tbeta*pow2(MuSUSY)*pow2(s2t)*(
666  -1 + pow2(Sbeta)) + 12*Mt*(4744*MuSUSY*s2t + 12729*Mt*Tbeta)*pow2(
667  Sbeta)) + s2t*(Mt*(90723*MuSUSY*s2t - 164264*Mt*Tbeta) + Mst1*s2t*(
668  50485*MuSUSY*s2t - 80628*Mt*Tbeta))*pow2(Sbeta)*pow3(Mst1) + 86*Tbeta*(
669  -1719*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 874*pow2(Sbeta)*pow3(
670  s2t)*pow5(Mst1))) - pow2(Mst1)*(16*Mst1*MuSUSY*pow2(Mt)*(20297*MuSUSY*
671  s2t*Tbeta*(-1 + pow2(Sbeta)) + 76009*Mt*pow2(Sbeta)) - Mt*pow2(Mst1)*(
672  37669*Tbeta*pow2(MuSUSY)*pow2(s2t)*(-1 + pow2(Sbeta)) - 4*Mt*(-32783*
673  MuSUSY*s2t + 33933*Mt*Tbeta)*pow2(Sbeta)) + 4*s2t*(Mt*(33783*MuSUSY*s2t
674  - 23402*Mt*Tbeta) + 18*Mst1*s2t*(143*MuSUSY*s2t - 37*Mt*Tbeta))*pow2(
675  Sbeta)*pow3(Mst1) + Tbeta*(197889*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(
676  Mt) + 24096*pow2(Sbeta)*pow3(s2t)*pow5(Mst1))) - 8*pow2(Dmglst1)*(-30*
677  Mst1*MuSUSY*pow2(Mt)*(7027*MuSUSY*s2t*Tbeta*(-1 + pow2(Sbeta)) + 49842*
678  Mt*pow2(Sbeta)) + 2*Mt*pow2(Mst1)*(74815*Tbeta*pow2(MuSUSY)*pow2(s2t)*(
679  -1 + pow2(Sbeta)) + 3*Mt*(127094*MuSUSY*s2t - 68199*Mt*Tbeta)*pow2(
680  Sbeta)) + s2t*(-6*Mt*s2t*(105405*MuSUSY + 22484*Mst1*Tbeta) + 2271672*
681  Tbeta*pow2(Mt) + 74815*Mst1*MuSUSY*pow2(s2t))*pow2(Sbeta)*pow3(Mst1) +
682  Tbeta*(-184068*pow2(MuSUSY)*(-1 + pow2(Sbeta))*pow3(Mt) + 105405*pow2(
683  Sbeta)*pow3(s2t)*pow5(Mst1))))))))/(1.90512e8*Tbeta*pow2(Sbeta)*pow4(
684  Msq)*pow5(Mst1)*pow6(Mst2));
685 }
686 
687 /**
688  * @return The off-diagonal (1, 2) = (2, 1) matrix element of the Higgs mass matrix as a double for the hierarchy 'H3q22g'
689  */
690 double H3q22g::getS12() const {
691  return -(MuSUSY*(40500*Dmst12*Mt*oneLoopFlag*s2t*pow4(Msq)*(-2*Dmst12*Mt*(
692  MuSUSY*s2t + 6*Mt*Tbeta)*pow2(Mst2) + pow2(Dmst12)*(2*Mt*MuSUSY*s2t +
693  8*Tbeta*pow2(Mt) - Tbeta*pow2(Mst1)*pow2(s2t)) + 24*Tbeta*pow2(Mt)*
694  pow4(Mst2))*pow5(Mst1) + 32*Al4p*pow2(Msq)*(4*xDmglst1*pow3(Dmglst1)*(
695  45*Mst1*twoLoopFlag*pow2(Msq)*pow2(Mt)*(80*Dmst12*Mt*MuSUSY*s2t*(pow2(
696  Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) - Tbeta*(pow2(Dmst12)*pow2(
697  Mst2)*(660*Mst1*Mt*s2t + (-181 + 90*lmMst1 - 90*lmMt)*pow2(Mt) - 60*
698  pow2(Mst1)*pow2(s2t)) + (-660*Mst1*Mt*s2t + (181 - 90*lmMst1 + 90*lmMt)
699  *pow2(Mt) + 120*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + Dmst12*Mt*((181 -
700  90*lmMst1 + 90*lmMt)*Mt - 660*Mst1*s2t)*pow4(Mst2) + 12*(16 - 15*lmMst1
701  + 15*lmMt)*pow2(Mt)*pow6(Mst2))) + Al4p*threeLoopFlag*(5*MuSUSY*pow2(
702  Msq)*pow3(Mt)*(Dmst12*Mst1*s2t*(14217821 + 161940*lmMst1 - 28800*pow2(
703  lmMst1))*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + 2*(2219399 +
704  9600*lmMst1 + 960*lmMt)*Mt*pow6(Mst2)) - (Mst1*Tbeta*pow2(Mt)*(pow2(
705  Msq)*(-15*pow2(Dmst12)*pow2(Mst2)*(980*Mst1*Mt*s2t*(520781 + 17172*
706  lmMst1 - 192*lmMt - 11520*pow2(lmMst1)) + (2834349878 - 802536*lmMt +
707  8*lmMst1*(4061413 + 37800*lmMt) - 6370560*pow2(lmMst1) + 423360*pow2(
708  lmMt))*pow2(Mt) + 49*(14217821 + 161940*lmMst1 - 28800*pow2(lmMst1))*
709  pow2(Mst1)*pow2(s2t)) + 30*((490*Mst1*Mt*s2t*(520781 + 17172*lmMst1 -
710  192*lmMt - 11520*pow2(lmMst1)) + (1417174939 - 401268*lmMt + 4*lmMst1*(
711  4061413 + 37800*lmMt) - 3185280*pow2(lmMst1) + 211680*pow2(lmMt))*pow2(
712  Mt) + 49*(14217821 + 161940*lmMst1 - 28800*pow2(lmMst1))*pow2(Mst1)*
713  pow2(s2t))*pow3(Dmst12) + Dmst12*Mt*(490*Mst1*s2t*(520781 + 17172*
714  lmMst1 - 192*lmMt - 11520*pow2(lmMst1)) + Mt*(1417174939 - 401268*lmMt
715  + 4*lmMst1*(4061413 + 37800*lmMt) - 3185280*pow2(lmMst1) + 211680*pow2(
716  lmMt)))*pow4(Mst2))) + 392*(60*Dmsqst1*(1541 - 420*lmMst1 + 420*lmMt) +
717  (54198467 + 43950*lmMt + 180*lmMst1*(6353 + 135*lmMt) - 272700*pow2(
718  lmMst1) + 32400*pow2(lmMt))*pow2(Msq))*pow2(Mt)*pow6(Mst2)))/196.)) +
719  45*Mt*twoLoopFlag*pow2(Msq)*pow2(Mst1)*(-25*s2t*pow3(Dmst12)*(8*(5 + 6*
720  lmMst1)*MuSUSY*pow2(Mst1)*pow2(Mt) + 2*pow2(Dmglst1)*(24*MuSUSY*pow2(
721  Mt) + (11 - 6*lmMst1)*Tbeta*pow2(s2t)*pow3(Mst1)) - 8*Dmglst1*((5 - 6*
722  lmMst1)*Mst1*MuSUSY*pow2(Mt) + (-1 + 3*lmMst1)*Tbeta*pow2(s2t)*pow4(
723  Mst1)) + 3*Tbeta*pow2(s2t)*pow5(Mst1)) + Mt*(50*Dmst12*MuSUSY*s2t*(
724  Mst1*s2t*pow2(Dmst12)*(16*Dmglst1*(1 - 3*lmMst1)*Mst1 + (44 - 24*
725  lmMst1)*pow2(Dmglst1) + (1 + 6*lmMst1)*pow2(Mst1)) + 4*Dmst12*Mt*(6*
726  pow2(Dmglst1) + pow2(Mst1))*pow2(Mst2) - Dmst12*Mst1*s2t*(8*Dmglst1*(1
727  - 3*lmMst1)*Mst1 + (22 - 12*lmMst1)*pow2(Dmglst1) + 3*pow2(Mst1))*pow2(
728  Mst2) - 4*Mt*(Dmglst1*(5 - 6*lmMst1)*Mst1 + 6*pow2(Dmglst1) - 3*(1 + 2*
729  lmMst1)*pow2(Mst1))*pow4(Mst2)) + 2*Tbeta*(-25*pow2(Mst1)*(-(pow2(
730  Dmst12)*pow2(Mst2)*(24*(1 - 3*lmMst1)*Mst1*Mt*s2t + 2*(5 + 6*lmMst1 +
731  6*lmMt)*pow2(Mt) + 9*(1 + 2*lmMst1)*pow2(Mst1)*pow2(s2t))) + (-48*(-1 +
732  lmMst1)*Mst1*Mt*s2t - 4*(5 + 6*lmMst1)*pow2(Mt) + 6*(1 + 3*lmMst1)*
733  pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + 4*Dmst12*Mt*(2*(5 + 6*lmMst1 + 3*
734  lmMt)*Mt - 9*(1 + 2*lmMst1)*Mst1*s2t)*pow4(Mst2) + 72*(1 + lmMst1 +
735  lmMt)*pow2(Mt)*pow6(Mst2)) + 2*pow2(Dmglst1)*(9*pow2(Dmst12)*pow2(Mst2)
736  *(5*Mst1*Mt*s2t + (-109 + 10*lmMst1 - 10*lmMt)*pow2(Mt) - 25*pow2(Mst1)
737  *pow2(s2t)) + (60*(-29 + 15*lmMst1)*Mst1*Mt*s2t + (981 - 90*lmMst1 +
738  90*lmMt)*pow2(Mt) + 450*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + 3*Dmst12*
739  Mt*((327 - 30*lmMst1 + 30*lmMt)*Mt + 50*(11 - 6*lmMst1)*Mst1*s2t)*pow4(
740  Mst2) + 100*(17 - 3*lmMst1 + 3*lmMt)*pow2(Mt)*pow6(Mst2)) + Dmglst1*
741  Mst1*(-3*pow2(Dmst12)*(pow2(Mst2)*(100*Mst1*Mt*s2t - 2*(-41 + 90*lmMst1
742  + 10*lmMt)*pow2(Mt) + 25*(5 - 6*lmMst1)*pow2(Mst1)*pow2(s2t)) + Dmst12*
743  (-1200*lmMst1*Mst1*Mt*s2t + 8*(17 - 30*lmMst1 + 5*lmMt)*pow2(Mt) + 25*(
744  -5 + 6*lmMst1)*pow2(Mst1)*pow2(s2t))) + 300*Dmst12*Mt*((3 - 6*lmMst1)*
745  Mt + 2*(1 - 6*lmMst1)*Mst1*s2t)*pow4(Mst2) - 200*(-7 + 15*lmMst1 + 3*
746  lmMt)*pow2(Mt)*pow6(Mst2)))))) + threeLoopFlag*pow2(Al4p)*(4*Mst1*pow2(
747  Msq)*((5*MuSUSY*pow2(Mt)*(16*pow2(Dmglst1)*pow2(Msq)*(pow2(Dmst12)*
748  pow2(Mst2)*(30*Mst1*Mt*s2t*(-954181 - 11256*lmMst1 + 11520*pow2(lmMst1)
749  ) + 4*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(lmMst1))*pow2(
750  Mt) + 15*(655743 + 5288*lmMst1 - 11820*pow2(lmMst1))*pow2(Mst1)*pow2(
751  s2t)) - 2*(15*Mst1*Mt*s2t*(-954181 - 11256*lmMst1 + 11520*pow2(lmMst1))
752  + 2*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(lmMst1))*pow2(Mt)
753  + 15*(655743 + 5288*lmMst1 - 11820*pow2(lmMst1))*pow2(Mst1)*pow2(s2t))*
754  pow3(Dmst12) + 2*Dmst12*Mt*(15*Mst1*s2t*(954181 + 11256*lmMst1 - 11520*
755  pow2(lmMst1)) - 2*Mt*(-6233611 - 58800*lmMst1 + 4560*lmMt + 14400*pow2(
756  lmMst1)))*pow4(Mst2) - 50*(-230009 - 4704*lmMst1 + 480*lmMt + 1152*
757  pow2(lmMst1))*pow2(Mt)*pow6(Mst2)) + 40*Dmglst1*Mst1*(-(pow2(Dmst12)*
758  pow2(Mst2)*(12000*Dmsqst1*Mst1*Mt*s2t + pow2(Msq)*(6*Mst1*Mt*s2t*(-
759  949861 + 1944*lmMst1 + 11520*pow2(lmMst1)) + 10*(403559 + 384*(lmMst1 +
760  lmMt) - 4608*pow2(lmMst1))*pow2(Mt) + 15*(-84209 - 1264*lmMst1 + 240*
761  pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))) + 2*(6000*Dmsqst1*Mst1*Mt*s2t +
762  pow2(Msq)*(6*Mst1*Mt*s2t*(-1282471 + 7264*lmMst1 + 18120*pow2(lmMst1))
763  + 5*(403559 + 384*(lmMst1 + lmMt) - 4608*pow2(lmMst1))*pow2(Mt) + 15*(-
764  84209 - 1264*lmMst1 + 240*pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))*pow3(
765  Dmst12) + 10*Dmst12*Mt*(1200*Dmsqst1*Mst1*s2t + (Mt*(403559 + 384*(
766  lmMst1 + lmMt) - 4608*pow2(lmMst1)) - 12*Mst1*s2t*(-33261 + 532*lmMst1
767  + 660*pow2(lmMst1)))*pow2(Msq))*pow4(Mst2) + 240*(9631 + 16*lmMst1 +
768  48*lmMt - 192*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*pow6(Mst2)) + 15*pow2(
769  Mst1)*(-(pow2(Dmst12)*pow2(Mst2)*(2400*Dmsqst1*Mst1*s2t*(56*Mt + Mst1*
770  s2t*(7 - 12*lmMst1*(-2 + shiftst1 + shiftst2) + 18*(shiftst1 +
771  shiftst2))) - pow2(Msq)*(-240*Mst1*Mt*s2t*(-10473 + 40*lmMst1 + 256*
772  pow2(lmMst1)) + (852541 + 9216*lmMst1 - 10240*lmMt + 6144*pow2(lmMst1))
773  *pow2(Mt) + 80*(1429 - 454*lmMst1 - 180*(shiftst1 + shiftst2) + 360*
774  lmMst1*(shiftst1 + shiftst2) - 126*shiftst3 + 108*lmMst1*shiftst3 + 24*
775  pow2(lmMst1))*pow2(Mst1)*pow2(s2t)))) + 2*(-2400*Dmsqst1*Mst1*s2t*(-28*
776  Mt + Mst1*s2t*(-7 + 24*lmMst1*(-1 + shiftst2) - 36*shiftst2)) + pow2(
777  Msq)*(80*Mst1*Mt*s2t*(-36863 + 80*lmMst1 + 552*pow2(lmMst1)) + (-
778  1763661 + 47104*lmMst1 - 5120*lmMt + 24576*pow2(lmMst1))*pow2(Mt) + (
779  350605 + 4320*shiftst1 + 2880*shiftst2 + 8352*shiftst3 - 96*lmMst1*(-
780  115 + 90*shiftst1 + 60*shiftst2 + 54*shiftst3) - 2160*pow2(lmMst1))*
781  pow2(Mst1)*pow2(s2t)))*pow3(Dmst12) + 160*Dmst12*(60*Dmsqst1*Mst1*s2t*(
782  14*Mt + 3*(3 - 2*lmMst1)*Mst1*s2t*(shiftst1 - shiftst2)) + pow2(Msq)*(
783  4*Mst1*Mt*s2t*(1361 + 10*lmMst1 + 54*pow2(lmMst1)) + (11389 - 704*
784  lmMst1 + 192*lmMt - 384*pow2(lmMst1))*pow2(Mt) + 18*(1 - 2*lmMst1)*(10*
785  shiftst1 - 10*shiftst2 + shiftst3)*pow2(Mst1)*pow2(s2t)))*pow4(Mst2) +
786  1920*(349 - 56*lmMst1 + 24*lmMt - 32*pow2(lmMst1))*pow2(Msq)*pow2(Mt)*
787  pow6(Mst2))))/2. - 2*Mst1*Tbeta*(225*Dmst12*Mst1*s2t*pow2(Mt)*(Dmst12*
788  Mst1*s2t*(2*pow2(Dmglst1)*(954181 + 11256*lmMst1 - 11520*pow2(lmMst1))*
789  pow2(Msq)*(2*Dmst12 - pow2(Mst2)) + 5*pow2(Mst1)*(1680*Dmsqst1*(2*
790  Dmst12 - pow2(Mst2)) + pow2(Msq)*(Dmst12*(-20531 + 200*lmMst1 + 1200*
791  pow2(lmMst1)) - 8*(1361 + 10*lmMst1 + 54*pow2(lmMst1))*pow2(Mst2))) +
792  Dmglst1*Mst1*(2000*Dmsqst1*(2*Dmst12 - pow2(Mst2)) - pow2(Msq)*(Dmst12*
793  (284641 + 8696*lmMst1 + 1680*pow2(lmMst1)) - 20*(-33261 + 532*lmMst1 +
794  660*pow2(lmMst1))*pow2(Mst2)))) + Mt*(4*pow2(Dmglst1)*(360*Dmsqst1*(
795  pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(Msq)*(4*pow2(
796  Dmst12)*(-274009 - 964*lmMst1 - 104*lmMt + 8310*pow2(lmMst1)) + Dmst12*
797  (515917 + 6972*lmMst1 - 192*lmMt - 11520*pow2(lmMst1))*pow2(Mst2) + 2*(
798  32101 - 5044*lmMst1 + 400*lmMt - 5100*pow2(lmMst1))*pow4(Mst2))) +
799  pow2(Mst1)*(100*Dmsqst1*((346 + 288*lmMst1)*pow2(Dmst12) - 161*Dmst12*
800  pow2(Mst2) - 24*(1 + 12*lmMst1)*pow4(Mst2)) + pow2(Msq)*(pow2(Dmst12)*(
801  102747 + 10592*lmMst1 - 640*lmMt - 13888*pow2(lmMst1)) + 20*Dmst12*(-
802  2071 - 296*lmMst1 + 96*lmMt + 600*pow2(lmMst1))*pow2(Mst2) - 160*(631 -
803  lmMst1 + 36*lmMt + 21*pow2(lmMst1))*pow4(Mst2))) + 16*Dmglst1*Mst1*(
804  4700*Dmsqst1*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(
805  Msq)*(pow2(Dmst12)*(-12383 - 4128*lmMst1 - 80*lmMt + 1260*pow2(lmMst1))
806  + Dmst12*(17539 + 574*lmMst1 + 160*lmMt - 1920*pow2(lmMst1))*pow2(Mst2)
807  + 5*(-4539 + 596*lmMst1 - 48*lmMt + 516*pow2(lmMst1))*pow4(Mst2))))) +
808  2*pow4(Mt)*(-3*pow2(Dmst12)*(4*pow2(Dmglst1)*(84334067 + 120*lmMst1*(
809  2843 - 120*lmMt) + 202200*lmMt - 828000*pow2(lmMst1) - 21600*pow2(lmMt)
810  )*pow2(Msq) + 2*Dmglst1*Mst1*(400*Dmsqst1*(557 + 120*lmMst1 - 120*lmMt)
811  + (-39474953 + lmMst1*(52780 - 87600*lmMt) - 3300*lmMt + 319200*pow2(
812  lmMst1) + 14400*pow2(lmMt))*pow2(Msq)) + 5*(4000*Dmsqst1*(65 - 6*lmMst1
813  + 6*lmMt) + (-2284899 + 49840*lmMt - 32*lmMst1*(-1793 + 555*lmMt) +
814  87360*pow2(lmMst1) + 28800*pow2(lmMt))*pow2(Msq))*pow2(Mst1))*pow2(
815  Mst2) + 2*(6*pow2(Dmglst1)*(84334067 + 120*lmMst1*(2843 - 120*lmMt) +
816  202200*lmMt - 828000*pow2(lmMst1) - 21600*pow2(lmMt))*pow2(Msq) + 2*
817  Dmglst1*Mst1*(600*Dmsqst1*(557 + 120*lmMst1 - 120*lmMt) + (-193364399 -
818  90000*lmMt + 300*lmMst1*(3781 + 582*lmMt) + 2005200*pow2(lmMst1) +
819  43200*pow2(lmMt))*pow2(Msq)) + 15*(2000*Dmsqst1*(65 - 6*lmMst1 + 6*
820  lmMt) + (-3454599 + 16840*lmMt + 48*lmMst1*(262 + 405*lmMt) + 46560*
821  pow2(lmMst1))*pow2(Msq))*pow2(Mst1))*pow3(Dmst12) + 40*Dmst12*((3*pow2(
822  Dmglst1)*(84334067 + 120*lmMst1*(2843 - 120*lmMt) + 202200*lmMt -
823  828000*pow2(lmMst1) - 21600*pow2(lmMt))*pow2(Msq))/10. + 2*Dmglst1*
824  Mst1*(30*Dmsqst1*(557 + 120*lmMst1 - 120*lmMt) + (3746977 + 4005*lmMt -
825  18*lmMst1*(2711 + 1215*lmMt) - 52380*pow2(lmMst1))*pow2(Msq)) + 75*(20*
826  Dmsqst1*(65 - 6*lmMst1 + 6*lmMt) + (11697 + 448*lmMst1 + 330*lmMt -
827  372*lmMst1*lmMt + 408*pow2(lmMst1) + 288*pow2(lmMt))*pow2(Msq))*pow2(
828  Mst1))*pow4(Mst2) + 160*(pow2(Dmglst1)*(90*Dmsqst1*(-206 + 15*lmMst1 -
829  15*lmMt) + (3030652 + 19305*lmMt - 6*lmMst1*(1183 + 645*lmMt) - 55530*
830  pow2(lmMst1) - 5400*pow2(lmMt))*pow2(Msq)) + 75*Dmglst1*Mst1*(10*
831  Dmsqst1*(26 + 3*lmMst1 - 3*lmMt) + (9961 + 66*lmMt - 10*lmMst1*(67 +
832  24*lmMt) + 42*pow2(lmMst1) + 72*pow2(lmMt))*pow2(Msq)) + 150*(15*
833  Dmsqst1*(20 - 3*lmMst1 + 3*lmMt) + (434 - 83*lmMst1 + 174*lmMt - 66*
834  lmMst1*lmMt + 183*pow2(lmMst1) + 108*pow2(lmMt))*pow2(Msq))*pow2(Mst1))
835  *pow6(Mst2)) + 750*pow3(Mst1)*(-36*shiftst2*(5*Mt*pow2(Dmst12)*(
836  Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*pow2(Mst1)*(2*Dmst12
837  - 3*pow2(Mst2))*pow3(s2t) + 120*s2t*pow3(Mt)*((1 - 2*lmMst1)*pow2(Msq)*
838  (Dmst12 + pow2(Mst2))*pow4(Mst2) + Dmsqst1*(-3 + 2*lmMst1)*(pow3(
839  Dmst12) - Dmst12*pow4(Mst2) - pow6(Mst2)))) + Mt*s2t*(-(((pow2(Dmglst1)
840  *(655743 + 5288*lmMst1 - 11820*pow2(lmMst1))*pow2(Msq))/5. + (Dmglst1*
841  Mst1*(84209 + 1264*lmMst1 - 240*pow2(lmMst1))*pow2(Msq))/2. - 30*
842  Dmsqst1*(7 + 24*lmMst1)*pow2(Mst1) + (1429 - 454*lmMst1 + 24*pow2(
843  lmMst1))*pow2(Msq)*pow2(Mst1))*pow2(s2t)*pow3(Dmst12)) + 180*shiftst1*(
844  Dmsqst1*(3 - 2*lmMst1) + (1 - 2*lmMst1)*pow2(Msq))*(pow2(Dmst12)*pow2(
845  Mst1)*(4*Dmst12 - 3*pow2(Mst2))*pow2(s2t) + 24*pow2(Mt)*pow6(Mst2)) -
846  18*shiftst3*pow2(Msq)*(pow2(Dmst12)*pow2(Mst1)*(Dmst12*(-13 + 14*
847  lmMst1) + 3*(1 - 2*lmMst1)*pow2(Mst2))*pow2(s2t) + 8*pow2(Mt)*(6*(-2 +
848  lmMst1)*pow2(Dmst12)*pow2(Mst2) + (14 - 6*lmMst1)*pow3(Dmst12) + 3*
849  Dmst12*(3 - 2*lmMst1)*pow4(Mst2) + 3*(-1 + 2*lmMst1)*pow6(Mst2))))))) +
850  75*(-576*z2*pow5(Mst1)*(-50*Mt*s2t*xDmsqst1*pow2(Dmsqst1)*(s2t*pow2(
851  Dmst12)*(2*Mt*MuSUSY*(shiftst1 + shiftst2) + 3*s2t*(-shiftst1 +
852  shiftst2)*Tbeta*pow2(Mst1))*pow2(Mst2) + (-8*Mt*MuSUSY*s2t*shiftst2 +
853  24*shiftst2*Tbeta*pow2(Mt) + 2*(2*shiftst1 - shiftst2)*Tbeta*pow2(Mst1)
854  *pow2(s2t))*pow3(Dmst12) - 4*Dmst12*Mt*(MuSUSY*s2t*(shiftst1 -
855  shiftst2) + 6*Mt*shiftst2*Tbeta)*pow4(Mst2) + 24*(shiftst1 - shiftst2)*
856  Tbeta*pow2(Mt)*pow6(Mst2)) + pow2(Msq)*(2*Dmst12*MuSUSY*pow2(Mt)*pow2(
857  s2t)*(2*pow2(Dmst12)*(100*Dmsqst1*shiftst2 + (15*shiftst1 + 10*shiftst2
858  + 9*shiftst3)*pow2(Msq)) - 5*Dmst12*(10*Dmsqst1*(shiftst1 + shiftst2) +
859  (10*(shiftst1 + shiftst2) + 3*shiftst3)*pow2(Msq))*pow2(Mst2) + 10*(10*
860  Dmsqst1*(shiftst1 - shiftst2) + (10*shiftst1 - 10*shiftst2 + shiftst3)*
861  pow2(Msq))*pow4(Mst2)) + 5*shiftst3*Tbeta*pow2(Msq)*(-(Mt*pow2(Dmst12)*
862  pow2(Mst1)*(7*Dmst12 - 3*pow2(Mst2))*pow3(s2t)) + 24*s2t*pow3(Mt)*(-(
863  pow2(Dmst12)*pow2(Mst2)) + pow3(Dmst12) + Dmst12*pow4(Mst2) - pow6(
864  Mst2))) - 50*Tbeta*(shiftst2*(-(Mt*pow2(Dmst12)*(Dmsqst1 + pow2(Msq))*
865  pow2(Mst1)*(2*Dmst12 - 3*pow2(Mst2))*pow3(s2t)) + 24*s2t*pow3(Mt)*(
866  Dmsqst1*pow3(Dmst12) - (Dmsqst1 + pow2(Msq))*(Dmst12 + pow2(Mst2))*
867  pow4(Mst2))) + shiftst1*(Dmsqst1 + pow2(Msq))*(Mt*pow2(Dmst12)*pow2(
868  Mst1)*(4*Dmst12 - 3*pow2(Mst2))*pow3(s2t) + 24*s2t*pow3(Mt)*pow6(Mst2))
869  ))) + Mt*((32*Mst1*xDmsqst1*pow2(Dmsqst1)*(-20*Dmglst1*Mt*pow2(Mst1)*(-
870  (pow2(Dmst12)*pow2(Mst2)*(-100*Mt*s2t*(5*MuSUSY - 141*Mst1*Tbeta) + 4*(
871  557 + 120*lmMst1 - 120*lmMt)*Tbeta*pow2(Mt) + 375*Tbeta*pow2(Mst1)*
872  pow2(s2t))) + 2*(-50*Mt*s2t*(5*MuSUSY - 141*Mst1*Tbeta) + 2*(557 + 120*
873  lmMst1 - 120*lmMt)*Tbeta*pow2(Mt) + 375*Tbeta*pow2(Mst1)*pow2(s2t))*
874  pow3(Dmst12) + 4*Dmst12*Mt*(-125*MuSUSY*s2t + (557 + 120*lmMst1 - 120*
875  lmMt)*Mt*Tbeta + 3525*Mst1*s2t*Tbeta)*pow4(Mst2) + 100*(44 + 3*lmMst1 -
876  3*lmMt)*Tbeta*pow2(Mt)*pow6(Mst2)) - 125*pow3(Mst1)*(-2*pow2(Dmst12)*
877  pow2(Mst2)*(-6*s2t*(28*MuSUSY + Mst1*(11 + 36*shiftst2)*Tbeta)*pow2(Mt)
878  - 3*Mst1*Mt*(MuSUSY*(7 + 30*shiftst1 + 6*shiftst2 - 12*lmMst1*(-2 +
879  shiftst1 + shiftst2)) - 42*Mst1*Tbeta)*pow2(s2t) + Tbeta*(8*(65 - 6*
880  lmMst1 + 6*lmMt)*pow3(Mt) - 54*(-2 + lmMst1)*(shiftst1 - shiftst2)*
881  pow3(Mst1)*pow3(s2t))) + pow3(Dmst12)*(-(s2t*(336*MuSUSY - 9*Mst1*(47 +
882  96*(lmMst1 + shiftst2 - lmMst1*shiftst2))*Tbeta)*pow2(Mt)) - 12*Mst1*
883  Mt*(MuSUSY*(7 + 6*shiftst1 - 24*lmMst1*(-1 + shiftst2) + 30*shiftst2) -
884  42*Mst1*Tbeta)*pow2(s2t) + Tbeta*(16*(65 - 6*lmMst1 + 6*lmMt)*pow3(Mt)
885  + 3*(7 + 108*shiftst1 - 72*shiftst2 + 24*lmMst1*(1 - 2*shiftst1 +
886  shiftst2))*pow3(Mst1)*pow3(s2t))) - Dmst12*Mt*(3*Mt*s2t*(112*MuSUSY +
887  Mst1*(229 - 288*lmMst1*(-1 + shiftst2) + 576*shiftst2)*Tbeta) + 16*(-65
888  + 6*lmMst1 - 6*lmMt)*Tbeta*pow2(Mt) - 144*(-2 + lmMst1)*Mst1*MuSUSY*(
889  shiftst1 - shiftst2)*pow2(s2t))*pow4(Mst2) + 16*((91 - 12*lmMst1 + 12*
890  lmMt)*Mt - 54*(-2 + lmMst1)*Mst1*s2t*(shiftst1 - shiftst2))*Tbeta*pow2(
891  Mt)*pow6(Mst2)) + 8*Tbeta*pow2(Dmglst1)*pow2(Mt)*(4*Dmglst1*(-1541 +
892  420*lmMst1 - 420*lmMt)*Mt*xDmglst1*pow6(Mst2) + 15*Mst1*(-45*Dmst12*
893  Mst1*s2t*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + 4*(206 - 15*
894  lmMst1 + 15*lmMt)*Mt*pow6(Mst2)))))/5. - 15*z3*(-5040*Dmst12*s2t*
895  xDmsqst1*pow2(Dmsqst1)*(2*Dmst12*Mt*MuSUSY*s2t*pow2(Mst2) + pow2(
896  Dmst12)*(-4*Mt*MuSUSY*s2t + 15*Tbeta*pow2(Mt) + Tbeta*pow2(Mst1)*pow2(
897  s2t)) - 15*Tbeta*pow2(Mt)*pow4(Mst2))*pow5(Mst1) + 32*Mt*xDmglst1*pow3(
898  Dmglst1)*pow4(Msq)*(44*Mt*MuSUSY*(4789*Dmst12*Mst1*s2t*(pow2(Dmst12) -
899  Dmst12*pow2(Mst2) + pow4(Mst2)) + 1503*Mt*pow6(Mst2)) - 3*Mst1*Tbeta*(
900  52679*pow2(Dmst12)*pow2(Mst1)*(2*Dmst12 - pow2(Mst2))*pow2(s2t) +
901  38236*Dmst12*Mst1*Mt*s2t*(pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)
902  ) + 3*pow2(Mt)*(-71189*pow2(Dmst12)*pow2(Mst2) + 71189*Dmst12*(pow2(
903  Dmst12) + pow4(Mst2)) + 35412*pow6(Mst2)))) + Mst1*pow2(Msq)*(Mt*
904  MuSUSY*(-2*pow2(Dmst12)*pow2(Mst1)*(20*Dmglst1*(29926*Dmglst1 + 10097*
905  Mst1)*pow2(Msq)*(2*Dmst12 - pow2(Mst2)) - pow2(Mst1)*(5040*Dmsqst1*(2*
906  Dmst12 - pow2(Mst2)) + pow2(Msq)*(37669*Dmst12 + 10296*pow2(Mst2))))*
907  pow2(s2t) + Mt*pow2(Msq)*(16*Dmst12*Mst1*s2t*(210810*pow2(Dmglst1)*(
908  pow2(Dmst12) - Dmst12*pow2(Mst2) + pow4(Mst2)) + pow2(Mst1)*(-40594*
909  pow2(Dmst12) + 17285*Dmst12*pow2(Mst2) + 6024*pow4(Mst2)) + Dmglst1*
910  Mst1*(-285974*pow2(Dmst12) + 105405*Dmst12*pow2(Mst2) + 75164*pow4(
911  Mst2))) + 3*Mt*(-(pow2(Mst1)*(-31963*pow2(Dmst12)*pow2(Mst2) + 131926*
912  pow3(Dmst12) - 68000*Dmst12*pow4(Mst2) - 24064*pow6(Mst2))) + 48*
913  Dmglst1*Mst1*(-8213*pow2(Dmst12)*pow2(Mst2) + 8213*Dmst12*(pow2(Dmst12)
914  + pow4(Mst2)) + 4664*pow6(Mst2)) + 96*pow2(Dmglst1)*(-10226*pow2(
915  Dmst12)*pow2(Mst2) + 10226*Dmst12*(pow2(Dmst12) + pow4(Mst2)) + 4715*
916  pow6(Mst2))))) - 4*Mst1*Tbeta*(-((5*Dmglst1*(29926*Dmglst1 + 10097*
917  Mst1)*pow2(Msq) - 18*(70*Dmsqst1 - 143*pow2(Msq))*pow2(Mst1))*pow3(
918  Dmst12)*pow3(Mst1)*pow3(s2t)) - Dmst12*Mst1*s2t*pow2(Mt)*(-(pow2(Mst1)*
919  (pow2(Msq)*(32783*pow2(Dmst12) - 14508*Dmst12*pow2(Mst2) - 15840*pow4(
920  Mst2)) + 1260*Dmsqst1*(22*pow2(Dmst12) - 7*Dmst12*pow2(Mst2) - 8*pow4(
921  Mst2)))) + 24*Dmglst1*pow2(Msq)*(11*Dmglst1*(5777*pow2(Dmst12) - 2607*
922  Dmst12*pow2(Mst2) - 563*pow4(Mst2)) + Mst1*(4744*pow2(Dmst12) - 3821*
923  Dmst12*pow2(Mst2) + 2898*pow4(Mst2)))) + pow2(Msq)*(3*Mt*pow2(Dmst12)*
924  pow2(Mst1)*(210810*pow2(Dmglst1)*(2*Dmst12 - pow2(Mst2)) - pow2(Mst1)*(
925  11261*Dmst12 + 6024*pow2(Mst2)) - Dmglst1*Mst1*(30241*Dmst12 + 75164*
926  pow2(Mst2)))*pow2(s2t) + 2*pow3(Mt)*(1404*pow2(Dmglst1)*(-1065*pow2(
927  Dmst12)*pow2(Mst2) + 1065*pow3(Dmst12) + 1065*Dmst12*pow4(Mst2) + 512*
928  pow6(Mst2)) + pow2(Mst1)*(51181*pow2(Dmst12)*pow2(Mst2) - 152018*pow3(
929  Dmst12) + 49656*Dmst12*pow4(Mst2) + 10176*pow6(Mst2)) + 4*Dmglst1*Mst1*
930  (86833*pow2(Dmst12)*pow2(Mst2) - 287078*pow3(Dmst12) + 113412*Dmst12*
931  pow4(Mst2) + 47112*pow6(Mst2))))))))))))/(3.888e6*Tbeta*pow4(Msq)*pow5(
932  Mst1)*pow6(Mst2));
933 }
934 
935 /**
936  * @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)
937  */
938 double H3q22g::calc_coef_at_as2_no_sm_logs_log0() const {
939 
940  const double result =
941  (8*pow3(Dmglst1)*(Mt*pow2(Dmst12)*pow2(Mst2)*(31305120*Mst1*Mt*s2t*pow2(
942  Dmsqst1) + 31305120*Dmsqst1*Mst1*Mt*s2t*pow2(Msq) + (Mst1*Mt*s2t*(-
943  21284082326 + 17749864125*z3) + (673066160 - 615195000*z3)*pow2(Mt) +
944  7350*(-520877 + 430155*z3)*pow2(Mst1)*pow2(s2t))*pow4(Msq)) + pow3(
945  Dmst12)*(-31305120*Mst1*s2t*pow2(Dmsqst1)*pow2(Mt) - 31305120*Dmsqst1*
946  Mst1*s2t*pow2(Msq)*pow2(Mt) + (2*Mst1*s2t*(31902674758 - 26534253375*
947  z3)*pow2(Mt) - 14700*Mt*(-520877 + 430155*z3)*pow2(Mst1)*pow2(s2t) +
948  40*(-16826654 + 15379875*z3)*pow3(Mt) + 245*(14217821 - 11852775*z3)*
949  pow3(Mst1)*pow3(s2t))*pow4(Msq)) + 4*Dmst12*pow2(Mt)*(-7826280*Mst1*
950  s2t*pow2(Dmsqst1) - 7826280*Dmsqst1*Mst1*s2t*pow2(Msq) + (10*Mt*(-
951  16826654 + 15379875*z3) + 49*Mst1*s2t*(-108352984 + 89636625*z3))*pow4(
952  Msq))*pow4(Mst2) - 392*pow3(Mt)*(154740*pow2(Dmsqst1) + 154740*Dmsqst1*
953  pow2(Msq) + (10583177 - 8913375*z3)*pow4(Msq))*pow6(Mst2)) + 49*pow3(
954  Mst1)*(-150*Dmsqst1*pow2(Msq)*(-4*Mt*pow2(Dmst12)*pow2(Mst2)*(19600*
955  Mst1*Mt*s2t + (-9122 + 14175*z3)*pow2(Mt) + 450*(2 - 21*z3)*pow2(Mst1)*
956  pow2(s2t)) + pow3(Dmst12)*(-20800*Mst1*s2t*pow2(Mt) - 75*Mt*(274 + 63*
957  z3)*pow2(Mst1)*pow2(s2t) + 16*(-6536 + 14175*z3)*pow3(Mt) - 8400*pow3(
958  Mst1)*pow3(s2t)) - 200*Dmst12*(-888*Mst1*s2t + Mt*(-158 + 567*z3))*
959  pow2(Mt)*pow4(Mst2) - 201600*pow3(Mt)*pow6(Mst2)) - 150*pow2(Dmsqst1)*(
960  -(Mt*pow2(Dmst12)*pow2(Mst2)*(36800*Mst1*Mt*s2t + 20112*pow2(Mt) + 75*(
961  458 - 945*z3)*pow2(Mst1)*pow2(s2t))) + 3*pow3(Dmst12)*(-20800*Mst1*s2t*
962  pow2(Mt) - 525*Mt*(-26 + 45*z3)*pow2(Mst1)*pow2(s2t) + 4*(-3998 +
963  14175*z3)*pow3(Mt) - 2800*pow3(Mst1)*pow3(s2t)) - 100*Dmst12*(-1360*
964  Mst1*s2t + 63*Mt*(-14 + 27*z3))*pow2(Mt)*pow4(Mst2) - 200*(410 + 567*
965  z3)*pow3(Mt)*pow6(Mst2)) - pow4(Msq)*(-20*Mt*pow2(Dmst12)*pow2(Mst2)*(
966  300*Mst1*Mt*s2t*(-17512 + 14805*z3) + (-375892 + 621675*z3)*pow2(Mt) -
967  900*(-1226 + 495*z3)*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-30*Mst1*
968  s2t*(-1117238 + 877575*z3)*pow2(Mt) - 2250*Mt*(-5570 + 333*z3)*pow2(
969  Mst1)*pow2(s2t) + (-41715182 + 38174625*z3)*pow3(Mt) + 3000*(-2722 +
970  2259*z3)*pow3(Mst1)*pow3(s2t)) - 6000*Dmst12*(81*Mt*(-406 + 285*z3) +
971  8*Mst1*s2t*(-694 + 477*z3))*pow2(Mt)*pow4(Mst2) + 48000*(623 + 963*z3)*
972  pow3(Mt)*pow6(Mst2))) + 196*Dmglst1*pow2(Mst1)*(600*pow2(Dmsqst1)*(-6*
973  Mt*pow2(Dmst12)*pow2(Mst2)*(-347*Mst1*Mt*s2t + 2538*pow2(Mt) + 1175*
974  pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(386*Mst1*s2t*pow2(Mt) + 14100*Mt*
975  pow2(Mst1)*pow2(s2t) + 15228*pow3(Mt) + 125*pow3(Mst1)*pow3(s2t)) + 2*
976  Dmst12*(7614*Mt - 2275*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 22600*pow3(Mt)*
977  pow6(Mst2)) + 600*Dmsqst1*pow2(Msq)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*(-
978  1516*Mst1*Mt*s2t + 7614*pow2(Mt) + 3525*pow2(Mst1)*pow2(s2t)) + pow3(
979  Dmst12)*(-564*Mst1*s2t*pow2(Mt) + 14100*Mt*pow2(Mst1)*pow2(s2t) +
980  15228*pow3(Mt) + 125*pow3(Mst1)*pow3(s2t)) + 4*Dmst12*(3807*Mt - 1375*
981  Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 29200*pow3(Mt)*pow6(Mst2)) + pow4(Msq)*
982  (Mt*pow2(Dmst12)*pow2(Mst2)*(40*Mst1*Mt*s2t*(-4511549 + 3729375*z3) + (
983  45149198 - 35285625*z3)*pow2(Mt) - 47250*(-430 + 207*z3)*pow2(Mst1)*
984  pow2(s2t)) + pow3(Dmst12)*(2*Mst1*s2t*(-28188929 + 23099625*z3)*pow2(
985  Mt) + 225*Mt*(-160136 + 100785*z3)*pow2(Mst1)*pow2(s2t) + (115176444 -
986  85920750*z3)*pow3(Mt) + 1125*(22174 - 18791*z3)*pow3(Mst1)*pow3(s2t)) +
987  40*Dmst12*(150*Mst1*s2t*(-19856 + 17667*z3) + Mt*(-5136871 + 3912300*
988  z3))*pow2(Mt)*pow4(Mst2) + 6000*(-31142 + 22653*z3)*pow3(Mt)*pow6(Mst2)
989  )) + 2*Mst1*pow2(Dmglst1)*(35280*Dmsqst1*Mt*pow2(Msq)*(-(pow2(Dmst12)*
990  pow2(Mst2)*(7940*Mst1*Mt*s2t + 914*pow2(Mt) + 225*pow2(Mst1)*pow2(s2t))
991  ) + (7940*Mst1*Mt*s2t + 914*pow2(Mt) + 450*pow2(Mst1)*pow2(s2t))*pow3(
992  Dmst12) + 2*Dmst12*Mt*(457*Mt + 3970*Mst1*s2t)*pow4(Mst2) + 4760*pow2(
993  Mt)*pow6(Mst2)) + 35280*Mt*pow2(Dmsqst1)*(-(pow2(Dmst12)*pow2(Mst2)*(
994  7940*Mst1*Mt*s2t + 914*pow2(Mt) + 225*pow2(Mst1)*pow2(s2t))) + (7940*
995  Mst1*Mt*s2t + 914*pow2(Mt) + 450*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) +
996  2*Dmst12*Mt*(457*Mt + 3970*Mst1*s2t)*pow4(Mst2) + 5022*pow2(Mt)*pow6(
997  Mst2)) + pow4(Msq)*(Mt*pow2(Dmst12)*pow2(Mst2)*(196*Mst1*Mt*s2t*(-
998  263717842 + 218365875*z3) + (27129768542 - 22522586625*z3)*pow2(Mt) +
999  22050*(-63802 + 92895*z3)*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-392*
1000  Mst1*s2t*(-384557822 + 319453875*z3)*pow2(Mt) + 66150*Mt*(-150737 +
1001  112420*z3)*pow2(Mst1)*pow2(s2t) + (-25287306692 + 22556819250*z3)*pow3(
1002  Mt) + 3675*(1908362 - 1581075*z3)*pow3(Mst1)*pow3(s2t)) + 392*Dmst12*(
1003  20*Mst1*s2t*(-6041999 + 5054400*z3) + Mt*(-73908751 + 57368250*z3))*
1004  pow2(Mt)*pow4(Mst2) + 3920*(-8223692 + 6125625*z3)*pow3(Mt)*pow6(Mst2))
1005  ))/(1.9845e6*pow3(Mst1)*pow3(Mt)*pow4(Msq)*pow6(Mst2));
1006 
1007  return result;
1008 }
1009 
1010 /**
1011  * @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)
1012  */
1013 double H3q22g::calc_coef_at_as2_no_sm_logs_log1() const {
1014 
1015  const double result =
1016  (4*(2*pow3(Dmglst1)*(-2*Mt*pow2(Dmst12)*pow2(Mst2)*(82320*Mst1*Mt*s2t*
1017  pow2(Dmsqst1) + 82320*Dmsqst1*Mst1*Mt*s2t*pow2(Msq) + (346639*Mst1*Mt*
1018  s2t + 555463*pow2(Mt) + 1051785*pow2(Mst1)*pow2(s2t))*pow4(Msq)) +
1019  pow3(Dmst12)*(164640*Mst1*s2t*pow2(Dmsqst1)*pow2(Mt) + 164640*Dmsqst1*
1020  Mst1*s2t*pow2(Msq)*pow2(Mt) + (8778304*Mst1*s2t*pow2(Mt) + 4207140*Mt*
1021  pow2(Mst1)*pow2(s2t) + 1110926*pow3(Mt) + 661255*pow3(Mst1)*pow3(s2t))*
1022  pow4(Msq)) + 2*Dmst12*pow2(Mt)*(82320*Mst1*s2t*pow2(Dmsqst1) + 82320*
1023  Dmsqst1*Mst1*s2t*pow2(Msq) + (555463*Mt - 3695874*Mst1*s2t)*pow4(Msq))*
1024  pow4(Mst2) - 4704*pow3(Mt)*(10*pow2(Dmsqst1) + 10*Dmsqst1*pow2(Msq) +
1025  141*pow4(Msq))*pow6(Mst2)) + Mst1*pow2(Dmglst1)*(17640*pow2(Dmsqst1)*
1026  pow2(Mt)*((-17*Mt + 10*Mst1*s2t)*pow2(Dmst12)*pow2(Mst2) + (17*Mt - 10*
1027  Mst1*s2t)*pow3(Dmst12) + Dmst12*(17*Mt - 10*Mst1*s2t)*pow4(Mst2) + 21*
1028  Mt*pow6(Mst2)) + 17640*Dmsqst1*pow2(Msq)*pow2(Mt)*((-17*Mt + 10*Mst1*
1029  s2t)*pow2(Dmst12)*pow2(Mst2) + (17*Mt - 10*Mst1*s2t)*pow3(Dmst12) +
1030  Dmst12*(17*Mt - 10*Mst1*s2t)*pow4(Mst2) + 35*Mt*pow6(Mst2)) + pow4(Msq)
1031  *(Mt*pow2(Dmst12)*pow2(Mst2)*(-4088560*Mst1*Mt*s2t + 968629*pow2(Mt) +
1032  1853670*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(7502488*Mst1*s2t*pow2(Mt)
1033  - 3134775*Mt*pow2(Mst1)*pow2(s2t) + 2019394*pow3(Mt) + 689430*pow3(
1034  Mst1)*pow3(s2t)) - 196*Dmst12*(20187*Mt - 3442*Mst1*s2t)*pow2(Mt)*pow4(
1035  Mst2) - 8199072*pow3(Mt)*pow6(Mst2))) - 98*Dmglst1*pow2(Mst1)*(600*
1036  Dmsqst1*pow2(Msq)*pow2(Mt)*(-((6*Mt + Mst1*s2t)*pow2(Dmst12)*pow2(Mst2)
1037  ) + (6*Mt - 3*Mst1*s2t)*pow3(Dmst12) + Dmst12*(6*Mt + 5*Mst1*s2t)*pow4(
1038  Mst2) + 25*Mt*pow6(Mst2)) + 300*pow2(Dmsqst1)*pow2(Mt)*(3*(-4*Mt +
1039  Mst1*s2t)*pow2(Dmst12)*pow2(Mst2) + (12*Mt - 11*Mst1*s2t)*pow3(Dmst12)
1040  + Dmst12*(12*Mt + 5*Mst1*s2t)*pow4(Mst2) + 35*Mt*pow6(Mst2)) + pow4(
1041  Msq)*(Mt*pow2(Dmst12)*pow2(Mst2)*(29758*Mst1*Mt*s2t + 6677*pow2(Mt) +
1042  22350*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(-34587*Mst1*s2t*pow2(Mt) -
1043  18045*Mt*pow2(Mst1)*pow2(s2t) + 22414*pow3(Mt) + 3325*pow3(Mst1)*pow3(
1044  s2t)) - 8*Dmst12*(4471*Mt + 6875*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 50800*
1045  pow3(Mt)*pow6(Mst2))) - 49*pow3(Mst1)*(-150*Dmsqst1*Mt*pow2(Msq)*(pow2(
1046  Dmst12)*pow2(Mst2)*(-80*Mst1*Mt*s2t - 17*pow2(Mt) + 180*pow2(Mst1)*
1047  pow2(s2t)) + 2*(20*Mst1*Mt*s2t + 187*pow2(Mt) - 90*pow2(Mst1)*pow2(s2t)
1048  )*pow3(Dmst12) - 20*Dmst12*Mt*(17*Mt - 6*Mst1*s2t)*pow4(Mst2) - 1080*
1049  pow2(Mt)*pow6(Mst2)) - 150*Mt*pow2(Dmsqst1)*(-4*pow2(Dmst12)*pow2(Mst2)
1050  *(10*Mst1*Mt*s2t + 3*pow2(Mt) - 45*pow2(Mst1)*pow2(s2t)) + 9*(41*pow2(
1051  Mt) - 20*pow2(Mst1)*pow2(s2t))*pow3(Dmst12) + 5*Dmst12*Mt*(-69*Mt + 16*
1052  Mst1*s2t)*pow4(Mst2) - 1010*pow2(Mt)*pow6(Mst2)) - pow4(Msq)*(-2*Mt*
1053  pow2(Dmst12)*pow2(Mst2)*(25850*Mst1*Mt*s2t + 21033*pow2(Mt) + 75*pow2(
1054  Mst1)*pow2(s2t)) + pow3(Dmst12)*(68264*Mst1*s2t*pow2(Mt) + 5700*Mt*
1055  pow2(Mst1)*pow2(s2t) + 36107*pow3(Mt) + 250*pow3(Mst1)*pow3(s2t)) +
1056  200*Dmst12*(131*Mt + 100*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 400*(-533 +
1057  54*z3)*pow3(Mt)*pow6(Mst2)))))/(33075.*pow3(Mst1)*pow3(Mt)*pow4(Msq)*
1058  pow6(Mst2));
1059 
1060  return result;
1061 }
1062 
1063 /**
1064  * @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)
1065  */
1066 double H3q22g::calc_coef_at_as2_no_sm_logs_log2() const {
1067 
1068  const double result =
1069  (8*(-7*Dmglst1*pow2(Mst1)*pow4(Msq)*(2*Mt*pow2(Dmst12)*pow2(Mst2)*(-1304*
1070  Mst1*Mt*s2t + 888*pow2(Mt) + 645*pow2(Mst1)*pow2(s2t)) + pow3(Dmst12)*(
1071  1544*Mst1*s2t*pow2(Mt) - 2250*Mt*pow2(Mst1)*pow2(s2t) - 1640*pow3(Mt) +
1072  275*pow3(Mst1)*pow3(s2t)) - 8*Dmst12*(239*Mt - 35*Mst1*s2t)*pow2(Mt)*
1073  pow4(Mst2) - 6520*pow3(Mt)*pow6(Mst2)) - 8*pow3(Dmglst1)*pow4(Msq)*(-3*
1074  Mt*pow2(Dmst12)*pow2(Mst2)*(-75*Mst1*Mt*s2t + 1122*pow2(Mt) + 560*pow2(
1075  Mst1)*pow2(s2t)) + pow3(Dmst12)*(1671*Mst1*s2t*pow2(Mt) + 3360*Mt*pow2(
1076  Mst1)*pow2(s2t) + 3366*pow3(Mt) + 140*pow3(Mst1)*pow3(s2t)) + 3*Dmst12*
1077  (1122*Mt - 707*Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 2163*pow3(Mt)*pow6(Mst2)
1078  ) + Mst1*pow2(Dmglst1)*pow4(Msq)*(Mt*pow2(Dmst12)*pow2(Mst2)*(4088*
1079  Mst1*Mt*s2t - 22884*pow2(Mt) + 8925*pow2(Mst1)*pow2(s2t)) + pow3(
1080  Dmst12)*(-42728*Mst1*s2t*pow2(Mt) + 1155*Mt*pow2(Mst1)*pow2(s2t) +
1081  56772*pow3(Mt) - 3360*pow3(Mst1)*pow3(s2t)) - 28*Dmst12*(393*Mt - 1234*
1082  Mst1*s2t)*pow2(Mt)*pow4(Mst2) + 3948*pow3(Mt)*pow6(Mst2)) + 7*pow3(
1083  Mst1)*(Mt*pow2(Dmst12)*pow2(Mst2)*(1760*Mst1*Mt*s2t + 538*pow2(Mt) +
1084  105*pow2(Mst1)*pow2(s2t))*pow4(Msq) - pow3(Dmst12)*(1032*Mst1*s2t*pow2(
1085  Mt) + 480*Mt*pow2(Mst1)*pow2(s2t) + 644*pow3(Mt) - 45*pow3(Mst1)*pow3(
1086  s2t))*pow4(Msq) - 40*Dmst12*(11*Mt + 61*Mst1*s2t)*pow2(Mt)*pow4(Msq)*
1087  pow4(Mst2) + 10*pow3(Mt)*(45*pow2(Dmsqst1) + 90*Dmsqst1*pow2(Msq) +
1088  442*pow4(Msq))*pow6(Mst2))))/(315.*pow3(Mst1)*pow3(Mt)*pow4(Msq)*pow6(
1089  Mst2));
1090 
1091  return result;
1092 }
1093 
1094 /**
1095  * @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)
1096  */
1097 double H3q22g::calc_coef_at_as2_no_sm_logs_log3() const {
1098 
1099  const double result =
1100  -298.6666666666667;
1101 
1102  return result;
1103 }
1104 
1105 } // namespace hierarchies
1106 } // namespace himalaya
himalaya::hierarchies::H3q22g::calc_coef_at_as2_no_sm_logs_log3
double calc_coef_at_as2_no_sm_logs_log3() const
Definition: H3q22g.cpp:1097
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
Definition: H3.cpp:14
himalaya::hierarchies::H3q22g::xDmsqst1
int xDmsqst1
MDR and truncation flags.
Definition: H3q22g.hpp:73
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::H3q22g::lmMt
double lmMt
common variables
Definition: H3q22g.hpp:72
himalaya::hierarchies::ExpansionDepth::Msq
truncate the expansion depth in the average squark mass by one order
Definition: HierarchyFlags.hpp:48
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::H3q22g::calc_coef_at_as2_no_sm_logs_log2
double calc_coef_at_as2_no_sm_logs_log2() const
Definition: H3q22g.cpp:1066
himalaya::hierarchies::H3q22g::H3q22g
H3q22g(const ExpansionFlags_t &expansionDepth, double Al4p, double beta, double Dmglst1, double Dmst12, double Dmsqst1, double lmMt, double lmMst1, double Mt, double Mst1, double Mst2, double Msq, double MuSUSY, double s2t, int mdrFlag, int oneLoopFlag, int twoLoopFlag, int threeLoopFlag)
Definition: H3q22g.cpp:38
himalaya::hierarchies::H3q22g::threeLoopFlag
int threeLoopFlag
loop flags
Definition: H3q22g.hpp:74
himalaya::hierarchies::H3q22g::calc_coef_at_as2_no_sm_logs_log1
double calc_coef_at_as2_no_sm_logs_log1() const
Definition: H3q22g.cpp:1013
himalaya::hierarchies::H3q22g::calc_coef_at_as2_no_sm_logs_log0
double calc_coef_at_as2_no_sm_logs_log0() const
Definition: H3q22g.cpp:938
Generated by   doxygen 1.8.13