MSSM: MSSM.nb

File MSSM.nb, 13.0 KB (added by Benjamin Fuks, 14 years ago)

Example notebook

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[ 13172, 405]
14NotebookOptionsPosition[ 11175, 330]
15NotebookOutlinePosition[ 11533, 346]
16CellTagsIndexPosition[ 11490, 343]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell["Setup", "Section"],
24
25Cell[BoxData[{
26 RowBox[{
27 RowBox[{"$OldDir", "=",
28 RowBox[{"Directory", "[", "]"}]}], ";"}], "\n",
29 RowBox[{
30 RowBox[{"$FeynRulesPath", "=",
31 RowBox[{
32 "SetDirectory", "[", "\"\<~/FeynRules/trunk/feynrules-development\>\"",
33 "]"}]}], ";"}], "\n",
34 RowBox[{"<<", "FeynRules`"}], "\n",
35 RowBox[{
36 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\n",
37 RowBox[{
38 RowBox[{"LoadModel", "[", "\"\<mssm.fr\>\"", "]"}], ";"}], "\n",
39 RowBox[{
40 RowBox[{"$CalcLag", "=", "False"}], ";"}]}], "Input",
41 CellChangeTimes->{{3.510511714359152*^9, 3.5105117232569723`*^9}}]
42}, Open ]],
43
44Cell[CellGroupData[{
45
46Cell["Lagrangian computation and update of the parameters", "Section",
47 CellChangeTimes->{{3.4912194701600723`*^9, 3.4912194715325117`*^9}, {
48 3.50686100488665*^9, 3.506861008141419*^9}, {3.506866546045188*^9,
49 3.50686654675209*^9}, 3.510511732200982*^9}],
50
51Cell[CellGroupData[{
52
53Cell["\<\
54Calculation of the Lagrangian itself. This can take ~10 minutes. Therefore, \
55it is saved in a file which can be read each time it is necessary.\
56\>", "Subsubsection",
57 CellChangeTimes->{{3.510511740551598*^9, 3.510511755990991*^9}}],
58
59Cell[BoxData[
60 RowBox[{
61 RowBox[{"If", "[",
62 RowBox[{"$CalcLag", ",",
63 RowBox[{
64 RowBox[{"lagr", "=", "Lag"}], ";", "\[IndentingNewLine]",
65 RowBox[{
66 RowBox[{"Definition", "[", "lagr", "]"}], ">>", "lagrangian.dat"}]}],
67 ",",
68 RowBox[{"<<", "lagrangian.dat"}]}], "]"}], ";"}]], "Input"]
69}, Open ]],
70
71Cell[CellGroupData[{
72
73Cell["Access to the different pieces of the Lagrangian", "Subsubsection",
74 CellChangeTimes->{{3.510511759479219*^9, 3.510511767623085*^9}}],
75
76Cell[BoxData[{
77 RowBox[{
78 RowBox[{"LVec", "=", "LVector"}], ";"}], "\[IndentingNewLine]",
79 RowBox[{
80 RowBox[{"LC", "=", "LChiral"}], ";"}], "\[IndentingNewLine]",
81 RowBox[{
82 RowBox[{"LW", "=", "LSuperW"}], ";"}], "\[IndentingNewLine]",
83 RowBox[{
84 RowBox[{"LSUSYBr", "=", "LSoft"}], ";"}], "\[IndentingNewLine]",
85 RowBox[{
86 RowBox[{"LGaugeFixing", "=", "LFeynmanGFix"}], ";"}]}], "Input",
87 CellChangeTimes->{
88 3.506861831832621*^9, {3.50686396987472*^9, 3.50686397956723*^9}, {
89 3.506864125221086*^9, 3.506864126560401*^9}, {3.506864198923688*^9,
90 3.506864211817175*^9}, {3.510511899307577*^9, 3.510511913390203*^9}}]
91}, Open ]],
92
93Cell[CellGroupData[{
94
95Cell["\<\
96Unitarity gauge (no ghosts and Goldstones, if the flag $FeynmanGauge was not \
97set correctly)\
98\>", "Subsubsection",
99 CellChangeTimes->{{3.508650830315112*^9, 3.508650832348791*^9}, {
100 3.510511936427951*^9, 3.510511947862441*^9}}],
101
102Cell[BoxData[
103 RowBox[{
104 RowBox[{"LagFeynman", "=",
105 RowBox[{"lagr", "/.",
106 RowBox[{"{",
107 RowBox[{
108 RowBox[{
109 RowBox[{"ghG", "[", "__", "]"}], "\[Rule]", "0"}], ",", " ",
110 RowBox[{
111 RowBox[{"ghGbar", "[", "__", "]"}], "\[Rule]", "0"}], ",",
112 RowBox[{"ghWp", "\[Rule]", "0"}], ",",
113 RowBox[{"ghWpbar", "\[Rule]", "0"}], ",",
114 RowBox[{"ghWmbar", "\[Rule]", "0"}], ",",
115 RowBox[{"ghWm", "\[Rule]", "0"}], ",",
116 RowBox[{"ghZ", "\[Rule]", "0"}], ",",
117 RowBox[{"ghZbar", "\[Rule]", "0"}], ",",
118 RowBox[{"ghA", "\[Rule]", "0"}], ",",
119 RowBox[{"ghAbar", "\[Rule]", "0"}], ",", " ",
120 RowBox[{"G0", "\[Rule]", "0"}], ",",
121 RowBox[{"GP", "\[Rule]", "0"}], ",",
122 RowBox[{"GPbar", "\[Rule]", "0"}]}], "}"}]}]}], ";"}]], "Input",
123 CellChangeTimes->{{3.507308829401106*^9, 3.5073089177178507`*^9},
124 3.507310603266782*^9}]
125}, Open ]],
126
127Cell[CellGroupData[{
128
129Cell["\<\
130Update of the parameters (not compatible with the Whizard interface, but \
131compatible with the CalcHep interface)\
132\>", "Subsubsection",
133 CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, {
134 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9,
135 3.510514606763693*^9}, {3.510514811212164*^9, 3.510514819639832*^9}}],
136
137Cell[BoxData[{
138 RowBox[{
139 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
140 RowBox[{
141 RowBox[{"ReadLHAFile", "[",
142 RowBox[{"Input", "\[Rule]", "\"\<sps1a.dat\>\""}], "]"}], ";"}]}], "Input",\
143
144 CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, {
145 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9,
146 3.510511927139632*^9}, 3.510511971795094*^9}]
147}, Open ]],
148
149Cell[CellGroupData[{
150
151Cell["\<\
152Update of the parameters (compatible the Whizard interface, but not \
153compatible with the Calchep interface)\
154\>", "Subsubsection",
155 CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, {
156 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9,
157 3.510514606763693*^9}, {3.5105147982031612`*^9, 3.510514828422709*^9}}],
158
159Cell[BoxData[{
160 RowBox[{
161 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
162 RowBox[{
163 RowBox[{"ReadLHAFile", "[",
164 RowBox[{"Input", "\[Rule]", "\"\<sps1a_wo.dat\>\""}], "]"}],
165 ";"}]}], "Input",
166 CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, {
167 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9,
168 3.510511927139632*^9}, 3.510511971795094*^9, {3.510514852406598*^9,
169 3.510514852821206*^9}}]
170}, Open ]]
171}, Open ]],
172
173Cell[CellGroupData[{
174
175Cell["Feynman rules calculation", "Section",
176 CellChangeTimes->{{3.508650859336681*^9, 3.5086508623247147`*^9}}],
177
178Cell[BoxData[
179 RowBox[{"rules", "=",
180 RowBox[{"FeynmanRules", "[",
181 RowBox[{"lagr", ",",
182 RowBox[{"ScreenOutput", "\[Rule]", "False"}], ",",
183 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}]}]], "Input",
184 CellChangeTimes->{{3.508061383518343*^9, 3.508061413431903*^9}, {
185 3.510571497855986*^9, 3.510571568813633*^9}, {3.510572173986182*^9,
186 3.5105721743854856`*^9}, 3.510572210161304*^9, {3.510572541184834*^9,
187 3.510572566899229*^9}}]
188}, Open ]],
189
190Cell[CellGroupData[{
191
192Cell[TextData[{
193 Cell[BoxData[
194 FormBox[
195 StyleBox[
196 RowBox[{"L",
197 StyleBox[
198 AdjustmentBox["A",
199 BoxBaselineShift->-0.4,
200 BoxMargins->{{-0.5, -0.3}, {0, 0}}],
201 FontSize->Smaller], "T",
202 AdjustmentBox["E",
203 BoxBaselineShift->0.5,
204 BoxMargins->{{-0.3, 0}, {0, 0}}], "X"}],
205 SingleLetterItalics->False], TraditionalForm]]],
206 " Output"
207}], "Section",
208 CellChangeTimes->{{3.510516585335*^9, 3.510516590423654*^9}}],
209
210Cell[BoxData[
211 RowBox[{"WriteLaTeXOutput", "[", "]"}]], "Input",
212 CellChangeTimes->{{3.510516592624229*^9, 3.510516605419121*^9}}]
213}, Open ]],
214
215Cell[CellGroupData[{
216
217Cell["Generating MC model files with optimized restrictions", "Section",
218 CellChangeTimes->{{3.51051202409053*^9, 3.5105120509193287`*^9}}],
219
220Cell[CellGroupData[{
221
222Cell["Optimization", "Subsubsection",
223 CellChangeTimes->{{3.508650844785081*^9, 3.508650846244998*^9}}],
224
225Cell[BoxData[
226 RowBox[{
227 RowBox[{"WriteRestrictionFile", "[", "]"}], ";", " ",
228 RowBox[{"LoadRestriction", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";", " ",
229 RowBox[{"DeleteFile", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";"}]], "Input",\
230
231 CellChangeTimes->{{3.5080579535629387`*^9, 3.508057961830206*^9}, {
232 3.508058109308621*^9, 3.508058122229617*^9}}]
233}, Open ]],
234
235Cell[CellGroupData[{
236
237Cell["UFO interface", "Subsubsection",
238 CellChangeTimes->{{3.508650849632307*^9, 3.508650853575732*^9}}],
239
240Cell[BoxData[{
241 RowBox[{
242 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
243 RowBox[{
244 RowBox[{"WriteUFO", "[",
245 RowBox[{"LagFeynman", ",", " ",
246 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
247
248 CellChangeTimes->{{3.5073089286898727`*^9, 3.507308937458375*^9}, {
249 3.507311910076086*^9, 3.507311939653351*^9}, {3.507398897853717*^9,
250 3.507398901626676*^9}, {3.507399273736373*^9, 3.507399274557005*^9}, {
251 3.510512090662549*^9, 3.510512091990087*^9}}]
252}, Open ]],
253
254Cell[CellGroupData[{
255
256Cell["CalcHEP model files", "Subsubsection",
257 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}}],
258
259Cell[BoxData[{
260 RowBox[{
261 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
262 RowBox[{
263 RowBox[{"WriteCHOutput", "[",
264 RowBox[{"lagr", ",", " ",
265 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}], ",",
266 RowBox[{"CHSimplify", "\[Rule]", "False"}]}], "]"}], ";"}]}], "Input",
267 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
268 3.5105127146047497`*^9, 3.510512717826407*^9}}]
269}, Open ]],
270
271Cell[CellGroupData[{
272
273Cell["MadGraph model files", "Subsubsection",
274 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
275 3.510513948444434*^9, 3.510513949301979*^9}}],
276
277Cell[BoxData[{
278 RowBox[{
279 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
280 RowBox[{
281 RowBox[{"WriteMGOutput", "[",
282 RowBox[{"lagr", ",", " ",
283 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
284
285 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
286 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
287 3.510513955605171*^9}}]
288}, Open ]],
289
290Cell[CellGroupData[{
291
292Cell["Whizard model files", "Subsubsection",
293 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
294 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
295 3.510514891766358*^9}}],
296
297Cell[BoxData[{
298 RowBox[{
299 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
300 RowBox[{
301 RowBox[{"WriteWOOutput", "[",
302 RowBox[{"lagr", ",", " ",
303 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
304
305 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
306 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
307 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
308 3.510572235119359*^9, 3.51057223551847*^9}, {3.510572377588068*^9,
309 3.510572378543983*^9}, {3.510572531362316*^9, 3.510572531820734*^9}}]
310}, Open ]],
311
312Cell[CellGroupData[{
313
314Cell["FeynArts model files", "Subsubsection",
315 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
316 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
317 3.510514891766358*^9}, {3.5105162477375813`*^9, 3.5105162494145308`*^9}}],
318
319Cell[BoxData[{
320 RowBox[{
321 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
322 RowBox[{
323 RowBox[{"WriteFeynArtsOutput", "[", "lagr", "]"}], ";"}]}], "Input",
324 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
325 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
326 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
327 3.5105162369063063`*^9, 3.510516243526369*^9}}]
328}, Open ]]
329}, Open ]]
330},
331WindowSize->{1425, 913},
332WindowMargins->{{3, Automatic}, {Automatic, 23}},
333ShowSelection->True,
334FrontEndVersion->"7.0 for Linux x86 (64-bit) (February 25, 2009)",
335StyleDefinitions->"Default.nb"
336]
337(* End of Notebook Content *)
338
339(* Internal cache information *)
340(*CellTagsOutline
341CellTagsIndex->{}
342*)
343(*CellTagsIndex
344CellTagsIndex->{}
345*)
346(*NotebookFileOutline
347Notebook[{
348Cell[CellGroupData[{
349Cell[567, 22, 24, 0, 73, "Section"],
350Cell[594, 24, 585, 16, 143, "Input"]
351}, Open ]],
352Cell[CellGroupData[{
353Cell[1216, 45, 260, 3, 73, "Section"],
354Cell[CellGroupData[{
355Cell[1501, 52, 243, 4, 29, "Subsubsection"],
356Cell[1747, 58, 316, 9, 55, "Input"]
357}, Open ]],
358Cell[CellGroupData[{
359Cell[2100, 72, 140, 1, 29, "Subsubsection"],
360Cell[2243, 75, 631, 14, 121, "Input"]
361}, Open ]],
362Cell[CellGroupData[{
363Cell[2911, 94, 241, 5, 29, "Subsubsection"],
364Cell[3155, 101, 911, 22, 32, "Input"]
365}, Open ]],
366Cell[CellGroupData[{
367Cell[4103, 128, 360, 6, 29, "Subsubsection"],
368Cell[4466, 136, 414, 9, 55, "Input"]
369}, Open ]],
370Cell[CellGroupData[{
371Cell[4917, 150, 357, 6, 29, "Subsubsection"],
372Cell[5277, 158, 468, 10, 55, "Input"]
373}, Open ]]
374}, Open ]],
375Cell[CellGroupData[{
376Cell[5794, 174, 112, 1, 73, "Section"],
377Cell[5909, 177, 463, 9, 32, "Input"]
378}, Open ]],
379Cell[CellGroupData[{
380Cell[6409, 191, 461, 16, 75, "Section"],
381Cell[6873, 209, 130, 2, 32, "Input"]
382}, Open ]],
383Cell[CellGroupData[{
384Cell[7040, 216, 139, 1, 73, "Section"],
385Cell[CellGroupData[{
386Cell[7204, 221, 103, 1, 29, "Subsubsection"],
387Cell[7310, 224, 356, 7, 32, "Input"]
388}, Open ]],
389Cell[CellGroupData[{
390Cell[7703, 236, 104, 1, 29, "Subsubsection"],
391Cell[7810, 239, 520, 11, 55, "Input"]
392}, Open ]],
393Cell[CellGroupData[{
394Cell[8367, 255, 110, 1, 29, "Subsubsection"],
395Cell[8480, 258, 429, 9, 55, "Input"]
396}, Open ]],
397Cell[CellGroupData[{
398Cell[8946, 272, 160, 2, 29, "Subsubsection"],
399Cell[9109, 276, 426, 10, 55, "Input"]
400}, Open ]],
401Cell[CellGroupData[{
402Cell[9572, 291, 208, 3, 29, "Subsubsection"],
403Cell[9783, 296, 615, 12, 55, "Input"]
404}, Open ]],
405Cell[CellGroupData[{
406Cell[10435, 313, 259, 3, 29, "Subsubsection"],
407Cell[10697, 318, 450, 8, 55, "Input"]
408}, Open ]]
409}, Open ]]
410}
411]
412*)
413
414(* End of internal cache information *)