# This file was automatically created by FeynRules 2.3.24
# Mathematica version: 10.3.0 for Mac OS X x86 (64-bit) (October 9, 2015)
# Date: Mon 15 Oct 2018 09:52:26


from .object_library import all_decays, Decay
from . import particles as P


Decay_b = Decay(name = 'Decay_b',
                particle = P.b,
                partial_widths = {(P.DP__minus__,P.t):'(((6*MB**2*msbarDownQuark3**2)/VPi**2 - (6*MDPc**2*msbarDownQuark3**2)/VPi**2 + (6*MB**2*msbarUpQuark3**2)/VPi**2 - (6*MDPc**2*msbarUpQuark3**2)/VPi**2 - (24*MB*msbarDownQuark3*msbarUpQuark3*MT)/VPi**2 + (6*msbarDownQuark3**2*MT**2)/VPi**2 + (6*msbarUpQuark3**2*MT**2)/VPi**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.rho0,P.d):'(((3*CKM1x1*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MD**2 + MD**4 - 2*MB**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MB)**3)',
                                  (P.rho0,P.s):'(((3*CKM1x2*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MB)**3)',
                                  (P.rho__minus__,P.c):'(((3*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*ee**4*MC**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*ee**4*MB**2*MC**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*ee**4*MC**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*Mro**2 - 2*MC**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MB)**3)',
                                  (P.rho__minus__,P.t):'(((3*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) + (3*CKM3x3*ee**4*MT**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x3*ee**4*MB**2*MT**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x3*ee**4*MT**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MB)**3)',
                                  (P.rho__minus__,P.u):'(((3*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*sw**4) + (3*CKM1x3*ee**4*MU**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*ee**4*MB**2*MU**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*ee**4*MU**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MU**2 - 2*Mro**2*MU**2 + MU**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):'(((6*MC**2*msbarDownQuark2**2)/VPi**2 - (6*MDPc**2*msbarDownQuark2**2)/VPi**2 + (6*MS**2*msbarDownQuark2**2)/VPi**2 - (24*MC*MS*msbarDownQuark2*msbarUpQuark2)/VPi**2 + (6*MC**2*msbarUpQuark2**2)/VPi**2 - (6*MDPc**2*msbarUpQuark2**2)/VPi**2 + (6*MS**2*msbarUpQuark2**2)/VPi**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.rho__plus__,P.b):'(((3*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*ee**4*MC**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*ee**4*MB**2*MC**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*ee**4*MC**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*Mro**2 - 2*MC**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MC)**3)',
                                  (P.rho__plus__,P.d):'(((3*CKM2x1*ee**4*MC**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*ee**4*MC**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*ee**4*MC**2*MD**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MC)**3)',
                                  (P.rho__plus__,P.s):'(((3*CKM2x2*ee**4*MC**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*ee**4*MC**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*sw**4) + (3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*ee**4*MC**2*MS**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*Mro**2 + Mro**4 - 2*MC**2*MS**2 - 2*Mro**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)',
                                  (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.DP__minus__,P.u):'(((6*MD**2*msbarDownQuark1**2)/VPi**2 - (6*MDPc**2*msbarDownQuark1**2)/VPi**2 + (6*MD**2*msbarUpQuark1**2)/VPi**2 - (6*MDPc**2*msbarUpQuark1**2)/VPi**2 - (24*MD*msbarDownQuark1*msbarUpQuark1*MU)/VPi**2 + (6*msbarDownQuark1**2*MU**2)/VPi**2 + (6*msbarUpQuark1**2*MU**2)/VPi**2)*cmath.sqrt(MD**4 - 2*MD**2*MDPc**2 + MDPc**4 - 2*MD**2*MU**2 - 2*MDPc**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (P.rho0,P.b):'(((3*CKM1x1*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MD**2 + MD**4 - 2*MB**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (P.rho0,P.s):'(((3*CKM1x1*CKM1x2*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (P.rho__minus__,P.c):'(((3*CKM2x1*ee**4*MC**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*ee**4*MC**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*ee**4*MC**2*MD**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (P.rho__minus__,P.t):'(((3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) + (3*CKM3x1*ee**4*MT**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*ee**4*MD**2*MT**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*ee**4*MT**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (P.rho__minus__,P.u):'(((3*CKM1x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*sw**4) + (3*CKM1x1*ee**4*MU**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*ee**4*MD**2*MU**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*ee**4*MU**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MU**2 - 2*Mro**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MD)**3)',
                                  (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_DP__plus__ = Decay(name = 'Decay_DP__plus__',
                         particle = P.DP__plus__,
                         partial_widths = {(P.c,P.s__tilde__):'(((-6*MC**2*msbarDownQuark2**2)/VPi**2 + (6*MDPc**2*msbarDownQuark2**2)/VPi**2 - (6*MS**2*msbarDownQuark2**2)/VPi**2 + (24*MC*MS*msbarDownQuark2*msbarUpQuark2)/VPi**2 - (6*MC**2*msbarUpQuark2**2)/VPi**2 + (6*MDPc**2*msbarUpQuark2**2)/VPi**2 - (6*MS**2*msbarUpQuark2**2)/VPi**2)*cmath.sqrt(MC**4 - 2*MC**2*MDPc**2 + MDPc**4 - 2*MC**2*MS**2 - 2*MDPc**2*MS**2 + MS**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.DP0,P.rho__plus__):'(((-2*CL**2*MDP**2)/NDark - (2*CL**2*MDPc**2)/NDark + (CL**2*MDP**4)/(Mro**2*NDark) - (2*CL**2*MDP**2*MDPc**2)/(Mro**2*NDark) + (CL**2*MDPc**4)/(Mro**2*NDark) + (CL**2*Mro**2)/NDark)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*Mro**2 - 2*MDPc**2*Mro**2 + Mro**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.DP0,P.W__plus__):'(((-2*ee**2*MDP**2)/sw**2 - (2*ee**2*MDPc**2)/sw**2 + (ee**2*MDP**4)/(MW**2*sw**2) - (2*ee**2*MDP**2*MDPc**2)/(MW**2*sw**2) + (ee**2*MDPc**4)/(MW**2*sw**2) + (ee**2*MW**2)/sw**2)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*MW**2 - 2*MDPc**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.H,P.W__plus__):'(((MDPc**4*Xi**2)/VPi**2 - (2*MDPc**2*MH**2*Xi**2)/VPi**2 + (MH**4*Xi**2)/VPi**2 - (2*MDPc**2*MW**2*Xi**2)/VPi**2 - (2*MH**2*MW**2*Xi**2)/VPi**2 + (MW**4*Xi**2)/VPi**2)*cmath.sqrt(MDPc**4 - 2*MDPc**2*MH**2 + MH**4 - 2*MDPc**2*MW**2 - 2*MH**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.t,P.b__tilde__):'(((-6*MB**2*msbarDownQuark3**2)/VPi**2 + (6*MDPc**2*msbarDownQuark3**2)/VPi**2 - (6*MB**2*msbarUpQuark3**2)/VPi**2 + (6*MDPc**2*msbarUpQuark3**2)/VPi**2 + (24*MB*msbarDownQuark3*msbarUpQuark3*MT)/VPi**2 - (6*msbarDownQuark3**2*MT**2)/VPi**2 - (6*msbarUpQuark3**2*MT**2)/VPi**2)*cmath.sqrt(MB**4 - 2*MB**2*MDPc**2 + MDPc**4 - 2*MB**2*MT**2 - 2*MDPc**2*MT**2 + MT**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.u,P.d__tilde__):'(((-6*MD**2*msbarDownQuark1**2)/VPi**2 + (6*MDPc**2*msbarDownQuark1**2)/VPi**2 - (6*MD**2*msbarUpQuark1**2)/VPi**2 + (6*MDPc**2*msbarUpQuark1**2)/VPi**2 + (24*MD*msbarDownQuark1*msbarUpQuark1*MU)/VPi**2 - (6*msbarDownQuark1**2*MU**2)/VPi**2 - (6*msbarUpQuark1**2*MU**2)/VPi**2)*cmath.sqrt(MD**4 - 2*MD**2*MDPc**2 + MDPc**4 - 2*MD**2*MU**2 - 2*MDPc**2*MU**2 + MU**4))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.ve,P.e__plus__):'((MDPc**2 - Me**2)*((2*MDPc**2*Me**2)/VPi**2 - (2*Me**4)/VPi**2))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.vm,P.mu__plus__):'((MDPc**2 - MMU**2)*((2*MDPc**2*MMU**2)/VPi**2 - (2*MMU**4)/VPi**2))/(16.*cmath.pi*abs(MDPc)**3)',
                                           (P.vt,P.ta__plus__):'((MDPc**2 - MTA**2)*((2*MDPc**2*MTA**2)/VPi**2 - (2*MTA**4)/VPi**2))/(16.*cmath.pi*abs(MDPc)**3)'})

Decay_DP0 = Decay(name = 'Decay_DP0',
                  particle = P.DP0,
                  partial_widths = {(P.b,P.b__tilde__):'(3*MDP**2*msbarDownQuark3**2*cmath.sqrt(-4*MB**2*MDP**2 + MDP**4))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.c,P.c__tilde__):'(3*MDP**2*msbarUpQuark2**2*cmath.sqrt(-4*MC**2*MDP**2 + MDP**4))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.d,P.d__tilde__):'(3*MDP**2*msbarDownQuark1**2*cmath.sqrt(-4*MD**2*MDP**2 + MDP**4))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.DP__plus__,P.rho__minus__):'(((-2*CL**2*MDP**2)/NDark - (2*CL**2*MDPc**2)/NDark + (CL**2*MDP**4)/(Mro**2*NDark) - (2*CL**2*MDP**2*MDPc**2)/(Mro**2*NDark) + (CL**2*MDPc**4)/(Mro**2*NDark) + (CL**2*Mro**2)/NDark)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*Mro**2 - 2*MDPc**2*Mro**2 + Mro**4))/(16.*cmath.pi*abs(MDP)**3)',
                                    (P.DP__plus__,P.W__minus__):'(((-2*ee**2*MDP**2)/sw**2 - (2*ee**2*MDPc**2)/sw**2 + (ee**2*MDP**4)/(MW**2*sw**2) - (2*ee**2*MDP**2*MDPc**2)/(MW**2*sw**2) + (ee**2*MDPc**4)/(MW**2*sw**2) + (ee**2*MW**2)/sw**2)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*MW**2 - 2*MDPc**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MDP)**3)',
                                    (P.DP__minus__,P.rho__plus__):'(((-2*CL**2*MDP**2)/NDark - (2*CL**2*MDPc**2)/NDark + (CL**2*MDP**4)/(Mro**2*NDark) - (2*CL**2*MDP**2*MDPc**2)/(Mro**2*NDark) + (CL**2*MDPc**4)/(Mro**2*NDark) + (CL**2*Mro**2)/NDark)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*Mro**2 - 2*MDPc**2*Mro**2 + Mro**4))/(16.*cmath.pi*abs(MDP)**3)',
                                    (P.DP__minus__,P.W__plus__):'(((-2*ee**2*MDP**2)/sw**2 - (2*ee**2*MDPc**2)/sw**2 + (ee**2*MDP**4)/(MW**2*sw**2) - (2*ee**2*MDP**2*MDPc**2)/(MW**2*sw**2) + (ee**2*MDPc**4)/(MW**2*sw**2) + (ee**2*MW**2)/sw**2)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*MW**2 - 2*MDPc**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MDP)**3)',
                                    (P.e__minus__,P.e__plus__):'(MDP**2*Me**2*cmath.sqrt(MDP**4 - 4*MDP**2*Me**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.H,P.Z):'(((-2*MDP**2*MW**2*Xi**2)/(cw**2*VPi**2) - (2*MH**2*MW**2*Xi**2)/(cw**2*VPi**2) + (MDP**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) - (2*MDP**2*MH**2*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MH**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MW**2*MZ**2*Xi**2)/(cw**2*VPi**2))*cmath.sqrt(MDP**4 - 2*MDP**2*MH**2 + MH**4 - 2*MDP**2*MZ**2 - 2*MH**2*MZ**2 + MZ**4))/(16.*cmath.pi*abs(MDP)**3)',
                                    (P.mu__minus__,P.mu__plus__):'(MDP**2*MMU**2*cmath.sqrt(MDP**4 - 4*MDP**2*MMU**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.s,P.s__tilde__):'(3*MDP**2*msbarDownQuark2**2*cmath.sqrt(MDP**4 - 4*MDP**2*MS**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.t,P.t__tilde__):'(3*MDP**2*msbarUpQuark3**2*cmath.sqrt(MDP**4 - 4*MDP**2*MT**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.ta__minus__,P.ta__plus__):'(MDP**2*MTA**2*cmath.sqrt(MDP**4 - 4*MDP**2*MTA**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)',
                                    (P.u,P.u__tilde__):'(3*MDP**2*msbarUpQuark1**2*cmath.sqrt(MDP**4 - 4*MDP**2*MU**2))/(8.*cmath.pi*VPi**2*abs(MDP)**3)'})

Decay_e__minus__ = Decay(name = 'Decay_e__minus__',
                         particle = P.e__minus__,
                         partial_widths = {(P.DP__minus__,P.ve):'((-MDPc**2 + Me**2)*((-2*MDPc**2*Me**2)/VPi**2 + (2*Me**4)/VPi**2))/(32.*cmath.pi*abs(Me)**3)',
                                           (P.rho__minus__,P.ve):'((Me**2 - Mro**2)*((ee**4*Me**2*NDark)/(32.*cmath.pi**2*sw**4) + (ee**4*Me**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4) - (ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4)))/(32.*cmath.pi*abs(Me)**3)',
                                           (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)'})

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.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.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.DP__plus__,P.W__minus__):'(((MDPc**4*Xi**2)/VPi**2 - (2*MDPc**2*MH**2*Xi**2)/VPi**2 + (MH**4*Xi**2)/VPi**2 - (2*MDPc**2*MW**2*Xi**2)/VPi**2 - (2*MH**2*MW**2*Xi**2)/VPi**2 + (MW**4*Xi**2)/VPi**2)*cmath.sqrt(MDPc**4 - 2*MDPc**2*MH**2 + MH**4 - 2*MDPc**2*MW**2 - 2*MH**2*MW**2 + MW**4))/(16.*cmath.pi*abs(MH)**3)',
                                  (P.DP0,P.Z):'(((-2*MDP**2*MW**2*Xi**2)/(cw**2*VPi**2) - (2*MH**2*MW**2*Xi**2)/(cw**2*VPi**2) + (MDP**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) - (2*MDP**2*MH**2*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MH**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MW**2*MZ**2*Xi**2)/(cw**2*VPi**2))*cmath.sqrt(MDP**4 - 2*MDP**2*MH**2 + MH**4 - 2*MDP**2*MZ**2 - 2*MH**2*MZ**2 + MZ**4))/(16.*cmath.pi*abs(MH)**3)',
                                  (P.DP__minus__,P.W__plus__):'(((MDPc**4*Xi**2)/VPi**2 - (2*MDPc**2*MH**2*Xi**2)/VPi**2 + (MH**4*Xi**2)/VPi**2 - (2*MDPc**2*MW**2*Xi**2)/VPi**2 - (2*MH**2*MW**2*Xi**2)/VPi**2 + (MW**4*Xi**2)/VPi**2)*cmath.sqrt(MDPc**4 - 2*MDPc**2*MH**2 + MH**4 - 2*MDPc**2*MW**2 - 2*MH**2*MW**2 + MW**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.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.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.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.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_mu__minus__ = Decay(name = 'Decay_mu__minus__',
                          particle = P.mu__minus__,
                          partial_widths = {(P.DP__minus__,P.vm):'((-MDPc**2 + MMU**2)*((-2*MDPc**2*MMU**2)/VPi**2 + (2*MMU**4)/VPi**2))/(32.*cmath.pi*abs(MMU)**3)',
                                            (P.rho__minus__,P.vm):'((MMU**2 - Mro**2)*((ee**4*MMU**2*NDark)/(32.*cmath.pi**2*sw**4) + (ee**4*MMU**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4) - (ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4)))/(32.*cmath.pi*abs(MMU)**3)',
                                            (P.W__minus__,P.vm):'((MMU**2 - MW**2)*((ee**2*MMU**2)/(2.*sw**2) + (ee**2*MMU**4)/(2.*MW**2*sw**2) - (ee**2*MW**2)/sw**2))/(32.*cmath.pi*abs(MMU)**3)'})

Decay_rho__plus__ = Decay(name = 'Decay_rho__plus__',
                          particle = P.rho__plus__,
                          partial_widths = {(P.c,P.b__tilde__):'(((-3*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*ee**4*MC**2*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*ee**4*MB**2*MC**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*ee**4*MC**4*NDark*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MC**2 + MC**4 - 2*MB**2*Mro**2 - 2*MC**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.c,P.d__tilde__):'(((-3*CKM2x1*ee**4*MC**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*ee**4*MC**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*ee**4*MC**2*MD**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1))/(16.*cmath.pi**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*MD**2 + MD**4 - 2*MC**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.c,P.s__tilde__):'(((-3*CKM2x2*ee**4*MC**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*ee**4*MC**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*sw**4) - (3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*ee**4*MC**2*MS**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*Mro**2 + Mro**4 - 2*MC**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.DP0,P.DP__plus__):'(((-2*CL**2*MDP**2)/NDark - (2*CL**2*MDPc**2)/NDark + (CL**2*MDP**4)/(Mro**2*NDark) - (2*CL**2*MDP**2*MDPc**2)/(Mro**2*NDark) + (CL**2*MDPc**4)/(Mro**2*NDark) + (CL**2*Mro**2)/NDark)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*Mro**2 - 2*MDPc**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.t,P.b__tilde__):'(((-3*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) - (3*CKM3x3*ee**4*MT**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x3*ee**4*MB**2*MT**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x3*ee**4*MT**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.t,P.d__tilde__):'(((-3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) - (3*CKM3x1*ee**4*MT**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*ee**4*MD**2*MT**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*ee**4*MT**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.t,P.s__tilde__):'(((3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) - (3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*ee**4*MT**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM3x2*ee**4*MS**2*MT**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*ee**4*MT**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MT**2 - 2*MS**2*MT**2 + MT**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.u,P.b__tilde__):'(((-3*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*sw**4) - (3*CKM1x3*ee**4*MU**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*ee**4*MB**2*MU**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*ee**4*MU**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MU**2 - 2*Mro**2*MU**2 + MU**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.u,P.d__tilde__):'(((-3*CKM1x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*sw**4) - (3*CKM1x1*ee**4*MU**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*ee**4*MD**2*MU**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*ee**4*MU**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MU**2 - 2*Mro**2*MU**2 + MU**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.u,P.s__tilde__):'(((3*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*sw**4) - (3*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*ee**4*MU**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*ee**4*MS**2*MU**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*ee**4*MU**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MU**2 - 2*MS**2*MU**2 + MU**4))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.ve,P.e__plus__):'((-Me**2 + Mro**2)*(-(ee**4*Me**2*NDark)/(32.*cmath.pi**2*sw**4) - (ee**4*Me**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4) + (ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4)))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.vm,P.mu__plus__):'((-MMU**2 + Mro**2)*(-(ee**4*MMU**2*NDark)/(32.*cmath.pi**2*sw**4) - (ee**4*MMU**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4) + (ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4)))/(48.*cmath.pi*abs(Mro)**3)',
                                            (P.vt,P.ta__plus__):'((Mro**2 - MTA**2)*((ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4) - (ee**4*MTA**2*NDark)/(32.*cmath.pi**2*sw**4) - (ee**4*MTA**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4)))/(48.*cmath.pi*abs(Mro)**3)'})

Decay_rho0 = Decay(name = 'Decay_rho0',
                   particle = P.rho0,
                   partial_widths = {(P.b,P.b__tilde__):'(((-3*CKM1x3**2*ee**4*MB**2*NDark*complexconjugate(CKM1x3)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM1x3**2*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x3))/(16.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x3))/(16.*cmath.pi**2*sw**4) - (3*CKM2x3**2*ee**4*MB**2*NDark*complexconjugate(CKM2x3)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM2x3**2*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) - (3*CKM3x3**2*ee**4*MB**2*NDark*complexconjugate(CKM3x3)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM3x3**2*ee**4*Mro**2*NDark*complexconjugate(CKM3x3)**2)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(-4*MB**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.b,P.d__tilde__):'(((-3*CKM1x1*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MD**2 + MD**4 - 2*MB**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.b,P.s__tilde__):'(((-3*CKM1x2*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.c,P.c__tilde__):'(((-3*ee**4*MC**2*NDark)/(32.*cmath.pi**2*sw**4) + (3*ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(-4*MC**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.d,P.b__tilde__):'(((-3*CKM1x1*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x3*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x3*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MB**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*MB**2*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x3*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*MD**2 + MD**4 - 2*MB**2*Mro**2 - 2*MD**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.d,P.d__tilde__):'(((-3*CKM1x1**2*ee**4*MD**2*NDark*complexconjugate(CKM1x1)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM1x1**2*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x1))/(16.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x1))/(16.*cmath.pi**2*sw**4) - (3*CKM2x1**2*ee**4*MD**2*NDark*complexconjugate(CKM2x1)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM2x1**2*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) - (3*CKM3x1**2*ee**4*MD**2*NDark*complexconjugate(CKM3x1)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM3x1**2*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)**2)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(-4*MD**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.d,P.s__tilde__):'(((-3*CKM1x1*CKM1x2*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.DP__minus__,P.DP__plus__):'(((-4*CL**2*MDPc**2)/NDark + (CL**2*Mro**2)/NDark)*cmath.sqrt(-4*MDPc**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.e__minus__,P.e__plus__):'((-(ee**4*Me**2*NDark)/(32.*cmath.pi**2*sw**4) + (ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(-4*Me**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.mu__minus__,P.mu__plus__):'((-(ee**4*MMU**2*NDark)/(32.*cmath.pi**2*sw**4) + (ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(-4*MMU**2*Mro**2 + Mro**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.s,P.b__tilde__):'(((-3*CKM1x2*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.s,P.d__tilde__):'(((-3*CKM1x1*CKM1x2*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.s,P.s__tilde__):'(((3*CKM1x2**2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM1x2**2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x2))/(16.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x2))/(16.*cmath.pi**2*sw**4) + (3*CKM2x2**2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM2x2**2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)**2)/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) + (3*CKM3x2**2*ee**4*Mro**2*NDark*complexconjugate(CKM3x2)**2)/(32.*cmath.pi**2*sw**4) - (3*CKM3x2**2*ee**4*MS**2*NDark*complexconjugate(CKM3x2)**2)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(Mro**4 - 4*Mro**2*MS**2))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.t,P.t__tilde__):'(((3*ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4) - (3*ee**4*MT**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(Mro**4 - 4*Mro**2*MT**2))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.ta__minus__,P.ta__plus__):'(((ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4) - (ee**4*MTA**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(Mro**4 - 4*Mro**2*MTA**2))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.u,P.u__tilde__):'(((3*ee**4*Mro**2*NDark)/(32.*cmath.pi**2*sw**4) - (3*ee**4*MU**2*NDark)/(32.*cmath.pi**2*sw**4))*cmath.sqrt(Mro**4 - 4*Mro**2*MU**2))/(48.*cmath.pi*abs(Mro)**3)',
                                     (P.ve,P.ve__tilde__):'(ee**4*Mro**4*NDark)/(1536.*cmath.pi**3*sw**4*abs(Mro)**3)',
                                     (P.vm,P.vm__tilde__):'(ee**4*Mro**4*NDark)/(1536.*cmath.pi**3*sw**4*abs(Mro)**3)',
                                     (P.vt,P.vt__tilde__):'(ee**4*Mro**4*NDark)/(1536.*cmath.pi**3*sw**4*abs(Mro)**3)'})

Decay_s = Decay(name = 'Decay_s',
                particle = P.s,
                partial_widths = {(P.DP__minus__,P.c):'(((6*MC**2*msbarDownQuark2**2)/VPi**2 - (6*MDPc**2*msbarDownQuark2**2)/VPi**2 + (6*MS**2*msbarDownQuark2**2)/VPi**2 - (24*MC*MS*msbarDownQuark2*msbarUpQuark2)/VPi**2 + (6*MC**2*msbarUpQuark2**2)/VPi**2 - (6*MDPc**2*msbarUpQuark2**2)/VPi**2 + (6*MS**2*msbarUpQuark2**2)/VPi**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(MS)**3)',
                                  (P.rho0,P.b):'(((3*CKM1x2*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM1x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM1x3*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM1x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM2x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM2x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM2x3*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM2x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM1x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM2x3)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM1x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MB**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM2x3*CKM3x2*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*sw**4) - (3*CKM3x2*CKM3x3*ee**4*MB**2*MS**2*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*CKM3x3*ee**4*MS**4*NDark*complexconjugate(CKM3x2)*complexconjugate(CKM3x3))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.rho0,P.d):'(((3*CKM1x1*CKM1x2*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM1x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM1x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM2x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM2x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM2x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM2x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM2x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM1x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*sw**4) - (3*CKM2x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2)*complexconjugate(CKM3x1))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM1x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM1x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM2x2*CKM3x1*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*CKM3x1*ee**4*MS**4*NDark*complexconjugate(CKM2x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MD**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*sw**4) - (3*CKM3x1*CKM3x2*ee**4*MD**2*MS**2*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x1)*complexconjugate(CKM3x2))/(64.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.rho__minus__,P.c):'(((3*CKM2x2*ee**4*MC**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) + (3*CKM2x2*ee**4*MC**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM2x2*ee**4*Mro**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*sw**4) + (3*CKM2x2*ee**4*MS**2*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*sw**4) - (3*CKM2x2*ee**4*MC**2*MS**2*NDark*complexconjugate(CKM2x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM2x2*ee**4*MS**4*NDark*complexconjugate(CKM2x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MC**4 - 2*MC**2*Mro**2 + Mro**4 - 2*MC**2*MS**2 - 2*Mro**2*MS**2 + MS**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.rho__minus__,P.t):'(((-3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) + (3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*ee**4*MT**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM3x2*ee**4*MS**2*MT**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*ee**4*MT**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MT**2 - 2*MS**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.rho__minus__,P.u):'(((-3*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*sw**4) + (3*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*ee**4*MU**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*ee**4*MS**2*MU**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*ee**4*MU**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MU**2 - 2*MS**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.W__minus__,P.c):'(((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(MS)**3)',
                                  (P.W__minus__,P.t):'(((3*CKM3x2*ee**2*MS**2*complexconjugate(CKM3x2))/(2.*sw**2) + (3*CKM3x2*ee**2*MT**2*complexconjugate(CKM3x2))/(2.*sw**2) + (3*CKM3x2*ee**2*MS**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) - (3*CKM3x2*ee**2*MS**2*MT**2*complexconjugate(CKM3x2))/(MW**2*sw**2) + (3*CKM3x2*ee**2*MT**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) - (3*CKM3x2*ee**2*MW**2*complexconjugate(CKM3x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MT**2 + MT**4 - 2*MS**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MS)**3)',
                                  (P.W__minus__,P.u):'(((3*CKM1x2*ee**2*MS**2*complexconjugate(CKM1x2))/(2.*sw**2) + (3*CKM1x2*ee**2*MU**2*complexconjugate(CKM1x2))/(2.*sw**2) + (3*CKM1x2*ee**2*MS**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) - (3*CKM1x2*ee**2*MS**2*MU**2*complexconjugate(CKM1x2))/(MW**2*sw**2) + (3*CKM1x2*ee**2*MU**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) - (3*CKM1x2*ee**2*MW**2*complexconjugate(CKM1x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MU**2 + MU**4 - 2*MS**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MS)**3)'})

Decay_t = Decay(name = 'Decay_t',
                particle = P.t,
                partial_widths = {(P.DP__plus__,P.b):'(((6*MB**2*msbarDownQuark3**2)/VPi**2 - (6*MDPc**2*msbarDownQuark3**2)/VPi**2 + (6*MB**2*msbarUpQuark3**2)/VPi**2 - (6*MDPc**2*msbarUpQuark3**2)/VPi**2 - (24*MB*msbarDownQuark3*msbarUpQuark3*MT)/VPi**2 + (6*msbarDownQuark3**2*MT**2)/VPi**2 + (6*msbarUpQuark3**2*MT**2)/VPi**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(MT)**3)',
                                  (P.rho__plus__,P.b):'(((3*CKM3x3*ee**4*MB**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) + (3*CKM3x3*ee**4*MB**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x3*ee**4*Mro**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*sw**4) + (3*CKM3x3*ee**4*MT**2*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*sw**4) - (3*CKM3x3*ee**4*MB**2*MT**2*NDark*complexconjugate(CKM3x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x3*ee**4*MT**4*NDark*complexconjugate(CKM3x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MT)**3)',
                                  (P.rho__plus__,P.d):'(((3*CKM3x1*ee**4*MD**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) + (3*CKM3x1*ee**4*MD**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM3x1*ee**4*Mro**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*sw**4) + (3*CKM3x1*ee**4*MT**2*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*sw**4) - (3*CKM3x1*ee**4*MD**2*MT**2*NDark*complexconjugate(CKM3x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x1*ee**4*MT**4*NDark*complexconjugate(CKM3x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MT**2 - 2*Mro**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MT)**3)',
                                  (P.rho__plus__,P.s):'(((-3*CKM3x2*ee**4*Mro**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*sw**4) + (3*CKM3x2*ee**4*MS**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) + (3*CKM3x2*ee**4*MS**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*ee**4*MT**2*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*sw**4) - (3*CKM3x2*ee**4*MS**2*MT**2*NDark*complexconjugate(CKM3x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM3x2*ee**4*MT**4*NDark*complexconjugate(CKM3x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MT**2 - 2*MS**2*MT**2 + MT**4))/(96.*cmath.pi*abs(MT)**3)',
                                  (P.W__plus__,P.b):'(((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(MT)**3)',
                                  (P.W__plus__,P.d):'(((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(MT)**3)',
                                  (P.W__plus__,P.s):'(((3*CKM3x2*ee**2*MS**2*complexconjugate(CKM3x2))/(2.*sw**2) + (3*CKM3x2*ee**2*MT**2*complexconjugate(CKM3x2))/(2.*sw**2) + (3*CKM3x2*ee**2*MS**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) - (3*CKM3x2*ee**2*MS**2*MT**2*complexconjugate(CKM3x2))/(MW**2*sw**2) + (3*CKM3x2*ee**2*MT**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) - (3*CKM3x2*ee**2*MW**2*complexconjugate(CKM3x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MT**2 + MT**4 - 2*MS**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.DP__minus__,P.vt):'((-MDPc**2 + MTA**2)*((-2*MDPc**2*MTA**2)/VPi**2 + (2*MTA**4)/VPi**2))/(32.*cmath.pi*abs(MTA)**3)',
                                            (P.rho__minus__,P.vt):'((-Mro**2 + MTA**2)*(-(ee**4*Mro**2*NDark)/(16.*cmath.pi**2*sw**4) + (ee**4*MTA**2*NDark)/(32.*cmath.pi**2*sw**4) + (ee**4*MTA**4*NDark)/(32.*cmath.pi**2*Mro**2*sw**4)))/(32.*cmath.pi*abs(MTA)**3)',
                                            (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_u = Decay(name = 'Decay_u',
                particle = P.u,
                partial_widths = {(P.DP__plus__,P.d):'(((6*MD**2*msbarDownQuark1**2)/VPi**2 - (6*MDPc**2*msbarDownQuark1**2)/VPi**2 + (6*MD**2*msbarUpQuark1**2)/VPi**2 - (6*MDPc**2*msbarUpQuark1**2)/VPi**2 - (24*MD*msbarDownQuark1*msbarUpQuark1*MU)/VPi**2 + (6*msbarDownQuark1**2*MU**2)/VPi**2 + (6*msbarUpQuark1**2*MU**2)/VPi**2)*cmath.sqrt(MD**4 - 2*MD**2*MDPc**2 + MDPc**4 - 2*MD**2*MU**2 - 2*MDPc**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MU)**3)',
                                  (P.rho__plus__,P.b):'(((3*CKM1x3*ee**4*MB**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) + (3*CKM1x3*ee**4*MB**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x3*ee**4*Mro**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*sw**4) + (3*CKM1x3*ee**4*MU**2*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*sw**4) - (3*CKM1x3*ee**4*MB**2*MU**2*NDark*complexconjugate(CKM1x3))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x3*ee**4*MU**4*NDark*complexconjugate(CKM1x3))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MB**4 - 2*MB**2*Mro**2 + Mro**4 - 2*MB**2*MU**2 - 2*Mro**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MU)**3)',
                                  (P.rho__plus__,P.d):'(((3*CKM1x1*ee**4*MD**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) + (3*CKM1x1*ee**4*MD**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4) - (3*CKM1x1*ee**4*Mro**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*sw**4) + (3*CKM1x1*ee**4*MU**2*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*sw**4) - (3*CKM1x1*ee**4*MD**2*MU**2*NDark*complexconjugate(CKM1x1))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x1*ee**4*MU**4*NDark*complexconjugate(CKM1x1))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(MD**4 - 2*MD**2*Mro**2 + Mro**4 - 2*MD**2*MU**2 - 2*Mro**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MU)**3)',
                                  (P.rho__plus__,P.s):'(((-3*CKM1x2*ee**4*Mro**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*sw**4) + (3*CKM1x2*ee**4*MS**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) + (3*CKM1x2*ee**4*MS**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*ee**4*MU**2*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*sw**4) - (3*CKM1x2*ee**4*MS**2*MU**2*NDark*complexconjugate(CKM1x2))/(16.*cmath.pi**2*Mro**2*sw**4) + (3*CKM1x2*ee**4*MU**4*NDark*complexconjugate(CKM1x2))/(32.*cmath.pi**2*Mro**2*sw**4))*cmath.sqrt(Mro**4 - 2*Mro**2*MS**2 + MS**4 - 2*Mro**2*MU**2 - 2*MS**2*MU**2 + MU**4))/(96.*cmath.pi*abs(MU)**3)',
                                  (P.W__plus__,P.b):'(((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(MU)**3)',
                                  (P.W__plus__,P.d):'(((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(MU)**3)',
                                  (P.W__plus__,P.s):'(((3*CKM1x2*ee**2*MS**2*complexconjugate(CKM1x2))/(2.*sw**2) + (3*CKM1x2*ee**2*MU**2*complexconjugate(CKM1x2))/(2.*sw**2) + (3*CKM1x2*ee**2*MS**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) - (3*CKM1x2*ee**2*MS**2*MU**2*complexconjugate(CKM1x2))/(MW**2*sw**2) + (3*CKM1x2*ee**2*MU**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) - (3*CKM1x2*ee**2*MW**2*complexconjugate(CKM1x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MU**2 + MU**4 - 2*MS**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(96.*cmath.pi*abs(MU)**3)'})

Decay_W__plus__ = Decay(name = 'Decay_W__plus__',
                        particle = P.W__plus__,
                        partial_widths = {(P.c,P.b__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.c,P.d__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.c,P.s__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.DP__plus__,P.H):'(((MDPc**4*Xi**2)/VPi**2 - (2*MDPc**2*MH**2*Xi**2)/VPi**2 + (MH**4*Xi**2)/VPi**2 - (2*MDPc**2*MW**2*Xi**2)/VPi**2 - (2*MH**2*MW**2*Xi**2)/VPi**2 + (MW**4*Xi**2)/VPi**2)*cmath.sqrt(MDPc**4 - 2*MDPc**2*MH**2 + MH**4 - 2*MDPc**2*MW**2 - 2*MH**2*MW**2 + MW**4))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.DP0,P.DP__plus__):'(((-2*ee**2*MDP**2)/sw**2 - (2*ee**2*MDPc**2)/sw**2 + (ee**2*MDP**4)/(MW**2*sw**2) - (2*ee**2*MDP**2*MDPc**2)/(MW**2*sw**2) + (ee**2*MDPc**4)/(MW**2*sw**2) + (ee**2*MW**2)/sw**2)*cmath.sqrt(MDP**4 - 2*MDP**2*MDPc**2 + MDPc**4 - 2*MDP**2*MW**2 - 2*MDPc**2*MW**2 + MW**4))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.t,P.b__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.t,P.d__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.t,P.s__tilde__):'(((-3*CKM3x2*ee**2*MS**2*complexconjugate(CKM3x2))/(2.*sw**2) - (3*CKM3x2*ee**2*MT**2*complexconjugate(CKM3x2))/(2.*sw**2) - (3*CKM3x2*ee**2*MS**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) + (3*CKM3x2*ee**2*MS**2*MT**2*complexconjugate(CKM3x2))/(MW**2*sw**2) - (3*CKM3x2*ee**2*MT**4*complexconjugate(CKM3x2))/(2.*MW**2*sw**2) + (3*CKM3x2*ee**2*MW**2*complexconjugate(CKM3x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MT**2 + MT**4 - 2*MS**2*MW**2 - 2*MT**2*MW**2 + MW**4))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.u,P.b__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.u,P.d__tilde__):'(((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.u,P.s__tilde__):'(((-3*CKM1x2*ee**2*MS**2*complexconjugate(CKM1x2))/(2.*sw**2) - (3*CKM1x2*ee**2*MU**2*complexconjugate(CKM1x2))/(2.*sw**2) - (3*CKM1x2*ee**2*MS**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) + (3*CKM1x2*ee**2*MS**2*MU**2*complexconjugate(CKM1x2))/(MW**2*sw**2) - (3*CKM1x2*ee**2*MU**4*complexconjugate(CKM1x2))/(2.*MW**2*sw**2) + (3*CKM1x2*ee**2*MW**2*complexconjugate(CKM1x2))/sw**2)*cmath.sqrt(MS**4 - 2*MS**2*MU**2 + MU**4 - 2*MS**2*MW**2 - 2*MU**2*MW**2 + MW**4))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.ve,P.e__plus__):'((-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))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.vm,P.mu__plus__):'((-MMU**2 + MW**2)*(-(ee**2*MMU**2)/(2.*sw**2) - (ee**2*MMU**4)/(2.*MW**2*sw**2) + (ee**2*MW**2)/sw**2))/(48.*cmath.pi*abs(MW)**3)',
                                          (P.vt,P.ta__plus__):'((-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))/(48.*cmath.pi*abs(MW)**3)'})

Decay_Z = Decay(name = 'Decay_Z',
                particle = P.Z,
                partial_widths = {(P.b,P.b__tilde__):'((-7*ee**2*MB**2 + ee**2*MZ**2 - (3*cw**2*ee**2*MB**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) - (17*ee**2*MB**2*sw**2)/(6.*cw**2) + (5*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MB**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.c,P.c__tilde__):'((-11*ee**2*MC**2 - ee**2*MZ**2 - (3*cw**2*ee**2*MC**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*MC**2*sw**2)/(6.*cw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MC**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.d,P.d__tilde__):'((-7*ee**2*MD**2 + ee**2*MZ**2 - (3*cw**2*ee**2*MD**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) - (17*ee**2*MD**2*sw**2)/(6.*cw**2) + (5*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MD**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.DP0,P.H):'(((-2*MDP**2*MW**2*Xi**2)/(cw**2*VPi**2) - (2*MH**2*MW**2*Xi**2)/(cw**2*VPi**2) + (MDP**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) - (2*MDP**2*MH**2*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MH**4*MW**2*Xi**2)/(cw**2*MZ**2*VPi**2) + (MW**2*MZ**2*Xi**2)/(cw**2*VPi**2))*cmath.sqrt(MDP**4 - 2*MDP**2*MH**2 + MH**4 - 2*MDP**2*MZ**2 - 2*MH**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.DP__minus__,P.DP__plus__):'(((-4*cw**2*ee**2*MDPc**2)/sw**2 + (cw**2*ee**2*MZ**2)/sw**2)*cmath.sqrt(-4*MDPc**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.e__minus__,P.e__plus__):'((-5*ee**2*Me**2 - ee**2*MZ**2 - (cw**2*ee**2*Me**2)/(2.*sw**2) + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*Me**2*sw**2)/(2.*cw**2) + (5*ee**2*MZ**2*sw**2)/(2.*cw**2))*cmath.sqrt(-4*Me**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.mu__minus__,P.mu__plus__):'((-5*ee**2*MMU**2 - ee**2*MZ**2 - (cw**2*ee**2*MMU**2)/(2.*sw**2) + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*MMU**2*sw**2)/(2.*cw**2) + (5*ee**2*MZ**2*sw**2)/(2.*cw**2))*cmath.sqrt(-4*MMU**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.s,P.s__tilde__):'((-7*ee**2*MS**2 + ee**2*MZ**2 - (3*cw**2*ee**2*MS**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) - (17*ee**2*MS**2*sw**2)/(6.*cw**2) + (5*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MS**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.t,P.t__tilde__):'((-11*ee**2*MT**2 - ee**2*MZ**2 - (3*cw**2*ee**2*MT**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*MT**2*sw**2)/(6.*cw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MT**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.ta__minus__,P.ta__plus__):'((-5*ee**2*MTA**2 - ee**2*MZ**2 - (cw**2*ee**2*MTA**2)/(2.*sw**2) + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*MTA**2*sw**2)/(2.*cw**2) + (5*ee**2*MZ**2*sw**2)/(2.*cw**2))*cmath.sqrt(-4*MTA**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.u,P.u__tilde__):'((-11*ee**2*MU**2 - ee**2*MZ**2 - (3*cw**2*ee**2*MU**2)/(2.*sw**2) + (3*cw**2*ee**2*MZ**2)/(2.*sw**2) + (7*ee**2*MU**2*sw**2)/(6.*cw**2) + (17*ee**2*MZ**2*sw**2)/(6.*cw**2))*cmath.sqrt(-4*MU**2*MZ**2 + MZ**4))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.ve,P.ve__tilde__):'(MZ**2*(ee**2*MZ**2 + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (ee**2*MZ**2*sw**2)/(2.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.vm,P.vm__tilde__):'(MZ**2*(ee**2*MZ**2 + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (ee**2*MZ**2*sw**2)/(2.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
                                  (P.vt,P.vt__tilde__):'(MZ**2*(ee**2*MZ**2 + (cw**2*ee**2*MZ**2)/(2.*sw**2) + (ee**2*MZ**2*sw**2)/(2.*cw**2)))/(48.*cmath.pi*abs(MZ)**3)',
                                  (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)'})

