HeavyN: heavyN.nb

File heavyN.nb, 7.6 KB (added by richardphysics, 20 months ago)

Type I Seesaw Model FR mathematica notebook that loads model file, performs basic sanity checks, and outputs to desired format, e.g., UFO.

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[      7586,        217]
14NotebookOptionsPosition[      6428,        175]
15NotebookOutlinePosition[      6866,        192]
16CellTagsIndexPosition[      6823,        189]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell["FeynRules Setup", "Title",
24 CellChangeTimes->{{3.596207285657744*^9, 3.596207286276866*^9}, {
25  3.596207330352247*^9, 3.59620733250506*^9}, {3.655035967648223*^9,
26  3.65503597269068*^9}}],
27
28Cell[BoxData[{
29 RowBox[{
30  RowBox[{"$CurrentPath", "=",
31   RowBox[{"Directory", "[", "]"}]}], ";"}], "\[IndentingNewLine]",
32 RowBox[{
33  RowBox[{"$FeynRulesPath", "=",
34   RowBox[{
35   "SetDirectory", "[", "\"\<~/Scripts/FeynRules/feynrules-2.0.x\>\"",
36    "]"}]}], ";"}], "\[IndentingNewLine]",
37 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
38 RowBox[{
39  RowBox[{"SetDirectory", "[", "$CurrentPath", "]"}], ";"}]}], "Input",
40 CellChangeTimes->{{3.41265862251538*^9, 3.412658649947229*^9}, {
41   3.423415585782702*^9, 3.423415597189939*^9}, {3.4234163173467493`*^9,
42   3.4234163227881193`*^9}, {3.572102427116931*^9, 3.572102456206238*^9}, {
43   3.572280529035656*^9, 3.572280529468644*^9}, {3.572366766720894*^9,
44   3.572366773370627*^9}, {3.596207212826934*^9, 3.5962072447909603`*^9}, {
45   3.596914045279478*^9, 3.5969140472383204`*^9}, 3.596914102616413*^9,
46   3.596915921742405*^9, {3.596915961175353*^9, 3.59691596259039*^9}, {
47   3.60744636894991*^9, 3.607446377507101*^9}, {3.6364757082312813`*^9,
48   3.636475708315627*^9}, {3.655035981353272*^9, 3.655035982244761*^9}, {
49   3.6550360677717857`*^9, 3.655036081259776*^9}, {3.655036117419112*^9,
50   3.655036143336364*^9}, {3.655037954145546*^9, 3.6550379644813023`*^9}}],
51
52Cell[BoxData[
53 RowBox[{"(*",
54  RowBox[{
55  "Must", " ", "load", " ", "default", " ", "SM", " ", "file", " ", "with",
56   " ", "Heavy", " ", "N", " ",
57   RowBox[{"extension", "."}]}], "*)"}]], "Input",
58 CellChangeTimes->{{3.6550381551475554`*^9, 3.655038155711425*^9}}],
59
60Cell[BoxData[
61 RowBox[{"LoadModel", "[",
62  RowBox[{"\"\<sm.fr\>\"", ",", "\"\<heavyN.fr\>\""}], "]"}]], "Input",
63 CellChangeTimes->{{3.6550379628078337`*^9, 3.6550379664979897`*^9}}],
64
65Cell[BoxData[
66 RowBox[{"LoadRestriction", "[",
67  RowBox[{"\"\<diagonalCKM.rst\>\"", ",", "\"\<massless.rst\>\""}],
68  "]"}]], "Input"],
69
70Cell[CellGroupData[{
71
72Cell["New Lagrangian Terms", "Section",
73 CellChangeTimes->{{3.411910065227421*^9, 3.411910071535137*^9}, {
74  3.572348211890565*^9, 3.572348234275403*^9}, {3.596207324453108*^9,
75  3.596207326711195*^9}, {3.655036212669145*^9, 3.6550362220422897`*^9}, {
76  3.655036282693001*^9, 3.655036287060457*^9}, {3.655037866940884*^9,
77  3.6550378926486187`*^9}}],
78
79Cell[BoxData[
80 RowBox[{
81  RowBox[{"FeynmanGauge", "=", "False"}], ";"}]], "Input",
82 CellChangeTimes->{{3.572348235565138*^9, 3.572348241872657*^9}}],
83
84Cell[BoxData[{"LNKin", "\[IndentingNewLine]", "LNCC", "\[IndentingNewLine]", \
85"LNNC", "\[IndentingNewLine]", "LNYukawa"}], "Input",
86 CellChangeTimes->{{3.5723482442775097`*^9, 3.572348252612296*^9}, {
87  3.596212546796714*^9, 3.596212547440735*^9}, {3.655037865036394*^9,
88  3.6550378653420067`*^9}}]
89}, Open  ]],
90
91Cell[CellGroupData[{
92
93Cell["Sanity Checks for New Lagrangian Terms", "Section",
94 CellChangeTimes->{{3.4234156496121893`*^9, 3.423415657027569*^9}, {
95  3.655037883884348*^9, 3.655037906883815*^9}}],
96
97Cell[BoxData[
98 RowBox[{"CheckHermiticity", "[", "LN", "]"}]], "Input",
99 CellChangeTimes->{3.607446407230465*^9}],
100
101Cell[BoxData[
102 RowBox[{"CheckMassSpectrum", "[", "LN", "]"}]], "Input",
103 CellChangeTimes->{{3.607446409988699*^9, 3.6074464110954323`*^9},
104   3.6550378392495413`*^9}],
105
106Cell[BoxData[
107 RowBox[{"CheckKineticTermNormalisation", "[", "LN", "]"}]], "Input",
108 CellChangeTimes->{{3.414763473002973*^9, 3.4147634845507174`*^9}, {
109   3.572348320182221*^9, 3.572348326613598*^9}, {3.572358912362153*^9,
110   3.572358917148378*^9}, {3.596207478919518*^9, 3.596207486580676*^9}, {
111   3.596212578276376*^9, 3.5962125890639143`*^9}, 3.607446414234477*^9}]
112}, Open  ]],
113
114Cell[CellGroupData[{
115
116Cell["\<\
117New Feynman Rules for Heavy N Interactions\
118\>", "Section",
119 CellChangeTimes->{{3.411910637332127*^9, 3.41191064360071*^9}, {
120  3.655037529059973*^9, 3.655037541428491*^9}, {3.6550379116719217`*^9,
121  3.655037913954034*^9}}],
122
123Cell[BoxData[
124 RowBox[{"(*", " ",
125  RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
126   RowBox[{"appear", ":", " ",
127    RowBox[{"\"\<\!\(\*
128StyleBox[\"QN\", \"MessageName\"]\)\!\(\*
129StyleBox[\"::\", \"MessageName\"]\)\!\(\*
130StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"", "Warning"}], ":", " ",
131    RowBox[{
132    "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
133     " ",
134     RowBox[{"encountered", "!"}],
135     "\"\<;\nThis is expected as Lepton Number is violated and FeynRules \
136checks for this.\>"}]}]}], "*)"}]], "Input",
137 CellChangeTimes->{{3.6550374096249533`*^9, 3.655037447663039*^9},
138   3.655037551040966*^9, {3.656940976048616*^9, 3.6569409786068287`*^9}}],
139
140Cell[BoxData[
141 RowBox[{"FeynmanRules", "[", "LN", "]"}]], "Input"]
142}, Open  ]],
143
144Cell[CellGroupData[{
145
146Cell["UFO Output", "Section",
147 CellChangeTimes->{
148  3.411910944409371*^9, {3.412185514054689*^9, 3.412185517014236*^9}, {
149   3.59620768137986*^9, 3.596207684212212*^9}}],
150
151Cell[BoxData[
152 RowBox[{"(*", " ",
153  RowBox[{"The", " ", "following", " ", "warnings", " ", "will", " ",
154   RowBox[{"appear", ":", " ",
155    RowBox[{"\"\<\!\(\*
156StyleBox[\"QN\", \"MessageName\"]\)\!\(\*
157StyleBox[\"::\", \"MessageName\"]\)\!\(\*
158StyleBox[\"NonConserv\", \"MessageName\"]\): \>\"", "Warning"}], ":", " ",
159    RowBox[{
160    "non", " ", "quantum", " ", "number", " ", "conserving", " ", "vertex",
161     " ",
162     RowBox[{"encountered", "!"}],
163     "\"\<;\nThis is expected as Lepton Number is violated and FeynRules \
164checks for this.\>"}]}]}], "*)"}]], "Input",
165 CellChangeTimes->{{3.656940935403554*^9, 3.656940952617087*^9}}],
166
167Cell[BoxData[
168 RowBox[{
169  RowBox[{"WriteUFO", "[",
170   RowBox[{
171   "LGauge", ",", "LHiggs", ",", "LFermions", ",", "LYukawa", ",", "LGhost",
172    ",", "LN"}], "]"}], ";"}]], "Input"]
173}, Open  ]]
174}, Open  ]]
175},
176WindowSize->{1196, 628},
177WindowMargins->{{192, Automatic}, {Automatic, 64}},
178PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
179ShowSelection->True,
180FrontEndVersion->"9.0 for Linux x86 (64-bit) (November 20, 2012)",
181StyleDefinitions->"Default.nb"
182]
183(* End of Notebook Content *)
184
185(* Internal cache information *)
186(*CellTagsOutline
187CellTagsIndex->{}
188*)
189(*CellTagsIndex
190CellTagsIndex->{}
191*)
192(*NotebookFileOutline
193Notebook[{
194Cell[CellGroupData[{
195Cell[567, 22, 194, 3, 93, "Title"],
196Cell[764, 27, 1239, 22, 99, "Input"],
197Cell[2006, 51, 269, 6, 32, "Input"],
198Cell[2278, 59, 183, 3, 32, "Input"],
199Cell[2464, 64, 135, 3, 32, "Input"],
200Cell[CellGroupData[{
201Cell[2624, 71, 351, 5, 81, "Section"],
202Cell[2978, 78, 148, 3, 32, "Input"],
203Cell[3129, 83, 300, 4, 99, "Input"]
204}, Open  ]],
205Cell[CellGroupData[{
206Cell[3466, 92, 174, 2, 81, "Section"],
207Cell[3643, 96, 112, 2, 32, "Input"],
208Cell[3758, 100, 167, 3, 32, "Input"],
209Cell[3928, 105, 371, 5, 32, "Input"]
210}, Open  ]],
211Cell[CellGroupData[{
212Cell[4336, 115, 234, 5, 81, "Section"],
213Cell[4573, 122, 720, 15, 77, "Input"],
214Cell[5296, 139, 66, 1, 32, "Input"]
215}, Open  ]],
216Cell[CellGroupData[{
217Cell[5399, 145, 169, 3, 81, "Section"],
218Cell[5571, 150, 644, 14, 77, "Input"],
219Cell[6218, 166, 182, 5, 32, "Input"]
220}, Open  ]]
221}, Open  ]]
222}
223]
224*)
225
226(* End of internal cache information *)