| 1 | (****** This is the FeynRules mod-file for the extension of the SM *********)
|
|---|
| 2 | (* where a new quark Q is added *)
|
|---|
| 3 | (* This model file needs the Standard Model model file SM.fr. *)
|
|---|
| 4 | (* *)
|
|---|
| 5 | (* We use unitary gauge (default in SM.fr) *)
|
|---|
| 6 | (* *)
|
|---|
| 7 | (****************************************************************************)
|
|---|
| 8 |
|
|---|
| 9 |
|
|---|
| 10 | M$ModelName = "topFCNC_JLOperators";
|
|---|
| 11 |
|
|---|
| 12 | (* We work in the unitary gauge (just to make sure we do it again) *)
|
|---|
| 13 | FeynmanGauge = False;
|
|---|
| 14 |
|
|---|
| 15 | M$Information = {Authors -> "A.Amorim",
|
|---|
| 16 | Institutions ->
|
|---|
| 17 | "University of Porto & LIP",
|
|---|
| 18 | Emails -> "artur@lip.pt",
|
|---|
| 19 | Date -> "25.September.2015",
|
|---|
| 20 | Version -> "1.0"};
|
|---|
| 21 |
|
|---|
| 22 | (* ************************** *)
|
|---|
| 23 | (* ***** Parameters ***** *)
|
|---|
| 24 | (* ************************** *)
|
|---|
| 25 |
|
|---|
| 26 | M$Parameters = {
|
|---|
| 27 |
|
|---|
| 28 | (* External parameters *)
|
|---|
| 29 |
|
|---|
| 30 | ReXLut == { TeX -> Subscript[ReX,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 1, InteractionOrder -> {NP,1} },
|
|---|
| 31 | ImXLut == { TeX -> Subscript[ImX,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 2, InteractionOrder -> {NP,1} },
|
|---|
| 32 | ReXRut == { TeX -> Subscript[ReX,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 3, InteractionOrder -> {NP,1} },
|
|---|
| 33 | ImXRut == { TeX -> Subscript[ImX,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 4, InteractionOrder -> {NP,1} },
|
|---|
| 34 | ReXLct == { TeX -> Subscript[ReX,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 5, InteractionOrder -> {NP,1} },
|
|---|
| 35 | ImXLct == { TeX -> Subscript[ImX,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 6, InteractionOrder -> {NP,1} },
|
|---|
| 36 | ReXRct == { TeX -> Subscript[ReX,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 7, InteractionOrder -> {NP,1} },
|
|---|
| 37 | ImXRct == { TeX -> Subscript[ImX,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 8, InteractionOrder -> {NP,1} },
|
|---|
| 38 | ReKLut == { TeX -> Subscript[ReK,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 9, InteractionOrder -> {NP,1} },
|
|---|
| 39 | ImKLut == { TeX -> Subscript[ImK,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 10, InteractionOrder -> {NP,1} },
|
|---|
| 40 | ReKRut == { TeX -> Subscript[ReK,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 11, InteractionOrder -> {NP,1} },
|
|---|
| 41 | ImKRut == { TeX -> Subscript[ImK,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 12, InteractionOrder -> {NP,1} },
|
|---|
| 42 | ReKLct == { TeX -> Subscript[ReK,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 13, InteractionOrder -> {NP,1} },
|
|---|
| 43 | ImKLct == { TeX -> Subscript[ImK,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 14, InteractionOrder -> {NP,1} },
|
|---|
| 44 | ReKRct == { TeX -> Subscript[ReK,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 15, InteractionOrder -> {NP,1} },
|
|---|
| 45 | ImKRct == { TeX -> Subscript[ImK,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 16, InteractionOrder -> {NP,1} },
|
|---|
| 46 | ReZetaLut == { TeX -> Subscript[ReZeta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 17, InteractionOrder -> {NP,1} },
|
|---|
| 47 | ImZetaLut == { TeX -> Subscript[ImZeta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 18, InteractionOrder -> {NP,1} },
|
|---|
| 48 | ReZetaRut == { TeX -> Subscript[ReZeta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 19, InteractionOrder -> {NP,1} },
|
|---|
| 49 | ImZetaRut == { TeX -> Subscript[ImZeta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 20, InteractionOrder -> {NP,1} },
|
|---|
| 50 | ReZetaLct == { TeX -> Subscript[ReZeta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 21, InteractionOrder -> {NP,1} },
|
|---|
| 51 | ImZetaLct == { TeX -> Subscript[ImZeta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 22, InteractionOrder -> {NP,1} },
|
|---|
| 52 | ReZetaRct == { TeX -> Subscript[ReZeta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 23, InteractionOrder -> {NP,1} },
|
|---|
| 53 | ImZetaRct == { TeX -> Subscript[ImZeta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 24, InteractionOrder -> {NP,1} },
|
|---|
| 54 | ReEtaLut == { TeX -> Subscript[ReEta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 25, InteractionOrder -> {NP,1} },
|
|---|
| 55 | ImEtaLut == { TeX -> Subscript[ImEta,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 26, InteractionOrder -> {NP,1} },
|
|---|
| 56 | ReEtaRut == { TeX -> Subscript[ReEta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 27, InteractionOrder -> {NP,1} },
|
|---|
| 57 | ImEtaRut == { TeX -> Subscript[ImEta,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 28, InteractionOrder -> {NP,1} },
|
|---|
| 58 | ReEtaLct == { TeX -> Subscript[ReEta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 29, InteractionOrder -> {NP,1} },
|
|---|
| 59 | ImEtaLct == { TeX -> Subscript[ImEta,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 30, InteractionOrder -> {NP,1} },
|
|---|
| 60 | ReEtaRct == { TeX -> Subscript[ReEta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 31, InteractionOrder -> {NP,1} },
|
|---|
| 61 | ImEtaRct == { TeX -> Subscript[ImEta,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 32, InteractionOrder -> {NP,1} },
|
|---|
| 62 | ReLambdaLut == { TeX -> Subscript[ReLambda,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 33, InteractionOrder -> {NP,1} },
|
|---|
| 63 | ImLambdaLut == { TeX -> Subscript[ImLambda,Lut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 34, InteractionOrder -> {NP,1} },
|
|---|
| 64 | ReLambdaRut == { TeX -> Subscript[ReLambda,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 35, InteractionOrder -> {NP,1} },
|
|---|
| 65 | ImLambdaRut == { TeX -> Subscript[ImLambda,Rut], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 36, InteractionOrder -> {NP,1} },
|
|---|
| 66 | ReLambdaLct == { TeX -> Subscript[ReLambda,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 37, InteractionOrder -> {NP,1} },
|
|---|
| 67 | ImLambdaLct == { TeX -> Subscript[ImLambda,Lct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 38, InteractionOrder -> {NP,1} },
|
|---|
| 68 | ReLambdaRct == { TeX -> Subscript[ReLambda,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 39, InteractionOrder -> {NP,1} },
|
|---|
| 69 | ImLambdaRct == { TeX -> Subscript[ImLambda,Rct], ParameterType -> External, Value -> 0.0, BlockName -> NEWCOUP, OrderBlock -> 40, InteractionOrder -> {NP,1} }
|
|---|
| 70 |
|
|---|
| 71 |
|
|---|
| 72 |
|
|---|
| 73 |
|
|---|
| 74 | (* Internal Parameters *)
|
|---|
| 75 |
|
|---|
| 76 |
|
|---|
| 77 | };
|
|---|
| 78 |
|
|---|
| 79 | (* ************************** *)
|
|---|
| 80 | (* ***** Lagrangian ***** *)
|
|---|
| 81 | (* ************************** *)
|
|---|
| 82 |
|
|---|
| 83 | (* H to glue glue *)
|
|---|
| 84 |
|
|---|
| 85 | LZFCNC := Block[{la},
|
|---|
| 86 | la = -(gw*(ReXLct+I*ImXLct)/(2*cw))*(cbar.Ga[mu].ProjM.t)*Z[mu]-(gw*(ReXLut+I*ImXLut)/(2*cw))*(ubar.Ga[mu].ProjM.t)*Z[mu]
|
|---|
| 87 | -(gw*(ReXRct+I*ImXRct)/(2*cw))*(cbar.Ga[mu].ProjP.t)*Z[mu]-(gw*(ReXRut+I*ImXRut)/(2*cw))*(ubar.Ga[mu].ProjP.t)*Z[mu]
|
|---|
| 88 | +(gw*(ReKLct+I*ImKLct)*I/(4*cw*MZ))*(cbar.Ga[mu].Ga[nu].ProjM.t)*del[Z[nu],mu] - (gw*(ReKLct+I*ImKLct)*I/(4*cw*MZ))*(cbar.Ga[nu].Ga[mu].ProjM.t)*del[Z[nu],mu]
|
|---|
| 89 | +(gw*(ReKRct+I*ImKRct)*I/(4*cw*MZ))*(cbar.Ga[mu].Ga[nu].ProjP.t)*del[Z[nu],mu] - (gw*(ReKRct+I*ImKRct)*I/(4*cw*MZ))*(cbar.Ga[nu].Ga[mu].ProjP.t)*del[Z[nu],mu]
|
|---|
| 90 | +(gw*(ReKLut+I*ImKLut)*I/(4*cw*MZ))*(ubar.Ga[mu].Ga[nu].ProjM.t)*del[Z[nu],mu] - (gw*(ReKLut+I*ImKLut)*I/(4*cw*MZ))*(ubar.Ga[nu].Ga[mu].ProjM.t)*del[Z[nu],mu]
|
|---|
| 91 | +(gw*(ReKRut+I*ImKRut)*I/(4*cw*MZ))*(ubar.Ga[mu].Ga[nu].ProjP.t)*del[Z[nu],mu] - (gw*(ReKRut+I*ImKRut)*I/(4*cw*MZ))*(ubar.Ga[nu].Ga[mu].ProjP.t)*del[Z[nu],mu];
|
|---|
| 92 | la + HC[la]];
|
|---|
| 93 |
|
|---|
| 94 | LphotonFCNC := Block[{la},
|
|---|
| 95 | la =((ReLambdaLct+I*ImLambdaLct)*I*ee/(2*MT))*(cbar.Ga[mu].Ga[nu].ProjM.t)*del[A[nu],mu]-((ReLambdaLct+I*ImLambdaLct)*I*ee/(2*MT))*(cbar.Ga[nu].Ga[mu].ProjM.t)*del[A[nu],mu]
|
|---|
| 96 | +((ReLambdaRct+I*ImLambdaRct)*I*ee/(2*MT))*(cbar.Ga[mu].Ga[nu].ProjP.t)*del[A[nu],mu]-((ReLambdaRct+I*ImLambdaRct)*I*ee/(2*MT))*(cbar.Ga[nu].Ga[mu].ProjP.t)*del[A[nu],mu]
|
|---|
| 97 | +((ReLambdaLut+I*ImLambdaLut)*I*ee/(2*MT))*(ubar.Ga[mu].Ga[nu].ProjM.t)*del[A[nu],mu]-((ReLambdaLut+I*ImLambdaLut)*I*ee/(2*MT))*(ubar.Ga[nu].Ga[mu].ProjM.t)*del[A[nu],mu]
|
|---|
| 98 | +((ReLambdaRut+I*ImLambdaRut)*I*ee/(2*MT))*(ubar.Ga[mu].Ga[nu].ProjP.t)*del[A[nu],mu]-((ReLambdaRut+I*ImLambdaRut)*I*ee/(2*MT))*(ubar.Ga[nu].Ga[mu].ProjP.t)*del[A[nu],mu];
|
|---|
| 99 | la + HC[la]];
|
|---|
| 100 |
|
|---|
| 101 | LgluonFCNC := Block[{la},
|
|---|
| 102 | la = (gs*(ReZetaLct+I*ImZetaLct)*I/(2*MT))*2*(cbar.T[a].Ga[mu].Ga[nu].ProjM.t)*del[G[nu,a],mu]-(gs*(ReZetaLct+I*ImZetaLct)*I/(2*MT))*2*(cbar.T[a].Ga[nu].Ga[mu].ProjM.t)*del[G[nu,a],mu]
|
|---|
| 103 | +(gs*(ReZetaRct+I*ImZetaRct)*I/(2*MT))*2*(cbar.T[a].Ga[mu].Ga[nu].ProjP.t)*del[G[nu,a],mu]-(gs*(ReZetaRct+I*ImZetaRct)*I/(2*MT))*2*(cbar.T[a].Ga[nu].Ga[mu].ProjP.t)*del[G[nu,a],mu]
|
|---|
| 104 | +(gs*(ReZetaLut+I*ImZetaLut)*I/(2*MT))*2*(ubar.T[a].Ga[mu].Ga[nu].ProjM.t)*del[G[nu,a],mu]-(gs*(ReZetaLut+I*ImZetaLut)*I/(2*MT))*2*(ubar.T[a].Ga[nu].Ga[mu].ProjM.t)*del[G[nu,a],mu]
|
|---|
| 105 | +(gs*(ReZetaRut+I*ImZetaRut)*I/(2*MT))*2*(ubar.T[a].Ga[mu].Ga[nu].ProjP.t)*del[G[nu,a],mu]-(gs*(ReZetaRut+I*ImZetaRut)*I/(2*MT))*2*(ubar.T[a].Ga[nu].Ga[mu].ProjP.t)*del[G[nu,a],mu];
|
|---|
| 106 | la + HC[la]];
|
|---|
| 107 |
|
|---|
| 108 | LHiggsFCNC := Block[{la},
|
|---|
| 109 | la = -((ReEtaLut+I*ImEtaLut)/Sqrt[2])*(ubar.ProjM.t)*H-((ReEtaRut+I*ImEtaRut)/Sqrt[2])*(ubar.ProjP.t)*H
|
|---|
| 110 | -((ReEtaLct+I*ImEtaLct)/Sqrt[2])*(cbar.ProjM.t)*H-((ReEtaRct+I*ImEtaRct)/Sqrt[2])*(cbar.ProjP.t)*H;
|
|---|
| 111 | la + HC[la]];
|
|---|
| 112 |
|
|---|
| 113 |
|
|---|
| 114 | Ltot := LSM + LZFCNC + LgluonFCNC + LHiggsFCNC + LphotonFCNC;
|
|---|
| 115 |
|
|---|
