ASperGe: 2HDM.nb

File 2HDM.nb, 31.1 KB (added by Karen De Causmaecker, 12 years ago)
Line 
1(* Content-type: application/vnd.wolfram.mathematica *)
2
3(*** Wolfram Notebook File ***)
4(* http://www.wolfram.com/nb *)
5
6(* CreatedBy='Mathematica 8.0' *)
7
8(*CacheID: 234*)
9(* Internal cache information:
10NotebookFileLineBreakTest
11NotebookFileLineBreakTest
12NotebookDataPosition[ 157, 7]
13NotebookDataLength[ 31607, 900]
14NotebookOptionsPosition[ 28108, 773]
15NotebookOutlinePosition[ 28462, 789]
16CellTagsIndexPosition[ 28419, 786]
17WindowFrame->Normal*)
18
19(* Beginning of Notebook Content *)
20Notebook[{
21
22Cell[CellGroupData[{
23Cell["\<\
24ASperGe:
25The Two-Higgs-Doublet Model\
26\>", "Title",
27 CellChangeTimes->{{3.5678552816101*^9, 3.5678553411857243`*^9}, {
28 3.567855585658388*^9, 3.567855586280493*^9}, {3.568094017817658*^9,
29 3.568094022181663*^9}}],
30
31Cell[BoxData[
32 RowBox[{"Quit", "[", "]"}]], "Input",
33 CellChangeTimes->{{3.567855348415655*^9, 3.567855350834606*^9}}],
34
35Cell[CellGroupData[{
36
37Cell["Set the path", "Section",
38 CellChangeTimes->{{3.56785531856068*^9, 3.567855344393437*^9}}],
39
40Cell[CellGroupData[{
41
42Cell[BoxData[{
43 RowBox[{
44 RowBox[{"$FeynRulesPath", " ", "=", " ",
45 RowBox[{
46 "SetDirectory", "[", "\"\<<< address of the package >>\>\"", "]"}]}],
47 ";"}], "\[IndentingNewLine]",
48 RowBox[{"<<", "FeynRules`"}], "\[IndentingNewLine]",
49 RowBox[{
50 RowBox[{
51 "SetDirectory", "[", "\"\<<< address of the model file >>\>\"", "]"}],
52 ";"}]}], "Input",
53 CellChangeTimes->{{3.567862017826684*^9, 3.567862020512813*^9}, {
54 3.567862072745294*^9, 3.567862099825725*^9}, {3.568293595866456*^9,
55 3.568293621441877*^9}}],
56
57Cell[CellGroupData[{
58
59Cell[BoxData["\<\" - FeynRules - \"\>"], "Print",
60 CellChangeTimes->{
61 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
62 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
63 3.568096511733856*^9, 3.5682935538776083`*^9}],
64
65Cell[BoxData[
66 InterpretationBox[
67 RowBox[{"\<\"Version: \"\>", "\[InvisibleSpace]", "\<\"1.7.134\"\>",
68 "\[InvisibleSpace]",
69 RowBox[{"\<\" (\"\>", " ", "\<\"27 January 2013\"\>"}],
70 "\[InvisibleSpace]", "\<\").\"\>"}],
71 SequenceForm["Version: ", "1.7.134", " (" "27 January 2013", ")."],
72 Editable->False]], "Print",
73 CellChangeTimes->{
74 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
75 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
76 3.568096511733856*^9, 3.5682935538810253`*^9}],
77
78Cell[BoxData["\<\"Authors: A. Alloul, N. Christensen, C. Degrande, C. Duhr, \
79B. Fuks\"\>"], "Print",
80 CellChangeTimes->{
81 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
82 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
83 3.568096511733856*^9, 3.5682935538841057`*^9}],
84
85Cell[BoxData["\<\" \"\>"], "Print",
86 CellChangeTimes->{
87 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
88 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
89 3.568096511733856*^9, 3.5682935538873873`*^9}],
90
91Cell[BoxData["\<\"Please cite: Comput.Phys.Commun.180:1614-1641,2009 \
92(arXiv:0806.4194).\"\>"], "Print",
93 CellChangeTimes->{
94 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
95 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
96 3.568096511733856*^9, 3.568293553891555*^9}],
97
98Cell[BoxData["\<\"http://feynrules.phys.ucl.ac.be\"\>"], "Print",
99 CellChangeTimes->{
100 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
101 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
102 3.568096511733856*^9, 3.568293553895336*^9}],
103
104Cell[BoxData["\<\" \"\>"], "Print",
105 CellChangeTimes->{
106 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
107 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
108 3.568096511733856*^9, 3.568293553898706*^9}],
109
110Cell[BoxData["\<\"The FeynRules palette can be opened using the command \
111FRPalette[].\"\>"], "Print",
112 CellChangeTimes->{
113 3.567862106393355*^9, 3.568032097277495*^9, 3.568093569924898*^9, {
114 3.5680938564749527`*^9, 3.5680938787686234`*^9}, 3.568095912005393*^9,
115 3.568096511733856*^9, 3.568293553901369*^9}]
116}, Open ]]
117}, Open ]]
118}, Open ]],
119
120Cell[CellGroupData[{
121
122Cell["Import the model", "Section",
123 CellChangeTimes->{{3.567855368640649*^9, 3.5678553714967747`*^9}}],
124
125Cell[CellGroupData[{
126
127Cell[BoxData[{
128 RowBox[{
129 RowBox[{"LoadModel", "[", "\"\<2HDM_mix.fr\>\"", "]"}],
130 ";"}], "\[IndentingNewLine]",
131 RowBox[{
132 RowBox[{"SetDirectory", "[", "$FeynRulesPath", "]"}], ";"}]}], "Input"],
133
134Cell[CellGroupData[{
135
136Cell[BoxData["\<\"This model implementation was created by\"\>"], "Print",
137 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
138 3.568096517242443*^9, 3.568293561865912*^9}],
139
140Cell[BoxData["\<\"K. De Causmaecker\"\>"], "Print",
141 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
142 3.568096517242443*^9, 3.568293561869576*^9}],
143
144Cell[BoxData[
145 InterpretationBox[
146 RowBox[{"\<\"Model Version: \"\>", "\[InvisibleSpace]", "\<\"1.0.0\"\>"}],
147 SequenceForm["Model Version: ", "1.0.0"],
148 Editable->False]], "Print",
149 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
150 3.568096517242443*^9, 3.56829356187389*^9}],
151
152Cell[BoxData["\<\"For more information, type ModelInformation[].\"\>"], \
153"Print",
154 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
155 3.568096517242443*^9, 3.568293561876544*^9}],
156
157Cell[BoxData["\<\"\"\>"], "Print",
158 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
159 3.568096517242443*^9, 3.5682935618804493`*^9}],
160
161Cell[BoxData["\<\" - Loading particle classes.\"\>"], "Print",
162 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
163 3.568096517242443*^9, 3.568293561885044*^9}],
164
165Cell[BoxData["\<\" - Loading gauge group classes.\"\>"], "Print",
166 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
167 3.568096517242443*^9, 3.5682935620747223`*^9}],
168
169Cell[BoxData["\<\" - Loading parameter classes.\"\>"], "Print",
170 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
171 3.568096517242443*^9, 3.5682935620996113`*^9}],
172
173Cell[BoxData[
174 InterpretationBox[
175 RowBox[{"\<\"\\nModel \"\>", "\[InvisibleSpace]", "\<\"2HDM_mix\"\>",
176 "\[InvisibleSpace]", "\<\" loaded.\"\>"}],
177 SequenceForm["\nModel ", "2HDM_mix", " loaded."],
178 Editable->False]], "Print",
179 CellChangeTimes->{3.568093884068632*^9, 3.568095917489497*^9,
180 3.568096517242443*^9, 3.568293562174679*^9}]
181}, Open ]]
182}, Open ]]
183}, Open ]],
184
185Cell[CellGroupData[{
186
187Cell["Compute the tree level mass matrices", "Section",
188 CellChangeTimes->{{3.567855382889579*^9, 3.567855388328291*^9}}],
189
190Cell[CellGroupData[{
191
192Cell[TextData[StyleBox["Compute all the mass matrices", "Subsection"]], \
193"Subsubsection",
194 CellChangeTimes->{{3.567855393721282*^9, 3.567855397792397*^9}}],
195
196Cell[CellGroupData[{
197
198Cell[BoxData[
199 RowBox[{"ComputeTreeLevelMassMatrix", "[", "L2HDM", "]"}]], "Input",
200 CellChangeTimes->{{3.568035094618717*^9, 3.568035125565981*^9}}],
201
202Cell[CellGroupData[{
203
204Cell[BoxData["\<\"Index expansion ...\"\>"], "Print",
205 CellChangeTimes->{{3.568035102003182*^9, 3.568035110486939*^9},
206 3.568093596316161*^9, {3.568093888110133*^9, 3.568093906582636*^9},
207 3.568093942054863*^9, 3.56809592057561*^9, 3.568096519641309*^9,
208 3.568293564042506*^9}],
209
210Cell[BoxData[
211 InterpretationBox[
212 RowBox[{"\<\" => done in \"\>", "\[InvisibleSpace]", "0.627605`",
213 "\[InvisibleSpace]", "\<\" seconds.\"\>"}],
214 SequenceForm[" => done in ", 0.627605, " seconds."],
215 Editable->False]], "Print",
216 CellChangeTimes->{{3.568035102003182*^9, 3.568035110486939*^9},
217 3.568093596316161*^9, {3.568093888110133*^9, 3.568093906582636*^9},
218 3.568093942054863*^9, 3.56809592057561*^9, 3.568096519641309*^9,
219 3.5682935646694183`*^9}],
220
221Cell[BoxData["\<\"Tree level mass matrix computations achieved. Only the mass \
222matrices necessary for the C++ code have been computed.\"\>"], "Print",
223 CellChangeTimes->{{3.568035102003182*^9, 3.568035110486939*^9},
224 3.568093596316161*^9, {3.568093888110133*^9, 3.568093906582636*^9},
225 3.568093942054863*^9, 3.56809592057561*^9, 3.568096519641309*^9,
226 3.568293565074601*^9}],
227
228Cell[BoxData["\<\"Check FR$MassMatrices for the list of the MixingIDs.\"\>"], \
229"Print",
230 CellChangeTimes->{{3.568035102003182*^9, 3.568035110486939*^9},
231 3.568093596316161*^9, {3.568093888110133*^9, 3.568093906582636*^9},
232 3.568093942054863*^9, 3.56809592057561*^9, 3.568096519641309*^9,
233 3.568293565076906*^9}],
234
235Cell[BoxData["\<\"Use MassMatrix[ MixingID ], GaugeBasis[ MixingID ], \
236MassBasis[ MixingID ], BlockName[ MixingID ] and MatrixSymbol[ MixingID ] for \
237more information on the results.\"\>"], "Print",
238 CellChangeTimes->{{3.568035102003182*^9, 3.568035110486939*^9},
239 3.568093596316161*^9, {3.568093888110133*^9, 3.568093906582636*^9},
240 3.568093942054863*^9, 3.56809592057561*^9, 3.568096519641309*^9,
241 3.5682935650794163`*^9}]
242}, Open ]]
243}, Open ]]
244}, Open ]],
245
246Cell[CellGroupData[{
247
248Cell["\<\
249Compute all the mass matrices without generation of the output about the mass \
250matrices\
251\>", "Subsection",
252 CellChangeTimes->{{3.5680352005369177`*^9, 3.568035215142868*^9}, {
253 3.56809361626239*^9, 3.568093620750354*^9}}],
254
255Cell[BoxData[
256 RowBox[{"ComputeTreeLevelMassMatrix", "[",
257 RowBox[{"L2HDM", ",", " ",
258 RowBox[{"ScreenOutput", "\[Rule]", " ", "False"}]}], "]"}]], "Input",
259 CellChangeTimes->{{3.568035137500168*^9, 3.568035158886353*^9}}],
260
261Cell[CellGroupData[{
262
263Cell[TextData[StyleBox["Compute only the mass matrix of mixing of the neutral \
264Higgses (1s)", "Subsection"]], "Subsubsection",
265 CellChangeTimes->{{3.5678554211273003`*^9, 3.567855435470812*^9}, {
266 3.568093674230916*^9, 3.568093681100155*^9}, {3.5680937152983913`*^9,
267 3.568093723634061*^9}}],
268
269Cell[CellGroupData[{
270
271Cell[BoxData[
272 RowBox[{"ComputeTreeLevelMassMatrix", "[",
273 RowBox[{"L2HDM", ",", " ",
274 RowBox[{"Mix", "\[Rule]", "\"\<1s\>\""}]}], "]"}]], "Input",
275 CellChangeTimes->{{3.568032145034424*^9, 3.568032146074194*^9}, {
276 3.568093642312999*^9, 3.5680936439442167`*^9}}],
277
278Cell[CellGroupData[{
279
280Cell[BoxData["\<\"Index expansion ...\"\>"], "Print",
281 CellChangeTimes->{3.568032146920103*^9, 3.568093699459681*^9,
282 3.5680939136317167`*^9, 3.5680939682239017`*^9, 3.568095927447856*^9,
283 3.568096567491313*^9, 3.568293571684527*^9}],
284
285Cell[BoxData[
286 InterpretationBox[
287 RowBox[{"\<\" => done in \"\>", "\[InvisibleSpace]", "0.6513780000000002`",
288 "\[InvisibleSpace]", "\<\" seconds.\"\>"}],
289 SequenceForm[" => done in ", 0.6513780000000002, " seconds."],
290 Editable->False]], "Print",
291 CellChangeTimes->{3.568032146920103*^9, 3.568093699459681*^9,
292 3.5680939136317167`*^9, 3.5680939682239017`*^9, 3.568095927447856*^9,
293 3.568096567491313*^9, 3.5682935723352137`*^9}],
294
295Cell[BoxData["\<\"Computation of the tree level mass matrix for the mixing id \
2961s achieved.\"\>"], "Print",
297 CellChangeTimes->{3.568032146920103*^9, 3.568093699459681*^9,
298 3.5680939136317167`*^9, 3.5680939682239017`*^9, 3.568095927447856*^9,
299 3.568096567491313*^9, 3.568293572505595*^9}],
300
301Cell[BoxData["\<\"For the results, use MassMatrix[ 1s ], GaugeBasis[ 1s ], \
302MassBasis[ 1s ], BlockName [ 1s ] and MatrixSymbol [ 1s ].\"\>"], "Print",
303 CellChangeTimes->{3.568032146920103*^9, 3.568093699459681*^9,
304 3.5680939136317167`*^9, 3.5680939682239017`*^9, 3.568095927447856*^9,
305 3.568096567491313*^9, 3.5682935725091248`*^9}],
306
307Cell[BoxData[
308 TagBox[
309 RowBox[{"(", "\[NoBreak]", GridBox[{
310 {
311 RowBox[{"mu1", "+",
312 RowBox[{"3", " ", "lambda1", " ",
313 SuperscriptBox["vev1", "2"]}]}],
314 RowBox[{
315 FractionBox["1", "4"], " ",
316 RowBox[{"(",
317 RowBox[{
318 RowBox[{"2", " ", "mu3"}], "+",
319 RowBox[{"3", " ", "lambda6", " ",
320 SuperscriptBox["vev1", "2"]}], "+",
321 RowBox[{"3", " ",
322 SuperscriptBox["vev1", "2"], " ",
323 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
324 RowBox[{"2", " ",
325 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
326 RowBox[{
327 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
328 RowBox[{"(",
329 RowBox[{
330 RowBox[{"2", " ", "mu3"}], "+",
331 RowBox[{"3", " ", "lambda6", " ",
332 SuperscriptBox["vev1", "2"]}], "-",
333 RowBox[{"3", " ",
334 SuperscriptBox["vev1", "2"], " ",
335 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
336 RowBox[{"2", " ",
337 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}]},
338 {
339 RowBox[{
340 FractionBox["1", "4"], " ",
341 RowBox[{"(",
342 RowBox[{
343 RowBox[{"2", " ", "mu3"}], "+",
344 RowBox[{"3", " ", "lambda6", " ",
345 SuperscriptBox["vev1", "2"]}], "+",
346 RowBox[{"3", " ",
347 SuperscriptBox["vev1", "2"], " ",
348 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
349 RowBox[{"2", " ",
350 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
351 RowBox[{"mu2", "+",
352 RowBox[{
353 FractionBox["1", "2"], " ",
354 RowBox[{"(",
355 RowBox[{"lambda3", "+", "lambda4", "+",
356 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
357 SuperscriptBox["vev1", "2"]}]}], "0"},
358 {
359 RowBox[{
360 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
361 RowBox[{"(",
362 RowBox[{
363 RowBox[{"2", " ", "mu3"}], "+",
364 RowBox[{"3", " ", "lambda6", " ",
365 SuperscriptBox["vev1", "2"]}], "-",
366 RowBox[{"3", " ",
367 SuperscriptBox["vev1", "2"], " ",
368 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
369 RowBox[{"2", " ",
370 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}], "0",
371 RowBox[{
372 FractionBox["1", "2"], " ",
373 RowBox[{"(",
374 RowBox[{
375 RowBox[{"2", " ", "mu2"}], "+",
376 RowBox[{
377 RowBox[{"(",
378 RowBox[{"lambda3", "+", "lambda4", "-",
379 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
380 SuperscriptBox["vev1", "2"]}]}], ")"}]}]}
381 },
382 GridBoxAlignment->{
383 "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
384 "RowsIndexed" -> {}},
385 GridBoxSpacings->{"Columns" -> {
386 Offset[0.27999999999999997`], {
387 Offset[0.7]},
388 Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
389 Offset[0.2], {
390 Offset[0.4]},
391 Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
392 Function[BoxForm`e$,
393 MatrixForm[BoxForm`e$]]]], "Print",
394 CellChangeTimes->{3.568032146920103*^9, 3.568093699459681*^9,
395 3.5680939136317167`*^9, 3.5680939682239017`*^9, 3.568095927447856*^9,
396 3.568096567491313*^9, 3.5682935725637627`*^9}]
397}, Open ]]
398}, Open ]]
399}, Open ]]
400}, Open ]],
401
402Cell[CellGroupData[{
403
404Cell["\<\
405Ask for all the information about the mixing of the neutral Higgses (1s)\
406\>", "Subsection",
407 CellChangeTimes->{{3.56785545031773*^9, 3.567855460702108*^9}, {
408 3.568093685467507*^9, 3.568093729338764*^9}}],
409
410Cell[CellGroupData[{
411
412Cell[BoxData[
413 RowBox[{"MixingSummary", "[", "\"\<1s\>\"", "]"}]], "Input"],
414
415Cell[CellGroupData[{
416
417Cell[BoxData[
418 InterpretationBox[
419 RowBox[{"\<\"Gauge basis = \"\>", "\[InvisibleSpace]",
420 RowBox[{"{",
421 RowBox[{"H0", ",", "R0", ",", "I0"}], "}"}]}],
422 SequenceForm[
423 "Gauge basis = ", {$CellContext`H0, $CellContext`R0, $CellContext`I0}],
424 Editable->False]], "Print",
425 CellChangeTimes->{
426 3.5680937028476753`*^9, {3.568093956653385*^9, 3.568093971205803*^9},
427 3.568095930253375*^9, {3.568096552576515*^9, 3.5680965701244593`*^9},
428 3.568293574230373*^9}],
429
430Cell[BoxData[
431 InterpretationBox[
432 RowBox[{"\<\"Mass basis = \"\>", "\[InvisibleSpace]",
433 RowBox[{"{",
434 RowBox[{"h1", ",", "h2", ",", "h3"}], "}"}],
435 "\[InvisibleSpace]", "\<\"\\nBlock = \"\>", "\[InvisibleSpace]", "NHMIX",
436 "\[InvisibleSpace]", "\<\"\\nSymbol = \"\>", "\[InvisibleSpace]", "NH"}],
437 SequenceForm[
438 "Mass basis = ", {$CellContext`h1, $CellContext`h2, $CellContext`h3},
439 "\nBlock = ", $CellContext`NHMIX, "\nSymbol = ", $CellContext`NH],
440 Editable->False]], "Print",
441 CellChangeTimes->{
442 3.5680937028476753`*^9, {3.568093956653385*^9, 3.568093971205803*^9},
443 3.568095930253375*^9, {3.568096552576515*^9, 3.5680965701244593`*^9},
444 3.5682935742319736`*^9}],
445
446Cell[BoxData[
447 InterpretationBox[
448 RowBox[{"\<\"Squared mass matrix = \"\>", "\[InvisibleSpace]",
449 TagBox[
450 RowBox[{"(", "\[NoBreak]", GridBox[{
451 {
452 RowBox[{"mu1", "+",
453 RowBox[{"3", " ", "lambda1", " ",
454 SuperscriptBox["vev1", "2"]}]}],
455 RowBox[{
456 FractionBox["1", "4"], " ",
457 RowBox[{"(",
458 RowBox[{
459 RowBox[{"2", " ", "mu3"}], "+",
460 RowBox[{"3", " ", "lambda6", " ",
461 SuperscriptBox["vev1", "2"]}], "+",
462 RowBox[{"3", " ",
463 SuperscriptBox["vev1", "2"], " ",
464 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
465 RowBox[{"2", " ",
466 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
467 RowBox[{
468 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
469 RowBox[{"(",
470 RowBox[{
471 RowBox[{"2", " ", "mu3"}], "+",
472 RowBox[{"3", " ", "lambda6", " ",
473 SuperscriptBox["vev1", "2"]}], "-",
474 RowBox[{"3", " ",
475 SuperscriptBox["vev1", "2"], " ",
476 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
477 RowBox[{"2", " ",
478 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}]},
479 {
480 RowBox[{
481 FractionBox["1", "4"], " ",
482 RowBox[{"(",
483 RowBox[{
484 RowBox[{"2", " ", "mu3"}], "+",
485 RowBox[{"3", " ", "lambda6", " ",
486 SuperscriptBox["vev1", "2"]}], "+",
487 RowBox[{"3", " ",
488 SuperscriptBox["vev1", "2"], " ",
489 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
490 RowBox[{"2", " ",
491 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
492 RowBox[{"mu2", "+",
493 RowBox[{
494 FractionBox["1", "2"], " ",
495 RowBox[{"(",
496 RowBox[{"lambda3", "+", "lambda4", "+",
497 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
498 SuperscriptBox["vev1", "2"]}]}], "0"},
499 {
500 RowBox[{
501 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
502 RowBox[{"(",
503 RowBox[{
504 RowBox[{"2", " ", "mu3"}], "+",
505 RowBox[{"3", " ", "lambda6", " ",
506 SuperscriptBox["vev1", "2"]}], "-",
507 RowBox[{"3", " ",
508 SuperscriptBox["vev1", "2"], " ",
509 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
510 RowBox[{"2", " ",
511 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}], "0",
512 RowBox[{
513 FractionBox["1", "2"], " ",
514 RowBox[{"(",
515 RowBox[{
516 RowBox[{"2", " ", "mu2"}], "+",
517 RowBox[{
518 RowBox[{"(",
519 RowBox[{"lambda3", "+", "lambda4", "-",
520 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
521 SuperscriptBox["vev1", "2"]}]}], ")"}]}]}
522 },
523 GridBoxAlignment->{
524 "Columns" -> {{Center}}, "ColumnsIndexed" -> {},
525 "Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
526 GridBoxSpacings->{"Columns" -> {
527 Offset[0.27999999999999997`], {
528 Offset[0.7]},
529 Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
530 Offset[0.2], {
531 Offset[0.4]},
532 Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
533 Function[BoxForm`e$,
534 MatrixForm[BoxForm`e$]]]}],
535 SequenceForm["Squared mass matrix = ",
536 MatrixForm[{{
537 FeynRules`mu1 + 3 $CellContext`lambda1 $CellContext`vev1^2,
538 Rational[1, 4] (2 $CellContext`mu3 +
539 3 $CellContext`lambda6 $CellContext`vev1^2 +
540 3 $CellContext`vev1^2 Conjugate[$CellContext`lambda6] +
541 2 Conjugate[$CellContext`mu3]), Complex[0,
542 Rational[1, 4]] (2 $CellContext`mu3 +
543 3 $CellContext`lambda6 $CellContext`vev1^2 - 3 $CellContext`vev1^2
544 Conjugate[$CellContext`lambda6] - 2 Conjugate[$CellContext`mu3])}, {
545 Rational[1, 4] (2 $CellContext`mu3 +
546 3 $CellContext`lambda6 $CellContext`vev1^2 +
547 3 $CellContext`vev1^2 Conjugate[$CellContext`lambda6] +
548 2 Conjugate[$CellContext`mu3]), FeynRules`mu2 +
549 Rational[1, 2] ($CellContext`lambda3 + $CellContext`lambda4 +
550 2 $CellContext`lambda5) $CellContext`vev1^2, 0}, {Complex[0,
551 Rational[1, 4]] (2 $CellContext`mu3 +
552 3 $CellContext`lambda6 $CellContext`vev1^2 - 3 $CellContext`vev1^2
553 Conjugate[$CellContext`lambda6] - 2 Conjugate[$CellContext`mu3]), 0,
554 Rational[1, 2] (
555 2 FeynRules`mu2 + ($CellContext`lambda3 + $CellContext`lambda4 -
556 2 $CellContext`lambda5) $CellContext`vev1^2)}}]],
557 Editable->False]], "Print",
558 CellChangeTimes->{
559 3.5680937028476753`*^9, {3.568093956653385*^9, 3.568093971205803*^9},
560 3.568095930253375*^9, {3.568096552576515*^9, 3.5680965701244593`*^9},
561 3.568293574234605*^9}],
562
563Cell[BoxData[
564 InterpretationBox[
565 RowBox[{"\<\"The mixing matrix = \"\>", "\[InvisibleSpace]",
566 TagBox["\<\"Please use the numerical code.\"\>",
567 Function[BoxForm`e$,
568 MatrixForm[BoxForm`e$]]]}],
569 SequenceForm["The mixing matrix = ",
570 MatrixForm["Please use the numerical code."]],
571 Editable->False]], "Print",
572 CellChangeTimes->{
573 3.5680937028476753`*^9, {3.568093956653385*^9, 3.568093971205803*^9},
574 3.568095930253375*^9, {3.568096552576515*^9, 3.5680965701244593`*^9},
575 3.568293574237418*^9}]
576}, Open ]]
577}, Open ]]
578}, Open ]]
579}, Open ]],
580
581Cell[CellGroupData[{
582
583Cell["Write the ASperGE package", "Section",
584 CellChangeTimes->{{3.567855490212146*^9, 3.5678554958361464`*^9}}],
585
586Cell[CellGroupData[{
587
588Cell["\<\
589Write the output for all mass matrices to the folder \
590\[OpenCurlyQuote]2HDM_mix_MD\[CloseCurlyQuote] in $FeynRulesPath\
591\>", "Subsection",
592 CellChangeTimes->{{3.568035248605268*^9, 3.5680352719368563`*^9}, {
593 3.568035342262802*^9, 3.568035362863311*^9}}],
594
595Cell[CellGroupData[{
596
597Cell[BoxData[
598 RowBox[{"WriteASperGe", "[", "L2HDM", "]"}]], "Input",
599 CellChangeTimes->{{3.568035279856433*^9, 3.568035287557234*^9}}],
600
601Cell[CellGroupData[{
602
603Cell[BoxData["\<\"Parameters.cpp and Parameters.hpp output is finished.\\n \
604They have been stored in /Users/Karen/trunk/feynrules-development/2HDM_mix_MD\
605\"\>"], "Print",
606 CellChangeTimes->{3.5682935784949207`*^9}],
607
608Cell[BoxData["\<\"Checking parameter list...\"\>"], "Print",
609 CellChangeTimes->{3.5682935785269203`*^9}],
610
611Cell[BoxData["\<\"Warning: Ghost fields will be ignored.\"\>"], "Print",
612 CellChangeTimes->{3.568293578529491*^9}],
613
614Cell[BoxData[
615 InterpretationBox[
616 RowBox[{"\<\"Output written on \"\>",
617 "\[InvisibleSpace]", "\<\"input/externals.dat\"\>"}],
618 SequenceForm["Output written on ", "input/externals.dat"],
619 Editable->False]], "Print",
620 CellChangeTimes->{3.568293578531893*^9}]
621}, Open ]]
622}, Open ]]
623}, Open ]],
624
625Cell[CellGroupData[{
626
627Cell["\<\
628Write the output to the folder \[OpenCurlyQuote]2HDM\[CloseCurlyQuote] in \
629$FeynRulesPath only for the mixing \[OpenCurlyQuote]1s\[CloseCurlyQuote]
630For mutliple mixings use e.g. Mix -> {\[OpenCurlyDoubleQuote]1s\
631\[CloseCurlyDoubleQuote], \[OpenCurlyDoubleQuote]2s\[CloseCurlyDoubleQuote]}\
632\>", "Subsection",
633 CellChangeTimes->{{3.568035291371396*^9, 3.56803531000459*^9}, {
634 3.5680353661825333`*^9, 3.568035413797114*^9}, {3.5680934954740963`*^9,
635 3.568093532566071*^9}, {3.568096444862878*^9, 3.568096445036433*^9}}],
636
637Cell[CellGroupData[{
638
639Cell[BoxData[
640 RowBox[{"WriteASperGe", "[",
641 RowBox[{"L2HDM", ",", " ",
642 RowBox[{"Output", "\[Rule]", "\"\<2HDM\>\""}], ",", " ",
643 RowBox[{"Mix", "\[Rule]", "\"\<1s\>\""}]}], "]"}]], "Input",
644 CellChangeTimes->{{3.568035418086753*^9, 3.5680354407144623`*^9}}],
645
646Cell[CellGroupData[{
647
648Cell[BoxData[
649 TagBox[
650 RowBox[{"(", "\[NoBreak]", GridBox[{
651 {
652 RowBox[{"mu1", "+",
653 RowBox[{"3", " ", "lambda1", " ",
654 SuperscriptBox["vev1", "2"]}]}],
655 RowBox[{
656 FractionBox["1", "4"], " ",
657 RowBox[{"(",
658 RowBox[{
659 RowBox[{"2", " ", "mu3"}], "+",
660 RowBox[{"3", " ", "lambda6", " ",
661 SuperscriptBox["vev1", "2"]}], "+",
662 RowBox[{"3", " ",
663 SuperscriptBox["vev1", "2"], " ",
664 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
665 RowBox[{"2", " ",
666 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
667 RowBox[{
668 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
669 RowBox[{"(",
670 RowBox[{
671 RowBox[{"2", " ", "mu3"}], "+",
672 RowBox[{"3", " ", "lambda6", " ",
673 SuperscriptBox["vev1", "2"]}], "-",
674 RowBox[{"3", " ",
675 SuperscriptBox["vev1", "2"], " ",
676 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
677 RowBox[{"2", " ",
678 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}]},
679 {
680 RowBox[{
681 FractionBox["1", "4"], " ",
682 RowBox[{"(",
683 RowBox[{
684 RowBox[{"2", " ", "mu3"}], "+",
685 RowBox[{"3", " ", "lambda6", " ",
686 SuperscriptBox["vev1", "2"]}], "+",
687 RowBox[{"3", " ",
688 SuperscriptBox["vev1", "2"], " ",
689 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "+",
690 RowBox[{"2", " ",
691 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}],
692 RowBox[{"mu2", "+",
693 RowBox[{
694 FractionBox["1", "2"], " ",
695 RowBox[{"(",
696 RowBox[{"lambda3", "+", "lambda4", "+",
697 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
698 SuperscriptBox["vev1", "2"]}]}], "0"},
699 {
700 RowBox[{
701 FractionBox["1", "4"], " ", "\[ImaginaryI]", " ",
702 RowBox[{"(",
703 RowBox[{
704 RowBox[{"2", " ", "mu3"}], "+",
705 RowBox[{"3", " ", "lambda6", " ",
706 SuperscriptBox["vev1", "2"]}], "-",
707 RowBox[{"3", " ",
708 SuperscriptBox["vev1", "2"], " ",
709 SuperscriptBox["lambda6", "\<\"*\"\>"]}], "-",
710 RowBox[{"2", " ",
711 SuperscriptBox["mu3", "\<\"*\"\>"]}]}], ")"}]}], "0",
712 RowBox[{
713 FractionBox["1", "2"], " ",
714 RowBox[{"(",
715 RowBox[{
716 RowBox[{"2", " ", "mu2"}], "+",
717 RowBox[{
718 RowBox[{"(",
719 RowBox[{"lambda3", "+", "lambda4", "-",
720 RowBox[{"2", " ", "lambda5"}]}], ")"}], " ",
721 SuperscriptBox["vev1", "2"]}]}], ")"}]}]}
722 },
723 GridBoxAlignment->{
724 "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
725 "RowsIndexed" -> {}},
726 GridBoxSpacings->{"Columns" -> {
727 Offset[0.27999999999999997`], {
728 Offset[0.7]},
729 Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
730 Offset[0.2], {
731 Offset[0.4]},
732 Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
733 Function[BoxForm`e$,
734 MatrixForm[BoxForm`e$]]]], "Print",
735 CellChangeTimes->{3.568293581864567*^9}],
736
737Cell[BoxData["\<\"Parameters.cpp and Parameters.hpp output is finished.\\n \
738They have been stored in /Users/Karen/trunk/feynrules-development/2HDM\"\>"], \
739"Print",
740 CellChangeTimes->{3.568293581978196*^9}],
741
742Cell[BoxData["\<\"Checking parameter list...\"\>"], "Print",
743 CellChangeTimes->{3.568293582003612*^9}],
744
745Cell[BoxData["\<\"Warning: Ghost fields will be ignored.\"\>"], "Print",
746 CellChangeTimes->{3.568293582005931*^9}],
747
748Cell[BoxData[
749 InterpretationBox[
750 RowBox[{"\<\"Output written on \"\>",
751 "\[InvisibleSpace]", "\<\"input/externals.dat\"\>"}],
752 SequenceForm["Output written on ", "input/externals.dat"],
753 Editable->False]], "Print",
754 CellChangeTimes->{3.568293582008362*^9}]
755}, Open ]]
756}, Open ]]
757}, Open ]]
758}, Open ]],
759
760Cell[CellGroupData[{
761
762Cell["\<\
763Compile and excecute ASperGe and
764load the output file back into FeynRules\
765\>", "Section",
766 CellChangeTimes->{{3.5678555212604*^9, 3.567855543609857*^9}}],
767
768Cell[BoxData[
769 RowBox[{"RunASperGe", "[", "]"}]], "Input",
770 CellChangeTimes->{{3.56785554597211*^9, 3.567855549491005*^9}}]
771}, Open ]]
772}, Open ]]
773},
774WindowSize->{740, 656},
775WindowMargins->{{Automatic, 26}, {Automatic, 0}},
776FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (August 26, \
7772012)",
778StyleDefinitions->"Default.nb"
779]
780(* End of Notebook Content *)
781
782(* Internal cache information *)
783(*CellTagsOutline
784CellTagsIndex->{}
785*)
786(*CellTagsIndex
787CellTagsIndex->{}
788*)
789(*NotebookFileOutline
790Notebook[{
791Cell[CellGroupData[{
792Cell[579, 22, 225, 6, 123, "Title"],
793Cell[807, 30, 118, 2, 27, "Input"],
794Cell[CellGroupData[{
795Cell[950, 36, 96, 1, 67, "Section"],
796Cell[CellGroupData[{
797Cell[1071, 41, 521, 13, 58, "Input"],
798Cell[CellGroupData[{
799Cell[1617, 58, 264, 4, 20, "Print"],
800Cell[1884, 64, 545, 11, 20, "Print"],
801Cell[2432, 77, 316, 5, 20, "Print"],
802Cell[2751, 84, 250, 4, 20, "Print"],
803Cell[3004, 90, 318, 5, 20, "Print"],
804Cell[3325, 97, 278, 4, 20, "Print"],
805Cell[3606, 103, 248, 4, 20, "Print"],
806Cell[3857, 109, 315, 5, 20, "Print"]
807}, Open ]]
808}, Open ]]
809}, Open ]],
810Cell[CellGroupData[{
811Cell[4233, 121, 103, 1, 67, "Section"],
812Cell[CellGroupData[{
813Cell[4361, 126, 200, 5, 43, "Input"],
814Cell[CellGroupData[{
815Cell[4586, 135, 185, 2, 20, "Print"],
816Cell[4774, 139, 162, 2, 20, "Print"],
817Cell[4939, 143, 294, 6, 20, "Print"],
818Cell[5236, 151, 193, 3, 20, "Print"],
819Cell[5432, 156, 147, 2, 20, "Print"],
820Cell[5582, 160, 175, 2, 20, "Print"],
821Cell[5760, 164, 180, 2, 20, "Print"],
822Cell[5943, 168, 178, 2, 20, "Print"],
823Cell[6124, 172, 345, 7, 36, "Print"]
824}, Open ]]
825}, Open ]]
826}, Open ]],
827Cell[CellGroupData[{
828Cell[6530, 186, 121, 1, 67, "Section"],
829Cell[CellGroupData[{
830Cell[6676, 191, 156, 2, 26, "Subsubsection"],
831Cell[CellGroupData[{
832Cell[6857, 197, 149, 2, 27, "Input"],
833Cell[CellGroupData[{
834Cell[7031, 203, 286, 4, 20, "Print"],
835Cell[7320, 209, 471, 9, 20, "Print"],
836Cell[7794, 220, 384, 5, 36, "Print"],
837Cell[8181, 227, 321, 5, 20, "Print"],
838Cell[8505, 234, 435, 6, 51, "Print"]
839}, Open ]]
840}, Open ]]
841}, Open ]],
842Cell[CellGroupData[{
843Cell[9001, 247, 234, 5, 34, "Subsection"],
844Cell[9238, 254, 227, 4, 27, "Input"],
845Cell[CellGroupData[{
846Cell[9490, 262, 295, 4, 26, "Subsubsection"],
847Cell[CellGroupData[{
848Cell[9810, 270, 269, 5, 27, "Input"],
849Cell[CellGroupData[{
850Cell[10104, 279, 237, 3, 20, "Print"],
851Cell[10344, 284, 442, 8, 20, "Print"],
852Cell[10789, 294, 292, 4, 20, "Print"],
853Cell[11084, 300, 337, 4, 36, "Print"],
854Cell[11424, 306, 3272, 89, 80, "Print"]
855}, Open ]]
856}, Open ]]
857}, Open ]]
858}, Open ]],
859Cell[CellGroupData[{
860Cell[14769, 403, 216, 4, 34, "Subsection"],
861Cell[CellGroupData[{
862Cell[15010, 411, 75, 1, 27, "Input"],
863Cell[CellGroupData[{
864Cell[15110, 416, 472, 11, 20, "Print"],
865Cell[15585, 429, 696, 14, 51, "Print"],
866Cell[16284, 445, 4803, 115, 80, "Print"],
867Cell[21090, 562, 521, 12, 20, "Print"]
868}, Open ]]
869}, Open ]]
870}, Open ]]
871}, Open ]],
872Cell[CellGroupData[{
873Cell[21684, 582, 112, 1, 67, "Section"],
874Cell[CellGroupData[{
875Cell[21821, 587, 266, 5, 35, "Subsection"],
876Cell[CellGroupData[{
877Cell[22112, 596, 135, 2, 27, "Input"],
878Cell[CellGroupData[{
879Cell[22272, 602, 216, 3, 36, "Print"],
880Cell[22491, 607, 104, 1, 20, "Print"],
881Cell[22598, 610, 114, 1, 20, "Print"],
882Cell[22715, 613, 264, 6, 20, "Print"]
883}, Open ]]
884}, Open ]]
885}, Open ]],
886Cell[CellGroupData[{
887Cell[23040, 626, 533, 8, 50, "Subsection"],
888Cell[CellGroupData[{
889Cell[23598, 638, 267, 5, 27, "Input"],
890Cell[CellGroupData[{
891Cell[23890, 647, 3128, 87, 80, "Print"],
892Cell[27021, 736, 207, 3, 36, "Print"],
893Cell[27231, 741, 102, 1, 20, "Print"],
894Cell[27336, 744, 114, 1, 20, "Print"],
895Cell[27453, 747, 264, 6, 20, "Print"]
896}, Open ]]
897}, Open ]]
898}, Open ]]
899}, Open ]],
900Cell[CellGroupData[{
901Cell[27790, 761, 164, 4, 89, "Section"],
902Cell[27957, 767, 123, 2, 27, "Input"]
903}, Open ]]
904}, Open ]]
905}
906]
907*)
908
909(* End of internal cache information *)