WZPrimeAtNLO: vPrimeNLO.fr

File vPrimeNLO.fr, 7.3 KB (added by richardphysics, 10 months ago)

W' Z' FeynRules model file in Feynman gauge

Line 
1(* ********************************************************* *)
2(* *****                                               ***** *)
3(* *****  FeynRules model file: SM + W_SSM + Z_SSM     ***** *)
4(* *****  Authors: B. Fuks, R. Ruiz                    ***** *)
5(* *****                                               ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* *****  Information   ***** *)
10(* ************************** *)
11M$ModelName = "vPrimeNLO";
12M$Information = { Authors->{"B. Fuks, R. Ruiz"},
13                  Emails->{"benjamin.fuks@iphc.cnrs.fr,rruiz@durham.ac.uk"},
14                  Institutions->{"IPHC Strasbourg / University of Strasbourg, IPPP/University of Durham"},
15                  Date->"2017 January 20",
16                  Version->"1.1",
17                  References->{"B. Fuks, R. Ruiz [arXiv:1701.YYYYY], and references therein."},
18                  URLs->{"feynrules.irmp.ucl.ac.be/"} };
19FeynmanGauge = True;
20
21
22
23
24
25(* ************************** *)
26(* *****     Fields     ***** *)
27(* ************************** *)
28M$ClassesDescription = {
29(* Sequential SM Z prime boson *)
30  V[32] == {
31        ClassName        -> Zp,
32        SelfConjugate    -> True,
33        Mass             -> {MZp,  3000.00},
34        Width            -> {WZp,    89.59},
35        ParticleName     -> "Zp",
36        PDG              -> 32,
37        PropagatorLabel  -> "Zp",
38        PropagatorType   -> Sine,
39        PropagatorArrow  -> None,
40        FullName         -> "Zp"
41  },
42(* Sequential SM W prime boson *)
43  V[34] == {
44        ClassName       ->Wp,
45        SelfConjugate   ->False, 
46        Mass            ->{MWp, 3000.00},
47        Width           ->{WWp,  101.27},
48        ParticleName    ->"Wp+",
49        AntiParticleName->"Wp-",
50        QuantumNumbers  ->{Q->1},
51        PDG             ->34,
52        PropagatorLabel ->"Wp",
53        PropagatorType  ->Sine,
54        PropagatorArrow ->Forward
55        FullName        ->"Wp"
56  }
57};
58
59(* ************************** *)
60(* *****   Parameters   ***** *)
61(* ************************** *)
62M$Parameters = {
63  CRq == { TeX->Subsuperscript[C,q,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRq,
64            Value-> {
65              CRq[1,1]->1.0,    CRq[1,2]->.225773, CRq[1,3]->0.,
66              CRq[2,1]->0,      CRq[2,2]->.97418,  CRq[2,3]->0.,
67              CRq[3,1]->0,      CRq[3,2]->0.,      CRq[3,3]->1. },
68            Description->"Right-handed W' couplings to quarks"},
69  CRl == { TeX->Subsuperscript[C,l,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRl,
70            Value-> {
71              CRl[1,1]->1.,  CRl[1,2]->0., CRl[1,3]->0.,
72              CRl[2,1]->0.,  CRl[2,2]->1., CRl[2,3]->0.,
73              CRl[3,1]->0.,  CRl[3,2]->0., CRl[3,3]->1. },
74            Description->"Right-handed W' couplings to leptons"},
75  CLq == { TeX->Subsuperscript[C,q,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLq,
76            Value-> {
77              CLq[1,1]->.97418,   CLq[1,2]->.225773, CLq[1,3]->0.,
78              CLq[2,1]->-.225773, CLq[2,2]->.97418,  CLq[2,3]->0.,
79              CLq[3,1]->0.,       CLq[3,2]->0.,      CLq[3,3]->1. },
80            Description->"Left-handed W' couplings to quarks"},
81  CLl == { TeX->Subsuperscript[C,l,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLl,
82            Value-> {
83              CLl[1,1]->1.,  CLl[1,2]->0., CLl[1,3]->0.,
84              CLl[2,1]->0.,  CLl[2,2]->1., CLl[2,3]->0.,
85              CLl[3,1]->0.,  CLl[3,2]->0., CLl[3,3]->1. },
86            Description->"Left-handed W' couplings to leptons"
87        },
88
89
90  kL == { ParameterType -> External,
91          Value -> 1.0,
92          BlockName->SSMCOUP,
93          OrderBlock->1,
94          TeX -> Subscript[kappa,L],
95          Description -> "Left-handed W' coupling constant scale factor"
96        },
97
98  kR == { ParameterType -> External,
99          Value -> 0.0,
100          TeX -> Subscript[kappa,R],
101          BlockName->SSMCOUP,
102          OrderBlock->2,
103          Description -> "Right-handed W' coupling constant scale factor"
104        },
105
106(* Z_SSM vector coupling to fermions: gZp,f_V =  (1/2)TL3 - Qf*sW2 *)
107  gZpuV == {
108        ParameterType -> Internal,
109        Value         -> (1/2)*( 1/2) - ( 2/3)*sw2,
110        Description   -> "Vector Zp-u-uBar Coupling"
111  },
112  gZpdV == {
113        ParameterType -> Internal,
114        Value         -> (1/2)*(-1/2) - (-1/3)*sw2,
115        Description   -> "Vector Zp-d-dBar Coupling"
116  },
117  gZpvV == {
118        ParameterType -> Internal,
119        Value         -> (1/2)*( 1/2) - (0)*sw2,
120        Description   -> "Vector Zp-v-vBar Coupling"
121  },
122  gZpeV == {
123        ParameterType -> Internal,
124        Value         -> (1/2)*(-1/2) - (-1)*sw2,
125        Description   -> "Vector Zp-e-eBar Coupling"
126  },
127
128(* Z_SSM axial-vector coupling to fermions: gZp,f_A = (-1/2)TL3 *)
129  gZpuA == {
130        ParameterType -> Internal,
131        Value         -> (-1/2)*(1/2),
132        Description   -> "Axial-vector Zp-u-uBar Coupling"
133  },
134  gZpdA == {
135        ParameterType -> Internal,
136        Value         -> (-1/2)*(-1/2),
137        Description   -> "Axial-vector Zp-d-dBar Coupling"
138  },
139  gZpvA == {
140        ParameterType -> Internal,
141        Value         -> (-1/2)*(1/2),
142        Description   -> "Axial-vector Zp-v-vBar Coupling"
143  },
144  gZpeA == {
145        ParameterType -> Internal,
146        Value         -> (-1/2)*(-1/2),
147        Description   -> "Axial-vector Zp-e-eBar Coupling"
148  },
149
150(* Z_SSM LH coupling to fermions: gZp,f_L = gZpfV - gZpfA  *)
151  gZpuL == {
152        ParameterType -> Internal,
153        Value         -> gZpuV - gZpuA,
154        Description   -> "LH Zp-u-uBar Coupling"
155  },
156  gZpdL == {
157        ParameterType -> Internal,
158        Value         -> gZpdV - gZpdA,
159        Description   -> "LH Zp-d-dBar Coupling"
160  },
161  gZpvL == {
162        ParameterType -> Internal,
163        Value         -> gZpvV - gZpvA,
164        Description   -> "LH Zp-v-vBar Coupling"
165  },
166  gZpeL == {
167        ParameterType -> Internal,
168        Value         -> gZpeV - gZpeA,
169        Description   -> "LH Zp-e-eBar Coupling"
170  },
171
172(* Z_SSM RH coupling to fermions: gZp,f_R = gZpfV + gZpfA  *)
173  gZpuR == {
174        ParameterType -> Internal,
175        Value         -> gZpuV + gZpuA,
176        Description   -> "LH Zp-u-uBar Coupling"
177  },
178  gZpdR == {
179        ParameterType -> Internal,
180        Value         -> gZpdV + gZpdA,
181        Description   -> "LH Zp-d-dBar Coupling"
182  },
183  gZpvR == {
184        ParameterType -> Internal,
185        Value         -> gZpvV + gZpvA,
186        Description   -> "LH Zp-v-vBar Coupling"
187  },
188  gZpeR == {
189        ParameterType -> Internal,
190        Value         -> gZpeV + gZpeA,
191        Description   -> "LH Zp-e-eBar Coupling"
192  }
193
194};
195
196(* ************************** *)
197(* *** Interaction orders *** *)
198(* ************************** *)
199
200(* ZSSM Currents *)
201LZSSMTmpU := -kL*ee/sw/cw *Zp[mu]*(gZpuR* uqbar.Ga[mu].ProjP.uq \
202                           + gZpuL* uqbar.Ga[mu].ProjM.uq );
203LZSSMTmpD := -kL*ee/sw/cw *Zp[mu]*(gZpdR* dqbar.Ga[mu].ProjP.dq \
204                           + gZpdL* dqbar.Ga[mu].ProjM.dq );
205LZSSMTmpV := -kL*ee/sw/cw *Zp[mu]*(gZpvR* vlbar.Ga[mu].ProjP.vl \
206                           + gZpvL* vlbar.Ga[mu].ProjM.vl );
207LZSSMTmpE := -kL*ee/sw/cw *Zp[mu]*(gZpeR*  lbar.Ga[mu].ProjP.l  \
208                           + gZpeL*  lbar.Ga[mu].ProjM.l  );
209LZSSM     := LZSSMTmpU + LZSSMTmpD + LZSSMTmpV + LZSSMTmpE;
210
211(* WSSM Currents *)
212LWSSMTmp := -kL*ee/sw/Sqrt[2]*(vlbar.Ga[mu].ProjM.CLl.l + uqbar.Ga[mu].ProjM.CLq.dq)*Wp[mu] \
213            -kR*ee/sw/Sqrt[2]*(vlbar.Ga[mu].ProjP.CRl.l + uqbar.Ga[mu].ProjP.CRq.dq)*Wp[mu];
214
215LWSSM := LWSSMTmp + HC[LWSSMTmp];
216
217(* Combine Everything *)
218LSSM  := LWSSM + LZSSM;
219LFull := LSM + LSSM;