TopHeavyRes: PseudoScalarSinglet.fr

File PseudoScalarSinglet.fr, 2.2 KB (added by LucDarme, 4 years ago)
Line 
1(* ****************************************************************** *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + singlet Top-philic pseudo-scalars ***** *)
4(* ***** Author: L. Darme, B. Fuks, ***** *)
5(* ***** ***** *)
6(* ****************************************************************** *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "PseudoScalarSingletTopPhilic";
12M$Information = {
13 Authors -> {"Luc Darme","Benjamin Fuks"},
14 Emails -> {"luc.darme@lnf.infn.it","fuks@lpthe.jussieu.fr"},
15 Institutions -> {"LPTHE / UPMC","LNF, INFN"},
16 Date -> "26.02.2021",
17 Version -> "1",
18 References -> "xxx"};
19
20
21(* ************************** *)
22(* ***** Orders ***** *)
23(* ************************** *)
24M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
25M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
26
27(* ************************** *)
28(* ***** Fields ***** *)
29(* ************************** *)
30M$ClassesDescription = {
31 S[4] == {
32 ClassName -> PS1,
33 SelfConjugate -> True,
34 Mass -> {mPS1,400},
35 Width -> {wPS1,10},
36 ParticleName -> "PS1",
37 PDG -> 9000001,
38 PropagatorLabel -> "PS1",
39 PropagatorType -> ScalarDash,
40 PropagatorArrow -> Forward
41 }
42};
43
44(* ************************** *)
45(* ***** Parameters ***** *)
46(* ************************** *)
47M$Parameters = {
48 yPS1 == {
49 TeX -> Subsuperscript[y,"1P"],
50 ComplexParameter -> False,
51 ParameterType -> External,
52 BlockName -> NPyPS,
53 Value -> 1,
54 InteractionOrder -> {NP,1},
55 Description -> "Pseudo-Scalar singlet coupling"
56 }
57};
58
59(* ************************** *)
60(* ***** Lagrangian ***** *)
61(* ************************** *)
62(* We keep yLu and yRu, set yLu=yRu for the standard pseudo-scalar case *)
63
64L1 := 1/2 DC[Xt,mu] DC[X,mu] - 1/2 m8^2 Xt Xt;
65
66L2a := Block[{sp,sp2,sp3,ff,gg,mm,nn},ExpandIndices[ PS1*yPS1*(-I ProjM[sp,sp2] + I ProjP[sp,sp2]) uqbar[sp,3,m].uq[sp2,3,m]]];
67LagNP := L1 + L2a;