ASperGe: 2HDM.nb

File 2HDM.nb, 31.1 KB (added by Karen De Causmaecker, 4 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 *)