Triplets: TripletDiquarks.fr

File TripletDiquarks.fr, 3.5 KB (added by BenjF, 5 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