StandardModel: SM.nb

File SM.nb, 28.8 KB (added by degrande, 5 years ago)
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[     29348,        888]
14NotebookOptionsPosition[     23926,        711]
15NotebookOutlinePosition[     24399,        730]
16CellTagsIndexPosition[     24356,        727]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"Quit", "[", "]"}]], "Input",
23 CellChangeTimes->{{3.4921467751527157`*^9, 3.492146776183146*^9}}],
24
25Cell[BoxData[
26 RowBox[{"$FeynRulesPath", "=",
27  RowBox[{
28  "SetDirectory", "[", "\"\<< address of the package >\>\"", "]"}]}]], "Input",\
29
30 CellChangeTimes->{{3.41265862251538*^9, 3.412658649947229*^9}, {
31  3.423415585782702*^9, 3.423415597189939*^9}, {3.4234163173467493`*^9,
32  3.4234163227881193`*^9}}],
33
34Cell[BoxData[
35 RowBox[{"<<", "FeynRules`"}]], "Input",
36 CellChangeTimes->{{3.547535564344927*^9, 3.547535564971527*^9}}],
37
38Cell[BoxData[
39 RowBox[{
40  RowBox[{"SetDirectory", "[",
41   RowBox[{"$FeynRulesPath", "<>", "\"\</Models/SM\>\""}], "]"}],
42  ";"}]], "Input",
43 CellChangeTimes->{{3.419073170860696*^9, 3.419073182827229*^9}}],
44
45Cell[CellGroupData[{
46
47Cell["The Standard Model", "Title"],
48
49Cell["We first load in the Standard Model model-file", "Text"],
50
51Cell[BoxData[
52 RowBox[{"LoadModel", "[", "\"\<SM.fr\>\"", "]"}]], "Input",
53 CellChangeTimes->{{3.4022069973481913`*^9, 3.402207011768662*^9},
54   3.4022081212072697`*^9, {3.402208250379383*^9, 3.402208254043104*^9},
55   3.4027466057480917`*^9, {3.403240270135737*^9, 3.403240277228945*^9},
56   3.403266503388291*^9, {3.403267649630335*^9, 3.40326765417397*^9}, {
57   3.403269919787421*^9, 3.403269921965273*^9}, {3.403347551273425*^9,
58   3.403347555049163*^9}, 3.4044490490588417`*^9, {3.411744339876704*^9,
59   3.411744340012457*^9}, {3.4121886924550533`*^9, 3.412188699157571*^9},
60   3.412188808811866*^9, 3.4121888580521603`*^9, {3.412450464077868*^9,
61   3.412450464378695*^9}, {3.413715097460478*^9, 3.41371509757642*^9}, {
62   3.41440825334604*^9, 3.414408254159686*^9}, {3.41862573831756*^9,
63   3.4186257392223186`*^9}, {3.4190731862389174`*^9, 3.419073187003003*^9}, {
64   3.419073336802393*^9, 3.4190733374513063`*^9}, {3.4190828041801767`*^9,
65   3.4190828048079023`*^9}, 3.542453089813714*^9, 3.547534567644828*^9}],
66
67Cell[BoxData[
68 RowBox[{"LoadRestriction", "[",
69  RowBox[{"\"\<Massless.rst\>\"", ",", "\"\<DiagonalCKM.rst\>\""}],
70  "]"}]], "Input",
71 CellChangeTimes->{{3.547535575308606*^9, 3.547535589236651*^9}}],
72
73Cell[CellGroupData[{
74
75Cell["The SM lagrangian", "Section",
76 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}}],
77
78Cell[CellGroupData[{
79
80Cell["Unitary Gauge", "Subsection",
81 CellChangeTimes->{{3.411910080775518*^9, 3.41191008414785*^9}}],
82
83Cell["\<\
84The full lagrangian in unitary gauge can be accessed immediatly via\
85\>", "Text",
86 CellChangeTimes->{{3.411910016771447*^9, 3.411910062087528*^9}, {
87  3.411910094087652*^9, 3.411910096590087*^9}, {3.411910215381351*^9,
88  3.4119102162371798`*^9}}],
89
90Cell[BoxData["LSM"], "Input",
91 CellChangeTimes->{{3.411910106011015*^9, 3.411910106351081*^9}}],
92
93Cell[CellGroupData[{
94
95Cell["Gauge sector", "Subsubsection",
96 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}}],
97
98Cell["\<\
99The part of the lagrangian representing the gauge sector can be accessed via\
100\>", "Text",
101 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
102  3.411910235813308*^9, 3.4119102509651337`*^9}}],
103
104Cell[BoxData["LGauge"], "Input",
105 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}}]
106}, Open  ]],
107
108Cell[CellGroupData[{
109
110Cell["Scalar sector", "Subsubsection",
111 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
112  3.411910274288354*^9, 3.41191027703699*^9}}],
113
114Cell["\<\
115The part of the lagrangian representing the scalar sector can be accessed via\
116\
117\>", "Text",
118 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
119  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.41191028690167*^9,
120  3.411910287709094*^9}}],
121
122Cell[BoxData["LHiggs"], "Input",
123 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
124  3.411910290186841*^9, 3.411910292005392*^9}}]
125}, Open  ]],
126
127Cell[CellGroupData[{
128
129Cell["Fermion sector", "Subsubsection",
130 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
131  3.4119103528868923`*^9, 3.411910357404933*^9}}],
132
133Cell["\<\
134The part of the lagrangian representing the fermion sector can be accessed via\
135\>", "Text",
136 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
137  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.4119103873434134`*^9,
138  3.411910388077118*^9}}],
139
140Cell[BoxData["LFermions"], "Input",
141 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
142  3.411910379677638*^9, 3.4119103807255*^9}}]
143}, Open  ]],
144
145Cell[CellGroupData[{
146
147Cell["Yukawa sector", "Subsubsection",
148 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
149  3.411910393742139*^9, 3.411910394852953*^9}}],
150
151Cell["\<\
152The part of the lagrangian representing the Yukawa sector can be accessed via\
153\
154\>", "Text",
155 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
156  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.411910399061659*^9,
157  3.411910400421185*^9}}],
158
159Cell[BoxData["LYukawa"], "Input",
160 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
161  3.4119104032372026`*^9, 3.411910404165523*^9}}]
162}, Open  ]]
163}, Open  ]],
164
165Cell[CellGroupData[{
166
167Cell["Feynman Gauge", "Subsection",
168 CellChangeTimes->{{3.411910080775518*^9, 3.41191008414785*^9}, {
169  3.411910431304595*^9, 3.41191043361451*^9}}],
170
171Cell["\<\
172By default, SM lagrangian is loaded in unitary gauge. However, Feynman gauge \
173can be obtained by putting\
174\>", "Text",
175 CellChangeTimes->{{3.411910440565042*^9, 3.411910480762147*^9}}],
176
177Cell[BoxData[
178 RowBox[{
179  RowBox[{"FeynmanGauge", "=", "True"}], ";"}]], "Input",
180 CellChangeTimes->{{3.411910487433539*^9, 3.4119104927114267`*^9}}],
181
182Cell["\<\
183The full lagrangian in Feynman gauge can be accessed immediatly via\
184\>", "Text",
185 CellChangeTimes->{{3.411910016771447*^9, 3.411910062087528*^9}, {
186  3.411910094087652*^9, 3.411910096590087*^9}, {3.411910215381351*^9,
187  3.4119102162371798`*^9}, {3.411910436148773*^9, 3.411910437013105*^9}}],
188
189Cell[BoxData["LSM"], "Input",
190 CellChangeTimes->{{3.411910106011015*^9, 3.411910106351081*^9}}],
191
192Cell[CellGroupData[{
193
194Cell["Gauge sector", "Subsubsection",
195 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}}],
196
197Cell["\<\
198The part of the lagrangian representing the gauge sector can be accessed via\
199\>", "Text",
200 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
201  3.411910235813308*^9, 3.4119102509651337`*^9}}],
202
203Cell[BoxData["LGauge"], "Input",
204 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}}]
205}, Open  ]],
206
207Cell[CellGroupData[{
208
209Cell["Scalar sector", "Subsubsection",
210 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
211  3.411910274288354*^9, 3.41191027703699*^9}}],
212
213Cell["\<\
214The part of the lagrangian representing the scalar sector can be accessed via\
215\
216\>", "Text",
217 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
218  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.41191028690167*^9,
219  3.411910287709094*^9}}],
220
221Cell[BoxData["LHiggs"], "Input",
222 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
223  3.411910290186841*^9, 3.411910292005392*^9}}]
224}, Open  ]],
225
226Cell[CellGroupData[{
227
228Cell["Fermion sector", "Subsubsection",
229 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
230  3.4119103528868923`*^9, 3.411910357404933*^9}}],
231
232Cell["\<\
233The part of the lagrangian representing the fermion sector can be accessed via\
234\>", "Text",
235 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
236  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.4119103873434134`*^9,
237  3.411910388077118*^9}}],
238
239Cell[BoxData["LFermions"], "Input",
240 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
241  3.411910379677638*^9, 3.4119103807255*^9}}]
242}, Open  ]],
243
244Cell[CellGroupData[{
245
246Cell["Yukawa sector", "Subsubsection",
247 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
248  3.411910393742139*^9, 3.411910394852953*^9}}],
249
250Cell["\<\
251The part of the lagrangian representing the Yukawa sector can be accessed via\
252\
253\>", "Text",
254 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
255  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.411910399061659*^9,
256  3.411910400421185*^9}}],
257
258Cell[BoxData["LYukawa"], "Input",
259 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
260  3.4119104032372026`*^9, 3.411910404165523*^9}}]
261}, Open  ]],
262
263Cell[CellGroupData[{
264
265Cell["Ghost sector", "Subsubsection",
266 CellChangeTimes->{{3.4119102216473637`*^9, 3.411910224608881*^9}, {
267  3.411910393742139*^9, 3.411910394852953*^9}, {3.4119106057673063`*^9,
268  3.411910606397141*^9}}],
269
270Cell["\<\
271The part of the lagrangian representing the ghost sector can be accessed via\
272\>", "Text",
273 CellChangeTimes->{{3.4119101255458603`*^9, 3.411910192225008*^9}, {
274  3.411910235813308*^9, 3.4119102509651337`*^9}, {3.411910399061659*^9,
275  3.411910400421185*^9}, {3.4119106110136833`*^9, 3.411910611661195*^9}}],
276
277Cell[BoxData["LGhost"], "Input",
278 CellChangeTimes->{{3.411910181169772*^9, 3.411910184667088*^9}, {
279  3.4119104032372026`*^9, 3.411910404165523*^9}, {3.411910615100849*^9,
280  3.411910615781582*^9}}]
281}, Open  ]]
282}, Open  ]]
283}, Closed]],
284
285Cell[CellGroupData[{
286
287Cell["Checking the Lagrangian", "Section",
288 CellChangeTimes->{{3.4234156496121893`*^9, 3.423415657027569*^9}}],
289
290Cell[BoxData[
291 RowBox[{
292  RowBox[{"FeynmanGauge", "=", "False"}], ";"}]], "Input",
293 CellChangeTimes->{{3.4234159280429783`*^9, 3.423415934798614*^9}, {
294  3.588417616520516*^9, 3.588417617062751*^9}}],
295
296Cell[CellGroupData[{
297
298Cell["Checking hermiticity", "Subsection",
299 CellChangeTimes->{{3.4234156831783743`*^9, 3.423415689540821*^9}}],
300
301Cell["The hermiticity of the Lagrangian can be checked via", "Text",
302 CellChangeTimes->{{3.423415694323595*^9, 3.423415707867804*^9}}],
303
304Cell[BoxData[
305 RowBox[{"CheckHermiticity", "[", "LSM", "]"}]], "Input",
306 CellChangeTimes->{{3.414763473002973*^9, 3.4147634845507174`*^9}, {
307  3.5884182415538187`*^9, 3.5884182416874647`*^9}}],
308
309Cell["\<\
310The same command can be applied to the flavor - expanded Lagrangian\
311\>", "Text",
312 CellChangeTimes->{{3.423415721582508*^9, 3.423415739787958*^9}}],
313
314Cell[BoxData[
315 RowBox[{"CheckHermiticity", "[",
316  RowBox[{"LSM", ",",
317   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]], "Input",
318 CellChangeTimes->{{3.414763604111637*^9, 3.4147636044973803`*^9}}]
319}, Open  ]],
320
321Cell[CellGroupData[{
322
323Cell["Checking the mass spectrum", "Subsection",
324 CellChangeTimes->{{3.423415831209465*^9, 3.423415835278742*^9}}],
325
326Cell["\<\
327FeynRules can check whether all the values of the masses given in the model \
328file ar ein agreement with the masses given in the Lagrangian:\
329\>", "Text",
330 CellChangeTimes->{{3.423415751305388*^9, 3.423415792594919*^9}}],
331
332Cell[BoxData[
333 RowBox[{"CheckMassSpectrum", "[", "LSM", "]"}]], "Input",
334 CellChangeTimes->{{3.414736325829801*^9, 3.414736326065997*^9},
335   3.4147528281460342`*^9}]
336}, Open  ]],
337
338Cell[CellGroupData[{
339
340Cell["Checking the kinetic terms", "Subsection",
341 CellChangeTimes->{{3.423415843532096*^9, 3.42341584850035*^9}}],
342
343Cell["\<\
344In the same way, FeynRules can check whether all kinetic terms are correctly \
345normalized :\
346\>", "Text",
347 CellChangeTimes->{{3.423415798800013*^9, 3.423415817964087*^9}}],
348
349Cell[BoxData[
350 RowBox[{"CheckKineticTermNormalisation", "[",
351  RowBox[{"LSM", ",",
352   RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}]}], "]"}]], "Input",
353 CellChangeTimes->{{3.4147636651636972`*^9, 3.414763667250185*^9}}],
354
355Cell["\<\
356The same command cen be used at the level of the flavor - expanded Lagrangian \
357:\
358\>", "Text",
359 CellChangeTimes->{{3.423415868906054*^9, 3.4234158864202538`*^9}}],
360
361Cell[BoxData[
362 RowBox[{"CheckKineticTermNormalisation", "[",
363  RowBox[{"LSM", ",",
364   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]], "Input"]
365}, Open  ]]
366}, Closed]],
367
368Cell[CellGroupData[{
369
370Cell["Calculation of the Feynman rules", "Section",
371 CellChangeTimes->{{3.411910637332127*^9, 3.41191064360071*^9}}],
372
373Cell["\<\
374We will calculate now explicitly the Feynman rules for the SM in Feynman \
375gauge.\
376\>", "Text",
377 CellChangeTimes->{{3.41191064804084*^9, 3.411910672936666*^9}}],
378
379Cell[BoxData[
380 RowBox[{
381  RowBox[{"FeynmanGauge", "=", "True"}], ";"}]], "Input",
382 CellChangeTimes->{{3.411910674121689*^9, 3.4119106799564953`*^9}}],
383
384Cell[CellGroupData[{
385
386Cell["Gauge sector", "Subsection",
387 CellChangeTimes->{{3.411910690443768*^9, 3.411910692233274*^9}}],
388
389Cell[BoxData[
390 RowBox[{
391  RowBox[{"vertsGauge", "=",
392   RowBox[{"FeynmanRules", "[",
393    RowBox[{"LGauge", ",",
394     RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}]}], "]"}]}], ";"}]], "Input",
395 CellChangeTimes->{
396  3.40274690437012*^9, {3.402746954727231*^9, 3.4027469691642714`*^9}, {
397   3.4032592998045692`*^9, 3.403259304564069*^9}, 3.403267976055098*^9, {
398   3.411910705212987*^9, 3.411910705553741*^9}, {3.423415993272615*^9,
399   3.42341600147605*^9}}]
400}, Open  ]],
401
402Cell[CellGroupData[{
403
404Cell["Higgs Sector", "Subsection",
405 CellChangeTimes->{3.404470070069872*^9}],
406
407Cell[BoxData[
408 RowBox[{
409  RowBox[{"vertsHiggs", "=",
410   RowBox[{"FeynmanRules", "[", "LHiggs", "]"}]}], ";"}]], "Input",
411 CellChangeTimes->{{3.411910728609757*^9, 3.41191072908109*^9}, {
412   3.423416009086747*^9, 3.423416012982971*^9}, 3.423416065626593*^9}]
413}, Open  ]],
414
415Cell[CellGroupData[{
416
417Cell["Matter sector", "Subsection"],
418
419Cell[BoxData[
420 RowBox[{
421  RowBox[{"vertsFermions", "=",
422   RowBox[{"FeynmanRules", "[", "LFermions", "]"}]}], ";"}]], "Input",
423 CellChangeTimes->{{3.402205118991502*^9, 3.4022051228508387`*^9},
424   3.402747846432675*^9, {3.4032613449872217`*^9, 3.403261349245253*^9}, {
425   3.403266977133403*^9, 3.4032669793125563`*^9}, {3.411910788644576*^9,
426   3.411910789233307*^9}, {3.423416034294303*^9, 3.423416042739814*^9}}]
427}, Open  ]],
428
429Cell[CellGroupData[{
430
431Cell["Yukawa sector", "Subsection"],
432
433Cell[BoxData[
434 RowBox[{
435  RowBox[{"vertsYukawa", "=",
436   RowBox[{"FeynmanRules", "[",
437    RowBox[{"LYukawa", ",",
438     RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]}], ";"}]], "Input",
439 CellChangeTimes->{{3.403266987858165*^9, 3.4032670026080303`*^9}, {
440  3.411910826105542*^9, 3.411910826777452*^9}, {3.423416089772278*^9,
441  3.423416095643961*^9}}]
442}, Open  ]],
443
444Cell[CellGroupData[{
445
446Cell["LGhost sector", "Subsection",
447 CellChangeTimes->{{3.411910849714159*^9, 3.411910851088325*^9}}],
448
449Cell[BoxData[
450 RowBox[{
451  RowBox[{"vertsGhosts", "=",
452   RowBox[{"FeynmanRules", "[",
453    RowBox[{"LGhost", ",",
454     RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}]}], "]"}]}], ";"}]], "Input",
455 CellChangeTimes->{{3.403266987858165*^9, 3.4032670026080303`*^9}, {
456   3.411910826105542*^9, 3.411910826777452*^9}, {3.411910857580064*^9,
457   3.411910870001601*^9}, 3.415256462464088*^9, {3.423416107236608*^9,
458   3.423416112835993*^9}}]
459}, Open  ]]
460}, Closed]],
461
462Cell[CellGroupData[{
463
464Cell["Outputs and interfaces", "Section",
465 CellChangeTimes->{
466  3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}}],
467
468Cell[CellGroupData[{
469
470Cell["FeynArts output", "Subsection",
471 CellChangeTimes->{{3.411911000132907*^9, 3.411911005135379*^9}}],
472
473Cell["The FeynArts output for the SM can be obtained via", "Text",
474 CellChangeTimes->{{3.411910955110715*^9, 3.4119109690877657`*^9}, {
475  3.411911062212926*^9, 3.4119110634765463`*^9}}],
476
477Cell[BoxData[{
478 RowBox[{
479  RowBox[{
480   RowBox[{"FeynmanGauge", " ", "=", " ", "False"}], ";"}],
481  "\[IndentingNewLine]"}], "\[IndentingNewLine]",
482 RowBox[{"WriteFeynArtsOutput", "[",
483  RowBox[{"LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",",
484   RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}]}], "]"}]}], "Input",
485 CellChangeTimes->{{3.411911018354619*^9, 3.4119110196967573`*^9}, {
486   3.412185470680571*^9, 3.412185480912915*^9}, {3.41245048911446*^9,
487   3.4124504892674522`*^9}, 3.414315634453483*^9, {3.423416180751266*^9,
488   3.423416182892613*^9}}],
489
490Cell["FeynArts also supports the Feynman gauge", "Text",
491 CellChangeTimes->{{3.4121854851672792`*^9, 3.412185494183435*^9}}],
492
493Cell[BoxData[{
494 RowBox[{
495  RowBox[{
496   RowBox[{"FeynmanGauge", " ", "=", " ", "True"}], ";"}],
497  "\[IndentingNewLine]"}], "\[IndentingNewLine]",
498 RowBox[{"WriteFeynArtsOutput", "[",
499  RowBox[{
500  "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LGhost",
501   ",",
502   RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}], ",",
503   RowBox[{"Output", "\[Rule]", "\"\<SM.mod\>\""}]}], "]"}]}], "Input",
504 CellChangeTimes->{{3.411911018354619*^9, 3.4119110196967573`*^9}, {
505  3.412185470680571*^9, 3.412185506240796*^9}, {3.4234161971513453`*^9,
506  3.423416209716672*^9}}]
507}, Open  ]],
508
509Cell[CellGroupData[{
510
511Cell["Sherpa output", "Subsection",
512 CellChangeTimes->{{3.4119110460271273`*^9, 3.411911049405038*^9}}],
513
514Cell["The Sherpa output for the SM can be obtained via", "Text",
515 CellChangeTimes->{{3.411910955110715*^9, 3.4119109690877657`*^9}, {
516  3.411911059590749*^9, 3.4119110601885567`*^9}}],
517
518Cell[BoxData[{
519 RowBox[{
520  RowBox[{
521   RowBox[{"FeynmanGauge", "=", "False"}], ";"}],
522  "\[IndentingNewLine]"}], "\[IndentingNewLine]",
523 RowBox[{"WriteSHOutput", "[",
524  RowBox[{"LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa"}],
525  "]"}]}], "Input",
526 CellChangeTimes->{{3.4119110943383636`*^9, 3.41191110054212*^9}, {
527  3.412185455856553*^9, 3.412185456628922*^9}, {3.4143156544792233`*^9,
528  3.414315668290971*^9}, {3.423416222700699*^9, 3.423416235756197*^9}}]
529}, Open  ]],
530
531Cell[CellGroupData[{
532
533Cell["CalcHep output", "Subsection",
534 CellChangeTimes->{{3.4119110460271273`*^9, 3.411911049405038*^9}, {
535  3.411911150323291*^9, 3.4119111515122623`*^9}}],
536
537Cell["CalcHep also supports the Feynman gauge", "Text",
538 CellChangeTimes->{{3.4121853722243843`*^9, 3.412185380407339*^9}}],
539
540Cell[BoxData[{
541 RowBox[{
542  RowBox[{
543   RowBox[{"FeynmanGauge", "=", "True"}], ";"}],
544  "\[IndentingNewLine]"}], "\[IndentingNewLine]",
545 RowBox[{
546  RowBox[{"WriteCHOutput", "[",
547   RowBox[{
548   "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LGhost"}],
549    "]"}], ";"}]}], "Input",
550 CellChangeTimes->{{3.4119110943383636`*^9, 3.41191110054212*^9}, {
551  3.412185409690489*^9, 3.4121854275198507`*^9}, {3.412450479942018*^9,
552  3.4124504817138367`*^9}}],
553
554Cell[BoxData[""], "Input",
555 CellChangeTimes->{3.412185367693318*^9, 3.412185407947097*^9}]
556}, Open  ]],
557
558Cell[CellGroupData[{
559
560Cell["Whizard output", "Subsection",
561 CellChangeTimes->{{3.4119110460271273`*^9, 3.411911049405038*^9}, {
562  3.411911150323291*^9, 3.4119111515122623`*^9}, {3.5074389822322817`*^9,
563  3.507438983565652*^9}}],
564
565Cell["Whizard also supports the Feynman gauge", "Text",
566 CellChangeTimes->{{3.4121853722243843`*^9, 3.412185380407339*^9}, {
567  3.507438990997637*^9, 3.507438992828805*^9}}],
568
569Cell[BoxData[{
570 RowBox[{
571  RowBox[{
572   RowBox[{"FeynmanGauge", "=", "True"}], ";"}],
573  "\[IndentingNewLine]"}], "\[IndentingNewLine]",
574 RowBox[{
575  RowBox[{"WriteWOOutput", "[",
576   RowBox[{
577   "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LGhost"}],
578    "]"}], ";"}]}], "Input",
579 CellChangeTimes->{{3.4119110943383636`*^9, 3.41191110054212*^9}, {
580  3.412185409690489*^9, 3.4121854275198507`*^9}, {3.412450479942018*^9,
581  3.4124504817138367`*^9}, {3.507438986168689*^9, 3.507438987383316*^9}}],
582
583Cell[BoxData[""], "Input",
584 CellChangeTimes->{3.412185367693318*^9, 3.412185407947097*^9}]
585}, Open  ]],
586
587Cell[CellGroupData[{
588
589Cell["UFO output", "Subsection",
590 CellChangeTimes->{{3.4119110460271273`*^9, 3.411911049405038*^9}, {
591  3.411911150323291*^9, 3.4119111515122623`*^9}, {3.5074389822322817`*^9,
592  3.507438983565652*^9}, {3.507439158477231*^9, 3.507439159035309*^9}}],
593
594Cell[BoxData[
595 RowBox[{
596  RowBox[{"WriteUFO", "[",
597   RowBox[{
598   "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LGhost"}],
599    "]"}], ";"}]], "Input",
600 CellChangeTimes->{
601  3.412185367693318*^9, 3.412185407947097*^9, {3.5074391727475023`*^9,
602   3.5074391804929533`*^9}}]
603}, Open  ]]
604}, Closed]],
605
606Cell[CellGroupData[{
607
608Cell["NLO (QCD)", "Subsubtitle",
609 CellChangeTimes->{{3.606568224887179*^9, 3.606568225228436*^9}, {
610  3.606568338418679*^9, 3.606568340840477*^9}}],
611
612Cell[BoxData[
613 RowBox[{"(*",
614  RowBox[{
615  "Renormalization", " ", "and", " ", "output", " ", "to", " ", "FA"}],
616  "*)"}]], "Input",
617 CellChangeTimes->{{3.606568280586434*^9, 3.606568289181669*^9}, {
618  3.606568347601158*^9, 3.606568356208818*^9}}],
619
620Cell[BoxData[{
621 RowBox[{
622  RowBox[{"Lren", "=",
623   RowBox[{"OnShellRenormalization", "[",
624    RowBox[{"LSM", ",",
625     RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
626     RowBox[{"FlavorMixing", "\[Rule]", "False"}]}], "]"}]}],
627  ";"}], "\[IndentingNewLine]",
628 RowBox[{
629  RowBox[{"SetDirectory", "[", "\"\<~/celine/FeynArts-3.7/Models\>\"", "]"}],
630  ";"}], "\[IndentingNewLine]",
631 RowBox[{"WriteFeynArtsOutput", "[",
632  RowBox[{"Lren", ",",
633   RowBox[{"Output", "\[Rule]", "\"\<SMQCDrenoL\>\""}], ",",
634   RowBox[{"GenericFile", "\[Rule]", "False"}]}], "]"}]}], "Input",
635 CellChangeTimes->{{3.606568250159355*^9, 3.6065682552288523`*^9}, {
636  3.6065683087252502`*^9, 3.6065683093453608`*^9}, {3.6065684035728207`*^9,
637  3.606568404147846*^9}}],
638
639Cell[BoxData[
640 RowBox[{"(*",
641  RowBox[{"Computation", " ", "of", " ", "the", " ", "Counterterms"}],
642  "*)"}]], "Input",
643 CellChangeTimes->{{3.606568413917823*^9, 3.6065684261433983`*^9}}],
644
645Cell[BoxData[
646 RowBox[{"Quit", "[", "]"}]], "Input",
647 CellChangeTimes->{{3.606568433162343*^9, 3.606568438694559*^9}}],
648
649Cell[BoxData[{
650 RowBox[{
651  RowBox[{"SetDirectory", "[", "\"\<< FeynArts directory >\>\"", "]"}],
652  ";"}], "\[IndentingNewLine]",
653 RowBox[{"<<", " ", "FeynArts`"}], "\[IndentingNewLine]",
654 RowBox[{
655  RowBox[{"SetDirectory", "[", "\"\<< NLOCT directory >\>\"", "]"}],
656  ";"}], "\[IndentingNewLine]",
657 RowBox[{"<<", " ", "NLOCT`"}]}], "Input",
658 CellChangeTimes->{{3.6065684720815067`*^9, 3.6065685038224916`*^9}}],
659
660Cell[BoxData[
661 RowBox[{
662  RowBox[{"WriteCT", "[",
663   RowBox[{
664   "\"\<SMQCDrenoL/SMQCDrenoL\>\"", ",", "\"\<Lorentz\>\"", ",",
665    "\"\<SMQCDreno\>\"", ",",
666    RowBox[{"LabelInternal", "\[Rule]", "True"}], ",",
667    RowBox[{"QCDOnly", "\[Rule]", "True"}], ",",
668    RowBox[{"KeptIndices", "\[Rule]",
669     RowBox[{"{", "}"}]}], ",",
670    RowBox[{"ZeroMom", "\[Rule]",
671     RowBox[{"{",
672      RowBox[{"{",
673       RowBox[{"aS", ",",
674        RowBox[{"{",
675         RowBox[{
676          RowBox[{"F", "[", "7", "]"}], ",",
677          RowBox[{"V", "[", "4", "]"}], ",",
678          RowBox[{"-",
679           RowBox[{"F", "[", "7", "]"}]}]}], "}"}], ",", "0"}], "}"}],
680      "}"}]}], ",",
681    RowBox[{"ComplexMass", "\[Rule]", "False"}]}], "]"}], "//", "Timing",
682  " "}]], "Input",
683 CellChangeTimes->{{3.606569782037321*^9, 3.606569798425932*^9}}],
684
685Cell[BoxData[
686 RowBox[{"(*",
687  RowBox[{
688  "Quit", " ", "and", " ", "reload", " ", "FeynRules", " ", "and", " ", "the",
689    " ", "model", " ", "with", " ", "the", " ", "restrictions"}],
690  "*)"}]], "Input",
691 CellChangeTimes->{{3.6065685476941853`*^9, 3.6065685843784647`*^9}}],
692
693Cell[BoxData[{
694 RowBox[{
695  RowBox[{
696  "SetDirectory", "[",
697   "\"\<~/celine/feynrules/trunk/feynrules-development/R2\>\"", "]"}],
698  ";"}], "\[IndentingNewLine]",
699 RowBox[{
700  RowBox[{"Get", "[", "\"\<SMQCDreno.nlo\>\"", "]"}], ";"}]}], "Input"],
701
702Cell[BoxData[
703 RowBox[{"WriteUFO", "[",
704  RowBox[{"LSM", ",",
705   RowBox[{"UVCounterterms", "\[Rule]", "UV$vertlist"}], ",",
706   RowBox[{"R2Vertices", "\[Rule]", "R2$vertlist"}], ",",
707   RowBox[{"Output", "\[Rule]", "\"\<SM_NLO\>\""}]}], "]"}]], "Input",
708 CellChangeTimes->{{3.60656864069383*^9, 3.606568704839243*^9}}]
709}, Open  ]]
710}, Open  ]]
711},
712WindowSize->{1276, 756},
713WindowMargins->{{0, Automatic}, {Automatic, 0}},
714ShowSelection->True,
715Magnification:>FEPrivate`If[
716  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 2., 2. Inherited],
717FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 6, \
7182010)",
719StyleDefinitions->"Default.nb"
720]
721(* End of Notebook Content *)
722
723(* Internal cache information *)
724(*CellTagsOutline
725CellTagsIndex->{}
726*)
727(*CellTagsIndex
728CellTagsIndex->{}
729*)
730(*NotebookFileOutline
731Notebook[{
732Cell[545, 20, 120, 2, 54, "Input"],
733Cell[668, 24, 306, 7, 54, "Input"],
734Cell[977, 33, 120, 2, 54, "Input"],
735Cell[1100, 37, 207, 5, 54, "Input"],
736Cell[CellGroupData[{
737Cell[1332, 46, 35, 0, 152, "Title"],
738Cell[1370, 48, 62, 0, 52, "Text"],
739Cell[1435, 50, 1026, 14, 54, "Input"],
740Cell[2464, 66, 201, 4, 54, "Input"],
741Cell[CellGroupData[{
742Cell[2690, 74, 102, 1, 134, "Section"],
743Cell[CellGroupData[{
744Cell[2817, 79, 100, 1, 136, "Subsection"],
745Cell[2920, 82, 257, 5, 104, "Text"],
746Cell[3180, 89, 95, 1, 108, "Input"],
747Cell[CellGroupData[{
748Cell[3300, 94, 105, 1, 96, "Subsubsection"],
749Cell[3408, 97, 219, 4, 104, "Text"],
750Cell[3630, 103, 98, 1, 108, "Input"]
751}, Open  ]],
752Cell[CellGroupData[{
753Cell[3765, 109, 154, 2, 96, "Subsubsection"],
754Cell[3922, 113, 270, 6, 104, "Text"],
755Cell[4195, 121, 147, 2, 108, "Input"]
756}, Open  ]],
757Cell[CellGroupData[{
758Cell[4379, 128, 158, 2, 96, "Subsubsection"],
759Cell[4540, 132, 272, 5, 104, "Text"],
760Cell[4815, 139, 148, 2, 108, "Input"]
761}, Open  ]],
762Cell[CellGroupData[{
763Cell[5000, 146, 155, 2, 96, "Subsubsection"],
764Cell[5158, 150, 271, 6, 104, "Text"],
765Cell[5432, 158, 150, 2, 108, "Input"]
766}, Open  ]]
767}, Open  ]],
768Cell[CellGroupData[{
769Cell[5631, 166, 148, 2, 136, "Subsection"],
770Cell[5782, 170, 195, 4, 104, "Text"],
771Cell[5980, 176, 149, 3, 108, "Input"],
772Cell[6132, 181, 303, 5, 104, "Text"],
773Cell[6438, 188, 95, 1, 108, "Input"],
774Cell[CellGroupData[{
775Cell[6558, 193, 105, 1, 96, "Subsubsection"],
776Cell[6666, 196, 219, 4, 104, "Text"],
777Cell[6888, 202, 98, 1, 108, "Input"]
778}, Open  ]],
779Cell[CellGroupData[{
780Cell[7023, 208, 154, 2, 96, "Subsubsection"],
781Cell[7180, 212, 270, 6, 104, "Text"],
782Cell[7453, 220, 147, 2, 108, "Input"]
783}, Open  ]],
784Cell[CellGroupData[{
785Cell[7637, 227, 158, 2, 96, "Subsubsection"],
786Cell[7798, 231, 272, 5, 104, "Text"],
787Cell[8073, 238, 148, 2, 108, "Input"]
788}, Open  ]],
789Cell[CellGroupData[{
790Cell[8258, 245, 155, 2, 96, "Subsubsection"],
791Cell[8416, 249, 271, 6, 104, "Text"],
792Cell[8690, 257, 150, 2, 108, "Input"]
793}, Open  ]],
794Cell[CellGroupData[{
795Cell[8877, 264, 205, 3, 96, "Subsubsection"],
796Cell[9085, 269, 316, 5, 104, "Text"],
797Cell[9404, 276, 198, 3, 108, "Input"]
798}, Open  ]]
799}, Open  ]]
800}, Closed]],
801Cell[CellGroupData[{
802Cell[9663, 286, 110, 1, 70, "Section"],
803Cell[9776, 289, 199, 4, 108, "Input"],
804Cell[CellGroupData[{
805Cell[10000, 297, 110, 1, 136, "Subsection"],
806Cell[10113, 300, 134, 1, 104, "Text"],
807Cell[10250, 303, 192, 3, 108, "Input"],
808Cell[10445, 308, 157, 3, 104, "Text"],
809Cell[10605, 313, 206, 4, 108, "Input"]
810}, Open  ]],
811Cell[CellGroupData[{
812Cell[10848, 322, 114, 1, 136, "Subsection"],
813Cell[10965, 325, 230, 4, 164, "Text"],
814Cell[11198, 331, 166, 3, 108, "Input"]
815}, Open  ]],
816Cell[CellGroupData[{
817Cell[11401, 339, 113, 1, 136, "Subsection"],
818Cell[11517, 342, 181, 4, 104, "Text"],
819Cell[11701, 348, 219, 4, 108, "Input"],
820Cell[11923, 354, 173, 4, 104, "Text"],
821Cell[12099, 360, 151, 3, 108, "Input"]
822}, Open  ]]
823}, Closed]],
824Cell[CellGroupData[{
825Cell[12299, 369, 116, 1, 70, "Section"],
826Cell[12418, 372, 170, 4, 104, "Text"],
827Cell[12591, 378, 149, 3, 108, "Input"],
828Cell[CellGroupData[{
829Cell[12765, 385, 100, 1, 136, "Subsection"],
830Cell[12868, 388, 457, 10, 108, "Input"]
831}, Open  ]],
832Cell[CellGroupData[{
833Cell[13362, 403, 76, 1, 136, "Subsection"],
834Cell[13441, 406, 258, 5, 108, "Input"]
835}, Open  ]],
836Cell[CellGroupData[{
837Cell[13736, 416, 35, 0, 136, "Subsection"],
838Cell[13774, 418, 417, 7, 108, "Input"]
839}, Open  ]],
840Cell[CellGroupData[{
841Cell[14228, 430, 35, 0, 136, "Subsection"],
842Cell[14266, 432, 360, 8, 108, "Input"]
843}, Open  ]],
844Cell[CellGroupData[{
845Cell[14663, 445, 101, 1, 136, "Subsection"],
846Cell[14767, 448, 433, 9, 108, "Input"]
847}, Open  ]]
848}, Closed]],
849Cell[CellGroupData[{
850Cell[15249, 463, 132, 2, 70, "Section"],
851Cell[CellGroupData[{
852Cell[15406, 469, 103, 1, 136, "Subsection"],
853Cell[15512, 472, 185, 2, 104, "Text"],
854Cell[15700, 476, 569, 11, 290, "Input"],
855Cell[16272, 489, 124, 1, 104, "Text"],
856Cell[16399, 492, 575, 13, 290, "Input"]
857}, Open  ]],
858Cell[CellGroupData[{
859Cell[17011, 510, 103, 1, 136, "Subsection"],
860Cell[17117, 513, 183, 2, 104, "Text"],
861Cell[17303, 517, 474, 10, 230, "Input"]
862}, Open  ]],
863Cell[CellGroupData[{
864Cell[17814, 532, 155, 2, 272, "Subsection"],
865Cell[17972, 536, 123, 1, 208, "Text"],
866Cell[18098, 539, 468, 12, 460, "Input"],
867Cell[18569, 553, 90, 1, 216, "Input"]
868}, Open  ]],
869Cell[CellGroupData[{
870Cell[18696, 559, 206, 3, 272, "Subsection"],
871Cell[18905, 564, 172, 2, 208, "Text"],
872Cell[19080, 568, 514, 12, 460, "Input"],
873Cell[19597, 582, 90, 1, 216, "Input"]
874}, Open  ]],
875Cell[CellGroupData[{
876Cell[19724, 588, 248, 3, 272, "Subsection"],
877Cell[19975, 593, 291, 8, 216, "Input"]
878}, Open  ]]
879}, Closed]],
880Cell[CellGroupData[{
881Cell[20315, 607, 147, 2, 52, "Subsubtitle"],
882Cell[20465, 611, 248, 6, 54, "Input"],
883Cell[20716, 619, 749, 17, 115, "Input"],
884Cell[21468, 638, 189, 4, 54, "Input"],
885Cell[21660, 644, 118, 2, 54, "Input"],
886Cell[21781, 648, 415, 9, 145, "Input"],
887Cell[22199, 659, 847, 23, 145, "Input"],
888Cell[23049, 684, 276, 6, 54, "Input"],
889Cell[23328, 692, 246, 7, 85, "Input"],
890Cell[23577, 701, 321, 6, 85, "Input"]
891}, Open  ]]
892}, Open  ]]
893}
894]
895*)
896
897(* End of internal cache information *)