AnomalousGaugeCoupling: quartic_ls01.nb

File quartic_ls01.nb, 10.1 KB (added by eboli, 5 years ago)

Notebook used to generate the SM, LS0, and LS1 interactions.

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[     10156,        277]
14NotebookOptionsPosition[      8876,        235]
15NotebookOutlinePosition[      9273,        253]
16CellTagsIndexPosition[      9230,        250]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"\[IndentingNewLine]",
23  FrameBox[
24   RowBox[{
25   "\[SelectionPlaceholder]", " ", "FeynRules", " ",
26    "initialization"}]]}]], "Input",
27 CellChangeTimes->{{3.5542065470381117`*^9, 3.554206578929636*^9}}],
28
29Cell[BoxData[
30 RowBox[{"$FeynRulesPath", "=",
31  RowBox[{
32  "SetDirectory", "[",
33   "\"\</Users/oscar/Trabalho/WorkInProgress/FeynRules/feynrules-development\>\
34\"", "]"}]}]], "Input",
35 CellChangeTimes->{{3.41265862251538*^9, 3.412658649947229*^9}, {
36  3.423415585782702*^9, 3.423415597189939*^9}, {3.4234163173467493`*^9,
37  3.4234163227881193`*^9}, {3.549259989465773*^9, 3.549259999606881*^9}}],
38
39Cell[BoxData[
40 RowBox[{"<<", "FeynRules`"}]], "Input"],
41
42Cell[BoxData[
43 RowBox[{
44  RowBox[{"SetDirectory", "[",
45   RowBox[{"$FeynRulesPath", "<>", "\"\</Models/Quartic_linear\>\""}], "]"}],
46  ";"}]], "Input",
47 CellChangeTimes->{{3.419073170860696*^9, 3.419073182827229*^9}, {
48  3.549260010957612*^9, 3.549260018221593*^9}, {3.5527968095887547`*^9,
49  3.552796815826297*^9}}],
50
51Cell[BoxData[
52 RowBox[{"\[IndentingNewLine]",
53  FrameBox[
54   RowBox[{
55    RowBox[{
56    "\[SelectionPlaceholder]", " ", "Loading", " ", "the", " ", "file", " ",
57     "containing", " ", "the", " ", "Standard", " ", "Model", " ", "and", " ",
58      "L_S0"}], " ", "+", " ", "L_S1"}]], " "}]], "Input",
59 CellChangeTimes->{{3.554206586417383*^9, 3.554206720572693*^9}}],
60
61Cell[BoxData[
62 RowBox[{"LoadModel", "[", "\"\<quartic.fr\>\"", "]"}]], "Input",
63 CellChangeTimes->{{3.4022069973481913`*^9, 3.402207011768662*^9},
64   3.4022081212072697`*^9, {3.402208250379383*^9, 3.402208254043104*^9},
65   3.4027466057480917`*^9, {3.403240270135737*^9, 3.403240277228945*^9},
66   3.403266503388291*^9, {3.403267649630335*^9, 3.40326765417397*^9}, {
67   3.403269919787421*^9, 3.403269921965273*^9}, {3.403347551273425*^9,
68   3.403347555049163*^9}, 3.4044490490588417`*^9, {3.411744339876704*^9,
69   3.411744340012457*^9}, {3.4121886924550533`*^9, 3.412188699157571*^9},
70   3.412188808811866*^9, 3.4121888580521603`*^9, {3.412450464077868*^9,
71   3.412450464378695*^9}, {3.413715097460478*^9, 3.41371509757642*^9}, {
72   3.41440825334604*^9, 3.414408254159686*^9}, {3.41862573831756*^9,
73   3.4186257392223186`*^9}, {3.4190731862389174`*^9, 3.419073187003003*^9}, {
74   3.419073336802393*^9, 3.4190733374513063`*^9}, {3.4190828041801767`*^9,
75   3.4190828048079023`*^9}, {3.549356567929675*^9, 3.549356569284753*^9}, {
76   3.552796819881452*^9, 3.552796823177125*^9}}],
77
78Cell[BoxData[
79 RowBox[{"LoadRestriction", "[",
80  RowBox[{"\"\<Cabibbo.rst\>\"", ",", "\"\<Massless.rst\>\""}],
81  "]"}]], "Input",
82 CellChangeTimes->{{3.453790588503598*^9, 3.4537906008856297`*^9},
83   3.5529941847635183`*^9}],
84
85Cell[BoxData[
86 RowBox[{
87  FrameBox[
88   RowBox[{
89   "\[SelectionPlaceholder]", " ", "Checking", " ", "the", " ",
90    "lagrangian"}]], "\[IndentingNewLine]"}]], "Input",
91 CellChangeTimes->{{3.554207149082877*^9, 3.5542071697416487`*^9}}],
92
93Cell[BoxData[
94 RowBox[{
95  RowBox[{"FeynmanGauge", "=", "True"}], ";"}]], "Input",
96 CellChangeTimes->{{3.4234159280429783`*^9, 3.423415934798614*^9},
97   3.554210133339581*^9}],
98
99Cell["The hermiticity of the Lagrangian can be checked via", "Text",
100 CellChangeTimes->{{3.423415694323595*^9, 3.423415707867804*^9}}],
101
102Cell[BoxData[
103 RowBox[{"CheckHermiticity", "[",
104  RowBox[{"LQS", ",",
105   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]], "Input",
106 CellChangeTimes->{{3.414763604111637*^9, 3.4147636044973803`*^9}, {
107  3.554207230374851*^9, 3.554207231973926*^9}}],
108
109Cell["\<\
110FeynRules can check whether all the values of the masses given in the model \
111file are in agreement with the masses given in the Lagrangian:\
112\>", "Text",
113 CellChangeTimes->{{3.423415751305388*^9, 3.423415792594919*^9}, {
114  3.5542072544005003`*^9, 3.554207255988462*^9}}],
115
116Cell[BoxData[
117 RowBox[{"CheckMassSpectrum", "[",
118  RowBox[{"LSM", "+", "LQS"}], "]"}]], "Input",
119 CellChangeTimes->{{3.414736325829801*^9, 3.414736326065997*^9},
120   3.4147528281460342`*^9, {3.5542072604055557`*^9, 3.554207262590489*^9}, {
121   3.554290298409437*^9, 3.554290300070582*^9}}],
122
123Cell["\<\
124In the same way, FeynRules can check whether all kinetic terms are correctly \
125normalized :\
126\>", "Text",
127 CellChangeTimes->{{3.423415798800013*^9, 3.423415817964087*^9}}],
128
129Cell[BoxData[
130 RowBox[{"CheckKineticTermNormalisation", "[",
131  RowBox[{
132   RowBox[{"LSM", "+", "LQS"}], ",",
133   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]], "Input",
134 CellChangeTimes->{{3.554290748499344*^9, 3.554290750913568*^9}}],
135
136Cell[BoxData[
137 FrameBox[
138  RowBox[{
139  "\[SelectionPlaceholder]", " ", "Feynman", " ", "rules", " ", "originating",
140    " ", "from", " ", "L_S0", " ", "and", " ", "L_S1", " ", "using", " ",
141   "the", " ", "unitary", " ", "gauge"}]]], "Input",
142 CellChangeTimes->{{3.554206797674385*^9, 3.554206865976857*^9}, {
143  3.5542069112258472`*^9, 3.554206917083228*^9}}],
144
145Cell[BoxData[{
146 RowBox[{
147  RowBox[{"FeynmanGauge", "=", "False"}], ";"}], "\[IndentingNewLine]",
148 RowBox[{"FeynmanRules", "[",
149  RowBox[{"LS0", " ", ",",
150   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]}], "Input",
151 CellChangeTimes->{{3.5542147888623667`*^9, 3.5542147904036713`*^9},
152   3.5542148289281187`*^9}],
153
154Cell[BoxData[{
155 RowBox[{
156  RowBox[{"FeynmanGauge", "=", "False"}], ";"}], "\[IndentingNewLine]",
157 RowBox[{"FeynmanRules", "[",
158  RowBox[{"LS1", ",",
159   RowBox[{"FlavorExpand", "\[Rule]", "True"}]}], "]"}]}], "Input",
160 CellChangeTimes->{{3.55421481319798*^9, 3.5542148158116817`*^9}}],
161
162Cell[BoxData[
163 FrameBox[
164  RowBox[{
165  "\[SelectionPlaceholder]", " ", "Generating", " ", "the", " ", "UFO", " ",
166   "output", " ", "used", " ", "by", " ", "MadGraph5"}]]], "Input",
167 CellChangeTimes->{{3.554207020764146*^9, 3.554207046179934*^9}, {
168  3.554207083712941*^9, 3.55420709708432*^9}}],
169
170Cell[BoxData[{
171 RowBox[{
172  RowBox[{"FeynmanGauge", "=", "False"}], ";"}], "\[IndentingNewLine]",
173 RowBox[{"WriteUFO", "[",
174  RowBox[{
175  "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LS0", ",",
176   " ", "LS1", ",", " ",
177   RowBox[{"Output", " ", "\[Rule]", " ", "\"\<SM_LS0_LS1_UFO\>\""}]}],
178  "]"}]}], "Input"],
179
180Cell[BoxData[
181 RowBox[{
182  FrameBox[
183   RowBox[{"\[SelectionPlaceholder]", " ", "FeynArts", " ", "output"}]],
184  " "}]], "Input",
185 CellChangeTimes->{{3.554210222303082*^9, 3.554210236874051*^9}}],
186
187Cell["FeynArts also supports the Feynman gauge", "Text",
188 CellChangeTimes->{{3.4121854851672792`*^9, 3.412185494183435*^9}}],
189
190Cell[BoxData[{
191 RowBox[{
192  RowBox[{"FeynmanGauge", " ", "=", " ", "False"}],
193  ";"}], "\[IndentingNewLine]",
194 RowBox[{"WriteFeynArtsOutput", "[",
195  RowBox[{
196  "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LS0", ",",
197   "LS1", ",", " ",
198   RowBox[{"FlavorExpand", "\[Rule]", "SU2W"}], ",",
199   RowBox[{"Output", "\[Rule]", "\"\<SM_LS0_LS1_FA\>\""}]}], "]"}]}], "Input",\
200
201 CellChangeTimes->{{3.411911018354619*^9, 3.4119110196967573`*^9}, {
202   3.412185470680571*^9, 3.412185506240796*^9}, {3.4234161971513453`*^9,
203   3.423416209716672*^9}, {3.5542103799698553`*^9, 3.554210404862816*^9}, {
204   3.554213142736615*^9, 3.554213144826766*^9}, {3.554214531350573*^9,
205   3.554214532452024*^9}, {3.5542147194543047`*^9, 3.554214722231299*^9},
206   3.5542907766028643`*^9}],
207
208Cell[BoxData[
209 FrameBox[
210  RowBox[{
211  "\[SelectionPlaceholder]", " ", "CalcHep", " ", "output"}]]], "Input",
212 CellChangeTimes->{{3.554210492982093*^9, 3.554210505875268*^9}}],
213
214Cell["CalcHep also supports the Feynman gauge", "Text",
215 CellChangeTimes->{{3.4121853722243843`*^9, 3.412185380407339*^9}}],
216
217Cell[BoxData[{
218 RowBox[{
219  RowBox[{"FeynmanGauge", "=", "False"}], ";"}], "\[IndentingNewLine]",
220 RowBox[{
221  RowBox[{"WriteCHOutput", "[",
222   RowBox[{
223   "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", " ", "LS0",
224     ",", " ", "LS1", ",",
225    RowBox[{"Output", "\[Rule]", "\"\<SM_LS0_LS1_CH\>\""}]}], "]"}],
226  ";"}]}], "Input",
227 CellChangeTimes->{{3.4119110943383636`*^9, 3.41191110054212*^9}, {
228   3.412185409690489*^9, 3.4121854275198507`*^9}, {3.412450479942018*^9,
229   3.4124504817138367`*^9}, {3.554210518056747*^9, 3.554210525845241*^9}, {
230   3.554212177238345*^9, 3.554212186985408*^9}, {3.554214738990678*^9,
231   3.5542147440692997`*^9}, 3.5542907733944674`*^9}],
232
233Cell[BoxData[""], "Input",
234 CellChangeTimes->{3.412185367693318*^9, 3.412185407947097*^9}]
235},
236WindowSize->{1024, 656},
237WindowMargins->{{Automatic, -2}, {Automatic, 0}},
238ShowSelection->True,
239Magnification->1.5,
240FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 6, \
2412010)",
242StyleDefinitions->"Default.nb"
243]
244(* End of Notebook Content *)
245
246(* Internal cache information *)
247(*CellTagsOutline
248CellTagsIndex->{}
249*)
250(*CellTagsIndex
251CellTagsIndex->{}
252*)
253(*NotebookFileOutline
254Notebook[{
255Cell[545, 20, 229, 6, 82, "Input"],
256Cell[777, 28, 398, 8, 87, "Input"],
257Cell[1178, 38, 54, 1, 40, "Input"],
258Cell[1235, 41, 319, 7, 42, "Input"],
259Cell[1557, 50, 365, 8, 83, "Input"],
260Cell[1925, 60, 1083, 15, 40, "Input"],
261Cell[3011, 77, 227, 5, 40, "Input"],
262Cell[3241, 84, 236, 6, 83, "Input"],
263Cell[3480, 92, 175, 4, 40, "Input"],
264Cell[3658, 98, 134, 1, 39, "Text"],
265Cell[3795, 101, 255, 5, 40, "Input"],
266Cell[4053, 108, 281, 5, 62, "Text"],
267Cell[4337, 115, 289, 5, 40, "Input"],
268Cell[4629, 122, 181, 4, 39, "Text"],
269Cell[4813, 128, 243, 5, 40, "Input"],
270Cell[5059, 135, 359, 7, 59, "Input"],
271Cell[5421, 144, 322, 7, 64, "Input"],
272Cell[5746, 153, 286, 6, 64, "Input"],
273Cell[6035, 161, 295, 6, 58, "Input"],
274Cell[6333, 169, 333, 8, 65, "Input"],
275Cell[6669, 179, 194, 5, 58, "Input"],
276Cell[6866, 186, 124, 1, 39, "Text"],
277Cell[6993, 189, 788, 16, 88, "Input"],
278Cell[7784, 207, 174, 4, 58, "Input"],
279Cell[7961, 213, 123, 1, 39, "Text"],
280Cell[8087, 216, 692, 14, 88, "Input"],
281Cell[8782, 232, 90, 1, 40, "Input"]
282}
283]
284*)
285
286(* End of internal cache information *)