DMEffFerm: dm_eff_fermionic_2607.fr

File dm_eff_fermionic_2607.fr, 4.6 KB (added by Ayse Kuday, 3 years ago)
Line 
1(* ****************************************************************** *)
2(* ***** ***** *)
3(* ***** FeynRules model file: Effective Fermionic DM Model ***** *)
4(* ***** Author: A. E. Kuday, F. Ozok, U. Saka ***** *)
5(* ***** ***** *)
6(* ***** ***** *)
7(* ****************************************************************** *)
8
9
10(* ************************** *)
11(* ***** Information ***** *)
12(* ************************** *)
13M$ModelName = "dm_eff_fermionic.fr";
14
15M$Information = {
16Authors -> {"A. E. Kuday F. Ozok, U. Saka"},
17Institutions -> {"Mimar Sinan Fine Arts University, Istanbul University"},
18Emails -> {"aelciboga@gmail.com"},
19URLs -> "http://feynrules.irmp.ucl.ac.be/wiki/DMEFT/",
20References -> {"H. Zhang et al., arXiv:0912.4511v2",
21 "A. E. Kuday et al.},
22Version -> "1.0",
23Date -> "2018.07.26"
24};
25
26(* ************************** *)
27(* ***** Change log ***** *)
28(* ************************** *)
29
30(* 2018.xx.xx v1.0 - release version. *)
31
32(***** Setting for interaction order (as e.g. used by MadGraph 5) ******)
33
34M$InteractionOrderLimit = {
35{DMEFT, 2}
36};
37
38M$InteractionOrderHierarchy = {
39{DMEFT, 1}
40};
41
42
43(* ************************** *)
44(* ***** Fields ***** *)
45(* ************************** *)
46
47M$ClassesDescription = {
48
49F[21] == {
50ClassName -> Xf,
51SelfConjugate -> False,
52Mass -> {MXf, 70},
53Width -> {MXf, 0},
54PropagatorLabel -> "Xf",
55ParticleName -> "Xf",
56FullName -> "Xf"}
57};
58
59
60(* ************************** *)
61(* ***** Parameters ***** *)
62(* ************************** *)
63M$Parameters = {
64
65(* External parameters *)
66
67Lambda == {
68ParameterType -> External,
69ParameterName -> Lambda,
70Value -> 2000,
71TeX -> \[CapitalLambda],
72Description -> "Scale of the new physics"
73},
74
75AlphauX == {
76ParameterType -> External,
77InteractionOrder -> {DMEFT, 1},
78BlockName -> DMINPUTS,
79TeX -> Subscript[Alpha,uX],
80Description -> "Alpha-uX coupling"
81},
82
83AlphadX == {
84ParameterType -> External,
85InteractionOrder -> {DMEFT, 1},
86BlockName -> DMINPUTS,
87TeX -> Subscript[Alpha,dX],
88Description -> "Alpha-dX coupling"
89},
90
91AlphaeX == {
92ParameterType -> External,
93InteractionOrder -> {DMEFT, 1},
94BlockName -> DMINPUTS,
95TeX -> Subscript[Alpha,eX],
96Description -> "Alpha-eX coupling"
97},
98
99AlphaqX == {
100ParameterType -> External,
101InteractionOrder -> {DMEFT, 1},
102BlockName -> DMINPUTS,
103TeX -> Subscript[Alpha,qX],
104Description -> "Alpha-qX coupling"
105},
106
107AlphalX == {
108ParameterType -> External,
109InteractionOrder -> {DMEFT, 1},
110BlockName -> DMINPUTS,
111TeX -> Subscript[Alpha,lX],
112Description -> "Alpha-lX coupling"
113},
114
115AlphaPhiX == {
116ParameterType -> External,
117InteractionOrder -> {DMEFT, 1},
118BlockName -> DMINPUTS,
119TeX -> Subscript[Alpha,PhiX],
120Description -> "Alpha-PhiX coupling"
121},
122
123(* Internal parameters *)
124
125guR == {
126ParameterType -> Internal,
127InteractionOrder -> {DMEFT, 1},
128TeX -> Subscript[g,uR],
129Value -> AlphauX/2,
130Description -> "Xf-uR coupling"},
131
132gdR == {
133ParameterType -> Internal,
134InteractionOrder -> {DMEFT, 1},
135TeX -> Subscript[g,dR],
136Value -> AlphadX/2,
137Description -> "Xf-dR coupling"},
138
139
140geR == {
141ParameterType -> Internal,
142InteractionOrder -> {DMEFT, 1},
143TeX -> Subscript[g,eR],
144Value -> AlphaeX/2,
145Description -> "Xf-eR coupling"},
146
147
148gQl == {
149ParameterType -> Internal,
150InteractionOrder -> {DMEFT, 1},
151TeX -> Subscript[g,QL],
152Value -> -AlphaqX/2,
153Description -> "Xf-Ql coupling"},
154
155gLl == {
156ParameterType -> Internal,
157InteractionOrder -> {DMEFT, 1},
158TeX -> Subscript[g,Ll],
159Value -> -AlphalX/2,
160Description -> "Xf-Ll coupling"}
161};
162
163
164(* ************************** *)
165(* ***** Lagrangian ***** *)
166(* ************************** *)
167
168LEff4f :=( AlphauX/4 Block[{sp1,sp2,sp3,sp4,mu,ff,cc}, (uRbar[sp1,ff,cc].uR[sp2,ff,cc] Ga[mu,sp1,sp2]) Xfbar[sp3].Xf[sp4] Ga[mu,sp3,sp4]] +
169AlphadX/4 Block[{sp1,sp2,sp3,sp4,mu,ff,cc}, (dRbar[sp1,ff,cc].dR[sp2,ff,cc] Ga[mu,sp1,sp2]) Xfbar[sp3].Xf[sp4] Ga[mu,sp3,sp4]] +
170AlphaeX/4 Block[{sp1,sp2,sp3,sp4,mu,ff}, (eRbar[sp1].eR[sp2] Ga[mu,sp1,sp2]) Xfbar[sp3].Xf[sp4] Ga[mu,sp3,sp4]] +
171Block[{sp1,sp3,ii,ff,cc}, ExpandIndices[-AlphaqX/2 QLbar[sp1,ii,ff,cc].Xf[sp1] Xfbar[sp3].QL[sp3,ii,ff,cc] , FlavorExpand->{SU2D}];] +
172Block[{sp1,sp3,ii,ff}, ExpandIndices[-AlphalX/2 LLbar[sp1,ii,ff].Xf[sp1] Xfbar[sp3].LL[sp3,ii,ff], FlavorExpand->{SU2D,SU2W}]] ) / Lambda^2;
173
174
175LEffXfPhi := Block[{ii}, ExpandIndices[
176I/Lambda^2 AlphaPhiX*(Phibar[ii] DC[Phi[ii],mu] - DC[Phibar[ii],mu] Phi[ii]) (Xfbar.Ga[mu].Xf),
177FlavorExpand->{SU2D,SU2W}]];
178
179LEffFDM := LEff4f + LEffXfPhi ;