DMsimp: dm_s_spin0_ew.fr

File dm_s_spin0_ew.fr, 2.7 KB (added by mawatari, 5 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