Triplets: TripletDiquarks.fr

File TripletDiquarks.fr, 3.5 KB (added by Benjamin Fuks, 13 years ago)

Model file

Line 
1
2M$ModelName = "Triplet_Diquarks";
3
4(*
5
6 The convention and notations follow 0909.2666
7 We also allow for non intergeneration couplings between quarks.
8 The mixing matrices are however implemented in general, and put diagonal via the independent
9 restriction file
10
11 MFV.rst
12
13 The new particles are
14
15 trip1 = (6, 1, -1/3)
16
17*)
18
19M$Information = {Authors -> {"C. Duhr"},
20 Version -> "1.0",
21 Date -> "27. 10. 2010",
22 Institutions -> {"IPPP, Durham"},
23 Emails -> {"claude.duhr@durham.ac.uk"}};
24
25
26(* Coupling matrices are symmetric *)
27
28
29SetAttributes[LQQR, Orderless];
30SetAttributes[LUDL, Orderless];
31SetAttributes[LUUL, Orderless];
32SetAttributes[LDDL, Orderless];
33
34M$Parameters = {
35
36 LQQRR == {Indices -> {Index[Generation], Index[Generation]},
37 Value -> {LQQRR[1,1] -> 0.1,
38 LQQRR[2,2] -> 0.1,
39 LQQRR[3,3] -> 0.1,
40 LQQRR[i_, j_] :> 0 /; NumericQ[i] && NumericQ[j] && (i !=j)},
41 InteractionOrder -> {QCD, 1},
42 ParameterType -> External,
43 ComplexParameter -> False
44 },
45
46 LQQRI == {Indices -> {Index[Generation], Index[Generation]},
47 Value -> {LQQRI[_,_] -> 0},
48 InteractionOrder -> {QCD, 1},
49 ParameterType -> External,
50 ComplexParameter -> False
51 },
52
53 LUDLR == {Indices -> {Index[Generation], Index[Generation]},
54 Value -> {LUDLR[1,1] -> 0.1,
55 LUDLR[2,2] -> 0.1,
56 LUDLR[3,3] -> 0.1,
57 LUDLR[i_, j_] :> 0 /; NumericQ[i] && NumericQ[j] && (i !=j)},
58 InteractionOrder -> {QCD, 1},
59 ParameterType -> External,
60 ComplexParameter -> False
61 },
62
63 LUDLI == {Indices -> {Index[Generation], Index[Generation]},
64 Value -> {LUDLI[_,_] -> 0},
65 InteractionOrder -> {QCD, 1},
66 ParameterType -> External,
67 ComplexParameter -> False
68 },
69
70
71
72 LHS1 == {InteractionOrder -> {QED, 2},
73 Value -> 1,
74 ParameterType -> External
75 },
76
77 LSS11 == {InteractionOrder -> {QCD, 2},
78 Value -> 1,
79 ParameterType -> External
80 },
81
82
83(* Internal parameters *)
84
85 LQQR == {Indices -> {Index[Generation], Index[Generation]},
86 Value -> {LQQR[i_,j_] :> LQQRR[i,j] + I LQQRI[i,j]},
87 InteractionOrder -> {QCD, 1},
88 ComplexParameter -> True
89 },
90
91 LUDL == {Indices -> {Index[Generation], Index[Generation]},
92 Value -> {LUDL[i_,j_] :> LUDLR[i,j] + I LUDLI[i,j]},
93 InteractionOrder -> {QCD, 1},
94 ComplexParameter -> True
95 }
96};
97
98M$ClassesDescription = {
99
100 S[100] == {
101 ClassName -> trip1,
102 SelfConjugate -> False,
103 Indices -> {Index[Colour]},
104 Mass -> {Mtrip1, 500},
105 Width -> {Wtrip1, 4.4108},
106 QuantumNumbers -> {Q -> -1/3, Y -> -1/3}
107 }
108};
109
110
111(* the Lagrangian *)
112
113
114LTripKin := DC[trip1bar[k], mu]DC[trip1[k],mu] - Mtrip1^2 trip1bar[k]trip1[k];
115
116LD11 := 2 (Eps[k,i,j] trip1bar[k] LQQR[n,m] ProjP[s,r] dqbar[s,n,i].CC[uq][r,m,j] +
117 Eps[k,i,j] trip1bar[k] LUDL[n,m] ProjM[s,r] dqbar[s,n,i].CC[uq][r,m,j]);
118
119LD1 := LD11 + HC[LD11];
120
121LPot := ExpandIndices[LHS1 Phibar[ii] Phi[ii] trip1bar[k]trip1[k] +
122 LSS11 trip1bar[k1]trip1[k1]trip1bar[k2]trip1[k2], FlavorExpand->{SU2W,SU2D}];
123
124LTrip := LTripKin + LD1 + LPot;
125
126
127
128
129
130