ZeeBabu: zeebabu_NLO_public.nb

File zeebabu_NLO_public.nb, 21.9 KB (added by Richard Ruiz, 2 years ago)

Mathematica notebook to generate UFO files from Zee-Babu FeynRules model file.

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[ 22236, 545]
14NotebookOptionsPosition[ 18643, 484]
15NotebookOutlinePosition[ 19116, 502]
16CellTagsIndexPosition[ 19073, 499]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"(*", " ",
23 RowBox[{
24 RowBox[{"R", ".", " ", "Ruiz"}], ";", " ",
25 RowBox[{"2022", " ", "June"}]}], " ", "*)"}]], "Input",
26 CellChangeTimes->{{3.665142881510211*^9, 3.665142884049021*^9}, {
27 3.665142942303646*^9, 3.6651429529879627`*^9}, {3.8578849922141743`*^9,
28 3.857884996676725*^9}, {3.865578290201123*^9,
29 3.865578290708436*^9}},ExpressionUUID->"cee4e897-cbc1-4bf1-a304-\
3044bbfcd5aed3"],
31
32Cell[BoxData[
33 RowBox[{
34 RowBox[{"(*", " ",
35 RowBox[{
36 RowBox[{"To", " ", "generate", " ", "LO", " ",
37 RowBox[{"UFO", ":", " ",
38 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
39 RowBox[{"load", " ", "model"}], ",", " ",
40 RowBox[{"write", " ", "UFO"}]}], " ", "*)"}], "\[IndentingNewLine]",
41 RowBox[{"(*", " ",
42 RowBox[{
43 RowBox[{"To", " ", "generate", " ", "NLO", " ",
44 RowBox[{"UFO", ":", " ",
45 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
46 RowBox[{"load", " ", "model"}], ",", " ",
47 RowBox[{
48 RowBox[{"extract", " ", "on"}], "-",
49 RowBox[{"shell", " ", "renormalization", " ", "terms"}]}], ",", " ",
50 RowBox[{"write", " ", "CT", " ", "terms"}], ",", " ",
51 RowBox[{"reset", " ", "paths"}], ",", " ",
52 RowBox[{"reload", " ", "model"}], ",", " ",
53 RowBox[{"write", " ", "everything", " ", "to", " ",
54 RowBox[{"file", "."}]}]}], " ", "*)"}]}]], "Input",
55 CellChangeTimes->{{3.6648662671878157`*^9, 3.664866418502534*^9}, {
56 3.664866852892112*^9,
57 3.664866879416404*^9}},ExpressionUUID->"4149d1d2-f5f9-4692-ab77-\
58ed3366c4e0f7"],
59
60Cell[BoxData[
61 RowBox[{
62 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
63 CellChangeTimes->{{3.800203881213688*^9, 3.800203885091346*^9}, {
64 3.814192661966413*^9, 3.81419266729278*^9}, {3.814196498380603*^9,
65 3.814196499409697*^9}, {3.814196533345257*^9, 3.814196534014574*^9}, {
66 3.814537179057211*^9, 3.814537179988576*^9}, {3.817278556475902*^9,
67 3.817278557221746*^9}, {3.8308506728039427`*^9, 3.83085067396842*^9}},
68 CellLabel->"In[1]:=",ExpressionUUID->"2dd3093a-c7ea-41c7-a007-d589dc397211"],
69
70Cell[CellGroupData[{
71
72Cell["FeynRules Setup", "Title",
73 CellChangeTimes->{{3.596207285657744*^9, 3.596207286276866*^9}, {
74 3.596207330352247*^9, 3.59620733250506*^9}, {3.655035967648223*^9,
75 3.65503597269068*^9}},ExpressionUUID->"36a8985f-3b71-4bb1-acae-\
766fbe5493c081"],
77
78Cell[BoxData[{
79 RowBox[{
80 RowBox[{"$CurrentPath", "=",
81 RowBox[{"Directory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
82 RowBox[{
83 RowBox[{"$FeynRulesPath", "=",
84 RowBox[{
85 "SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}]}],
86 ";"}], "\[IndentingNewLine]",
87 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
88 RowBox[{
89 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",Express\
90ionUUID->"c02dce38-3048-4633-94e9-e77e8ef0c800"],
91
92Cell[BoxData[
93 RowBox[{"(*",
94 RowBox[{"Load", " ", "default", " ", "SM", " ", "file", " ", "v1", ".4",
95 RowBox[{".7", "."}]}], "*)"}]], "Input",
96 CellChangeTimes->{{3.6550381551475554`*^9, 3.655038155711425*^9}, {
97 3.6622141225281467`*^9, 3.662214139672209*^9}, {3.730294820431093*^9,
98 3.730294820982045*^9},
99 3.8655804703528767`*^9},ExpressionUUID->"d74de1e8-93f7-4afe-980e-\
1008f4e6a416a1f"],
101
102Cell[BoxData[{
103 RowBox[{
104 RowBox[{
105 RowBox[{"isLeptonMassive", "=", "True"}], ";"}], " ",
106 RowBox[{"(*", " ",
107 RowBox[{"switch", " ", "for", " ",
108 RowBox[{"massless", " ", "/", " ", "massive"}], " ", "charged", " ",
109 "leptons"}], " ", "*)"}]}], "\[IndentingNewLine]",
110 RowBox[{
111 RowBox[{"LoadModel", "[",
112 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<zeebabu_NLO_public.fr\>\""}], "]"}],
113 ";"}], "\[IndentingNewLine]",
114 RowBox[{
115 RowBox[{"If", "[",
116 RowBox[{"isLeptonMassive", ",", "\[IndentingNewLine]",
117 RowBox[{"LoadRestriction", "[",
118 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<massless_quarks.rst\>\""}],
119 "]"}], ",", "\[IndentingNewLine]",
120 RowBox[{"LoadRestriction", "[",
121 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<massless.rst\>\""}],
122 "]"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]",
123 RowBox[{
124 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}], "Input",ExpressionUUID->\
125"5ea88f54-35d3-4b35-bd98-245cd541651b"],
126
127Cell[CellGroupData[{
128
129Cell["New Feynman Rules for ZB Interactions", "Section",
130 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
131 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
132 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
133 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
134 3.6550378926486187`*^9}, {3.662214296374611*^9, 3.662214296674329*^9}, {
135 3.865578449133216*^9,
136 3.8655784503336477`*^9}},ExpressionUUID->"430f6c12-ab98-447b-8fa4-\
1378a13c875b6c3"],
138
139Cell[BoxData[
140 RowBox[{"(*", " ",
141 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
142 RowBox[{"appear", ":", " ",
143 RowBox[{
144 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
145\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
146 "Warning"}], ":", " ",
147 RowBox[{
148 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
149 " ",
150 RowBox[{"encountered", "!"}],
151 "\"\<;\nThis is expected as Lepton Number is violated and FeynRules \
152checks for this.\>"}]}]}], "*)"}]], "Input",ExpressionUUID->"ace2f1a4-cdb5-\
1534121-9963-39b4d1feecad"],
154
155Cell[BoxData[{"LZBKin", "\[IndentingNewLine]",
156 RowBox[{"FeynmanRules", "[", "LZBKin", "]"}]}], "Input",
157 CellChangeTimes->{{3.5723482442775097`*^9, 3.572348252612296*^9}, {
158 3.596212546796714*^9, 3.596212547440735*^9}, {3.655037865036394*^9,
159 3.6550378653420067`*^9}, {3.6622141995703773`*^9, 3.6622142141549997`*^9}, {
160 3.857877763779788*^9, 3.857877766487311*^9}, {3.865578361212719*^9,
161 3.86557836149465*^9}},ExpressionUUID->"4b15b6dc-b821-4c12-bd21-\
162ceb678b9a346"],
163
164Cell[BoxData[{"Lkll", "\[IndentingNewLine]",
165 RowBox[{"ExpandIndices", "[", "Lkll", "]"}]}], "Input",
166 CellChangeTimes->{{3.6622142327067413`*^9, 3.662214233288702*^9}, {
167 3.857877874835691*^9, 3.857877880139495*^9}, {3.857878029397729*^9,
168 3.857878052035183*^9}, {3.85787818070706*^9, 3.8578781868335533`*^9}, {
169 3.857878442957178*^9, 3.857878443624498*^9}, {3.857889837254807*^9,
170 3.8578898455728083`*^9},
171 3.865578418659614*^9},ExpressionUUID->"e6072a7f-a1db-426d-b53f-\
1720b3a0d059746"],
173
174Cell[BoxData[{"Lhll", "\[IndentingNewLine]",
175 RowBox[{"ExpandIndices", "[", "Lhll", "]"}]}], "Input",
176 CellChangeTimes->{{3.6622142357666893`*^9, 3.6622142376786947`*^9},
177 3.662826298532291*^9, {3.857877985083164*^9, 3.857877986765386*^9}, {
178 3.865578409051498*^9,
179 3.865578412702682*^9}},ExpressionUUID->"e4b33362-cb94-4ba2-a771-\
180443401bf60bf"],
181
182Cell[BoxData[
183 RowBox[{"FeynmanRules", "[", "Llnv", "]"}]], "Input",ExpressionUUID->"d12be5a0-8a03-4271-bb2f-33d3a6709310"]
184}, Open ]],
185
186Cell[CellGroupData[{
187
188Cell["Basic checks for new Lagrangian terms", "Section",
189 CellChangeTimes->{{3.4234156496121893`*^9, 3.423415657027569*^9}, {
190 3.655037883884348*^9, 3.655037906883815*^9}, {3.865293916945572*^9,
191 3.8652939246393557`*^9}},ExpressionUUID->"963bd5da-1517-491b-8768-\
1923e81fb288ca2"],
193
194Cell[BoxData[{
195 RowBox[{"CheckHermiticity", "[", "LZB", "]"}], "\[IndentingNewLine]",
196 RowBox[{"CheckMassSpectrum", "[", "LZB", "]"}], "\[IndentingNewLine]",
197 RowBox[{"CheckKineticTermNormalisation", "[", "LZB", "]"}]}], "Input",Express\
198ionUUID->"bf30a8c1-c27e-4d9e-896e-422310148051"]
199}, Open ]],
200
201Cell[CellGroupData[{
202
203Cell["UFO Output", "Section",
204 CellChangeTimes->{
205 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
206 3.59620768137986*^9,
207 3.596207684212212*^9}},ExpressionUUID->"47f009e8-0603-4d37-95af-\
208884018211353"],
209
210Cell[BoxData[
211 RowBox[{"(*", " ",
212 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
213 RowBox[{"appear", ":", " ",
214 RowBox[{
215 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
216\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
217 "Warning"}], ":", " ",
218 RowBox[{
219 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
220 " ",
221 RowBox[{"encountered", "!"}],
222 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
223checks for this.\>"}]}]}], "*)"}]], "Input",
224 CellChangeTimes->{{3.656940935403554*^9, 3.656940952617087*^9}, {
225 3.6648662530247602`*^9,
226 3.664866256652151*^9}},ExpressionUUID->"95ec60e1-8ba3-42ae-912b-\
227c6b75ece5652"],
228
229Cell[BoxData[{
230 RowBox[{
231 RowBox[{"LpQCD", " ", "=", " ",
232 RowBox[{"ExpandIndices", "/@",
233 RowBox[{"Expand", "[",
234 RowBox[{"LSM", "+", "LZB"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
235 RowBox[{
236 RowBox[{"If", "[",
237 RowBox[{"isLeptonMassive", ",", "\[IndentingNewLine]",
238 RowBox[{"WriteUFO", "[",
239 RowBox[{"LpQCD", ",",
240 RowBox[{
241 "Output", "\[Rule]", "\"\<SM_ZeeBabu_MassiveLeptons_XLO\>\""}]}], "]"}],
242 ",", "\[IndentingNewLine]",
243 RowBox[{"WriteUFO", "[",
244 RowBox[{"LpQCD", ",",
245 RowBox[{"Output", "\[Rule]", "\"\<SM_ZeeBabu_XLO\>\""}]}], "]"}]}],
246 "\[IndentingNewLine]", "]"}], ";"}]}], "Input",ExpressionUUID->"34c50153-\
2471c7b-45c0-b793-17a5f54e6210"]
248}, Open ]],
249
250Cell[CellGroupData[{
251
252Cell["On-Shell Renormalization with FeynArts", "Section",
253 CellChangeTimes->{
254 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
255 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199318805483*^9,
256 3.6621993287719097`*^9}, {3.662216149623213*^9,
257 3.662216170908692*^9}},ExpressionUUID->"f35c5a82-2f44-4096-9399-\
258e82f5b97ff8d"],
259
260Cell[BoxData[
261 RowBox[{"(*", " ",
262 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
263 RowBox[{"appear", ":", " ",
264 RowBox[{
265 "\"\<\!\(\*StyleBox[\"QN\", \"MessageName\"]\)\!\(\*StyleBox[\"::\", \
266\"MessageName\"]\)\!\(\*StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"",
267 "Warning"}], ":", " ",
268 RowBox[{
269 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
270 " ",
271 RowBox[{"encountered", "!"}],
272 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
273checks for this.\>"}]}]}], "*)"}]], "Input",ExpressionUUID->"eafcd19a-c48f-\
27446e1-8302-30202d5af7c6"],
275
276Cell[BoxData[{
277 RowBox[{
278 RowBox[{"LRen", " ", "=", " ",
279 RowBox[{"OnShellRenormalization", "[",
280 RowBox[{
281 RowBox[{"LSM", "+", "LZB"}], ",",
282 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
283 RowBox[{"FlavorMixing", "\[Rule]", " ", "False"}], ",",
284 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}]}],
285 ";"}], "\[IndentingNewLine]",
286 RowBox[{"SetDirectory", "[",
287 "\"\<~/Scripts/FeynArts/FeynArts-3.11/Models\>\"",
288 "]"}], "\[IndentingNewLine]",
289 RowBox[{
290 RowBox[{"WriteFeynArtsOutput", "[",
291 RowBox[{"LRen", ",",
292 RowBox[{"GenericFile", "\[Rule]", "True"}], ",",
293 RowBox[{"FlavorExpand", "\[Rule]", "True"}], ",",
294 RowBox[{"Output", "\[Rule]", "\"\<SM_ZeeBabu_FA\>\""}]}], "]"}],
295 ";"}]}], "Input",ExpressionUUID->"02f82583-a5c0-4ac9-8b08-6f2cc230ee4b"]
296}, Open ]],
297
298Cell[CellGroupData[{
299
300Cell["Write Counter Terms to BSM.nlo", "Section",
301 CellChangeTimes->{
302 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
303 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199461128889*^9,
304 3.662199467953211*^9}, {3.662216541485611*^9,
305 3.662216544705964*^9}},ExpressionUUID->"637e3231-e751-489a-ad59-\
306470d98c11b8b"],
307
308Cell[BoxData[
309 RowBox[{"(*", " ",
310 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
311 "*)"}]], "Input",
312 CellChangeTimes->{{3.664866502336069*^9,
313 3.664866536918194*^9}},ExpressionUUID->"558b9e27-1057-416f-a429-\
314231889bb6a86"],
315
316Cell[BoxData[
317 RowBox[{
318 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
319 CellLabel->"In[1]:=",ExpressionUUID->"3640fc82-001c-4a97-8551-70d99a9b683d"],
320
321Cell[BoxData[{
322 RowBox[{
323 RowBox[{
324 "SetDirectory", "[", "\"\<~/Scripts/FeynArts/FeynArts-3.11/\>\"", "]"}],
325 ";"}], "\[IndentingNewLine]",
326 RowBox[{"<<", "FeynArts`"}], "\[IndentingNewLine]",
327 RowBox[{
328 RowBox[{"SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}],
329 ";"}], "\[IndentingNewLine]",
330 RowBox[{"<<", "NLOCT`"}], "\[IndentingNewLine]",
331 RowBox[{
332 RowBox[{"$CurrentPath", "=",
333 RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
334 RowBox[{
335 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",Express\
336ionUUID->"415ceffa-5035-4fe0-bcb8-4cc2ac3a3a63"],
337
338Cell[BoxData[
339 RowBox[{
340 RowBox[{"(*",
341 RowBox[{
342 RowBox[{"Note", ":", " ",
343 RowBox[{"Grab", " ", "some", " ", "coffee"}]}], ";", " ",
344 RowBox[{"this", " ", "will", " ", "take", " ", "a", " ", "few", " ",
345 RowBox[{"minutes", ".", " ", "Output"}], " ", "is", " ", "located", " ",
346 "in", " ", "FeynRules", " ",
347 RowBox[{"directory", "."}]}]}], "*)"}], "\[IndentingNewLine]",
348 RowBox[{"WriteCT", "[",
349 RowBox[{
350 "\"\<SM_ZeeBabu_FA/SM_ZeeBabu_FA\>\"", ",",
351 "\"\<SM_ZeeBabu_FA/SM_ZeeBabu_FA\>\"", ",",
352 RowBox[{"Output", "\[Rule]", "\"\<SM_ZeeBabu\>\""}], ",",
353 "\[IndentingNewLine]",
354 RowBox[{"ZeroMom", "\[Rule]",
355 RowBox[{"{",
356 RowBox[{"{",
357 RowBox[{"aS", ",",
358 RowBox[{"{",
359 RowBox[{
360 RowBox[{"F", "[", "7", "]"}], ",",
361 RowBox[{"V", "[", "4", "]"}], ",",
362 RowBox[{"-",
363 RowBox[{"F", "[", "7", "]"}]}]}], "}"}]}], "}"}], "}"}]}], ",",
364 RowBox[{"ComplexMass", "\[Rule]", "False"}], ",",
365 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
366 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}]}]], "Input",
367 CellChangeTimes->{{3.865584086094201*^9,
368 3.865584090328697*^9}},ExpressionUUID->"7744ad0f-d54b-4fc9-a3e5-\
3696ec91e11b757"]
370}, Open ]],
371
372Cell[CellGroupData[{
373
374Cell["Load BSM.nlo and Generate UFO at NLO", "Section",
375 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
376 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
377 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
378 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
379 3.6550378926486187`*^9}, {3.660583469669636*^9,
380 3.6605834886500196`*^9}},ExpressionUUID->"493a4b49-0255-47ae-87b0-\
381254220da35ea"],
382
383Cell[BoxData[
384 RowBox[{"(*", " ",
385 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
386 "*)"}]], "Input",ExpressionUUID->"6c2b66dd-ba7d-4d84-96e6-433ecd080754"],
387
388Cell[BoxData[
389 RowBox[{
390 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
391 CellChangeTimes->{{3.660583502811754*^9, 3.6605835046945877`*^9},
392 3.66212390294204*^9},
393 CellLabel->"In[1]:=",ExpressionUUID->"8506311a-f288-4488-903d-47a5940847e6"],
394
395Cell[BoxData[
396 RowBox[{
397 RowBox[{"(*", " ",
398 RowBox[{"load", " ", "FeynRules"}], " ", "*)"}], "\[IndentingNewLine]",
399 RowBox[{
400 RowBox[{
401 RowBox[{"$CurrentPath", "=",
402 RowBox[{"NotebookDirectory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
403
404 RowBox[{
405 RowBox[{"$FeynRulesPath", "=",
406 RowBox[{
407 "SetDirectory", "[", "\"\<~/Scripts/FeynRules/FeynRules\>\"", "]"}]}],
408 ";"}], "\[IndentingNewLine]",
409 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
410 RowBox[{
411 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}]}]], "Input",
412 CellChangeTimes->{{3.865584286630831*^9,
413 3.8655842870432568`*^9}},ExpressionUUID->"84058f34-7b2a-47b2-a076-\
414a7a271d6cba9"],
415
416Cell[BoxData[
417 RowBox[{
418 RowBox[{"(*", " ",
419 RowBox[{"load", " ", "model", " ", "files"}], " ", "*)"}],
420 "\[IndentingNewLine]",
421 RowBox[{
422 RowBox[{
423 RowBox[{"isLeptonMassive", "=", "True"}], ";", " ",
424 RowBox[{"(*", " ",
425 RowBox[{
426 RowBox[{"switch", " ", "for", " ",
427 RowBox[{"massless", " ", "/", " ", "massive"}], " ", "charged", " ",
428 "leptons"}], ";", " ",
429 RowBox[{"must", " ", "match", " ", "above"}]}], " ", "*)"}],
430 "\[IndentingNewLine]",
431 RowBox[{"LoadModel", "[",
432 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<zeebabu_NLO_public.fr\>\""}], "]"}],
433 ";"}], "\[IndentingNewLine]",
434 RowBox[{
435 RowBox[{"If", "[",
436 RowBox[{"isLeptonMassive", ",", "\[IndentingNewLine]",
437 RowBox[{"LoadRestriction", "[",
438 RowBox[{
439 "\"\<diagonalCKM.rst\>\"", ",", "\"\<massless_quarks.rst\>\""}], "]"}],
440 ",", "\[IndentingNewLine]",
441 RowBox[{"LoadRestriction", "[",
442 RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<massless.rst\>\""}],
443 "]"}]}], "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]",
444 RowBox[{
445 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}]}]], "Input",
446 CellChangeTimes->{{3.865584246879266*^9,
447 3.865584252863028*^9}},ExpressionUUID->"06c39969-fb4d-4bdb-9dae-\
448a2d392b421e8"],
449
450Cell[BoxData[
451 RowBox[{
452 RowBox[{"(*", " ",
453 RowBox[{"generate", " ", "UFO"}], " ", "*)"}], "\[IndentingNewLine]",
454 RowBox[{
455 RowBox[{
456 RowBox[{"Get", "[", "\"\<./SM_ZeeBabu.nlo\>\"", "]"}], ";"}],
457 "\[IndentingNewLine]",
458 RowBox[{
459 RowBox[{"LpQCD", " ", "=", " ",
460 RowBox[{"ExpandIndices", "/@",
461 RowBox[{"Expand", "[",
462 RowBox[{"LSM", "+", "LZB"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
463 RowBox[{
464 RowBox[{"If", "[",
465 RowBox[{"isLeptonMassive", ",", "\[IndentingNewLine]",
466 RowBox[{"WriteUFO", "[",
467 RowBox[{"LpQCD", ",",
468 RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
469 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
470 RowBox[{
471 "Output", "\[Rule]", "\"\<SM_ZeeBabu_MassiveLeptons_NLO\>\""}]}],
472 "]"}], ",", "\[IndentingNewLine]",
473 RowBox[{"WriteUFO", "[",
474 RowBox[{"LpQCD", ",",
475 RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
476 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
477 RowBox[{"Output", "\[Rule]", "\"\<SM_ZeeBabu_NLO\>\""}]}], "]"}]}],
478 "\[IndentingNewLine]", "]"}], ";"}]}]}]], "Input",
479 CellChangeTimes->{{3.865584131474765*^9,
480 3.865584136767414*^9}},ExpressionUUID->"335b5ec7-5fa2-4add-be4c-\
4816a63111321ae"]
482}, Open ]]
483}, Open ]]
484},
485WindowSize->{1330, 918},
486WindowMargins->{{Automatic, 18}, {0, Automatic}},
487PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
488ShowSelection->True,
489TaggingRules->{"TryRealOnly" -> False},
490FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
491StyleDefinitions->"Default.nb"
492]
493(* End of Notebook Content *)
494
495(* Internal cache information *)
496(*CellTagsOutline
497CellTagsIndex->{}
498*)
499(*CellTagsIndex
500CellTagsIndex->{}
501*)
502(*NotebookFileOutline
503Notebook[{
504Cell[545, 20, 426, 9, 31, "Input",ExpressionUUID->"cee4e897-cbc1-4bf1-a304-44bbfcd5aed3"],
505Cell[974, 31, 1108, 26, 78, "Input",ExpressionUUID->"4149d1d2-f5f9-4692-ab77-ed3366c4e0f7"],
506Cell[2085, 59, 502, 8, 31, "Input",ExpressionUUID->"2dd3093a-c7ea-41c7-a007-d589dc397211"],
507Cell[CellGroupData[{
508Cell[2612, 71, 251, 4, 99, "Title",ExpressionUUID->"36a8985f-3b71-4bb1-acae-6fbe5493c081"],
509Cell[2866, 77, 483, 12, 101, "Input",ExpressionUUID->"c02dce38-3048-4633-94e9-e77e8ef0c800"],
510Cell[3352, 91, 404, 8, 31, "Input",ExpressionUUID->"d74de1e8-93f7-4afe-980e-8f4e6a416a1f"],
511Cell[3759, 101, 990, 23, 170, "Input",ExpressionUUID->"5ea88f54-35d3-4b35-bd98-245cd541651b"],
512Cell[CellGroupData[{
513Cell[4774, 128, 525, 8, 68, "Section",ExpressionUUID->"430f6c12-ab98-447b-8fa4-8a13c875b6c3"],
514Cell[5302, 138, 645, 14, 55, "Input",ExpressionUUID->"ace2f1a4-cdb5-4121-9963-39b4d1feecad"],
515Cell[5950, 154, 477, 7, 55, "Input",ExpressionUUID->"4b15b6dc-b821-4c12-bd21-ceb678b9a346"],
516Cell[6430, 163, 502, 8, 55, "Input",ExpressionUUID->"e6072a7f-a1db-426d-b53f-0b3a0d059746"],
517Cell[6935, 173, 355, 6, 55, "Input",ExpressionUUID->"e4b33362-cb94-4ba2-a771-443401bf60bf"],
518Cell[7293, 181, 123, 1, 31, "Input",ExpressionUUID->"d12be5a0-8a03-4271-bb2f-33d3a6709310"]
519}, Open ]],
520Cell[CellGroupData[{
521Cell[7453, 187, 281, 4, 68, "Section",ExpressionUUID->"963bd5da-1517-491b-8768-3e81fb288ca2"],
522Cell[7737, 193, 288, 4, 78, "Input",ExpressionUUID->"bf30a8c1-c27e-4d9e-896e-422310148051"]
523}, Open ]],
524Cell[CellGroupData[{
525Cell[8062, 202, 230, 5, 68, "Section",ExpressionUUID->"47f009e8-0603-4d37-95af-884018211353"],
526Cell[8295, 209, 765, 17, 55, "Input",ExpressionUUID->"95ec60e1-8ba3-42ae-912b-c6b75ece5652"],
527Cell[9063, 228, 716, 18, 124, "Input",ExpressionUUID->"34c50153-1c7b-45c0-b793-17a5f54e6210"]
528}, Open ]],
529Cell[CellGroupData[{
530Cell[9816, 251, 356, 6, 68, "Section",ExpressionUUID->"f35c5a82-2f44-4096-9399-e82f5b97ff8d"],
531Cell[10175, 259, 645, 14, 55, "Input",ExpressionUUID->"eafcd19a-c48f-46e1-8302-30202d5af7c6"],
532Cell[10823, 275, 813, 19, 78, "Input",ExpressionUUID->"02f82583-a5c0-4ac9-8b08-6f2cc230ee4b"]
533}, Open ]],
534Cell[CellGroupData[{
535Cell[11673, 299, 346, 6, 68, "Section",ExpressionUUID->"637e3231-e751-489a-ad59-470d98c11b8b"],
536Cell[12022, 307, 255, 6, 31, "Input",ExpressionUUID->"558b9e27-1057-416f-a429-231889bb6a86"],
537Cell[12280, 315, 148, 3, 31, "Input",ExpressionUUID->"3640fc82-001c-4a97-8551-70d99a9b683d"],
538Cell[12431, 320, 630, 15, 147, "Input",ExpressionUUID->"415ceffa-5035-4fe0-bcb8-4cc2ac3a3a63"],
539Cell[13064, 337, 1270, 31, 78, "Input",ExpressionUUID->"7744ad0f-d54b-4fc9-a3e5-6ec91e11b757"]
540}, Open ]],
541Cell[CellGroupData[{
542Cell[14371, 373, 475, 7, 68, "Section",ExpressionUUID->"493a4b49-0255-47ae-87b0-254220da35ea"],
543Cell[14849, 382, 184, 3, 31, "Input",ExpressionUUID->"6c2b66dd-ba7d-4d84-96e6-433ecd080754"],
544Cell[15036, 387, 241, 5, 31, "Input",ExpressionUUID->"8506311a-f288-4488-903d-47a5940847e6"],
545Cell[15280, 394, 709, 19, 124, "Input",ExpressionUUID->"84058f34-7b2a-47b2-a076-a7a271d6cba9"],
546Cell[15992, 415, 1301, 32, 193, "Input",ExpressionUUID->"06c39969-fb4d-4bdb-9dae-a2d392b421e8"],
547Cell[17296, 449, 1319, 31, 170, "Input",ExpressionUUID->"335b5ec7-5fa2-4add-be4c-6a63111321ae"]
548}, Open ]]
549}, Open ]]
550}
551]
552*)
553