GeneralFCNTop: topFCNC.fr

File topFCNC.fr, 10.1 KB (added by arturamorim, 21 months ago)
Line 
1(****** This is the FeynRules mod-file for the extension of the SM  *********)
2(*      where a new quark Q is added                                        *)
3(*      This model file needs the Standard Model model file SM.fr.          *)
4(*                                                                          *)
5(*      We use unitary gauge (default in SM.fr)                             *)
6(*                                                                          *)
7(****************************************************************************)
8
9
10M$ModelName = "topFCNC_JLOperators";
11
12(* We work in the unitary gauge (just to make sure we do it again) *)
13FeynmanGauge = False;
14
15M$Information = {Authors -> "A.Amorim",                   
16   Institutions ->
17    "University of Porto & LIP",
18                    Emails -> "artur@lip.pt",
19                    Date -> "25.September.2015",
20                    Version  ->  "1.0"};
21
22(* ************************** *)
23(* *****   Parameters   ***** *)
24(* ************************** *)
25
26M$Parameters = {
27 
28(* External parameters *)
29 
30ReXLut  == { TeX -> Subscript[ReX,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 1, InteractionOrder -> {NP,1} },
31ImXLut  == { TeX -> Subscript[ImX,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 2, InteractionOrder -> {NP,1} },
32ReXRut  == { TeX -> Subscript[ReX,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 3, InteractionOrder -> {NP,1} },
33ImXRut  == { TeX -> Subscript[ImX,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 4, InteractionOrder -> {NP,1} },
34ReXLct  == { TeX -> Subscript[ReX,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 5, InteractionOrder -> {NP,1} },
35ImXLct  == { TeX -> Subscript[ImX,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 6, InteractionOrder -> {NP,1} },
36ReXRct  == { TeX -> Subscript[ReX,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 7, InteractionOrder -> {NP,1} },
37ImXRct  == { TeX -> Subscript[ImX,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 8, InteractionOrder -> {NP,1} },
38ReKLut  == { TeX -> Subscript[ReK,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 9, InteractionOrder -> {NP,1} },
39ImKLut  == { TeX -> Subscript[ImK,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 10, InteractionOrder -> {NP,1} },
40ReKRut  == { TeX -> Subscript[ReK,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 11, InteractionOrder -> {NP,1} },
41ImKRut  == { TeX -> Subscript[ImK,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 12, InteractionOrder -> {NP,1} },
42ReKLct  == { TeX -> Subscript[ReK,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 13, InteractionOrder -> {NP,1} },
43ImKLct  == { TeX -> Subscript[ImK,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 14, InteractionOrder -> {NP,1} },
44ReKRct  == { TeX -> Subscript[ReK,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 15, InteractionOrder -> {NP,1} },
45ImKRct  == { TeX -> Subscript[ImK,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 16, InteractionOrder -> {NP,1} },
46ReZetaLut  == { TeX -> Subscript[ReZeta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 17, InteractionOrder -> {NP,1} },
47ImZetaLut  == { TeX -> Subscript[ImZeta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 18, InteractionOrder -> {NP,1} },
48ReZetaRut  == { TeX -> Subscript[ReZeta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 19, InteractionOrder -> {NP,1} },
49ImZetaRut  == { TeX -> Subscript[ImZeta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 20, InteractionOrder -> {NP,1} },
50ReZetaLct  == { TeX -> Subscript[ReZeta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 21, InteractionOrder -> {NP,1} },
51ImZetaLct  == { TeX -> Subscript[ImZeta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 22, InteractionOrder -> {NP,1} },
52ReZetaRct  == { TeX -> Subscript[ReZeta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 23, InteractionOrder -> {NP,1} },
53ImZetaRct  == { TeX -> Subscript[ImZeta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 24, InteractionOrder -> {NP,1} },
54ReEtaLut  == { TeX -> Subscript[ReEta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 25, InteractionOrder -> {NP,1} },
55ImEtaLut  == { TeX -> Subscript[ImEta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 26, InteractionOrder -> {NP,1} },
56ReEtaRut  == { TeX -> Subscript[ReEta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 27, InteractionOrder -> {NP,1} },
57ImEtaRut  == { TeX -> Subscript[ImEta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 28, InteractionOrder -> {NP,1} },
58ReEtaLct  == { TeX -> Subscript[ReEta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 29, InteractionOrder -> {NP,1} },
59ImEtaLct  == { TeX -> Subscript[ImEta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 30, InteractionOrder -> {NP,1} },
60ReEtaRct  == { TeX -> Subscript[ReEta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 31, InteractionOrder -> {NP,1} },
61ImEtaRct  == { TeX -> Subscript[ImEta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 32, InteractionOrder -> {NP,1} },
62ReLambdaLut  == { TeX -> Subscript[ReLambda,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 33, InteractionOrder -> {NP,1} },
63ImLambdaLut  == { TeX -> Subscript[ImLambda,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 34, InteractionOrder -> {NP,1} },
64ReLambdaRut  == { TeX -> Subscript[ReLambda,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 35, InteractionOrder -> {NP,1} },
65ImLambdaRut  == { TeX -> Subscript[ImLambda,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 36, InteractionOrder -> {NP,1} },
66ReLambdaLct  == { TeX -> Subscript[ReLambda,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 37, InteractionOrder -> {NP,1} },
67ImLambdaLct  == { TeX -> Subscript[ImLambda,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 38, InteractionOrder -> {NP,1} },
68ReLambdaRct  == { TeX -> Subscript[ReLambda,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 39, InteractionOrder -> {NP,1} },
69ImLambdaRct  == { TeX -> Subscript[ImLambda,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 40, InteractionOrder -> {NP,1} }
70
71
72
73
74(* Internal Parameters *)
75
76
77};
78
79(* ************************** *)
80(* *****   Lagrangian   ***** *)
81(* ************************** *)
82
83(* H to glue glue *)
84
85LZFCNC := Block[{la},
86        la = -(gw*(ReXLct+I*ImXLct)/(2*cw))*(cbar.Ga[mu].ProjM.t)*Z[mu]-(gw*(ReXLut+I*ImXLut)/(2*cw))*(ubar.Ga[mu].ProjM.t)*Z[mu]
87        -(gw*(ReXRct+I*ImXRct)/(2*cw))*(cbar.Ga[mu].ProjP.t)*Z[mu]-(gw*(ReXRut+I*ImXRut)/(2*cw))*(ubar.Ga[mu].ProjP.t)*Z[mu]
88        +(gw*(ReKLct+I*ImKLct)*I/(4*cw*MZ))*(cbar.Ga[mu].Ga[nu].ProjM.t)*del[Z[nu],mu] - (gw*(ReKLct+I*ImKLct)*I/(4*cw*MZ))*(cbar.Ga[nu].Ga[mu].ProjM.t)*del[Z[nu],mu]
89        +(gw*(ReKRct+I*ImKRct)*I/(4*cw*MZ))*(cbar.Ga[mu].Ga[nu].ProjP.t)*del[Z[nu],mu] - (gw*(ReKRct+I*ImKRct)*I/(4*cw*MZ))*(cbar.Ga[nu].Ga[mu].ProjP.t)*del[Z[nu],mu]
90        +(gw*(ReKLut+I*ImKLut)*I/(4*cw*MZ))*(ubar.Ga[mu].Ga[nu].ProjM.t)*del[Z[nu],mu] - (gw*(ReKLut+I*ImKLut)*I/(4*cw*MZ))*(ubar.Ga[nu].Ga[mu].ProjM.t)*del[Z[nu],mu]
91        +(gw*(ReKRut+I*ImKRut)*I/(4*cw*MZ))*(ubar.Ga[mu].Ga[nu].ProjP.t)*del[Z[nu],mu] - (gw*(ReKRut+I*ImKRut)*I/(4*cw*MZ))*(ubar.Ga[nu].Ga[mu].ProjP.t)*del[Z[nu],mu];
92        la + HC[la]];
93
94LphotonFCNC := Block[{la},
95        la =((ReLambdaLct+I*ImLambdaLct)*I*ee/(2*MT))*(cbar.Ga[mu].Ga[nu].ProjM.t)*del[A[nu],mu]-((ReLambdaLct+I*ImLambdaLct)*I*ee/(2*MT))*(cbar.Ga[nu].Ga[mu].ProjM.t)*del[A[nu],mu]
96        +((ReLambdaRct+I*ImLambdaRct)*I*ee/(2*MT))*(cbar.Ga[mu].Ga[nu].ProjP.t)*del[A[nu],mu]-((ReLambdaRct+I*ImLambdaRct)*I*ee/(2*MT))*(cbar.Ga[nu].Ga[mu].ProjP.t)*del[A[nu],mu]
97        +((ReLambdaLut+I*ImLambdaLut)*I*ee/(2*MT))*(ubar.Ga[mu].Ga[nu].ProjM.t)*del[A[nu],mu]-((ReLambdaLut+I*ImLambdaLut)*I*ee/(2*MT))*(ubar.Ga[nu].Ga[mu].ProjM.t)*del[A[nu],mu]
98        +((ReLambdaRut+I*ImLambdaRut)*I*ee/(2*MT))*(ubar.Ga[mu].Ga[nu].ProjP.t)*del[A[nu],mu]-((ReLambdaRut+I*ImLambdaRut)*I*ee/(2*MT))*(ubar.Ga[nu].Ga[mu].ProjP.t)*del[A[nu],mu];
99        la + HC[la]];
100
101LgluonFCNC := Block[{la},
102        la = (gs*(ReZetaLct+I*ImZetaLct)*I/(2*MT))*2*(cbar.T[a].Ga[mu].Ga[nu].ProjM.t)*del[G[nu,a],mu]-(gs*(ReZetaLct+I*ImZetaLct)*I/(2*MT))*2*(cbar.T[a].Ga[nu].Ga[mu].ProjM.t)*del[G[nu,a],mu]
103        +(gs*(ReZetaRct+I*ImZetaRct)*I/(2*MT))*2*(cbar.T[a].Ga[mu].Ga[nu].ProjP.t)*del[G[nu,a],mu]-(gs*(ReZetaRct+I*ImZetaRct)*I/(2*MT))*2*(cbar.T[a].Ga[nu].Ga[mu].ProjP.t)*del[G[nu,a],mu]
104        +(gs*(ReZetaLut+I*ImZetaLut)*I/(2*MT))*2*(ubar.T[a].Ga[mu].Ga[nu].ProjM.t)*del[G[nu,a],mu]-(gs*(ReZetaLut+I*ImZetaLut)*I/(2*MT))*2*(ubar.T[a].Ga[nu].Ga[mu].ProjM.t)*del[G[nu,a],mu]
105        +(gs*(ReZetaRut+I*ImZetaRut)*I/(2*MT))*2*(ubar.T[a].Ga[mu].Ga[nu].ProjP.t)*del[G[nu,a],mu]-(gs*(ReZetaRut+I*ImZetaRut)*I/(2*MT))*2*(ubar.T[a].Ga[nu].Ga[mu].ProjP.t)*del[G[nu,a],mu];
106        la + HC[la]];
107
108LHiggsFCNC := Block[{la},
109        la = -((ReEtaLut+I*ImEtaLut)/Sqrt[2])*(ubar.ProjM.t)*H-((ReEtaRut+I*ImEtaRut)/Sqrt[2])*(ubar.ProjP.t)*H
110        -((ReEtaLct+I*ImEtaLct)/Sqrt[2])*(cbar.ProjM.t)*H-((ReEtaRct+I*ImEtaRct)/Sqrt[2])*(cbar.ProjP.t)*H;
111        la + HC[la]];
112
113
114Ltot := LSM + LZFCNC + LgluonFCNC + LHiggsFCNC + LphotonFCNC;
115