WAsymm: DYNarrowWidth.nb

File DYNarrowWidth.nb, 14.6 KB (added by (none), 13 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[ 14721, 497]
14NotebookOptionsPosition[ 12913, 429]
15NotebookOutlinePosition[ 13290, 446]
16CellTagsIndexPosition[ 13247, 443]
17WindowFrame->Normal
18ContainsDynamic->False*)
19
20(* Beginning of Notebook Content *)
21Notebook[{
22
23Cell[CellGroupData[{
24Cell["\<\
25The Drell-Yan cross-section in the narrow width approximation\
26\>", "Title",
27 CellChangeTimes->{{3.42892672627883*^9, 3.428926824323577*^9}, {
28 3.428937703319875*^9, 3.428937712432496*^9}, {3.428941361957245*^9,
29 3.4289413752252607`*^9}}],
30
31Cell[CellGroupData[{
32
33Cell["Input FeynCalc", "Subsection"],
34
35Cell[BoxData[
36 RowBox[{
37 RowBox[{"<<", "HighEnergyPhysics`fc`"}], ";"}]], "Input",
38 CellChangeTimes->{{3.428941446581498*^9, 3.428941449045033*^9}}]
39}, Open ]],
40
41Cell[CellGroupData[{
42
43Cell["The partonic cross - section", "Section",
44 CellChangeTimes->{{3.428941607370674*^9, 3.428941611186956*^9}}],
45
46Cell[BoxData[
47 RowBox[{
48 RowBox[{"Near", " ", "the", " ", "W", " ", "resonance"}], ",", " ",
49 RowBox[{
50 RowBox[{
51 "the", " ", "partonic", " ", "cross", " ", "section", " ", "\[Sigma]p",
52 " ", "for", " ", "u", " ", "dbar"}], " ", "\[RightArrow]", " ",
53 SuperscriptBox["W", "+"], "\[RightArrow]", " ",
54 RowBox[{
55 SuperscriptBox["e", "+"], " ",
56 SubscriptBox["\[Nu]",
57 RowBox[{"e", " "}]], "can", " ", "be", " ", "written"}]}]}]], "Text",
58 CellChangeTimes->{{3.428941462807769*^9, 3.428941586342557*^9}}],
59
60Cell[CellGroupData[{
61
62Cell[BoxData[
63 RowBox[{"\[Sigma]p", "=",
64 RowBox[{"4",
65 RowBox[{"Pi", "/", "3"}], " ", "\[CapitalGamma]ud", " ",
66 RowBox[{"\[CapitalGamma]ev", " ", "/",
67 RowBox[{"(",
68 RowBox[{
69 RowBox[{
70 RowBox[{"(",
71 RowBox[{"sp", " ", "-", " ",
72 RowBox[{"mw", "^", "2"}]}], ")"}], "^", "2"}], " ", "+",
73 RowBox[{
74 RowBox[{"mw", "^", "2"}],
75 RowBox[{"\[CapitalGamma]tot", "^", "2"}]}]}], ")"}]}]}]}]], "Input",
76 CellChangeTimes->{{3.428941591115911*^9, 3.428941670219965*^9}, {
77 3.428941861583827*^9, 3.428941864412772*^9}}],
78
79Cell[BoxData[
80 FormBox[
81 FractionBox[
82 RowBox[{
83 "4", " ", "\[Pi]", " ", "\[CapitalGamma]ev", " ", "\[CapitalGamma]ud"}],
84 RowBox[{"3", " ",
85 RowBox[{"(",
86 RowBox[{
87 SuperscriptBox[
88 RowBox[{"(",
89 RowBox[{"sp", "-",
90 SuperscriptBox["mw", "2"]}], ")"}], "2"], "+",
91 RowBox[{
92 SuperscriptBox["mw", "2"], " ",
93 SuperscriptBox["\[CapitalGamma]tot", "2"]}]}], ")"}]}]],
94 TraditionalForm]], "Output",
95 CellChangeTimes->{3.428941671372353*^9, 3.4289418646862993`*^9,
96 3.428942151070161*^9, 3.428942597586359*^9}]
97}, Open ]],
98
99Cell["\<\
100where
101
102* \[CapitalGamma]ud is the partial width of the W into u d
103* \[CapitalGamma]ev is the partial width of the W into e \[Nu]
104*\[CapitalGamma]tot is the totel width
105
106At tree-level, we can write approximatively
107
108\[CapitalGamma]ev = \[CapitalGamma]
109\[CapitalGamma]ud = 3 \[CapitalGamma]
110\[CapitalGamma]tot = 9 \[CapitalGamma]
111\
112\>", "Text",
113 CellChangeTimes->{{3.428941673641809*^9, 3.428941779762138*^9}, {
114 3.4289418146980762`*^9, 3.428941823801876*^9}}],
115
116Cell[CellGroupData[{
117
118Cell[BoxData[
119 RowBox[{"\[Sigma]p", "=",
120 RowBox[{"\[Sigma]p", "/.",
121 RowBox[{"{",
122 RowBox[{
123 RowBox[{"\[CapitalGamma]ud", "\[Rule]",
124 RowBox[{"3", "\[CapitalGamma]"}]}], ",",
125 RowBox[{"\[CapitalGamma]ev", "->", "\[CapitalGamma]"}], ",",
126 RowBox[{"\[CapitalGamma]tot", "\[Rule]",
127 RowBox[{"9", "\[CapitalGamma]"}]}]}], "}"}]}]}]], "Input",
128 CellChangeTimes->{{3.4289418266514587`*^9, 3.4289418551565723`*^9}}],
129
130Cell[BoxData[
131 FormBox[
132 FractionBox[
133 RowBox[{"4", " ", "\[Pi]", " ",
134 SuperscriptBox["\[CapitalGamma]", "2"]}],
135 RowBox[{
136 SuperscriptBox[
137 RowBox[{"(",
138 RowBox[{"sp", "-",
139 SuperscriptBox["mw", "2"]}], ")"}], "2"], "+",
140 RowBox[{"81", " ",
141 SuperscriptBox["mw", "2"], " ",
142 SuperscriptBox["\[CapitalGamma]", "2"]}]}]], TraditionalForm]], "Output",\
143
144 CellChangeTimes->{{3.428941855531948*^9, 3.428941867874544*^9},
145 3.4289421541207457`*^9, 3.42894260090372*^9}]
146}, Open ]],
147
148Cell[TextData[{
149 "In the narrow width approximation, using\n\n\[Delta](x) = lim ",
150 Cell[BoxData[
151 FormBox[
152 FractionBox["1", "\[Pi]"], TraditionalForm]]],
153 " ",
154 Cell[BoxData[
155 FormBox[
156 FractionBox["\[Epsilon]",
157 RowBox[{
158 SuperscriptBox["x", "2"], "+",
159 SuperscriptBox["\[Epsilon]", "2"]}]], TraditionalForm]]],
160 "\n \[Epsilon]\[RightArrow]0\nwe have"
161}], "Text",
162 CellChangeTimes->{{3.428941885053216*^9, 3.4289420421627827`*^9}}],
163
164Cell[CellGroupData[{
165
166Cell[BoxData[
167 RowBox[{"\[Sigma]pNW", "=",
168 RowBox[{"9", "mw", " ", "\[CapitalGamma]", " ", "\[Sigma]pw", " ",
169 RowBox[{"DiracDelta", "[",
170 RowBox[{"sp", " ", "-", " ",
171 SuperscriptBox["mw", "2"]}], "]"}]}]}]], "Input",
172 CellChangeTimes->{{3.428942043694252*^9, 3.4289421458530283`*^9}, {
173 3.4289423809661303`*^9, 3.428942388280025*^9}, {3.4289426098473177`*^9,
174 3.428942609966979*^9}}],
175
176Cell[BoxData[
177 FormBox[
178 RowBox[{"9", " ", "mw", " ", "\[CapitalGamma]", " ", "\[Sigma]pw", " ",
179 RowBox[{
180 InterpretationBox["\[Delta]",
181 DiracDelta,
182 Editable->False,
183 Selectable->False], "(",
184 RowBox[{"sp", "-",
185 SuperscriptBox["mw", "2"]}], ")"}]}], TraditionalForm]], "Output",
186 CellChangeTimes->{{3.42894212937996*^9, 3.4289421574093657`*^9},
187 3.428942389290112*^9, 3.4289426110242777`*^9}]
188}, Open ]],
189
190Cell["\<\
191where \[Sigma]pw is \[Sigma]p evaluated at sp = mw^2\
192\>", "Text",
193 CellChangeTimes->{{3.4289421593269777`*^9, 3.428942176694747*^9}}]
194}, Open ]],
195
196Cell[CellGroupData[{
197
198Cell["The hadronic cross - section @ Tevatron", "Section",
199 CellChangeTimes->{{3.4289421851120443`*^9, 3.4289421961269503`*^9}}],
200
201Cell["We now fold the partonic cross - section with the pdf' s", "Text",
202 CellChangeTimes->{{3.428942204656125*^9, 3.4289422175994797`*^9}}],
203
204Cell[BoxData[{
205 RowBox[{
206 RowBox[{
207 RowBox[{"u", "[", "x_", "]"}], ":=",
208 RowBox[{"6",
209 RowBox[{
210 RowBox[{"(",
211 RowBox[{"1", "-", "x"}], ")"}], "^", "2"}]}]}],
212 ";"}], "\[IndentingNewLine]",
213 RowBox[{
214 RowBox[{
215 RowBox[{"d", "[", "x_", "]"}], ":=",
216 RowBox[{"3",
217 RowBox[{
218 RowBox[{"(",
219 RowBox[{"1", "-", "x"}], ")"}], "^", "2"}]}]}], ";"}]}], "Input",
220 CellChangeTimes->{{3.4289422228675957`*^9, 3.428942241154451*^9}}],
221
222Cell[TextData[{
223 "The limits of the integration are \n\nx1: [0, 1]\nx2: ",
224 Cell[BoxData[
225 FormBox["[", TraditionalForm]]],
226 " ",
227 Cell[BoxData[
228 FormBox[
229 RowBox[{
230 RowBox[{
231 RowBox[{
232 SuperscriptBox["mw", "2"], "/", "s"}], ",", " ", "1"}], "]"}],
233 TraditionalForm]]]
234}], "Text",
235 CellChangeTimes->{{3.428942264303419*^9, 3.428942324108612*^9}}],
236
237Cell["Integration in x1:", "Text",
238 CellChangeTimes->{{3.42894232586976*^9, 3.428942341099662*^9}}],
239
240Cell[CellGroupData[{
241
242Cell[BoxData[
243 RowBox[{
244 RowBox[{
245 RowBox[{"u", "[", "x1", "]"}],
246 RowBox[{"d", "[", "x2", "]"}], "\[Sigma]pNW"}], "/.",
247 RowBox[{"{",
248 RowBox[{"sp", "\[Rule]",
249 RowBox[{"s", " ", "x1", " ", "x2"}]}], "}"}]}]], "Input",
250 CellChangeTimes->{{3.4289431681668463`*^9, 3.4289431694285183`*^9}}],
251
252Cell[BoxData[
253 FormBox[
254 RowBox[{"162", " ", "mw", " ",
255 SuperscriptBox[
256 RowBox[{"(",
257 RowBox[{"1", "-", "x1"}], ")"}], "2"], " ",
258 SuperscriptBox[
259 RowBox[{"(",
260 RowBox[{"1", "-", "x2"}], ")"}], "2"], " ", "\[CapitalGamma]", " ",
261 "\[Sigma]pw", " ",
262 RowBox[{
263 InterpretationBox["\[Delta]",
264 DiracDelta,
265 Editable->False,
266 Selectable->False], "(",
267 RowBox[{
268 SuperscriptBox["mw", "2"], "-",
269 FractionBox["s",
270 RowBox[{"x1", " ", "x2"}]]}], ")"}]}], TraditionalForm]], "Output",
271 CellChangeTimes->{3.428943155912224*^9}]
272}, Open ]],
273
274Cell[CellGroupData[{
275
276Cell[BoxData[
277 RowBox[{"\[Sigma]", "=",
278 RowBox[{"Integrate", "[",
279 RowBox[{
280 RowBox[{
281 RowBox[{
282 RowBox[{"u", "[", "x1", "]"}],
283 RowBox[{"d", "[", "x2", "]"}], "\[Sigma]pNW"}], "/.",
284 RowBox[{"{",
285 RowBox[{"sp", "\[Rule]",
286 RowBox[{"s", " ", "x1", " ", "x2"}]}], "}"}]}], ",",
287 RowBox[{"{",
288 RowBox[{"x1", ",", "0", ",", "1"}], "}"}], ",",
289 RowBox[{"Assumptions", "\[Rule]",
290 RowBox[{
291 RowBox[{"(",
292 RowBox[{
293 RowBox[{
294 RowBox[{"mw", "^", "2"}], "/", "s"}], "<", "x2", "<", "1"}], ")"}], "&&",
295
296 RowBox[{"(",
297 RowBox[{"mw", ">", "0"}], ")"}], "&&",
298 RowBox[{"(",
299 RowBox[{"s", ">", "0"}], ")"}]}]}]}], "]"}]}]], "Input",
300 CellChangeTimes->{{3.428942342103321*^9, 3.428942375261568*^9}, {
301 3.428942406883278*^9, 3.428942462485826*^9}, {3.428942620176518*^9,
302 3.428942620422737*^9}, {3.428943124604053*^9, 3.42894321367962*^9}}],
303
304Cell[BoxData[
305 FormBox[
306 FractionBox[
307 RowBox[{"162", " ", "mw", " ",
308 SuperscriptBox[
309 RowBox[{"(",
310 RowBox[{"x2", "-", "1"}], ")"}], "2"], " ",
311 SuperscriptBox[
312 RowBox[{"(",
313 RowBox[{
314 SuperscriptBox["mw", "2"], "-",
315 RowBox[{"s", " ", "x2"}]}], ")"}], "2"], " ", "\[CapitalGamma]", " ",
316 "\[Sigma]pw"}],
317 RowBox[{
318 SuperscriptBox["s", "3"], " ",
319 SuperscriptBox["x2", "3"]}]], TraditionalForm]], "Output",
320 CellChangeTimes->{
321 3.428942400145671*^9, {3.4289424466850643`*^9, 3.428942465656705*^9},
322 3.4289425295180407`*^9, 3.428942621257402*^9, 3.428942755459882*^9,
323 3.4289429048806868`*^9, 3.428943215227071*^9}]
324}, Open ]],
325
326Cell["Integration in x2:", "Text",
327 CellChangeTimes->{{3.42894232586976*^9, 3.428942341099662*^9}, {
328 3.428942481054343*^9, 3.428942481204427*^9}}],
329
330Cell[CellGroupData[{
331
332Cell[BoxData[
333 RowBox[{"\[Sigma]", "=",
334 RowBox[{
335 RowBox[{
336 RowBox[{"Integrate", "[",
337 RowBox[{"\[Sigma]", ",",
338 RowBox[{"{",
339 RowBox[{"x2", ",",
340 RowBox[{
341 RowBox[{"mw", "^", "2"}], "/", "s"}], ",", "1"}], "}"}]}], "]"}], "//",
342 "Expand"}], "//", "Simplify"}]}]], "Input",
343 CellChangeTimes->{{3.428942483394573*^9, 3.42894252477667*^9}, {
344 3.4289427584498053`*^9, 3.4289427585792522`*^9}, {3.428942966307247*^9,
345 3.428942997247208*^9}, {3.428943223022503*^9, 3.428943273216826*^9}}],
346
347Cell[BoxData[
348 FormBox[
349 RowBox[{"-",
350 FractionBox[
351 RowBox[{
352 "162", " ", "mw", " ", "\[CapitalGamma]", " ", "\[Sigma]pw", " ",
353 RowBox[{"(",
354 RowBox[{
355 RowBox[{
356 RowBox[{"-", "3"}], " ",
357 SuperscriptBox["mw", "4"]}], "+",
358 RowBox[{"3", " ",
359 SuperscriptBox["s", "2"]}], "+",
360 RowBox[{
361 RowBox[{"(",
362 RowBox[{
363 SuperscriptBox["mw", "4"], "+",
364 RowBox[{"4", " ", "s", " ",
365 SuperscriptBox["mw", "2"]}], "+",
366 SuperscriptBox["s", "2"]}], ")"}], " ",
367 RowBox[{"log", "(",
368 FractionBox[
369 SuperscriptBox["mw", "2"], "s"], ")"}]}]}], ")"}]}],
370 SuperscriptBox["s", "3"]]}], TraditionalForm]], "Output",
371 CellChangeTimes->{{3.428942519914587*^9, 3.428942553721361*^9},
372 3.428942628541798*^9, 3.428942765391282*^9, 3.428942937050954*^9,
373 3.428942998707487*^9, {3.428943229928903*^9, 3.428943278718606*^9}}]
374}, Open ]],
375
376Cell[CellGroupData[{
377
378Cell[BoxData[
379 RowBox[{"\[Sigma]pw", "=",
380 RowBox[{"\[Sigma]p", "/.",
381 RowBox[{"sp", "\[Rule]",
382 RowBox[{"mw", "^", "2"}]}]}]}]], "Input",
383 CellChangeTimes->{{3.428943280891795*^9, 3.428943298023265*^9}}],
384
385Cell[BoxData[
386 FormBox[
387 FractionBox[
388 RowBox[{"4", " ", "\[Pi]"}],
389 RowBox[{"81", " ",
390 SuperscriptBox["mw", "2"]}]], TraditionalForm]], "Output",
391 CellChangeTimes->{{3.428943298212599*^9, 3.4289433089075317`*^9}}]
392}, Open ]],
393
394Cell[CellGroupData[{
395
396Cell[BoxData[
397 RowBox[{
398 RowBox[{"\[Sigma]", "//", "Expand"}], "//", "Simplify"}]], "Input",
399 CellChangeTimes->{{3.42894331467972*^9, 3.428943330039212*^9}}],
400
401Cell[BoxData[
402 FormBox[
403 RowBox[{"-",
404 FractionBox[
405 RowBox[{"8", " ", "\[Pi]", " ", "\[CapitalGamma]", " ",
406 RowBox[{"(",
407 RowBox[{
408 RowBox[{
409 RowBox[{"-", "3"}], " ",
410 SuperscriptBox["mw", "4"]}], "+",
411 RowBox[{"3", " ",
412 SuperscriptBox["s", "2"]}], "+",
413 RowBox[{
414 RowBox[{"(",
415 RowBox[{
416 SuperscriptBox["mw", "4"], "+",
417 RowBox[{"4", " ", "s", " ",
418 SuperscriptBox["mw", "2"]}], "+",
419 SuperscriptBox["s", "2"]}], ")"}], " ",
420 RowBox[{"log", "(",
421 FractionBox[
422 SuperscriptBox["mw", "2"], "s"], ")"}]}]}], ")"}]}],
423 RowBox[{"mw", " ",
424 SuperscriptBox["s", "3"]}]]}], TraditionalForm]], "Output",
425 CellChangeTimes->{{3.428943318260807*^9, 3.42894333060226*^9}}]
426}, Open ]]
427}, Open ]]
428}, Open ]]
429},
430WindowSize->{1280, 683},
431WindowMargins->{{0, Automatic}, {0, Automatic}},
432ShowSelection->True,
433Magnification->1.5,
434FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (April 20, 2007)",
435StyleDefinitions->"Default.nb"
436]
437(* End of Notebook Content *)
438
439(* Internal cache information *)
440(*CellTagsOutline
441CellTagsIndex->{}
442*)
443(*CellTagsIndex
444CellTagsIndex->{}
445*)
446(*NotebookFileOutline
447Notebook[{
448Cell[CellGroupData[{
449Cell[590, 23, 251, 5, 185, "Title"],
450Cell[CellGroupData[{
451Cell[866, 32, 36, 0, 51, "Subsection"],
452Cell[905, 34, 149, 3, 40, "Input"]
453}, Open ]],
454Cell[CellGroupData[{
455Cell[1091, 42, 113, 1, 100, "Section"],
456Cell[1207, 45, 531, 12, 39, "Text"],
457Cell[CellGroupData[{
458Cell[1763, 61, 574, 15, 40, "Input"],
459Cell[2340, 78, 573, 17, 75, "Output"]
460}, Open ]],
461Cell[2928, 98, 469, 15, 287, "Text"],
462Cell[CellGroupData[{
463Cell[3422, 117, 444, 10, 40, "Input"],
464Cell[3869, 129, 511, 15, 79, "Output"]
465}, Open ]],
466Cell[4395, 147, 462, 14, 144, "Text"],
467Cell[CellGroupData[{
468Cell[4882, 165, 403, 8, 49, "Input"],
469Cell[5288, 175, 427, 11, 49, "Output"]
470}, Open ]],
471Cell[5730, 189, 144, 3, 39, "Text"]
472}, Open ]],
473Cell[CellGroupData[{
474Cell[5911, 197, 128, 1, 100, "Section"],
475Cell[6042, 200, 140, 1, 39, "Text"],
476Cell[6185, 203, 464, 16, 66, "Input"],
477Cell[6652, 221, 365, 13, 115, "Text"],
478Cell[7020, 236, 99, 1, 39, "Text"],
479Cell[CellGroupData[{
480Cell[7144, 241, 304, 8, 40, "Input"],
481Cell[7451, 251, 584, 19, 59, "Output"]
482}, Open ]],
483Cell[CellGroupData[{
484Cell[8072, 275, 948, 26, 64, "Input"],
485Cell[9023, 303, 682, 19, 75, "Output"]
486}, Open ]],
487Cell[9720, 325, 148, 2, 39, "Text"],
488Cell[CellGroupData[{
489Cell[9893, 331, 532, 13, 40, "Input"],
490Cell[10428, 346, 954, 26, 83, "Output"]
491}, Open ]],
492Cell[CellGroupData[{
493Cell[11419, 377, 213, 5, 40, "Input"],
494Cell[11635, 384, 225, 6, 67, "Output"]
495}, Open ]],
496Cell[CellGroupData[{
497Cell[11897, 395, 158, 3, 40, "Input"],
498Cell[12058, 400, 815, 24, 83, "Output"]
499}, Open ]]
500}, Open ]]
501}, Open ]]
502}
503]
504*)
505
506(* End of internal cache information *)