1 | (* Content-type: application/vnd.wolfram.mathematica *)
|
---|
2 |
|
---|
3 | (*** Wolfram Notebook File ***)
|
---|
4 | (* http://www.wolfram.com/nb *)
|
---|
5 |
|
---|
6 | (* CreatedBy='Mathematica 10.4' *)
|
---|
7 |
|
---|
8 | (*CacheID: 234*)
|
---|
9 | (* Internal cache information:
|
---|
10 | NotebookFileLineBreakTest
|
---|
11 | NotebookFileLineBreakTest
|
---|
12 | NotebookDataPosition[ 158, 7]
|
---|
13 | NotebookDataLength[ 20744, 719]
|
---|
14 | NotebookOptionsPosition[ 18991, 651]
|
---|
15 | NotebookOutlinePosition[ 19348, 667]
|
---|
16 | CellTagsIndexPosition[ 19305, 664]
|
---|
17 | WindowFrame->Normal*)
|
---|
18 |
|
---|
19 | (* Beginning of Notebook Content *)
|
---|
20 | Notebook[{
|
---|
21 |
|
---|
22 | Cell[CellGroupData[{
|
---|
23 | Cell[BoxData[
|
---|
24 | RowBox[{"<<",
|
---|
25 | "\"\</Users/marcozaro/Physics/feyncalc-master/FeynCalc/fc.m\>\""}]], "Input"],
|
---|
26 |
|
---|
27 | Cell[CellGroupData[{
|
---|
28 |
|
---|
29 | Cell[BoxData[
|
---|
30 | FormBox[
|
---|
31 | InterpretationBox[
|
---|
32 | RowBox[{
|
---|
33 | StyleBox["\<\"FeynCalc \"\>", "Text",
|
---|
34 | StripOnInput->False,
|
---|
35 | FontWeight->Bold], "\[InvisibleSpace]",
|
---|
36 | StyleBox["\<\"9.1.0 (development version). For help, use the \"\>", "Text",
|
---|
37 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
38 | StyleBox[
|
---|
39 | TagBox[
|
---|
40 | ButtonBox[
|
---|
41 | RowBox[{"documentation", " ", "center"}],
|
---|
42 | BaseStyle->"Link",
|
---|
43 | ButtonData:>"paclet:FeynCalc/",
|
---|
44 | ButtonNote->"paclet:FeynCalc/"],
|
---|
45 | DisplayForm], "Text",
|
---|
46 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
47 | StyleBox["\<\", check out the \"\>", "Text",
|
---|
48 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
49 | StyleBox[
|
---|
50 | TagBox[
|
---|
51 | ButtonBox["wiki",
|
---|
52 | BaseStyle->"Hyperlink",
|
---|
53 | ButtonData:>{
|
---|
54 | URL["https://github.com/FeynCalc/feyncalc/wiki"], None},
|
---|
55 | ButtonNote->"https://github.com/FeynCalc/feyncalc/wiki"],
|
---|
56 | DisplayForm], "Text",
|
---|
57 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
58 | StyleBox["\<\" or write to the \"\>", "Text",
|
---|
59 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
60 | StyleBox[
|
---|
61 | TagBox[
|
---|
62 | ButtonBox[
|
---|
63 | RowBox[{"mailing", " ",
|
---|
64 | RowBox[{"list", "."}]}],
|
---|
65 | BaseStyle->"Hyperlink",
|
---|
66 | ButtonData:>{
|
---|
67 | URL["http://www.feyncalc.org/forum/"], None},
|
---|
68 | ButtonNote->"http://www.feyncalc.org/forum/"],
|
---|
69 | DisplayForm], "Text",
|
---|
70 | StripOnInput->False]}],
|
---|
71 | SequenceForm[
|
---|
72 | Style["FeynCalc ", "Text", Bold],
|
---|
73 | Style["9.1.0 (development version). For help, use the ", "Text"],
|
---|
74 | Style[
|
---|
75 | DisplayForm[
|
---|
76 | ButtonBox[
|
---|
77 | "documentation center", BaseStyle -> "Link", ButtonData :>
|
---|
78 | "paclet:FeynCalc/", ButtonNote -> "paclet:FeynCalc/"]], "Text"],
|
---|
79 | Style[", check out the ", "Text"],
|
---|
80 | Style[
|
---|
81 | DisplayForm[
|
---|
82 | ButtonBox["wiki", ButtonData :> {
|
---|
83 | URL["https://github.com/FeynCalc/feyncalc/wiki"], None}, BaseStyle ->
|
---|
84 | "Hyperlink", ButtonNote ->
|
---|
85 | "https://github.com/FeynCalc/feyncalc/wiki"]], "Text"],
|
---|
86 | Style[" or write to the ", "Text"],
|
---|
87 | Style[
|
---|
88 | DisplayForm[
|
---|
89 | ButtonBox["mailing list.", ButtonData :> {
|
---|
90 | URL["http://www.feyncalc.org/forum/"], None}, BaseStyle ->
|
---|
91 | "Hyperlink", ButtonNote -> "http://www.feyncalc.org/forum/"]], "Text"]],
|
---|
92 | Editable->False], TraditionalForm]], "Print",
|
---|
93 | CellChangeTimes->{3.69303902949454*^9, 3.6930399100242023`*^9}],
|
---|
94 |
|
---|
95 | Cell[BoxData[
|
---|
96 | FormBox[
|
---|
97 | InterpretationBox[
|
---|
98 | RowBox[{
|
---|
99 | StyleBox["\<\"See also the supplied \"\>", "Text",
|
---|
100 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
101 | StyleBox[
|
---|
102 | TagBox[
|
---|
103 | ButtonBox[
|
---|
104 | RowBox[{"examples", "."}],
|
---|
105 | BaseStyle->"Hyperlink",
|
---|
106 | ButtonFunction:>SystemOpen[
|
---|
107 | FileNameJoin[{FeynCalc`$FeynCalcDirectory, "Examples"}]],
|
---|
108 | Evaluator->Automatic,
|
---|
109 | Method->"Preemptive"],
|
---|
110 | DisplayForm], "Text",
|
---|
111 | StripOnInput->False], "\[InvisibleSpace]",
|
---|
112 | StyleBox["\<\" If you use FeynCalc in your research, please cite\"\>",
|
---|
113 | "Text",
|
---|
114 | StripOnInput->False]}],
|
---|
115 | SequenceForm[
|
---|
116 | Style["See also the supplied ", "Text"],
|
---|
117 | Style[
|
---|
118 | DisplayForm[
|
---|
119 | ButtonBox[
|
---|
120 | "examples.", BaseStyle -> "Hyperlink", ButtonFunction :> SystemOpen[
|
---|
121 | FileNameJoin[{FeynCalc`$FeynCalcDirectory, "Examples"}]], Evaluator ->
|
---|
122 | Automatic, Method -> "Preemptive"]], "Text"],
|
---|
123 | Style[" If you use FeynCalc in your research, please cite", "Text"]],
|
---|
124 | Editable->False], TraditionalForm]], "Print",
|
---|
125 | CellChangeTimes->{3.69303902949454*^9, 3.693039910029832*^9}],
|
---|
126 |
|
---|
127 | Cell[BoxData[
|
---|
128 | FormBox[
|
---|
129 | StyleBox["\<\" \[Bullet] V. Shtabovenko, R. Mertig and F. Orellana, TUM-EFT \
|
---|
130 | 71/15, arXiv:1601.01167\"\>", "Text",
|
---|
131 | StripOnInput->False], TraditionalForm]], "Print",
|
---|
132 | CellChangeTimes->{3.69303902949454*^9, 3.693039910034374*^9}],
|
---|
133 |
|
---|
134 | Cell[BoxData[
|
---|
135 | FormBox[
|
---|
136 | StyleBox["\<\" \[Bullet] R. Mertig, M. B\[ODoubleDot]hm, and A. Denner, \
|
---|
137 | Comput. Phys. Commun., 64, 345-359, 1991.\"\>", "Text",
|
---|
138 | StripOnInput->False], TraditionalForm]], "Print",
|
---|
139 | CellChangeTimes->{3.69303902949454*^9, 3.69303991003861*^9}]
|
---|
140 | }, Open ]]
|
---|
141 | }, Open ]],
|
---|
142 |
|
---|
143 | Cell[CellGroupData[{
|
---|
144 |
|
---|
145 | Cell["Spin-1/2 vs scalar quarks", "Title",
|
---|
146 | CellChangeTimes->{{3.693039048153269*^9, 3.693039055120265*^9}}],
|
---|
147 |
|
---|
148 | Cell[CellGroupData[{
|
---|
149 |
|
---|
150 | Cell["we compute the cross section e+(p1) e-(p2) > q(p3) q~(p4)", "Section",
|
---|
151 | CellChangeTimes->{{3.693039062567734*^9, 3.693039094965802*^9}}],
|
---|
152 |
|
---|
153 | Cell[BoxData[{
|
---|
154 | RowBox[{
|
---|
155 | RowBox[{
|
---|
156 | RowBox[{"ScalarProduct", "[",
|
---|
157 | RowBox[{"p1", ",", "p1"}], "]"}], " ", "=", " ", "0"}],
|
---|
158 | ";"}], "\[IndentingNewLine]",
|
---|
159 | RowBox[{
|
---|
160 | RowBox[{
|
---|
161 | RowBox[{"ScalarProduct", "[",
|
---|
162 | RowBox[{"p2", ",", "p2"}], "]"}], " ", "=", " ", "0"}],
|
---|
163 | ";"}], "\[IndentingNewLine]",
|
---|
164 | RowBox[{
|
---|
165 | RowBox[{
|
---|
166 | RowBox[{"ScalarProduct", "[",
|
---|
167 | RowBox[{"p3", ",", "p3"}], "]"}], " ", "=", " ", "0"}],
|
---|
168 | ";"}], "\[IndentingNewLine]",
|
---|
169 | RowBox[{
|
---|
170 | RowBox[{
|
---|
171 | RowBox[{"ScalarProduct", "[",
|
---|
172 | RowBox[{"p4", ",", "p4"}], "]"}], " ", "=", " ", "0"}],
|
---|
173 | ";"}], "\[IndentingNewLine]",
|
---|
174 | RowBox[{
|
---|
175 | RowBox[{
|
---|
176 | RowBox[{"ScalarProduct", "[",
|
---|
177 | RowBox[{"p1", ",", "p2"}], "]"}], " ", "=", " ",
|
---|
178 | RowBox[{"s", "/", "2"}]}], ";"}], "\[IndentingNewLine]",
|
---|
179 | RowBox[{
|
---|
180 | RowBox[{
|
---|
181 | RowBox[{"ScalarProduct", "[",
|
---|
182 | RowBox[{"p3", ",", "p4"}], "]"}], " ", "=", " ",
|
---|
183 | RowBox[{"s", "/", "2"}]}], ";"}], "\[IndentingNewLine]",
|
---|
184 | RowBox[{
|
---|
185 | RowBox[{
|
---|
186 | RowBox[{"ScalarProduct", "[",
|
---|
187 | RowBox[{"p1", ",", "p3"}], "]"}], " ", "=", " ",
|
---|
188 | RowBox[{
|
---|
189 | RowBox[{"-", "t"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]",
|
---|
190 | RowBox[{
|
---|
191 | RowBox[{
|
---|
192 | RowBox[{"ScalarProduct", "[",
|
---|
193 | RowBox[{"p2", ",", "p4"}], "]"}], " ", "=", " ",
|
---|
194 | RowBox[{
|
---|
195 | RowBox[{"-", "t"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]",
|
---|
196 | RowBox[{
|
---|
197 | RowBox[{
|
---|
198 | RowBox[{"ScalarProduct", "[",
|
---|
199 | RowBox[{"p1", ",", "p4"}], "]"}], " ", "=", " ",
|
---|
200 | RowBox[{
|
---|
201 | RowBox[{"-", "u"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]",
|
---|
202 | RowBox[{
|
---|
203 | RowBox[{
|
---|
204 | RowBox[{"ScalarProduct", "[",
|
---|
205 | RowBox[{"p2", ",", "p3"}], "]"}], " ", "=", " ",
|
---|
206 | RowBox[{
|
---|
207 | RowBox[{"-", "u"}], "/", "2"}]}], ";"}]}], "Input",
|
---|
208 | CellChangeTimes->{{3.6930391053897133`*^9, 3.6930391058309107`*^9}}]
|
---|
209 | }, Open ]],
|
---|
210 |
|
---|
211 | Cell[CellGroupData[{
|
---|
212 |
|
---|
213 | Cell["The leptonic part of the cross section", "Section",
|
---|
214 | CellChangeTimes->{{3.693039062567734*^9, 3.693039121124161*^9}}],
|
---|
215 |
|
---|
216 | Cell[CellGroupData[{
|
---|
217 |
|
---|
218 | Cell[BoxData[
|
---|
219 | RowBox[{"Lmunu", " ", "=", " ",
|
---|
220 | RowBox[{"Tr", "[",
|
---|
221 | RowBox[{
|
---|
222 | RowBox[{"GS", "[", "p1", "]"}], ".",
|
---|
223 | RowBox[{"GA", "[", "mu", "]"}], ".",
|
---|
224 | RowBox[{"GS", "[", "p2", "]"}], ".",
|
---|
225 | RowBox[{"GA", "[", "nu", "]"}]}], "]"}]}]], "Input",
|
---|
226 | CellChangeTimes->{{3.693039125452209*^9, 3.6930391752723217`*^9}}],
|
---|
227 |
|
---|
228 | Cell[BoxData[
|
---|
229 | FormBox[
|
---|
230 | RowBox[{
|
---|
231 | RowBox[{"-", "2"}], " ",
|
---|
232 | RowBox[{"(",
|
---|
233 | RowBox[{
|
---|
234 | RowBox[{"s", " ",
|
---|
235 | FormBox[
|
---|
236 | SuperscriptBox[
|
---|
237 | OverscriptBox["g", "_"],
|
---|
238 | RowBox[{
|
---|
239 | FormBox[
|
---|
240 | FormBox[
|
---|
241 | FormBox["mu",
|
---|
242 | TraditionalForm],
|
---|
243 | TraditionalForm],
|
---|
244 | TraditionalForm],
|
---|
245 | FormBox[
|
---|
246 | FormBox[
|
---|
247 | FormBox["nu",
|
---|
248 | TraditionalForm],
|
---|
249 | TraditionalForm],
|
---|
250 | TraditionalForm]}]],
|
---|
251 | TraditionalForm]}], "-",
|
---|
252 | RowBox[{"2", " ",
|
---|
253 | FormBox[
|
---|
254 | SuperscriptBox[
|
---|
255 | FormBox[
|
---|
256 | OverscriptBox[
|
---|
257 | FormBox["p2",
|
---|
258 | TraditionalForm], "_"],
|
---|
259 | TraditionalForm],
|
---|
260 | FormBox[
|
---|
261 | FormBox[
|
---|
262 | FormBox["mu",
|
---|
263 | TraditionalForm],
|
---|
264 | TraditionalForm],
|
---|
265 | TraditionalForm]],
|
---|
266 | TraditionalForm], " ",
|
---|
267 | FormBox[
|
---|
268 | SuperscriptBox[
|
---|
269 | FormBox[
|
---|
270 | OverscriptBox[
|
---|
271 | FormBox["p1",
|
---|
272 | TraditionalForm], "_"],
|
---|
273 | TraditionalForm],
|
---|
274 | FormBox[
|
---|
275 | FormBox[
|
---|
276 | FormBox["nu",
|
---|
277 | TraditionalForm],
|
---|
278 | TraditionalForm],
|
---|
279 | TraditionalForm]],
|
---|
280 | TraditionalForm]}], "-",
|
---|
281 | RowBox[{"2", " ",
|
---|
282 | FormBox[
|
---|
283 | SuperscriptBox[
|
---|
284 | FormBox[
|
---|
285 | OverscriptBox[
|
---|
286 | FormBox["p1",
|
---|
287 | TraditionalForm], "_"],
|
---|
288 | TraditionalForm],
|
---|
289 | FormBox[
|
---|
290 | FormBox[
|
---|
291 | FormBox["mu",
|
---|
292 | TraditionalForm],
|
---|
293 | TraditionalForm],
|
---|
294 | TraditionalForm]],
|
---|
295 | TraditionalForm], " ",
|
---|
296 | FormBox[
|
---|
297 | SuperscriptBox[
|
---|
298 | FormBox[
|
---|
299 | OverscriptBox[
|
---|
300 | FormBox["p2",
|
---|
301 | TraditionalForm], "_"],
|
---|
302 | TraditionalForm],
|
---|
303 | FormBox[
|
---|
304 | FormBox[
|
---|
305 | FormBox["nu",
|
---|
306 | TraditionalForm],
|
---|
307 | TraditionalForm],
|
---|
308 | TraditionalForm]],
|
---|
309 | TraditionalForm]}]}], ")"}]}], TraditionalForm]], "Output",
|
---|
310 | CellChangeTimes->{{3.6930391756863317`*^9, 3.693039188289091*^9},
|
---|
311 | 3.693039910482806*^9}]
|
---|
312 | }, Open ]]
|
---|
313 | }, Open ]],
|
---|
314 |
|
---|
315 | Cell[CellGroupData[{
|
---|
316 |
|
---|
317 | Cell["\<\
|
---|
318 | The hadronic part of the cross section for spin-1/2, and the total matrix \
|
---|
319 | element\
|
---|
320 | \>", "Section",
|
---|
321 | CellChangeTimes->{{3.693039062567734*^9, 3.693039121124161*^9}, {
|
---|
322 | 3.693039203772478*^9, 3.693039252982416*^9}}],
|
---|
323 |
|
---|
324 | Cell[CellGroupData[{
|
---|
325 |
|
---|
326 | Cell[BoxData[
|
---|
327 | RowBox[{"WmunuFermion", " ", "=", " ",
|
---|
328 | RowBox[{"Tr", "[",
|
---|
329 | RowBox[{
|
---|
330 | RowBox[{"GS", "[", "p3", "]"}], ".",
|
---|
331 | RowBox[{"GA", "[", "mu", "]"}], ".",
|
---|
332 | RowBox[{"GS", "[", "p4", "]"}], ".",
|
---|
333 | RowBox[{"GA", "[", "nu", "]"}]}], "]"}]}]], "Input",
|
---|
334 | CellChangeTimes->{{3.693039125452209*^9, 3.6930391752723217`*^9}, {
|
---|
335 | 3.6930392149801083`*^9, 3.693039230919681*^9}, {3.693127070874337*^9,
|
---|
336 | 3.6931271027239437`*^9}}],
|
---|
337 |
|
---|
338 | Cell[BoxData[
|
---|
339 | FormBox[
|
---|
340 | RowBox[{
|
---|
341 | RowBox[{"-", "2"}], " ",
|
---|
342 | RowBox[{"(",
|
---|
343 | RowBox[{
|
---|
344 | RowBox[{"s", " ",
|
---|
345 | FormBox[
|
---|
346 | SuperscriptBox[
|
---|
347 | OverscriptBox["g", "_"],
|
---|
348 | RowBox[{
|
---|
349 | FormBox[
|
---|
350 | FormBox[
|
---|
351 | FormBox["mu",
|
---|
352 | TraditionalForm],
|
---|
353 | TraditionalForm],
|
---|
354 | TraditionalForm],
|
---|
355 | FormBox[
|
---|
356 | FormBox[
|
---|
357 | FormBox["nu",
|
---|
358 | TraditionalForm],
|
---|
359 | TraditionalForm],
|
---|
360 | TraditionalForm]}]],
|
---|
361 | TraditionalForm]}], "-",
|
---|
362 | RowBox[{"2", " ",
|
---|
363 | FormBox[
|
---|
364 | SuperscriptBox[
|
---|
365 | FormBox[
|
---|
366 | OverscriptBox[
|
---|
367 | FormBox["p4",
|
---|
368 | TraditionalForm], "_"],
|
---|
369 | TraditionalForm],
|
---|
370 | FormBox[
|
---|
371 | FormBox[
|
---|
372 | FormBox["mu",
|
---|
373 | TraditionalForm],
|
---|
374 | TraditionalForm],
|
---|
375 | TraditionalForm]],
|
---|
376 | TraditionalForm], " ",
|
---|
377 | FormBox[
|
---|
378 | SuperscriptBox[
|
---|
379 | FormBox[
|
---|
380 | OverscriptBox[
|
---|
381 | FormBox["p3",
|
---|
382 | TraditionalForm], "_"],
|
---|
383 | TraditionalForm],
|
---|
384 | FormBox[
|
---|
385 | FormBox[
|
---|
386 | FormBox["nu",
|
---|
387 | TraditionalForm],
|
---|
388 | TraditionalForm],
|
---|
389 | TraditionalForm]],
|
---|
390 | TraditionalForm]}], "-",
|
---|
391 | RowBox[{"2", " ",
|
---|
392 | FormBox[
|
---|
393 | SuperscriptBox[
|
---|
394 | FormBox[
|
---|
395 | OverscriptBox[
|
---|
396 | FormBox["p3",
|
---|
397 | TraditionalForm], "_"],
|
---|
398 | TraditionalForm],
|
---|
399 | FormBox[
|
---|
400 | FormBox[
|
---|
401 | FormBox["mu",
|
---|
402 | TraditionalForm],
|
---|
403 | TraditionalForm],
|
---|
404 | TraditionalForm]],
|
---|
405 | TraditionalForm], " ",
|
---|
406 | FormBox[
|
---|
407 | SuperscriptBox[
|
---|
408 | FormBox[
|
---|
409 | OverscriptBox[
|
---|
410 | FormBox["p4",
|
---|
411 | TraditionalForm], "_"],
|
---|
412 | TraditionalForm],
|
---|
413 | FormBox[
|
---|
414 | FormBox[
|
---|
415 | FormBox["nu",
|
---|
416 | TraditionalForm],
|
---|
417 | TraditionalForm],
|
---|
418 | TraditionalForm]],
|
---|
419 | TraditionalForm]}]}], ")"}]}], TraditionalForm]], "Output",
|
---|
420 | CellChangeTimes->{3.693039242861771*^9, 3.693039910541121*^9}]
|
---|
421 | }, Open ]],
|
---|
422 |
|
---|
423 | Cell[CellGroupData[{
|
---|
424 |
|
---|
425 | Cell[BoxData[
|
---|
426 | RowBox[{"AmpFermion", " ", "=", " ",
|
---|
427 | RowBox[{
|
---|
428 | RowBox[{
|
---|
429 | RowBox[{"Contract", " ", "[",
|
---|
430 | RowBox[{"WmunuFermion", " ", "Lmunu"}], "]"}], "/", " ",
|
---|
431 | RowBox[{
|
---|
432 | RowBox[{"SP", "[",
|
---|
433 | RowBox[{"p1", "+", "p2"}], "]"}], "^", "2"}]}], " ", "//",
|
---|
434 | "ExpandScalarProduct"}]}]], "Input",
|
---|
435 | CellChangeTimes->{{3.693039245802354*^9, 3.693039312937352*^9}}],
|
---|
436 |
|
---|
437 | Cell[BoxData[
|
---|
438 | FormBox[
|
---|
439 | FractionBox[
|
---|
440 | RowBox[{
|
---|
441 | RowBox[{"8", " ",
|
---|
442 | SuperscriptBox["t", "2"]}], "+",
|
---|
443 | RowBox[{"8", " ",
|
---|
444 | SuperscriptBox["u", "2"]}]}],
|
---|
445 | SuperscriptBox["s", "2"]], TraditionalForm]], "Output",
|
---|
446 | CellChangeTimes->{
|
---|
447 | 3.693039276340749*^9, {3.693039308093046*^9, 3.6930393131384697`*^9},
|
---|
448 | 3.693039910587453*^9}]
|
---|
449 | }, Open ]]
|
---|
450 | }, Open ]],
|
---|
451 |
|
---|
452 | Cell[CellGroupData[{
|
---|
453 |
|
---|
454 | Cell["\<\
|
---|
455 | The hadronic part of the cross section for a scalar, and the total matrix \
|
---|
456 | element. The S(k)S(k\[CloseCurlyQuote])V(-k-k\[CloseCurlyQuote], mu) vertex \
|
---|
457 | is ie(k-k\[CloseCurlyQuote])^mu\
|
---|
458 | \>", "Section",
|
---|
459 | CellChangeTimes->{{3.693039062567734*^9, 3.693039121124161*^9}, {
|
---|
460 | 3.693039203772478*^9, 3.693039252982416*^9}, {3.693039320241497*^9,
|
---|
461 | 3.693039408731784*^9}}],
|
---|
462 |
|
---|
463 | Cell[CellGroupData[{
|
---|
464 |
|
---|
465 | Cell[BoxData[
|
---|
466 | RowBox[{"WmunuScalar", "=",
|
---|
467 | RowBox[{
|
---|
468 | RowBox[{"FV", "[",
|
---|
469 | RowBox[{
|
---|
470 | RowBox[{"p3", "-", "p4"}], ",", "mu"}], "]"}],
|
---|
471 | RowBox[{"FV", "[",
|
---|
472 | RowBox[{
|
---|
473 | RowBox[{"p3", "-", "p4"}], ",", "nu"}], "]"}]}]}]], "Input",
|
---|
474 | CellChangeTimes->{{3.693039323525537*^9, 3.6930393334056673`*^9}, {
|
---|
475 | 3.693039420533457*^9, 3.6930394353167067`*^9}}],
|
---|
476 |
|
---|
477 | Cell[BoxData[
|
---|
478 | FormBox[
|
---|
479 | RowBox[{
|
---|
480 | FormBox[
|
---|
481 | FormBox[
|
---|
482 | SuperscriptBox[
|
---|
483 | RowBox[{"(",
|
---|
484 | FormBox[
|
---|
485 | FormBox[
|
---|
486 | RowBox[{
|
---|
487 | OverscriptBox[
|
---|
488 | FormBox["p3",
|
---|
489 | TraditionalForm], "_"], "-",
|
---|
490 | OverscriptBox[
|
---|
491 | FormBox["p4",
|
---|
492 | TraditionalForm], "_"]}],
|
---|
493 | TraditionalForm],
|
---|
494 | TraditionalForm], ")"}],
|
---|
495 | FormBox[
|
---|
496 | FormBox[
|
---|
497 | FormBox["mu",
|
---|
498 | TraditionalForm],
|
---|
499 | TraditionalForm],
|
---|
500 | TraditionalForm]],
|
---|
501 | TraditionalForm],
|
---|
502 | TraditionalForm], " ",
|
---|
503 | FormBox[
|
---|
504 | FormBox[
|
---|
505 | SuperscriptBox[
|
---|
506 | RowBox[{"(",
|
---|
507 | FormBox[
|
---|
508 | FormBox[
|
---|
509 | RowBox[{
|
---|
510 | OverscriptBox[
|
---|
511 | FormBox["p3",
|
---|
512 | TraditionalForm], "_"], "-",
|
---|
513 | OverscriptBox[
|
---|
514 | FormBox["p4",
|
---|
515 | TraditionalForm], "_"]}],
|
---|
516 | TraditionalForm],
|
---|
517 | TraditionalForm], ")"}],
|
---|
518 | FormBox[
|
---|
519 | FormBox[
|
---|
520 | FormBox["nu",
|
---|
521 | TraditionalForm],
|
---|
522 | TraditionalForm],
|
---|
523 | TraditionalForm]],
|
---|
524 | TraditionalForm],
|
---|
525 | TraditionalForm]}], TraditionalForm]], "Output",
|
---|
526 | CellChangeTimes->{3.693039458151698*^9, 3.693039910631071*^9}]
|
---|
527 | }, Open ]],
|
---|
528 |
|
---|
529 | Cell[CellGroupData[{
|
---|
530 |
|
---|
531 | Cell[BoxData[
|
---|
532 | RowBox[{"AmpScalar", " ", "=", " ",
|
---|
533 | RowBox[{
|
---|
534 | RowBox[{
|
---|
535 | RowBox[{
|
---|
536 | RowBox[{"Contract", " ", "[",
|
---|
537 | RowBox[{"WmunuScalar", " ", "Lmunu"}], "]"}], "/", " ",
|
---|
538 | RowBox[{
|
---|
539 | RowBox[{"SP", "[",
|
---|
540 | RowBox[{"p1", "+", "p2"}], "]"}], "^", "2"}]}], " ", "//",
|
---|
541 | "ExpandScalarProduct"}], "//", "Simplify"}]}]], "Input",
|
---|
542 | CellChangeTimes->{{3.693039447693351*^9, 3.693039466254558*^9}}],
|
---|
543 |
|
---|
544 | Cell[BoxData[
|
---|
545 | FormBox[
|
---|
546 | FractionBox[
|
---|
547 | RowBox[{"2", " ",
|
---|
548 | RowBox[{"(",
|
---|
549 | RowBox[{
|
---|
550 | SuperscriptBox["s", "2"], "-",
|
---|
551 | SuperscriptBox[
|
---|
552 | RowBox[{"(",
|
---|
553 | RowBox[{"t", "-", "u"}], ")"}], "2"]}], ")"}]}],
|
---|
554 | SuperscriptBox["s", "2"]], TraditionalForm]], "Output",
|
---|
555 | CellChangeTimes->{{3.693039454137134*^9, 3.693039466695468*^9},
|
---|
556 | 3.693039910690022*^9}]
|
---|
557 | }, Open ]]
|
---|
558 | }, Open ]],
|
---|
559 |
|
---|
560 | Cell[CellGroupData[{
|
---|
561 |
|
---|
562 | Cell["In DIS, we need to replace s->T, t->S, u->U", "Section",
|
---|
563 | CellChangeTimes->{{3.693039062567734*^9, 3.693039121124161*^9}, {
|
---|
564 | 3.693039203772478*^9, 3.693039252982416*^9}, {3.693039320241497*^9,
|
---|
565 | 3.693039408731784*^9}, {3.693039486354397*^9, 3.6930395019281816`*^9}}],
|
---|
566 |
|
---|
567 | Cell[BoxData[{
|
---|
568 | RowBox[{
|
---|
569 | RowBox[{"DISreplace", " ", "=", " ",
|
---|
570 | RowBox[{"{",
|
---|
571 | RowBox[{
|
---|
572 | RowBox[{"s", "\[Rule]", "T"}], ",", " ",
|
---|
573 | RowBox[{"t", "\[Rule]", "S"}], ",", " ",
|
---|
574 | RowBox[{"u", "\[Rule]", "U"}]}], "}"}]}], ";"}], "\[IndentingNewLine]",
|
---|
575 | RowBox[{
|
---|
576 | RowBox[{"DISreplace2", " ", "=", " ",
|
---|
577 | RowBox[{"{",
|
---|
578 | RowBox[{
|
---|
579 | RowBox[{"T", "\[Rule]",
|
---|
580 | RowBox[{"-", "Q2"}]}], ",", " ",
|
---|
581 | RowBox[{"U", "\[Rule]", " ",
|
---|
582 | RowBox[{"S",
|
---|
583 | RowBox[{"(",
|
---|
584 | RowBox[{"y", "-", "1"}], ")"}]}]}]}], "}"}]}], ";"}]}], "Input",
|
---|
585 | CellChangeTimes->{{3.693039511539997*^9, 3.693039521223592*^9}, {
|
---|
586 | 3.693039662301601*^9, 3.69303966997472*^9}, {3.6930398660399427`*^9,
|
---|
587 | 3.6930399277280912`*^9}}],
|
---|
588 |
|
---|
589 | Cell[BoxData[{
|
---|
590 | RowBox[{"AmpFermionDIS", " ", "=", " ",
|
---|
591 | RowBox[{
|
---|
592 | RowBox[{"AmpFermion", " ", "/.", "DISreplace"}], " ", "/.",
|
---|
593 | "DISreplace2"}]}], "\[IndentingNewLine]",
|
---|
594 | RowBox[{"AmpScalarDIS", " ", "=", " ",
|
---|
595 | RowBox[{
|
---|
596 | RowBox[{
|
---|
597 | RowBox[{
|
---|
598 | RowBox[{
|
---|
599 | RowBox[{"AmpScalar", " ", "/.", "DISreplace"}], "/.", "DISreplace2"}],
|
---|
600 | " ", "/.",
|
---|
601 | RowBox[{"{",
|
---|
602 | RowBox[{"Q2", "\[Rule]", " ",
|
---|
603 | RowBox[{
|
---|
604 | RowBox[{"-", "S"}], " ", "x", " ", "y"}]}], "}"}]}], " ", "/.", " ",
|
---|
605 |
|
---|
606 | RowBox[{"x", "\[Rule]", "1"}]}], "//", "Simplify", " "}]}]}], "Input",
|
---|
607 | CellChangeTimes->{{3.693039562590953*^9, 3.693039609767437*^9}, {
|
---|
608 | 3.69303988989855*^9, 3.6930399001178102`*^9}, {3.693039933415969*^9,
|
---|
609 | 3.693039948918563*^9}, {3.6930412006969213`*^9, 3.6930412010710793`*^9}, {
|
---|
610 | 3.6930419330796423`*^9, 3.693041934379951*^9}, {3.693041982923629*^9,
|
---|
611 | 3.693041984041925*^9}, {3.6930423392707987`*^9, 3.693042370214587*^9}, {
|
---|
612 | 3.693042462398498*^9, 3.6930425194736233`*^9}}],
|
---|
613 |
|
---|
614 | Cell[CellGroupData[{
|
---|
615 |
|
---|
616 | Cell[BoxData[
|
---|
617 | RowBox[{"(*",
|
---|
618 | RowBox[{
|
---|
619 | RowBox[{
|
---|
620 | RowBox[{"at", " ", "LO"}], " ", "+", " ",
|
---|
621 | RowBox[{
|
---|
622 | "with", " ", "our", " ", "choice", " ", "of", " ", "momenta", " ", "one",
|
---|
623 | " ", "has", " ", "x"}]}], "=", "1"}]}]], "Input",
|
---|
624 | CellChangeTimes->{{3.693042525358066*^9, 3.693042539751954*^9}}],
|
---|
625 |
|
---|
626 | Cell[BoxData[
|
---|
627 | FormBox[
|
---|
628 | FractionBox[
|
---|
629 | RowBox[{
|
---|
630 | RowBox[{"8", " ",
|
---|
631 | SuperscriptBox["S", "2"], " ",
|
---|
632 | SuperscriptBox[
|
---|
633 | RowBox[{"(",
|
---|
634 | RowBox[{"y", "-", "1"}], ")"}], "2"]}], "+",
|
---|
635 | RowBox[{"8", " ",
|
---|
636 | SuperscriptBox["S", "2"]}]}],
|
---|
637 | SuperscriptBox["Q2", "2"]], TraditionalForm]], "Output",
|
---|
638 | CellChangeTimes->{{3.693042470856336*^9, 3.6930424885340147`*^9}}],
|
---|
639 |
|
---|
640 | Cell[BoxData[
|
---|
641 | FormBox[
|
---|
642 | FractionBox[
|
---|
643 | RowBox[{"8", " ",
|
---|
644 | RowBox[{"(",
|
---|
645 | RowBox[{"y", "-", "1"}], ")"}]}],
|
---|
646 | SuperscriptBox["y", "2"]], TraditionalForm]], "Output",
|
---|
647 | CellChangeTimes->{{3.693042470856336*^9, 3.693042488546664*^9}}]
|
---|
648 | }, Open ]]
|
---|
649 | }, Open ]]
|
---|
650 | }, Open ]]
|
---|
651 | },
|
---|
652 | WindowSize->{808, 596},
|
---|
653 | WindowMargins->{{12, Automatic}, {Automatic, 0}},
|
---|
654 | FrontEndVersion->"10.4 for Mac OS X x86 (32-bit, 64-bit Kernel) (February 25, \
|
---|
655 | 2016)",
|
---|
656 | StyleDefinitions->"Default.nb"
|
---|
657 | ]
|
---|
658 | (* End of Notebook Content *)
|
---|
659 |
|
---|
660 | (* Internal cache information *)
|
---|
661 | (*CellTagsOutline
|
---|
662 | CellTagsIndex->{}
|
---|
663 | *)
|
---|
664 | (*CellTagsIndex
|
---|
665 | CellTagsIndex->{}
|
---|
666 | *)
|
---|
667 | (*NotebookFileOutline
|
---|
668 | Notebook[{
|
---|
669 | Cell[CellGroupData[{
|
---|
670 | Cell[580, 22, 109, 2, 28, "Input"],
|
---|
671 | Cell[CellGroupData[{
|
---|
672 | Cell[714, 28, 2391, 64, 44, "Print"],
|
---|
673 | Cell[3108, 94, 1130, 30, 25, "Print"],
|
---|
674 | Cell[4241, 126, 257, 5, 26, "Print"],
|
---|
675 | Cell[4501, 133, 270, 5, 26, "Print"]
|
---|
676 | }, Open ]]
|
---|
677 | }, Open ]],
|
---|
678 | Cell[CellGroupData[{
|
---|
679 | Cell[4820, 144, 108, 1, 92, "Title"],
|
---|
680 | Cell[CellGroupData[{
|
---|
681 | Cell[4953, 149, 142, 1, 64, "Section"],
|
---|
682 | Cell[5098, 152, 1767, 55, 182, "Input"]
|
---|
683 | }, Open ]],
|
---|
684 | Cell[CellGroupData[{
|
---|
685 | Cell[6902, 212, 123, 1, 64, "Section"],
|
---|
686 | Cell[CellGroupData[{
|
---|
687 | Cell[7050, 217, 332, 8, 28, "Input"],
|
---|
688 | Cell[7385, 227, 2085, 83, 34, "Output"]
|
---|
689 | }, Open ]]
|
---|
690 | }, Open ]],
|
---|
691 | Cell[CellGroupData[{
|
---|
692 | Cell[9519, 316, 225, 5, 99, "Section"],
|
---|
693 | Cell[CellGroupData[{
|
---|
694 | Cell[9769, 325, 441, 10, 28, "Input"],
|
---|
695 | Cell[10213, 337, 2055, 82, 34, "Output"]
|
---|
696 | }, Open ]],
|
---|
697 | Cell[CellGroupData[{
|
---|
698 | Cell[12305, 424, 383, 10, 28, "Input"],
|
---|
699 | Cell[12691, 436, 349, 11, 53, "Output"]
|
---|
700 | }, Open ]]
|
---|
701 | }, Open ]],
|
---|
702 | Cell[CellGroupData[{
|
---|
703 | Cell[13089, 453, 376, 7, 99, "Section"],
|
---|
704 | Cell[CellGroupData[{
|
---|
705 | Cell[13490, 464, 364, 10, 28, "Input"],
|
---|
706 | Cell[13857, 476, 1205, 49, 34, "Output"]
|
---|
707 | }, Open ]],
|
---|
708 | Cell[CellGroupData[{
|
---|
709 | Cell[15099, 530, 420, 11, 28, "Input"],
|
---|
710 | Cell[15522, 543, 382, 12, 54, "Output"]
|
---|
711 | }, Open ]]
|
---|
712 | }, Open ]],
|
---|
713 | Cell[CellGroupData[{
|
---|
714 | Cell[15953, 561, 274, 3, 64, "Section"],
|
---|
715 | Cell[16230, 566, 733, 20, 46, "Input"],
|
---|
716 | Cell[16966, 588, 1008, 23, 46, "Input"],
|
---|
717 | Cell[CellGroupData[{
|
---|
718 | Cell[17999, 615, 312, 8, 28, "Input"],
|
---|
719 | Cell[18314, 625, 391, 12, 55, "Output"],
|
---|
720 | Cell[18708, 639, 243, 7, 52, "Output"]
|
---|
721 | }, Open ]]
|
---|
722 | }, Open ]]
|
---|
723 | }, Open ]]
|
---|
724 | }
|
---|
725 | ]
|
---|
726 | *)
|
---|
727 |
|
---|