RPVMSSM: RPV.nb

File RPV.nb, 15.1 KB (added by Benjamin Fuks, 14 years ago)

Example notebook for the RPV-MSSM

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[ 15241, 447]
14NotebookOptionsPosition[ 13235, 372]
15NotebookOutlinePosition[ 13594, 388]
16CellTagsIndexPosition[ 13551, 385]
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", "[", "\"\<rpvmssm.fr\>\"", "]"}], ";"}], "\n",
39 RowBox[{
40 RowBox[{"$CalcLag", "=", "False"}], ";"}]}], "Input",
41 CellChangeTimes->{{3.510511714359152*^9, 3.5105117232569723`*^9}, {
42 3.51167265039277*^9, 3.511672653856179*^9}, {3.511680648528541*^9,
43 3.5116806490432587`*^9}, {3.511686519419437*^9, 3.511686519752368*^9}, {
44 3.511690090879174*^9, 3.511690091381048*^9}, {3.5117562062782717`*^9,
45 3.511756206602289*^9}, {3.511757059875284*^9, 3.5117570603649693`*^9}}]
46}, Open ]],
47
48Cell[CellGroupData[{
49
50Cell["Lagrangian computation and update of the parameters", "Section",
51 CellChangeTimes->{{3.4912194701600723`*^9, 3.4912194715325117`*^9}, {
52 3.50686100488665*^9, 3.506861008141419*^9}, {3.506866546045188*^9,
53 3.50686654675209*^9}, 3.510511732200982*^9}],
54
55Cell[CellGroupData[{
56
57Cell["\<\
58Calculation of the Lagrangian itself. This can take ~10 minutes. Therefore, \
59it is saved in a file which can be read each time it is necessary.\
60\>", "Subsubsection",
61 CellChangeTimes->{{3.510511740551598*^9, 3.510511755990991*^9}}],
62
63Cell[BoxData[
64 RowBox[{
65 RowBox[{"If", "[",
66 RowBox[{"$CalcLag", ",",
67 RowBox[{
68 RowBox[{"lagr", "=", "Lag"}], ";", "\[IndentingNewLine]",
69 RowBox[{
70 RowBox[{"Definition", "[", "lagr", "]"}], ">>", "lagrangian.dat"}]}],
71 ",",
72 RowBox[{"<<", "lagrangian.dat"}]}], "]"}], ";"}]], "Input"]
73}, Open ]],
74
75Cell[CellGroupData[{
76
77Cell["Access to the different pieces of the Lagrangian", "Subsubsection",
78 CellChangeTimes->{{3.510511759479219*^9, 3.510511767623085*^9}}],
79
80Cell[BoxData[{
81 RowBox[{
82 RowBox[{"LVec", "=", "LVector"}], ";"}], "\[IndentingNewLine]",
83 RowBox[{
84 RowBox[{"LC", "=", "LChiral"}], ";"}], "\[IndentingNewLine]",
85 RowBox[{
86 RowBox[{"LW", "=", "LSuperW"}], ";"}], "\[IndentingNewLine]",
87 RowBox[{
88 RowBox[{"LSUSYBr", "=", "LSoft"}], ";"}], "\[IndentingNewLine]",
89 RowBox[{
90 RowBox[{"LGaugeFixing", "=", "LFeynmanGFix"}], ";"}]}], "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.511767566647472*^9, 3.51176757020861*^9}, {3.5117687354968033`*^9,
96 3.511768736624299*^9}}]
97}, Open ]],
98
99Cell[CellGroupData[{
100
101Cell["\<\
102Unitarity gauge (no ghosts and Goldstones, if the flag $FeynmanGauge was not \
103set correctly)\
104\>", "Subsubsection",
105 CellChangeTimes->{{3.508650830315112*^9, 3.508650832348791*^9}, {
106 3.510511936427951*^9, 3.510511947862441*^9}}],
107
108Cell[BoxData[
109 RowBox[{
110 RowBox[{"LagFeynman", "=",
111 RowBox[{"lagr", "/.",
112 RowBox[{"{",
113 RowBox[{
114 RowBox[{
115 RowBox[{"ghG", "[", "__", "]"}], "\[Rule]", "0"}], ",", " ",
116 RowBox[{
117 RowBox[{"ghGbar", "[", "__", "]"}], "\[Rule]", "0"}], ",",
118 RowBox[{"ghWp", "\[Rule]", "0"}], ",",
119 RowBox[{"ghWpbar", "\[Rule]", "0"}], ",",
120 RowBox[{"ghWmbar", "\[Rule]", "0"}], ",",
121 RowBox[{"ghWm", "\[Rule]", "0"}], ",",
122 RowBox[{"ghZ", "\[Rule]", "0"}], ",",
123 RowBox[{"ghZbar", "\[Rule]", "0"}], ",",
124 RowBox[{"ghA", "\[Rule]", "0"}], ",",
125 RowBox[{"ghAbar", "\[Rule]", "0"}], ",", " ",
126 RowBox[{"G0", "\[Rule]", "0"}], ",",
127 RowBox[{"GP", "\[Rule]", "0"}], ",",
128 RowBox[{"GPbar", "\[Rule]", "0"}]}], "}"}]}]}], ";"}]], "Input",
129 CellChangeTimes->{{3.507308829401106*^9, 3.5073089177178507`*^9},
130 3.507310603266782*^9}]
131}, Open ]],
132
133Cell[CellGroupData[{
134
135Cell["\<\
136Update of the parameters (not compatible with the Whizard interface, but \
137compatible with the CalcHep interface)\
138\>", "Subsubsection",
139 CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, {
140 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9,
141 3.510514606763693*^9}, {3.510514811212164*^9, 3.510514819639832*^9}}],
142
143Cell[BoxData[{
144 RowBox[{
145 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
146 RowBox[{
147 RowBox[{"ReadLHAFile", "[",
148 RowBox[{"Input", "\[Rule]", "\"\<af1.dat\>\""}], "]"}], ";"}]}], "Input",
149 CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, {
150 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9,
151 3.510511927139632*^9}, 3.510511971795094*^9, {3.51168065590665*^9,
152 3.511680656147051*^9}}]
153}, Open ]],
154
155Cell[CellGroupData[{
156
157Cell["\<\
158Update of the parameters (compatible the Whizard interface, but not \
159compatible with the Calchep interface)\
160\>", "Subsubsection",
161 CellChangeTimes->{{3.508650840502183*^9, 3.508650841544434*^9}, {
162 3.510511953402149*^9, 3.5105119585185213`*^9}, {3.5105145827435713`*^9,
163 3.510514606763693*^9}, {3.5105147982031612`*^9, 3.510514828422709*^9}}],
164
165Cell[BoxData[{
166 RowBox[{
167 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
168 RowBox[{
169 RowBox[{"ReadLHAFile", "[",
170 RowBox[{"Input", "\[Rule]", "\"\<af1_wo.dat\>\""}], "]"}], ";"}]}], "Input",\
171
172 CellChangeTimes->{{3.507310458304277*^9, 3.507310473976326*^9}, {
173 3.507310607848939*^9, 3.507310608185873*^9}, {3.510511925987886*^9,
174 3.510511927139632*^9}, 3.510511971795094*^9, {3.510514852406598*^9,
175 3.510514852821206*^9}, {3.511691760941457*^9, 3.511691761294777*^9}, {
176 3.511767574101495*^9, 3.511767577751384*^9}, {3.511768731486969*^9,
177 3.511768733085568*^9}}]
178}, Open ]]
179}, Open ]],
180
181Cell[CellGroupData[{
182
183Cell["Feynman rules calculation", "Section",
184 CellChangeTimes->{{3.508650859336681*^9, 3.5086508623247147`*^9}}],
185
186Cell[BoxData[
187 RowBox[{"rules", "=",
188 RowBox[{"FeynmanRules", "[",
189 RowBox[{"lagr", ",",
190 RowBox[{"ScreenOutput", "\[Rule]", "False"}], ",",
191 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}]}]], "Input",
192 CellChangeTimes->{{3.508061383518343*^9, 3.508061413431903*^9}, {
193 3.510571497855986*^9, 3.510571568813633*^9}, {3.510572173986182*^9,
194 3.5105721743854856`*^9}, 3.510572210161304*^9, {3.510572541184834*^9,
195 3.510572566899229*^9}, {3.511767581887206*^9, 3.5117675831212263`*^9}, {
196 3.511768728775072*^9, 3.511768729278853*^9}}]
197}, Open ]],
198
199Cell[CellGroupData[{
200
201Cell[TextData[{
202 Cell[BoxData[
203 FormBox[
204 StyleBox[
205 RowBox[{"L",
206 StyleBox[
207 AdjustmentBox["A",
208 BoxBaselineShift->-0.4,
209 BoxMargins->{{-0.5, -0.3}, {0, 0}}],
210 FontSize->Smaller], "T",
211 AdjustmentBox["E",
212 BoxBaselineShift->0.5,
213 BoxMargins->{{-0.3, 0}, {0, 0}}], "X"}],
214 SingleLetterItalics->False], TraditionalForm]]],
215 " Output"
216}], "Section",
217 CellChangeTimes->{{3.510516585335*^9, 3.510516590423654*^9}}],
218
219Cell[BoxData[
220 RowBox[{"WriteLaTeXOutput", "[", "]"}]], "Input",
221 CellChangeTimes->{{3.510516592624229*^9, 3.510516605419121*^9}}]
222}, Open ]],
223
224Cell[CellGroupData[{
225
226Cell["Generating MC model files with optimized restrictions", "Section",
227 CellChangeTimes->{{3.51051202409053*^9, 3.5105120509193287`*^9}}],
228
229Cell[CellGroupData[{
230
231Cell["Optimization", "Subsubsection",
232 CellChangeTimes->{{3.508650844785081*^9, 3.508650846244998*^9}}],
233
234Cell[BoxData[
235 RowBox[{
236 RowBox[{"WriteRestrictionFile", "[", "]"}], ";", " ",
237 RowBox[{"LoadRestriction", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";", " ",
238 RowBox[{"DeleteFile", "[", "\"\<ZeroValues.rst\>\"", "]"}], ";"}]], "Input",\
239
240 CellChangeTimes->{{3.5080579535629387`*^9, 3.508057961830206*^9}, {
241 3.508058109308621*^9, 3.508058122229617*^9}}]
242}, Open ]],
243
244Cell[CellGroupData[{
245
246Cell["UFO interface", "Subsubsection",
247 CellChangeTimes->{{3.508650849632307*^9, 3.508650853575732*^9}}],
248
249Cell[BoxData[{
250 RowBox[{
251 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
252 RowBox[{
253 RowBox[{"WriteUFO", "[",
254 RowBox[{"LagFeynman", ",", " ",
255 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
256
257 CellChangeTimes->{{3.5073089286898727`*^9, 3.507308937458375*^9}, {
258 3.507311910076086*^9, 3.507311939653351*^9}, {3.507398897853717*^9,
259 3.507398901626676*^9}, {3.507399273736373*^9, 3.507399274557005*^9}, {
260 3.510512090662549*^9, 3.510512091990087*^9}, {3.511707532536199*^9,
261 3.5117075657661657`*^9}, 3.511708094750585*^9, {3.511708971080892*^9,
262 3.511708974317479*^9}, {3.511757080212324*^9, 3.511757081138858*^9}, {
263 3.511759587557466*^9, 3.5117595883076143`*^9}, {3.511767591006217*^9,
264 3.511767591302368*^9}, {3.511861703077114*^9, 3.511861704343335*^9}}]
265}, Open ]],
266
267Cell[CellGroupData[{
268
269Cell["CalcHEP model files", "Subsubsection",
270 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}}],
271
272Cell[BoxData[{
273 RowBox[{
274 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
275 RowBox[{
276 RowBox[{"WriteCHOutput", "[",
277 RowBox[{
278 RowBox[{"lagr", "/.",
279 RowBox[{"{",
280 RowBox[{
281 RowBox[{
282 RowBox[{"LUDD", "[", "__", "]"}], "\[Rule]", "0"}], ",",
283 RowBox[{
284 RowBox[{"TUDD", "[", "__", "]"}], "\[Rule]", "0"}]}], "}"}]}], ",",
285 " ",
286 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}], ",",
287 RowBox[{"CHSimplify", "\[Rule]", "False"}]}], "]"}], ";"}]}], "Input",
288 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
289 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.511681238205276*^9,
290 3.511681251773086*^9}, {3.5116991968628283`*^9, 3.5116992038702507`*^9}}]
291}, Open ]],
292
293Cell[CellGroupData[{
294
295Cell["MadGraph model files", "Subsubsection",
296 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
297 3.510513948444434*^9, 3.510513949301979*^9}}],
298
299Cell[BoxData[{
300 RowBox[{
301 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
302 RowBox[{
303 RowBox[{"WriteMGOutput", "[",
304 RowBox[{
305 RowBox[{"lagr", "/.",
306 RowBox[{"{",
307 RowBox[{
308 RowBox[{
309 RowBox[{"LUDD", "[", "__", "]"}], "\[Rule]", "0"}], ",",
310 RowBox[{
311 RowBox[{"TUDD", "[", "__", "]"}], "\[Rule]", "0"}]}], "}"}]}], ",",
312 " ",
313 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
314
315 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
316 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
317 3.510513955605171*^9}, 3.511681261085073*^9, {3.511699210333851*^9,
318 3.511699214431435*^9}, 3.511767650012632*^9}]
319}, Open ]],
320
321Cell[CellGroupData[{
322
323Cell["Whizard model files", "Subsubsection",
324 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
325 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
326 3.510514891766358*^9}}],
327
328Cell[BoxData[{
329 RowBox[{
330 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
331 RowBox[{
332 RowBox[{"WriteWOOutput", "[",
333 RowBox[{
334 RowBox[{"lagr", "/.",
335 RowBox[{"{",
336 RowBox[{
337 RowBox[{
338 RowBox[{"LUDD", "[", "__", "]"}], "\[Rule]", "0"}], ",",
339 RowBox[{
340 RowBox[{"TUDD", "[", "__", "]"}], "\[Rule]", "0"}]}], "}"}]}], ",",
341 " ",
342 RowBox[{"Exclude4Scalars", "\[Rule]", "True"}]}], "]"}], ";"}]}], "Input",\
343
344 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
345 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
346 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
347 3.510572235119359*^9, 3.51057223551847*^9}, {3.510572377588068*^9,
348 3.510572378543983*^9}, {3.510572531362316*^9, 3.510572531820734*^9},
349 3.511681267132571*^9, {3.5116992155985603`*^9, 3.511699218031599*^9}, {
350 3.511767595410138*^9, 3.511767597874593*^9}, {3.511768742819017*^9,
351 3.51176874413035*^9}}]
352}, Open ]],
353
354Cell[CellGroupData[{
355
356Cell["FeynArts model files", "Subsubsection",
357 CellChangeTimes->{{3.510512057239231*^9, 3.510512065782227*^9}, {
358 3.510513948444434*^9, 3.510513949301979*^9}, {3.510514890712902*^9,
359 3.510514891766358*^9}, {3.5105162477375813`*^9, 3.5105162494145308`*^9}}],
360
361Cell[BoxData[{
362 RowBox[{
363 RowBox[{"SetDirectory", "[", "$OldDir", "]"}], ";"}], "\[IndentingNewLine]",
364 RowBox[{
365 RowBox[{"WriteFeynArtsOutput", "[", "lagr", "]"}], ";"}]}], "Input",
366 CellChangeTimes->{{3.5105120994316807`*^9, 3.510512105814425*^9}, {
367 3.5105127146047497`*^9, 3.510512717826407*^9}, {3.510513952348884*^9,
368 3.510513955605171*^9}, {3.510514894507822*^9, 3.510514894666054*^9}, {
369 3.5105162369063063`*^9, 3.510516243526369*^9}}]
370}, Open ]]
371}, Open ]]
372},
373WindowSize->{1914, 1123},
374WindowMargins->{{3, Automatic}, {Automatic, 23}},
375ShowSelection->True,
376FrontEndVersion->"7.0 for Linux x86 (64-bit) (February 25, 2009)",
377StyleDefinitions->"Default.nb"
378]
379(* End of Notebook Content *)
380
381(* Internal cache information *)
382(*CellTagsOutline
383CellTagsIndex->{}
384*)
385(*CellTagsIndex
386CellTagsIndex->{}
387*)
388(*NotebookFileOutline
389Notebook[{
390Cell[CellGroupData[{
391Cell[567, 22, 24, 0, 73, "Section"],
392Cell[594, 24, 881, 20, 143, "Input"]
393}, Open ]],
394Cell[CellGroupData[{
395Cell[1512, 49, 260, 3, 73, "Section"],
396Cell[CellGroupData[{
397Cell[1797, 56, 243, 4, 29, "Subsubsection"],
398Cell[2043, 62, 316, 9, 55, "Input"]
399}, Open ]],
400Cell[CellGroupData[{
401Cell[2396, 76, 140, 1, 29, "Subsubsection"],
402Cell[2539, 79, 732, 16, 121, "Input"]
403}, Open ]],
404Cell[CellGroupData[{
405Cell[3308, 100, 241, 5, 29, "Subsubsection"],
406Cell[3552, 107, 911, 22, 32, "Input"]
407}, Open ]],
408Cell[CellGroupData[{
409Cell[4500, 134, 360, 6, 29, "Subsubsection"],
410Cell[4863, 142, 459, 9, 55, "Input"]
411}, Open ]],
412Cell[CellGroupData[{
413Cell[5359, 156, 357, 6, 29, "Subsubsection"],
414Cell[5719, 164, 611, 12, 55, "Input"]
415}, Open ]]
416}, Open ]],
417Cell[CellGroupData[{
418Cell[6379, 182, 112, 1, 73, "Section"],
419Cell[6494, 185, 561, 10, 32, "Input"]
420}, Open ]],
421Cell[CellGroupData[{
422Cell[7092, 200, 461, 16, 75, "Section"],
423Cell[7556, 218, 130, 2, 32, "Input"]
424}, Open ]],
425Cell[CellGroupData[{
426Cell[7723, 225, 139, 1, 73, "Section"],
427Cell[CellGroupData[{
428Cell[7887, 230, 103, 1, 29, "Subsubsection"],
429Cell[7993, 233, 356, 7, 32, "Input"]
430}, Open ]],
431Cell[CellGroupData[{
432Cell[8386, 245, 104, 1, 29, "Subsubsection"],
433Cell[8493, 248, 841, 15, 55, "Input"]
434}, Open ]],
435Cell[CellGroupData[{
436Cell[9371, 268, 110, 1, 29, "Subsubsection"],
437Cell[9484, 271, 756, 18, 55, "Input"]
438}, Open ]],
439Cell[CellGroupData[{
440Cell[10277, 294, 160, 2, 29, "Subsubsection"],
441Cell[10440, 298, 750, 19, 55, "Input"]
442}, Open ]],
443Cell[CellGroupData[{
444Cell[11227, 322, 208, 3, 29, "Subsubsection"],
445Cell[11438, 327, 1020, 23, 55, "Input"]
446}, Open ]],
447Cell[CellGroupData[{
448Cell[12495, 355, 259, 3, 29, "Subsubsection"],
449Cell[12757, 360, 450, 8, 55, "Input"]
450}, Open ]]
451}, Open ]]
452}
453]
454*)
455
456(* End of internal cache information *)