MSSMatNLO: use-MoGRe.nb

File use-MoGRe.nb, 47.2 KB (added by Benjamin Fuks, 7 years ago)

Illustrative Mathematica motebook using MoGRe

Line 
1(* Content-type: application/mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 7.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 145, 7]
13NotebookDataLength[ 48186, 1232]
14NotebookOptionsPosition[ 45995, 1158]
15NotebookOutlinePosition[ 46372, 1175]
16CellTagsIndexPosition[ 46329, 1172]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"Quit", "[", "]"}]], "Input",
23 CellChangeTimes->{{3.642999466487681*^9, 3.6429994678496943`*^9}, {
24 3.6477455431974373`*^9, 3.647745544864357*^9}, {3.648269677005701*^9,
25 3.648269677954204*^9}, {3.744954871093508*^9, 3.7449548719161787`*^9}, {
26 3.74584274317243*^9, 3.74584274398195*^9}, {3.745857688966254*^9,
27 3.7458576899132833`*^9}, {3.757993454379569*^9, 3.757993455232729*^9}, {
28 3.758250980508205*^9, 3.7582509813882504`*^9}}],
29
30Cell[CellGroupData[{
31
32Cell["Setup", "Subsection",
33 CellChangeTimes->{{3.518892288721245*^9, 3.518892291440427*^9}, {
34 3.648272586359296*^9, 3.648272589589151*^9}, {3.744951163067635*^9,
35 3.744951164882929*^9}}],
36
37Cell[BoxData[{
38 RowBox[{
39 RowBox[{"$FeynRulesPath", "=",
40 RowBox[{"SetDirectory", "[", "\"\<<FeynRules path\>\"", "]"}]}],
41 ";"}], "\[IndentingNewLine]",
42 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
43 RowBox[{
44 RowBox[{"LoadModel", "[", "\"\<<Path to mssm-nlo.fr>\>\"", "]"}],
45 ";"}], "\[IndentingNewLine]",
46 RowBox[{
47 RowBox[{"LoadRestriction", "[", "\"\<<Path to Massless.rst>\>\"", "]"}],
48 ";"}], "\[IndentingNewLine]",
49 RowBox[{
50 RowBox[{"LoadRestriction", "[", "\"\<<Path to DiagonalCKM.rst\>\"", "]"}],
51 ";"}], "\[IndentingNewLine]",
52 RowBox[{
53 RowBox[{
54 RowBox[{"SetDirectory", "[",
55 RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], " ",
56 RowBox[{"(*",
57 RowBox[{"Where", " ", "MoGRe", " ", "is", " ", "located"}], " ",
58 "*)"}]}], "\[IndentingNewLine]",
59 RowBox[{
60 RowBox[{"Begin", "[", "\"\<MoGRe`\>\"", "]"}], ";"}], "\[IndentingNewLine]",
61 RowBox[{
62 RowBox[{"<<", "MoGRe.m"}], ";"}], "\[IndentingNewLine]",
63 RowBox[{
64 RowBox[{"End", "[", "]"}], ";"}]}], "Input",
65 CellChangeTimes->{{3.500793438723139*^9, 3.500793441143552*^9}, {
66 3.503666867066519*^9, 3.503666869077112*^9}, {3.518892141094482*^9,
67 3.518892168024193*^9}, {3.57172073024898*^9, 3.5717207668526573`*^9}, {
68 3.623126068798016*^9, 3.623126077634707*^9}, {3.6429991439343348`*^9,
69 3.642999184233075*^9}, {3.642999745375638*^9, 3.642999772385662*^9}, {
70 3.6431183608038282`*^9, 3.6431183609567623`*^9}, {3.643434363202712*^9,
71 3.6434343707164593`*^9}, {3.643452164235368*^9, 3.643452168865634*^9}, {
72 3.643452224521432*^9, 3.643452265717387*^9}, {3.643453069623658*^9,
73 3.6434530698976583`*^9}, {3.643453150320168*^9, 3.643453155861721*^9}, {
74 3.643454581300136*^9, 3.6434545849965267`*^9}, 3.6434805724133167`*^9, {
75 3.643483032880064*^9, 3.643483037703815*^9}, {3.643505734766756*^9,
76 3.64350573544672*^9}, 3.643510615279973*^9, {3.644042120992399*^9,
77 3.64404213031954*^9}, {3.644042232974127*^9, 3.644042235500444*^9}, {
78 3.644053058481842*^9, 3.6440530586177*^9}, {3.6441215850229073`*^9,
79 3.644121586348798*^9}, {3.6463093381530237`*^9, 3.646309341741721*^9},
80 3.646309402062256*^9, {3.648271608223464*^9, 3.648271614048018*^9}, {
81 3.6482717399947*^9, 3.648271742624647*^9}, {3.65271887449249*^9,
82 3.652718881699897*^9}, {3.652758115172617*^9, 3.652758138006729*^9}, {
83 3.65278445775749*^9, 3.652784457934544*^9}, {3.687700274467225*^9,
84 3.687700276625421*^9}, {3.687700349669055*^9, 3.687700365740831*^9}, {
85 3.6877005443628674`*^9, 3.6877005473043213`*^9}, {3.687700830581005*^9,
86 3.6877008370565042`*^9}, {3.688211032762238*^9, 3.688211033010005*^9}, {
87 3.688211128176875*^9, 3.688211128964767*^9}, {3.688211218473372*^9,
88 3.688211258461564*^9}, 3.688211613614009*^9, {3.688234821681107*^9,
89 3.688234828127536*^9}, {3.688236820715692*^9, 3.688236820899551*^9}, {
90 3.688914299821698*^9, 3.688914305788663*^9}, {3.688915189003973*^9,
91 3.688915189162331*^9}, {3.688930995258505*^9, 3.68893100017045*^9}, {
92 3.68898319860035*^9, 3.6889831987631693`*^9}, {3.688984936179865*^9,
93 3.6889849376475973`*^9}, {3.689307114415225*^9, 3.689307116753755*^9}, {
94 3.7136993130469503`*^9, 3.713699321159403*^9}, {3.7137204429773073`*^9,
95 3.713720445250545*^9}, {3.745153119299178*^9, 3.745153119658555*^9}, {
96 3.758247811625256*^9, 3.758247827461388*^9}, {3.75825834893437*^9,
97 3.758258349643984*^9}, {3.771263107602776*^9, 3.771263121935561*^9}, {
98 3.7715515410537024`*^9, 3.771551604765225*^9}}]
99}, Open ]],
100
101Cell[CellGroupData[{
102
103Cell["Computation of the Lagrangian", "Subsection",
104 CellChangeTimes->{{3.652758944615611*^9, 3.652758945925817*^9}, {
105 3.652788971327319*^9, 3.652788972429942*^9}, {3.717927237042218*^9,
106 3.7179272380337563`*^9}, {3.771263145046253*^9, 3.771263159779201*^9}}],
107
108Cell[TextData[{
109 " - The first time the notebook is run, the Lagrangian is computed. This can \
110take some time (pease be patient).\n - The results are store in a file \
111(called ",
112 StyleBox["lagrangian.dat",
113 FontSlant->"Italic"],
114 ") so that this part can be skiped in the future (the file with the \
115Lagrangian will instead be loaded).\n - After calculation or loading of the \
116file, the Lagrangian is stored into a variabl called \"lagr\""
117}], "Text",
118 CellChangeTimes->{{3.771263163734034*^9, 3.771263178118503*^9}, {
119 3.771263277883705*^9, 3.7712633902566853`*^9}}],
120
121Cell[BoxData[{
122 RowBox[{
123 RowBox[{"SetDirectory", "[",
124 RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
125 ";"}], "\[IndentingNewLine]",
126 RowBox[{
127 RowBox[{"If", "[",
128 RowBox[{
129 RowBox[{"!",
130 RowBox[{"FileExistsQ", "[", "\"\<lagrangian.dat\>\"", "]"}]}], ",",
131 "\[IndentingNewLine]",
132 RowBox[{
133 RowBox[{"lagr", "=", "LSUSY"}], ";", "\[IndentingNewLine]",
134 RowBox[{
135 RowBox[{"Definition", "[", "lagr", "]"}], " ", ">>", " ",
136 "lagrangian.dat"}]}], ",", "\[IndentingNewLine]",
137 RowBox[{
138 RowBox[{"<<", "lagrangian.dat"}], ";"}]}], "\[IndentingNewLine]", "]"}],
139 ";"}]}], "Input",
140 CellChangeTimes->{{3.65277982379673*^9, 3.6527798242214746`*^9}, {
141 3.65278011347823*^9, 3.652780125570778*^9}, {3.652780822315733*^9,
142 3.6527808271772423`*^9}, 3.652780876249255*^9, {3.65278103371167*^9,
143 3.652781036007701*^9}, {3.652784471052888*^9, 3.6527844714268417`*^9}, {
144 3.771263259581438*^9, 3.771263273327767*^9}, {3.7712634003294086`*^9,
145 3.771263400637439*^9}}]
146}, Open ]],
147
148Cell[CellGroupData[{
149
150Cell["\<\
151Using the MoGRe package to fix the renormalisation conditions\
152\>", "Subsection",
153 CellChangeTimes->{{3.744954096662758*^9, 3.744954097340728*^9}, {
154 3.745151786089308*^9, 3.745151794439547*^9}, {3.745151839727529*^9,
155 3.745151841215969*^9}, {3.745663645010592*^9, 3.7456636463056517`*^9}, {
156 3.745687290761694*^9, 3.745687291593329*^9}, {3.745841616116417*^9,
157 3.745841617247401*^9}, {3.74584378174969*^9, 3.745843800667904*^9}, {
158 3.74589615607681*^9, 3.74589617695035*^9}, 3.757417549664226*^9, {
159 3.7579096806644897`*^9, 3.7579096841682577`*^9}, {3.758159046212031*^9,
160 3.758159082144144*^9}, {3.758168073340328*^9, 3.7581680742810583`*^9}, {
161 3.758250674754568*^9, 3.758250677493423*^9}, {3.771263448171097*^9,
162 3.771263462210712*^9}, {3.771264763367551*^9, 3.771264765778284*^9}}],
163
164Cell[CellGroupData[{
165
166Cell["Setup", "Subsubsection",
167 CellChangeTimes->{{3.7712636194863377`*^9, 3.771263622434663*^9}}],
168
169Cell[BoxData[
170 RowBox[{
171 RowBox[{"(*", " ",
172 RowBox[{
173 RowBox[{"To", " ", "get", " ", "more", " ", "print"}], "-",
174 RowBox[{"outs", " ", "to", " ", "the", " ", "screen"}]}], " ", "*)"}],
175 "\[IndentingNewLine]",
176 RowBox[{
177 RowBox[{"MoGRe$Debug", "=", "True"}], ";"}]}]], "Input",
178 CellChangeTimes->{{3.771263627141055*^9, 3.7712636335707293`*^9}}],
179
180Cell[BoxData[
181 RowBox[{
182 RowBox[{"(*", " ",
183 RowBox[{
184 "Setting", " ", "the", " ", "options", " ", "of", " ", "the", " ", "main",
185 " ", "function"}], " ", "*)"}], "\[IndentingNewLine]",
186 RowBox[{
187 RowBox[{"SetOptions", "[",
188 RowBox[{"MoGRe$Renormalize", ",", " ",
189 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}], ",",
190 RowBox[{"FlavorMixing", "\[Rule]",
191 RowBox[{"{",
192 RowBox[{
193 RowBox[{"{",
194 RowBox[{"st1", ",", "st2"}], "}"}], ",",
195 RowBox[{"{",
196 RowBox[{"sb1", ",", "sb2"}], "}"}]}], "}"}]}], ",",
197 RowBox[{"CouplingOrders", "\[Rule]",
198 RowBox[{"{", "QCD", "}"}]}]}], "]"}], ";"}]}]], "Input",
199 CellChangeTimes->{{3.771263641132224*^9, 3.7712636559080553`*^9}}]
200}, Open ]],
201
202Cell[CellGroupData[{
203
204Cell["Naive on-shell (OS) scheme", "Subsubsection",
205 CellChangeTimes->{{3.771264549679812*^9, 3.771264557709045*^9}, {
206 3.771264725891897*^9, 3.771264729628981*^9}}],
207
208Cell[BoxData[
209 RowBox[{
210 RowBox[{"(*", " ",
211 RowBox[{"Unrenormalized", " ", "parameters"}], " ", "*)"}],
212 "\[IndentingNewLine]",
213 RowBox[{
214 RowBox[{"MoGRe`DefineUnrenormalizedParameters", "[",
215 RowBox[{"{", "\[IndentingNewLine]", " ",
216 RowBox[{
217 "Gf", ",", "aEWM1", ",", "MZ", ",", " ", "MUH", ",", "alp", ",", "tb",
218 ",", "\[IndentingNewLine]", "Mx1", ",", "Mx2", ",", "mHu2", ",", "mHd2",
219 ",", "meL", ",", "mmuL", ",", "mtauL", ",", "meR", ",", "mmuR", ",",
220 "mtauR", ",", "\[IndentingNewLine]",
221 RowBox[{"Sequence", "@@",
222 RowBox[{"Flatten", "[",
223 RowBox[{"Table", "[",
224 RowBox[{
225 RowBox[{"{",
226 RowBox[{"ae", "[",
227 RowBox[{"i", ",", "i"}], "]"}], "}"}], ",",
228 RowBox[{"{",
229 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], ",",
230 "\[IndentingNewLine]",
231 RowBox[{"Sequence", "@@",
232 RowBox[{"Flatten", "[",
233 RowBox[{"Table", "[",
234 RowBox[{
235 RowBox[{"{",
236 RowBox[{
237 RowBox[{"au", "[",
238 RowBox[{"i", ",", "i"}], "]"}], ",",
239 RowBox[{"ad", "[",
240 RowBox[{"i", ",", "i"}], "]"}]}], "}"}], ",",
241 RowBox[{"{",
242 RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}], "]"}]}], ",",
243 "\[IndentingNewLine]",
244 RowBox[{"Sequence", "@@",
245 RowBox[{"Flatten", "[",
246 RowBox[{"Table", "[",
247 RowBox[{
248 RowBox[{"{",
249 RowBox[{
250 RowBox[{"VV", "[",
251 RowBox[{"i", ",", "j"}], "]"}], ",",
252 RowBox[{"UU", "[",
253 RowBox[{"i", ",", "j"}], "]"}]}], "}"}], ",",
254 RowBox[{"{",
255 RowBox[{"i", ",", "1", ",", "2"}], "}"}], ",",
256 RowBox[{"{",
257 RowBox[{"j", ",", "1", ",", "2"}], "}"}]}], "]"}], "]"}]}], ",",
258 "\[IndentingNewLine]",
259 RowBox[{"Sequence", "@@",
260 RowBox[{"Flatten", "[",
261 RowBox[{"Table", "[",
262 RowBox[{
263 RowBox[{"{",
264 RowBox[{"NN", "[",
265 RowBox[{"i", ",", "j"}], "]"}], "}"}], ",",
266 RowBox[{"{",
267 RowBox[{"i", ",", "1", ",", "4"}], "}"}], ",",
268 RowBox[{"{",
269 RowBox[{"j", ",", "1", ",", "4"}], "}"}]}], "]"}], "]"}]}]}],
270 "\[IndentingNewLine]", "}"}], "]"}], ";"}]}]], "Input",
271 CellChangeTimes->{{3.771264570309547*^9, 3.771264574380278*^9}}],
272
273Cell[BoxData[
274 RowBox[{
275 RowBox[{"(*", " ", "Fields", " ", "*)"}], "\[IndentingNewLine]",
276 RowBox[{
277 RowBox[{
278 RowBox[{"MoGRe`DeclareUnrenormalizedFields", "[",
279 RowBox[{"W", ",", "A", ",", "Z"}], "]"}], ";"}], "\[IndentingNewLine]",
280 RowBox[{
281 RowBox[{"MoGRe`RealFieldRenormalization", "[", "]"}], " ",
282 ";"}]}]}]], "Input",
283 CellChangeTimes->{{3.771264582333989*^9, 3.771264608943572*^9}}],
284
285Cell[BoxData[
286 RowBox[{
287 RowBox[{"(*", " ",
288 RowBox[{"Removing", " ", "some", " ", "parameters", " ",
289 RowBox[{"(",
290 RowBox[{
291 RowBox[{"cf", ".", " ", "the"}], " ", "way", " ", "NLOCT", " ",
292 "works"}], ")"}]}], " ", "*)"}], "\[IndentingNewLine]",
293 RowBox[{
294 RowBox[{
295 RowBox[{
296 RowBox[{"MoGRe`RemovingInternalCst", "[", "#", "]"}], "&"}], "/@",
297 RowBox[{"{", "\[IndentingNewLine]",
298 RowBox[{"gs", ",", "\[IndentingNewLine]",
299 RowBox[{"Sequence", "@@",
300 RowBox[{"Flatten", "[",
301 RowBox[{"Table", "[",
302 RowBox[{
303 RowBox[{"{",
304 RowBox[{
305 RowBox[{"yu", "[",
306 RowBox[{"i", ",", "i"}], "]"}], ",",
307 RowBox[{"yd", "[",
308 RowBox[{"i", ",", "i"}], "]"}], ",",
309 RowBox[{"tu", "[",
310 RowBox[{"i", ",", "i"}], "]"}], ",",
311 RowBox[{"td", "[",
312 RowBox[{"i", ",", "i"}], "]"}]}], "}"}], ",",
313 RowBox[{"{",
314 RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], ",",
315 "\[IndentingNewLine]",
316 RowBox[{"Sequence", "@@",
317 RowBox[{"Flatten", "[",
318 RowBox[{"Table", "[",
319 RowBox[{
320 RowBox[{"{",
321 RowBox[{"CKM", "[",
322 RowBox[{"i", ",", "j"}], "]"}], "}"}], ",",
323 RowBox[{"{",
324 RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",",
325 RowBox[{"{",
326 RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}]}],
327 "\[IndentingNewLine]", "}"}]}], ";"}]}]], "Input",
328 CellChangeTimes->{{3.771264625204327*^9, 3.771264640138379*^9}}],
329
330Cell[BoxData[
331 RowBox[{
332 RowBox[{"(*", " ",
333 RowBox[{
334 RowBox[{
335 "Implementing", " ", "the", " ", "naive", " ", "OS", " ", "scheme"}], " ",
336 "-", " ",
337 RowBox[{
338 RowBox[{"A", ".", " ", "stop"}], " ", "mixings"}]}], " ", "*)"}],
339 "\[IndentingNewLine]",
340 RowBox[{
341 RowBox[{
342 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
343 RowBox[{
344 RowBox[{"FR$delta", "[",
345 RowBox[{
346 RowBox[{"{",
347 RowBox[{"Rtop", "[",
348 RowBox[{"1", ",", "1"}], "]"}], "}"}], ",",
349 RowBox[{"{", "}"}]}], "]"}], ",",
350 RowBox[{
351 RowBox[{"1", "/", "4"}], " ",
352 RowBox[{"Rtop", "[",
353 RowBox[{"2", ",", "1"}], "]"}], "*",
354 RowBox[{"(",
355 RowBox[{
356 RowBox[{"FR$deltaZ", "[",
357 RowBox[{
358 RowBox[{"{",
359 RowBox[{"st1", ",", "st2"}], "}"}], ",",
360 RowBox[{"{",
361 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
362 RowBox[{"Conjugate", "[",
363 RowBox[{"FR$deltaZ", "[",
364 RowBox[{
365 RowBox[{"{",
366 RowBox[{"st2", ",", "st1"}], "}"}], ",",
367 RowBox[{"{",
368 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
369 ";"}], "\[IndentingNewLine]",
370 RowBox[{
371 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
372 RowBox[{
373 RowBox[{"FR$delta", "[",
374 RowBox[{
375 RowBox[{"{",
376 RowBox[{"Rtop", "[",
377 RowBox[{"1", ",", "2"}], "]"}], "}"}], ",",
378 RowBox[{"{", "}"}]}], "]"}], ",",
379 RowBox[{
380 RowBox[{"1", "/", "4"}], " ",
381 RowBox[{"Rtop", "[",
382 RowBox[{"1", ",", "1"}], "]"}], "*",
383 RowBox[{"(",
384 RowBox[{
385 RowBox[{"FR$deltaZ", "[",
386 RowBox[{
387 RowBox[{"{",
388 RowBox[{"st1", ",", "st2"}], "}"}], ",",
389 RowBox[{"{",
390 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
391 RowBox[{"Conjugate", "[",
392 RowBox[{"FR$deltaZ", "[",
393 RowBox[{
394 RowBox[{"{",
395 RowBox[{"st2", ",", "st1"}], "}"}], ",",
396 RowBox[{"{",
397 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
398 ";"}], "\[IndentingNewLine]",
399 RowBox[{
400 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
401 RowBox[{
402 RowBox[{"FR$delta", "[",
403 RowBox[{
404 RowBox[{"{",
405 RowBox[{"Rtop", "[",
406 RowBox[{"2", ",", "1"}], "]"}], "}"}], ",",
407 RowBox[{"{", "}"}]}], "]"}], ",",
408 RowBox[{
409 RowBox[{"1", "/", "4"}], " ",
410 RowBox[{"Rtop", "[",
411 RowBox[{"1", ",", "1"}], "]"}], "*",
412 RowBox[{"(",
413 RowBox[{
414 RowBox[{"FR$deltaZ", "[",
415 RowBox[{
416 RowBox[{"{",
417 RowBox[{"st1", ",", "st2"}], "}"}], ",",
418 RowBox[{"{",
419 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
420 RowBox[{"Conjugate", "[",
421 RowBox[{"FR$deltaZ", "[",
422 RowBox[{
423 RowBox[{"{",
424 RowBox[{"st2", ",", "st1"}], "}"}], ",",
425 RowBox[{"{",
426 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
427 ";"}], "\[IndentingNewLine]",
428 RowBox[{
429 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
430 RowBox[{
431 RowBox[{"FR$delta", "[",
432 RowBox[{
433 RowBox[{"{",
434 RowBox[{"Rtop", "[",
435 RowBox[{"2", ",", "2"}], "]"}], "}"}], ",",
436 RowBox[{"{", "}"}]}], "]"}], ",",
437 RowBox[{
438 RowBox[{"1", "/", "4"}], " ",
439 RowBox[{"Rtop", "[",
440 RowBox[{"2", ",", "1"}], "]"}], "*",
441 RowBox[{"(",
442 RowBox[{
443 RowBox[{"FR$deltaZ", "[",
444 RowBox[{
445 RowBox[{"{",
446 RowBox[{"st1", ",", "st2"}], "}"}], ",",
447 RowBox[{"{",
448 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
449 RowBox[{"Conjugate", "[",
450 RowBox[{"FR$deltaZ", "[",
451 RowBox[{
452 RowBox[{"{",
453 RowBox[{"st2", ",", "st1"}], "}"}], ",",
454 RowBox[{"{",
455 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
456 ";"}]}]}]], "Input",
457 CellChangeTimes->{{3.771264647562903*^9, 3.771264712956895*^9}}],
458
459Cell[BoxData[
460 RowBox[{
461 RowBox[{"(*", " ",
462 RowBox[{
463 RowBox[{
464 "Implementing", " ", "the", " ", "naive", " ", "OS", " ", "scheme"}], " ",
465 "-", " ",
466 RowBox[{
467 RowBox[{"B", ".", " ", "sbottom"}], " ", "mixings"}]}], " ", "*)"}],
468 "\[IndentingNewLine]",
469 RowBox[{
470 RowBox[{
471 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
472 RowBox[{
473 RowBox[{"FR$delta", "[",
474 RowBox[{
475 RowBox[{"{",
476 RowBox[{"Rbot", "[",
477 RowBox[{"1", ",", "1"}], "]"}], "}"}], ",",
478 RowBox[{"{", "}"}]}], "]"}], ",",
479 RowBox[{
480 RowBox[{"1", "/", "4"}], " ",
481 RowBox[{"Rbot", "[",
482 RowBox[{"2", ",", "1"}], "]"}], "*",
483 RowBox[{"(",
484 RowBox[{
485 RowBox[{"FR$deltaZ", "[",
486 RowBox[{
487 RowBox[{"{",
488 RowBox[{"sb1", ",", "sb2"}], "}"}], ",",
489 RowBox[{"{",
490 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
491 RowBox[{"Conjugate", "[",
492 RowBox[{"FR$deltaZ", "[",
493 RowBox[{
494 RowBox[{"{",
495 RowBox[{"sb2", ",", "sb1"}], "}"}], ",",
496 RowBox[{"{",
497 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
498 ";"}], "\[IndentingNewLine]",
499 RowBox[{
500 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
501 RowBox[{
502 RowBox[{"FR$delta", "[",
503 RowBox[{
504 RowBox[{"{",
505 RowBox[{"Rbot", "[",
506 RowBox[{"1", ",", "2"}], "]"}], "}"}], ",",
507 RowBox[{"{", "}"}]}], "]"}], ",",
508 RowBox[{
509 RowBox[{"1", "/", "4"}], " ",
510 RowBox[{"Rbot", "[",
511 RowBox[{"1", ",", "1"}], "]"}], "*",
512 RowBox[{"(",
513 RowBox[{
514 RowBox[{"FR$deltaZ", "[",
515 RowBox[{
516 RowBox[{"{",
517 RowBox[{"sb1", ",", "sb2"}], "}"}], ",",
518 RowBox[{"{",
519 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
520 RowBox[{"Conjugate", "[",
521 RowBox[{"FR$deltaZ", "[",
522 RowBox[{
523 RowBox[{"{",
524 RowBox[{"sb2", ",", "sb1"}], "}"}], ",",
525 RowBox[{"{",
526 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
527 ";"}], "\[IndentingNewLine]",
528 RowBox[{
529 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
530 RowBox[{
531 RowBox[{"FR$delta", "[",
532 RowBox[{
533 RowBox[{"{",
534 RowBox[{"Rbot", "[",
535 RowBox[{"2", ",", "1"}], "]"}], "}"}], ",",
536 RowBox[{"{", "}"}]}], "]"}], ",",
537 RowBox[{
538 RowBox[{"1", "/", "4"}], " ",
539 RowBox[{"Rbot", "[",
540 RowBox[{"1", ",", "1"}], "]"}], "*",
541 RowBox[{"(",
542 RowBox[{
543 RowBox[{"FR$deltaZ", "[",
544 RowBox[{
545 RowBox[{"{",
546 RowBox[{"sb1", ",", "sb2"}], "}"}], ",",
547 RowBox[{"{",
548 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
549 RowBox[{"Conjugate", "[",
550 RowBox[{"FR$deltaZ", "[",
551 RowBox[{
552 RowBox[{"{",
553 RowBox[{"sb2", ",", "sb1"}], "}"}], ",",
554 RowBox[{"{",
555 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
556 ";"}], "\[IndentingNewLine]",
557 RowBox[{
558 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
559 RowBox[{
560 RowBox[{"FR$delta", "[",
561 RowBox[{
562 RowBox[{"{",
563 RowBox[{"Rbot", "[",
564 RowBox[{"2", ",", "2"}], "]"}], "}"}], ",",
565 RowBox[{"{", "}"}]}], "]"}], ",",
566 RowBox[{
567 RowBox[{"1", "/", "4"}], " ",
568 RowBox[{"Rbot", "[",
569 RowBox[{"2", ",", "1"}], "]"}], "*",
570 RowBox[{"(",
571 RowBox[{
572 RowBox[{"FR$deltaZ", "[",
573 RowBox[{
574 RowBox[{"{",
575 RowBox[{"sb1", ",", "sb2"}], "}"}], ",",
576 RowBox[{"{",
577 RowBox[{"{", "}"}], "}"}]}], "]"}], "-",
578 RowBox[{"Conjugate", "[",
579 RowBox[{"FR$deltaZ", "[",
580 RowBox[{
581 RowBox[{"{",
582 RowBox[{"sb2", ",", "sb1"}], "}"}], ",",
583 RowBox[{"{",
584 RowBox[{"{", "}"}], "}"}]}], "]"}], "]"}]}], ")"}]}]}], "]"}],
585 ";"}]}]}]], "Input",
586 CellChangeTimes->{{3.7712646812949867`*^9, 3.771264714690206*^9}}],
587
588Cell[BoxData[
589 RowBox[{
590 RowBox[{"(*", " ",
591 RowBox[{
592 RowBox[{
593 "Implementing", " ", "the", " ", "naive", " ", "OS", " ", "scheme"}], " ",
594 "-", " ",
595 RowBox[{
596 RowBox[{"C", ".", " ", "trilinear"}], " ", "couplings"}]}], " ", "*)"}],
597 "\[IndentingNewLine]",
598 RowBox[{
599 RowBox[{
600 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
601 RowBox[{
602 RowBox[{"FR$delta", "[",
603 RowBox[{
604 RowBox[{"{",
605 RowBox[{"Au", "[",
606 RowBox[{"3", ",", "3"}], "]"}], "}"}], ",",
607 RowBox[{"{", "}"}]}], "]"}], ",",
608 RowBox[{
609 RowBox[{"1", "/", "MT"}], "*",
610 RowBox[{"(",
611 RowBox[{
612 RowBox[{"(",
613 RowBox[{
614 RowBox[{"Mst1", "*",
615 RowBox[{"FR$delta", "[",
616 RowBox[{
617 RowBox[{"{", "Mst1", "}"}], ",",
618 RowBox[{"{", "}"}]}], "]"}]}], "-",
619 RowBox[{"Mst2", "*",
620 RowBox[{"FR$delta", "[",
621 RowBox[{
622 RowBox[{"{", "Mst2", "}"}], ",",
623 RowBox[{"{", "}"}]}], "]"}]}]}], ")"}], "+",
624 RowBox[{
625 RowBox[{"1", "/", "MT"}], "*",
626 RowBox[{"(",
627 RowBox[{
628 RowBox[{"Mst1", "^", "2"}], "-",
629 RowBox[{"Mst2", "^", "2"}]}], ")"}], "*",
630 RowBox[{"Rtop", "[",
631 RowBox[{"1", ",", "1"}], "]"}], "*",
632 RowBox[{"Rtop", "[",
633 RowBox[{"1", ",", "2"}], "]"}], "*",
634 RowBox[{"FR$delta", "[",
635 RowBox[{
636 RowBox[{"{", "MT", "}"}], ",",
637 RowBox[{"{", "}"}]}], "]"}]}], "+",
638 RowBox[{
639 RowBox[{"(",
640 RowBox[{
641 RowBox[{"Mst1", "^", "2"}], "-",
642 RowBox[{"Mst2", "^", "2"}]}], ")"}], "*",
643 RowBox[{"(",
644 RowBox[{
645 RowBox[{
646 RowBox[{"Rtop", "[",
647 RowBox[{"1", ",", "1"}], "]"}], "*",
648 RowBox[{"FR$delta", "[",
649 RowBox[{
650 RowBox[{"{",
651 RowBox[{"Rtop", "[",
652 RowBox[{"1", ",", "2"}], "]"}], "}"}], ",",
653 RowBox[{"{", "}"}]}], "]"}]}], "+",
654 RowBox[{
655 RowBox[{"Rtop", "[",
656 RowBox[{"1", ",", "2"}], "]"}], "*",
657 RowBox[{"FR$delta", "[",
658 RowBox[{
659 RowBox[{"{",
660 RowBox[{"Rtop", "[",
661 RowBox[{"1", ",", "1"}], "]"}], "}"}], ",",
662 RowBox[{"{", "}"}]}], "]"}]}]}], ")"}]}]}], ")"}]}]}], "]"}],
663 ";"}], "\[IndentingNewLine]",
664 RowBox[{
665 RowBox[{"MoGRe`AddRenormalizationCondition", "[",
666 RowBox[{
667 RowBox[{"FR$delta", "[",
668 RowBox[{
669 RowBox[{"{",
670 RowBox[{"Ad", "[",
671 RowBox[{"3", ",", "3"}], "]"}], "}"}], ",",
672 RowBox[{"{", "}"}]}], "]"}], ",",
673 RowBox[{
674 RowBox[{"1", "/", "MB"}], "*",
675 RowBox[{"(",
676 RowBox[{
677 RowBox[{"(",
678 RowBox[{
679 RowBox[{"Msb1", "*",
680 RowBox[{"FR$delta", "[",
681 RowBox[{
682 RowBox[{"{", "Msb1", "}"}], ",",
683 RowBox[{"{", "}"}]}], "]"}]}], "-",
684 RowBox[{"Msb2", "*",
685 RowBox[{"FR$delta", "[",
686 RowBox[{
687 RowBox[{"{", "Msb2", "}"}], ",",
688 RowBox[{"{", "}"}]}], "]"}]}]}], ")"}], "+",
689 RowBox[{
690 RowBox[{"1", "/", "MB"}], "*",
691 RowBox[{"(",
692 RowBox[{
693 RowBox[{"Msb1", "^", "2"}], "-",
694 RowBox[{"Msb2", "^", "2"}]}], ")"}], "*",
695 RowBox[{"Rbot", "[",
696 RowBox[{"1", ",", "1"}], "]"}], "*",
697 RowBox[{"Rbot", "[",
698 RowBox[{"1", ",", "2"}], "]"}], "*",
699 RowBox[{"FR$delta", "[",
700 RowBox[{
701 RowBox[{"{", "MB", "}"}], ",",
702 RowBox[{"{", "}"}]}], "]"}]}], "+",
703 RowBox[{
704 RowBox[{"(",
705 RowBox[{
706 RowBox[{"Msb1", "^", "2"}], "-",
707 RowBox[{"Msb2", "^", "2"}]}], ")"}], "*",
708 RowBox[{"(",
709 RowBox[{
710 RowBox[{
711 RowBox[{"Rbot", "[",
712 RowBox[{"1", ",", "1"}], "]"}], "*",
713 RowBox[{"FR$delta", "[",
714 RowBox[{
715 RowBox[{"{",
716 RowBox[{"Rbot", "[",
717 RowBox[{"1", ",", "2"}], "]"}], "}"}], ",",
718 RowBox[{"{", "}"}]}], "]"}]}], "+",
719 RowBox[{
720 RowBox[{"Rbot", "[",
721 RowBox[{"1", ",", "2"}], "]"}], "*",
722 RowBox[{"FR$delta", "[",
723 RowBox[{
724 RowBox[{"{",
725 RowBox[{"Rbot", "[",
726 RowBox[{"1", ",", "1"}], "]"}], "}"}], ",",
727 RowBox[{"{", "}"}]}], "]"}]}]}], ")"}]}]}], "\n", ")"}]}]}],
728 "]"}], ";"}]}]}]], "Input",
729 CellChangeTimes->{{3.7712646987797956`*^9, 3.771264709832576*^9}}]
730}, Open ]]
731}, Open ]],
732
733Cell[CellGroupData[{
734
735Cell["Renormalisation and check of the output", "Subsection",
736 CellChangeTimes->{{3.771264772072339*^9, 3.771264780115871*^9}}],
737
738Cell[CellGroupData[{
739
740Cell["Computation of the renormalised Lagrangian", "Subsubsection",
741 CellChangeTimes->{{3.771264736607092*^9, 3.771264744086411*^9}}],
742
743Cell[BoxData[
744 RowBox[{
745 RowBox[{"LReno", " ", "=",
746 RowBox[{"MoGRe$Renormalize", "[", "lagr", "]"}]}], ";"}]], "Input",
747 CellChangeTimes->{{3.771264750774624*^9, 3.771264754345985*^9}}]
748}, Open ]],
749
750Cell[CellGroupData[{
751
752Cell["Check of the output", "Subsubsection",
753 CellChangeTimes->{{3.771265709950657*^9, 3.771265712926222*^9}}],
754
755Cell[BoxData[
756 RowBox[{
757 RowBox[{"(*",
758 RowBox[{"The", " ", "renormalization", " ", "conditions", " ",
759 RowBox[{"(",
760 RowBox[{"this", " ", "is", " ", "an", " ", "input"}], ")"}]}], " ",
761 "*)"}], "\[IndentingNewLine]", "MoGRe`SchemeRules"}]], "Input",
762 CellChangeTimes->{{3.771265755382235*^9, 3.771265786258514*^9}}],
763
764Cell[BoxData[
765 RowBox[{
766 RowBox[{"(*", " ", "Masses", " ", "*)"}], "\[IndentingNewLine]",
767 RowBox[{"RenormalizationRules", "[", "\"\<Masses\>\"", "]"}]}]], "Input",
768 CellChangeTimes->{{3.7451518713843946`*^9, 3.745151922162116*^9}, {
769 3.745164095936433*^9, 3.745164096200345*^9}, {3.745164380017013*^9,
770 3.745164380432691*^9}, {3.7451660736095877`*^9, 3.7451660750228863`*^9}, {
771 3.745213815312695*^9, 3.745213818856022*^9}, {3.74523151148385*^9,
772 3.7452315119531517`*^9}, {3.745303047472493*^9, 3.7453030478622513`*^9}, {
773 3.745303157313552*^9, 3.745303157889686*^9}, 3.745303406136817*^9, {
774 3.745574845804298*^9, 3.745574846145966*^9}, {3.745575407401681*^9,
775 3.7455754079116373`*^9}, {3.745578284288937*^9, 3.745578290465178*^9}, {
776 3.745582587034992*^9, 3.745582592977302*^9}, {3.7455826657703257`*^9,
777 3.7455826676585293`*^9}, {3.745582880576159*^9, 3.7455828928842697`*^9}, {
778 3.745583654476617*^9, 3.74558365461695*^9}, {3.7455837983158484`*^9,
779 3.745583798571127*^9}, 3.74564688658142*^9, {3.745647697834746*^9,
780 3.745647728311248*^9}, {3.745647806746401*^9, 3.745647810713956*^9}, {
781 3.7456480115435867`*^9, 3.7456481102599573`*^9}, {3.745663218654118*^9,
782 3.74566325795298*^9}, {3.745663711108547*^9, 3.745663711921837*^9}, {
783 3.745664612538536*^9, 3.7456646408908463`*^9}, 3.7456654013484077`*^9, {
784 3.745665650589293*^9, 3.7456656519560213`*^9}, {3.745665936303863*^9,
785 3.74566593728722*^9}, {3.7456688354046164`*^9, 3.7456688378674917`*^9}, {
786 3.745669024959618*^9, 3.745669030940509*^9}, {3.745669163428668*^9,
787 3.7456691987440863`*^9}, {3.745686812472076*^9, 3.74568681368325*^9}, {
788 3.745686880794334*^9, 3.74568688213133*^9}, {3.7456871775241747`*^9,
789 3.745687178863327*^9}, {3.745857725037485*^9, 3.745857735208394*^9}, {
790 3.7458586075203047`*^9, 3.745858613138812*^9}, {3.7458627781148767`*^9,
791 3.745862788627599*^9}, 3.745863379016315*^9, {3.745895772747778*^9,
792 3.7458957948440533`*^9}, {3.7458958293849087`*^9, 3.745895832308837*^9}, {
793 3.7463392933769217`*^9, 3.746339293953823*^9}, {3.746340958247019*^9,
794 3.746340980757111*^9}, {3.746446967049057*^9, 3.7464469672316236`*^9}, {
795 3.757309129229148*^9, 3.757309129577918*^9}, {3.757311847604082*^9,
796 3.757311850951976*^9}, {3.7573119060012207`*^9, 3.7573119133483257`*^9}, {
797 3.757400119485875*^9, 3.7574001368021793`*^9}, {3.757402948489932*^9,
798 3.757402960729457*^9}, {3.757403682464596*^9, 3.757403685895455*^9}, {
799 3.757406015321528*^9, 3.757406016888089*^9}, {3.75741526102887*^9,
800 3.757415265769416*^9}, {3.7574152965611763`*^9, 3.757415298314205*^9}, {
801 3.757416463439416*^9, 3.757416466010469*^9}, {3.757417500086125*^9,
802 3.757417523855345*^9}, {3.7574201860450373`*^9, 3.757420186670478*^9}, {
803 3.757448602343779*^9, 3.757448611302442*^9}, 3.757449335625184*^9, {
804 3.757449874383556*^9, 3.757449880587248*^9}, {3.757912547323681*^9,
805 3.757912548642529*^9}, {3.7580005002291107`*^9, 3.758000500760256*^9}, {
806 3.758005087249868*^9, 3.758005090807074*^9}, {3.758006435804907*^9,
807 3.758006444604958*^9}, {3.758161578806641*^9, 3.758161582999338*^9}, {
808 3.758164979635907*^9, 3.758164982381565*^9}, {3.758165013510289*^9,
809 3.758165029147532*^9}, {3.758165070730723*^9, 3.7581650858708963`*^9}, {
810 3.758165462675087*^9, 3.758165464347835*^9}, {3.758166197918043*^9,
811 3.7581662072139597`*^9}, {3.7581671637742367`*^9, 3.75816718202489*^9}, {
812 3.758177963109878*^9, 3.758177981847004*^9}, {3.758178134359646*^9,
813 3.758178159140918*^9}, {3.7582313015106*^9, 3.758231315264636*^9}, {
814 3.758248921123448*^9, 3.758248931762648*^9}, {3.758250734518733*^9,
815 3.758250783798354*^9}, {3.758250835228381*^9, 3.758250845550149*^9}, {
816 3.7582511056739693`*^9, 3.7582511090483103`*^9}, {3.75825131359096*^9,
817 3.758251325445635*^9}, {3.758251381280675*^9, 3.758251384116638*^9},
818 3.758251422640869*^9, {3.75825319671157*^9, 3.758253224083932*^9}, {
819 3.758257134438253*^9, 3.758257136759152*^9}, {3.758257209158524*^9,
820 3.7582572123266*^9}, {3.7582574013488417`*^9, 3.758257401749174*^9}, {
821 3.7582575392174883`*^9, 3.758257543960134*^9}, {3.758257710440304*^9,
822 3.758257731760181*^9}, {3.758257835439458*^9, 3.7582578864437447`*^9}, {
823 3.758258210240507*^9, 3.758258254261258*^9}, {3.758259393453458*^9,
824 3.7582594699190207`*^9}, {3.758259849130109*^9, 3.758259849496891*^9},
825 3.760938752558837*^9, {3.760938791593569*^9, 3.76093879283946*^9}, {
826 3.76094336376334*^9, 3.760943369090857*^9}, {3.760943562738961*^9,
827 3.760943588064082*^9}, 3.760944298912347*^9, {3.760952987831949*^9,
828 3.760952987979011*^9}, {3.760953226892173*^9, 3.760953236996063*^9}, {
829 3.76095348865825*^9, 3.760953489028057*^9}, {3.763457602937347*^9,
830 3.763457611310288*^9}, 3.763982907134116*^9, {3.763983088407049*^9,
831 3.763983089444828*^9}, {3.763983130719948*^9, 3.763983131117784*^9}, {
832 3.763983964034717*^9, 3.763983965670892*^9}, {3.763987056564139*^9,
833 3.763987059126238*^9}, {3.764005252103383*^9, 3.764005252336042*^9}, {
834 3.7640062913058977`*^9, 3.7640062939055023`*^9}, {3.7640073431045647`*^9,
835 3.764007410995483*^9}, {3.76452231518419*^9, 3.764522315471777*^9},
836 3.764522354075388*^9, {3.7645232675315247`*^9, 3.764523275689969*^9}, {
837 3.7645238003142233`*^9, 3.764523801367539*^9}, {3.764523872479026*^9,
838 3.764523881589555*^9}, {3.771263472262384*^9, 3.771263478300324*^9}, {
839 3.771263615037691*^9, 3.77126364922724*^9}, {3.771264545348506*^9,
840 3.771264606885642*^9}, 3.771264660161273*^9, 3.771264747418335*^9, {
841 3.771265717192766*^9, 3.771265723035942*^9}}],
842
843Cell[BoxData[
844 RowBox[{
845 RowBox[{"(*", " ",
846 RowBox[{"External", " ", "parameters"}], " ", "*)"}],
847 "\[IndentingNewLine]",
848 RowBox[{"RenormalizationRules", "[", "\"\<Externals\>\"", "]"}]}]], "Input",\
849
850 CellChangeTimes->{{3.771265730383052*^9, 3.771265734792964*^9}}],
851
852Cell[BoxData[
853 RowBox[{
854 RowBox[{"(*", " ",
855 RowBox[{"Internal", " ", "parameters"}], " ", "*)"}],
856 "\[IndentingNewLine]",
857 RowBox[{"RenormalizationRules", "[", "\"\<Internals\>\"", "]"}]}]], "Input",\
858
859 CellChangeTimes->{{3.758258382884212*^9, 3.758258384402162*^9}, {
860 3.7712657927013063`*^9, 3.77126579615121*^9}}]
861}, Open ]]
862}, Open ]],
863
864Cell[CellGroupData[{
865
866Cell["Writing the FeynArts model file (to run NLOCT with it)", "Subsection",
867 CellChangeTimes->{{3.77126584929543*^9, 3.771265858635962*^9}}],
868
869Cell[BoxData[{
870 RowBox[{
871 RowBox[{"SetDirectory", "[",
872 RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
873 ";"}], "\[IndentingNewLine]",
874 RowBox[{
875 RowBox[{"WriteFeynArtsOutput", "[",
876 RowBox[{"LReno", ",",
877 RowBox[{"FlavorExpand", "\[Rule]", "True"}], ",",
878 RowBox[{"Output", "\[Rule]", "\"\<MSSMatnlo_FA\>\""}]}], "]"}],
879 ";"}]}], "Input",
880 CellChangeTimes->{{3.7712658774901943`*^9, 3.771265895113875*^9}}]
881}, Open ]],
882
883Cell[CellGroupData[{
884
885Cell["Running NLOCT", "Subsection",
886 CellChangeTimes->{{3.7712672044251623`*^9, 3.771267206204439*^9}}],
887
888Cell[BoxData[
889 RowBox[{"Quit", "[", "]"}]], "Input",
890 CellChangeTimes->{{3.77126739224209*^9, 3.77126739298962*^9}}],
891
892Cell[BoxData[
893 RowBox[{
894 RowBox[{"(*", " ",
895 RowBox[{"Loading", " ", "NLOCTR", " ", "and", " ", "FeynArts"}], " ",
896 "*)"}], "\[IndentingNewLine]",
897 RowBox[{
898 RowBox[{"<<", " ", "FeynArts`"}], "\[IndentingNewLine]",
899 RowBox[{
900 RowBox[{"SetDirectory", "[", "\"\<<Path to NLOCT>\>\"", "]"}], ";"}],
901 "\[IndentingNewLine]",
902 RowBox[{"<<", " ", "NLOCT`"}], "\[IndentingNewLine]",
903 RowBox[{
904 RowBox[{"SetDirectory", "[",
905 RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]}]}]], "Input",
906 CellChangeTimes->{{3.77126720952652*^9, 3.7712672144737177`*^9}, {
907 3.771551625248365*^9, 3.771551629053658*^9}}],
908
909Cell[BoxData[
910 RowBox[{
911 RowBox[{"(*", " ",
912 RowBox[{"Running", " ", "NLOCT"}], " ", "*)"}], "\[IndentingNewLine]",
913 RowBox[{
914 RowBox[{"WriteCT", "[",
915 RowBox[{
916 "\"\<MSSMatnlo_FA/MSSMatnlo_FA\>\"", ",",
917 "\"\<MSSMatnlo_FA/MSSMatnlo_FA\>\"", ",",
918 RowBox[{"Output", "->", "\"\<MSSMatnlo\>\""}], ",",
919 RowBox[{"LabelInternal", "\[Rule]", "True"}], ",",
920 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
921 RowBox[{"KeptIndices", "\[Rule]",
922 RowBox[{"{", "}"}]}], ",",
923 RowBox[{"ZeroMom", "\[Rule]",
924 RowBox[{"{",
925 RowBox[{"{",
926 RowBox[{"aS", ",",
927 RowBox[{"{",
928 RowBox[{
929 RowBox[{"F", "[", "14", "]"}], ",",
930 RowBox[{"V", "[", "4", "]"}], ",",
931 RowBox[{"-",
932 RowBox[{"F", "[", "14", "]"}]}]}], "}"}], ",", "0"}], "}"}],
933 "}"}]}], ",",
934 RowBox[{"Assumptions", "->",
935 RowBox[{"{",
936 RowBox[{
937 RowBox[{"MT", ">", "0"}], ",",
938 RowBox[{"MsuL", ">", "0"}], ",",
939 RowBox[{"MscL", ">", "0"}], ",",
940 RowBox[{"Mst1", ">", "0"}], ",",
941 RowBox[{"MsuR", ">", "0"}], ",",
942 RowBox[{"MscR", ">", "0"}], ",",
943 RowBox[{"Mst2", ">", "0"}], ",",
944 RowBox[{"MsdL", ">", "0"}], ",",
945 RowBox[{"MssL", ">", "0"}], ",",
946 RowBox[{"Msb1", ">", "0"}], ",",
947 RowBox[{"MsdR", ">", "0"}], ",",
948 RowBox[{"MssR", ">", "0"}], ",",
949 RowBox[{"Msb2", ">", "0"}], ",",
950 RowBox[{"Mgo", ">", "0"}]}], "}"}]}], ",",
951 RowBox[{"CTparameters", "\[Rule]", "True"}], ",",
952 RowBox[{"MaxDim", "\[Rule]", "4"}], ",",
953 RowBox[{"EvenOnly", "\[Rule]", "False"}], ",",
954 RowBox[{"IsFeynmanGauge", "\[Rule]", "True"}], ",",
955 RowBox[{"GenericVertexList", "\[Rule]",
956 RowBox[{"{",
957 RowBox[{
958 RowBox[{"{", "V", "}"}], ",",
959 RowBox[{"{", "S", "}"}], ",",
960 RowBox[{"{",
961 RowBox[{"F", ",", "F"}], "}"}], ",",
962 RowBox[{"{",
963 RowBox[{"S", ",", "S"}], "}"}], ",",
964 RowBox[{"{",
965 RowBox[{"V", ",", "V"}], "}"}], ",",
966 RowBox[{"{",
967 RowBox[{"V", ",", "S"}], "}"}], ",",
968 RowBox[{"{",
969 RowBox[{"F", ",", "F", ",", "S"}], "}"}], ",",
970 RowBox[{"{",
971 RowBox[{"F", ",", "F", ",", "V"}], "}"}], ",",
972 RowBox[{"{",
973 RowBox[{"V", ",", "V", ",", "V"}], "}"}], ",",
974 RowBox[{"{",
975 RowBox[{"V", ",", "V", ",", "S"}], "}"}], ",",
976 RowBox[{"{",
977 RowBox[{"V", ",", "S", ",", "S"}], "}"}], ",",
978 RowBox[{"{",
979 RowBox[{"S", ",", "S", ",", "S"}], "}"}], ",",
980 RowBox[{"{",
981 RowBox[{"V", ",", "V", ",", "V", ",", "V"}], "}"}], ",",
982 RowBox[{"{",
983 RowBox[{"V", ",", "V", ",", "V", ",", "S"}], "}"}], ",",
984 RowBox[{"{",
985 RowBox[{"V", ",", "V", ",", "S", ",", "S"}], "}"}], ",",
986 RowBox[{"{",
987 RowBox[{"V", ",", "S", ",", "S", ",", "S"}], "}"}]}],
988 RowBox[{"(*",
989 RowBox[{",",
990 RowBox[{"{",
991 RowBox[{"S", ",", "S", ",", "S", ",", "S"}], "}"}]}], "*)"}],
992 "}"}]}]}], "]"}], ";"}], " ",
993 RowBox[{"(*",
994 RowBox[{
995 "No", " ", "need", " ", "of", " ", "the", " ", "4", " ", "scalar", " ",
996 "R2"}], " ", "*)"}]}]], "Input",
997 CellChangeTimes->{{3.7712672183822308`*^9, 3.7712672677743483`*^9}, {
998 3.7712673246381702`*^9, 3.771267400685985*^9}, {3.7715516440984097`*^9,
999 3.771551654020927*^9}}]
1000}, Open ]],
1001
1002Cell[CellGroupData[{
1003
1004Cell["Generating the UFO model files", "Subsection",
1005 CellChangeTimes->{{3.518892302672626*^9, 3.518892305688608*^9}, {
1006 3.6485269255019417`*^9, 3.648526926658614*^9}, {3.648820158462558*^9,
1007 3.648820160261292*^9}}],
1008
1009Cell[BoxData[
1010 RowBox[{"Quit", "[", "]"}]], "Input",
1011 CellChangeTimes->{{3.643510536070684*^9, 3.643510536921255*^9}, {
1012 3.648885571488657*^9, 3.648885572262779*^9}, {3.6490120454401617`*^9,
1013 3.649012045967602*^9}}],
1014
1015Cell[BoxData[
1016 RowBox[{
1017 RowBox[{"(*",
1018 RowBox[{"Loading", " ", "FeynRules"}], " ", "*)"}], "\[IndentingNewLine]",
1019
1020 RowBox[{
1021 RowBox[{
1022 RowBox[{"$FeynRulesPath", "=",
1023 RowBox[{"SetDirectory", "[", "\"\<<FeynRules path\>\"", "]"}]}], ";"}],
1024 "\[IndentingNewLine]",
1025 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
1026 RowBox[{
1027 RowBox[{"LoadModel", "[", "\"\<<Path to mssm-nlo.fr>\>\"", "]"}], ";"}],
1028 "\[IndentingNewLine]",
1029 RowBox[{
1030 RowBox[{"LoadRestriction", "[", "\"\<<Path to Massless.rst>\>\"", "]"}],
1031 ";"}], "\[IndentingNewLine]",
1032 RowBox[{
1033 RowBox[{"LoadRestriction", "[", "\"\<<Path to DiagonalCKM.rst\>\"", "]"}],
1034 ";"}], "\[IndentingNewLine]",
1035 RowBox[{
1036 RowBox[{"SetDirectory", "[",
1037 RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]}]}]], "Input",
1038 CellChangeTimes->{{3.771899971227213*^9, 3.7718999799282093`*^9}}],
1039
1040Cell[BoxData[
1041 RowBox[{
1042 RowBox[{"(*",
1043 RowBox[{"Getting", " ", "the", " ", "NLO", " ", "file"}], "*)"}],
1044 "\[IndentingNewLine]",
1045 RowBox[{
1046 RowBox[{"Get", "[", "\"\<MSSMatnlo.nlo\>\"", "]"}], ";"}]}]], "Input"],
1047
1048Cell[BoxData[
1049 RowBox[{
1050 RowBox[{"(*", " ",
1051 RowBox[{
1052 "Computing", " ", "the", " ", "Lagrangian", " ", "if", " ", "needed"}],
1053 " ", "*)"}], "\[IndentingNewLine]",
1054 RowBox[{
1055 RowBox[{
1056 RowBox[{"SetDirectory", "[",
1057 RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}],
1058 "\[IndentingNewLine]",
1059 RowBox[{
1060 RowBox[{"If", "[",
1061 RowBox[{
1062 RowBox[{"!",
1063 RowBox[{"FileExistsQ", "[", "\"\<lagrangian.dat\>\"", "]"}]}], ",",
1064 "\[IndentingNewLine]",
1065 RowBox[{
1066 RowBox[{"lagr", "=", "LSUSY"}], ";", "\[IndentingNewLine]",
1067 RowBox[{
1068 RowBox[{"Definition", "[", "lagr", "]"}], " ", ">>", " ",
1069 "lagrangian.dat"}]}], ",", "\[IndentingNewLine]",
1070 RowBox[{
1071 RowBox[{"<<", "lagrangian.dat"}], ";"}]}], "\[IndentingNewLine]",
1072 "]"}], ";"}]}]}]], "Input"],
1073
1074Cell[BoxData[
1075 RowBox[{
1076 RowBox[{"(*", " ",
1077 RowBox[{"SUSY", " ", "restoring", " ", "counterterms"}], " ", "*)"}],
1078 "\[IndentingNewLine]",
1079 RowBox[{
1080 RowBox[{"CTrules", "=",
1081 RowBox[{"FeynmanRules", "[", "LCT", "]"}]}], ";"}]}]], "Input",
1082 CellChangeTimes->{{3.616474426822394*^9, 3.616474431300906*^9}, {
1083 3.616474476497079*^9, 3.616474500189231*^9}, 3.616474613066862*^9, {
1084 3.616494779207423*^9, 3.616494783828376*^9}, {3.616580249965948*^9,
1085 3.616580268353311*^9}, {3.616580334915408*^9, 3.616580354074359*^9}, {
1086 3.616730859450924*^9, 3.6167308627148323`*^9}, {3.618221520121941*^9,
1087 3.618221562907246*^9}, {3.618575869806182*^9, 3.6185758768210573`*^9}, {
1088 3.622267142817052*^9, 3.622267158950058*^9}, {3.622517903231851*^9,
1089 3.6225179034762793`*^9}, {3.643130455506323*^9, 3.643130457362434*^9}, {
1090 3.6434346190475397`*^9, 3.643434619424014*^9}, {3.64348057990516*^9,
1091 3.643480585441132*^9}, {3.643510546826219*^9, 3.643510547036009*^9},
1092 3.643510612235943*^9, {3.644167903839278*^9, 3.6441679129268084`*^9}, {
1093 3.644168079035636*^9, 3.644168086658132*^9}, 3.644168122851153*^9, {
1094 3.644168157533206*^9, 3.644168163660873*^9}, {3.64428642875541*^9,
1095 3.644286429840966*^9}, {3.644286474789322*^9, 3.6442864751012707`*^9}, {
1096 3.6442907764504147`*^9, 3.644290780850091*^9}, {3.644295335498075*^9,
1097 3.644295420355406*^9}, {3.644313480427791*^9, 3.644313483330103*^9}, {
1098 3.6443187949698143`*^9, 3.6443187975911818`*^9}, {3.644322102713567*^9,
1099 3.64432210330866*^9}, {3.644323739975934*^9, 3.644323741558518*^9}, {
1100 3.64455865826359*^9, 3.644558694344214*^9}, {3.644560355751863*^9,
1101 3.644560360988126*^9}, {3.644666766031262*^9, 3.6446667688757353`*^9},
1102 3.648465060332839*^9, 3.6484702356962557`*^9, 3.648475630789858*^9,
1103 3.648488825342585*^9, 3.6485491231721897`*^9, 3.648559751598288*^9, {
1104 3.648744806616398*^9, 3.648744811285426*^9}, {3.648829765484006*^9,
1105 3.648829767668827*^9}, {3.648838329122327*^9, 3.6488383400295486`*^9},
1106 3.648841075097126*^9, {3.648868134543252*^9, 3.648868143331616*^9}, {
1107 3.648878569185281*^9, 3.64887857015184*^9}, {3.648878724046915*^9,
1108 3.648878732013793*^9}, {3.648879167019032*^9, 3.648879202669491*^9},
1109 3.6488800001412497`*^9, {3.648915215173975*^9, 3.648915216873856*^9}, {
1110 3.648915859118635*^9, 3.648915860842957*^9}, {3.648919484393269*^9,
1111 3.648919485853923*^9}, {3.648921190392002*^9, 3.648921192470207*^9}, {
1112 3.648984982197277*^9, 3.648984984755639*^9}, {3.6489853031837873`*^9,
1113 3.648985307937234*^9}, {3.648985489148918*^9, 3.648985491780357*^9}, {
1114 3.648985592485896*^9, 3.648985611274666*^9}, 3.6489872177669*^9, {
1115 3.648987289747376*^9, 3.648987291297887*^9}, 3.648988818570665*^9, {
1116 3.648988996197752*^9, 3.648989001188278*^9}, {3.649012040072267*^9,
1117 3.649012042569594*^9}, {3.649015200388323*^9, 3.6490152033349524`*^9}, {
1118 3.6490792754392233`*^9, 3.6490793002540894`*^9}, {3.6490793377647142`*^9,
1119 3.649079339041106*^9}, {3.649089280363574*^9, 3.649089286072534*^9}, {
1120 3.6491015811843557`*^9, 3.649101589381792*^9}, {3.649130696400343*^9,
1121 3.6491307322830257`*^9}, {3.7718999879465647`*^9, 3.771900126473785*^9}}],
1122
1123Cell[BoxData[
1124 RowBox[{
1125 RowBox[{"(*", " ",
1126 RowBox[{"Writing", " ", "the", " ", "UFO"}], " ", "*)"}],
1127 "\[IndentingNewLine]",
1128 RowBox[{
1129 RowBox[{
1130 RowBox[{"WriteUFO", "[",
1131 RowBox[{"lagr", ",",
1132 RowBox[{"UVCounterterms", "\[Rule]",
1133 RowBox[{"(",
1134 RowBox[{"UV$vertlist", "/.",
1135 RowBox[{"{",
1136 RowBox[{
1137 RowBox[{"FR$IR", "->", "1"}], ",",
1138 RowBox[{"FR$IRLog", "->", "1"}]}], "}"}]}], ")"}]}], " ", ",",
1139 RowBox[{
1140 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], "/.",
1141 RowBox[{"{",
1142 RowBox[{
1143 RowBox[{"FR$IR", "->", "1"}], ",",
1144 RowBox[{"FR$IRLog", "->", "1"}]}], "}"}]}], ",",
1145 RowBox[{"UVLoopCounterterms", "\[Rule]", "CTrules"}], ",", " ",
1146 RowBox[{"Output", "\[Rule]", "\"\<MSSMatNLO_UFO\>\""}], ",",
1147 RowBox[{"AddDecays", "\[Rule]", "False"}], ",", " ",
1148 RowBox[{
1149 RowBox[{"CTParameters", "\[Rule]", "FR$CTparam"}], "/.",
1150 RowBox[{"{",
1151 RowBox[{
1152 RowBox[{"FR$IR", "->", "1"}], ",",
1153 RowBox[{"FR$IRLog", "->", "1"}]}], "}"}]}]}], "]"}], ";"}],
1154 "\[IndentingNewLine]"}]}]], "Input",
1155 CellChangeTimes->{{3.7719001200292463`*^9, 3.771900125275017*^9},
1156 3.772272367240961*^9}]
1157}, Open ]]
1158},
1159WindowSize->{1354, 800},
1160WindowMargins->{{Automatic, 0}, {-30, Automatic}},
1161ShowSelection->True,
1162FrontEndVersion->"10.4 for Mac OS X x86 (32-bit, 64-bit Kernel) (April 11, \
11632016)",
1164StyleDefinitions->"Default.nb"
1165]
1166(* End of Notebook Content *)
1167
1168(* Internal cache information *)
1169(*CellTagsOutline
1170CellTagsIndex->{}
1171*)
1172(*CellTagsIndex
1173CellTagsIndex->{}
1174*)
1175(*NotebookFileOutline
1176Notebook[{
1177Cell[545, 20, 463, 7, 28, "Input"],
1178Cell[CellGroupData[{
1179Cell[1033, 31, 191, 3, 44, "Subsection"],
1180Cell[1227, 36, 3498, 61, 165, "Input"]
1181}, Open ]],
1182Cell[CellGroupData[{
1183Cell[4762, 102, 263, 3, 44, "Subsection"],
1184Cell[5028, 107, 572, 11, 68, "Text"],
1185Cell[5603, 120, 1021, 24, 114, "Input"]
1186}, Open ]],
1187Cell[CellGroupData[{
1188Cell[6661, 149, 818, 12, 44, "Subsection"],
1189Cell[CellGroupData[{
1190Cell[7504, 165, 98, 1, 35, "Subsubsection"],
1191Cell[7605, 168, 363, 9, 46, "Input"],
1192Cell[7971, 179, 749, 19, 46, "Input"]
1193}, Open ]],
1194Cell[CellGroupData[{
1195Cell[8757, 203, 166, 2, 35, "Subsubsection"],
1196Cell[8926, 207, 2442, 63, 165, "Input"],
1197Cell[11371, 272, 414, 10, 63, "Input"],
1198Cell[11788, 284, 1630, 43, 114, "Input"],
1199Cell[13421, 329, 4227, 127, 97, "Input"],
1200Cell[17651, 458, 4232, 127, 97, "Input"],
1201Cell[21886, 587, 4941, 141, 148, "Input"]
1202}, Open ]]
1203}, Open ]],
1204Cell[CellGroupData[{
1205Cell[26876, 734, 127, 1, 44, "Subsection"],
1206Cell[CellGroupData[{
1207Cell[27028, 739, 133, 1, 35, "Subsubsection"],
1208Cell[27164, 742, 190, 4, 28, "Input"]
1209}, Open ]],
1210Cell[CellGroupData[{
1211Cell[27391, 751, 110, 1, 35, "Subsubsection"],
1212Cell[27504, 754, 332, 7, 46, "Input"],
1213Cell[27839, 763, 5597, 77, 46, "Input"],
1214Cell[33439, 842, 276, 7, 46, "Input"],
1215Cell[33718, 851, 326, 8, 46, "Input"]
1216}, Open ]]
1217}, Open ]],
1218Cell[CellGroupData[{
1219Cell[34093, 865, 141, 1, 44, "Subsection"],
1220Cell[34237, 868, 428, 11, 46, "Input"]
1221}, Open ]],
1222Cell[CellGroupData[{
1223Cell[34702, 884, 103, 1, 44, "Subsection"],
1224Cell[34808, 887, 116, 2, 28, "Input"],
1225Cell[34927, 891, 634, 15, 97, "Input"],
1226Cell[35564, 908, 3506, 90, 148, "Input"]
1227}, Open ]],
1228Cell[CellGroupData[{
1229Cell[39107, 1003, 218, 3, 44, "Subsection"],
1230Cell[39328, 1008, 218, 4, 28, "Input"],
1231Cell[39549, 1014, 900, 23, 131, "Input"],
1232Cell[40452, 1039, 221, 6, 46, "Input"],
1233Cell[40676, 1047, 841, 24, 131, "Input"],
1234Cell[41520, 1073, 3189, 47, 46, "Input"],
1235Cell[44712, 1122, 1267, 33, 80, "Input"]
1236}, Open ]]
1237}
1238]
1239*)
1240