# This file was automatically created by FeynRules 2.4.55
# Mathematica version: 10.1.0 for Mac OS X x86 (64-bit) (March 24, 2015)
# Date: Tue 31 Jan 2017 12:52:07
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
# This file was automatically created by FeynRules 2.1
# Mathematica version: 9.0 for Mac OS X x86 (64-bit) (January 24, 2013)
# Date: Thu 16 Oct 2014 15:27:30
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.b,P.b__tilde__):'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__,P.ta__plus__):'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t,P.t__tilde__):'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_sc = Decay(name = 'Decay_sc',
particle = P.sc,
partial_widths = {(P.a,P.a):'(Msc**2*((ee**4*gu**2*Msc**4)/(18.*cmath.pi**4*vev**2) + (ee**4*gu**2*Msc**16)/(5.6448e6*cmath.pi**4*MT**12*vev**2) + (ee**4*gu**2*Msc**14)/(453600.*cmath.pi**4*MT**10*vev**2) + (191*ee**4*gu**2*Msc**12)/(8.1648e6*cmath.pi**4*MT**8*vev**2) + (41*ee**4*gu**2*Msc**10)/(136080.*cmath.pi**4*MT**6*vev**2) + (7*ee**4*gu**2*Msc**8)/(4320.*cmath.pi**4*MT**4*vev**2) + (ee**4*gu**2*Msc**6)/(108.*cmath.pi**4*MT**2*vev**2)))/(32.*cmath.pi*abs(Msc)**3)',
(P.b,P.b__tilde__):'(3*gd**2*MB**2*Msc**2*cmath.sqrt(-4*MB**2*Msc**2 + Msc**4))/(4.*cmath.pi*vev**2*abs(Msc)**3)',
(P.chi,P.chi__tilde__):'((-8*gchi**2*Mchi**2 + 2*gchi**2*Msc**2)*cmath.sqrt(-4*Mchi**2*Msc**2 + Msc**4))/(16.*cmath.pi*abs(Msc)**3)',
(P.g,P.g):'(Msc**2*((G**4*gu**2*Msc**4)/(16.*cmath.pi**4*vev**2) + (G**4*gu**2*Msc**16)/(5.0176e6*cmath.pi**4*MT**12*vev**2) + (G**4*gu**2*Msc**14)/(403200.*cmath.pi**4*MT**10*vev**2) + (191*G**4*gu**2*Msc**12)/(7.2576e6*cmath.pi**4*MT**8*vev**2) + (41*G**4*gu**2*Msc**10)/(120960.*cmath.pi**4*MT**6*vev**2) + (7*G**4*gu**2*Msc**8)/(3840.*cmath.pi**4*MT**4*vev**2) + (G**4*gu**2*Msc**6)/(96.*cmath.pi**4*MT**2*vev**2)))/(32.*cmath.pi*abs(Msc)**3)',
(P.ta__minus__,P.ta__plus__):'(gl**2*Msc**2*MTA**2*cmath.sqrt(Msc**4 - 4*Msc**2*MTA**2))/(4.*cmath.pi*vev**2*abs(Msc)**3)',
(P.t,P.t__tilde__):'(3*gu**2*Msc**2*MT**2*cmath.sqrt(Msc**4 - 4*Msc**2*MT**2))/(4.*cmath.pi*vev**2*abs(Msc)**3)'})
Decay_t = Decay(name = 'Decay_t',
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 10.0 for Linux x86 (64-bit) (September 9, 2014)
# Date: Mon 28 Aug 2017 23:20:25
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_N1 = Decay(
name='Decay_N1',
particle=P.N1,
partial_widths={
(P.WR__plus__, P.e__minus__):
'((mN1**2 - MWR**2)*((ee**2*kRl**2*mN1**2)/(2.*sw**2) + (ee**2*kRl**2*mN1**4)/(2.*MWR**2*sw**2) - (ee**2*kRl**2*MWR**2)/sw**2))/(32.*cmath.pi*abs(mN1)**3)',
(P.WR__minus__, P.e__plus__):
'((mN1**2 - MWR**2)*((ee**2*kRl**2*mN1**2)/(2.*sw**2) + (ee**2*kRl**2*mN1**4)/(2.*MWR**2*sw**2) - (ee**2*kRl**2*MWR**2)/sw**2))/(32.*cmath.pi*abs(mN1)**3)'
})
Decay_N2 = Decay(
name='Decay_N2',
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.0.0 for Mac OS X x86 (64-bit) (July 28, 2016)
# Date: Mon 28 Aug 2017 15:47:16
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_Z = Decay(
name='Decay_Z',
particle=P.Z,
partial_widths={
(P.W__minus__, P.W__plus__):
'(((-12*cw**2*ee**2*MW**2)/sw**2 - (17*cw**2*ee**2*MZ**2)/sw**2 + (4*cw**2*ee**2*MZ**4)/(MW**2*sw**2) + (cw**2*ee**2*MZ**6)/(4.*MW**4*sw**2))*cmath.sqrt(-4*MW**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
(P.u, P.u__tilde__):
'(MZ**2*(-(ee**2*MZ**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
# This file was automatically created by FeynRules 2.0.30
# Mathematica version: 9.0 for Linux x86 (64-bit) (January 24, 2013)
# Date: Fri 14 Nov 2014 11:30:01
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.W__minus__,P.c):'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__,P.t):'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__,P.u):'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_c = Decay(name = 'Decay_c',
particle = P.c,
partial_widths = {(P.W__plus__,P.b):'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__,P.d):'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__,P.s):'(((3*CKM2x2*ee**2*MC**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MS**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MC**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MC**2*MS**2*complexconjugate(CKM2x2))/(MW**2*sw**2) + (3*CKM2x2*ee**2*MS**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MW**2*complexconjugate(CKM2x2))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MS**2 + MS**4 - 2*MC**2*MW**2 - 2*MS**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)'})
Decay_d = Decay(name = 'Decay_d',
particle = P.d,
partial_widths = {(P.W__minus__,P.c):'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MD)**3)',
(P.W__minus__,P.t):'(((3*CKM3x1*ee**2*MD**2*complexconjugate(CKM3x1))/(2.*sw**2) + (3*CKM3x1*ee**2*MT**2*complexconjugate(CKM3x1))/(2.*sw**2) + (3*CKM3x1*ee**2*MD**4*complexconjugate(CKM3x1))/(2.*MW**2*sw**2) - (3*CKM3x1*ee**2*MD**2*MT**2*complexconjugate(CKM3x1))/(MW**2*sw**2) + (3*CKM3x1*ee**2*MT**4*complexconjugate(CKM3x1))/(2.*MW**2*sw**2) - (3*CKM3x1*ee**2*MW**2*complexconjugate(CKM3x1))/sw**2)*cmath.sqrt(MD**4 - 2*MD**2*MT**2 + MT**4 - 2*MD**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MD)**3)',
(P.W__minus__,P.u):'(((3*CKM1x1*ee**2*MD**2*complexconjugate(CKM1x1))/(2.*sw**2) + (3*CKM1x1*ee**2*MU**2*complexconjugate(CKM1x1))/(2.*sw**2) + (3*CKM1x1*ee**2*MD**4*complexconjugate(CKM1x1))/(2.*MW**2*sw**2) - (3*CKM1x1*ee**2*MD**2*MU**2*complexconjugate(CKM1x1))/(MW**2*sw**2) + (3*CKM1x1*ee**2*MU**4*complexconjugate(CKM1x1))/(2.*MW**2*sw**2) - (3*CKM1x1*ee**2*MW**2*complexconjugate(CKM1x1))/sw**2)*cmath.sqrt(MD**4 - 2*MD**2*MU**2 + MU**4 - 2*MD**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MD)**3)'})
Decay_e__minus__ = Decay(name = 'Decay_e__minus__',
particle = P.e__minus__,
partial_widths = {(P.W__minus__,P.ve):'((Me**2 - MW**2)*((ee**2*Me**2)/(2.*sw**2) + (ee**2*Me**4)/(2.*MW**2*sw**2) - (ee**2*MW**2)/sw**2))/(32.*cmath.pi*abs(Me)**3)'})
# This file was automatically created by FeynRules 2.3.24
# Mathematica version: 12.0.0 for Mac OS X x86 (64-bit) (April 7, 2019)
# Date: Tue 5 May 2020 11:41:20
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.DP__minus__, P.t):
'(((24*FermionEpsilon**2*FermionY**4*MB**2*msbarDownQuark3**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (24*FermionEpsilon**2*FermionY**4*MDPc**2*msbarDownQuark3**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*MB**2*msbarUpQuark3**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (24*FermionEpsilon**2*FermionY**4*MDPc**2*msbarUpQuark3**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (96*FermionEpsilon**2*FermionY**4*MB*msbarDownQuark3*msbarUpQuark3*MT)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*msbarDownQuark3**2*MT**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*msbarUpQuark3**2*MT**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2))*cmath.sqrt(MB**4 - 2*MB**2*MDPc**2 + MDPc**4 - 2*MB**2*MT**2 - 2*MDPc**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
particle=P.c,
partial_widths={
(P.DP__plus__, P.s):
'(((24*FermionEpsilon**2*FermionY**4*MC**2*msbarDownQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (24*FermionEpsilon**2*FermionY**4*MDPc**2*msbarDownQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*MS**2*msbarDownQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (96*FermionEpsilon**2*FermionY**4*MC*MS*msbarDownQuark2*msbarUpQuark2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*MC**2*msbarUpQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) - (24*FermionEpsilon**2*FermionY**4*MDPc**2*msbarUpQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2) + (24*FermionEpsilon**2*FermionY**4*MS**2*msbarUpQuark2**2)/(FermionA**2*(2 + FermionA)**2*FermionEta**6*cmath.pi**2*PionMass**2))*cmath.sqrt(MC**4 - 2*MC**2*MDPc**2 + MDPc**4 - 2*MC**2*MS**2 - 2*MDPc**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.b):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.d):
'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.a, P.a):
'(MH**2*((2*eAASM**2*kAA**2*MH**4)/vF**2 + (2*eAASM**2*lAACP**2*MH**4)/vF**2))/(32.*cmath.pi*abs(MH)**3)',
(P.g, P.g):
'(eggSM**2*MH**6)/(2.*cmath.pi*vF**2*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((2*eWW**2*MH**4)/vF**2 + (2*eWWCP**2*MH**4)/vF**2 + (kWW**2*MH**4)/vF**2 - (8*eWW**2*MH**2*MW**2)/vF**2 - (8*eWWCP**2*MH**2*MW**2)/vF**2 + (12*eWW*kWW*MH**2*MW**2)/vF**2 - (4*kWW**2*MH**2*MW**2)/vF**2 + (12*eWW**2*MW**4)/vF**2 - (24*eWW*kWW*MW**4)/vF**2 + (12*kWW**2*MW**4)/vF**2 + (2*kWW*MH**4*noise)/vF**2 + (12*eWW*MH**2*MW**2*noise)/vF**2 - (8*kWW*MH**2*MW**2*noise)/vF**2 - (24*eWW*MW**4*noise)/vF**2 + (24*kWW*MW**4*noise)/vF**2 + (MH**4*noise**2)/vF**2 - (4*MH**2*MW**2*noise**2)/vF**2 + (12*MW**4*noise**2)/vF**2)*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.a, P.Z):
'((MH**2 - MZ**2)*((2*eZASM**2*kZA**2*MH**4)/vF**2 + (2*eZASM**2*lZACP**2*MH**4)/vF**2 - (4*eZASM**2*kZA**2*MH**2*MZ**2)/vF**2 - (4*eZASM**2*lZACP**2*MH**2*MZ**2)/vF**2 + (2*eZASM**2*kZA**2*MZ**4)/vF**2 + (2*eZASM**2*lZACP**2*MZ**4)/vF**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((2*eZZ**2*MH**4)/vF**2 + (2*eZZCP**2*MH**4)/vF**2 + (kZZ**2*MH**4)/vF**2 - (8*eZZ**2*MH**2*MZ**2)/vF**2 - (8*eZZCP**2*MH**2*MZ**2)/vF**2 + (12*eZZ*kZZ*MH**2*MZ**2)/vF**2 - (4*kZZ**2*MH**2*MZ**2)/vF**2 + (12*eZZ**2*MZ**4)/vF**2 - (24*eZZ*kZZ*MZ**4)/vF**2 + (12*kZZ**2*MZ**4)/vF**2 + (2*kZZ*MH**4*noise)/vF**2 + (12*eZZ*MH**2*MZ**2*noise)/vF**2 - (8*kZZ*MH**2*MZ**2*noise)/vF**2 - (24*eZZ*MZ**4*noise)/vF**2 + (24*kZZ*MZ**4*noise)/vF**2 + (MH**4*noise**2)/vF**2 - (4*MH**2*MZ**2*noise**2)/vF**2 + (12*MZ**4*noise**2)/vF**2)*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)',
(P.b, P.b__tilde__):
'((-12*kb**2*MB**2*ybeff**2 + 3*kb**2*MH**2*ybeff**2 + 3*lb**2*MH**2*ybeff**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.c, P.c__tilde__):
'((-12*kc**2*MC**2*yceff**2 + 3*kc**2*MH**2*yceff**2 + 3*lc**2*MH**2*yceff**2)*cmath.sqrt(-4*MC**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.mu__minus__, P.mu__plus__):
'(MH**2*(kmu**2*MH**2*ymueff**2 + lmu**2*MH**2*ymueff**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((ktau**2*MH**2*ytaueff**2 + ltau**2*MH**2*ytaueff**2 - 4*ktau**2*MTA**2*ytaueff**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)'
})
Decay_Z = Decay(
name='Decay_Z',
particle=P.Z,
# This file was automatically created by FeynRules 2.3.36
# Mathematica version: 12.0.0 for Mac OS X x86 (64-bit) (April 7, 2019)
# Date: Sat 28 Mar 2020 01:43:25
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
particle=P.c,
partial_widths={
(P.W__plus__, P.b):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.d):
'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.s):
'(((3*CKM2x2*ee**2*MC**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MS**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MC**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MC**2*MS**2*complexconjugate(CKM2x2))/(MW**2*sw**2) + (3*CKM2x2*ee**2*MS**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MW**2*complexconjugate(CKM2x2))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MS**2 + MS**4 - 2*MC**2*MW**2 - 2*MS**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)'
})
# This file was automatically created by FeynRules 2.3.13
# Mathematica version: 11.0.1 for Mac OS X x86 (64-bit) (September 21, 2016)
# Date: Wed 23 May 2018 22:20:22
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.SS1, P.SS2):
'(gHSS**2*cmath.sqrt(MH**4 - 2*MH**2*MSS1**2 + MSS1**4 - 2*MH**2*MSS2**2 - 2*MSS1**2*MSS2**2 + MSS2**4))/(16.*cmath.pi*abs(MH)**3)',
(P.g, P.g):
'(GH**2*MH**6)/(8.*cmath.pi*abs(MH)**3)',
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_SS1 = Decay(
name='Decay_SS1',
particle=P.SS1,
partial_widths={
(P.H, P.SS2):
'(gHSS**2*cmath.sqrt(MH**4 - 2*MH**2*MSS1**2 + MSS1**4 - 2*MH**2*MSS2**2 - 2*MSS1**2*MSS2**2 + MSS2**4))/(16.*cmath.pi*abs(MSS1)**3)',
# This file was automatically created by FeynRules 2.3.22
# Mathematica version: 10.4.1 for Linux x86 (64-bit) (April 11, 2016)
# Date: Thu 8 Sep 2016 01:36:08
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
particle=P.c,
partial_widths={
(P.W__plus__, P.b):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.d):
'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.s):
'(((3*CKM2x2*ee**2*MC**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MS**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MC**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MC**2*MS**2*complexconjugate(CKM2x2))/(MW**2*sw**2) + (3*CKM2x2*ee**2*MS**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MW**2*complexconjugate(CKM2x2))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MS**2 + MS**4 - 2*MC**2*MW**2 - 2*MS**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)'
})
# This file was automatically created by FeynRules 2.1.48
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
# Date: Wed 2 Apr 2014 11:19:40
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.h__minus__, P.t):
'((3*I3a33*MB**2*complexconjugate(I3a33) - 3*I3a33*mhc**2*complexconjugate(I3a33) - 6*I4a33*MB*MT*complexconjugate(I3a33) + 3*I3a33*MT**2*complexconjugate(I3a33) + 3*I4a33*MB**2*complexconjugate(I4a33) - 3*I4a33*mhc**2*complexconjugate(I4a33) - 6*I3a33*MB*MT*complexconjugate(I4a33) + 3*I4a33*MT**2*complexconjugate(I4a33))*cmath.sqrt(MB**4 - 2*MB**2*mhc**2 + mhc**4 - 2*MB**2*MT**2 - 2*mhc**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*gw**2*MB**2)/2. + (3*gw**2*MT**2)/2. + (3*gw**2*MB**4)/(2.*MW**2) - (3*gw**2*MB**2*MT**2)/MW**2 + (3*gw**2*MT**4)/(2.*MW**2) - 3*gw**2*MW**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_h1 = Decay(
name='Decay_h1',
particle=P.h1,
partial_widths={
(P.a, P.a):
'(mh1**2*(cw**4*g1**4*TH1x1**2*vev**2 - 4*cw**3*g1**3*gw*sw*TH1x1**2*vev**2 + 6*cw**2*g1**2*gw**2*sw**2*TH1x1**2*vev**2 - 4*cw*g1*gw**3*sw**3*TH1x1**2*vev**2 + gw**4*sw**4*TH1x1**2*vev**2))/(32.*cmath.pi*abs(mh1)**3)',
(P.a, P.h3):
'((mh1**2 - mh3**2)*(-(cw**2*g1**2*mh1**2*TH2x1**2*TH3x3**2)/2. - (cw**2*g1**2*mh3**2*TH2x1**2*TH3x3**2)/2. + cw*g1*gw*mh1**2*sw*TH2x1**2*TH3x3**2 + cw*g1*gw*mh3**2*sw*TH2x1**2*TH3x3**2 - (gw**2*mh1**2*sw**2*TH2x1**2*TH3x3**2)/2. - (gw**2*mh3**2*sw**2*TH2x1**2*TH3x3**2)/2.))/(16.*cmath.pi*abs(mh1)**3)',
(P.a, P.Z):
'((mh1**2 - MZ**2)*((3*cw**4*g1**2*gw**2*TH1x1**2*vev**2)/4. + (3*cw**3*g1**3*gw*sw*TH1x1**2*vev**2)/2. - (3*cw**3*g1*gw**3*sw*TH1x1**2*vev**2)/2. + (3*cw**2*g1**4*sw**2*TH1x1**2*vev**2)/4. - 3*cw**2*g1**2*gw**2*sw**2*TH1x1**2*vev**2 + (3*cw**2*gw**4*sw**2*TH1x1**2*vev**2)/4. - (3*cw*g1**3*gw*sw**3*TH1x1**2*vev**2)/2. + (3*cw*g1*gw**3*sw**3*TH1x1**2*vev**2)/2. + (3*g1**2*gw**2*sw**4*TH1x1**2*vev**2)/4.))/(16.*cmath.pi*abs(mh1)**3)',
(P.b, P.b__tilde__):
'((-3*GD3x3**2*MB**2*TH2x1**2 - 12*GD3x3*MB**2*TH1x1*TH2x1*yb + 3*GD3x3*mh1**2*TH1x1*TH2x1*yb - 12*MB**2*TH1x1**2*yb**2 + 3*mh1**2*TH1x1**2*yb**2 - 6*GD3x3*MB**2*TH2x1**2*complexconjugate(GD3x3) + 3*GD3x3*mh1**2*TH2x1**2*complexconjugate(GD3x3) - 12*MB**2*TH1x1*TH2x1*yb*complexconjugate(GD3x3) + 3*mh1**2*TH1x1*TH2x1*yb*complexconjugate(GD3x3) - 3*MB**2*TH2x1**2*complexconjugate(GD3x3)**2)*cmath.sqrt(-4*MB**2*mh1**2 + mh1**4))/(16.*cmath.pi*abs(mh1)**3)',
(P.h2, P.h2):
'((36*l1**2*TH1x1**2*TH1x2**4*vev**2 + 18*l1*l6*TH1x1*TH1x2**4*TH2x1*vev**2 + 36*l1*l6*TH1x1**2*TH1x2**3*TH2x2*vev**2 + 24*l1*l3*TH1x1*TH1x2**3*TH2x1*TH2x2*vev**2 + 24*l1*l4*TH1x1*TH1x2**3*TH2x1*TH2x2*vev**2 + 24*l1*l5*TH1x1*TH1x2**3*TH2x1*TH2x2*vev**2 + 6*l3*l6*TH1x2**3*TH2x1**2*TH2x2*vev**2 + 6*l4*l6*TH1x2**3*TH2x1**2*TH2x2*vev**2 + 12*l1*l3*TH1x1**2*TH1x2**2*TH2x2**2*vev**2 + 12*l1*l4*TH1x1**2*TH1x2**2*TH2x2**2*vev**2 + 12*l1*l5*TH1x1**2*TH1x2**2*TH2x2**2*vev**2 + 15*l3*l6*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2 + 15*l4*l6*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2 + 18*l1*l7*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2 + 4*l3**2*TH1x2**2*TH2x1**2*TH2x2**2*vev**2 + 8*l3*l4*TH1x2**2*TH2x1**2*TH2x2**2*vev**2 + 4*l4**2*TH1x2**2*TH2x1**2*TH2x2**2*vev**2 + 8*l3*l5*TH1x2**2*TH2x1**2*TH2x2**2*vev**2 + 8*l4*l5*TH1x2**2*TH2x1**2*TH2x2**2*vev**2 + 6*l3*l6*TH1x1**2*TH1x2*TH2x2**3*vev**2 + 6*l4*l6*TH1x1**2*TH1x2*TH2x2**3*vev**2 + 4*l3**2*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2 + 8*l3*l4*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2 + 4*l4**2*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2 + 8*l3*l5*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2 + 8*l4*l5*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2 + 6*l3*l7*TH1x2*TH2x1**2*TH2x2**3*vev**2 + 6*l4*l7*TH1x2*TH2x1**2*TH2x2**3*vev**2 + l3**2*TH1x1**2*TH2x2**4*vev**2 + 2*l3*l4*TH1x1**2*TH2x2**4*vev**2 + l4**2*TH1x1**2*TH2x2**4*vev**2 + 2*l3*l5*TH1x1**2*TH2x2**4*vev**2 + 2*l4*l5*TH1x1**2*TH2x2**4*vev**2 + 3*l3*l7*TH1x1*TH2x1*TH2x2**4*vev**2 + 3*l4*l7*TH1x1*TH2x1*TH2x2**4*vev**2 + 24*l1*TH1x1*TH1x2**3*TH2x1*TH2x2*vev**2*complexconjugate(l5) + 12*l6*TH1x2**3*TH2x1**2*TH2x2*vev**2*complexconjugate(l5) + 12*l1*TH1x1**2*TH1x2**2*TH2x2**2*vev**2*complexconjugate(l5) + 30*l6*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2*complexconjugate(l5) + 8*l3*TH1x2**2*TH2x1**2*TH2x2**2*vev**2*complexconjugate(l5) + 8*l4*TH1x2**2*TH2x1**2*TH2x2**2*vev**2*complexconjugate(l5) + 16*l5*TH1x2**2*TH2x1**2*TH2x2**2*vev**2*complexconjugate(l5) + 12*l6*TH1x1**2*TH1x2*TH2x2**3*vev**2*complexconjugate(l5) + 8*l3*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2*complexconjugate(l5) + 8*l4*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2*complexconjugate(l5) + 16*l5*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2*complexconjugate(l5) + 12*l7*TH1x2*TH2x1**2*TH2x2**3*vev**2*complexconjugate(l5) + 2*l3*TH1x1**2*TH2x2**4*vev**2*complexconjugate(l5) + 2*l4*TH1x1**2*TH2x2**4*vev**2*complexconjugate(l5) + 4*l5*TH1x1**2*TH2x2**4*vev**2*complexconjugate(l5) + 6*l7*TH1x1*TH2x1*TH2x2**4*vev**2*complexconjugate(l5) + 18*l1*TH1x1*TH1x2**4*TH2x1*vev**2*complexconjugate(l6) + 9*l6*TH1x2**4*TH2x1**2*vev**2*complexconjugate(l6) + 36*l1*TH1x1**2*TH1x2**3*TH2x2*vev**2*complexconjugate(l6) + 36*l6*TH1x1*TH1x2**3*TH2x1*TH2x2*vev**2*complexconjugate(l6) + 6*l3*TH1x2**3*TH2x1**2*TH2x2*vev**2*complexconjugate(l6) + 6*l4*TH1x2**3*TH2x1**2*TH2x2*vev**2*complexconjugate(l6) + 12*l5*TH1x2**3*TH2x1**2*TH2x2*vev**2*complexconjugate(l6) + 36*l6*TH1x1**2*TH1x2**2*TH2x2**2*vev**2*complexconjugate(l6) + 15*l3*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2*complexconjugate(l6) + 15*l4*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2*complexconjugate(l6) + 30*l5*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2*complexconjugate(l6) + 9*l7*TH1x2**2*TH2x1**2*TH2x2**2*vev**2*complexconjugate(l6) + 6*l3*TH1x1**2*TH1x2*TH2x2**3*vev**2*complexconjugate(l6) + 6*l4*TH1x1**2*TH1x2*TH2x2**3*vev**2*complexconjugate(l6) + 12*l5*TH1x1**2*TH1x2*TH2x2**3*vev**2*complexconjugate(l6) + 18*l7*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2*complexconjugate(l6) + 18*l1*TH1x1*TH1x2**2*TH2x1*TH2x2**2*vev**2*complexconjugate(l7) + 9*l6*TH1x2**2*TH2x1**2*TH2x2**2*vev**2*complexconjugate(l7) + 18*l6*TH1x1*TH1x2*TH2x1*TH2x2**3*vev**2*complexconjugate(l7) + 6*l3*TH1x2*TH2x1**2*TH2x2**3*vev**2*complexconjugate(l7) + 6*l4*TH1x2*TH2x1**2*TH2x2**3*vev**2*complexconjugate(l7) + 12*l5*TH1x2*TH2x1**2*TH2x2**3*vev**2*complexconjugate(l7) + 3*l3*TH1x1*TH2x1*TH2x2**4*vev**2*complexconjugate(l7) + 3*l4*TH1x1*TH2x1*TH2x2**4*vev**2*complexconjugate(l7) + 6*l5*TH1x1*TH2x1*TH2x2**4*vev**2*complexconjugate(l7) + 9*l7*TH2x1**2*TH2x2**4*vev**2*complexconjugate(l7))*cmath.sqrt(mh1**4 - 4*mh1**2*mh2**2))/(32.*cmath.pi*abs(mh1)**3)',
# This file was automatically created by FeynRules 2.3.2
# Mathematica version: 12.1.0 for Mac OS X x86 (64-bit) (March 14, 2020)
# Date: Fri 28 May 2021 15:01:09
from object_library import all_decays, Decay
import particles as P
Decay_S0 = Decay(
name='Decay_S0',
particle=P.S0,
partial_widths={
(P.g, P.g):
'(2*gSg**2*MS0**6)/(cmath.pi*Lambda**2*abs(MS0)**3)',
(P.Vq, P.Vq__tilde__):
'((2*gSQv**2*MS0**2 - 8*gSQv**2*MVq**2)*cmath.sqrt(MS0**4 - 4*MS0**2*MVq**2))/(16.*cmath.pi*abs(MS0)**3)'
})
Decay_A0 = Decay(
name='Decay_A0',
particle=P.A0,
partial_widths={
(P.g, P.g):
'(2*gPg**2*MA0**6)/(cmath.pi*Lambda**2*abs(MA0)**3)',
(P.Vq, P.Vq__tilde__):
'(gPQv**2*MA0**2*cmath.sqrt(MA0**4 - 4*MA0**2*MVq**2))/(8.*cmath.pi*abs(MA0)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
# This file was automatically created by FeynRules 2.0.6
# Mathematica version: 8.0 for Linux x86 (64-bit) (October 10, 2011)
# Date: Tue 7 Apr 2015 19:45:59
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.X1__tilde__,P.d__tilde__):'((MB**2 - MX1**2)*(6*Lam11x1x3*lam1Norm**2*MB**2*complexconjugate(Lam11x1x3) - 6*Lam11x3x1*lam1Norm**2*MB**2*complexconjugate(Lam11x1x3) - 6*Lam11x1x3*lam1Norm**2*MX1**2*complexconjugate(Lam11x1x3) + 6*Lam11x3x1*lam1Norm**2*MX1**2*complexconjugate(Lam11x1x3) - 6*Lam11x1x3*lam1Norm**2*MB**2*complexconjugate(Lam11x3x1) + 6*Lam11x3x1*lam1Norm**2*MB**2*complexconjugate(Lam11x3x1) + 6*Lam11x1x3*lam1Norm**2*MX1**2*complexconjugate(Lam11x3x1) - 6*Lam11x3x1*lam1Norm**2*MX1**2*complexconjugate(Lam11x3x1)))/(96.*cmath.pi*abs(MB)**3)',
(P.X1__tilde__,P.s__tilde__):'((MB**2 - MX1**2)*(6*Lam11x2x3*lam1Norm**2*MB**2*complexconjugate(Lam11x2x3) - 6*Lam11x3x2*lam1Norm**2*MB**2*complexconjugate(Lam11x2x3) - 6*Lam11x2x3*lam1Norm**2*MX1**2*complexconjugate(Lam11x2x3) + 6*Lam11x3x2*lam1Norm**2*MX1**2*complexconjugate(Lam11x2x3) - 6*Lam11x2x3*lam1Norm**2*MB**2*complexconjugate(Lam11x3x2) + 6*Lam11x3x2*lam1Norm**2*MB**2*complexconjugate(Lam11x3x2) + 6*Lam11x2x3*lam1Norm**2*MX1**2*complexconjugate(Lam11x3x2) - 6*Lam11x3x2*lam1Norm**2*MX1**2*complexconjugate(Lam11x3x2)))/(96.*cmath.pi*abs(MB)**3)',
(P.X2__tilde__,P.d__tilde__):'((MB**2 - MX2**2)*(6*Lam12x1x3*lam1Norm**2*MB**2*complexconjugate(Lam12x1x3) - 6*Lam12x3x1*lam1Norm**2*MB**2*complexconjugate(Lam12x1x3) - 6*Lam12x1x3*lam1Norm**2*MX2**2*complexconjugate(Lam12x1x3) + 6*Lam12x3x1*lam1Norm**2*MX2**2*complexconjugate(Lam12x1x3) - 6*Lam12x1x3*lam1Norm**2*MB**2*complexconjugate(Lam12x3x1) + 6*Lam12x3x1*lam1Norm**2*MB**2*complexconjugate(Lam12x3x1) + 6*Lam12x1x3*lam1Norm**2*MX2**2*complexconjugate(Lam12x3x1) - 6*Lam12x3x1*lam1Norm**2*MX2**2*complexconjugate(Lam12x3x1)))/(96.*cmath.pi*abs(MB)**3)',
(P.X2__tilde__,P.s__tilde__):'((MB**2 - MX2**2)*(6*Lam12x2x3*lam1Norm**2*MB**2*complexconjugate(Lam12x2x3) - 6*Lam12x3x2*lam1Norm**2*MB**2*complexconjugate(Lam12x2x3) - 6*Lam12x2x3*lam1Norm**2*MX2**2*complexconjugate(Lam12x2x3) + 6*Lam12x3x2*lam1Norm**2*MX2**2*complexconjugate(Lam12x2x3) - 6*Lam12x2x3*lam1Norm**2*MB**2*complexconjugate(Lam12x3x2) + 6*Lam12x3x2*lam1Norm**2*MB**2*complexconjugate(Lam12x3x2) + 6*Lam12x2x3*lam1Norm**2*MX2**2*complexconjugate(Lam12x3x2) - 6*Lam12x3x2*lam1Norm**2*MX2**2*complexconjugate(Lam12x3x2)))/(96.*cmath.pi*abs(MB)**3)'})
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.b,P.b__tilde__):'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t,P.t__tilde__):'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__,P.ta__plus__):'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_n = Decay(name = 'Decay_n',
particle = P.n,
partial_widths = {(P.X1,P.c__tilde__):'((Mn**2 - MX1**2)*(3*Lam21x2*lam2Left**2*Mn**2*complexconjugate(Lam21x2) + 3*Lam21x2*lam2Norm**2*Mn**2*complexconjugate(Lam21x2) - 3*Lam21x2*lam2Left**2*MX1**2*complexconjugate(Lam21x2) - 3*Lam21x2*lam2Norm**2*MX1**2*complexconjugate(Lam21x2)))/(32.*cmath.pi*abs(Mn)**3)',
(P.X1,P.t__tilde__):'((3*Lam21x3*lam2Left**2*Mn**2*complexconjugate(Lam21x3) + 3*Lam21x3*lam2Norm**2*Mn**2*complexconjugate(Lam21x3) + 12*Lam21x3*lam2Left*lam2Norm*Mn*MT*complexconjugate(Lam21x3) + 3*Lam21x3*lam2Left**2*MT**2*complexconjugate(Lam21x3) + 3*Lam21x3*lam2Norm**2*MT**2*complexconjugate(Lam21x3) - 3*Lam21x3*lam2Left**2*MX1**2*complexconjugate(Lam21x3) - 3*Lam21x3*lam2Norm**2*MX1**2*complexconjugate(Lam21x3))*cmath.sqrt(Mn**4 - 2*Mn**2*MT**2 + MT**4 - 2*Mn**2*MX1**2 - 2*MT**2*MX1**2 + MX1**4))/(32.*cmath.pi*abs(Mn)**3)',
(P.X1,P.u__tilde__):'((Mn**2 - MX1**2)*(3*Lam21x1*lam2Left**2*Mn**2*complexconjugate(Lam21x1) + 3*Lam21x1*lam2Norm**2*Mn**2*complexconjugate(Lam21x1) - 3*Lam21x1*lam2Left**2*MX1**2*complexconjugate(Lam21x1) - 3*Lam21x1*lam2Norm**2*MX1**2*complexconjugate(Lam21x1)))/(32.*cmath.pi*abs(Mn)**3)',
(P.X1__tilde__,P.c):'((Mn**2 - MX1**2)*(3*Lam21x2*lam2Left**2*Mn**2*complexconjugate(Lam21x2) + 3*Lam21x2*lam2Norm**2*Mn**2*complexconjugate(Lam21x2) - 3*Lam21x2*lam2Left**2*MX1**2*complexconjugate(Lam21x2) - 3*Lam21x2*lam2Norm**2*MX1**2*complexconjugate(Lam21x2)))/(32.*cmath.pi*abs(Mn)**3)',
# This file was automatically created by FeynRules 2.3.26
# Mathematica version: 10.3.1 for Mac OS X x86 (64-bit) (December 9, 2015)
# Date: Sun 3 Sep 2017 23:30:31
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.b,P.b__tilde__):'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__,P.ta__plus__):'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t,P.t__tilde__):'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_t = Decay(name = 'Decay_t',
particle = P.t,
partial_widths = {(P.W__plus__,P.b):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MT)**3)'})
Decay_ta__minus__ = Decay(name = 'Decay_ta__minus__',
particle = P.ta__minus__,
partial_widths = {(P.W__minus__,P.vt):'((MTA**2 - MW**2)*((ee**2*MTA**2)/(2.*sw**2) + (ee**2*MTA**4)/(2.*MW**2*sw**2) - (ee**2*MW**2)/sw**2))/(32.*cmath.pi*abs(MTA)**3)'})
Decay_W__plus__ = Decay(name = 'Decay_W__plus__',
particle = P.W__plus__,
# This file was automatically created by FeynRules 2.3.28
# Mathematica version: 10.2.0 for Mac OS X x86 (64-bit) (July 7, 2015)
# Date: Tue 7 Nov 2017 10:51:55
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.b,P.b__tilde__):'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__,P.ta__plus__):'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t,P.t__tilde__):'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_t = Decay(name = 'Decay_t',
particle = P.t,
partial_widths = {(P.W__plus__,P.b):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MT)**3)'})
Decay_ta__minus__ = Decay(name = 'Decay_ta__minus__',
particle = P.ta__minus__,
partial_widths = {(P.W__minus__,P.vt):'((MTA**2 - MW**2)*((ee**2*MTA**2)/(2.*sw**2) + (ee**2*MTA**4)/(2.*MW**2*sw**2) - (ee**2*MW**2)/sw**2))/(32.*cmath.pi*abs(MTA)**3)'})
Decay_W__plus__ = Decay(name = 'Decay_W__plus__',
particle = P.W__plus__,
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.2.0 for Mac OS X x86 (64-bit) (September 11, 2017)
# Date: Tue 12 Mar 2019 11:12:13
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.Xr, P.F3dL3):
'((3*gF3L3x3*MB**2*complexconjugate(gF3L3x3) + 3*gF3L3x3*MF3dL3**2*complexconjugate(gF3L3x3) - 3*gF3L3x3*MXr**2*complexconjugate(gF3L3x3))*cmath.sqrt(MB**4 - 2*MB**2*MF3dL3**2 + MF3dL3**4 - 2*MB**2*MXr**2 - 2*MF3dL3**2*MXr**2 + MXr**4))/(96.*cmath.pi*abs(MB)**3)',
(P.Xr, P.F3dR3):
'((3*gF3dR3x3*MB**2*complexconjugate(gF3dR3x3) + 3*gF3dR3x3*MF3dR3**2*complexconjugate(gF3dR3x3) - 3*gF3dR3x3*MXr**2*complexconjugate(gF3dR3x3))*cmath.sqrt(MB**4 - 2*MB**2*MF3dR3**2 + MF3dR3**4 - 2*MB**2*MXr**2 - 2*MF3dR3**2*MXr**2 + MXr**4))/(96.*cmath.pi*abs(MB)**3)',
(P.Xv, P.F3dL3):
'((3*gF3L3x3*MB**2*complexconjugate(gF3L3x3) + 3*gF3L3x3*MF3dL3**2*complexconjugate(gF3L3x3) + (3*gF3L3x3*MB**4*complexconjugate(gF3L3x3))/MXv**2 - (6*gF3L3x3*MB**2*MF3dL3**2*complexconjugate(gF3L3x3))/MXv**2 + (3*gF3L3x3*MF3dL3**4*complexconjugate(gF3L3x3))/MXv**2 - 6*gF3L3x3*MXv**2*complexconjugate(gF3L3x3))*cmath.sqrt(MB**4 - 2*MB**2*MF3dL3**2 + MF3dL3**4 - 2*MB**2*MXv**2 - 2*MF3dL3**2*MXv**2 + MXv**4))/(96.*cmath.pi*abs(MB)**3)',
(P.Xv, P.F3dR3):
'((3*gF3dR3x3*MB**2*complexconjugate(gF3dR3x3) + 3*gF3dR3x3*MF3dR3**2*complexconjugate(gF3dR3x3) + (3*gF3dR3x3*MB**4*complexconjugate(gF3dR3x3))/MXv**2 - (6*gF3dR3x3*MB**2*MF3dR3**2*complexconjugate(gF3dR3x3))/MXv**2 + (3*gF3dR3x3*MF3dR3**4*complexconjugate(gF3dR3x3))/MXv**2 - 6*gF3dR3x3*MXv**2*complexconjugate(gF3dR3x3))*cmath.sqrt(MB**4 - 2*MB**2*MF3dR3**2 + MF3dR3**4 - 2*MB**2*MXv**2 - 2*MF3dR3**2*MXv**2 + MXv**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_F3dL1 = Decay(
name='Decay_F3dL1',
particle=P.F3dL1,
partial_widths={
(P.W__minus__, P.F3uL1):
'(((3*ee**2*MF3dL1**2)/sw**2 - (18*ee**2*MF3dL1*MF3uL1)/sw**2 + (3*ee**2*MF3uL1**2)/sw**2 + (3*ee**2*MF3dL1**4)/(MW**2*sw**2) - (6*ee**2*MF3dL1**2*MF3uL1**2)/(MW**2*sw**2) + (3*ee**2*MF3uL1**4)/(MW**2*sw**2) - (6*ee**2*MW**2)/sw**2)*cmath.sqrt(MF3dL1**4 - 2*MF3dL1**2*MF3uL1**2 + MF3uL1**4 - 2*MF3dL1**2*MW**2 - 2*MF3uL1**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MF3dL1)**3)',
(P.Xr, P.d):
'((MF3dL1**2 - MXr**2)*(3*gF3L1x1*MF3dL1**2*complexconjugate(gF3L1x1) - 3*gF3L1x1*MXr**2*complexconjugate(gF3L1x1)))/(96.*cmath.pi*abs(MF3dL1)**3)',
# This file was automatically created by FeynRules 2.3.22
# Mathematica version: 11.2.0 for Mac OS X x86 (64-bit) (September 11, 2017)
# Date: Thu 19 Jul 2018 13:38:02
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_Z = Decay(
name='Decay_Z',
particle=P.Z,
partial_widths={
(P.W__minus__, P.W__plus__):
'(((-12*cw**2*ee**2*MW**2)/sw**2 - (17*cw**2*ee**2*MZ**2)/sw**2 + (4*cw**2*ee**2*MZ**4)/(MW**2*sw**2) + (cw**2*ee**2*MZ**6)/(4.*MW**4*sw**2))*cmath.sqrt(-4*MW**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
(P.u, P.u__tilde__):
'(MZ**2*(-(ee**2*MZ**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
(P.c, P.c__tilde__):
'(MZ**2*(-(ee**2*MZ**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
# This file was automatically created by FeynRules 2.1
# Mathematica version: 9.0 for Linux x86 (64-bit) (February 7, 2013)
# Date: Wed 7 Jan 2015 14:05:50
from object_library import all_decays, Decay
import particles as P
Decay_A0 = Decay(
name='Decay_A0',
particle=P.A0,
partial_widths={
(P.g, P.g):
'(G**4*gA0tt**2*ggA0*MA0**6*complexconjugate(ggA0))/(512.*cmath.pi**5*vev**2*abs(MA0)**3)',
(P.t, P.t__tilde__):
'(3*gA0tt**2*MA0**2*MT**2*cmath.sqrt(MA0**4 - 4*MA0**2*MT**2))/(8.*cmath.pi*vev**2*abs(MA0)**3)'
})
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
# This file was automatically created by FeynRules 2.4.1
# Mathematica version: 10.1.0 for Mac OS X x86 (64-bit) (March 24, 2015)
# Date: Thu 18 Jun 2015 16:26:34
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2 + (3*C3phiq*ee**2*MB**2*vev**2)/(2.*Lambda**2*sw**2) + (3*C3phiq*ee**2*MT**2*vev**2)/(2.*Lambda**2*sw**2) + (3*C3phiq*ee**2*MB**4*vev**2)/(2.*Lambda**2*MW**2*sw**2) - (3*C3phiq*ee**2*MB**2*MT**2*vev**2)/(Lambda**2*MW**2*sw**2) + (3*C3phiq*ee**2*MT**4*vev**2)/(2.*Lambda**2*MW**2*sw**2) - (3*C3phiq*ee**2*MW**2*vev**2)/(Lambda**2*sw**2) + (3*ee**2*MB**2*vev**2*complexconjugate(C3phiq))/(2.*Lambda**2*sw**2) + (3*ee**2*MT**2*vev**2*complexconjugate(C3phiq))/(2.*Lambda**2*sw**2) + (3*ee**2*MB**4*vev**2*complexconjugate(C3phiq))/(2.*Lambda**2*MW**2*sw**2) - (3*ee**2*MB**2*MT**2*vev**2*complexconjugate(C3phiq))/(Lambda**2*MW**2*sw**2) + (3*ee**2*MT**4*vev**2*complexconjugate(C3phiq))/(2.*Lambda**2*MW**2*sw**2) - (3*ee**2*MW**2*vev**2*complexconjugate(C3phiq))/(Lambda**2*sw**2) + (3*C3phiq*ee**2*MB**2*vev**4*complexconjugate(C3phiq))/(2.*Lambda**4*sw**2) + (3*C3phiq*ee**2*MT**2*vev**4*complexconjugate(C3phiq))/(2.*Lambda**4*sw**2) + (3*C3phiq*ee**2*MB**4*vev**4*complexconjugate(C3phiq))/(2.*Lambda**4*MW**2*sw**2) - (3*C3phiq*ee**2*MB**2*MT**2*vev**4*complexconjugate(C3phiq))/(Lambda**4*MW**2*sw**2) + (3*C3phiq*ee**2*MT**4*vev**4*complexconjugate(C3phiq))/(2.*Lambda**4*MW**2*sw**2) - (3*C3phiq*ee**2*MW**2*vev**4*complexconjugate(C3phiq))/(Lambda**4*sw**2) + (24*CtW*MB**4*vev**2*complexconjugate(CtW))/Lambda**4 - (48*CtW*MB**2*MT**2*vev**2*complexconjugate(CtW))/Lambda**4 + (24*CtW*MT**4*vev**2*complexconjugate(CtW))/Lambda**4 - (12*CtW*MB**2*MW**2*vev**2*complexconjugate(CtW))/Lambda**4 - (12*CtW*MT**2*MW**2*vev**2*complexconjugate(CtW))/Lambda**4 - (12*CtW*MW**4*vev**2*complexconjugate(CtW))/Lambda**4 - (9*CtW*ee*MB**2*MT*vev*cmath.sqrt(2))/(Lambda**2*sw) + (9*CtW*ee*MT**3*vev*cmath.sqrt(2))/(Lambda**2*sw) - (9*CtW*ee*MT*MW**2*vev*cmath.sqrt(2))/(Lambda**2*sw) - (9*C3phiq*CtW*ee*MB**2*MT*vev**3*cmath.sqrt(2))/(Lambda**4*sw) + (9*C3phiq*CtW*ee*MT**3*vev**3*cmath.sqrt(2))/(Lambda**4*sw) - (9*C3phiq*CtW*ee*MT*MW**2*vev**3*cmath.sqrt(2))/(Lambda**4*sw) - (9*ee*MB**2*MT*vev*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**2*sw) + (9*ee*MT**3*vev*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**2*sw) - (9*ee*MT*MW**2*vev*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**2*sw) - (9*ee*MB**2*MT*vev**3*complexconjugate(C3phiq)*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**4*sw) + (9*ee*MT**3*vev**3*complexconjugate(C3phiq)*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**4*sw) - (9*ee*MT*MW**2*vev**3*complexconjugate(C3phiq)*complexconjugate(CtW)*cmath.sqrt(2))/(Lambda**4*sw))*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'(((-3*C3phiq**2*MB**2*MH**2*vev**2)/(2.*Lambda**4) - 12*MB**2*yb**2 + 3*MH**2*yb**2 + (3*C3phiq*MB**2*MH**2*vev**2*complexconjugate(C3phiq))/Lambda**4 - (3*MB**2*MH**2*vev**2*complexconjugate(C3phiq)**2)/(2.*Lambda**4))*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.g, P.g):
'(CphiG**2*MH**6*vev**2)/(2.*cmath.pi*Lambda**4*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'(((-3*C3phiq**2*MH**2*MT**2*vev**2)/(2.*Lambda**4) + 3*MH**2*yt**2 - 12*MT**2*yt**2 + (3*C3phiq*MH**2*MT**2*vev**2*complexconjugate(C3phiq))/Lambda**4 - (3*MH**2*MT**2*vev**2*complexconjugate(C3phiq)**2)/(2.*Lambda**4))*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
# This file was automatically created by FeynRules 2.3.31
# Mathematica version: 11.0.0 for Mac OS X x86 (64-bit) (July 28, 2016)
# Date: Tue 9 Apr 2019 14:23:53
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
particle=P.c,
partial_widths={
(P.W__plus__, P.b):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.d):
'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.s):
'(((3*CKM2x2*ee**2*MC**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MS**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MC**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MC**2*MS**2*complexconjugate(CKM2x2))/(MW**2*sw**2) + (3*CKM2x2*ee**2*MS**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MW**2*complexconjugate(CKM2x2))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MS**2 + MS**4 - 2*MC**2*MW**2 - 2*MS**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)'
})
# This file was automatically created by FeynRules 2.3.24
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
# Date: Thu 25 Oct 2018 14:19:25
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.S1, P.vm__tilde__):
'((MB**2 - MS1**2)*(3*gL32**2*MB**2 - 3*gL32**2*MS1**2))/(96.*cmath.pi*abs(MB)**3)',
(P.S1, P.vt__tilde__):
'((MB**2 - MS1**2)*(3*gL3**2*MB**2 - 3*gL3**2*MS1**2))/(96.*cmath.pi*abs(MB)**3)',
(P.U1, P.ta__minus__):
'((3*hL33**2*MB**2 + 3*hL33**2*MTA**2 + (3*hL33**2*MB**4)/MU1**2 - (6*hL33**2*MB**2*MTA**2)/MU1**2 + (3*hL33**2*MTA**4)/MU1**2 - 6*hL33**2*MU1**2)*cmath.sqrt(MB**4 - 2*MB**2*MTA**2 + MTA**4 - 2*MB**2*MU1**2 - 2*MTA**2*MU1**2 + MU1**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
# This file was automatically created by FeynRules 2.3.21
# Mathematica version: 10.3.1 for Mac OS X x86 (64-bit) (December 9, 2015)
# Date: Thu 20 Oct 2016 00:50:38
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.g, P.omega):
'((MB**2 - Momega**2)*((coef**2*I2b3*MB**4*complexconjugate(I2b3))/(8.*cmath.pi**4*moct**2) - (coef**2*I2b3*MB**2*Momega**2*complexconjugate(I2b3))/(4.*cmath.pi**4*moct**2) + (coef**2*I2b3*Momega**4*complexconjugate(I2b3))/(8.*cmath.pi**4*moct**2)))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.Zp, P.omega):
'(((3*gz**2*I2b3*MB**2*zomega**2*complexconjugate(I2b3))/4. + (3*gz**2*I2b3*Momega**2*zomega**2*complexconjugate(I2b3))/4. + (3*gz**2*I2b3*MB**4*zomega**2*complexconjugate(I2b3))/(4.*MZp**2) - (3*gz**2*I2b3*MB**2*Momega**2*zomega**2*complexconjugate(I2b3))/(2.*MZp**2) + (3*gz**2*I2b3*Momega**4*zomega**2*complexconjugate(I2b3))/(4.*MZp**2) - (3*gz**2*I2b3*MZp**2*zomega**2*complexconjugate(I2b3))/2.)*cmath.sqrt(MB**4 - 2*MB**2*Momega**2 + Momega**4 - 2*MB**2*MZp**2 - 2*Momega**2*MZp**2 + MZp**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
# This file was automatically created by FeynRules 2.3.2
# Mathematica version: 8.0 for Linux x86 (64-bit) (February 23, 2011)
# Date: Sun 7 Aug 2016 10:05:39
from object_library import all_decays, Decay
import particles as P
Decay_P__tilde__A0 = Decay(
name='Decay_P__tilde__A0',
particle=P.P__tilde__A0,
partial_widths={
(P.P__tilde__H0, P.Z):
'((-(ee**2*MA0**2) - ee**2*MH0**2 + (ee**2*MA0**4)/(2.*MZ**2) - (ee**2*MA0**2*MH0**2)/MZ**2 + (ee**2*MH0**4)/(2.*MZ**2) + (ee**2*MZ**2)/2. - (CW**2*ee**2*MA0**2)/(2.*SW**2) - (CW**2*ee**2*MH0**2)/(2.*SW**2) + (CW**2*ee**2*MA0**4)/(4.*MZ**2*SW**2) - (CW**2*ee**2*MA0**2*MH0**2)/(2.*MZ**2*SW**2) + (CW**2*ee**2*MH0**4)/(4.*MZ**2*SW**2) + (CW**2*ee**2*MZ**2)/(4.*SW**2) - (ee**2*MA0**2*SW**2)/(2.*CW**2) - (ee**2*MH0**2*SW**2)/(2.*CW**2) + (ee**2*MA0**4*SW**2)/(4.*CW**2*MZ**2) - (ee**2*MA0**2*MH0**2*SW**2)/(2.*CW**2*MZ**2) + (ee**2*MH0**4*SW**2)/(4.*CW**2*MZ**2) + (ee**2*MZ**2*SW**2)/(4.*CW**2))*cmath.sqrt(MA0**4 - 2*MA0**2*MH0**2 + MH0**4 - 2*MA0**2*MZ**2 - 2*MH0**2*MZ**2 + MZ**4))/(16.*cmath.pi*abs(MA0)**3)',
(P.P__tilde__H__plus__, P.W__minus__):
'((-(ee**2*MA0**2)/(2.*SW**2) - (ee**2*MHch**2)/(2.*SW**2) + (ee**2*MA0**4)/(4.*MW**2*SW**2) - (ee**2*MA0**2*MHch**2)/(2.*MW**2*SW**2) + (ee**2*MHch**4)/(4.*MW**2*SW**2) + (ee**2*MW**2)/(4.*SW**2))*cmath.sqrt(MA0**4 - 2*MA0**2*MHch**2 + MHch**4 - 2*MA0**2*MW**2 - 2*MHch**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MA0)**3)',
(P.P__tilde__H__minus__, P.W__plus__):
'((-(ee**2*MA0**2)/(2.*SW**2) - (ee**2*MHch**2)/(2.*SW**2) + (ee**2*MA0**4)/(4.*MW**2*SW**2) - (ee**2*MA0**2*MHch**2)/(2.*MW**2*SW**2) + (ee**2*MHch**4)/(4.*MW**2*SW**2) + (ee**2*MW**2)/(4.*SW**2))*cmath.sqrt(MA0**4 - 2*MA0**2*MHch**2 + MHch**4 - 2*MA0**2*MW**2 - 2*MHch**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MA0)**3)'
})
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*SW**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*SW**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*SW**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*SW**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*SW**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/SW**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*SW**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*SW**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*SW**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*SW**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*SW**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/SW**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*SW**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*SW**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*SW**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*SW**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*SW**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/SW**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.0.1 for Linux x86 (64-bit) (September 21, 2016)
# Date: Tue 27 Feb 2018 15:06:58
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.bsq,P.chi):'((3*gDM**2*MB**2 - 3*gDM**2*MBSQ**2 + 3*gDM**2*Mchi**2)*cmath.sqrt(MB**4 - 2*MB**2*MBSQ**2 + MBSQ**4 - 2*MB**2*Mchi**2 - 2*MBSQ**2*Mchi**2 + Mchi**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_bsq = Decay(name = 'Decay_bsq',
particle = P.bsq,
partial_widths = {(P.b,P.chi):'((-3*gDM**2*MB**2 + 3*gDM**2*MBSQ**2 - 3*gDM**2*Mchi**2)*cmath.sqrt(MB**4 - 2*MB**2*MBSQ**2 + MBSQ**4 - 2*MB**2*Mchi**2 - 2*MBSQ**2*Mchi**2 + Mchi**4))/(48.*cmath.pi*abs(MBSQ)**3)'})
Decay_chi = Decay(name = 'Decay_chi',
particle = P.chi,
partial_widths = {(P.bsq,P.b__tilde__):'((3*gDM**2*MB**2 - 3*gDM**2*MBSQ**2 + 3*gDM**2*Mchi**2)*cmath.sqrt(MB**4 - 2*MB**2*MBSQ**2 + MBSQ**4 - 2*MB**2*Mchi**2 - 2*MBSQ**2*Mchi**2 + Mchi**4))/(32.*cmath.pi*abs(Mchi)**3)',
(P.bsq__tilde__,P.b):'((3*gDM**2*MB**2 - 3*gDM**2*MBSQ**2 + 3*gDM**2*Mchi**2)*cmath.sqrt(MB**4 - 2*MB**2*MBSQ**2 + MBSQ**4 - 2*MB**2*Mchi**2 - 2*MBSQ**2*Mchi**2 + Mchi**4))/(32.*cmath.pi*abs(Mchi)**3)',
(P.csq,P.c__tilde__):'((Mchi**2 - MCSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MCSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.csq__tilde__,P.c):'((Mchi**2 - MCSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MCSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.dsq,P.d__tilde__):'((Mchi**2 - MDSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MDSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.dsq__tilde__,P.d):'((Mchi**2 - MDSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MDSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.ssq,P.s__tilde__):'((Mchi**2 - MSSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MSSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.ssq__tilde__,P.s):'((Mchi**2 - MSSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MSSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
(P.tsq,P.t__tilde__):'((3*gDM**2*Mchi**2 + 3*gDM**2*MT**2 - 3*gDM**2*MTSQ**2)*cmath.sqrt(Mchi**4 - 2*Mchi**2*MT**2 + MT**4 - 2*Mchi**2*MTSQ**2 - 2*MT**2*MTSQ**2 + MTSQ**4))/(32.*cmath.pi*abs(Mchi)**3)',
(P.tsq__tilde__,P.t):'((3*gDM**2*Mchi**2 + 3*gDM**2*MT**2 - 3*gDM**2*MTSQ**2)*cmath.sqrt(Mchi**4 - 2*Mchi**2*MT**2 + MT**4 - 2*Mchi**2*MTSQ**2 - 2*MT**2*MTSQ**2 + MTSQ**4))/(32.*cmath.pi*abs(Mchi)**3)',
(P.usq,P.u__tilde__):'((Mchi**2 - MUSQ**2)*(3*gDM**2*Mchi**2 - 3*gDM**2*MUSQ**2))/(32.*cmath.pi*abs(Mchi)**3)',
# This file was automatically created by FeynRules 2.4.68
# Mathematica version: 10.4.1 for Mac OS X x86 (64-bit) (April 11, 2016)
# Date: Thu 6 Jun 2019 21:52:44
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_t = Decay(
name='Decay_t',
particle=P.t,
partial_widths={
(P.W__plus__, P.b):
'((MT**2 - MW**2)*((3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2))/(96.*cmath.pi*abs(MT)**3)',
(P.Xs, P.YF3Qu3):
'((3*lamF3Q3x3**2*MT**2 - 3*lamF3Q3x3**2*MXs**2 + 3*lamF3Q3x3**2*MYF3Qu3**2)*cmath.sqrt(MT**4 - 2*MT**2*MXs**2 + MXs**4 - 2*MT**2*MYF3Qu3**2 - 2*MXs**2*MYF3Qu3**2 + MYF3Qu3**4))/(96.*cmath.pi*abs(MT)**3)',
(P.Xs, P.YF3u3):
'((3*lamF3u3x3**2*MT**2 - 3*lamF3u3x3**2*MXs**2 + 3*lamF3u3x3**2*MYF3u3**2)*cmath.sqrt(MT**4 - 2*MT**2*MXs**2 + MXs**4 - 2*MT**2*MYF3u3**2 - 2*MXs**2*MYF3u3**2 + MYF3u3**4))/(96.*cmath.pi*abs(MT)**3)',
# This file was automatically created by FeynRules 2.1.78
# Mathematica version: 9.0 for Mac OS X x86 (64-bit) (January 24, 2013)
# Date: Fri 17 Oct 2014 08:59:37
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(name = 'Decay_b',
particle = P.b,
partial_widths = {(P.sig8,P.d):'((-m8**2 + MB**2)*(-4*g8Ld1x3**2*m8**2 - 4*g8Ld3x1**2*m8**2 - 8*g8Ld3x1*g8Rd1x3*m8**2 - 4*g8Rd1x3**2*m8**2 - 8*g8Ld1x3*g8Rd3x1*m8**2 - 4*g8Rd3x1**2*m8**2 + 4*g8Ld1x3**2*MB**2 + 4*g8Ld3x1**2*MB**2 + 8*g8Ld3x1*g8Rd1x3*MB**2 + 4*g8Rd1x3**2*MB**2 + 8*g8Ld1x3*g8Rd3x1*MB**2 + 4*g8Rd3x1**2*MB**2))/(96.*cmath.pi*abs(MB)**3)',
(P.sig8,P.s):'((-m8**2 + MB**2)*(-4*g8Ld2x3**2*m8**2 - 4*g8Ld3x2**2*m8**2 - 8*g8Ld3x2*g8Rd2x3*m8**2 - 4*g8Rd2x3**2*m8**2 - 8*g8Ld2x3*g8Rd3x2*m8**2 - 4*g8Rd3x2**2*m8**2 + 4*g8Ld2x3**2*MB**2 + 4*g8Ld3x2**2*MB**2 + 8*g8Ld3x2*g8Rd2x3*MB**2 + 4*g8Rd2x3**2*MB**2 + 8*g8Ld2x3*g8Rd3x2*MB**2 + 4*g8Rd3x2**2*MB**2))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__,P.t):'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'})
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.b,P.b__tilde__):'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t,P.t__tilde__):'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__,P.ta__plus__):'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_sig8 = Decay(name = 'Decay_sig8',
particle = P.sig8,
partial_widths = {(P.b,P.b__tilde__):'((8*g8Ld3x3**2*m8**2 + 16*g8Ld3x3*g8Rd3x3*m8**2 + 8*g8Rd3x3**2*m8**2 - 32*g8Ld3x3**2*MB**2 - 64*g8Ld3x3*g8Rd3x3*MB**2 - 32*g8Rd3x3**2*MB**2)*cmath.sqrt(m8**4 - 4*m8**2*MB**2))/(128.*cmath.pi*abs(m8)**3)',
(P.b,P.d__tilde__):'((m8**2 - MB**2)*(4*g8Ld1x3**2*m8**2 + 4*g8Ld3x1**2*m8**2 + 8*g8Ld3x1*g8Rd1x3*m8**2 + 4*g8Rd1x3**2*m8**2 + 8*g8Ld1x3*g8Rd3x1*m8**2 + 4*g8Rd3x1**2*m8**2 - 4*g8Ld1x3**2*MB**2 - 4*g8Ld3x1**2*MB**2 - 8*g8Ld3x1*g8Rd1x3*MB**2 - 4*g8Rd1x3**2*MB**2 - 8*g8Ld1x3*g8Rd3x1*MB**2 - 4*g8Rd3x1**2*MB**2))/(128.*cmath.pi*abs(m8)**3)',
(P.b,P.s__tilde__):'((m8**2 - MB**2)*(4*g8Ld2x3**2*m8**2 + 4*g8Ld3x2**2*m8**2 + 8*g8Ld3x2*g8Rd2x3*m8**2 + 4*g8Rd2x3**2*m8**2 + 8*g8Ld2x3*g8Rd3x2*m8**2 + 4*g8Rd3x2**2*m8**2 - 4*g8Ld2x3**2*MB**2 - 4*g8Ld3x2**2*MB**2 - 8*g8Ld3x2*g8Rd2x3*MB**2 - 4*g8Rd2x3**2*MB**2 - 8*g8Ld2x3*g8Rd3x2*MB**2 - 4*g8Rd3x2**2*MB**2))/(128.*cmath.pi*abs(m8)**3)',
(P.c,P.c__tilde__):'(m8**2*(8*g8Lu2x2**2*m8**2 + 16*g8Lu2x2*g8Ru2x2*m8**2 + 8*g8Ru2x2**2*m8**2))/(128.*cmath.pi*abs(m8)**3)',
(P.c,P.t__tilde__):'((m8**2 - MT**2)*(4*g8Lu2x3**2*m8**2 + 4*g8Lu3x2**2*m8**2 + 8*g8Lu3x2*g8Ru2x3*m8**2 + 4*g8Ru2x3**2*m8**2 + 8*g8Lu2x3*g8Ru3x2*m8**2 + 4*g8Ru3x2**2*m8**2 - 4*g8Lu2x3**2*MT**2 - 4*g8Lu3x2**2*MT**2 - 8*g8Lu3x2*g8Ru2x3*MT**2 - 4*g8Ru2x3**2*MT**2 - 8*g8Lu2x3*g8Ru3x2*MT**2 - 4*g8Ru3x2**2*MT**2))/(128.*cmath.pi*abs(m8)**3)',
(P.c,P.u__tilde__):'(m8**2*(4*g8Lu1x2**2*m8**2 + 4*g8Lu2x1**2*m8**2 + 8*g8Lu2x1*g8Ru1x2*m8**2 + 4*g8Ru1x2**2*m8**2 + 8*g8Lu1x2*g8Ru2x1*m8**2 + 4*g8Ru2x1**2*m8**2))/(128.*cmath.pi*abs(m8)**3)',
import particles as P
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.d, P.d__tilde__):
'((-12*MD**2*ydo**2 + 3*MH**2*ydo**2)*cmath.sqrt(-4*MD**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.s, P.s__tilde__):
'((3*MH**2*ys**2 - 12*MS**2*ys**2)*cmath.sqrt(MH**4 - 4*MH**2*MS**2))/(16.*cmath.pi*abs(MH)**3)',
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.e__minus__, P.e__plus__):
'((-4*Me**2*ye**2 + MH**2*ye**2)*cmath.sqrt(-4*Me**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.mu__minus__, P.mu__plus__):
'((MH**2*ym**2 - 4*MMU**2*ym**2)*cmath.sqrt(MH**4 - 4*MH**2*MMU**2))/(16.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.u, P.u__tilde__):
'((3*MH**2*yup**2 - 12*MU**2*yup**2)*cmath.sqrt(MH**4 - 4*MH**2*MU**2))/(16.*cmath.pi*abs(MH)**3)',
(P.c, P.c__tilde__):
'((-12*MC**2*yc**2 + 3*MH**2*yc**2)*cmath.sqrt(-4*MC**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
(P.W__minus__, P.W__plus__):
'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z, P.Z):
'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'
})
Decay_so = Decay(name='Decay_so',
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 10.0 for Linux x86 (64-bit) (September 9, 2014)
# Date: Mon 2 Oct 2017 06:25:15
from object_library import all_decays, Decay
import particles as P
Decay_H = Decay(name = 'Decay_H',
particle = P.H,
partial_widths = {(P.W__minus__,P.W__plus__):'(((3*ee**4*vev**2)/(4.*sw**4) + (ee**4*MH**4*vev**2)/(16.*MW**4*sw**4) - (ee**4*MH**2*vev**2)/(4.*MW**2*sw**4))*cmath.sqrt(MH**4 - 4*MH**2*MW**2))/(16.*cmath.pi*abs(MH)**3)',
(P.Z,P.Z):'(((9*ee**4*vev**2)/2. + (3*ee**4*MH**4*vev**2)/(8.*MZ**4) - (3*ee**4*MH**2*vev**2)/(2.*MZ**2) + (3*cw**4*ee**4*vev**2)/(4.*sw**4) + (cw**4*ee**4*MH**4*vev**2)/(16.*MZ**4*sw**4) - (cw**4*ee**4*MH**2*vev**2)/(4.*MZ**2*sw**4) + (3*cw**2*ee**4*vev**2)/sw**2 + (cw**2*ee**4*MH**4*vev**2)/(4.*MZ**4*sw**2) - (cw**2*ee**4*MH**2*vev**2)/(MZ**2*sw**2) + (3*ee**4*sw**2*vev**2)/cw**2 + (ee**4*MH**4*sw**2*vev**2)/(4.*cw**2*MZ**4) - (ee**4*MH**2*sw**2*vev**2)/(cw**2*MZ**2) + (3*ee**4*sw**4*vev**2)/(4.*cw**4) + (ee**4*MH**4*sw**4*vev**2)/(16.*cw**4*MZ**4) - (ee**4*MH**2*sw**4*vev**2)/(4.*cw**4*MZ**2))*cmath.sqrt(MH**4 - 4*MH**2*MZ**2))/(32.*cmath.pi*abs(MH)**3)'})
Decay_W__plus__ = Decay(name = 'Decay_W__plus__',
particle = P.W__plus__,
partial_widths = {(P.c,P.b__tilde__):'(CKM2x3*ee**2*MW**4*complexconjugate(CKM2x3))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.c,P.d__tilde__):'(CKM2x1*ee**2*MW**4*complexconjugate(CKM2x1))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.c,P.s__tilde__):'(CKM2x2*ee**2*MW**4*complexconjugate(CKM2x2))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.t,P.b__tilde__):'(CKM3x3*ee**2*MW**4*complexconjugate(CKM3x3))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.t,P.d__tilde__):'(CKM3x1*ee**2*MW**4*complexconjugate(CKM3x1))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.t,P.s__tilde__):'(CKM3x2*ee**2*MW**4*complexconjugate(CKM3x2))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.u,P.b__tilde__):'(CKM1x3*ee**2*MW**4*complexconjugate(CKM1x3))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.u,P.d__tilde__):'(CKM1x1*ee**2*MW**4*complexconjugate(CKM1x1))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.u,P.s__tilde__):'(CKM1x2*ee**2*MW**4*complexconjugate(CKM1x2))/(16.*cmath.pi*sw**2*abs(MW)**3)',
(P.ve,P.e__plus__):'(ee**2*MW**4)/(48.*cmath.pi*sw**2*abs(MW)**3)',
(P.vm,P.mu__plus__):'(ee**2*MW**4)/(48.*cmath.pi*sw**2*abs(MW)**3)',
(P.vt,P.ta__plus__):'(ee**2*MW**4)/(48.*cmath.pi*sw**2*abs(MW)**3)'})
Decay_WR__plus__ = Decay(name = 'Decay_WR__plus__',
particle = P.WR__plus__,
# This file was automatically created by FeynRules 2.3.36
# Mathematica version: 12.0.0 for Mac OS X x86 (64-bit) (April 7, 2019)
# Date: Mon 27 Jan 2020 17:55:07
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.c):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.t):
'(((3*CKM3x3*ee**2*MB**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MT**2*complexconjugate(CKM3x3))/(2.*sw**2) + (3*CKM3x3*ee**2*MB**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MB**2*MT**2*complexconjugate(CKM3x3))/(MW**2*sw**2) + (3*CKM3x3*ee**2*MT**4*complexconjugate(CKM3x3))/(2.*MW**2*sw**2) - (3*CKM3x3*ee**2*MW**2*complexconjugate(CKM3x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)',
(P.W__minus__, P.u):
'(((3*CKM1x3*ee**2*MB**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MU**2*complexconjugate(CKM1x3))/(2.*sw**2) + (3*CKM1x3*ee**2*MB**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MB**2*MU**2*complexconjugate(CKM1x3))/(MW**2*sw**2) + (3*CKM1x3*ee**2*MU**4*complexconjugate(CKM1x3))/(2.*MW**2*sw**2) - (3*CKM1x3*ee**2*MW**2*complexconjugate(CKM1x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MU**2 + MU**4 - 2*MB**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_c = Decay(
name='Decay_c',
particle=P.c,
partial_widths={
(P.W__plus__, P.b):
'(((3*CKM2x3*ee**2*MB**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MC**2*complexconjugate(CKM2x3))/(2.*sw**2) + (3*CKM2x3*ee**2*MB**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MB**2*MC**2*complexconjugate(CKM2x3))/(MW**2*sw**2) + (3*CKM2x3*ee**2*MC**4*complexconjugate(CKM2x3))/(2.*MW**2*sw**2) - (3*CKM2x3*ee**2*MW**2*complexconjugate(CKM2x3))/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*MW**2 - 2*MC**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.d):
'(((3*CKM2x1*ee**2*MC**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MD**2*complexconjugate(CKM2x1))/(2.*sw**2) + (3*CKM2x1*ee**2*MC**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MC**2*MD**2*complexconjugate(CKM2x1))/(MW**2*sw**2) + (3*CKM2x1*ee**2*MD**4*complexconjugate(CKM2x1))/(2.*MW**2*sw**2) - (3*CKM2x1*ee**2*MW**2*complexconjugate(CKM2x1))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*MW**2 - 2*MD**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)',
(P.W__plus__, P.s):
'(((3*CKM2x2*ee**2*MC**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MS**2*complexconjugate(CKM2x2))/(2.*sw**2) + (3*CKM2x2*ee**2*MC**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MC**2*MS**2*complexconjugate(CKM2x2))/(MW**2*sw**2) + (3*CKM2x2*ee**2*MS**4*complexconjugate(CKM2x2))/(2.*MW**2*sw**2) - (3*CKM2x2*ee**2*MW**2*complexconjugate(CKM2x2))/sw**2)*cmath.sqrt(MC**4 - 2*MC**2*MS**2 + MS**4 - 2*MC**2*MW**2 - 2*MS**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MC)**3)'
})
# This file was automatically created by FeynRules =.2.2
# Mathematica version: 9.0 for Linux x86 (64-bit) (February 7, 2013)
# Date: Fri 29 May 2015 17:02:46
from object_library import all_decays, Decay
import particles as P
Decay_b = Decay(
name='Decay_b',
particle=P.b,
partial_widths={
(P.W__minus__, P.t):
'(((3*ee**2*MB**2)/(2.*sw**2) + (3*ee**2*MT**2)/(2.*sw**2) + (3*ee**2*MB**4)/(2.*MW**2*sw**2) - (3*ee**2*MB**2*MT**2)/(MW**2*sw**2) + (3*ee**2*MT**4)/(2.*MW**2*sw**2) - (3*ee**2*MW**2)/sw**2)*cmath.sqrt(MB**4 - 2*MB**2*MT**2 + MT**4 - 2*MB**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MB)**3)'
})
Decay_H = Decay(
name='Decay_H',
particle=P.H,
partial_widths={
(P.a, P.a):
'(kHaa**2*MH**6)/(64.*cmath.pi*abs(MH)**3)',
(P.a, P.Z):
'((MH**2 - MZ**2)*((kHZa**2*MH**4)/2. - kHZa**2*MH**2*MZ**2 + (kHZa**2*MZ**4)/2.))/(16.*cmath.pi*abs(MH)**3)',
(P.b, P.b__tilde__):
'((-12*MB**2*yb**2 + 3*MH**2*yb**2)*cmath.sqrt(-4*MB**2*MH**2 + MH**4))/(16.*cmath.pi*abs(MH)**3)',
(P.g, P.g):
'(kHgg**2*MH**6)/(8.*cmath.pi*abs(MH)**3)',
(P.ta__minus__, P.ta__plus__):
'((MH**2*ytau**2 - 4*MTA**2*ytau**2)*cmath.sqrt(MH**4 - 4*MH**2*MTA**2))/(16.*cmath.pi*abs(MH)**3)',
(P.t, P.t__tilde__):
'((3*MH**2*yt**2 - 12*MT**2*yt**2)*cmath.sqrt(MH**4 - 4*MH**2*MT**2))/(16.*cmath.pi*abs(MH)**3)',
