SMWeinberg: smWeinberg_NLO_public.nb

File smWeinberg_NLO_public.nb, 25.8 KB (added by Richard Ruiz, 4 years ago)

Mathematica notebook for UFO generation

Line 
1(* Content-type: application/mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 6.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 145, 7]
13NotebookDataLength[ 26276, 621]
14NotebookOptionsPosition[ 22119, 553]
15NotebookOutlinePosition[ 22555, 570]
16CellTagsIndexPosition[ 22512, 567]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"(*", " ",
23 RowBox[{
24 RowBox[{"R", ".", " ", "Ruiz"}], ";", " ",
25 RowBox[{"2020", " ", "December", " ", "18"}]}], " ", "*)"}]], "Input",
26 CellChangeTimes->{{3.665142881510211*^9, 3.665142884049021*^9}, {
27 3.665142942303646*^9, 3.6651429529879627`*^9}, {3.8002851169494963`*^9,
28 3.8002851236251*^9}, {3.814537185371134*^9, 3.814537193009823*^9}, {
29 3.817283510000004*^9,
30 3.817283510271615*^9}},ExpressionUUID->"cee4e897-cbc1-4bf1-a304-\
3144bbfcd5aed3"],
32
33Cell[BoxData[
34 RowBox[{
35 RowBox[{"(*", " ",
36 RowBox[{
37 RowBox[{"To", " ", "generate", " ", "LO", " ",
38 RowBox[{"UFO", ":", " ",
39 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
40 RowBox[{"load", " ", "model"}], ",", " ",
41 RowBox[{"write", " ", "UFO"}]}], " ", "*)"}], "\[IndentingNewLine]",
42 RowBox[{"(*", " ",
43 RowBox[{
44 RowBox[{"To", " ", "generate", " ", "NLO", " ",
45 RowBox[{"UFO", ":", " ",
46 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
47 RowBox[{"load", " ", "model"}], ",", " ",
48 RowBox[{
49 RowBox[{"extract", " ", "on"}], "-",
50 RowBox[{"shell", " ", "renormalization", " ", "terms"}]}], ",", " ",
51 RowBox[{"write", " ", "CT", " ", "terms"}], ",", " ",
52 RowBox[{"reset", " ", "paths"}], ",", " ",
53 RowBox[{"reload", " ", "model"}], ",", " ",
54 RowBox[{"write", " ", "everything", " ", "to", " ",
55 RowBox[{"file", "."}]}]}], " ", "*)"}]}]], "Input",
56 CellChangeTimes->{{3.6648662671878157`*^9, 3.664866418502534*^9}, {
57 3.664866852892112*^9,
58 3.664866879416404*^9}},ExpressionUUID->"4149d1d2-f5f9-4692-ab77-\
59ed3366c4e0f7"],
60
61Cell[BoxData[
62 RowBox[{
63 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
64 CellChangeTimes->{{3.800203881213688*^9, 3.800203885091346*^9}, {
65 3.814192661966413*^9, 3.81419266729278*^9}, {3.814196498380603*^9,
66 3.814196499409697*^9}, {3.814196533345257*^9, 3.814196534014574*^9}, {
67 3.814537179057211*^9, 3.814537179988576*^9}, {3.817278556475902*^9,
68 3.817278557221746*^9}},
69 CellLabel->"In[5]:=",ExpressionUUID->"69969515-0b1d-48d5-8594-bd13de83a811"],
70
71Cell[CellGroupData[{
72
73Cell["FeynRules Setup", "Title",
74 CellChangeTimes->{{3.596207285657744*^9, 3.596207286276866*^9}, {
75 3.596207330352247*^9, 3.59620733250506*^9}, {3.655035967648223*^9,
76 3.65503597269068*^9}},ExpressionUUID->"36a8985f-3b71-4bb1-acae-\
776fbe5493c081"],
78
79Cell[BoxData[{
80 RowBox[{
81 RowBox[{"$CurrentPath", "=",
82 RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
83 RowBox[{
84 RowBox[{"$FeynRulesPath", "=",
85 RowBox[{
86 "SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}]}],
87 ";"}], "\[IndentingNewLine]",
88 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
89 RowBox[{
90 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
91 CellChangeTimes->{{3.41265862251538*^9, 3.412658649947229*^9}, {
92 3.423415585782702*^9, 3.423415597189939*^9}, {3.4234163173467493`*^9,
93 3.4234163227881193`*^9}, {3.572102427116931*^9, 3.572102456206238*^9}, {
94 3.572280529035656*^9, 3.572280529468644*^9}, {3.572366766720894*^9,
95 3.572366773370627*^9}, {3.596207212826934*^9, 3.5962072447909603`*^9}, {
96 3.596914045279478*^9, 3.5969140472383204`*^9}, 3.596914102616413*^9,
97 3.596915921742405*^9, {3.596915961175353*^9, 3.59691596259039*^9}, {
98 3.60744636894991*^9, 3.607446377507101*^9}, {3.6364757082312813`*^9,
99 3.636475708315627*^9}, {3.655035981353272*^9, 3.655035982244761*^9}, {
100 3.6550360677717857`*^9, 3.655036081259776*^9}, {3.655036117419112*^9,
101 3.655036143336364*^9}, {3.655037954145546*^9, 3.6550379644813023`*^9}, {
102 3.6605209855658293`*^9, 3.6605209855964212`*^9}, {3.6605772969976788`*^9,
103 3.660577300733718*^9}, {3.6621928886588717`*^9, 3.662192889375732*^9}, {
104 3.6628197402754908`*^9, 3.662819763199905*^9}, 3.800201454262025*^9},
105 CellLabel->"In[1]:=",ExpressionUUID->"2b566197-5488-4f54-9e47-2d15cdf0dd38"],
106
107Cell[BoxData[
108 RowBox[{"(*",
109 RowBox[{
110 "Load", " ", "default", " ", "SM", " ", "file", " ", "v1", ".4", ".7", " ",
111 "with", " ", "Weinberg", " ",
112 RowBox[{"extension", "."}]}], "*)"}]], "Input",
113 CellChangeTimes->{{3.6550381551475554`*^9, 3.655038155711425*^9}, {
114 3.6622141225281467`*^9, 3.662214139672209*^9}, {3.730294820431093*^9,
115 3.730294820982045*^9}, {3.800287658406439*^9, 3.800287674074195*^9}, {
116 3.814539926940324*^9,
117 3.8145399303602047`*^9}},ExpressionUUID->"d74de1e8-93f7-4afe-980e-\
1188f4e6a416a1f"],
119
120Cell[BoxData[{
121 RowBox[{
122 RowBox[{"LoadModel", "[",
123 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<weinberg_NLO.fr\>\""}], "]"}],
124 ";"}], "\[IndentingNewLine]",
125 RowBox[{
126 RowBox[{"LoadRestriction", "[",
127 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<Massless_4f.rst\>\""}], "]"}],
128 ";", " ",
129 RowBox[{"(*", " ",
130 RowBox[{
131 RowBox[{"this", " ", "will", " ", "setup", " ", "nf"}], "=", "5"}], " ",
132 "*)"}], "\[IndentingNewLine]",
133 RowBox[{"(*",
134 RowBox[{
135 RowBox[{"LoadRestriction", "[",
136 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<Massless_5f.rst\>\""}],
137 "]"}], ";"}], "*)"}], "\[IndentingNewLine]",
138 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}], "Input",
139 CellChangeTimes->{
140 3.800287989386373*^9, {3.800288102386517*^9, 3.800288107761572*^9}, {
141 3.814192695990925*^9, 3.814192696173706*^9}, {3.814199115854692*^9,
142 3.8141991169590063`*^9}, {3.814199841138897*^9, 3.814199842278466*^9}, {
143 3.814537872771997*^9, 3.814537879456297*^9}, {3.814538589414817*^9,
144 3.8145386022711287`*^9}, {3.814539911494363*^9, 3.814539912985259*^9}, {
145 3.814543467392374*^9, 3.814543467577437*^9}, {3.8145445432764997`*^9,
146 3.814544543350997*^9}, {3.817279878725089*^9, 3.817279878887182*^9}, {
147 3.8172835245449867`*^9,
148 3.817283558366296*^9}},ExpressionUUID->"a4b9bd52-5e8a-4f48-8fba-\
149852bc4d59da3"],
150
151Cell[CellGroupData[{
152
153Cell["Lagrangian Terms", "Section",
154 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
155 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
156 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
157 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
158 3.6550378926486187`*^9}, {3.662214296374611*^9,
159 3.662214296674329*^9}},ExpressionUUID->"430f6c12-ab98-447b-8fa4-\
1608a13c875b6c3"],
161
162Cell[BoxData["LD5"], "Input",
163 CellChangeTimes->{{3.664866736738257*^9, 3.664866737292602*^9}, {
164 3.814539883683156*^9,
165 3.814539887169508*^9}},ExpressionUUID->"01776d56-05d4-475b-b025-\
1662f0ebd0c8fc8"],
167
168Cell[BoxData["LNKin"], "Input",
169 CellChangeTimes->{{3.5723482442775097`*^9, 3.572348252612296*^9}, {
170 3.596212546796714*^9, 3.596212547440735*^9}, {3.655037865036394*^9,
171 3.6550378653420067`*^9}, {3.6622141995703773`*^9,
172 3.6622142141549997`*^9}},ExpressionUUID->"4b15b6dc-b821-4c12-bd21-\
173ceb678b9a346"],
174
175Cell[BoxData["LNCC"], "Input",
176 CellChangeTimes->{{3.6622142176950073`*^9, 3.6622142303960037`*^9}, {
177 3.8141990155618563`*^9, 3.814199025276359*^9}},
178 CellLabel->"In[8]:=",ExpressionUUID->"cc1f0f3b-1d9f-45c8-a8b8-49b2fdacb74f"],
179
180Cell[BoxData["LNNC"], "Input",
181 CellChangeTimes->{{3.81453989606852*^9,
182 3.814539901449139*^9}},ExpressionUUID->"f0517d40-71c0-4541-a5c1-\
183c1cbc531173e"],
184
185Cell[BoxData["LNHX"], "Input",
186 CellChangeTimes->{{3.6622142357666893`*^9, 3.6622142376786947`*^9},
187 3.662826298532291*^9,
188 3.8141999256497183`*^9},ExpressionUUID->"e4b33362-cb94-4ba2-a771-\
189443401bf60bf"],
190
191Cell[BoxData["LNGX"], "Input",
192 CellChangeTimes->{{3.6622142239436293`*^9, 3.6622142242517548`*^9},
193 3.66282630175613*^9,
194 3.8141999279705963`*^9},ExpressionUUID->"f71c29b1-311a-4c25-824e-\
195bfe62f098b98"]
196}, Open ]],
197
198Cell[CellGroupData[{
199
200Cell["Sanity Checks for New Lagrangian Terms", "Section",
201 CellChangeTimes->{{3.4234156496121893`*^9, 3.423415657027569*^9}, {
202 3.655037883884348*^9,
203 3.655037906883815*^9}},ExpressionUUID->"963bd5da-1517-491b-8768-\
2043e81fb288ca2"],
205
206Cell[BoxData[{
207 RowBox[{"CheckHermiticity", "[", "LD5", "]"}], "\[IndentingNewLine]",
208 RowBox[{"CheckMassSpectrum", "[", "LD5", "]"}], "\[IndentingNewLine]",
209 RowBox[{"CheckKineticTermNormalisation", "[", "LD5", "]"}]}], "Input",
210 CellChangeTimes->{
211 3.607446407230465*^9, {3.800204083442677*^9, 3.8002040843202953`*^9}, {
212 3.814539858341589*^9, 3.814539866624523*^9}, {3.817278596818945*^9,
213 3.817278648403672*^9}},
214 CellLabel->"In[8]:=",ExpressionUUID->"ad90ef2a-90e9-4499-9f66-242849c1984c"]
215}, Open ]],
216
217Cell[CellGroupData[{
218
219Cell["\<\
220New Feynman Rules for Effective Field Theory (Weinberg, d=5 operator) \
221Interactions\
222\>", "Section",
223 CellChangeTimes->{{3.411910637332127*^9, 3.41191064360071*^9}, {
224 3.655037529059973*^9, 3.655037541428491*^9}, {3.6550379116719217`*^9,
225 3.655037913954034*^9}, {3.800256492276526*^9, 3.800256495041098*^9}, {
226 3.8145446909104767`*^9, 3.814544708256122*^9},
227 3.817283587962694*^9},ExpressionUUID->"4f7e70e8-ac63-452b-bc39-\
2282e283d37d5b9"],
229
230Cell[BoxData[
231 RowBox[{"(*", " ",
232 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
233 RowBox[{"appear", ":", " ",
234 RowBox[{
235 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
236\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
237 "Warning"}], ":", " ",
238 RowBox[{
239 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
240 " ",
241 RowBox[{"encountered", "!"}],
242 "\"\<;\nThis is expected as Lepton Number is violated and FeynRules \
243checks for this.\>"}]}]}], "*)"}]], "Input",
244 CellChangeTimes->{{3.6550374096249533`*^9, 3.655037447663039*^9},
245 3.655037551040966*^9, {3.656940976048616*^9,
246 3.6569409786068287`*^9}},ExpressionUUID->"f33d5026-da81-408c-99ad-\
247a2da2ac609d7"],
248
249Cell[BoxData[
250 RowBox[{"FeynmanRules", "[", "LD5", "]"}]], "Input",ExpressionUUID->"0fd184c5-e009-4777-a8c1-f1520332de5e"]
251}, Open ]],
252
253Cell[CellGroupData[{
254
255Cell["UFO Output", "Section",
256 CellChangeTimes->{
257 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
258 3.59620768137986*^9,
259 3.596207684212212*^9}},ExpressionUUID->"47f009e8-0603-4d37-95af-\
260884018211353"],
261
262Cell[BoxData[
263 RowBox[{"(*", " ",
264 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
265 RowBox[{"appear", ":", " ",
266 RowBox[{
267 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
268\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
269 "Warning"}], ":", " ",
270 RowBox[{
271 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
272 " ",
273 RowBox[{"encountered", "!"}],
274 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
275checks for this.\>"}]}]}], "*)"}]], "Input",
276 CellChangeTimes->{{3.656940935403554*^9, 3.656940952617087*^9}, {
277 3.6648662530247602`*^9,
278 3.664866256652151*^9}},ExpressionUUID->"95ec60e1-8ba3-42ae-912b-\
279c6b75ece5652"],
280
281Cell[BoxData[
282 RowBox[{"(*",
283 RowBox[{
284 RowBox[{"WriteUFO", "[", "LSM", "]"}], ";"}], "*)"}]], "Input",
285 CellChangeTimes->{{3.814539824786786*^9,
286 3.814539832318253*^9}},ExpressionUUID->"0a6453a6-cc72-4868-bf62-\
287a55b89adae23"],
288
289Cell[BoxData[{
290 RowBox[{
291 RowBox[{"LWeinberg", " ", "=", " ",
292 RowBox[{"ExpandIndices", "/@",
293 RowBox[{"Expand", "[",
294 RowBox[{"LSM", "+", "LD5"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
295 RowBox[{
296 RowBox[{"WriteUFO", "[",
297 RowBox[{"LWeinberg", ",",
298 RowBox[{"Output", "\[Rule]", "\"\<SMWeinbergXLO\>\""}]}], "]"}],
299 ";"}], "\[IndentingNewLine]",
300 RowBox[{
301 RowBox[{"(*",
302 RowBox[{
303 RowBox[{"WriteUFO", "[",
304 RowBox[{"LWeinberg", ",",
305 RowBox[{"Output", "\[Rule]", "\"\<SMWeinbergXLO_4FS\>\""}]}], "]"}],
306 ";"}], "*)"}]}]}], "Input",
307 CellChangeTimes->{{3.817283600637013*^9,
308 3.817283601298911*^9}},ExpressionUUID->"dbed1df9-509a-4a3f-981c-\
3099d0520eb5336"]
310}, Open ]],
311
312Cell[CellGroupData[{
313
314Cell["On-Shell Renormalization with FeynArts", "Section",
315 CellChangeTimes->{
316 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
317 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199318805483*^9,
318 3.6621993287719097`*^9}, {3.662216149623213*^9,
319 3.662216170908692*^9}},ExpressionUUID->"f35c5a82-2f44-4096-9399-\
320e82f5b97ff8d"],
321
322Cell[BoxData[
323 RowBox[{"(*", " ",
324 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
325 RowBox[{"appear", ":", " ",
326 RowBox[{
327 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
328\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
329 "Warning"}], ":", " ",
330 RowBox[{
331 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
332 " ",
333 RowBox[{"encountered", "!"}],
334 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
335checks for this.\>"}]}]}], "*)"}]], "Input",ExpressionUUID->"eafcd19a-c48f-\
33646e1-8302-30202d5af7c6"],
337
338Cell[BoxData[{
339 RowBox[{
340 RowBox[{"LRen", " ", "=", " ",
341 RowBox[{"OnShellRenormalization", "[",
342 RowBox[{"LWeinberg", ",",
343 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
344 RowBox[{"FlavorMixing", "\[Rule]", " ", "False"}], ",",
345 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}]}],
346 ";"}], "\[IndentingNewLine]",
347 RowBox[{"SetDirectory", "[",
348 "\"\<~/Scripts/FeynArts/FeynArts-3.11/Models\>\"",
349 "]"}], "\[IndentingNewLine]",
350 RowBox[{
351 RowBox[{"WriteFeynArtsOutput", "[",
352 RowBox[{"LRen", ",",
353 RowBox[{"FlavorExpand", "\[Rule]", "True"}], ",",
354 RowBox[{"Output", "\[Rule]", "\"\<SMWeinberg_FA\>\""}]}], "]"}],
355 ";"}]}], "Input",ExpressionUUID->"241b18a3-6df1-40c1-9540-ebe4d334b067"]
356}, Open ]],
357
358Cell[CellGroupData[{
359
360Cell["Write Counter Terms to BSM.nlo", "Section",
361 CellChangeTimes->{
362 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
363 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199461128889*^9,
364 3.662199467953211*^9}, {3.662216541485611*^9,
365 3.662216544705964*^9}},ExpressionUUID->"637e3231-e751-489a-ad59-\
366470d98c11b8b"],
367
368Cell[BoxData[
369 RowBox[{"(*", " ",
370 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
371 "*)"}]], "Input",
372 CellChangeTimes->{{3.664866502336069*^9,
373 3.664866536918194*^9}},ExpressionUUID->"558b9e27-1057-416f-a429-\
374231889bb6a86"],
375
376Cell[BoxData[
377 RowBox[{
378 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
379 CellLabel->"In[1]:=",ExpressionUUID->"3640fc82-001c-4a97-8551-70d99a9b683d"],
380
381Cell[BoxData[{
382 RowBox[{
383 RowBox[{
384 "SetDirectory", "[", "\"\<~/Scripts/FeynArts/FeynArts-3.11/\>\"", "]"}],
385 ";"}], "\[IndentingNewLine]",
386 RowBox[{"<<", "FeynArts`"}], "\[IndentingNewLine]",
387 RowBox[{
388 RowBox[{"SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}],
389 ";"}], "\[IndentingNewLine]",
390 RowBox[{"<<", "NLOCT`"}], "\[IndentingNewLine]",
391 RowBox[{
392 RowBox[{"$CurrentPath", "=",
393 RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
394 RowBox[{
395 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
396 CellChangeTimes->{{3.6605244352671556`*^9, 3.660524496529204*^9}, {
397 3.660524555733032*^9, 3.660524556119355*^9}, 3.660528429455056*^9, {
398 3.6605286020387897`*^9, 3.660528602100609*^9}, {3.660530319063929*^9,
399 3.660530319160571*^9}, {3.660530366154963*^9, 3.66053037265125*^9}, {
400 3.660567620251087*^9, 3.6605676243565598`*^9}, {3.660568912072307*^9,
401 3.660568912133321*^9}, {3.66056999189439*^9, 3.660569993076742*^9}, {
402 3.660570552697504*^9, 3.660570554442918*^9}, {3.6605786705389338`*^9,
403 3.660578696714662*^9}, {3.660579471210005*^9, 3.660579471258161*^9}, {
404 3.660579558812049*^9, 3.6605795588460417`*^9}, {3.6605809812710743`*^9,
405 3.660580981737924*^9}, {3.662198567340746*^9, 3.662198568029338*^9}, {
406 3.664866430517086*^9, 3.664866431258605*^9}, {3.664866469872073*^9,
407 3.664866472738167*^9}, {3.730294872498015*^9, 3.730294872583467*^9},
408 3.800257031453598*^9, {3.800257186336883*^9, 3.8002571995691557`*^9}, {
409 3.80025737377853*^9, 3.800257382190775*^9}, {3.800258132600809*^9,
410 3.800258156657604*^9}, {3.800258308009911*^9, 3.800258308377743*^9}, {
411 3.800260676659759*^9, 3.800260677117632*^9}},
412 CellLabel->"In[1]:=",ExpressionUUID->"cb2da490-5600-44da-b7e4-f7f288533aed"],
413
414Cell[BoxData[
415 RowBox[{"(*",
416 RowBox[{
417 RowBox[{"Note", ":", " ",
418 RowBox[{
419 "Grab", " ", "some", " ", "coffee", " ", "or", " ", "go", " ", "for", " ",
420 "lunch"}]}], ";", " ",
421 RowBox[{
422 "this", " ", "will", " ", "take", " ", "an", " ", "hour", " ", "for", " ",
423
424 RowBox[{"SMWeinberg", ".", " ", "Output"}], " ", "is", " ", "located",
425 " ", "in", " ", "FeynRules", " ",
426 RowBox[{"directory", "."}]}]}], "*)"}]], "Input",
427 CellChangeTimes->{{3.660583519626089*^9, 3.660583556727778*^9}, {
428 3.817283663764142*^9, 3.817283696526361*^9}, {3.817283782132594*^9,
429 3.817283786109083*^9}},ExpressionUUID->"ee8e0629-2939-4697-a0c4-\
430dba81be5b6fc"],
431
432Cell[BoxData[
433 RowBox[{
434 RowBox[{"WriteCT", "[",
435 RowBox[{
436 "\"\<SMWeinberg_FA/SMWeinberg_FA\>\"", ",",
437 "\"\<SMWeinberg_FA/SMWeinberg_FA\>\"", ",",
438 RowBox[{"Output", "\[Rule]", "\"\<SMWeinbergNLO\>\""}], ",",
439 "\[IndentingNewLine]",
440 RowBox[{"ZeroMom", "\[Rule]",
441 RowBox[{"{",
442 RowBox[{"{",
443 RowBox[{"aS", ",",
444 RowBox[{"{",
445 RowBox[{
446 RowBox[{"F", "[", "7", "]"}], ",",
447 RowBox[{"V", "[", "4", "]"}], ",",
448 RowBox[{"-",
449 RowBox[{"F", "[", "7", "]"}]}]}], "}"}]}], "}"}], "}"}]}], ",",
450 RowBox[{"ComplexMass", "\[Rule]", "False"}], ",", "\[IndentingNewLine]",
451 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
452 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}], ";"}]], "Input",\
453ExpressionUUID->"33fe1c0b-2e33-4218-af6d-a78907ba4306"]
454}, Open ]],
455
456Cell[CellGroupData[{
457
458Cell["Load BSM.nlo and Generate UFO at NLO", "Section",
459 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
460 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
461 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
462 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
463 3.6550378926486187`*^9}, {3.660583469669636*^9,
464 3.6605834886500196`*^9}},ExpressionUUID->"493a4b49-0255-47ae-87b0-\
465254220da35ea"],
466
467Cell[BoxData[
468 RowBox[{"(*", " ",
469 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
470 "*)"}]], "Input",ExpressionUUID->"6c2b66dd-ba7d-4d84-96e6-433ecd080754"],
471
472Cell[BoxData[
473 RowBox[{
474 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
475 CellChangeTimes->{{3.660583502811754*^9, 3.6605835046945877`*^9},
476 3.66212390294204*^9},
477 CellLabel->"In[1]:=",ExpressionUUID->"8506311a-f288-4488-903d-47a5940847e6"],
478
479Cell[BoxData[{
480 RowBox[{
481 RowBox[{"$CurrentPath", "=",
482 RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
483 RowBox[{
484 RowBox[{"$FeynRulesPath", "=",
485 RowBox[{
486 "SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}]}],
487 ";"}], "\[IndentingNewLine]",
488 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
489 RowBox[{
490 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
491 CellChangeTimes->{{3.6648666749712143`*^9, 3.664866676189664*^9}, {
492 3.800257218044265*^9, 3.80025723793865*^9}, {3.814539767536695*^9,
493 3.8145397709138327`*^9}},
494 CellLabel->"In[1]:=",ExpressionUUID->"95f917ed-a9c9-4c6c-b687-25a83f1532e1"],
495
496Cell[BoxData[{
497 RowBox[{
498 RowBox[{"LoadModel", "[",
499 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<weinbergD5_NLO.fr\>\""}], "]"}],
500 ";"}], "\[IndentingNewLine]",
501 RowBox[{
502 RowBox[{
503 RowBox[{"LoadRestriction", "[",
504 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<Massless_5f.rst\>\""}],
505 "]"}], ";"}], "\[IndentingNewLine]",
506 RowBox[{"(*",
507 RowBox[{
508 RowBox[{"LoadRestriction", "[",
509 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<Massless_4f.rst\>\""}],
510 "]"}], ";"}], "*)"}]}], "\[IndentingNewLine]",
511 RowBox[{
512 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}], "Input",
513 CellChangeTimes->{{3.662827026079749*^9, 3.662827026713024*^9},
514 3.8002572683646812`*^9, 3.800263542486582*^9, 3.8002860389151506`*^9, {
515 3.814192829834181*^9, 3.814192829893564*^9}, {3.814537900663053*^9,
516 3.8145379067232513`*^9}, {3.814538967510811*^9, 3.8145389754832497`*^9}, {
517 3.814539959130856*^9, 3.814539965058557*^9}, {3.8145434878319597`*^9,
518 3.814543487952228*^9}, {3.814544592237762*^9,
519 3.814544592291259*^9}},ExpressionUUID->"39c2e736-287a-43b4-afd0-\
520e71ce7e12993"],
521
522Cell[BoxData[{
523 RowBox[{
524 RowBox[{"Get", "[", "\"\<./SMWeinbergNLO.nlo\>\"", "]"}],
525 ";"}], "\[IndentingNewLine]",
526 RowBox[{
527 RowBox[{"LWeinberg", " ", "=", " ",
528 RowBox[{"ExpandIndices", "/@",
529 RowBox[{"Expand", "[",
530 RowBox[{"LSM", "+", "LD5"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
531 RowBox[{
532 RowBox[{"WriteUFO", "[",
533 RowBox[{"LWeinberg", ",",
534 RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
535 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
536 RowBox[{"Output", "\[Rule]", "\"\<SMWeinbergNLO\>\""}]}], "]"}],
537 ";"}], "\[IndentingNewLine]",
538 RowBox[{
539 RowBox[{"(*",
540 RowBox[{
541 RowBox[{"WriteUFO", "[",
542 RowBox[{"LWeinberg", ",",
543 RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
544 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
545 RowBox[{"Output", "\[Rule]", "\"\<SMWeinbergNLO_4FS\>\""}]}], "]"}],
546 ";"}], "*)"}]}]}], "Input",
547 CellChangeTimes->{{3.814544632491803*^9, 3.814544647510696*^9}, {
548 3.817283721438352*^9, 3.817283721686819*^9}, {3.817283771800737*^9,
549 3.8172837727700253`*^9}},ExpressionUUID->"e7724719-7b0c-4900-bb8e-\
55089b7365b3b21"]
551}, Open ]]
552}, Open ]]
553},
554WindowSize->{1330, 754},
555WindowMargins->{{44, Automatic}, {9, Automatic}},
556PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
557ShowSelection->True,
558FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)",
559StyleDefinitions->"Default.nb"
560]
561(* End of Notebook Content *)
562
563(* Internal cache information *)
564(*CellTagsOutline
565CellTagsIndex->{}
566*)
567(*CellTagsIndex
568CellTagsIndex->{}
569*)
570(*NotebookFileOutline
571Notebook[{
572Cell[545, 20, 487, 10, 30, "Input",ExpressionUUID->"cee4e897-cbc1-4bf1-a304-44bbfcd5aed3"],
573Cell[1035, 32, 1108, 26, 73, "Input",ExpressionUUID->"4149d1d2-f5f9-4692-ab77-ed3366c4e0f7"],
574Cell[2146, 60, 455, 8, 30, "Input",ExpressionUUID->"69969515-0b1d-48d5-8594-bd13de83a811"],
575Cell[CellGroupData[{
576Cell[2626, 72, 251, 4, 98, "Title",ExpressionUUID->"36a8985f-3b71-4bb1-acae-6fbe5493c081"],
577Cell[2880, 78, 1545, 26, 94, "Input",ExpressionUUID->"2b566197-5488-4f54-9e47-2d15cdf0dd38"],
578Cell[4428, 106, 529, 11, 30, "Input",ExpressionUUID->"d74de1e8-93f7-4afe-980e-8f4e6a416a1f"],
579Cell[4960, 119, 1340, 29, 94, "Input",ExpressionUUID->"a4b9bd52-5e8a-4f48-8fba-852bc4d59da3"],
580Cell[CellGroupData[{
581Cell[6325, 152, 453, 7, 67, "Section",ExpressionUUID->"430f6c12-ab98-447b-8fa4-8a13c875b6c3"],
582Cell[6781, 161, 204, 4, 30, "Input",ExpressionUUID->"01776d56-05d4-475b-b025-2f0ebd0c8fc8"],
583Cell[6988, 167, 309, 5, 30, "Input",ExpressionUUID->"4b15b6dc-b821-4c12-bd21-ceb678b9a346"],
584Cell[7300, 174, 229, 3, 30, "Input",ExpressionUUID->"cc1f0f3b-1d9f-45c8-a8b8-49b2fdacb74f"],
585Cell[7532, 179, 155, 3, 30, "Input",ExpressionUUID->"f0517d40-71c0-4541-a5c1-c1cbc531173e"],
586Cell[7690, 184, 211, 4, 30, "Input",ExpressionUUID->"e4b33362-cb94-4ba2-a771-443401bf60bf"],
587Cell[7904, 190, 210, 4, 30, "Input",ExpressionUUID->"f71c29b1-311a-4c25-824e-bfe62f098b98"]
588}, Open ]],
589Cell[CellGroupData[{
590Cell[8151, 199, 234, 4, 67, "Section",ExpressionUUID->"963bd5da-1517-491b-8768-3e81fb288ca2"],
591Cell[8388, 205, 502, 8, 73, "Input",ExpressionUUID->"ad90ef2a-90e9-4499-9f66-242849c1984c"]
592}, Open ]],
593Cell[CellGroupData[{
594Cell[8927, 218, 459, 9, 67, "Section",ExpressionUUID->"4f7e70e8-ac63-452b-bc39-2e283d37d5b9"],
595Cell[9389, 229, 791, 17, 52, "Input",ExpressionUUID->"f33d5026-da81-408c-99ad-a2da2ac609d7"],
596Cell[10183, 248, 122, 1, 30, "Input",ExpressionUUID->"0fd184c5-e009-4777-a8c1-f1520332de5e"]
597}, Open ]],
598Cell[CellGroupData[{
599Cell[10342, 254, 230, 5, 67, "Section",ExpressionUUID->"47f009e8-0603-4d37-95af-884018211353"],
600Cell[10575, 261, 765, 17, 52, "Input",ExpressionUUID->"95ec60e1-8ba3-42ae-912b-c6b75ece5652"],
601Cell[11343, 280, 233, 6, 30, "Input",ExpressionUUID->"0a6453a6-cc72-4868-bf62-a55b89adae23"],
602Cell[11579, 288, 708, 20, 73, "Input",ExpressionUUID->"dbed1df9-509a-4a3f-981c-9d0520eb5336"]
603}, Open ]],
604Cell[CellGroupData[{
605Cell[12324, 313, 356, 6, 67, "Section",ExpressionUUID->"f35c5a82-2f44-4096-9399-e82f5b97ff8d"],
606Cell[12683, 321, 645, 14, 52, "Input",ExpressionUUID->"eafcd19a-c48f-46e1-8302-30202d5af7c6"],
607Cell[13331, 337, 737, 17, 73, "Input",ExpressionUUID->"241b18a3-6df1-40c1-9540-ebe4d334b067"]
608}, Open ]],
609Cell[CellGroupData[{
610Cell[14105, 359, 346, 6, 67, "Section",ExpressionUUID->"637e3231-e751-489a-ad59-470d98c11b8b"],
611Cell[14454, 367, 255, 6, 30, "Input",ExpressionUUID->"558b9e27-1057-416f-a429-231889bb6a86"],
612Cell[14712, 375, 148, 3, 30, "Input",ExpressionUUID->"3640fc82-001c-4a97-8551-70d99a9b683d"],
613Cell[14863, 380, 1798, 31, 136, "Input",ExpressionUUID->"cb2da490-5600-44da-b7e4-f7f288533aed"],
614Cell[16664, 413, 678, 16, 30, "Input",ExpressionUUID->"ee8e0629-2939-4697-a0c4-dba81be5b6fc"],
615Cell[17345, 431, 851, 21, 73, "Input",ExpressionUUID->"33fe1c0b-2e33-4218-af6d-a78907ba4306"]
616}, Open ]],
617Cell[CellGroupData[{
618Cell[18233, 457, 475, 7, 67, "Section",ExpressionUUID->"493a4b49-0255-47ae-87b0-254220da35ea"],
619Cell[18711, 466, 184, 3, 30, "Input",ExpressionUUID->"6c2b66dd-ba7d-4d84-96e6-433ecd080754"],
620Cell[18898, 471, 241, 5, 30, "Input",ExpressionUUID->"8506311a-f288-4488-903d-47a5940847e6"],
621Cell[19142, 478, 679, 15, 94, "Input",ExpressionUUID->"95f917ed-a9c9-4c6c-b687-25a83f1532e1"],
622Cell[19824, 495, 1093, 24, 94, "Input",ExpressionUUID->"39c2e736-287a-43b4-afd0-e71ce7e12993"],
623Cell[20920, 521, 1171, 28, 94, "Input",ExpressionUUID->"e7724719-7b0c-4900-bb8e-89b7365b3b21"]
624}, Open ]]
625}, Open ]]
626}
627]
628*)
629