1 | (***************************************************************************************************************)
|
---|
2 | (****** Flavor changing Higgs current feynman-rules model ******)
|
---|
3 | (****** ******)
|
---|
4 | (****** Authors: S. Krastanov ******)
|
---|
5 | (****** ******)
|
---|
6 | (****** Requires the SM.fr module. ******)
|
---|
7 | (***************************************************************************************************************)
|
---|
8 |
|
---|
9 | (* ************************** *)
|
---|
10 | (* ***** Information ***** *)
|
---|
11 | (* ************************** *)
|
---|
12 | M$ModelName = "thu";
|
---|
13 |
|
---|
14 | M$Information = {
|
---|
15 | Authors -> {"Stefan Krastanov"},
|
---|
16 | Version -> "0.1",
|
---|
17 | Date -> "August 1, 2013",
|
---|
18 | Institutions -> {"ENS Lyon", "UCL Belgium"}.
|
---|
19 | Emails -> {"stefan.krastanov@ens-lyon.fr"},
|
---|
20 | URLs -> "http://feynrules.irmp.ucl.ac.be/"
|
---|
21 | };
|
---|
22 |
|
---|
23 | FeynmanGauge = False;
|
---|
24 |
|
---|
25 | (* ************************** *)
|
---|
26 | (* ***** Change log ***** *)
|
---|
27 | (* ************************** *)
|
---|
28 |
|
---|
29 | (* v0.1: Based on SM v1.4.2. *)
|
---|
30 |
|
---|
31 | (* ************************** *)
|
---|
32 | (* ***** Gauge groups ***** *)
|
---|
33 | (* ************************** *)
|
---|
34 |
|
---|
35 |
|
---|
36 | (* ************************** *)
|
---|
37 | (* ***** Indices ***** *)
|
---|
38 | (* ************************** *)
|
---|
39 |
|
---|
40 |
|
---|
41 | (* ************************** *)
|
---|
42 | (* *** Interaction orders *** *)
|
---|
43 | (* *** (as used by mg5) *** *)
|
---|
44 | (* ************************** *)
|
---|
45 |
|
---|
46 |
|
---|
47 | (* ************************** *)
|
---|
48 | (* **** Particle classes **** *)
|
---|
49 | (* ************************** *)
|
---|
50 |
|
---|
51 |
|
---|
52 | (* ************************** *)
|
---|
53 | (* ***** Gauge ***** *)
|
---|
54 | (* ***** Parameters ***** *)
|
---|
55 | (* ***** (FeynArts) ***** *)
|
---|
56 | (* ************************** *)
|
---|
57 |
|
---|
58 |
|
---|
59 | (* ************************** *)
|
---|
60 | (* ***** Parameters ***** *)
|
---|
61 | (* ************************** *)
|
---|
62 | M$Parameters = {
|
---|
63 |
|
---|
64 | LambdaS == {
|
---|
65 | ParameterType -> External,
|
---|
66 | ParameterName -> LambdaS,
|
---|
67 | BlockName -> DIM6,
|
---|
68 | InteractionOrder -> {NP,-1},
|
---|
69 | Value -> 1000,
|
---|
70 | TeX -> \[CapitalLambda],
|
---|
71 | Description -> "Scale of the new physics"},
|
---|
72 | O13uphi == {
|
---|
73 | ParameterType -> External,
|
---|
74 | ParameterName -> O13uphi,
|
---|
75 | BlockName -> DIM6,
|
---|
76 | InteractionOrder -> {QED,2},
|
---|
77 | Value -> 1,
|
---|
78 | TeX -> Subscript[o,13uphi],
|
---|
79 | Description -> "Coupling of the eff operator"},
|
---|
80 | O23uphi == {
|
---|
81 | ParameterType -> External,
|
---|
82 | ParameterName -> O23uphi,
|
---|
83 | BlockName -> DIM6,
|
---|
84 | InteractionOrder -> {QED,2},
|
---|
85 | Value -> 1,
|
---|
86 | TeX -> Subscript[o,23uphi],
|
---|
87 | Description -> "Coupling of the eff operator"},
|
---|
88 | O31uphi == {
|
---|
89 | ParameterType -> External,
|
---|
90 | ParameterName -> O31uphi,
|
---|
91 | BlockName -> DIM6,
|
---|
92 | InteractionOrder -> {QED,2},
|
---|
93 | Value -> 1,
|
---|
94 | TeX -> Subscript[o,31uphi],
|
---|
95 | Description -> "Coupling of the eff operator"},
|
---|
96 | O32uphi == {
|
---|
97 | ParameterType -> External,
|
---|
98 | ParameterName -> O32uphi,
|
---|
99 | BlockName -> DIM6,
|
---|
100 | InteractionOrder -> {QED,2},
|
---|
101 | Value -> 1,
|
---|
102 | TeX -> Subscript[o,32uphi],
|
---|
103 | Description -> "Coupling of the eff operator"},
|
---|
104 | O13ug == {
|
---|
105 | ParameterType -> External,
|
---|
106 | ParameterName -> O13ug,
|
---|
107 | BlockName -> DIM6,
|
---|
108 | InteractionOrder -> {QED,2},
|
---|
109 | Value -> 1,
|
---|
110 | TeX -> Subscript[o,13ug],
|
---|
111 | Description -> "Coupling of the eff operator"},
|
---|
112 | O23ug == {
|
---|
113 | ParameterType -> External,
|
---|
114 | ParameterName -> O23ug,
|
---|
115 | BlockName -> DIM6,
|
---|
116 | InteractionOrder -> {QED,2},
|
---|
117 | Value -> 1,
|
---|
118 | TeX -> Subscript[o,23ug],
|
---|
119 | Description -> "Coupling of the eff operator"},
|
---|
120 | O31ug == {
|
---|
121 | ParameterType -> External,
|
---|
122 | ParameterName -> O31ug,
|
---|
123 | BlockName -> DIM6,
|
---|
124 | InteractionOrder -> {QED,2},
|
---|
125 | Value -> 1,
|
---|
126 | TeX -> Subscript[o,31ug],
|
---|
127 | Description -> "Coupling of the eff operator"},
|
---|
128 | O32ug == {
|
---|
129 | ParameterType -> External,
|
---|
130 | ParameterName -> O32ug,
|
---|
131 | BlockName -> DIM6,
|
---|
132 | InteractionOrder -> {QED,2},
|
---|
133 | Value -> 1,
|
---|
134 | TeX -> Subscript[o,32ug],
|
---|
135 | Description -> "Coupling of the eff operator"}
|
---|
136 | };
|
---|
137 |
|
---|
138 | (* ************************** *)
|
---|
139 | (* ***** Lagrangian ***** *)
|
---|
140 | (* ************************** *)
|
---|
141 | Lthu := Block[{ii, jj, ll, sp, sp1, sp2, cc, cc1, cc2, aa, mu, nu},
|
---|
142 | LO13uphi = QLbar[sp, ii, 1, cc].uR[sp, 3, cc] Phibar[jj] Eps[ii, jj] Phi[ll] Phibar[ll] O13uphi;
|
---|
143 | LO31uphi = QLbar[sp, ii, 3, cc].uR[sp, 1, cc] Phibar[jj] Eps[ii, jj] Phi[ll] Phibar[ll] O31uphi;
|
---|
144 | LO23uphi = QLbar[sp, ii, 2, cc].uR[sp, 3, cc] Phibar[jj] Eps[ii, jj] Phi[ll] Phibar[ll] O23uphi;
|
---|
145 | LO32uphi = QLbar[sp, ii, 3, cc].uR[sp, 2, cc] Phibar[jj] Eps[ii, jj] Phi[ll] Phibar[ll] O32uphi;
|
---|
146 |
|
---|
147 |
|
---|
148 | sigma := I/2 (Ga[mu, sp1, sp] Ga[nu, sp, sp2] - Ga[nu, sp1, sp] Ga[mu, sp, sp2]);
|
---|
149 |
|
---|
150 | LO13ug = QLbar[sp1, ii, 1, cc1] sigma T[aa, cc1, cc2] uR[sp2, 3, cc2] Eps[ii, jj] Phibar[jj] FS[G, mu, nu, aa] O13ug;
|
---|
151 | LO31ug = QLbar[sp1, ii, 3, cc1] sigma T[aa, cc1, cc2] uR[sp2, 1, cc2] Eps[ii, jj] Phibar[jj] FS[G, mu, nu, aa] O31ug;
|
---|
152 | LO23ug = QLbar[sp1, ii, 2, cc1] sigma T[aa, cc1, cc2] uR[sp2, 3, cc2] Eps[ii, jj] Phibar[jj] FS[G, mu, nu, aa] O23ug;
|
---|
153 | LO32ug = QLbar[sp1, ii, 3, cc1] sigma T[aa, cc1, cc2] uR[sp2, 2, cc2] Eps[ii, jj] Phibar[jj] FS[G, mu, nu, aa] O32ug;
|
---|
154 |
|
---|
155 |
|
---|
156 | lthu = ExpandIndices[
|
---|
157 | -(LO13uphi + LO31uphi + LO23uphi + LO32uphi) (Sqrt[2] ymt/vev)^3/LambdaS^2
|
---|
158 | +(LO13ug + LO31ug + LO23ug + LO32ug ) gs Sqrt[2] ymt/vev /LambdaS^2,
|
---|
159 | FlavorExpand -> {SU2D}];
|
---|
160 | lthu + HC[lthu]];
|
---|