NLOModels: VLQ_v3.fr

File VLQ_v3.fr, 8.2 KB (added by BenjF, 11 months ago)
Line 
1(***************************************************************************************************************)
2(******                       FeynRules mod-file for Model Independent searches of top partners           ******)
3(******                         X(5/3), T(2/3), B(-1/3) & Y(-4/3) with arbitrary couplings                ******)
4(******                                                                                                   ******)
5(******     Authors: B. Fuks, L. Panizzi                                                                  ******)
6(******     Compatibility with the new SM.fr of FR 2.0                                                    ******)
7(***************************************************************************************************************)
8
9M$ModelName = "VLQ";
10
11M$Information = {
12  Authors -> {"L. Panizzi, B. Fuks"},
13  Version -> "3.0",
14  Date -> "08.10.2016"
15};
16
17FeynmanGauge =False;
18
19M$InteractionOrderHierarchy = { {QCD, 1},  {QED, 2},  {VLQ,1}};
20M$InteractionOrderLimit     = { {QCD, 99}, {QED, 99}, {VLQ,99}};
21
22
23
24(* ************************** *)
25(* *****   Parameters   ***** *)
26(* ************************** *)
27M$Parameters = {
28  KTLw == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "TuW coupling (left-handed)",
29    BlockName -> KTLW, Value -> {KTLw[1] -> 0.120112, KTLw[2] -> 0.120112, KTLw[3] -> 0.160149} },
30  KTRw == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "TuW coupling (right-handed)",
31    BlockName -> KTRW, Value -> {KTRw[1] -> 0., KTRw[2] -> 0., KTRw[3] -> 0.} },
32  KTLz == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "TuZ coupling (left-handed)" ,
33    BlockName -> KTLZ, Value -> {KTLz[1] -> 0.180284, KTLz[2] -> 0.180284, KTLz[3] -> 0.240379} },
34  KTRz == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "TuZ coupling (right-handed)",
35    BlockName -> KTRZ, Value -> {KTRz[1] -> 0., KTRz[2] -> 0., KTRz[3] -> 0.} },
36  KTLh == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "TuH coupling (left-handed)" ,
37    BlockName -> KTLH, Value -> {KTLh[1] -> 0.479336, KTLh[2] -> 0.479336, KTLh[3] -> 0.639117} },
38  KTRh == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "TuH coupling (right-handed)",
39    BlockName -> KTRH, Value -> {KTRh[1] -> 0., KTRh[2] -> 0., KTRh[3] -> 0.} },
40
41  KBLw == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "BuW coupling (left-handed)" ,
42    BlockName -> KBLW, Value -> {KBLw[1] -> 0.120112, KBLw[2] -> 0.120112, KBLw[3] -> 0.160149} },
43  KBRw == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "BuW coupling (right-handed)",
44    BlockName -> KBRW, Value -> {KBRw[1] -> 0., KBRw[2] -> 0., KBRw[3] -> 0.} },
45  KBLz == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "BuZ coupling (left-handed)" ,
46    BlockName -> KBLZ, Value -> {KBLz[1] -> 0.180284, KBLz[2] -> 0.180284, KBLz[3] -> 0.240379} },
47  KBRz == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {{VLQ, 1}, {QED, -1}}, Description -> "BuZ coupling (right-handed)",
48    BlockName -> KBRZ, Value -> {KBRz[1] -> 0., KBRz[2] -> 0., KBRz[3] -> 0.} },
49  KBLh == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "BuH coupling (left-handed)" ,
50    BlockName -> KBLH, Value -> {KBLh[1] -> 0.479336, KBLh[2] -> 0.479336, KBLh[3] -> 0.639117} },
51  KBRh == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "BuH coupling (right-handed)",
52    BlockName -> KBRH, Value -> {KBRh[1] -> 0., KBRh[2] -> 0., KBRh[3] -> 0.} },
53
54  KXL  == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "XuW coupling (left-handed)" ,
55    BlockName -> KXLW, Value -> {KXL[1] -> 0.300279, KXL[2] -> 0.300279, KXL[3] -> 0.400379} },
56  KXR  == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "XuW coupling (right-handed)",
57    BlockName -> KXRW, Value -> {KXR[1] -> 0., KXR[2] -> 0., KXR[3] -> 0.} },
58
59  KYL  == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "YdW coupling (left-handed)" ,
60    BlockName -> KYLW, Value -> {KYL[1] -> 0.300279, KYL[2] -> 0.300279, KYL[3] -> 0.400379} },
61  KYR  == { ParameterType -> External, ComplexParameter -> False, Indices->{Index[Generation]}, InteractionOrder -> {VLQ, 1}, Description -> "YdW coupling (right-handed)",
62    BlockName -> KYRW, Value -> {KYR[1] -> 0., KYR[2] -> 0., KYR[3] -> 0.} }
63};
64
65
66
67(* ************************** *)
68(* *****   Fields       ***** *)
69(* ************************** *)
70
71M$ClassesDescription = {
72  (* VLQ Quarks X, Q=5/3*)
73  F[5] == { ClassName -> x,  SelfConjugate -> False, Indices -> {Index[Colour]}, Mass -> {MX,600},  Width -> {WX, 5.}, QuantumNumbers -> {Q -> 5/3},
74    PropagatorLabel -> "x",  PropagatorType -> Straight, PropagatorArrow -> Forward, PDG -> 6000005, FullName -> "X-quark"},
75
76  (* VLQ Quarks T, Q=2/3 *)
77  F[6] == { ClassName -> tp, SelfConjugate -> False, Indices -> {Index[Colour]}, Mass -> {MTP,600}, Width -> {WTP,5.}, QuantumNumbers -> {Q -> 2/3},
78    PropagatorLabel -> "tp", PropagatorType -> Straight, PropagatorArrow -> Forward, PDG -> 6000006, FullName -> "T-quark"},
79
80  (* VLQ Quarks B, Q=-1/3 *)
81  F[7] == { ClassName -> bp, SelfConjugate -> False, Indices -> {Index[Colour]}, Mass -> {MBP,600}, Width -> {WBP, 5.}, QuantumNumbers -> {Q -> -1/3},
82    PropagatorLabel -> "bp", PropagatorType -> Straight, PropagatorArrow -> Forward, PDG -> 6000007, FullName -> "B-quark"},
83
84  (* VLQ Quarks Y, Q=-4/3 *)
85  F[8] == { ClassName -> y,  SelfConjugate -> False, Indices -> {Index[Colour]}, Mass -> {MY,600},  Width -> {WY, 5.}, QuantumNumbers -> {Q -> -4/3},
86    PropagatorLabel -> "y", PropagatorType -> Straight, PropagatorArrow -> Forward, PDG -> 6000008, FullName -> "Y-quark"}
87};
88
89
90(* ************************** *)
91(* *****   Lagrangian   ***** *)
92(* ************************** *)
93
94(* Kinetic and mass terms *)
95LTK := I tpbar.Ga[mu].DC[tp, mu] - MTP tpbar.tp + 2*ee/3 A[mu] tpbar.Ga[mu].tp;
96LBK := I bpbar.Ga[mu].DC[bp, mu] - MBP bpbar.bp -   ee/3 A[mu] bpbar.Ga[mu].bp;
97LXK := I xbar.Ga[mu].DC[x, mu]   - MX xbar.x + 5*ee/3 A[mu] xbar.Ga[mu].x;
98LYK := I ybar.Ga[mu].DC[y, mu]   - MY ybar.y - 4*ee/3 A[mu] ybar.Ga[mu].y;
99
100(* T interactions *)
101LTW := gw/Sqrt[2] W[mu] Ga[mu,sp1,sp2] tpbar[sp1,cc].dq[sp3,ff,cc] (ProjM[sp2,sp3] KTLw[ff] + ProjP[sp2,sp3] KTRw[ff]);
102LTZ := gw/(2 cw)  Z[mu] Ga[mu,sp1,sp2] tpbar[sp1,cc].uq[sp3,ff,cc] (ProjM[sp2,sp3] KTLz[ff] + ProjP[sp2,sp3] KTRz[ff]);
103LTH :=            H     tpbar[sp1,cc].uq[sp3,ff,cc] (ProjP[sp1,sp3] KTLh[ff] + ProjM[sp1,sp3] KTRh[ff]);
104
105
106(* B interactions *)
107LBW := gw/Sqrt[2] Wbar[mu] Ga[mu,sp1,sp2] bpbar[sp1,cc].uq[sp3,ff,cc] (ProjM[sp2,sp3] KBLw[ff] + ProjP[sp2,sp3] KBRw[ff]);
108LBZ := gw/(2 cw)  Z[mu]    Ga[mu,sp1,sp2] bpbar[sp1,cc].dq[sp3,ff,cc] (ProjM[sp2,sp3] KBLz[ff] + ProjP[sp2,sp3] KBRz[ff]);
109LBH :=            H        bpbar[sp1,cc].dq[sp3,ff,cc] (ProjP[sp1,sp3] KBLh[ff] + ProjM[sp1,sp3] KBRh[ff]);
110
111(* X interactions *)
112LXW := gw/Sqrt[2] W[mu] Ga[mu,sp1,sp2] xbar[sp1,cc].uq[sp3,ff,cc] (ProjM[sp2,sp3] KXL[ff] + ProjP[sp2,sp3] KXR[ff]);
113
114(* Y interactions *)
115LYW := gw/Sqrt[2] Wbar[mu] Ga[mu,sp1,sp2] ybar[sp1,cc].dq[sp3,ff,cc] (ProjM[sp2,sp3] KYL[ff] + ProjP[sp2,sp3] KYR[ff]);
116
117LVLQ := LTW + LTZ + LTH + LBW + LBZ + LBH + LXW + LYW + HC[LTW + LTZ + LTH + LBW + LBZ + LBH + LXW + LYW] + LTK + LBK + LXK + LYK + LSM ;