DMsimp: dm_s_spin1.fr

File dm_s_spin1.fr, 10.1 KB (added by mawatari, 4 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;