NMSSM: NMSSM.nb

File NMSSM.nb, 13.0 KB (added by Benjamin Fuks, 14 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 *)