(***************** Model file for the tutorial at MC4BSM***********) M$ModelName = "MC4BSM"; M$Information = {Authors -> {"C. Degrande"}, Institutions -> {"UCL"}, Emails -> {"celine.degrande@uclouvain.be"}, Date -> "15. 04. 2010", References -> {"G. F. Giudice, C. Grojean, A. Pomarol and R. Rattazzi, JHEP 0706 2007) 045 [arXiv:hep-ph/0703164]"}, URLs -> {"http://mc4bsm.nbi.dk"}, Version -> "1.0" }; (***** Parameter list ******) M$Parameters = { (***** External parameters *****) MZ == { ParameterType -> External, TeX -> MZ, Value -> 91.188, Description -> "mass of the Z"}, mrho == { ParameterType -> External, TeX -> Subscript[m,\[Rho]], Value -> 800, Description -> "mass of the Z prime"}, grho == { ParameterType -> External, TeX -> Subscript[g,\[Rho]], Value -> 1, InteractionOrder -> {NP,1}, Description -> "coupling constant of the Z prime"}, (***** Internal parameters *****) MZ1 == { ParameterType -> Internal, Value -> Sqrt[1/2 (mrho^2 + MZ^2 - Sqrt[mrho^4 - 2 MZ^2 mrho^2 + MZ^4 + 16 grho^2 MZ^4 ])], Description -> "mass of Z1" }, cz == { ParameterType -> Internal, Value -> 2 grho MZ^2/Sqrt[(MZ^2 - MZ1^2)^2 + (2 grho MZ^2 )^2], Description -> "cosine of the mixing angle between Z and Z prime" }, sz == { ParameterType -> Internal, Value -> Sqrt[1-cz^2], Description -> "sine of the mixing angle between Z and Z prime" }, MZ2 == { ParameterType -> Internal, Value -> Sqrt[mrho^2+MZ^2-MZ1^2], Description -> "mass of Z2" } } (***** Particle classes list ******) M$ClassesDescription = { V[7] == { ClassName -> Zprime, SelfConjugate -> True, Definitions -> {Zprime[mu_] -> -sz Z1[mu] + cz Z2[mu]}, Indices -> {}, Mass -> 0, Unphysical -> True}, V[8] == { ClassName -> Z1, SelfConjugate -> True, Indices -> {}, Mass -> {MZ1, Internal}, Width -> {WZ1, 2.41140351}, PropagatorLabel -> "z1", PropagatorType -> W, PropagatorArrow -> None, FullName -> "Z1" }, V[9] == { ClassName -> Z2, SelfConjugate -> True, Indices -> {}, Mass -> {MZ2, Internal}, Width -> {WZ2, 1}, PropagatorLabel -> "z2", PropagatorType -> W, PropagatorArrow -> None, FullName -> "Z2" }} (***** Lagrangian *****) LNP:=-1/4 FS[Zprime,mu,nu] FS[Zprime,mu,nu] + 1/2 mrho^2 Zprime[mu]^2 + 1/2 grho v^2 Zprime[mu] (gw Wi[mu,3]-g1 B[mu])*Sqrt[gw^2+g1^2];