DMsimp: dm_s_spin1.fr

File dm_s_spin1.fr, 10.1 KB (added by mawatari, 9 years ago)

Main FR model file for spin1 v2.1

Line 
1(* ****************************************************************** *)
2(* ***** ***** *)
3(* ***** FeynRules model file: Simplified DM models ***** *)
4(* ***** s-channel spin-1 mediator ***** *)
5(* ***** Author: A. Martini, K. Mawatari ***** *)
6(* ***** J. Wang, C.Zhang (EW) ***** *)
7(* ***** B. Zaldivar (lepton) ***** *)
8(* ***** B. Fuks (monotop) ***** *)
9(* ***** ***** *)
10(* ****************************************************************** *)
11
12
13(* ************************** *)
14(* ***** Information ***** *)
15(* ************************** *)
16M$ModelName = "DMsimp_s_spin1";
17
18M$Information = {
19 Authors -> {"A. Martini, K. Mawatari, J. Wang, C. Zhang, B. Zaldivar, B. Fuks"},
20 Institutions -> {"Universite catholique de Louvain, LPSC Grenoble, Johnnas Gutenberg University of Mainz, Brookhaven National Laboratory, LAPTh Annecy, LPTHE / UPMC"},
21 Emails -> {"kentarou.mawatari@lpsc.in2p3.fr","cenzhang@bnl.gov"},
22 URLs -> "http://feynrules.irmp.ucl.ac.be/wiki/DMsimp/",
23 References -> {"O. Mattelaer et al., arXiv:1508.00564",
24 "M. Backovic et al., arXiv:1508.05327",
25 "M. Neubert et al., arXiv:1509.05785",
26 "J. Andrea et al., arXiv:1106.6199"},
27 Version -> "2.1",
28 Date -> "2016.10.27"
29};
30
31(* ************************** *)
32(* ***** Change log ***** *)
33(* ************************** *)
34
35(* 2015.08.24 v1.0 - release version. *)
36(* 2016.06.02 v1.1 - changed the PDG numbers for parton-shower programs. *)
37(* 2016.09.23 v2.0 - included the Y1-lepton interactions (B.Zaldivar) *)
38(* 2016.10.27 v2.1 - included the monotop interactions (B. Fuks) *)
39
40
41(***** Setting for interaction order (as e.g. used by MadGraph 5) ******)
42
43M$InteractionOrderLimit = {
44 {DMV, 2}
45};
46
47M$InteractionOrderHierarchy = {
48 {QCD, 1}, {DMV, 2}, {QED, 2}
49};
50
51
52(* ************************** *)
53(* ***** Fields ***** *)
54(* ************************** *)
55M$ClassesDescription = {
56
57S[7] == { ClassName -> Xr,
58 SelfConjugate -> True,
59 Mass -> {MXr, 10.},
60 Width -> 0,
61 PDG -> 5000511,
62 TeX -> Subscript[X,r],
63 FullName -> "Real scalar DM" },
64
65S[8] == { ClassName -> Xc,
66 SelfConjugate -> False,
67 Mass -> {MXc, 10.},
68 Width -> 0,
69 PDG -> 5000512,
70 TeX -> Subscript[X,c],
71 FullName -> "Complex scalar DM" },
72
73F[7] == { ClassName -> Xd,
74 SelfConjugate -> False,
75 Mass -> {MXd, 10.},
76 Width -> 0,
77 PDG -> 5000521,
78 TeX -> Subscript[X,d],
79 FullName -> "Dirac DM" },
80
81V[7] == { ClassName -> Y1,
82 SelfConjugate -> True,
83 Mass -> {MY1, 1000.},
84 Width -> {WY1, 10.},
85 PDG -> 5000001,
86 TeX -> Subscript[Y,1],
87 FullName -> "spin-1 mediator" }
88
89};
90
91
92(* ************************** *)
93(* ***** Parameters ***** *)
94(* ************************** *)
95M$Parameters = {
96
97 gVXc == {
98 ParameterType -> External,
99 InteractionOrder -> {DMV, 1},
100 BlockName -> DMINPUTS,
101 TeX -> Subscript[g,VXc],
102 Description -> "Xc-Y1 vector coupling",
103 Value -> 0. },
104
105 gVXd == {
106 ParameterType -> External,
107 InteractionOrder -> {DMV, 1},
108 BlockName -> DMINPUTS,
109 TeX -> Subscript[g,VXd],
110 Description -> "Xd-Y1 vector coupling",
111 Value -> 1. },
112
113 gAXd == {
114 ParameterType -> External,
115 InteractionOrder -> {DMV, 1},
116 BlockName -> DMINPUTS,
117 TeX -> Subscript[g,AXd],
118 Description -> "Xd-Y1 axial-vector coupling",
119 Value -> 0. },
120
121 gVd11 == {
122 ParameterType -> External,
123 InteractionOrder -> {DMV, 1},
124 BlockName -> DMINPUTS,
125 TeX -> Subscript[g,Vd11],
126 Description -> "d-Y1 vector coupling",
127 Value -> 0.25 },
128
129 gVu11 == {
130 ParameterType -> External,
131 InteractionOrder -> {DMV, 1},
132 BlockName -> DMINPUTS,
133 TeX -> Subscript[g,Vu11],
134 Description -> "u-Y1 vector coupling",
135 Value -> 0.25 },
136
137 gVd22 == {
138 ParameterType -> External,
139 InteractionOrder -> {DMV, 1},
140 BlockName -> DMINPUTS,
141 TeX -> Subscript[g,Vd22],
142 Description -> "s-Y1 vector coupling",
143 Value -> 0.25 },
144
145 gVu22 == {
146 ParameterType -> External,
147 InteractionOrder -> {DMV, 1},
148 BlockName -> DMINPUTS,
149 TeX -> Subscript[g,Vu22],
150 Description -> "c-Y1 vector coupling",
151 Value -> 0.25 },
152
153 gVd33 == {
154 ParameterType -> External,
155 InteractionOrder -> {DMV, 1},
156 BlockName -> DMINPUTS,
157 TeX -> Subscript[g,Vd33],
158 Description -> "b-Y1 vector coupling",
159 Value -> 0.25 },
160
161 gVu33 == {
162 ParameterType -> External,
163 InteractionOrder -> {DMV, 1},
164 BlockName -> DMINPUTS,
165 TeX -> Subscript[g,Vu33],
166 Description -> "t-Y1 vector coupling",
167 Value -> 0.25 },
168
169 gVl11 == {
170 ParameterType -> External,
171 InteractionOrder -> {DMV, 1},
172 BlockName -> DMINPUTS,
173 TeX -> Subscript[g,Ve],
174 Description -> "e-Y1 vector coupling",
175 Value -> 0. },
176
177 gVl22 == {
178 ParameterType -> External,
179 InteractionOrder -> {DMV, 1},
180 BlockName -> DMINPUTS,
181 TeX -> Subscript[g,Vmu],
182 Description -> "mu-Y1 vector coupling",
183 Value -> 0. },
184
185 gVl33 == {
186 ParameterType -> External,
187 InteractionOrder -> {DMV, 1},
188 BlockName -> DMINPUTS,
189 TeX -> Subscript[g,Vta],
190 Description -> "tau-Y1 vector coupling",
191 Value -> 0. },
192
193 gAd11 == {
194 ParameterType -> External,
195 InteractionOrder -> {DMV, 1},
196 BlockName -> DMINPUTS,
197 TeX -> Subscript[g,Ad11],
198 Description -> "d-Y1 axial-vector coupling",
199 Value -> 0. },
200
201 gAu11 == {
202 ParameterType -> External,
203 InteractionOrder -> {DMV, 1},
204 BlockName -> DMINPUTS,
205 TeX -> Subscript[g,Au11],
206 Description -> "u-Y1 axial-vector coupling",
207 Value -> 0. },
208
209 gAd22 == {
210 ParameterType -> External,
211 InteractionOrder -> {DMV, 1},
212 BlockName -> DMINPUTS,
213 TeX -> Subscript[g,Ad22],
214 Description -> "s-Y1 axial-vector coupling",
215 Value -> 0. },
216
217 gAu22 == {
218 ParameterType -> External,
219 InteractionOrder -> {DMV, 1},
220 BlockName -> DMINPUTS,
221 TeX -> Subscript[g,Au22],
222 Description -> "c-Y1 axial-vector coupling",
223 Value -> 0. },
224
225 gAd33 == {
226 ParameterType -> External,
227 InteractionOrder -> {DMV, 1},
228 BlockName -> DMINPUTS,
229 TeX -> Subscript[g,Ad33],
230 Description -> "b-Y1 axial-vector coupling",
231 Value -> 0. },
232
233 gAu33 == {
234 ParameterType -> External,
235 InteractionOrder -> {DMV, 1},
236 BlockName -> DMINPUTS,
237 TeX -> Subscript[g,Au33],
238 Description -> "t-Y1 axial-vector coupling",
239 Value -> 0. },
240
241 gAl11 == {
242 ParameterType -> External,
243 InteractionOrder -> {DMV, 1},
244 BlockName -> DMINPUTS,
245 TeX -> Subscript[g,Ae],
246 Description -> "e-Y1 axial-vector coupling",
247 Value -> 0. },
248
249 gAl22 == {
250 ParameterType -> External,
251 InteractionOrder -> {DMV, 1},
252 BlockName -> DMINPUTS,
253 TeX -> Subscript[g,Amu],
254 Description -> "mu-Y1 axial-vector coupling",
255 Value -> 0. },
256
257 gAl33 == {
258 ParameterType -> External,
259 InteractionOrder -> {DMV, 1},
260 BlockName -> DMINPUTS,
261 TeX -> Subscript[g,Ata],
262 Description -> "tau-Y1 axial-vector coupling",
263 Value -> 0. },
264
265 gnu11 == {
266 ParameterType -> External,
267 InteractionOrder -> {DMV, 1},
268 BlockName -> DMINPUTS,
269 TeX -> Subscript[g,nue],
270 Description -> "nu_e-Y1 coupling",
271 Value -> 0. },
272
273 gnu22 == {
274 ParameterType -> External,
275 InteractionOrder -> {DMV, 1},
276 BlockName -> DMINPUTS,
277 TeX -> Subscript[g,num],
278 Description -> "nu_mu-Y1 coupling",
279 Value -> 0. },
280
281 gnu33 == {
282 ParameterType -> External,
283 InteractionOrder -> {DMV, 1},
284 BlockName -> DMINPUTS,
285 TeX -> Subscript[g,nut],
286 Description -> "nu_tau-Y1 coupling",
287 Value -> 0. },
288
289 gVu31 == {
290 ParameterType -> External,
291 InteractionOrder -> {DMV, 1},
292 BlockName -> DMINPUTS,
293 TeX -> Subscript[g,Vu31],
294 Description -> "FV t-u-Y1 vector coupling",
295 Value -> 0.25 },
296
297 gAu31 == {
298 ParameterType -> External,
299 InteractionOrder -> {DMV, 1},
300 BlockName -> DMINPUTS,
301 TeX -> Subscript[g,Au31],
302 Description -> "FV t-u-Y1 axial-vector coupling",
303 Value -> 0.25 },
304
305 gVd31 == {
306 ParameterType -> External,
307 InteractionOrder -> {DMV, 1},
308 BlockName -> DMINPUTS,
309 TeX -> Subscript[g,Vd31],
310 Description -> "FV b-d-Y1 vector coupling",
311 Value -> 0.25 },
312
313 gAd31 == {
314 ParameterType -> External,
315 InteractionOrder -> {DMV, 1},
316 BlockName -> DMINPUTS,
317 TeX -> Subscript[g,Ad31],
318 Description -> "FV b-d-Y1 axial-vector coupling",
319 Value -> 0.25 }
320
321
322};
323
324
325(* ************************** *)
326(* ***** Lagrangian ***** *)
327(* ************************** *)
328
329L1X := I gVXc/2 (Xcbar del[Xc,mu]-del[Xcbar,mu] Xc) Y1[mu] +
330 Xdbar.Ga[mu].(gVXd + gAXd Ga[5]).Xd Y1[mu];
331
332L1SM := dbar.Ga[mu].(gVd11 + gAd11 Ga[5]).d Y1[mu] +
333 ubar.Ga[mu].(gVu11 + gAu11 Ga[5]).u Y1[mu] +
334 sbar.Ga[mu].(gVd22 + gAd22 Ga[5]).s Y1[mu] +
335 cbar.Ga[mu].(gVu22 + gAu22 Ga[5]).c Y1[mu] +
336 bbar.Ga[mu].(gVd33 + gAd33 Ga[5]).b Y1[mu] +
337 tbar.Ga[mu].(gVu33 + gAu33 Ga[5]).t Y1[mu] +
338 tbar.Ga[mu].(gVu31 + gAu31 Ga[5]).u Y1[mu] +
339 ubar.Ga[mu].(gVu31 + gAu31 Ga[5]).t Y1[mu] +
340 bbar.Ga[mu].(gVd31 + gAd31 Ga[5]).d Y1[mu] +
341 dbar.Ga[mu].(gVd31 + gAd31 Ga[5]).b Y1[mu];
342
343L1SMlep := ebar. Ga[mu].(gVl11 + gAl11 Ga[5]).e Y1[mu] +
344 mubar.Ga[mu].(gVl22 + gAl22 Ga[5]).mu Y1[mu] +
345 tabar.Ga[mu].(gVl33 + gAl33 Ga[5]).ta Y1[mu] +
346 gnu11 vebar.Ga[mu].ProjM.ve Y1[mu] +
347 gnu22 vmbar.Ga[mu].ProjM.vm Y1[mu] +
348 gnu33 vtbar.Ga[mu].ProjM.vt Y1[mu] ;
349
350
351L1DM := L1X + L1SM + L1SMlep;