MSSM: MSSM.nb

File MSSM.nb, 13.0 KB (added by BenjF, 8 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 *)