EffLRSM: effLRSMnlo.nb

File effLRSMnlo.nb, 27.3 KB (added by richardphysics, 4 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 *)