DMsimp: dm_s_spin0_ew.fr

File dm_s_spin0_ew.fr, 2.7 KB (added by mawatari, 9 years ago)

Main FR model file for spin0(EW) v1.0

Line 
1(* ****************************************************************** *)
2(* ***** ***** *)
3(* ***** FeynRules model file: Simplified s-channel DM models ***** *)
4(* ***** Author: A. Martini, K. Mawatari ***** *)
5(* ***** J. Wang, C.Zhang ***** *)
6(* ***** ***** *)
7(* ****************************************************************** *)
8
9(* ************************** *)
10(* ***** Change log ***** *)
11(* ************************** *)
12
13(* 2015.08.24 v1.0 - release version. *)
14
15
16(* ************************** *)
17(* ***** Parameters ***** *)
18(* ************************** *)
19M$Parameters = {
20
21(* EW couplings with the spin-0 mediator *)
22 gSh1 == {
23 ParameterType -> External,
24 InteractionOrder -> {DMS, 1},
25 BlockName -> DMINPUTS,
26 TeX -> Subscript[g,Sh1],
27 Description -> "h-Y0 scalar coupling 1",
28 Value -> 1. },
29
30 gSh2 == {
31 ParameterType -> External,
32 InteractionOrder -> {DMS, 1},
33 BlockName -> DMINPUTS,
34 TeX -> Subscript[g,Sh2],
35 Description -> "h-Y0 scalar coupling 2",
36 Value -> 0. },
37
38 gSb == {
39 ParameterType -> External,
40 InteractionOrder -> {DMS, 1},
41 BlockName -> DMINPUTS,
42 TeX -> Subscript[g,Sb],
43 Description -> "B-Y0 scalar coupling",
44 Value -> 1. },
45
46 gPb == {
47 ParameterType -> External,
48 InteractionOrder -> {DMS, 1},
49 BlockName -> DMINPUTS,
50 TeX -> Subscript[g,Pb],
51 Description -> "B-Y0 pseudoscalar coupling",
52 Value -> 0. },
53
54 gSw == {
55 ParameterType -> External,
56 InteractionOrder -> {DMS, 1},
57 BlockName -> DMINPUTS,
58 TeX -> Subscript[g,Sw],
59 Description -> "W-Y0 scalar coupling",
60 Value -> 0. },
61
62 gPw == {
63 ParameterType -> External,
64 InteractionOrder -> {DMS, 1},
65 BlockName -> DMINPUTS,
66 TeX -> Subscript[g,Pw],
67 Description -> "W-Y0 pseudoscalar coupling",
68 Value -> 0. }
69
70};
71
72(* ************************** *)
73(* ***** Lagrangian ***** *)
74(* ************************** *)
75
76L0SMEW := Module[{ii,jj}, ExpandIndices[
77 gSh1/Lambda * DC[Phibar[ii],mu] DC[Phi[ii],mu] Y0 +
78 gSh2/Lambda * MH^2 * (Phibar[jj] Phi[jj] - vev^2 / 2) Y0,
79 FlavorExpand->{SU2D,SU2W}]] +
80 Module[{mu,nu}, ExpandIndices[
81 1/Lambda * FS[B,mu,nu] (gSb FS[B,mu,nu] +
82 gPb Dual[FS][B,mu,nu]) Y0,
83 FlavorExpand->{SU2D,SU2W}]] +
84 Module[{mu,nu,ii}, ExpandIndices[
85 1/Lambda * FS[Wi,mu,nu,ii] (gSw FS[Wi,mu,nu,ii] +
86 gPw Dual[FS][Wi,mu,nu,ii]) Y0,
87 FlavorExpand->{SU2D,SU2W}]] ;
88