1 | (***************************************************************)
|
---|
2 | (****** This is the FR file for Top-philic resonance ******)
|
---|
3 | (****** ******)
|
---|
4 | (****** Authors: G.Mohlabeng ******)
|
---|
5 | (****** ******)
|
---|
6 | (****** Choose whether Feynman gauge is desired. ******)
|
---|
7 | (****** If set to False, unitary gauge is assumed. ****)
|
---|
8 | (****** Feynman gauge is especially useful for ******)
|
---|
9 | (****** CalcHEP/CompHEP where the calculation is 10-100 ******)
|
---|
10 | (****** times faster. ******)
|
---|
11 | (****** Feynman gauge is not supported in MadGraph and Sherpa.**)
|
---|
12 | (***************************************************************)
|
---|
13 |
|
---|
14 | (* ************************** *)
|
---|
15 | (* ***** Information ***** *)
|
---|
16 | (* ************************** *)
|
---|
17 | M$ModelName = "Top_Philic";
|
---|
18 |
|
---|
19 | M$Information = {
|
---|
20 | Authors -> {"Jeong Han Kim","G. Mohlabeng"},
|
---|
21 | Version -> "1.0.0",
|
---|
22 | Date -> "04. 01. 2016",
|
---|
23 | Institutions -> {"University of Kansas"},
|
---|
24 | Emails -> {"jeonghan.kim@ku.edu", "gopolang.mohlabeng@ku.edu"}
|
---|
25 | };
|
---|
26 |
|
---|
27 | FeynmanGauge = False;
|
---|
28 |
|
---|
29 |
|
---|
30 | (* ************************** *)
|
---|
31 | (* ***** Change log ***** *)
|
---|
32 | (* ************************** *)
|
---|
33 |
|
---|
34 | (* v1.0: initial version *)
|
---|
35 |
|
---|
36 |
|
---|
37 | (* ************************** *)
|
---|
38 | (* *** Interaction orders *** *)
|
---|
39 | (* *** (as used by mg5) *** *)
|
---|
40 | (* ************************** *)
|
---|
41 |
|
---|
42 | M$InteractionOrderHierarchy = {
|
---|
43 | {Qv1, 8}
|
---|
44 |
|
---|
45 | };
|
---|
46 |
|
---|
47 | (* ************************** *)
|
---|
48 | (* ***** Indices ***** *)
|
---|
49 | (* ************************** *)
|
---|
50 |
|
---|
51 |
|
---|
52 | (* ************************** *)
|
---|
53 | (* **** Particle classes **** *)
|
---|
54 | (* ************************** *)
|
---|
55 | M$ClassesDescription = {
|
---|
56 |
|
---|
57 | (* Vector Fields *)
|
---|
58 |
|
---|
59 | V[7] == {
|
---|
60 | ClassName -> v1,
|
---|
61 | SelfConjugate -> True,
|
---|
62 | Mass -> {Mv1, 1000},
|
---|
63 | Width -> {Wv1, 1.0},
|
---|
64 | ParticleName -> "v1",
|
---|
65 | PDG -> 6000055,
|
---|
66 | PropagatorLabel -> "v1",
|
---|
67 | PropagatorType -> Sine,
|
---|
68 | PropagatorArrow -> None,
|
---|
69 | FullName -> "v1 vector singlet"
|
---|
70 | }
|
---|
71 | };
|
---|
72 |
|
---|
73 |
|
---|
74 | (* ************************** *)
|
---|
75 | (* ***** Gauge ***** *)
|
---|
76 | (* ***** Parameters ***** *)
|
---|
77 | (* ***** (FeynArts) ***** *)
|
---|
78 | (* ************************** *)
|
---|
79 |
|
---|
80 |
|
---|
81 |
|
---|
82 | (* ************************** *)
|
---|
83 | (* ***** Parameters ***** *)
|
---|
84 | (* ************************** *)
|
---|
85 | M$Parameters = {
|
---|
86 |
|
---|
87 |
|
---|
88 | (* Couplings with the vector color singlet *)
|
---|
89 |
|
---|
90 | ct1 == {
|
---|
91 | ParameterType -> External,
|
---|
92 | ParameterName -> ct1,
|
---|
93 | BlockName -> v0PARAMS,
|
---|
94 | Value -> 1.0,
|
---|
95 | InteractionOrder -> {Qv1, 1},
|
---|
96 | TeX -> Subscript[C,t1],
|
---|
97 | Description -> "v1 singlet vector-Fermion coupling"
|
---|
98 | },
|
---|
99 | theta1 == {
|
---|
100 | ParameterType -> External,
|
---|
101 | ParameterName -> theta1,
|
---|
102 | BlockName -> v1PARAMS,
|
---|
103 | Value -> 1.5708,
|
---|
104 | InteractionOrder -> {Qv1, 1},
|
---|
105 | TeX -> Subscript[\[Theta], 1],
|
---|
106 | Description -> "Vector Singlet chirality angle"
|
---|
107 | }
|
---|
108 |
|
---|
109 | };
|
---|
110 |
|
---|
111 | (* ************************** *)
|
---|
112 | (* ***** Lagrangian ***** *)
|
---|
113 | (* ************************** *)
|
---|
114 |
|
---|
115 | (*----------------------*)
|
---|
116 | (* Vectors *)
|
---|
117 | (*----------------------*)
|
---|
118 | (* PL = ProjM, PR = ProjP *)
|
---|
119 | (* Color Singlet vector *)
|
---|
120 |
|
---|
121 | Lv1 := ct1 tbar.Ga[mu].(Cos[theta1] ProjM + Sin[theta1] ProjP).t v1[mu];
|
---|
122 |
|
---|
123 |
|
---|
124 |
|
---|
125 |
|
---|
126 |
|
---|
127 |
|
---|
128 |
|
---|
129 |
|
---|
130 |
|
---|
131 |
|
---|