NLOModels: VLL_v1.fr

File VLL_v1.fr, 9.3 KB (added by Benjamin Fuks, 3 years ago)
Line 
1(********************************************************************************)
2(****** FeynRules mod-file for Vector-like Leptons *)
3(****** *)
4(****** Authors: B. Fuks *)
5(********************************************************************************)
6
7M$ModelName = "VLL";
8
9M$Information = {
10 Authors -> {"B. Fuks"},
11 Institutions -> {"LPTHE / Sorbonne U."},
12 Emails -> {"fuks@lpthe.jussieu.fr"},
13 Version -> "1.0",
14 Date -> "13.10.2022"
15};
16
17FeynmanGauge = False;
18M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {VLL,1}};
19M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {VLL,99}};
20
21(* ************************** *)
22(* ***** Fields ***** *)
23(* ************************** *)
24
25M$ClassesDescription = {
26 (* VLL doublet [gauge eigenstates] *)
27 F[21] == { ClassName -> VLL0, Unphysical->True, SelfConjugate->False, Indices -> {Index[SU2D]}, FlavorIndex -> SU2D, QuantumNumbers -> {Y ->-1/2},
28 Definitions -> { VLL0[sp1_,1] :> VLLN[sp1], VLL0[sp1_,2] :> VLLE[sp1] } },
29
30 (* VLL doublet [mass eigenstates] *)
31 F[31] == { ClassName -> VLLE, SelfConjugate->False, Mass->{MVLLE,601}, Width -> {WVLLE,5.}, QuantumNumbers -> {Q->-1, LeptonNumber->1}, PDG -> 9000001, FullName -> "VLL doublet - charged component"},
32 F[32] == { ClassName -> VLLN, SelfConjugate->False, Mass->{MVLLN,602}, Width -> {WVLLN,5.}, QuantumNumbers -> {LeptonNumber->1}, PDG -> 9000002, FullName -> "VLL doublet - neutral component"},
33
34 (* VLL singlets [gauge eigenstates] *)
35 F[22] == { ClassName -> VLE0, Unphysical->True, SelfConjugate->False, QuantumNumbers -> {Y->-1}, Definitions->{VLE0[sp_]->VLE[sp]} },
36 F[23] == { ClassName -> VLN0, Unphysical->True, SelfConjugate->False, Definitions->{VLN0[sp_]->VLN[sp]} },
37
38 (* VLL singlets [mass eigenstates] *)
39 F[33] == { ClassName -> VLE, SelfConjugate->False, Mass->{MVLE,603}, Width -> {WVLE,5.}, QuantumNumbers -> {Q->-1, LeptonNumber->1}, PDG -> 9000003, FullName -> "Charged VLL singlet"},
40 F[34] == { ClassName -> VLN, SelfConjugate->False, Mass->{MVLN,604}, Width -> {WVLN,5.}, QuantumNumbers -> {LeptonNumber->1}, PDG -> 9000004, FullName -> "Neutral VLL singlet"}
41};
42
43(* ************************** *)
44(* ***** Parameters ***** *)
45(* ************************** *)
46
47M$Parameters = {
48
49 (* Higgs boson interactions *)
50 KLLEH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
51 Description->"Doublet charged VLL - Higgs coupling (left-handed)" ,
52 BlockName -> KLLEH, Value -> {KLLEH[1] -> 0.001, KLLEH[2] -> 0.002, KLLEH[3] -> 0.003} },
53 KRLEH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
54 Description->"Doublet charged VLL - Higgs coupling (right-handed)" ,
55 BlockName -> KRLEH, Value -> {KRLEH[1] -> 0.004, KRLEH[2] -> 0.005, KRLEH[3] -> 0.006} },
56 KLEH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
57 Description->"Singlet charged VLL - Higgs coupling (left-handed)" ,
58 BlockName -> KLEH, Value -> {KLEH[1] -> 0.007, KLEH[2] -> 0.008, KLEH[3] -> 0.009} },
59 KREH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
60 Description->"Singlet charged VLL - Higgs coupling (right-handed)" ,
61 BlockName -> KREH, Value -> {KREH[1] -> 0.010, KREH[2] -> 0.011, KREH[3] -> 0.012} },
62 KLLNH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
63 Description->"Doublet neutral VLL - Higgs coupling (left-handed)" ,
64 BlockName -> KLLNH, Value -> {KLLNH[1] -> 0.013, KLLNH[2] -> 0.014, KLLNH[3] -> 0.015} },
65 KLNH == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
66 Description->"Singlet neutral VLL - Higgs coupling (left-handed)" ,
67 BlockName -> KLNH, Value -> {KLNH[1] -> 0.016, KLNH[2] -> 0.017, KLNH[3] -> 0.018} },
68
69 (* Z boson interactions *)
70 KLLEZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
71 Description->"Doublet charged VLL - Z coupling (left-handed)" ,
72 BlockName -> KLLEZ, Value -> {KLLEZ[1] -> 0.021, KLLEZ[2] -> 0.022, KLLEZ[3] -> 0.023} },
73 KRLEZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
74 Description->"Doublet charged VLL - Z coupling (right-handed)" ,
75 BlockName -> KRLEZ, Value -> {KRLEZ[1] -> 0.024, KRLEZ[2] -> 0.025, KRLEZ[3] -> 0.026} },
76 KLEZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
77 Description->"Singlet charged VLL - Z coupling (left-handed)" ,
78 BlockName -> KLEZ, Value -> {KLEZ[1] -> 0.027, KLEZ[2] -> 0.028, KLEZ[3] -> 0.029} },
79 KREZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
80 Description->"Singlet charged VLL - Z coupling (right-handed)" ,
81 BlockName -> KREZ, Value -> {KREZ[1] -> 0.030, KREZ[2] -> 0.031, KREZ[3] -> 0.032} },
82 KLLNZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
83 Description->"Doublet neutral VLL - Z coupling (left-handed)" ,
84 BlockName -> KLLNZ, Value -> {KLLNZ[1] -> 0.033, KLLNZ[2] -> 0.034, KLLNZ[3] -> 0.035} },
85 KLNZ == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
86 Description->"Singlet neutral VLL - Z coupling (left-handed)" ,
87 BlockName -> KLNZ, Value -> {KLNZ[1] -> 0.036, KLNZ[2] -> 0.037, KLNZ[3] -> 0.038} },
88
89 (* W boson interactions *)
90 KLLEW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
91 Description->"Doublet charged VLL - W coupling (left-handed)" ,
92 BlockName -> KLLEW, Value -> {KLLEW[1] -> 0.041, KLLEW[2] -> 0.042, KLLEW[3] -> 0.043} },
93 KLEW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
94 Description->"Singlet charged VLL - W coupling (left-handed)" ,
95 BlockName -> KLEW, Value -> {KLEW[1] -> 0.044, KLEW[2] -> 0.045, KLEW[3] -> 0.046} },
96 KLLNW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
97 Description->"Doublet neutral VLL - W coupling (left-handed)" ,
98 BlockName -> KLLNW, Value -> {KLLNW[1] -> 0.047, KLLNW[2] -> 0.048, KLLNW[3] -> 0.049} },
99 KRLNW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
100 Description->"Doublet neutral VLL - W coupling (right-handed)" ,
101 BlockName -> KRLNW, Value -> {KRLNW[1] -> 0.050, KRLNW[2] -> 0.051, KRLNW[3] -> 0.052} },
102 KLNW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
103 Description->"Singlet neutral VLL - W coupling (left-handed)" ,
104 BlockName -> KLNW, Value -> {KLNW[1] -> 0.053, KLNW[2] -> 0.054, KLNW[3] -> 0.055} },
105 KRNW == { ParameterType->External, ComplexParameter->False, Indices->{Index[Generation]}, InteractionOrder->{VLL,1},
106 Description->"Doublet neutral VLL - W coupling (right-handed)" ,
107 BlockName -> KRNW, Value -> {KRNW[1] -> 0.056, KRNW[2] -> 0.057, KRNW[3] -> 0.058} }
108};
109
110
111(* ************************** *)
112(* ***** Lagrangian ***** *)
113(* ************************** *)
114
115(* Kinetic and mass terms *)
116Lkin := I VLL0bar.Ga[mu].DC[VLL0,mu] - MVLLE VLLEbar.VLLE - MVLLN VLLNbar.VLLN \
117 + I VLE0bar.Ga[mu].DC[VLE0,mu] - MVLE VLEbar.VLE \
118 + I VLN0bar.Ga[mu].DC[VLN0,mu] - MVLN VLNbar.VLN;
119
120(* Interactions with the Higgs *)
121LVLH := H VLLEbar[sp1] . l[sp3,ff] (ProjM[sp1,sp3] KLLEH[ff] + ProjP[sp1,sp3] KRLEH[ff]) \
122 + H VLEbar[sp1] . l[sp3,ff] (ProjM[sp1,sp3] KLEH[ff] + ProjP[sp1,sp3] KREH[ff]) \
123 + H VLLNbar[sp1] . vl[sp3,ff] (ProjM[sp1,sp3] KLLNH[ff] ) \
124 + H VLNbar[sp1] . vl[sp3,ff] (ProjM[sp1,sp3] KLNH[ff] );
125
126(* Interactions with the Z boson*)
127LVLZ := gw/(2 cw) * Z[mu] * Ga[mu,sp1,sp2] * VLLEbar[sp1] . l[sp3,ff] (ProjM[sp2,sp3] KLLEZ[ff] + ProjP[sp2,sp3] KRLEZ[ff]) \
128 + gw/(2 cw) * Z[mu] * Ga[mu,sp1,sp2] * VLEbar[sp1] . l[sp3,ff] (ProjM[sp2,sp3] KLEZ[ff] + ProjP[sp2,sp3] KREZ[ff]) \
129 + gw/(2 cw) * Z[mu] * Ga[mu,sp1,sp2] * VLLNbar[sp1] . vl[sp3,ff] (ProjM[sp2,sp3] KLLNZ[ff] ) \
130 + gw/(2 cw) * Z[mu] * Ga[mu,sp1,sp2] * VLNbar[sp1] . vl[sp3,ff] (ProjM[sp2,sp3] KLNZ[ff] );
131
132(* Interactions with the W boson*)
133LVLW := gw/Sqrt[2] * Wbar[mu] * Ga[mu,sp1,sp2] * VLLEbar[sp1] . vl[sp3,ff] (ProjM[sp2,sp3] KLLEW[ff] ) \
134 + gw/Sqrt[2] * Wbar[mu] * Ga[mu,sp1,sp2] * VLEbar[sp1] . vl[sp3,ff] (ProjM[sp2,sp3] KLEW[ff] ) \
135 + gw/Sqrt[2] * W[mu] * Ga[mu,sp1,sp2] * VLLNbar[sp1] . l[sp3,ff] (ProjM[sp2,sp3] KLLNW[ff] + ProjP[sp2,sp3] KRLNW[ff] ) \
136 + gw/Sqrt[2] * W[mu] * Ga[mu,sp1,sp2] * VLNbar[sp1] . l[sp3,ff] (ProjM[sp2,sp3] KLNW[ff] + ProjP[sp2,sp3] KRNW[ff] );
137
138LVLL := { LSM, Lkin, LVLH, LVLZ, LVLW, HC[LVLH], HC[LVLZ], HC[LVLW] };