NMSSM: NMSSM.nb

File NMSSM.nb, 13.0 KB (added by BenjF, 8 years ago)

Example notebook for the NMSSM

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