StandardModelScalars: SMScalars.fr

File SMScalars.fr, 2.8 KB (added by claudeduhr, 9 years ago)

Model file

Line 
1(**************** This is the FeynRules mod-file for the extension of the SM where additional scalars are added **************)
2(*                                                                                                                           *)
3(*        This model file needs the Standard Model model file SM.fr.                                                         *) 
4(*                                                                                                                           *)
5(*        This model can be used in both the unitary and the Feynman gauge.                                                  *)   
6(*                                                                                                                           *)
7(*        The number of scalars is by default set to 4. It can by changed at wish just by changing the variable Nf below.    *)
8(*****************************************************************************************************************************)
9
10Nf = 4;
11
12M$ModelName = "SM_Plus_Scalars";
13
14M$Information = {Authors -> "C. Duhr",
15                 Institutions -> "Universite catholique de Louvain (CP3).",
16                 Emails -> "claude.duhr@uclouvain.be",
17                 Date -> "14. 06. 2009",
18                 Version  ->  "1.0",
19                 URLs  -> "http://feynrules.phys.ucl.ac.be/view/Main/StandardModelScalars",
20                 References -> "\"The minimal non-minimal Standard Model\", J.J. van der Bij, Phys.Lett.B636:56-59,2006, hep-ph/0603082"};
21
22
23
24(******* Index definitions ********)
25
26(* This is the declaration for the generation index of the new scalars *)
27
28IndexRange[ Index[SGen] ] = Range[Nf]
29
30
31IndexStyle[SGen, k]
32
33
34(****************  Parameters *************)
35
36M$Parameters = {
37
38(*  The parameters of the new sector *)
39
40  \[Lambda]S == {
41       InteractionOrder -> {QED, 2},
42       ParameterName -> lS,
43       Value -> 0.5,
44       BlockName -> INVSCALAR},
45
46  \[Omega] == {
47       InteractionOrder -> {QED, 2},
48       ParameterName -> om,
49       Value -> 0.5,
50       BlockName -> INVSCALAR}
51}
52
53
54
55(********* Particle Classes **********)
56
57M$ClassesDescription = {
58   (* The new scalar sector *)
59  S[4] == {
60        ClassName -> Sk,
61        ClassMembers :> Table[Symbol["S" <> ToString[k]], {k, Nf}],
62        SelfConjugate -> True,
63        Indices -> {Index[SGen]},
64        FlavorIndex -> SGen,
65        Mass -> {MSk, 40},
66        Width -> 0}
67}
68
69(*****************************************************************************************)
70
71(* The lagrangian for the new sector *)
72
73 LScalar := 1/2 del[Sk[k], mu]^2 - 1/2 MSk^2 Sk[k]^2   - \[Lambda]S/(8 Nf) Sk[k1]^2 Sk[k2]^2 - \[Omega]/(4 Sqrt[Nf]) Sk[k]^2 HC[Phi].Phi;
74   
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92   
93