LeptoQuark: vector_LQ.fr

File vector_LQ.fr, 7.1 KB (added by JavierFuentes, 2 years ago)

Changed PDG code for vector LQ

Line 
1(* ************************************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: vector_LQ ***** *)
4(* ***** Authors: M. Baker, J. Fuentes-Martin, M. König, G. Isidori ***** *)
5(* ***** ***** *)
6(* ************************************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11
12M$ModelName = "vector_LQ";
13
14M$Information = {Authors -> {"M. Baker", "J. Fuentes-Martin", "M. Koenig", "G. Isidori"},
15 Version -> "1.1",
16 Date -> "07.04.2021",
17 Institutions -> {"U. Zurich", "Johannes Gutenberg University Mainz"},
18 Emails -> {"baker@physik.uzh.ch", "jfuentes@uni-mainz.de", "matthias.koenig@tum.de", "isidori@physik.uzh.ch"}
19 };
20
21(* ************************** *)
22(* ***** Change log ***** *)
23(* ************************** *)
24
25(* v1.1: Fixed bug with vector leptoquark covariant derivative. *)
26
27(* ************************** *)
28(* *** Interaction orders *** *)
29(* ************************** *)
30
31M$InteractionOrderHierarchy = { {QCD,1}, {NP,1}, {QED,2} };
32
33
34(* ************************** *)
35(* **** Particle classes **** *)
36(* ************************** *)
37
38M$ClassesDescription = {
39 (* Physical leptoquark field *)
40 V[5] == {
41 ClassName -> VLQ,
42 SelfConjugate -> False,
43 QuantumNumbers -> {Q -> 2/3, LeptonNumber -> -1},
44 Indices -> {Index[Colour]},
45 Mass -> {MVLQ, 3000},
46 Width -> {WVLQ, 600},
47 ParticleName -> "VLQ",
48 PDG -> 42,
49 PropagatorLabel -> "VLQ",
50 PropagatorType -> Sine,
51 PropagatorArrow -> None,
52 FullName -> "Vector leptoquark U1"
53 },
54 (* Unphysical leptoquark field (needed to get the right covariant derivative) *)
55 V[13] == {
56 ClassName -> U,
57 Unphysical -> True,
58 SelfConjugate -> False,
59 QuantumNumbers -> {Y -> 2/3},
60 Indices -> {Index[Colour]},
61 Definitions -> {U[mu_,cc_] -> VLQ[mu,cc]}
62 }
63};
64
65(* ************************** *)
66(* ***** Parameters ***** *)
67(* ************************** *)
68
69M$Parameters = {
70(* External Parameters *)
71 gU == {
72 ParameterType -> External,
73 Value -> 3.0,
74 InteractionOrder -> {NP,1},
75 TeX -> Subscript[g,U],
76 BlockName -> NPLQCOUP,
77 OrderBlock -> 1,
78 Description -> "Overall LQ coupling strength"
79 },
80
81 betaL33 == {
82 ParameterType -> External,
83 Value -> 1.0,
84 InteractionOrder -> {NP,0},
85 TeX -> Subscript[Superscript[\[Beta],L],33],
86 BlockName -> NPLQCOUP,
87 OrderBlock -> 2,
88 Description -> "LQ LH b tau coupling"
89 },
90
91 betaRd33 == {
92 ParameterType -> External,
93 Value -> 1.0,
94 InteractionOrder -> {NP,0},
95 TeX -> Subscript[Superscript[\[Beta],R],33],
96 BlockName -> NPLQCOUP,
97 OrderBlock -> 3,
98 Description -> "LQ RH b tau coupling"
99 },
100
101 betaL23 == {
102 ParameterType -> External,
103 Value -> 0.0,
104 InteractionOrder -> {NP,0},
105 TeX -> Subscript[Superscript[\[Beta],L],23],
106 BlockName -> NPLQCOUP,
107 OrderBlock -> 4,
108 Description -> "LQ LH s tau coupling"
109 },
110
111 betaL32 == {
112 ParameterType -> External,
113 Value -> 0.0,
114 InteractionOrder -> {NP,0},
115 TeX -> Subscript[Superscript[\[Beta],L],32],
116 BlockName -> NPLQCOUP,
117 OrderBlock -> 5,
118 Description -> "LQ LH b mu coupling"
119 },
120
121 kappaU == {
122 ParameterType -> External,
123 Value -> 0.0,
124 InteractionOrder -> {NP,0},
125 TeX -> Subscript[\[Kappa],U],
126 BlockName -> NPLQCOUP,
127 OrderBlock -> 6,
128 Description -> "Non-minimal LQ coupling to gluons"
129 },
130
131 kappaUtilde == {
132 ParameterType -> External,
133 Value -> 0.0,
134 InteractionOrder -> {NP,0},
135 TeX -> Subscript[OverTilde[\[Kappa]],U],
136 BlockName -> NPLQCOUP,
137 OrderBlock -> 7,
138 Description -> "Non-minimal LQ coupling to B gauge boson"
139 },
140
141 (* Internal Parameters *)
142 betaL == {
143 ParameterType -> Internal,
144 Indices -> {Index[Generation], Index[Generation]},
145 Definitions -> {betaL[1, 2] :> 0, betaL[i_?NumericQ, 1] :> 0, betaL[2,2] -> 0},
146 Value -> {
147 betaL[1,3] -> Conjugate[newCKM[3,1]/newCKM[3,2]]*betaL23,
148 betaL[2,3] -> betaL23,
149 betaL[3,2] -> betaL32, betaL[3,3] -> betaL33
150 },
151 InteractionOrder -> {NP,0},
152 TeX -> Subscript[\[Beta],L],
153 Description -> "LQ LH coupling matrix"
154 },
155
156 betaRd == {
157 ParameterType -> Internal,
158 Indices -> {Index[Generation], Index[Generation]},
159 Definitions -> {betaRd[1, j_?NumericQ] :> 0, betaRd[i_?NumericQ, 1] :> 0, betaRd[2,i_?NumericQ] :> 0, betaRd[i_?NumericQ,2] :> 0 },
160 Value -> {betaRd[3,3] -> betaRd33},
161 InteractionOrder -> {NP,0},
162 TeX -> Subscript[\[Beta],Rd],
163 Description -> "LQ RH down-type coupling matrix"
164 },
165
166 (* N. B. : Inputs from UTFit NP fit 2018 *)
167 newCKM == {
168 ParameterType -> Internal,
169 Unitary -> True,
170 Indices -> {Index[Generation], Index[Generation]},
171 Value -> {newCKM[1,1] -> 0.97431, newCKM[1,2] -> 0.22514, newCKM[1,3] -> 0.00371 * Exp[-68.8*I*Pi/180],
172 newCKM[2,1] -> -0.22504 * Exp[0.0341*I*Pi/180], newCKM[2,2] -> 0.97353 * Exp[-0.00182*I*Pi/180], newCKM[2,3] -> 0.0397,
173 newCKM[3,1] -> 0.00837 * Exp[-23.6*I*Pi/180], newCKM[3,2] -> -0.0390 * Exp[1.137*I*Pi/180], newCKM[3,3] -> 0.999200},
174 TeX -> Superscript[V,CKM'],
175 Description -> "CKM-Matrix with 1-3 and 2-3 mixing"
176 }
177};
178
179
180(* ************************** *)
181(* ***** Lagrangian ***** *)
182(* ************************** *)
183
184LVLQKin := Block[{mu,nu,cc1}, - 1/2 * (DC[Ubar[mu,cc1], nu] - DC[Ubar[nu,cc1], mu]) * (DC[U[mu,cc1], nu] - DC[U[nu,cc1], mu]) ];
185
186LVLQF := Block[{ff1,ff2,ff3,s1,s2,s3,cc1,mu},
187 gU/(Sqrt[2]) * VLQ[mu,cc1] *
188 (
189 betaL[ff1,ff2] * newCKM[ff3,ff1] * uqbar[s1,ff3,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * vl[s3,ff2] +
190 betaL[ff1,ff2] * dqbar[s1,ff1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * l[s3,ff2] +
191 betaRd[ff1,ff2] * dqbar[s1,ff1,cc1] * Ga[mu,s1,s2] * ProjP[s2,s3] * l[s3,ff2]
192 )
193 ];
194
195LVLQG := Block[{cc1,cc2,aa1,mu,nu},
196 - I * gs * (1 - kappaU) * Ubar[mu,cc1] * T[aa1,cc1,cc2] * U[nu,cc2] * FS[G,mu,nu,aa1] - I * 2/3 * g1 * (1 - kappaUtilde) * Ubar[mu,cc1] * U[nu,cc1] * FS[B,mu,nu]
197 ];
198
199LVLQ := LVLQKin + LVLQF + HC[LVLQF] + LVLQG;