TypeIISeesaw: use-MoGre.nb

File use-MoGre.nb, 26.5 KB (added by BenjF, 8 months ago)
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[     27030,        677]
14NotebookOptionsPosition[     25101,        610]
15NotebookOutlinePosition[     25478,        627]
16CellTagsIndexPosition[     25435,        624]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"Quit", "[", "]"}]], "Input",
23 CellChangeTimes->{{3.7638865422974358`*^9, 3.763886542860983*^9}, {
24  3.767372147789777*^9, 3.7673721483716993`*^9}, {3.780996846995326*^9,
25  3.780996847734174*^9}}],
26
27Cell[CellGroupData[{
28
29Cell["Setup", "Section",
30 CellChangeTimes->{{3.518892288721245*^9, 3.518892291440427*^9}, {
31  3.648272586359296*^9, 3.648272589589151*^9}, {3.744951163067635*^9,
32  3.744951164882929*^9}}],
33
34Cell[BoxData[{
35 RowBox[{
36  RowBox[{"$FeynRulesPath", "=",
37   RowBox[{"SetDirectory", "[", "\"\<<path-to-FR>\>\"", "]"}]}],
38  ";"}], "\[IndentingNewLine]",
39 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
40 RowBox[{
41  RowBox[{"SetDirectory", "[",
42   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
43  ";"}], "\[IndentingNewLine]",
44 RowBox[{
45  RowBox[{"LoadModel", "[",
46   RowBox[{"\"\<sm_reduced.fr\>\"", ",", "\"\<type_ii.fr\>\""}], "]"}],
47  ";"}], "\[IndentingNewLine]",
48 RowBox[{
49  RowBox[{"LoadRestriction", "[",
50   RowBox[{
51    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
52    "\"\</Massless_4f.rst\>\""}], "]"}], ";"}], "\[IndentingNewLine]",
53 RowBox[{
54  RowBox[{"LoadRestriction", "[",
55   RowBox[{
56    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
57    "\"\</DiagonalCKM.rst\>\""}], "]"}], ";"}]}], "Input",
58 CellChangeTimes->{{3.715342142509688*^9, 3.715342163031452*^9}, {
59   3.715344057238492*^9, 3.715344058291609*^9}, {3.7178506900799713`*^9,
60   3.717850690973333*^9}, {3.7180002809816017`*^9, 3.718000300883692*^9}, {
61   3.718000394980835*^9, 3.7180003951009827`*^9}, {3.745153113604156*^9,
62   3.745153113886346*^9}, {3.757403422964182*^9, 3.757403423697225*^9}, {
63   3.758251145793989*^9, 3.7582511468706083`*^9}, 3.7582513470889874`*^9, {
64   3.76335111930016*^9, 3.763351139093457*^9}, {3.7633511693483973`*^9,
65   3.763351173687544*^9}, {3.763886545633397*^9, 3.763886571907737*^9},
66   3.767099220292862*^9, {3.767374514400957*^9, 3.767374515454331*^9}, {
67   3.767375022418717*^9, 3.76737502255937*^9}, {3.7674580186709423`*^9,
68   3.767458019546134*^9}, {3.767610478088133*^9, 3.767610507175609*^9}, {
69   3.767610538533186*^9, 3.76761053951677*^9}, {3.767610603703107*^9,
70   3.7676106297846937`*^9}, {3.7676900131323767`*^9,
71   3.7676900151887627`*^9}, {3.76858524470041*^9, 3.768585245204138*^9}, {
72   3.77027747721585*^9, 3.7702775056436033`*^9}, {3.7711262488766813`*^9,
73   3.771126253973398*^9}, {3.771126559639118*^9, 3.771126602632146*^9},
74   3.771139675031245*^9, {3.771139748794265*^9, 3.771139751639614*^9}, {
75   3.7711399153197002`*^9, 3.7711399456001577`*^9}, {3.771139981523079*^9,
76   3.771139981644162*^9}, {3.77857888321222*^9, 3.778578897476334*^9}, {
77   3.77857892985773*^9, 3.7785789375119877`*^9}, 3.7785789770568542`*^9, {
78   3.77857907445473*^9, 3.7785791018908577`*^9}, {3.779597422325054*^9,
79   3.779597426820354*^9}, {3.780921597734684*^9, 3.780921601148004*^9},
80   3.780942736968416*^9, {3.780998835004087*^9, 3.780998837598987*^9}}]
81}, Open  ]],
82
83Cell[CellGroupData[{
84
85Cell["NLO", "Section",
86 CellChangeTimes->{{3.7795974887503777`*^9, 3.779597489081195*^9}}],
87
88Cell[CellGroupData[{
89
90Cell["Setup", "Subsubsection",
91 CellChangeTimes->{{3.7712636194863377`*^9, 3.771263622434663*^9}, {
92  3.780998867365963*^9, 3.780998885433375*^9}}],
93
94Cell[BoxData[{
95 RowBox[{
96  RowBox[{"SetDirectory", "[",
97   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
98  ";"}], "\[IndentingNewLine]",
99 RowBox[{
100  RowBox[{"Begin", "[", "\"\<MoGRe`\>\"", "]"}], ";"}], "\[IndentingNewLine]",
101 RowBox[{
102  RowBox[{"<<", "MoGRe_v1.1.m"}], ";"}], "\[IndentingNewLine]",
103 RowBox[{
104  RowBox[{"End", "[", "]"}], ";"}], "\[IndentingNewLine]",
105 RowBox[{
106  RowBox[{"MoGRe`ClearRenormalizationScheme", "[", "]"}], ";"}]}], "Input",
107 CellChangeTimes->{{3.779621739750206*^9, 3.7796217508767843`*^9}, {
108  3.7796273631453953`*^9, 3.77962740535894*^9}, {3.78099943900846*^9,
109  3.780999439837109*^9}}],
110
111Cell[BoxData[
112 RowBox[{
113  RowBox[{"(*", " ",
114   RowBox[{"Calculation", " ", "of", " ", "the", " ", "Lagrangian", " ",
115    RowBox[{"(",
116     RowBox[{"once", " ", "and", " ", "for", " ", "all"}], ")"}], " ", "if",
117    " ", "not", " ", "done", " ", "already"}], " ", "*)"}], " ",
118  "\[IndentingNewLine]",
119  RowBox[{
120   RowBox[{"If", "[",
121    RowBox[{
122     RowBox[{"!",
123      RowBox[{"FileExistsQ", "[", "\"\<lagrangian.dat\>\"", "]"}]}], ",",
124     "\[IndentingNewLine]",
125     RowBox[{
126      RowBox[{"lags", "=", " ",
127       RowBox[{"ExpandIndices", "/@", "LType2"}]}], ";",
128      "\[IndentingNewLine]",
129      RowBox[{
130       RowBox[{"Definition", "[", "lags", "]"}], " ", ">>", " ",
131       "lagrangian.dat"}]}], ",", "\[IndentingNewLine]",
132     RowBox[{
133      RowBox[{"<<", "lagrangian.dat"}], ";"}]}], "\[IndentingNewLine]", "]"}],
134    ";"}]}]], "Input",
135 CellChangeTimes->{{3.779597667101206*^9, 3.7795976737764177`*^9}, {
136  3.779621768389236*^9, 3.7796217817142982`*^9}, {3.780998896348439*^9,
137  3.7809989075559053`*^9}}],
138
139Cell[BoxData[
140 RowBox[{
141  RowBox[{"(*", " ",
142   RowBox[{
143    RowBox[{"To", " ", "get", " ", "more", " ", "print"}], "-",
144    RowBox[{"outs", " ", "to", " ", "the", " ", "screen"}]}], " ", "*)"}],
145  "\[IndentingNewLine]",
146  RowBox[{
147   RowBox[{"MoGRe$Debug", "=", "True"}], ";"}]}]], "Input",
148 CellChangeTimes->{{3.771263627141055*^9, 3.7712636335707293`*^9}}],
149
150Cell[BoxData[
151 RowBox[{
152  RowBox[{"(*", " ",
153   RowBox[{
154   "Setting", " ", "the", " ", "options", " ", "of", " ", "the", " ", "main",
155    " ", "function"}], " ", "*)"}], "\[IndentingNewLine]",
156  RowBox[{
157   RowBox[{"SetOptions", "[",
158    RowBox[{"MoGRe$Renormalize", ",", " ",
159     RowBox[{"Exclude4Scalars", "\[Rule]", "False"}], ",",
160     RowBox[{"FlavorMixing", "\[Rule]", "False"}], ",",
161     RowBox[{"CouplingOrders", "\[Rule]",
162      RowBox[{"{", "QCD", "}"}]}], ",",
163     RowBox[{"CanonicalTwoPoints", "\[Rule]", "True"}]}], "]"}],
164   ";"}]}]], "Input",
165 CellChangeTimes->{{3.771263641132224*^9, 3.7712636559080553`*^9}, {
166  3.7796218121813602`*^9, 3.7796218127058563`*^9}, {3.7796229843052807`*^9,
167  3.779622984874784*^9}, {3.780979808313796*^9, 3.7809798160041103`*^9}}]
168}, Open  ]],
169
170Cell[CellGroupData[{
171
172Cell["OS renormalization", "Subsubsection",
173 CellChangeTimes->{{3.771264549679812*^9, 3.771264557709045*^9}, {
174  3.771264725891897*^9, 3.771264729628981*^9}, {3.779622421220426*^9,
175  3.779622424683118*^9}}],
176
177Cell[BoxData[
178 RowBox[{
179  RowBox[{"(*", " ",
180   RowBox[{"Unrenormalized", " ", "parameters"}], " ", "*)"}],
181  "\[IndentingNewLine]",
182  RowBox[{
183   RowBox[{"MoGRe`DefineUnrenormalizedParameters", "[",
184    RowBox[{"{", "\[IndentingNewLine]",
185     RowBox[{
186     "Gf", ",", "aEWM1", ",", "MZ", ",", "delCP", ",", "dmsq31", ",",
187      "dmsq21", ",", "th12", ",", "th13", ",", "th23", ",", "phiM1", ",",
188      "phiM2", ",", "vevD", ",", "lamHD1", ",", "lamD1"}], "}"}], "]"}],
189   ";"}]}]], "Input",
190 CellChangeTimes->{{3.771264570309547*^9, 3.771264574380278*^9}, {
191  3.779622328795891*^9, 3.779622332835644*^9}, {3.779623239502658*^9,
192  3.7796232588028*^9}, {3.779623369643303*^9, 3.7796233703313828`*^9}, {
193  3.779625801787382*^9, 3.779625812383876*^9}, {3.779625921812018*^9,
194  3.7796259254903584`*^9}, {3.779626038445567*^9, 3.779626041360046*^9}, {
195  3.779626348010645*^9, 3.779626365854025*^9}, {3.779626431993538*^9,
196  3.7796264326080313`*^9}, {3.7796265337645187`*^9, 3.779626563396529*^9}, {
197  3.780942696192215*^9, 3.780942696568851*^9}}],
198
199Cell[BoxData[
200 RowBox[{
201  RowBox[{"(*", " ", "Fields", " ", "*)"}], "\[IndentingNewLine]",
202  RowBox[{
203   RowBox[{
204    RowBox[{"MoGRe`DeclareUnrenormalizedFields", "[",
205     RowBox[{"W", ",", "A", ",", "Z"}], "]"}], ";"}], "\[IndentingNewLine]",
206   RowBox[{
207    RowBox[{"MoGRe`RealFieldRenormalization", "[", "]"}], " ",
208    ";"}]}]}]], "Input",
209 CellChangeTimes->{{3.771264582333989*^9, 3.771264608943572*^9}}],
210
211Cell[BoxData[
212 RowBox[{
213  RowBox[{"(*", " ",
214   RowBox[{"Removing", " ", "some", " ", "parameters", " ",
215    RowBox[{"(",
216     RowBox[{
217      RowBox[{"cf", ".", " ", "the"}], " ", "way", " ", "NLOCT", " ",
218      "works"}], ")"}]}], " ", "*)"}], "\[IndentingNewLine]",
219  RowBox[{
220   RowBox[{
221    RowBox[{
222     RowBox[{"MoGRe`RemovingInternalCst", "[", "#", "]"}], "&"}], "/@",
223    RowBox[{"{", "\[IndentingNewLine]",
224     RowBox[{"gs", ",", "\[IndentingNewLine]",
225      RowBox[{"Sequence", "@@",
226       RowBox[{"Flatten", "[",
227        RowBox[{"Table", "[",
228         RowBox[{
229          RowBox[{"{",
230           RowBox[{
231            RowBox[{"yu", "[",
232             RowBox[{"i", ",", "i"}], "]"}], ",",
233            RowBox[{"yd", "[",
234             RowBox[{"i", ",", "i"}], "]"}]}], "}"}], ",",
235          RowBox[{"{",
236           RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], ",",
237      "\[IndentingNewLine]",
238      RowBox[{"Sequence", "@@",
239       RowBox[{"Flatten", "[",
240        RowBox[{"Table", "[",
241         RowBox[{
242          RowBox[{"{",
243           RowBox[{"CKM", "[",
244            RowBox[{"i", ",", "j"}], "]"}], "}"}], ",",
245          RowBox[{"{",
246           RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",",
247          RowBox[{"{",
248           RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}]}],
249     "\[IndentingNewLine]", "}"}]}], ";"}]}]], "Input",
250 CellChangeTimes->{{3.771264625204327*^9, 3.771264640138379*^9},
251   3.779622343993898*^9}]
252}, Open  ]],
253
254Cell[CellGroupData[{
255
256Cell["Computation of the renormalised Lagrangian", "Subsubsection",
257 CellChangeTimes->{{3.771264736607092*^9, 3.771264744086411*^9}}],
258
259Cell[BoxData[
260 RowBox[{
261  RowBox[{"LReno", " ", "=",
262   RowBox[{"MoGRe$Renormalize", "[",
263    RowBox[{"Plus", "@@", "lags"}], "]"}]}], ";"}]], "Input",
264 CellChangeTimes->{{3.771264750774624*^9, 3.771264754345985*^9}, {
265   3.779622437680065*^9, 3.7796224387360353`*^9}, 3.779624783600335*^9, {
266   3.7796274220329533`*^9, 3.779627429379524*^9}, {3.7796294251206827`*^9,
267   3.779629444949071*^9}, {3.779629548413604*^9, 3.7796295518915253`*^9}}]
268}, Open  ]],
269
270Cell[CellGroupData[{
271
272Cell["Check of the output", "Subsubsection",
273 CellChangeTimes->{{3.771265709950657*^9, 3.771265712926222*^9}}],
274
275Cell[BoxData[
276 RowBox[{
277  RowBox[{"(*",
278   RowBox[{"The", " ", "renormalization", " ", "conditions", " ",
279    RowBox[{"(",
280     RowBox[{"this", " ", "is", " ", "an", " ", "input"}], ")"}]}], " ",
281   "*)"}], "\[IndentingNewLine]", "MoGRe`SchemeRules"}]], "Input",
282 CellChangeTimes->{{3.771265755382235*^9, 3.771265786258514*^9}}],
283
284Cell[BoxData[
285 RowBox[{
286  RowBox[{"(*", " ", "Masses", " ", "*)"}], "\[IndentingNewLine]",
287  RowBox[{"RenormalizationRules", "[", "\"\<Masses\>\"", "]"}]}]], "Input",
288 CellChangeTimes->{{3.7451518713843946`*^9, 3.745151922162116*^9}, {
289   3.745164095936433*^9, 3.745164096200345*^9}, {3.745164380017013*^9,
290   3.745164380432691*^9}, {3.7451660736095877`*^9, 3.7451660750228863`*^9}, {
291   3.745213815312695*^9, 3.745213818856022*^9}, {3.74523151148385*^9,
292   3.7452315119531517`*^9}, {3.745303047472493*^9, 3.7453030478622513`*^9}, {
293   3.745303157313552*^9, 3.745303157889686*^9}, 3.745303406136817*^9, {
294   3.745574845804298*^9, 3.745574846145966*^9}, {3.745575407401681*^9,
295   3.7455754079116373`*^9}, {3.745578284288937*^9, 3.745578290465178*^9}, {
296   3.745582587034992*^9, 3.745582592977302*^9}, {3.7455826657703257`*^9,
297   3.7455826676585293`*^9}, {3.745582880576159*^9, 3.7455828928842697`*^9}, {
298   3.745583654476617*^9, 3.74558365461695*^9}, {3.7455837983158484`*^9,
299   3.745583798571127*^9}, 3.74564688658142*^9, {3.745647697834746*^9,
300   3.745647728311248*^9}, {3.745647806746401*^9, 3.745647810713956*^9}, {
301   3.7456480115435867`*^9, 3.7456481102599573`*^9}, {3.745663218654118*^9,
302   3.74566325795298*^9}, {3.745663711108547*^9, 3.745663711921837*^9}, {
303   3.745664612538536*^9, 3.7456646408908463`*^9}, 3.7456654013484077`*^9, {
304   3.745665650589293*^9, 3.7456656519560213`*^9}, {3.745665936303863*^9,
305   3.74566593728722*^9}, {3.7456688354046164`*^9, 3.7456688378674917`*^9}, {
306   3.745669024959618*^9, 3.745669030940509*^9}, {3.745669163428668*^9,
307   3.7456691987440863`*^9}, {3.745686812472076*^9, 3.74568681368325*^9}, {
308   3.745686880794334*^9, 3.74568688213133*^9}, {3.7456871775241747`*^9,
309   3.745687178863327*^9}, {3.745857725037485*^9, 3.745857735208394*^9}, {
310   3.7458586075203047`*^9, 3.745858613138812*^9}, {3.7458627781148767`*^9,
311   3.745862788627599*^9}, 3.745863379016315*^9, {3.745895772747778*^9,
312   3.7458957948440533`*^9}, {3.7458958293849087`*^9, 3.745895832308837*^9}, {
313   3.7463392933769217`*^9, 3.746339293953823*^9}, {3.746340958247019*^9,
314   3.746340980757111*^9}, {3.746446967049057*^9, 3.7464469672316236`*^9}, {
315   3.757309129229148*^9, 3.757309129577918*^9}, {3.757311847604082*^9,
316   3.757311850951976*^9}, {3.7573119060012207`*^9, 3.7573119133483257`*^9}, {
317   3.757400119485875*^9, 3.7574001368021793`*^9}, {3.757402948489932*^9,
318   3.757402960729457*^9}, {3.757403682464596*^9, 3.757403685895455*^9}, {
319   3.757406015321528*^9, 3.757406016888089*^9}, {3.75741526102887*^9,
320   3.757415265769416*^9}, {3.7574152965611763`*^9, 3.757415298314205*^9}, {
321   3.757416463439416*^9, 3.757416466010469*^9}, {3.757417500086125*^9,
322   3.757417523855345*^9}, {3.7574201860450373`*^9, 3.757420186670478*^9}, {
323   3.757448602343779*^9, 3.757448611302442*^9}, 3.757449335625184*^9, {
324   3.757449874383556*^9, 3.757449880587248*^9}, {3.757912547323681*^9,
325   3.757912548642529*^9}, {3.7580005002291107`*^9, 3.758000500760256*^9}, {
326   3.758005087249868*^9, 3.758005090807074*^9}, {3.758006435804907*^9,
327   3.758006444604958*^9}, {3.758161578806641*^9, 3.758161582999338*^9}, {
328   3.758164979635907*^9, 3.758164982381565*^9}, {3.758165013510289*^9,
329   3.758165029147532*^9}, {3.758165070730723*^9, 3.7581650858708963`*^9}, {
330   3.758165462675087*^9, 3.758165464347835*^9}, {3.758166197918043*^9,
331   3.7581662072139597`*^9}, {3.7581671637742367`*^9, 3.75816718202489*^9}, {
332   3.758177963109878*^9, 3.758177981847004*^9}, {3.758178134359646*^9,
333   3.758178159140918*^9}, {3.7582313015106*^9, 3.758231315264636*^9}, {
334   3.758248921123448*^9, 3.758248931762648*^9}, {3.758250734518733*^9,
335   3.758250783798354*^9}, {3.758250835228381*^9, 3.758250845550149*^9}, {
336   3.7582511056739693`*^9, 3.7582511090483103`*^9}, {3.75825131359096*^9,
337   3.758251325445635*^9}, {3.758251381280675*^9, 3.758251384116638*^9},
338   3.758251422640869*^9, {3.75825319671157*^9, 3.758253224083932*^9}, {
339   3.758257134438253*^9, 3.758257136759152*^9}, {3.758257209158524*^9,
340   3.7582572123266*^9}, {3.7582574013488417`*^9, 3.758257401749174*^9}, {
341   3.7582575392174883`*^9, 3.758257543960134*^9}, {3.758257710440304*^9,
342   3.758257731760181*^9}, {3.758257835439458*^9, 3.7582578864437447`*^9}, {
343   3.758258210240507*^9, 3.758258254261258*^9}, {3.758259393453458*^9,
344   3.7582594699190207`*^9}, {3.758259849130109*^9, 3.758259849496891*^9},
345   3.760938752558837*^9, {3.760938791593569*^9, 3.76093879283946*^9}, {
346   3.76094336376334*^9, 3.760943369090857*^9}, {3.760943562738961*^9,
347   3.760943588064082*^9}, 3.760944298912347*^9, {3.760952987831949*^9,
348   3.760952987979011*^9}, {3.760953226892173*^9, 3.760953236996063*^9}, {
349   3.76095348865825*^9, 3.760953489028057*^9}, {3.763457602937347*^9,
350   3.763457611310288*^9}, 3.763982907134116*^9, {3.763983088407049*^9,
351   3.763983089444828*^9}, {3.763983130719948*^9, 3.763983131117784*^9}, {
352   3.763983964034717*^9, 3.763983965670892*^9}, {3.763987056564139*^9,
353   3.763987059126238*^9}, {3.764005252103383*^9, 3.764005252336042*^9}, {
354   3.7640062913058977`*^9, 3.7640062939055023`*^9}, {3.7640073431045647`*^9,
355   3.764007410995483*^9}, {3.76452231518419*^9, 3.764522315471777*^9},
356   3.764522354075388*^9, {3.7645232675315247`*^9, 3.764523275689969*^9}, {
357   3.7645238003142233`*^9, 3.764523801367539*^9}, {3.764523872479026*^9,
358   3.764523881589555*^9}, {3.771263472262384*^9, 3.771263478300324*^9}, {
359   3.771263615037691*^9, 3.77126364922724*^9}, {3.771264545348506*^9,
360   3.771264606885642*^9}, 3.771264660161273*^9, 3.771264747418335*^9, {
361   3.771265717192766*^9, 3.771265723035942*^9}}],
362
363Cell[BoxData[
364 RowBox[{
365  RowBox[{"(*", " ",
366   RowBox[{"External", " ", "parameters"}], " ", "*)"}],
367  "\[IndentingNewLine]",
368  RowBox[{"RenormalizationRules", "[", "\"\<Externals\>\"", "]"}]}]], "Input",\
369
370 CellChangeTimes->{{3.771265730383052*^9, 3.771265734792964*^9}}],
371
372Cell[BoxData[
373 RowBox[{
374  RowBox[{"(*", " ",
375   RowBox[{"Internal", " ", "parameters"}], " ", "*)"}],
376  "\[IndentingNewLine]",
377  RowBox[{"RenormalizationRules", "[", "\"\<Internals\>\"", "]"}]}]], "Input",\
378
379 CellChangeTimes->{{3.758258382884212*^9, 3.758258384402162*^9}, {
380  3.7712657927013063`*^9, 3.77126579615121*^9}}],
381
382Cell[BoxData[
383 RowBox[{
384  RowBox[{"(*", " ", "Fields", " ", "*)"}], "\[IndentingNewLine]",
385  RowBox[{"RenormalizationRules", "[", "\"\<Fields\>\"", "]"}]}]], "Input"]
386}, Open  ]],
387
388Cell[CellGroupData[{
389
390Cell["Writing the FA model", "Subsubsection",
391 CellChangeTimes->{{3.7796363600275784`*^9, 3.779636361893107*^9}}],
392
393Cell[BoxData[{
394 RowBox[{
395  RowBox[{"SetDirectory", "[",
396   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
397  ";"}], "\[IndentingNewLine]",
398 RowBox[{
399  RowBox[{"WriteFeynArtsOutput", "[",
400   RowBox[{"LReno", ",",
401    RowBox[{"Output", "\[Rule]", "\"\<typeII_v1_1_FA\>\""}]}], "]"}],
402  ";"}]}], "Input",
403 CellChangeTimes->{{3.779597577798992*^9, 3.779597584199922*^9}, {
404  3.779597680344357*^9, 3.779597685565913*^9}, {3.77959779313251*^9,
405  3.779597803082788*^9}, {3.779632511026588*^9, 3.77963252396973*^9}, {
406  3.779632589533596*^9, 3.779632594866434*^9}, {3.780944772665889*^9,
407  3.780944773643894*^9}, {3.7809989929665813`*^9, 3.780998995191955*^9}}]
408}, Open  ]],
409
410Cell[CellGroupData[{
411
412Cell["NLOCT run", "Subsubsection",
413 CellChangeTimes->{{3.779636351041431*^9, 3.779636352006412*^9}}],
414
415Cell[BoxData[
416 RowBox[{
417  RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
418 CellChangeTimes->{{3.779636370018443*^9, 3.779636371598173*^9}}],
419
420Cell[BoxData[{
421 RowBox[{"<<", " ", "FeynArts`"}], "\[IndentingNewLine]",
422 RowBox[{
423  RowBox[{
424  "SetDirectory", "[",
425   "\"\<~/Work/tools/FeynRules/trunk/feynrules-development\>\"", "]"}],
426  ";"}], "\[IndentingNewLine]",
427 RowBox[{"<<", " ", "NLOCT`"}], "\[IndentingNewLine]",
428 RowBox[{
429  RowBox[{"SetDirectory", "[",
430   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
431  ";"}], "\[IndentingNewLine]",
432 RowBox[{
433  RowBox[{"WriteCT", "[",
434   RowBox[{
435   "\"\<typeII_FA/typeII_v1_1_FA\>\"", ",",
436    "\"\<typeII_FA/typeII_v1_1_FA\>\"", ",", " ",
437    RowBox[{"Output", "->", "\"\<typeII_4FNS_v1_1\>\""}], ",",
438    "\[IndentingNewLine]",
439    RowBox[{"LabelInternal", "\[Rule]", "True"}], ",",
440    RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
441    RowBox[{"KeptIndices", "\[Rule]",
442     RowBox[{"{", "}"}]}], ",",
443    RowBox[{"ZeroMom", "\[Rule]",
444     RowBox[{"{",
445      RowBox[{"{",
446       RowBox[{"aS", ",",
447        RowBox[{"{",
448         RowBox[{
449          RowBox[{"F", "[", "7", "]"}], ",",
450          RowBox[{"V", "[", "4", "]"}], ",",
451          RowBox[{"-",
452           RowBox[{"F", "[", "7", "]"}]}]}], "}"}], ",", "0"}], "}"}],
453      "}"}]}], ",", "\[IndentingNewLine]",
454    RowBox[{"Assumptions", "->",
455     RowBox[{"{",
456      RowBox[{
457       RowBox[{"MT", ">", "0"}], ",",
458       RowBox[{"MB", ">", "0"}]}], "}"}]}], ",", "\[IndentingNewLine]",
459    RowBox[{"CTparameters", "\[Rule]", "True"}], ",",
460    RowBox[{"MaxDim", "\[Rule]", "4"}], ",",
461    RowBox[{"EvenOnly", "\[Rule]", "False"}], ",",
462    RowBox[{"IsFeynmanGauge", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",
463 CellChangeTimes->{{3.7796363829653883`*^9, 3.779636429912023*^9}, {
464  3.780999051176557*^9, 3.780999062807684*^9}, {3.7809991888797283`*^9,
465  3.7809991906985292`*^9}}]
466}, Open  ]],
467
468Cell[CellGroupData[{
469
470Cell["NLO UFO", "Subsubsection",
471 CellChangeTimes->{{3.7796397724027576`*^9, 3.7796397744779*^9}}],
472
473Cell[BoxData[
474 RowBox[{"Quit", "[", "]"}]], "Input",
475 CellChangeTimes->{{3.779639826390726*^9, 3.7796398270234957`*^9}}],
476
477Cell[BoxData[{
478 RowBox[{
479  RowBox[{"$FeynRulesPath", "=",
480   RowBox[{
481   "SetDirectory", "[",
482    "\"\<~/Work/tools/FeynRules/trunk/feynrules-development\>\"", "]"}]}],
483  ";"}], "\[IndentingNewLine]",
484 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
485 RowBox[{
486  RowBox[{"SetDirectory", "[",
487   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
488  ";"}], "\[IndentingNewLine]",
489 RowBox[{
490  RowBox[{"LoadModel", "[",
491   RowBox[{"\"\<sm_reduced.fr\>\"", ",", "\"\<type_ii.fr\>\""}], "]"}],
492  ";"}], "\[IndentingNewLine]",
493 RowBox[{
494  RowBox[{"LoadRestriction", "[",
495   RowBox[{
496    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
497    "\"\</Massless_4f.rst\>\""}], "]"}], ";"}], "\[IndentingNewLine]",
498 RowBox[{
499  RowBox[{"LoadRestriction", "[",
500   RowBox[{
501    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
502    "\"\</DiagonalCKM.rst\>\""}], "]"}], ";"}], "\[IndentingNewLine]",
503 RowBox[{
504  RowBox[{"Get", "[", "\"\<typeII_4FNS_v1_1.nlo\>\"", "]"}],
505  ";"}], "\[IndentingNewLine]",
506 RowBox[{
507  RowBox[{"If", "[",
508   RowBox[{
509    RowBox[{"!",
510     RowBox[{"FileExistsQ", "[", "\"\<lagrangian.dat\>\"", "]"}]}], ",",
511    "\[IndentingNewLine]",
512    RowBox[{
513     RowBox[{"lags", "=", " ",
514      RowBox[{"ExpandIndices", "/@", "LType2"}]}], ";", "\[IndentingNewLine]",
515     RowBox[{
516      RowBox[{"Definition", "[", "lags", "]"}], " ", ">>", " ",
517      "lagrangian.dat"}]}], ",", "\[IndentingNewLine]",
518    RowBox[{
519     RowBox[{"<<", "lagrangian.dat"}], ";"}]}], "\[IndentingNewLine]", "]"}],
520  ";"}], "\[IndentingNewLine]"}], "Input",
521 CellChangeTimes->{{3.7796397898851013`*^9, 3.779639824714785*^9}, {
522  3.780999092284133*^9, 3.780999094749385*^9}}],
523
524Cell[BoxData[
525 RowBox[{
526  RowBox[{"WriteUFO", "[",
527   RowBox[{"lags", ",",
528    RowBox[{"UVCounterterms", "\[Rule]",
529     RowBox[{"(",
530      RowBox[{"UV$vertlist", "/.",
531       RowBox[{"{",
532        RowBox[{
533         RowBox[{"FR$IR", "\[Rule]", "1"}], ",",
534         RowBox[{"FR$IRLog", "\[Rule]", "1"}]}], "}"}]}], ")"}]}], " ", ",",
535    RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
536    RowBox[{"CTParameters", "\[Rule]",
537     RowBox[{"(",
538      RowBox[{"FR$CTparam", "/.",
539       RowBox[{"{",
540        RowBox[{
541         RowBox[{"FR$IR", "\[Rule]", "1"}], ",",
542         RowBox[{"FR$IRLog", "\[Rule]", "1"}]}], "}"}]}], ")"}]}], ",",
543    RowBox[{"Output", "\[Rule]", "\"\<TypeII_NLO_v_1_1_UFO\>\""}]}], "]"}],
544  ";"}]], "Input",
545 CellChangeTimes->{{3.779639836439653*^9, 3.779639854005334*^9}, {
546  3.7809990847792883`*^9, 3.780999086952511*^9}}]
547}, Open  ]]
548}, Open  ]],
549
550Cell[CellGroupData[{
551
552Cell["LO ", "Section",
553 CellChangeTimes->{{3.7796397724027576`*^9, 3.7796397744779*^9},
554   3.7809756787787046`*^9, 3.780999031575129*^9}],
555
556Cell[BoxData[
557 RowBox[{"Quit", "[", "]"}]], "Input",
558 CellChangeTimes->{{3.779639826390726*^9, 3.7796398270234957`*^9}}],
559
560Cell[BoxData[{
561 RowBox[{
562  RowBox[{"$FeynRulesPath", "=",
563   RowBox[{
564   "SetDirectory", "[",
565    "\"\<~/Work/tools/FeynRules/trunk/feynrules-development\>\"", "]"}]}],
566  ";"}], "\[IndentingNewLine]",
567 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
568 RowBox[{
569  RowBox[{"SetDirectory", "[",
570   RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
571  ";"}], "\[IndentingNewLine]",
572 RowBox[{
573  RowBox[{"LoadModel", "[",
574   RowBox[{"\"\<sm_reduced.fr\>\"", ",", "\"\<type_ii.fr\>\""}], "]"}],
575  ";"}], "\[IndentingNewLine]",
576 RowBox[{
577  RowBox[{"LoadRestriction", "[",
578   RowBox[{
579    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
580    "\"\</Massless_4f.rst\>\""}], "]"}], ";"}], "\[IndentingNewLine]",
581 RowBox[{
582  RowBox[{"LoadRestriction", "[",
583   RowBox[{
584    RowBox[{"NotebookDirectory", "[", "]"}], "<>",
585    "\"\</DiagonalCKM.rst\>\""}], "]"}], ";"}], "\[IndentingNewLine]",
586 RowBox[{
587  RowBox[{"If", "[",
588   RowBox[{
589    RowBox[{"!",
590     RowBox[{"FileExistsQ", "[", "\"\<lagrangian.dat\>\"", "]"}]}], ",",
591    "\[IndentingNewLine]",
592    RowBox[{
593     RowBox[{"lags", "=", " ",
594      RowBox[{"ExpandIndices", "/@", "LType2"}]}], ";", "\[IndentingNewLine]",
595     RowBox[{
596      RowBox[{"Definition", "[", "lags", "]"}], " ", ">>", " ",
597      "lagrangian.dat"}]}], ",", "\[IndentingNewLine]",
598    RowBox[{
599     RowBox[{"<<", "lagrangian.dat"}], ";"}]}], "\[IndentingNewLine]", "]"}],
600  ";"}], "\[IndentingNewLine]",
601 RowBox[{
602  RowBox[{"WriteUFO", "[",
603   RowBox[{"lags", ",",
604    RowBox[{"Output", "\[Rule]", "\"\<TypeII_LO_v1_1_UFO\>\""}]}], "]"}],
605  ";"}], "\[IndentingNewLine]"}], "Input",
606 CellChangeTimes->{{3.7796397898851013`*^9, 3.779639824714785*^9}, {
607  3.780975707166958*^9, 3.780975707338049*^9}, {3.780999012384077*^9,
608  3.780999012577855*^9}}]
609}, Open  ]]
610},
611WindowSize->{1292, 1109},
612WindowMargins->{{Automatic, 3}, {Automatic, 31}},
613ShowSelection->True,
614FrontEndVersion->"10.4 for Mac OS X x86 (32-bit, 64-bit Kernel) (April 11, \
6152016)",
616StyleDefinitions->"Default.nb"
617]
618(* End of Notebook Content *)
619
620(* Internal cache information *)
621(*CellTagsOutline
622CellTagsIndex->{}
623*)
624(*CellTagsIndex
625CellTagsIndex->{}
626*)
627(*NotebookFileOutline
628Notebook[{
629Cell[545, 20, 220, 4, 28, "Input"],
630Cell[CellGroupData[{
631Cell[790, 28, 188, 3, 64, "Section"],
632Cell[981, 33, 2503, 46, 114, "Input"]
633}, Open  ]],
634Cell[CellGroupData[{
635Cell[3521, 84, 90, 1, 64, "Section"],
636Cell[CellGroupData[{
637Cell[3636, 89, 147, 2, 35, "Subsubsection"],
638Cell[3786, 93, 628, 15, 97, "Input"],
639Cell[4417, 110, 1037, 26, 114, "Input"],
640Cell[5457, 138, 363, 9, 46, "Input"],
641Cell[5823, 149, 791, 17, 46, "Input"]
642}, Open  ]],
643Cell[CellGroupData[{
644Cell[6651, 171, 207, 3, 35, "Subsubsection"],
645Cell[6861, 176, 1053, 20, 63, "Input"],
646Cell[7917, 198, 414, 10, 63, "Input"],
647Cell[8334, 210, 1490, 40, 114, "Input"]
648}, Open  ]],
649Cell[CellGroupData[{
650Cell[9861, 255, 133, 1, 35, "Subsubsection"],
651Cell[9997, 258, 445, 8, 28, "Input"]
652}, Open  ]],
653Cell[CellGroupData[{
654Cell[10479, 271, 110, 1, 35, "Subsubsection"],
655Cell[10592, 274, 332, 7, 46, "Input"],
656Cell[10927, 283, 5597, 77, 46, "Input"],
657Cell[16527, 362, 276, 7, 46, "Input"],
658Cell[16806, 371, 326, 8, 46, "Input"],
659Cell[17135, 381, 167, 3, 46, "Input"]
660}, Open  ]],
661Cell[CellGroupData[{
662Cell[17339, 389, 113, 1, 35, "Subsubsection"],
663Cell[17455, 392, 661, 14, 46, "Input"]
664}, Open  ]],
665Cell[CellGroupData[{
666Cell[18153, 411, 100, 1, 35, "Subsubsection"],
667Cell[18256, 414, 136, 3, 28, "Input"],
668Cell[18395, 419, 1776, 45, 148, "Input"]
669}, Open  ]],
670Cell[CellGroupData[{
671Cell[20208, 469, 98, 1, 35, "Subsubsection"],
672Cell[20309, 472, 120, 2, 28, "Input"],
673Cell[20432, 476, 1675, 45, 233, "Input"],
674Cell[22110, 523, 872, 22, 46, "Input"]
675}, Open  ]]
676}, Open  ]],
677Cell[CellGroupData[{
678Cell[23031, 551, 138, 2, 64, "Section"],
679Cell[23172, 555, 120, 2, 28, "Input"],
680Cell[23295, 559, 1790, 48, 233, "Input"]
681}, Open  ]]
682}
683]
684*)
685