topBSM: thu.fr

File thu.fr, 5.2 KB (added by stefankrastanov, 11 years ago)
Line 
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(* ************************** *)
12M$ModelName = "thu";
13
14M$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
23FeynmanGauge = 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(* ************************** *)
62M$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(* ************************** *)
141Lthu := 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]];