RMSSM: MRSSM.nb

File MRSSM.nb, 11.1 KB (added by BenjF, 5 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 *)