MSSMatNLO: use-MoGRe.nb

File use-MoGRe.nb, 47.2 KB (added by BenjF, 4 months 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