EffLRSM: effLRSMnlo.nb

File effLRSMnlo.nb, 27.3 KB (added by Richard Ruiz, 8 years ago)

Mathematica notebook file that generates UFO file from FeynRules model files. Allows user to also run quick sanity checks (optional) on model.

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[ 27774, 733]
14NotebookOptionsPosition[ 24630, 628]
15NotebookOutlinePosition[ 25067, 645]
16CellTagsIndexPosition[ 25024, 642]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"(*", " ",
23 RowBox[{
24 RowBox[{"R", ".", " ", "Ruiz"}], ";", " ",
25 RowBox[{"2016", " ", "October", " ", "30"}]}], " ", "*)"}]], "Input",
26 CellChangeTimes->{{3.665142881510211*^9, 3.665142884049021*^9}, {
27 3.665142942303646*^9, 3.6651429529879627`*^9}, {3.686847218103582*^9,
28 3.686847221267969*^9}}],
29
30Cell[BoxData[
31 RowBox[{
32 RowBox[{"(*", " ",
33 RowBox[{
34 RowBox[{"To", " ", "generate", " ", "LO", " ",
35 RowBox[{"UFO", ":", " ",
36 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
37 RowBox[{"load", " ", "model"}], ",", " ",
38 RowBox[{"write", " ", "UFO"}]}], " ", "*)"}], "\[IndentingNewLine]",
39 RowBox[{"(*", " ",
40 RowBox[{
41 RowBox[{"To", " ", "generate", " ", "NLO", " ",
42 RowBox[{"UFO", ":", " ",
43 RowBox[{"Set", " ", "paths"}]}]}], ",", " ",
44 RowBox[{"load", " ", "model"}], ",", " ",
45 RowBox[{
46 RowBox[{"extract", " ", "on"}], "-",
47 RowBox[{"shell", " ", "renormalization", " ", "terms"}]}], ",", " ",
48 RowBox[{"write", " ", "CT", " ", "terms"}], ",", " ",
49 RowBox[{"reset", " ", "paths"}], ",", " ",
50 RowBox[{"reload", " ", "model"}], ",", " ",
51 RowBox[{"write", " ", "everything", " ", "to", " ",
52 RowBox[{"file", "."}]}]}], " ", "*)"}]}]], "Input",
53 CellChangeTimes->{{3.6648662671878157`*^9, 3.664866418502534*^9}, {
54 3.664866852892112*^9, 3.664866879416404*^9}}],
55
56Cell[BoxData[
57 RowBox[{
58 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
59 CellChangeTimes->{{3.680432474319689*^9, 3.680432479018909*^9}}],
60
61Cell[CellGroupData[{
62
63Cell["FeynRules Setup", "Title",
64 CellChangeTimes->{{3.596207285657744*^9, 3.596207286276866*^9}, {
65 3.596207330352247*^9, 3.59620733250506*^9}, {3.655035967648223*^9,
66 3.65503597269068*^9}}],
67
68Cell[BoxData[{
69 RowBox[{
70 RowBox[{"$CurrentPath", "=",
71 RowBox[{"Directory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
72 RowBox[{
73 RowBox[{"$FeynRulesPath", "=",
74 RowBox[{
75 "SetDirectory", "[", "\"\<~/.../FeynRules/feynrules-2_x\>\"", "]"}]}],
76 ";"}], "\[IndentingNewLine]",
77 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
78 RowBox[{
79 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
80 CellChangeTimes->{{3.41265862251538*^9, 3.412658649947229*^9}, {
81 3.423415585782702*^9, 3.423415597189939*^9}, {3.4234163173467493`*^9,
82 3.4234163227881193`*^9}, {3.572102427116931*^9, 3.572102456206238*^9}, {
83 3.572280529035656*^9, 3.572280529468644*^9}, {3.572366766720894*^9,
84 3.572366773370627*^9}, {3.596207212826934*^9, 3.5962072447909603`*^9}, {
85 3.596914045279478*^9, 3.5969140472383204`*^9}, 3.596914102616413*^9,
86 3.596915921742405*^9, {3.596915961175353*^9, 3.59691596259039*^9}, {
87 3.60744636894991*^9, 3.607446377507101*^9}, {3.6364757082312813`*^9,
88 3.636475708315627*^9}, {3.655035981353272*^9, 3.655035982244761*^9}, {
89 3.6550360677717857`*^9, 3.655036081259776*^9}, {3.655036117419112*^9,
90 3.655036143336364*^9}, {3.655037954145546*^9, 3.6550379644813023`*^9}, {
91 3.6605209855658293`*^9, 3.6605209855964212`*^9}, {3.6605772969976788`*^9,
92 3.660577300733718*^9}, {3.6621928886588717`*^9, 3.662192889375732*^9}, {
93 3.6628197402754908`*^9, 3.662819763199905*^9}, {3.6868516060928926`*^9,
94 3.686851606385346*^9}}],
95
96Cell[BoxData[
97 RowBox[{"(*",
98 RowBox[{
99 "Load", " ", "default", " ", "SM", " ", "file", " ", "v1", ".4", ".6", " ",
100 "with", " ", "Eff", " ", "LRSM", " ",
101 RowBox[{"extension", "."}]}], "*)"}]], "Input",
102 CellChangeTimes->{{3.6550381551475554`*^9, 3.655038155711425*^9}, {
103 3.6622141225281467`*^9, 3.662214139672209*^9}, {3.682240952486841*^9,
104 3.6822409556512823`*^9}}],
105
106Cell[BoxData[{
107 RowBox[{
108 RowBox[{"LoadModel", "[",
109 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<effLRSM.fr\>\""}], "]"}],
110 ";"}], "\[IndentingNewLine]",
111 RowBox[{
112 RowBox[{"LoadRestriction", "[",
113 RowBox[{"\"\<massless.rst\>\"", ",", "\"\<diagonalMixing.rst\>\""}], "]"}],
114 ";"}], "\[IndentingNewLine]",
115 RowBox[{
116 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}], "Input",
117 CellChangeTimes->{{3.680432699128376*^9, 3.6804326995380583`*^9}, {
118 3.680432926277862*^9, 3.6804329279472427`*^9}, {3.68045424820811*^9,
119 3.680454257097041*^9}, {3.686847264199445*^9, 3.686847272821652*^9}}],
120
121Cell[CellGroupData[{
122
123Cell["Lagrangian Terms", "Section",
124 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
125 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
126 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
127 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
128 3.6550378926486187`*^9}, {3.662214296374611*^9, 3.662214296674329*^9}}],
129
130Cell[BoxData["LWRTmp"], "Input",
131 CellChangeTimes->{{3.664866736738257*^9, 3.664866737292602*^9}, {
132 3.6804327091032457`*^9, 3.680432711330409*^9}}],
133
134Cell[BoxData["LZRTmp"], "Input",
135 CellChangeTimes->{{3.680432715205942*^9, 3.6804327154717617`*^9}}],
136
137Cell[BoxData["LNKin"], "Input",
138 CellChangeTimes->{{3.5723482442775097`*^9, 3.572348252612296*^9}, {
139 3.596212546796714*^9, 3.596212547440735*^9}, {3.655037865036394*^9,
140 3.6550378653420067`*^9}, {3.6622141995703773`*^9, 3.6622142141549997`*^9}}],
141
142Cell[BoxData[
143 RowBox[{"FeynmanRules", "[", "LagLRSM", "]"}]], "Input",
144 CellChangeTimes->{
145 3.664866752982409*^9, {3.680432965086788*^9, 3.680432967514791*^9}}],
146
147Cell[BoxData[
148 RowBox[{
149 RowBox[{"(*", " ",
150 RowBox[{"LFull", " ", "=", " ",
151 RowBox[{"LSM", " ", "+", " ", "LagLRSM"}]}], " ", "*)"}],
152 "\[IndentingNewLine]",
153 RowBox[{"FeynmanRules", "[", "LagFull", "]"}]}]], "Input",
154 CellChangeTimes->{{3.662214278069661*^9, 3.66221429036182*^9}, {
155 3.6804329724299297`*^9, 3.680432984734503*^9}}]
156}, Open ]],
157
158Cell[CellGroupData[{
159
160Cell["Sanity Checks for New Lagrangian Terms", "Section",
161 CellChangeTimes->{{3.4234156496121893`*^9, 3.423415657027569*^9}, {
162 3.655037883884348*^9, 3.655037906883815*^9}}],
163
164Cell[CellGroupData[{
165
166Cell[BoxData[
167 RowBox[{"CheckHermiticity", "[", "LagLRSM", "]"}]], "Input",
168 CellChangeTimes->{
169 3.607446407230465*^9, {3.68043298839073*^9, 3.680432990821744*^9}}],
170
171Cell[CellGroupData[{
172
173Cell[BoxData["\<\"Checking for hermiticity by calculating the Feynman rules \
174contained in L-HC[L].\"\>"], "Print",
175 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
176 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
177 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
178 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
179 3.682346058917761*^9, 3.6868474512273293`*^9, 3.686849867990106*^9}],
180
181Cell[BoxData["\<\"If the lagrangian is hermitian, then the number of vertices \
182should be zero.\"\>"], "Print",
183 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
184 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
185 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
186 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
187 3.682346058917761*^9, 3.6868474512273293`*^9, 3.6868498679921913`*^9}],
188
189Cell[BoxData[
190 StyleBox["\<\"Starting Feynman rule calculation.\"\>",
191 StripOnInput->False,
192 LineColor->RGBColor[1, 0.5, 0],
193 FrontFaceColor->RGBColor[1, 0.5, 0],
194 BackFaceColor->RGBColor[1, 0.5, 0],
195 GraphicsColor->RGBColor[1, 0.5, 0],
196 FontWeight->Bold,
197 FontColor->RGBColor[1, 0.5, 0]]], "Print",
198 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
199 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
200 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
201 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
202 3.682346058917761*^9, 3.6868474512273293`*^9, 3.686849871094857*^9}],
203
204Cell[BoxData["\<\"Expanding the Lagrangian...\"\>"], "Print",
205 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
206 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
207 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
208 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
209 3.682346058917761*^9, 3.6868474512273293`*^9, 3.6868498710972843`*^9}],
210
211Cell[BoxData["\<\"No vertices found.\"\>"], "Print",
212 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
213 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
214 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
215 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
216 3.682346058917761*^9, 3.6868474512273293`*^9, 3.686849871098793*^9}],
217
218Cell[BoxData[
219 InterpretationBox[
220 RowBox[{"0", "\[InvisibleSpace]", "\<\" vertices obtained.\"\>"}],
221 SequenceForm[0, " vertices obtained."],
222 Editable->False]], "Print",
223 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
224 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
225 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
226 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
227 3.682346058917761*^9, 3.6868474512273293`*^9, 3.686849871100977*^9}],
228
229Cell[BoxData["\<\"The lagrangian is hermitian.\"\>"], "Print",
230 CellChangeTimes->{3.6760964157649927`*^9, 3.6804329942285748`*^9,
231 3.68224069823027*^9, 3.682240988832697*^9, 3.682241062538535*^9,
232 3.682242153019915*^9, 3.682249412718433*^9, 3.682330147646923*^9,
233 3.68233194987496*^9, 3.682332402523697*^9, 3.682333082927475*^9,
234 3.682346058917761*^9, 3.6868474512273293`*^9, 3.686849871103045*^9}]
235}, Open ]],
236
237Cell[BoxData[
238 RowBox[{"{", "}"}]], "Output",
239 CellChangeTimes->{3.676096421299799*^9, 3.6804329969955187`*^9,
240 3.682240703168305*^9, 3.6822409930052433`*^9, 3.6822410649946947`*^9,
241 3.682242155903184*^9, 3.682249415419217*^9, 3.682330150500434*^9,
242 3.6823319527628393`*^9, 3.682332405473526*^9, 3.682333086086486*^9,
243 3.68234606197195*^9, 3.686847454237342*^9, 3.686849871104382*^9}]
244}, Open ]],
245
246Cell[CellGroupData[{
247
248Cell[BoxData[
249 RowBox[{"CheckDiagonalMassTerms", "[", "LagLRSM", "]"}]], "Input",
250 CellChangeTimes->{{3.682346065096884*^9, 3.6823460776896048`*^9}}],
251
252Cell[CellGroupData[{
253
254Cell[BoxData[
255 InterpretationBox[
256 RowBox[{"\<\"Neglecting all terms with more than \"\>",
257 "\[InvisibleSpace]", "\<\"2\"\>",
258 "\[InvisibleSpace]", "\<\" particles.\"\>"}],
259 SequenceForm["Neglecting all terms with more than ", "2", " particles."],
260 Editable->False]], "Print",
261 CellChangeTimes->{3.682346079161729*^9, 3.686847456016287*^9,
262 3.6868498722964783`*^9}],
263
264Cell[BoxData["\<\"All mass terms are diagonal.\"\>"], "Print",
265 CellChangeTimes->{3.682346079161729*^9, 3.686847456016287*^9,
266 3.686849872336046*^9}]
267}, Open ]],
268
269Cell[BoxData["True"], "Output",
270 CellChangeTimes->{3.682346079207398*^9, 3.686847456038209*^9,
271 3.686849872337475*^9}]
272}, Open ]],
273
274Cell[CellGroupData[{
275
276Cell[BoxData[
277 RowBox[{"CheckMassSpectrum", "[", "LagLRSM", "]"}]], "Input",
278 CellChangeTimes->{{3.607446409988699*^9, 3.6074464110954323`*^9},
279 3.6550378392495413`*^9, 3.680432999605236*^9}],
280
281Cell[CellGroupData[{
282
283Cell[BoxData[
284 InterpretationBox[
285 RowBox[{"\<\"Neglecting all terms with more than \"\>",
286 "\[InvisibleSpace]", "\<\"2\"\>",
287 "\[InvisibleSpace]", "\<\" particles.\"\>"}],
288 SequenceForm["Neglecting all terms with more than ", "2", " particles."],
289 Editable->False]], "Print",
290 CellChangeTimes->{3.6760964240892897`*^9, 3.6804330019810057`*^9,
291 3.682241068968843*^9, 3.682346100456819*^9, 3.686847459186689*^9,
292 3.6868498761499863`*^9}],
293
294Cell[BoxData["\<\"All mass terms are diagonal.\"\>"], "Print",
295 CellChangeTimes->{3.6760964240892897`*^9, 3.6804330019810057`*^9,
296 3.682241068968843*^9, 3.682346100456819*^9, 3.686847459186689*^9,
297 3.686849876186857*^9}],
298
299Cell[BoxData["\<\"Getting mass spectrum.\"\>"], "Print",
300 CellChangeTimes->{3.6760964240892897`*^9, 3.6804330019810057`*^9,
301 3.682241068968843*^9, 3.682346100456819*^9, 3.686847459186689*^9,
302 3.686849877052225*^9}],
303
304Cell[BoxData["\<\"Checking for less then 0.1% agreement with model file \
305values.\"\>"], "Print",
306 CellChangeTimes->{3.6760964240892897`*^9, 3.6804330019810057`*^9,
307 3.682241068968843*^9, 3.682346100456819*^9, 3.686847459186689*^9,
308 3.686849877054696*^9}]
309}, Open ]],
310
311Cell[BoxData[
312 TagBox[GridBox[{
313 {"\<\"Particle\"\>", "\<\"Analytic value\"\>", "\<\"Numerical value\"\>", \
314"\<\"Model-file value\"\>"},
315 {"N1", "mN1", "173.3`", "173.3`"},
316 {"N2", "mN2", "1.`*^12", "1.`*^12"},
317 {"N3", "mN3", "1.`*^14", "1.`*^14"}
318 },
319 GridBoxAlignment->{
320 "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
321 "RowsIndexed" -> {}},
322 GridBoxSpacings->{"Columns" -> {
323 Offset[0.27999999999999997`], {
324 Offset[2.0999999999999996`]},
325 Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
326 Offset[0.2], {
327 Offset[0.4]},
328 Offset[0.2]}, "RowsIndexed" -> {}}],
329 Function[BoxForm`e$,
330 TableForm[BoxForm`e$]]]], "Output",
331 CellChangeTimes->{3.6760964255329437`*^9, 3.680433002777938*^9,
332 3.682241069641775*^9, 3.682346101314436*^9, 3.6868474605209923`*^9,
333 3.686849877057846*^9}]
334}, Open ]],
335
336Cell[CellGroupData[{
337
338Cell[BoxData[
339 RowBox[{"CheckKineticTermNormalisation", "[", "LagLRSM", "]"}]], "Input",
340 CellChangeTimes->{{3.414763473002973*^9, 3.4147634845507174`*^9}, {
341 3.572348320182221*^9, 3.572348326613598*^9}, {3.572358912362153*^9,
342 3.572358917148378*^9}, {3.596207478919518*^9, 3.596207486580676*^9}, {
343 3.596212578276376*^9, 3.5962125890639143`*^9}, 3.607446414234477*^9,
344 3.680433031289057*^9}],
345
346Cell[CellGroupData[{
347
348Cell[BoxData[
349 InterpretationBox[
350 RowBox[{"\<\"Neglecting all terms with more than \"\>",
351 "\[InvisibleSpace]", "\<\"2\"\>",
352 "\[InvisibleSpace]", "\<\" particles.\"\>"}],
353 SequenceForm["Neglecting all terms with more than ", "2", " particles."],
354 Editable->False]], "Print",
355 CellChangeTimes->{3.676096428280446*^9, 3.68043303458121*^9,
356 3.682241076298921*^9, 3.682346107202651*^9, 3.686847464684845*^9,
357 3.686849881961215*^9}],
358
359Cell[BoxData["\<\"All kinetic terms are diagonal.\"\>"], "Print",
360 CellChangeTimes->{3.676096428280446*^9, 3.68043303458121*^9,
361 3.682241076298921*^9, 3.682346107202651*^9, 3.686847464684845*^9,
362 3.6868498819967613`*^9}],
363
364Cell[BoxData["\<\"All kinetic terms are correctly normalized.\"\>"], "Print",
365 CellChangeTimes->{3.676096428280446*^9, 3.68043303458121*^9,
366 3.682241076298921*^9, 3.682346107202651*^9, 3.686847464684845*^9,
367 3.686849882859892*^9}]
368}, Open ]],
369
370Cell[BoxData["True"], "Output",
371 CellChangeTimes->{3.676096429594046*^9, 3.6804330353281107`*^9,
372 3.682241076962926*^9, 3.682346108233881*^9, 3.686847465513378*^9,
373 3.686849882862053*^9}]
374}, Open ]]
375}, Open ]],
376
377Cell[CellGroupData[{
378
379Cell["\<\
380New Feynman Rules for LRSM N Interactions\
381\>", "Section",
382 CellChangeTimes->{{3.411910637332127*^9, 3.41191064360071*^9}, {
383 3.655037529059973*^9, 3.655037541428491*^9}, {3.6550379116719217`*^9,
384 3.655037913954034*^9}, {3.6804330876628723`*^9, 3.6804330880539503`*^9}}],
385
386Cell[BoxData[
387 RowBox[{"(*", " ",
388 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
389 RowBox[{"appear", ":", " ",
390 RowBox[{"\"\<\!\(\*
391StyleBox[\"QN\", \"MessageName\"]\)\!\(\*
392StyleBox[\"::\", \"MessageName\"]\)\!\(\*
393StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"", "Warning"}], ":", " ",
394 RowBox[{
395 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
396 " ",
397 RowBox[{"encountered", "!"}],
398 "\"\<;\nThis is expected as Lepton Number is violated and FeynRules \
399checks for this.\>"}]}]}], "*)"}]], "Input",
400 CellChangeTimes->{{3.6550374096249533`*^9, 3.655037447663039*^9},
401 3.655037551040966*^9, {3.656940976048616*^9, 3.6569409786068287`*^9}}],
402
403Cell[BoxData[
404 RowBox[{"FeynmanRules", "[", "LagLRSM", "]"}]], "Input",
405 CellChangeTimes->{{3.680433078792603*^9, 3.6804330829390717`*^9}}]
406}, Open ]],
407
408Cell[CellGroupData[{
409
410Cell["UFO Output", "Section",
411 CellChangeTimes->{
412 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
413 3.59620768137986*^9, 3.596207684212212*^9}}],
414
415Cell[BoxData[
416 RowBox[{"(*", " ",
417 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
418 RowBox[{"appear", ":", " ",
419 RowBox[{"\"\<\!\(\*
420StyleBox[\"QN\", \"MessageName\"]\)\!\(\*
421StyleBox[\"::\", \"MessageName\"]\)\!\(\*
422StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"", "Warning"}], ":", " ",
423 RowBox[{
424 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
425 " ",
426 RowBox[{"encountered", "!"}],
427 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
428checks for this.\>"}]}]}], "*)"}]], "Input",
429 CellChangeTimes->{{3.656940935403554*^9, 3.656940952617087*^9}, {
430 3.6648662530247602`*^9, 3.664866256652151*^9}}],
431
432Cell[BoxData[
433 RowBox[{
434 RowBox[{"WriteUFO", "[",
435 RowBox[{"LagFull", ",", " ",
436 RowBox[{"Output", "\[Rule]", "\"\<EffLRSM_UFO\>\""}]}], "]"}],
437 ";"}]], "Input"]
438}, Open ]],
439
440Cell[CellGroupData[{
441
442Cell["On-Shell Renormalization with FeynArts", "Section",
443 CellChangeTimes->{
444 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
445 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199318805483*^9,
446 3.6621993287719097`*^9}, {3.662216149623213*^9, 3.662216170908692*^9}}],
447
448Cell[BoxData[
449 RowBox[{"(*", " ",
450 RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
451 RowBox[{"appear", ":", " ",
452 RowBox[{"\"\<\!\(\*
453StyleBox[\"QN\", \"MessageName\"]\)\!\(\*
454StyleBox[\"::\", \"MessageName\"]\)\!\(\*
455StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"", "Warning"}], ":", " ",
456 RowBox[{
457 "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
458 " ",
459 RowBox[{"encountered", "!"}],
460 "\"\<;\nThis is expected as lepton number is violated and FeynRules \
461checks for this.\>"}]}]}], "*)"}]], "Input"],
462
463Cell[BoxData[{
464 RowBox[{
465 RowBox[{"LRen", " ", "=", " ",
466 RowBox[{"OnShellRenormalization", "[",
467 RowBox[{"LagFull", ",",
468 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
469 RowBox[{"FlavorMixing", "\[Rule]", " ", "False"}], ",",
470 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}]}],
471 ";"}], "\[IndentingNewLine]",
472 RowBox[{"SetDirectory", "[", "\"\<~/.../FeynArts/FeynArts-3.8/Models\>\"",
473 "]"}], "\[IndentingNewLine]",
474 RowBox[{
475 RowBox[{"WriteFeynArtsOutput", "[",
476 RowBox[{"LRen", ",",
477 RowBox[{"GenericFile", "\[Rule]", "False"}], ",",
478 RowBox[{"FlavorExpand", "\[Rule]", "True"}], ",",
479 RowBox[{"Output", "\[Rule]", "\"\<EffLRSM_FA\>\""}]}], "]"}],
480 ";"}]}], "Input",
481 CellChangeTimes->{{3.686851598866076*^9, 3.686851599171583*^9}}]
482}, Open ]],
483
484Cell[CellGroupData[{
485
486Cell["Write Counter Terms to BSM.nlo", "Section",
487 CellChangeTimes->{
488 3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
489 3.59620768137986*^9, 3.596207684212212*^9}, {3.662199461128889*^9,
490 3.662199467953211*^9}, {3.662216541485611*^9, 3.662216544705964*^9}}],
491
492Cell[BoxData[
493 RowBox[{"(*", " ",
494 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
495 "*)"}]], "Input",
496 CellChangeTimes->{{3.664866502336069*^9, 3.664866536918194*^9}}],
497
498Cell[BoxData[
499 RowBox[{
500 RowBox[{"Quit", "[", "]"}], ";"}]], "Input"],
501
502Cell[BoxData[{
503 RowBox[{
504 RowBox[{"SetDirectory", "[", "\"\<~/.../FeynArts/FeynArts-3.8\>\"", "]"}],
505 ";"}], "\[IndentingNewLine]",
506 RowBox[{"<<", "FeynArts`"}], "\[IndentingNewLine]",
507 RowBox[{
508 RowBox[{"SetDirectory", "[", "\"\<~/.../FeynRules/feynrules-2_x\>\"", "]"}],
509 ";"}], "\[IndentingNewLine]",
510 RowBox[{"<<", "NLOCT`"}], "\[IndentingNewLine]",
511 RowBox[{
512 RowBox[{"$CurrentPath", "=",
513 RowBox[{"Directory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
514 RowBox[{
515 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
516 CellChangeTimes->{{3.6605244352671556`*^9, 3.660524496529204*^9}, {
517 3.660524555733032*^9, 3.660524556119355*^9}, 3.660528429455056*^9, {
518 3.6605286020387897`*^9, 3.660528602100609*^9}, {3.660530319063929*^9,
519 3.660530319160571*^9}, {3.660530366154963*^9, 3.66053037265125*^9}, {
520 3.660567620251087*^9, 3.6605676243565598`*^9}, {3.660568912072307*^9,
521 3.660568912133321*^9}, {3.66056999189439*^9, 3.660569993076742*^9}, {
522 3.660570552697504*^9, 3.660570554442918*^9}, {3.6605786705389338`*^9,
523 3.660578696714662*^9}, {3.660579471210005*^9, 3.660579471258161*^9}, {
524 3.660579558812049*^9, 3.6605795588460417`*^9}, {3.6605809812710743`*^9,
525 3.660580981737924*^9}, {3.662198567340746*^9, 3.662198568029338*^9}, {
526 3.664866430517086*^9, 3.664866431258605*^9}, {3.664866469872073*^9,
527 3.664866472738167*^9}, {3.6868516217227917`*^9, 3.686851625319178*^9}}],
528
529Cell[BoxData[
530 RowBox[{"(*",
531 RowBox[{
532 RowBox[{"Note", ":", " ",
533 RowBox[{"Grab", " ", "some", " ", "coffee"}]}], ";", " ",
534 RowBox[{"this", " ", "will", " ", "take", " ", "some", " ",
535 RowBox[{"minutes", ".", " ", "Output"}], " ", "is", " ", "located", " ",
536 "in", " ", "FeynRules", " ",
537 RowBox[{"directory", "."}]}]}], "*)"}]], "Input",
538 CellChangeTimes->{{3.660583519626089*^9, 3.660583556727778*^9}, {
539 3.6868513541484337`*^9, 3.686851354998522*^9}}],
540
541Cell[BoxData[
542 RowBox[{
543 RowBox[{"WriteCT", "[",
544 RowBox[{"\"\<EffLRSM_FA/EffLRSM_FA\>\"", ",", "\"\<Lorentz\>\"", ",",
545 RowBox[{"Output", "\[Rule]", "\"\<EffLRSMct\>\""}], ",",
546 "\[IndentingNewLine]",
547 RowBox[{"ZeroMom", "\[Rule]",
548 RowBox[{"{",
549 RowBox[{"{",
550 RowBox[{"aS", ",",
551 RowBox[{"{",
552 RowBox[{
553 RowBox[{"F", "[", "7", "]"}], ",",
554 RowBox[{"V", "[", "4", "]"}], ",",
555 RowBox[{"-",
556 RowBox[{"F", "[", "7", "]"}]}]}], "}"}]}], "}"}], "}"}]}], ",",
557 "\[IndentingNewLine]",
558 RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
559 RowBox[{"Exclude4ScalarsCT", "\[Rule]", "True"}]}], "]"}], ";"}]], "Input"]
560}, Open ]],
561
562Cell[CellGroupData[{
563
564Cell["\<\
565Load BSM.nlo and Generate UFO at NLO\
566\>", "Section",
567 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
568 3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
569 3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
570 3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
571 3.6550378926486187`*^9}, {3.660583469669636*^9, 3.6605834886500196`*^9}}],
572
573Cell[BoxData[
574 RowBox[{"(*", " ",
575 RowBox[{"quit", " ", "kernel", " ", "before", " ", "continuing"}], " ",
576 "*)"}]], "Input"],
577
578Cell[BoxData[
579 RowBox[{
580 RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
581 CellChangeTimes->{{3.660583502811754*^9, 3.6605835046945877`*^9},
582 3.66212390294204*^9}],
583
584Cell[BoxData[{
585 RowBox[{
586 RowBox[{"$CurrentPath", "=",
587 RowBox[{"Directory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
588 RowBox[{
589 RowBox[{"$FeynRulesPath", "=",
590 RowBox[{
591 "SetDirectory", "[", "\"\<~/.../FeynRules/feynrules-2_x\>\"", "]"}]}],
592 ";"}], "\[IndentingNewLine]",
593 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
594 RowBox[{
595 RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
596 CellChangeTimes->{{3.6648666749712143`*^9, 3.664866676189664*^9}, {
597 3.6868522686358957`*^9, 3.686852268940915*^9}}],
598
599Cell[BoxData[{
600 RowBox[{
601 RowBox[{"LoadModel", "[",
602 RowBox[{"\"\<sm.fr\>\"", ",", "\"\<effLRSM.fr\>\""}], "]"}],
603 ";"}], "\[IndentingNewLine]",
604 RowBox[{
605 RowBox[{"LoadRestriction", "[",
606 RowBox[{"\"\<massless.rst\>\"", ",", "\"\<diagonalMixing.rst\>\""}], "]"}],
607 ";"}], "\[IndentingNewLine]",
608 RowBox[{
609 RowBox[{"FeynmanGauge", "=", "True"}], ";"}]}], "Input",
610 CellChangeTimes->{{3.662827026079749*^9, 3.662827026713024*^9},
611 3.680433296257472*^9, {3.680455037101138*^9, 3.680455043996207*^9},
612 3.686848266936306*^9}],
613
614Cell[BoxData[{
615 RowBox[{
616 RowBox[{
617 "Get", "[", "\"\<~/.../FeynRules/feynrules-2_x/EffLRSMct.nlo\>\"", "]"}],
618 ";"}], "\[IndentingNewLine]",
619 RowBox[{
620 RowBox[{"WriteUFO", "[",
621 RowBox[{"LagFull", ",",
622 RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
623 RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",", " ",
624 RowBox[{"Output", "\[Rule]", "\"\<EffLRSM_NLO\>\""}]}], "]"}],
625 ";"}]}], "Input"]
626}, Open ]]
627}, Open ]]
628},
629WindowSize->{1330, 781},
630WindowMargins->{{Automatic, 23}, {Automatic, -8}},
631PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
632ShowSelection->True,
633FrontEndVersion->"9.0 for Linux x86 (64-bit) (November 20, 2012)",
634StyleDefinitions->"Default.nb"
635]
636(* End of Notebook Content *)
637
638(* Internal cache information *)
639(*CellTagsOutline
640CellTagsIndex->{}
641*)
642(*CellTagsIndex
643CellTagsIndex->{}
644*)
645(*NotebookFileOutline
646Notebook[{
647Cell[545, 20, 331, 7, 32, "Input"],
648Cell[879, 29, 1048, 24, 99, "Input"],
649Cell[1930, 55, 136, 3, 32, "Input"],
650Cell[CellGroupData[{
651Cell[2091, 62, 194, 3, 90, "Title"],
652Cell[2288, 67, 1489, 26, 99, "Input"],
653Cell[3780, 95, 382, 8, 32, "Input"],
654Cell[4165, 105, 591, 13, 77, "Input"],
655Cell[CellGroupData[{
656Cell[4781, 122, 393, 5, 80, "Section"],
657Cell[5177, 129, 149, 2, 32, "Input"],
658Cell[5329, 133, 100, 1, 32, "Input"],
659Cell[5432, 136, 249, 3, 32, "Input"],
660Cell[5684, 141, 162, 3, 32, "Input"],
661Cell[5849, 146, 347, 8, 55, "Input"]
662}, Open ]],
663Cell[CellGroupData[{
664Cell[6233, 159, 174, 2, 80, "Section"],
665Cell[CellGroupData[{
666Cell[6432, 165, 165, 3, 32, "Input"],
667Cell[CellGroupData[{
668Cell[6622, 172, 459, 6, 23, "Print"],
669Cell[7084, 180, 457, 6, 23, "Print"],
670Cell[7544, 188, 650, 13, 23, "Print"],
671Cell[8197, 203, 407, 5, 23, "Print"],
672Cell[8607, 210, 396, 5, 23, "Print"],
673Cell[9006, 217, 518, 9, 23, "Print"],
674Cell[9527, 228, 406, 5, 23, "Print"]
675}, Open ]],
676Cell[9948, 236, 392, 6, 32, "Output"]
677}, Open ]],
678Cell[CellGroupData[{
679Cell[10377, 247, 149, 2, 32, "Input"],
680Cell[CellGroupData[{
681Cell[10551, 253, 376, 8, 23, "Print"],
682Cell[10930, 263, 151, 2, 23, "Print"]
683}, Open ]],
684Cell[11096, 268, 120, 2, 32, "Output"]
685}, Open ]],
686Cell[CellGroupData[{
687Cell[11253, 275, 194, 3, 32, "Input"],
688Cell[CellGroupData[{
689Cell[11472, 282, 449, 9, 23, "Print"],
690Cell[11924, 293, 224, 3, 23, "Print"],
691Cell[12151, 298, 218, 3, 23, "Print"],
692Cell[12372, 303, 259, 4, 23, "Print"]
693}, Open ]],
694Cell[12646, 310, 894, 22, 102, "Output"]
695}, Open ]],
696Cell[CellGroupData[{
697Cell[13577, 337, 402, 6, 32, "Input"],
698Cell[CellGroupData[{
699Cell[14004, 347, 442, 9, 23, "Print"],
700Cell[14449, 358, 224, 3, 23, "Print"],
701Cell[14676, 363, 234, 3, 23, "Print"]
702}, Open ]],
703Cell[14925, 369, 191, 3, 32, "Output"]
704}, Open ]]
705}, Open ]],
706Cell[CellGroupData[{
707Cell[15165, 378, 283, 5, 80, "Section"],
708Cell[15451, 385, 720, 15, 55, "Input"],
709Cell[16174, 402, 139, 2, 32, "Input"]
710}, Open ]],
711Cell[CellGroupData[{
712Cell[16350, 409, 169, 3, 80, "Section"],
713Cell[16522, 414, 695, 15, 55, "Input"],
714Cell[17220, 431, 171, 5, 32, "Input"]
715}, Open ]],
716Cell[CellGroupData[{
717Cell[17428, 441, 295, 4, 80, "Section"],
718Cell[17726, 447, 578, 13, 55, "Input"],
719Cell[18307, 462, 790, 18, 77, "Input"]
720}, Open ]],
721Cell[CellGroupData[{
722Cell[19134, 485, 285, 4, 80, "Section"],
723Cell[19422, 491, 195, 4, 32, "Input"],
724Cell[19620, 497, 70, 2, 32, "Input"],
725Cell[19693, 501, 1436, 25, 143, "Input"],
726Cell[21132, 528, 480, 10, 32, "Input"],
727Cell[21615, 540, 705, 18, 77, "Input"]
728}, Open ]],
729Cell[CellGroupData[{
730Cell[22357, 563, 423, 7, 80, "Section"],
731Cell[22783, 572, 129, 3, 32, "Input"],
732Cell[22915, 577, 163, 4, 32, "Input"],
733Cell[23081, 583, 545, 13, 99, "Input"],
734Cell[23629, 598, 542, 13, 77, "Input"],
735Cell[24174, 613, 428, 11, 55, "Input"]
736}, Open ]]
737}, Open ]]
738}
739]
740*)
741
742(* End of internal cache information *)