DMsimp: edm_s_spin1.fr

File edm_s_spin1.fr, 3.8 KB (added by gianmarco, 2 days ago)

Spin-1 mediated DM-e interaction, FeynRules file

Line 
1(* ************************************************************************** *)
2(* ***** FeynRules model file: Simplified DM models ***** *)
3(* ***** DM-e scattering real spin-0 mediator ***** *)
4(* ***** ***** *)
5(* ***** Author: Gian Marco Lucchetti ***** *)
6(* ***** Date: 11/12/2023 ***** *)
7(* ************************************************************************** *)
8
9
10(* ************************** *)
11(* ***** Information ***** *)
12(* ************************** *)
13M$ModelName = "eDMsimp_s_spin1";
14
15M$Information = {
16 Authors -> {"Gian Marco Lucchetti"},
17 Institutions -> {"University of Bologna"},
18 Emails -> {"gianmarco.lucchetti@unibo.it"},
19 URLs -> "http://feynrules.irmp.ucl.ac.be/wiki/",
20 References -> {"Haider Alhazmi et al., arXiv:2006.16252v2"},
21 Version -> "1.0",
22 Date -> "2025.05.29"
23};
24
25(* ************************** *)
26(* ***** Change log ***** *)
27(* ************************** *)
28
29(* none *)
30
31
32(***** Setting for interaction order (as e.g. used by MadGraph 5) ******)
33
34M$InteractionOrderLimit = {
35 {DMC, 2}
36};
37
38M$InteractionOrderHierarchy = {
39 {QCD, 1}, {DMC, 2}, {QED, 2}
40};
41
42
43(* ************************** *)
44(* ***** Fields ***** *)
45(* ************************** *)
46M$ClassesDescription = {
47
48 S[8] == {
49 ClassName -> Xc,
50 SelfConjugate -> False,
51 Mass -> {MXc, 1.},
52 Width -> 0,
53 PDG -> 52,
54 TeX -> Subscript[X,c],
55 FullName -> "Complex scalar DM" },
56
57 F[7] == {
58 ClassName -> Xd,
59 SelfConjugate -> False,
60 Mass -> {MXd, 1.0},
61 Width -> 0,
62 PDG -> 53,
63 TeX -> Subscript[X,d],
64 FullName -> "Dirac DM" },
65
66 V[7] == {
67 ClassName -> Y1,
68 SelfConjugate -> True,
69 Mass -> {MY1, 100.},
70 Width -> 0,
71 PDG -> 54,
72 TeX -> Subscript[Y,1],
73 FullName -> "Spin-1 mediator" }
74};
75
76
77(* ************************** *)
78(* ***** Parameters ***** *)
79(* ************************** *)
80M$Parameters = {
81
82 gVXc =={
83 ParameterType -> External,
84 InteractionOrder -> {DMC, 1},
85 BlockName -> DMINPUTS,
86 TeX -> Subscript[g,VXc],
87 Description -> "Xc-Y1 vector coupling",
88 Value -> 0. },
89
90 gVXd == {
91 ParameterType -> External,
92 InteractionOrder -> {DMC, 1},
93 BlockName -> DMINPUTS,
94 TeX -> Subscript[g,VXd],
95 Description -> "Xd-Y1 vector coupling",
96 Value -> 1. },
97
98 gVe == {
99 ParameterType -> External,
100 InteractionOrder -> {DMC, 1},
101 BlockName -> DMINPUTS,
102 TeX -> Subscript[g,Ve],
103 Description -> "e-Y1 vector coupling",
104 Value -> 1. },
105
106 gAXc =={
107 ParameterType -> External,
108 InteractionOrder -> {DMC, 1},
109 BlockName -> DMINPUTS,
110 TeX -> Subscript[g,AXc],
111 Description -> "Xc-Y1 axial coupling",
112 Value -> 0. },
113
114 gAXd == {
115 ParameterType -> External,
116 InteractionOrder -> {DMC, 1},
117 BlockName -> DMINPUTS,
118 TeX -> Subscript[g,AXd],
119 Description -> "Xd-Y1 axial coupling",
120 Value -> 1. },
121
122 gAe == {
123 ParameterType -> External,
124 InteractionOrder -> {DMC, 1},
125 BlockName -> DMINPUTS,
126 TeX -> Subscript[g,Ae],
127 Description -> "e-Y1 axial coupling",
128 Value -> 1. }
129};
130
131
132(* ************************** *)
133(* ***** Lagrangian ***** *)
134(* ************************** *)
135
136LVDM := - I gVXc/2 (Xcbar del[Xc,mu]-del[Xcbar,mu] Xc) Y1[mu] - gVXd Xdbar.Ga[mu].Xd Y1[mu] - gAXd Xdbar.Ga[mu].Ga[5].Xd Y1[mu];
137
138LVSM := - gVe ebar.Ga[mu].e Y1[mu] - gAe ebar.Ga[mu].Ga[5].e Y1[mu];
139
140LVNew := LVDM + LVSM;