topBSM: thu.fr

File thu.fr, 5.2 KB (added by stefankrastanov, 7 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]];