RPVMSSM: RPV.nb

File RPV.nb, 15.1 KB (added by BenjF, 8 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 *)