| 1 | (* ************************************************************************* *)
|
|---|
| 2 | (* ***** ***** *)
|
|---|
| 3 | (* ***** FeynRules model file: coloron ***** *)
|
|---|
| 4 | (* ***** Authors: M. Baker, J. Fuentes-Martin, M. König, G. Isidori ***** *)
|
|---|
| 5 | (* ***** ***** *)
|
|---|
| 6 | (* ************************************************************************* *)
|
|---|
| 7 |
|
|---|
| 8 | (* ************************** *)
|
|---|
| 9 | (* ***** Information ***** *)
|
|---|
| 10 | (* ************************** *)
|
|---|
| 11 |
|
|---|
| 12 | M$ModelName = "coloron";
|
|---|
| 13 |
|
|---|
| 14 | M$Information = {Authors -> {"M. Baker", "J. Fuentes-Martin", "M. Koenig", "G. Isidori"},
|
|---|
| 15 | Version -> "1.1",
|
|---|
| 16 | Date -> "07.04.2021",
|
|---|
| 17 | Institutions -> {"U. Zurich", "Johannes Gutenberg University Mainz"},
|
|---|
| 18 | Emails -> {"baker@physik.uzh.ch", "jfuentes@uni-mainz.de", "matthias.koenig@tum.de", "isidori@physik.uzh.ch"}
|
|---|
| 19 | };
|
|---|
| 20 |
|
|---|
| 21 | (* ************************** *)
|
|---|
| 22 | (* *** Interaction orders *** *)
|
|---|
| 23 | (* ************************** *)
|
|---|
| 24 |
|
|---|
| 25 | M$InteractionOrderHierarchy = { {QCD,1}, {NP,1}, {QED,2} };
|
|---|
| 26 |
|
|---|
| 27 |
|
|---|
| 28 | (* ************************** *)
|
|---|
| 29 | (* **** Particle classes **** *)
|
|---|
| 30 | (* ************************** *)
|
|---|
| 31 |
|
|---|
| 32 | M$ClassesDescription = {
|
|---|
| 33 | V[7] == {
|
|---|
| 34 | ClassName -> Gp,
|
|---|
| 35 | SelfConjugate -> True,
|
|---|
| 36 | Indices -> {Index[Gluon]},
|
|---|
| 37 | Mass -> {MGp, 4000},
|
|---|
| 38 | Width -> {WGp, 800},
|
|---|
| 39 | ParticleName -> "Gp",
|
|---|
| 40 | PDG -> 44,
|
|---|
| 41 | PropagatorLabel -> "Gp",
|
|---|
| 42 | PropagatorType -> C,
|
|---|
| 43 | PropagatorArrow -> None,
|
|---|
| 44 | FullName -> "Color-octect (coloron)"
|
|---|
| 45 | }
|
|---|
| 46 | };
|
|---|
| 47 |
|
|---|
| 48 | (* ************************** *)
|
|---|
| 49 | (* ***** Parameters ***** *)
|
|---|
| 50 | (* ************************** *)
|
|---|
| 51 |
|
|---|
| 52 | M$Parameters = {
|
|---|
| 53 | (* External Parameters *)
|
|---|
| 54 | gGp == {
|
|---|
| 55 | ParameterType -> External,
|
|---|
| 56 | Value -> 3.0,
|
|---|
| 57 | InteractionOrder -> {NP,1},
|
|---|
| 58 | TeX -> Subscript[g,Gp],
|
|---|
| 59 | BlockName -> NPGPCOUP,
|
|---|
| 60 | OrderBlock -> 1,
|
|---|
| 61 | Description -> "Overall Gp coupling strength"
|
|---|
| 62 | },
|
|---|
| 63 |
|
|---|
| 64 | kappaq33 == {
|
|---|
| 65 | ParameterType -> External,
|
|---|
| 66 | Value -> 1.0,
|
|---|
| 67 | InteractionOrder -> {NP,0},
|
|---|
| 68 | TeX -> Subscript[Superscript[\[Kappa],L],33],
|
|---|
| 69 | BlockName -> NPGPCOUP,
|
|---|
| 70 | OrderBlock -> 2,
|
|---|
| 71 | Description -> "Gp LH coupling third generation"
|
|---|
| 72 | },
|
|---|
| 73 |
|
|---|
| 74 | kappaRu33 == {
|
|---|
| 75 | ParameterType -> External,
|
|---|
| 76 | Value -> 1.0,
|
|---|
| 77 | InteractionOrder -> {NP,0},
|
|---|
| 78 | TeX -> Subscript[Superscript[\[Kappa],Ru],33],
|
|---|
| 79 | BlockName -> NPGPCOUP,
|
|---|
| 80 | OrderBlock -> 3,
|
|---|
| 81 | Description -> "Gp RH tt coupling"
|
|---|
| 82 | },
|
|---|
| 83 |
|
|---|
| 84 | kappaRd33 == {
|
|---|
| 85 | ParameterType -> External,
|
|---|
| 86 | Value -> 1.0,
|
|---|
| 87 | InteractionOrder -> {NP,0},
|
|---|
| 88 | TeX -> Subscript[Superscript[\[Kappa],Rd],33],
|
|---|
| 89 | BlockName -> NPGPCOUP,
|
|---|
| 90 | OrderBlock -> 4,
|
|---|
| 91 | Description -> "Gp RH bb coupling"
|
|---|
| 92 | },
|
|---|
| 93 |
|
|---|
| 94 | kappaqll == {
|
|---|
| 95 | ParameterType -> External,
|
|---|
| 96 | Value -> 0.0,
|
|---|
| 97 | InteractionOrder -> {NP,0},
|
|---|
| 98 | TeX -> Subscript[Superscript[\[Kappa],L],ll],
|
|---|
| 99 | BlockName -> NPGPCOUP,
|
|---|
| 100 | OrderBlock -> 5,
|
|---|
| 101 | Description -> "Gp LH coupling light-generations"
|
|---|
| 102 | },
|
|---|
| 103 |
|
|---|
| 104 | kappaRull == {
|
|---|
| 105 | ParameterType -> External,
|
|---|
| 106 | Value -> 0.0,
|
|---|
| 107 | InteractionOrder -> {NP,0},
|
|---|
| 108 | TeX -> Subscript[Superscript[\[Kappa],Ru],ll],
|
|---|
| 109 | BlockName -> NPGPCOUP,
|
|---|
| 110 | OrderBlock -> 6,
|
|---|
| 111 | Description -> "Gp RH up-type light-quark coupling"
|
|---|
| 112 | },
|
|---|
| 113 |
|
|---|
| 114 | kappaRdll == {
|
|---|
| 115 | ParameterType -> External,
|
|---|
| 116 | Value -> 0.0,
|
|---|
| 117 | InteractionOrder -> {NP,0},
|
|---|
| 118 | TeX -> Subscript[Superscript[\[Kappa],Rd],ll],
|
|---|
| 119 | BlockName -> NPGPCOUP,
|
|---|
| 120 | OrderBlock -> 7,
|
|---|
| 121 | Description -> "Gp RH light down-quark coupling"
|
|---|
| 122 | },
|
|---|
| 123 |
|
|---|
| 124 | kappaG1 == {
|
|---|
| 125 | ParameterType -> External,
|
|---|
| 126 | Value -> 0.0,
|
|---|
| 127 | InteractionOrder -> {NP,0},
|
|---|
| 128 | TeX -> Subscript[\[Kappa],G2],
|
|---|
| 129 | BlockName -> NPGPCOUP,
|
|---|
| 130 | OrderBlock -> 8,
|
|---|
| 131 | Description -> "Non-minimal Gp coupling to gluons 1"
|
|---|
| 132 | },
|
|---|
| 133 |
|
|---|
| 134 | kappaG2 == {
|
|---|
| 135 | ParameterType -> External,
|
|---|
| 136 | Value -> 0.0,
|
|---|
| 137 | InteractionOrder -> {NP,0},
|
|---|
| 138 | TeX -> Subscript[\[Kappa],G1],
|
|---|
| 139 | BlockName -> NPGPCOUP,
|
|---|
| 140 | OrderBlock -> 9,
|
|---|
| 141 | Description -> "Non-minimal Gp coupling to gluons 2"
|
|---|
| 142 | },
|
|---|
| 143 |
|
|---|
| 144 | (* Internal Parameters *)
|
|---|
| 145 | kappaL == {
|
|---|
| 146 | ParameterType -> Internal,
|
|---|
| 147 | Indices -> {Index[Generation], Index[Generation]},
|
|---|
| 148 | Hermitian -> True,
|
|---|
| 149 | Definitions -> {kappaL[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
|
|---|
| 150 | Value -> {kappaL[1,1] -> kappaqll, kappaL[2,2] -> kappaqll, kappaL[3,3] -> kappaq33},
|
|---|
| 151 | InteractionOrder -> {NP,0},
|
|---|
| 152 | TeX -> Subscript[\[Kappa],L],
|
|---|
| 153 | Description -> "Gp LH coupling matrix"
|
|---|
| 154 | },
|
|---|
| 155 |
|
|---|
| 156 | kappaRu == {
|
|---|
| 157 | ParameterType -> Internal,
|
|---|
| 158 | Indices -> {Index[Generation], Index[Generation]},
|
|---|
| 159 | Hermitian -> True,
|
|---|
| 160 | Definitions -> {kappaRu[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
|
|---|
| 161 | Value -> {kappaRu[1,1] -> kappaRull, kappaRu[2,2] -> kappaRull, kappaRu[3,3] -> kappaRu33},
|
|---|
| 162 | InteractionOrder -> {NP,0},
|
|---|
| 163 | TeX -> Subscript[\[Kappa],Ru],
|
|---|
| 164 | Description -> "Gp RH up-type coupling matrix"
|
|---|
| 165 | },
|
|---|
| 166 |
|
|---|
| 167 | kappaRd == {
|
|---|
| 168 | ParameterType -> Internal,
|
|---|
| 169 | Indices -> {Index[Generation], Index[Generation]},
|
|---|
| 170 | Hermitian -> True,
|
|---|
| 171 | Definitions -> {kappaRd[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
|
|---|
| 172 | Value -> {kappaRd[1,1] -> kappaRdll, kappaRd[2,2] -> kappaRdll, kappaRd[3,3] -> kappaRd33},
|
|---|
| 173 | InteractionOrder -> {NP,0},
|
|---|
| 174 | TeX -> Subscript[\[Kappa],Rd],
|
|---|
| 175 | Description -> "Gp RH down-type coupling matrix"
|
|---|
| 176 | }
|
|---|
| 177 | };
|
|---|
| 178 |
|
|---|
| 179 |
|
|---|
| 180 |
|
|---|
| 181 | (* ************************** *)
|
|---|
| 182 | (* ***** Lagrangian ***** *)
|
|---|
| 183 | (* ************************** *)
|
|---|
| 184 |
|
|---|
| 185 | LGpKin := Block[{mu,nu,aa1}, - 1/4 * (DC[Gp[mu,aa1], nu] - DC[Gp[nu,aa1], mu]) * (DC[Gp[mu,aa1], nu] - DC[Gp[nu,aa1], mu]) ];
|
|---|
| 186 |
|
|---|
| 187 | LGpF := Block[{ff1,ff2,s1,s2,s3,aa1,cc1,cc2,mu},
|
|---|
| 188 | gGp * Gp[mu,aa1] *
|
|---|
| 189 | (
|
|---|
| 190 | kappaL[ff1,ff2] * newCKM[ff3,ff1] * uqbar[s1,ff3,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * Conjugate[newCKM[ff4,ff2]] * uq[s3,ff4,cc2] +
|
|---|
| 191 | kappaL[ff1,ff2] * dqbar[s1,ff1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * dq[s3,ff2,cc2] +
|
|---|
| 192 | kappaRu[ff1,ff2] * uqbar[s1,ff1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjP[s2,s3] * uq[s3,ff2,cc2] +
|
|---|
| 193 | kappaRd[ff1,ff2] * dqbar[s1,ff1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjP[s2,s3] * dq[s3,ff2,cc2]
|
|---|
| 194 | )
|
|---|
| 195 | ];
|
|---|
| 196 |
|
|---|
| 197 | LGpG := Block[{aa1,aa2,aa3,mu,nu},
|
|---|
| 198 | 1/2 * kappaG1 * (DC[Gp[nu,aa1], mu] - DC[Gp[mu,aa1], nu]) * FS[G,mu,nu,aa1] +
|
|---|
| 199 | gs * kappaG2 * f[aa1,aa2,aa3] * Gp[mu,aa1] * Gp[nu,aa2] * FS[G,mu,nu,aa3]
|
|---|
| 200 | ];
|
|---|
| 201 |
|
|---|
| 202 | LGp := LGpKin + LGpF + LGpG;
|
|---|