RMSSM: MRSSM.nb

File MRSSM.nb, 11.1 KB (added by Benjamin Fuks, 12 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[ 11173, 347]
14NotebookOptionsPosition[ 9378, 279]
15NotebookOutlinePosition[ 9736, 295]
16CellTagsIndexPosition[ 9693, 292]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21Cell[BoxData[
22 RowBox[{"Quit", "[", "]"}]], "Input",
23 CellChangeTimes->{{3.511874428548209*^9, 3.511874430667358*^9}}],
24
25Cell[CellGroupData[{
26
27Cell["Setup", "Section"],
28
29Cell[BoxData[{
30 RowBox[{
31 RowBox[{"$OldDir", "=",
32 RowBox[{"Directory", "[", "]"}]}], ";"}], "\n",
33 RowBox[{
34 RowBox[{"$FeynRulesPath", "=",
35 RowBox[{
36 "SetDirectory", "[", "\"\<~/FeynRules/trunk/feynrules-development\>\"",
37 "]"}]}], ";"}], "\n",
38 RowBox[{"<<", "FeynRules`"}], "\n",
39 RowBox[{
40 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\n",
41 RowBox[{
42 RowBox[{"LoadModel", "[", "\"\<mrssm.fr\>\"", "]"}], ";"}], "\n",
43 RowBox[{
44 RowBox[{"$CalcLag", "=", "False"}], ";"}]}], "Input",
45 CellChangeTimes->{{3.510511714359152*^9, 3.5105117232569723`*^9}, {
46 3.511873979058074*^9, 3.511873982075786*^9}, {3.537793110077643*^9,
47 3.537793110571989*^9}}]
48}, Open ]],
49
50Cell[CellGroupData[{
51
52Cell["Lagrangian computation and update of the parameters", "Section",
53 CellChangeTimes->{{3.4912194701600723`*^9, 3.4912194715325117`*^9}, {
54 3.50686100488665*^9, 3.506861008141419*^9}, {3.506866546045188*^9,
55 3.50686654675209*^9}, 3.510511732200982*^9}],
56
57Cell[CellGroupData[{
58
59Cell["\<\
60Calculation of the Lagrangian itself. This can take ~10 minutes. Therefore, \
61it is saved in a file which can be read each time it is necessary.\
62\>", "Subsubsection",
63 CellChangeTimes->{{3.510511740551598*^9, 3.510511755990991*^9}}],
64
65Cell[BoxData[
66 RowBox[{
67 RowBox[{"If", "[",
68 RowBox[{"$CalcLag", ",",
69 RowBox[{
70 RowBox[{"lagr", "=", "Lag"}], ";", "\[IndentingNewLine]",
71 RowBox[{
72 RowBox[{"Definition", "[", "lagr", "]"}], ">>", "lagrangian.dat"}]}],
73 ",",
74 RowBox[{"<<", "lagrangian.dat"}]}], "]"}], ";"}]], "Input"]
75}, Open ]],
76
77Cell[CellGroupData[{
78
79Cell["Access to the different pieces of the Lagrangian", "Subsubsection",
80 CellChangeTimes->{{3.510511759479219*^9, 3.510511767623085*^9}}],
81
82Cell[BoxData[{
83 RowBox[{
84 RowBox[{"LVec", "=", "LVector"}], ";"}], "\[IndentingNewLine]",
85 RowBox[{
86 RowBox[{"LC", "=", "LChiral"}], ";"}], "\[IndentingNewLine]",
87 RowBox[{
88 RowBox[{"LW", "=", "LSuperW"}], ";"}], "\[IndentingNewLine]",
89 RowBox[{
90 RowBox[{"LSUSYBr", "=", "LSoft"}], ";"}]}], "Input",
91 CellChangeTimes->{
92 3.506861831832621*^9, {3.50686396987472*^9, 3.50686397956723*^9}, {
93 3.506864125221086*^9, 3.506864126560401*^9}, {3.506864198923688*^9,
94 3.506864211817175*^9}, {3.510511899307577*^9, 3.510511913390203*^9},
95 3.5121224264764137`*^9, {3.512122465053269*^9, 3.512122467303591*^9}}]
96}, Open ]],
97
98Cell[CellGroupData[{
99
100Cell["\<\
101Update of the parameters (not compatible with the Whizard interface, but \
102compatible with the CalcHep interface)\
103\>", "Subsubsection",
104 CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, {
105 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9,
106 3.510514606763693*^9}, {3.510514811212164*^9, 3.510514819639832*^9}}],
107
108Cell[BoxData[{
109 RowBox[{
110 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
111 RowBox[{
112 RowBox[{"ReadLHAFile", "[",
113 RowBox[{"Input", "\[Rule]", "\"\<spectrum.dat\>\""}], "]"}],
114 ";"}]}], "Input",
115 CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, {
116 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9,
117 3.510511927139632*^9}, 3.510511971795094*^9, {3.5439086684386797`*^9,
118 3.543908669966134*^9}}]
119}, Open ]]
120}, Open ]],
121
122Cell[CellGroupData[{
123
124Cell["Feynman rules calculation", "Section",
125 CellChangeTimes->{{3.508650859336681*^9, 3.5086508623247147`*^9}}],
126
127Cell[BoxData[
128 RowBox[{"rules", "=",
129 RowBox[{"FeynmanRules", "[",
130 RowBox[{"lagr", ",",
131 RowBox[{"ScreenOutput", "\[Rule]", "False"}], ",",
132 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}]}]], "Input",
133 CellChangeTimes->{{3.508061383518343*^9, 3.508061413431903*^9}, {
134 3.510571497855986*^9, 3.510571568813633*^9}, {3.510572173986182*^9,
135 3.5105721743854856`*^9}, 3.510572210161304*^9, {3.510572541184834*^9,
136 3.510572566899229*^9}}]
137}, Open ]],
138
139Cell[CellGroupData[{
140
141Cell[TextData[{
142 Cell[BoxData[
143 FormBox[
144 StyleBox[
145 RowBox[{"L",
146 StyleBox[
147 AdjustmentBox["A",
148 BoxBaselineShift->-0.4,
149 BoxMargins->{{-0.5, -0.3}, {0, 0}}],
150 FontSize->Smaller], "T",
151 AdjustmentBox["E",
152 BoxBaselineShift->0.5,
153 BoxMargins->{{-0.3, 0}, {0, 0}}], "X"}],
154 SingleLetterItalics->False], TraditionalForm]]],
155 " Output"
156}], "Section",
157 CellChangeTimes->{{3.510516585335*^9, 3.510516590423654*^9}}],
158
159Cell[BoxData[
160 RowBox[{"WriteLaTeXOutput", "[", "]"}]], "Input",
161 CellChangeTimes->{{3.510516592624229*^9, 3.510516605419121*^9}}]
162}, Open ]],
163
164Cell[CellGroupData[{
165
166Cell["Generating MC model files with optimized restrictions", "Section",
167 CellChangeTimes->{{3.51051202409053*^9, 3.5105120509193287`*^9}}],
168
169Cell[CellGroupData[{
170
171Cell["Optimization", "Subsubsection",
172 CellChangeTimes->{{3.508650844785081*^9, 3.508650846244998*^9}}],
173
174Cell[BoxData[
175 RowBox[{
176 RowBox[{"WriteRestrictionFile", "[", "]"}], ";", " ",
177 RowBox[{"LoadRestriction", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";", " ",
178 RowBox[{"DeleteFile", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";"}]], "Input",\
179
180 CellChangeTimes->{{3.5080579535629387`*^9, 3.508057961830206*^9}, {
181 3.508058109308621*^9, 3.508058122229617*^9}}]
182}, Open ]],
183
184Cell[CellGroupData[{
185
186Cell["UFO interface", "Subsubsection",
187 CellChangeTimes->{{3.508650849632307*^9, 3.508650853575732*^9}}],
188
189Cell[BoxData[{
190 RowBox[{
191 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
192 RowBox[{
193 RowBox[{"WriteUFO", "[",
194 RowBox[{"LagFeynman", ",", " ",
195 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
196
197 CellChangeTimes->{{3.5073089286898727`*^9, 3.507308937458375*^9}, {
198 3.507311910076086*^9, 3.507311939653351*^9}, {3.507398897853717*^9,
199 3.507398901626676*^9}, {3.507399273736373*^9, 3.507399274557005*^9}, {
200 3.510512090662549*^9, 3.510512091990087*^9}}]
201}, Open ]],
202
203Cell[CellGroupData[{
204
205Cell["CalcHEP model files", "Subsubsection",
206 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}}],
207
208Cell[BoxData[{
209 RowBox[{
210 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
211 RowBox[{
212 RowBox[{"WriteCHOutput", "[",
213 RowBox[{"lagr", ",", " ",
214 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}], ",",
215 RowBox[{"CHSimplify", "\[Rule]", "False"}]}], "]"}], ";"}]}], "Input",
216 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
217 3.5105127146047497`*^9, 3.510512717826407*^9}}]
218}, Open ]],
219
220Cell[CellGroupData[{
221
222Cell["MadGraph model files", "Subsubsection",
223 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
224 3.510513948444434*^9, 3.510513949301979*^9}}],
225
226Cell[BoxData[{
227 RowBox[{
228 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
229 RowBox[{
230 RowBox[{"WriteMGOutput", "[",
231 RowBox[{"lagr", ",", " ",
232 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
233
234 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
235 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
236 3.510513955605171*^9}}]
237}, Open ]],
238
239Cell[CellGroupData[{
240
241Cell["Whizard model files", "Subsubsection",
242 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
243 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
244 3.510514891766358*^9}}],
245
246Cell[BoxData[{
247 RowBox[{
248 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
249 RowBox[{
250 RowBox[{"WriteWOOutput", "[",
251 RowBox[{"lagr", ",", " ",
252 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
253
254 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
255 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
256 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
257 3.510572235119359*^9, 3.51057223551847*^9}, {3.510572377588068*^9,
258 3.510572378543983*^9}, {3.510572531362316*^9, 3.510572531820734*^9}}]
259}, Open ]],
260
261Cell[CellGroupData[{
262
263Cell["FeynArts model files", "Subsubsection",
264 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
265 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
266 3.510514891766358*^9}, {3.5105162477375813`*^9, 3.5105162494145308`*^9}}],
267
268Cell[BoxData[{
269 RowBox[{
270 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
271 RowBox[{
272 RowBox[{"WriteFeynArtsOutput", "[", "lagr", "]"}], ";"}]}], "Input",
273 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
274 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
275 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
276 3.5105162369063063`*^9, 3.510516243526369*^9}}]
277}, Open ]]
278}, Open ]]
279},
280WindowSize->{1440, 825},
281WindowMargins->{{3, Automatic}, {Automatic, 23}},
282ShowSelection->True,
283FrontEndVersion->"7.0 for Linux x86 (64-bit) (February 25, 2009)",
284StyleDefinitions->"Default.nb"
285]
286(* End of Notebook Content *)
287
288(* Internal cache information *)
289(*CellTagsOutline
290CellTagsIndex->{}
291*)
292(*CellTagsIndex
293CellTagsIndex->{}
294*)
295(*NotebookFileOutline
296Notebook[{
297Cell[545, 20, 118, 2, 32, "Input"],
298Cell[CellGroupData[{
299Cell[688, 26, 24, 0, 73, "Section"],
300Cell[715, 28, 684, 18, 143, "Input"]
301}, Open ]],
302Cell[CellGroupData[{
303Cell[1436, 51, 260, 3, 73, "Section"],
304Cell[CellGroupData[{
305Cell[1721, 58, 243, 4, 29, "Subsubsection"],
306Cell[1967, 64, 316, 9, 55, "Input"]
307}, Open ]],
308Cell[CellGroupData[{
309Cell[2320, 78, 140, 1, 29, "Subsubsection"],
310Cell[2463, 81, 615, 13, 99, "Input"]
311}, Open ]],
312Cell[CellGroupData[{
313Cell[3115, 99, 360, 6, 29, "Subsubsection"],
314Cell[3478, 107, 470, 10, 55, "Input"]
315}, Open ]]
316}, Open ]],
317Cell[CellGroupData[{
318Cell[3997, 123, 112, 1, 73, "Section"],
319Cell[4112, 126, 463, 9, 32, "Input"]
320}, Open ]],
321Cell[CellGroupData[{
322Cell[4612, 140, 461, 16, 75, "Section"],
323Cell[5076, 158, 130, 2, 32, "Input"]
324}, Open ]],
325Cell[CellGroupData[{
326Cell[5243, 165, 139, 1, 73, "Section"],
327Cell[CellGroupData[{
328Cell[5407, 170, 103, 1, 29, "Subsubsection"],
329Cell[5513, 173, 356, 7, 32, "Input"]
330}, Open ]],
331Cell[CellGroupData[{
332Cell[5906, 185, 104, 1, 29, "Subsubsection"],
333Cell[6013, 188, 520, 11, 55, "Input"]
334}, Open ]],
335Cell[CellGroupData[{
336Cell[6570, 204, 110, 1, 29, "Subsubsection"],
337Cell[6683, 207, 429, 9, 55, "Input"]
338}, Open ]],
339Cell[CellGroupData[{
340Cell[7149, 221, 160, 2, 29, "Subsubsection"],
341Cell[7312, 225, 426, 10, 55, "Input"]
342}, Open ]],
343Cell[CellGroupData[{
344Cell[7775, 240, 208, 3, 29, "Subsubsection"],
345Cell[7986, 245, 615, 12, 55, "Input"]
346}, Open ]],
347Cell[CellGroupData[{
348Cell[8638, 262, 259, 3, 29, "Subsubsection"],
349Cell[8900, 267, 450, 8, 55, "Input"]
350}, Open ]]
351}, Open ]]
352}
353]
354*)
355
356(* End of internal cache information *)