TopHeavyRes: ScalarSinglet.fr

File ScalarSinglet.fr, 2.2 KB (added by LucDarme, 4 years ago)
Line 
1(* ********************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + singlet Top-philic scalars ***** *)
4(* ***** Author: B. Fuks, L. Darme ***** *)
5(* ***** ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "SingletTopPhilicScalar";
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(* Change log *)
21
22(* ************************** *)
23(* ***** Orders ***** *)
24(* ************************** *)
25M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
26M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
27
28(* ************************** *)
29(* ***** Fields ***** *)
30(* ************************** *)
31M$ClassesDescription = {
32 S[4] == {
33 ClassName -> sig1,
34 SelfConjugate -> True,
35 Mass -> {msig1,400},
36 Width -> {wsig1,10},
37 ParticleName -> "sig1",
38 PDG -> 9000001,
39 PropagatorLabel -> "sig1",
40 PropagatorType -> ScalarDash,
41 PropagatorArrow -> Forward
42 }
43};
44
45(* ************************** *)
46(* ***** Parameters ***** *)
47(* ************************** *)
48M$Parameters = {
49 ynpS1 == {
50 TeX -> Subsuperscript[y,"1S"],
51 ComplexParameter -> False,
52 ParameterType -> External,
53 BlockName -> NPyS1,
54 Value -> 0.1,
55 InteractionOrder -> {NP,1},
56 Description -> "Scalar top-philic singlet coupling"
57 }
58};
59
60(* ************************** *)
61(* ***** Lagrangian ***** *)
62(* ************************** *)
63L1 := 1/2 DC[sig1,mu] DC[sig1,mu] - 1/2 msig1^2 sig1 sig1;
64
65L2a := Block[{sp,sp2,sp3,ff,gg,mm,nn},ExpandIndices[ sig1 ynpS1*( ProjM[sp,sp2] + ProjP[sp,sp2]) uqbar[sp,3,m].uq[sp2,3,m]]];
66
67LagNP := L1 + L2a ;