anomalyfreeZprime: AnoFree_ZP_simple_for_NLO_16.fr

File AnoFree_ZP_simple_for_NLO_16.fr, 19.5 KB (added by martinbauer, 5 weeks ago)
Line 
1(***************************************************************************************************************)
2(******                       This is the FeynRules mod-file for an Anomaly free Z' model                 ******)
3(******                                                                                                   ******)
4(******     Author: Sascha Diefenbacher                                                                   ******)
5(******                                                                                                   ******)
6(****** Choose whether Feynman gauge is desired.                                                          ******)
7(****** If set to False, unitary gauge is assumed.                                                          ****)
8(****** Feynman gauge is especially useful for CalcHEP/CompHEP where the calculation is 10-100 times faster. ***)
9(****** Feynman gauge is not supported in MadGraph and Sherpa.                                              ****)
10(***************************************************************************************************************)
11
12(* ************************** *)
13(* *****  Information   ***** *)
14(* ************************** *)
15M$ModelName = "AnoFree_ZP";
16
17M$Information = {Authors -> {"Sascha D. Diefenbacher"},
18             Version -> "1.0",
19             Date -> "26. 8. 2017",
20             Institutions -> {"Uni-Heidelberg"},
21             Emails -> {""}
22             };
23
24(* ************************** *)
25(* *****    Indices     ***** *)
26(* ************************** *)
27
28
29
30
31(* Parameter list *)
32
33M$Parameters = {
34
35(****External Parameters****)
36   
37        MZptarget == {
38         ParameterType    -> External,
39         Value            -> 300,
40         InteractionOrder -> {Zp,1}, 
41         TeX              -> Subscript[M, zptarget],
42         Description      -> "mzp target"
43    }, 
44       
45        Mn1 == {
46         ParameterType    -> External,
47         Value            -> 250,
48         InteractionOrder -> {Mn1,1}, 
49         TeX              -> Subscript[M,n1],
50         Description      -> "DM mass"
51    },
52   
53    xi == {
54         ParameterType    -> External,
55         Value            -> 0.1,
56         TeX              -> "Chi",
57         Description      -> "Some mixing thing-y"
58    },
59   
60    QS0 == {
61         ParameterType    -> External,
62         Value            -> 2.0,
63         InteractionOrder -> {QED,1}, 
64         TeX              -> Subscript[Q,S0],
65         Description      -> "S0 Carge"
66    },
67   
68    QDM == {
69         ParameterType    -> External,
70         Value            -> 1.0,
71         TeX              -> Subscript[Q,DM],
72         Description      -> "DM Carge"
73    },
74   
75   
76    gzp == {
77         ParameterType    -> External,
78         Value            -> 1.0,
79         InteractionOrder -> {Zp,1}, 
80         TeX              -> Subscript[g, zp],
81         Description      -> "Z' gauge coupling"
82    },
83 
84        chi == {
85         ParameterType    -> External,
86         Value            -> Cos[0.0],
87         TeX              -> "cos(hi)",
88         Description      -> "asljdf"
89    },
90   
91        MStarget == {
92         ParameterType    -> External,
93         Value            -> 200,
94         TeX              -> "MStarget",
95         Description      -> "desired MS0"
96    },
97   
98    MHtarget == {
99         ParameterType    -> External,
100         Value            -> 125,
101         TeX              -> "MHtarget",
102         Description      -> "desred MH"
103    },
104   
105    lambHS == {
106         ParameterType    -> External,
107         Value            -> 0.01,
108         InteractionOrder -> {QED,2}, 
109         TeX              -> "lambdaHS",
110         Description      -> "lambda Mixing"
111    },
112 
113
114     
115
116(*************Internal Parameters************)
117       
118        sxi == {
119         ParameterType    -> Internal,
120         Value            -> Sin[xi],
121         TeX              -> "Sin(xi)",
122         Description      -> "Sin(xi)"
123    },
124   
125    cxi == {
126         ParameterType    -> Internal,
127         Value            -> Cos[xi],
128         TeX              -> "Cos(xi)",
129         Description      -> "Cos(xi)"
130    },
131   
132    txi == {
133         ParameterType    -> Internal,
134         Value            -> sxi/cxi,
135         TeX              -> "Tan(xi)",
136         Description      -> "Tan(xi)"
137    },
138   
139        VS == {
140         ParameterType    -> Internal,
141         Value            -> (Sqrt[2]*cxi*MZptarget*Sqrt[-MZ^2 + MZptarget^2 - MZ^2*sw^2*txi^2])/(gzp*Sqrt[-MZ^2 + MZptarget^2]*QS0),
142         InteractionOrder -> {QED,-1}, 
143         TeX              -> Subscript[vev, S0],
144         Description      -> "vev of S0"
145    },
146       
147        lambS == {
148         ParameterType    -> Internal,
149         Value            -> (MStarget^2 + MHtarget^2 + Sqrt[(MStarget^2-MHtarget^2)^2-(2*lambHS*vev*VS)^2])/(2*VS*VS),
150         InteractionOrder -> {QED,2}, 
151         TeX              -> Subscript[#lambda, S0],
152         Description      -> "lambda S0"
153    }, 
154   
155    lambH == {
156         ParameterType    -> Internal,
157         Value            -> (MStarget^2 + MHtarget^2)/(vev*vev) - lambS*VS*VS/vev/vev,
158                 InteractionOrder -> {QED,2}, 
159         TeX              -> Subscript[#lambda, H],
160         Description      -> "lambda H"
161    },
162       
163   
164    zi == {
165         ParameterType    -> Internal,
166         Value            -> (1/2)*ArcTan[2*MZ^2*txi*sw/(MZ^2(1-sw^2*txi^2) - xs*vev^2/(2*cxi^2*VS^2)*(gzp*QS0*VS)^2)],
167         TeX              -> "zi",
168         Description      -> "other mixing thing"
169    },
170   
171    szi == {
172         ParameterType    -> Internal,
173         Value            -> Sin[zi],
174         TeX              -> "Sin(zi)",
175         Description      -> "Sin(zi)"
176    },
177   
178    czi == {
179         ParameterType    -> Internal,
180         Value            -> Cos[zi],
181         TeX              -> "Cos(zi)",
182         Description      -> "Cos(zi)"
183    },
184     
185        alp == {
186         ParameterType    -> Internal,
187         Value            -> (1/2)*ArcTan[2*lambHS*vev*VS/(lambH*vev*vev-lambS*VS*VS)],
188         TeX              -> "alpha",
189         Description      -> "Higgs S0 Mixing"
190    },
191       
192        xs == {
193         ParameterType    -> Internal,
194         Value            -> (VS/vev)^2,
195         TeX              -> "xs",
196         Description      -> "xs"
197    },
198   
199    sal == {
200         ParameterType    -> Internal,
201         Value            -> Sin[alp],
202         TeX              -> "Sin(alpha)",
203         Description      -> "Sin(alpha)"
204    },
205   
206    cal == {
207         ParameterType    -> Internal,
208         Value            -> Cos[alp],
209         TeX              -> "Cos(alpha)",
210         Description      -> "Cos(alpha)"
211    },
212       
213        ozp == {
214         ParameterType    -> Internal,
215         Value            -> 1,
216         InteractionOrder -> {Zp,1}, 
217         TeX              -> "orderzp",
218         Description      -> "orderzp"
219    },
220   
221    MS0 == {
222         ParameterType    -> Internal,
223         Value            -> Sqrt[(1/2)*(lambH*vev*vev+lambS*VS*VS + Sqrt[(lambH*vev*vev-lambS*VS*VS)^2+(2*lambHS*vev*VS)^2])],
224         TeX              -> Subscript[M,S0],
225         Description      -> "Mass S0 after Mixing"
226    },
227       
228        MZp == {
229         ParameterType    -> Internal,
230         Value            -> Sqrt[(1/2)*(vev/vev)^2*(MZ^2*(1+txi^2*sw^2) + xs*vev^2/(2*cxi^2*VS^2)*(gzp*QS0*VS)^2
231                                                + Sqrt[(MZ^2*(1+txi^2*sw^2) + xs*vev^2/(2*cxi^2*VS^2)*(gzp*QS0*VS)^2)^2 - 2*xs*vev^2/(cxi^2*VS^2)*(gzp*QS0*VS)^2*MZ^2])],
232         TeX              -> Subscript[M,zp],
233         Description      -> "Z' Mass"
234    }
235};
236
237
238(*****************************************************************************)
239(*              New fields                                                   *)
240(*****************************************************************************)
241
242(************* New Quarks ***********)
243
244M$ClassesDescription = {
245
246(* Gauge bosons: physical vector fields *) 
247
248
249  S[4] == {
250        ClassName        -> S0,
251        SelfConjugate    -> True,
252        Indices          -> {},
253        Mass             -> {MS0, Internal},
254        Width            -> {WS0, 1.},
255        ParticleName     -> "S0",
256        PDG              -> 200002100,
257        PropagatorLabel  -> "S0",
258        PropagatorType   -> ScalarDash,
259        PropagatorArrow  -> None},
260       
261  F[7] == {
262    ClassName         -> dm,
263    SelfConjugate     -> False,
264    Mass             -> {Mn1, Internal},
265    Width            -> 0,
266    PDG              -> 200002200,
267    ParticleName     -> {"dm"},
268    AntiParticleName -> {"dm~"},
269    TeX               -> "dm",
270    FullName          -> "Dirac DM" },
271 
272  V[5] == {
273    ClassName       -> Zp,
274    SelfConjugate   -> True,
275    Mass            -> {MZp, Internal},
276    Width           -> {WZp, 1},
277    ParticleName    -> "Zp",
278    PDG             -> 23000,
279    PropagatorLabel -> "Zp",
280    PropagatorType  -> Sine,
281    PropagatorArrow -> None,
282    FullName        -> "Zp"
283  }
284};
285
286(*****************************************************************************)
287(*              New Lagrangian Terms                                         *)
288(*****************************************************************************)
289
290(*********************)
291(**** Kinetic terms***)
292(*********************)
293
294LZPkin :=   -1/4 FS[Zp,mu,nu] FS[Zp,mu,nu] + (MZp^2/2) Zp[mu].Zp[mu];
295
296LS0kin := 1/2 del[S0, mu] del[S0, mu] - 1/2 MS0^2 S0^2;
297
298LNewkin := LZPkin + LS0kin;
299
300(*******************************)
301(**** Z Z' H S0 interactions ***)
302(*******************************)
303
304(*LZHSMneg := -chi*(cw^2 + sw^2)^2*vev^2*ee^2/4 Z[mu].Z[mu]  H;*)
305LZHSMneg := -chi*(MZ^2/(vev)) Z[mu].Z[mu]  H;
306
307LHS0ZZ := Z[mu].Z[mu] ( (MZ^2/(vev))*(cal H + sal S0)*(czi+sw*szi*txi)^2 + ((gzp*QS0*VS/ozp)^2/(VS))*(cal S0 - sal H)*(szi^2/cxi^2) );
308
309LHS0ZPZP := Zp[mu].Zp[mu] ( (MZ^2/(vev))*(cal H + sal S0)*(-szi+sw*czi*txi)^2 + ((gzp*QS0*VS/ozp)^2/(VS))*(cal S0 - sal H)*(czi^2/cxi^2) );
310
311LHS0ZZP := 2*( (MZ^2/(vev))*(cal*(Zp[mu].Z[mu] H) + sal*(Zp[mu].Z[mu] S0))*(czi+sw*szi*txi)*(-szi+sw*czi*txi) + ((gzp*QS0*VS/ozp)^2/(VS))*(cal*(Zp[mu].Z[mu] S0) - sal*(Zp[mu].Z[mu] H))*(szi*czi/cxi^2) );
312
313LHS0ZZPtotal := LHS0ZZ + LHS0ZZP + LHS0ZPZP;
314
315
316(*****************************)
317(**** H S0 SM interactions ***)
318(*****************************)
319
320(*LHHHHneg := chi*lam/4*(vev + H)^4;*)
321
322LS0S0HH :=  S0 H H*(-3 cal^2 lambH sal vev + 2 cal^2 lambHS sal vev - lambHS sal^3 vev + cal^3 lambHS VS - 2 cal lambHS sal^2 VS + 3 cal lambS sal^2 VS)
323           +S0 H*ozp*(-3 cal lambH sal vev^2 + cal lambHS sal vev^2 + 2 cal^2 lambHS vev VS - 2 lambHS sal^2 vev VS - cal lambHS sal VS^2 + 3 cal lambS sal VS^2);
324
325LHSMneg := -chi*H*((1/Sqrt[2])*ydo dbar.d + (1/Sqrt[2])*yup ubar.u + (1/Sqrt[2])*ys sbar.s + (1/Sqrt[2])*yc cbar.c + (1/Sqrt[2])*yb bbar.b + (1/Sqrt[2])*yt tbar.t +  (1/Sqrt[2])*ye ebar.e + (1/Sqrt[2])*ym mubar.mu + (1/Sqrt[2])*ytau tabar.ta + ee^2*vev/(2*sw^2) W[mu].Wbar[mu]);
326
327LHS0SM := (cal*H + sal*S0)*((1/Sqrt[2])*ydo dbar.d + (1/Sqrt[2])*yup ubar.u + (1/Sqrt[2])*ys sbar.s + (1/Sqrt[2])*yc cbar.c + (1/Sqrt[2])*yb bbar.b + (1/Sqrt[2])*yt tbar.t +  (1/Sqrt[2])*ye ebar.e + (1/Sqrt[2])*ym mubar.mu + (1/Sqrt[2])*ytau tabar.ta + ee^2*vev/(2*sw^2) W[mu].Wbar[mu]);
328
329LHS0SMtotal :=LHSMneg + LHS0SM + LS0S0HH;
330
331(*****************************)
332(**** Z Z' SM interactions ***)
333(*****************************)
334
335LZJEM :=-cw*szi*txi*ee*( (-1)*ebar.Ga[mu].e Z[mu]   +   (2/3)*ubar.Ga[mu].u Z[mu] + (-1/3)*dbar.Ga[mu].d Z[mu] +
336                      (-1)*mubar.Ga[mu].mu Z[mu] +   (2/3)*cbar.Ga[mu].c Z[mu] + (-1/3)*sbar.Ga[mu].s Z[mu] +
337                      (-1)*tabar.Ga[mu].ta Z[mu] +   (2/3)*tbar.Ga[mu].t Z[mu] + (-1/3)*bbar.Ga[mu].b Z[mu] );
338
339LZPJEM :=-cw*czi*txi*ee*( (-1)*ebar.Ga[mu].e Zp[mu]   +   (2/3)*ubar.Ga[mu].u Zp[mu] + (-1/3)*dbar.Ga[mu].d Zp[mu] +
340                       (-1)*mubar.Ga[mu].mu Zp[mu] +   (2/3)*cbar.Ga[mu].c Zp[mu] + (-1/3)*sbar.Ga[mu].s Zp[mu] +
341                       (-1)*tabar.Ga[mu].ta Zp[mu] +   (2/3)*tbar.Ga[mu].t Zp[mu] + (-1/3)*bbar.Ga[mu].b Zp[mu] );
342
343JZJZSMneg := ee/(sw*cw)*(-chi)*( (1/2)*left[vlbar].Ga[mu].left[vl] Z[mu] + (-(1/2)+sw^2)*left[lbar].Ga[mu].left[l] Z[mu] + (sw^2)*right[lbar].Ga[mu].right[l] Z[mu]
344                                     +((1/2)-(2/3)*sw^2)*left[uqbar].Ga[mu].left[uq] Z[mu] + (-(2/3)*sw^2)*right[uqbar].Ga[mu].right[uq] Z[mu]
345                                     +(-(1/2)+(1/3)*sw^2)*left[dqbar].Ga[mu].left[dq] Z[mu] + ((1/3)*sw^2)*right[dqbar].Ga[mu].right[dq] Z[mu] );
346
347
348
349LZJZ := ee/(sw*cw)*(czi+sw*szi*txi)*( (1/2)*left[vebar].Ga[mu].left[ve] Z[mu] + (-(1/2)+sw^2)*left[ebar].Ga[mu].left[e] Z[mu] + (sw^2)*right[ebar].Ga[mu].right[e] Z[mu]
350                                     +((1/2)-(2/3)*sw^2)*left[ubar].Ga[mu].left[u] Z[mu] + (-(2/3)*sw^2)*right[ubar].Ga[mu].right[u] Z[mu]
351                                     +(-(1/2)+(1/3)*sw^2)*left[dbar].Ga[mu].left[d] Z[mu] + ((1/3)*sw^2)*right[dbar].Ga[mu].right[d] Z[mu]
352                                     +(1/2)*left[vmbar].Ga[mu].left[vm] Z[mu] + (-(1/2)+sw^2)*left[mubar].Ga[mu].left[mu] Z[mu] + (sw^2)*right[mubar].Ga[mu].right[mu] Z[mu]
353                                     +((1/2)-(2/3)*sw^2)*left[cbar].Ga[mu].left[c] Z[mu] + (-(2/3)*sw^2)*right[cbar].Ga[mu].right[c] Z[mu]
354                                     +(-(1/2)+(1/3)*sw^2)*left[sbar].Ga[mu].left[s] Z[mu] + ((1/3)*sw^2)*right[sbar].Ga[mu].right[s] Z[mu]
355                                     +(1/2)*left[vtbar].Ga[mu].left[vt] Z[mu] + (-(1/2)+sw^2)*left[tabar].Ga[mu].left[ta] Z[mu] + (sw^2)*right[tabar].Ga[mu].right[ta] Z[mu]
356                                     +((1/2)-(2/3)*sw^2)*left[tbar].Ga[mu].left[t] Z[mu] + (-(2/3)*sw^2)*right[tbar].Ga[mu].right[t] Z[mu]
357                                     +(-(1/2)+(1/3)*sw^2)*left[bbar].Ga[mu].left[b] Z[mu] + ((1/3)*sw^2)*right[bbar].Ga[mu].right[b] Z[mu] );
358
359LZPJZ := ee/(sw*cw)*(-szi+sw*czi*txi)*((1/2)*left[vebar].Ga[mu].left[ve] Zp[mu] + (-(1/2)+sw^2)*left[ebar].Ga[mu].left[e] Zp[mu] + (sw^2)*right[ebar].Ga[mu].right[e] Zp[mu]
360                                     +((1/2)-(2/3)*sw^2)*left[ubar].Ga[mu].left[u] Zp[mu] + (-(2/3)*sw^2)*right[ubar].Ga[mu].right[u] Zp[mu]
361                                     +(-(1/2)+(1/3)*sw^2)*left[dbar].Ga[mu].left[d] Zp[mu] + ((1/3)*sw^2)*right[dbar].Ga[mu].right[d] Zp[mu]
362                                     +(1/2)*left[vmbar].Ga[mu].left[vm] Zp[mu] + (-(1/2)+sw^2)*left[mubar].Ga[mu].left[mu] Zp[mu] + (sw^2)*right[mubar].Ga[mu].right[mu] Zp[mu]
363                                     +((1/2)-(2/3)*sw^2)*left[cbar].Ga[mu].left[c] Zp[mu] + (-(2/3)*sw^2)*right[cbar].Ga[mu].right[c] Zp[mu]
364                                     +(-(1/2)+(1/3)*sw^2)*left[sbar].Ga[mu].left[s] Zp[mu] + ((1/3)*sw^2)*right[sbar].Ga[mu].right[s] Zp[mu]
365                                     +(1/2)*left[vtbar].Ga[mu].left[vt] Zp[mu] + (-(1/2)+sw^2)*left[tabar].Ga[mu].left[ta] Zp[mu] + (sw^2)*right[tabar].Ga[mu].right[ta] Zp[mu]
366                                     +((1/2)-(2/3)*sw^2)*left[tbar].Ga[mu].left[t] Zp[mu] + (-(2/3)*sw^2)*right[tbar].Ga[mu].right[t] Zp[mu]
367                                     +(-(1/2)+(1/3)*sw^2)*left[bbar].Ga[mu].left[b] Zp[mu] + ((1/3)*sw^2)*right[bbar].Ga[mu].right[b] Zp[mu] );
368
369
370
371LZWWneg  = -chi*(1/2)*gw*((del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *cw*Z[nu]                     -
372                          (del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   cw*Z[mu]                     *(Wbar[nu] - W[nu])/Sqrt[2]/I +
373                          (del[cw*Z[nu],                     mu] -   del[cw*Z[mu],                     nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *(Wbar[nu] - W[nu])/Sqrt[2]/I -
374                                                  (del[cw*Z[nu],                     mu] -   del[cw*Z[mu],                     nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *(W[nu] + Wbar[nu])/Sqrt[2]   +
375                          (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   cw*Z[mu]                     *(W[nu] + Wbar[nu])/Sqrt[2]   -
376                          (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *cw*Z[nu]);                     
377                         
378
379LZWW =   ((czi+sw*szi*txi)  -cw*szi*txi*(sw/cw))*(1/2)*gw*((del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *cw*Z[nu]                    -
380                                                           (del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   cw*Z[mu]                     *(Wbar[nu] - W[nu])/Sqrt[2]/I +
381                                                           (del[cw*Z[nu],                     mu] -   del[cw*Z[mu],                     nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *(Wbar[nu] - W[nu])/Sqrt[2]/I -
382                                                                                   (del[cw*Z[nu],                     mu] -   del[cw*Z[mu],                     nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *(W[nu] + Wbar[nu])/Sqrt[2]   +
383                                                           (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   cw*Z[mu]                     *(W[nu] + Wbar[nu])/Sqrt[2]   -
384                                                           (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *cw*Z[nu]);                     
385                         
386LZPWW =  ((-szi+sw*czi*txi) -cw*czi*txi*(sw/cw))*(1/2)*gw*((del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *cw*Zp[nu]                  -
387                                                           (del[(W[nu] + Wbar[nu])/Sqrt[2],   mu] -   del[(W[mu] + Wbar[mu])/Sqrt[2],   nu])*   cw*Zp[mu]                    *(Wbar[nu] - W[nu])/Sqrt[2]/I +
388                                                           (del[cw*Zp[nu],                    mu] -   del[cw*Zp[mu],                    nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *(Wbar[nu] - W[nu])/Sqrt[2]/I -
389                                                                                   (del[cw*Zp[nu],                    mu] -   del[cw*Zp[mu],                    nu])*   (Wbar[mu] - W[mu])/Sqrt[2]/I *(W[nu] + Wbar[nu])/Sqrt[2]   +
390                                                           (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   cw*Zp[mu]                    *(W[nu] + Wbar[nu])/Sqrt[2]   -
391                                                           (del[(Wbar[nu] - W[nu])/Sqrt[2]/I, mu] -   del[(Wbar[mu] - W[mu])/Sqrt[2]/I, nu])*   (W[mu] + Wbar[mu])/Sqrt[2]   *cw*Zp[nu]);                     
392
393
394
395LZJBmL :=  gzp*(szi/cxi)*((-1)*vebar.Ga[mu].ve + (-1)*ebar.Ga[mu].e   + (1/3)*ubar.Ga[mu].u + (1/3)*dbar.Ga[mu].d+
396                          (-1)*vmbar.Ga[mu].vm + (-1)*mubar.Ga[mu].mu + (1/3)*cbar.Ga[mu].c + (1/3)*sbar.Ga[mu].s+
397                          (-1)*vtbar.Ga[mu].vt + (-1)*tabar.Ga[mu].ta + (1/3)*tbar.Ga[mu].t + (1/3)*bbar.Ga[mu].b)*Z[mu];
398
399LZPJBmL := gzp*(czi/cxi)*((-1)*vebar.Ga[mu].ve + (-1)*ebar.Ga[mu].e   + (1/3)*ubar.Ga[mu].u + (1/3)*dbar.Ga[mu].d+
400                          (-1)*vmbar.Ga[mu].vm + (-1)*mubar.Ga[mu].mu + (1/3)*cbar.Ga[mu].c + (1/3)*sbar.Ga[mu].s+
401                          (-1)*vtbar.Ga[mu].vt + (-1)*tabar.Ga[mu].ta + (1/3)*tbar.Ga[mu].t + (1/3)*bbar.Ga[mu].b)*Zp[mu];
402
403LZJMUTAU :=  gzp*(szi/cxi)*((1)*mubar.Ga[mu].mu + (-1)*tabar.Ga[mu].ta + (1)*vmbar.Ga[mu].vm + (-1)*vtbar.Ga[mu].vt)*Z[mu];
404
405LZPJMUTAU := gzp*(czi/cxi)*((1)*mubar.Ga[mu].mu + (-1)*tabar.Ga[mu].ta + (1)*vmbar.Ga[mu].vm + (-1)*vtbar.Ga[mu].vt)*Zp[mu];
406
407
408LZDM := gzp*(szi/cxi)*QDM*dmbar.Ga[mu].dm*Z[mu];
409
410LZPDM := gzp*(czi/cxi)*QDM*dmbar.Ga[mu].dm*Zp[mu];
411
412LZZPSMtotal := LZJEM + LZPJEM + JZJZSMneg + LZJZ + LZPJZ + LZPDM + LZDM + LZWWneg + LZWW + LZPWW;
413
414
415Ltotal := LSM + LNewkin + LHS0ZZPtotal + LHS0SMtotal + LZZPSMtotal + LZPJMUTAU + LZJMUTAU;