ChiChi2QQ_M2Radiation.F
No OneTemporary
Actions

Size

162 KB

Subscribers

None

ChiChi2QQ_M2Radiation.F
View Options

	************************************************

	double precision function radS()

	#include "ChiChi2QQ_Kinematics.h"
	#include "ChiChi2QQ_Model.h"
	#include "ChiChi2QQ_GenCouplings.h"

	integer i,j

	double precision temp1,temp2
	complex*16 DenSmZ, DenSmH(6)
	double precision DenSmZ1, DenSmH1(6)

	radS = 0d0

	***** widths of Z-boson and Higgs included in propagators

	if (icha.eq.1) then

	DenSmZ = (sqrtS*2 - MW2) + cIMW*wW
	DenSmZ1 = DenSmZdconjg(DenSmZ)/(sqrtS*2 - MW2)

	else

	DenSmZ = (sqrtS*2 - MZ2) + cIMZ*wZ
	DenSmZ1 = DenSmZdconjg(DenSmZ)/(sqrtS*2 - MZ2)

	endif

	***** itypeT,U follows the type of the sfermion exchanged in t,u-channel interactions
	***** att gives the type of the antiquark final state
	***** antitt is the opposite type to the final state quark type

	if (icha.eq.0) then
	att = tt
	itypeT = tt
	itypeU = tt
	antitt = tt
	else if (icha.eq.1) then
	* tt is defined as normally as up-type (quark) type in chargino-neutralino co-annihilations
	if (tt.eq.3) then
	att = tt+1
	itypeT = tt
	itypeU = tt+1
	antitt = tt+1
	else if (tt.eq.4) then
	att = tt-1
	itypeT = tt
	itypeU = tt-1
	antitt = tt-1
	endif
	else if (icha.eq.2) then
	if (tt.eq.3) then
	att = tt
	itypeT = tt-1
	itypeU = tt+1
	antitt = tt+1
	else if (tt.eq.4) then
	att = tt
	itypeT = tt-1
	itypeU = tt+1
	antitt = tt-1
	endif
	endif

	if (icha.eq.1) then
	c only charged Higgs exchange
	imin = 5
	imax = 6
	jmin = 5
	jmax = 6
	else
	c exchange of neutral Higgses
	imin = 1
	imax = 4
	jmin = 1
	jmax = 4
	endif

	do i = imin, imax
	DenSmH(i) = (sqrtS2 - MHk(i)2) + cIMHk(i)wHk(i)
	DenSmH1(i) = DenSmH(i)dconjg(DenSmH(i))/(sqrtS2 - MHk(i)*2)
	enddo


	** Zboson bremsstrahlung

	radS =

	** Zboson dipole contribution 312

	radS = radS - NCCF(64AL2GS*2pisqrtS2(6BLBR(1 - 6muu1*2)muu3**2 +
	- BL*2(1 - muu3*2 + 2muu1*2(-2 + 5muu32)) + BR2(1 - muu3*2 + 2muu1*2(-2 + 5muu32)))
	- ((-2 + 4muu32)/(-2 + x1 + x2) + (dsqrt((-1 + 4muu3*2)/(4muu32 - x22))*
	- (4muu34 + muu32(4 - 6x2) + (-1 + x2)(-1 + x1 + 2x2)))/((-1 + 2muu3*2)(-1 + x2))))/
	- (3d0(DenSmZdconjg(DenSmZ))*(1 - x2))

	** Zboson dipole contribution 321

	radS = radS - NCCF(64AL2GS*2pisqrtS2(6BLBR(1 - 6muu1*2)muu3**2 +
	- BL*2(1 - muu3*2 + 2muu1*2(-2 + 5muu32)) + BR2(1 - muu3*2 + 2muu1*2(-2 + 5muu32)))
	- ((-2 + 4muu32)/(-2 + x1 + x2) + (dsqrt((-1 + 4muu3*2)/(4muu32 - x12))*
	- (4muu34 + muu32(4 - 6x1) + (-1 + x1)(-1 + 2x1 + x2)))/((-1 + 2muu3*2)(-1 + x1))))/
	- (3d0(DenSmZdconjg(DenSmZ))*(1 - x1))


	do j = 1, 4
	do i = 1, 4

	** Higgs boson bremsstrahlung

	radS = radS + NCCF(-8GS2pisqrtS2(-(CL(i)CL(j))-CR(i)CR(j) + 2muu12(CL(i) + CR(i))*
	- (CL(j) + CR(j)))*
	- (DL(i)((4muu3*4(-2 + x1 + x2)*2 + (-1 + x1)(-1 + x2)(2 + x12 + 2x1(-1 + x2) + (-2 + x2)x2) +
	- 2muu32(-2 + x1*2(1 - 2x2) - 2x1(-2 + x2)x2 + x2*2))DL(j) +
	- 4muu32(muu3*2(-2 + x1 + x2)*2 - (-1 + x1)(-1 + x2)(-1 + x1 + x2))DR(j)) +
	- DR(i)(4muu3*2(muu3*2(-2 + x1 + x2)*2 - (-1 + x1)(-1 + x2)(-1 + x1 + x2))DL(j) +
	- (4muu34(-2 + x1 + x2)*2 + (-1 + x1)(-1 + x2)(2 + x12 + 2x1(-1 + x2) + (-2 + x2)x2) +
	- 2muu32(-2 + x1*2(1 - 2x2) - 2x1(-2 + x2)x2 + x2*2))DR(j))))/
	- ((-1 + x1)*2(-1 + x2)*2DenSmH(i)*dconjg(DenSmH(j)))

	** Hboson dipole contribution 312

	radS = radS - NCCF(-8GS2pisqrtS2((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x22))(4muu34 + muu32(4 - 6x2) +
	- (-1 + x2)(-1 + x1 + 2x2)))/((-1 + 2muu32)(-1 + x2)))*
	- (-(CL(i)CL(j)) - CR(i)CR(j) + 2muu12(CL(i) + CR(i))(CL(j) + CR(j)))(-(DL(i)*DL(j)) -
	- DR(i)DR(j) + 2muu3*2(DL(i) + DR(i))*(DL(j) + DR(j))))/
	- ((-1 + x2)DenSmH(i)dconjg(DenSmH(j)))

	** Hboson dipole contribution 321

	radS = radS - NCCF(-8GS2pisqrtS2((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x12))(4muu34 + muu32(4 - 6x1) +
	- (-1 + x1)(-1 + 2x1 + x2)))/((-1 + 2muu32)(-1 + x1)))*
	- (-(CL(i)CL(j)) - CR(i)CR(j) + 2muu12(CL(i) + CR(i))(CL(j) + CR(j)))(-(DL(i)*DL(j)) -
	- DR(i)DR(j) + 2muu3*2(DL(i) + DR(i))*(DL(j) + DR(j))))/
	- ((-1 + x1)DenSmH(i)dconjg(DenSmH(j)))

	enddo

	** Z-Higgs interference

	radS = radS + NCCF(32AL(BL - BR)GS2muu1muu3pisqrtS2(2muu32(-2 + x1 + x2)*2 - (-1 + x1)
	- (-1 + x2)(2 + x12 + 2x1(-1 + x2) - 2x2 + x2*2))(CL(j) - CR(j))*(DL(j) - DR(j)))/
	- (-dreal(DenSmZdconjg(DenSmH(j)))(-1 + x1)*2(-1 + x2)**2)

	** Z-Higgs dipole contribution 312

	radS = radS - 2d0NCCF(16AL(BL - BR)GS*2muu1muu3pisqrtS2((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x22))(4muu34 + muu32(4 - 6x2) + (-1 + x2)(-1 + x1 + 2x2)))/
	- ((-1 + 2muu32)(-1 + x2)))(CL(j) - CR(j))(DL(j) - DR(j)))/
	- (-dreal(DenSmZdconjg(DenSmH(j)))(-1 + x2))

	** Z-Higgs dipole contribution 321

	radS = radS - 2d0NCCF(16AL(BL - BR)GS*2muu1muu3pisqrtS2((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x12))(4muu34 + muu32(4 - 6x1) + (-1 + x1)(-1 + 2x1 + x2)))/
	- ((-1 + 2muu32)(-1 + x1)))(CL(j) - CR(j))(DL(j) - DR(j)))/
	- (-dreal(DenSmZdconjg(DenSmH(j)))(-1 + x1))

	enddo

	end

	************************************************

	double precision function radT()

	#include "ChiChi2QQ_Kinematics.h"
	#include "ChiChi2QQ_Model.h"
	#include "ChiChi2QQ_GenCouplings.h"

	integer i,j

	radT = 0d0

	do j = 1, 2
	do i = 1, 2

	** Tchannel bremsstrahlung 11

	radT = radT + NCCF(-4gs2(ELc(i,tt,gen)((2p1k1(-1 + x2) + 2p1k2(-1 + 2muu3**2 + x2) -
	- sqrtS*2(-1 + x2)*
	- (-1 + 2muu32 + x1 + x2))ELc(j,tt,gen) + 2muu1muu3sqrtS2(1 + 2muu32 - x2)ER(j,tt,gen)) +
	- ER(i,tt,gen)(2muu1muu3sqrtS*2(1 + 2muu32 - x2)ELc(j,tt,gen) + (2p1k1(-1 + x2) +
	- 2p1k2(-1 + 2muu32 + x2) - sqrtS2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))ER(j,tt,gen)))*
	- (FR(i,tt,gen)(muu1muu3sqrtS2FL(j,tt,gen) + p1k2FR(j,tt,gen)) + FL(i,tt,gen)(p1k2*FL(j,tt,gen) +
	- muu1muu3sqrtS*2FR(j,tt,gen))))/
	- (sqrtS*2(-1 + x2)*2(2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)*
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))

	** Tchannel bremsstrahlung 22

	radT = radT + NCCF(2gs2(ELc(i,tt,gen)((-2p1k1 + sqrtS*2x1)*ELc(j,tt,gen) +
	- 2muu1muu3sqrtS2ER(j,tt,gen)) + ER(i,tt,gen)(2muu1muu3sqrtS*2ELc(j,tt,gen) +
	- (-2p1k1 + sqrtS2x1) ER(j,tt,gen)))(FL(i,tt,gen)((2p1k2(-1 + x1) + 2p1k1*
	- (-1 + 2muu32 + x1) - sqrtS2(-1 + 2muu32 + x1))FL(j,tt,gen) + 2muu1muu3sqrtS2
	- (-1 - 2muu32 + x1)FR(j,tt,gen)) + FR(i,tt,gen)(-2muu1muu3sqrtS*2(1 + 2muu32 - x1)FL(j,tt,gen) +
	- (2p1k2(-1 + x1) + 2p1k1(-1 + 2muu32 + x1) - sqrtS2(-1 + 2muu32 + x1))FR(j,tt,gen))))/
	- (sqrtS*2(-1 + x1)*2(2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(j,tt,gen)**2))

	** Tchannel bremsstrahlung 33

	radT = radT + NCCF(4gs2(2p1k1 - 2p1k2 + sqrtS*2(2muu12 + 2muu3*2 - x1))
	- (ELc(i,tt,gen)((2p1k1 - sqrtS*2x1)ELc(j,tt,gen) - 2muu1muu3sqrtS*2ER(j,tt,gen)) -
	- ER(i,tt,gen)(2muu1muu3sqrtS*2ELc(j,tt,gen) + (-2p1k1 + sqrtS2x1)ER(j,tt,gen)))
	- (FR(i,tt,gen)(muu1muu3sqrtS2FL(j,tt,gen) + p1k2FR(j,tt,gen)) + FL(i,tt,gen)
	- (p1k2FL(j,tt,gen) + muu1muu3sqrtS2FR(j,tt,gen))))/
	- ((2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)*
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(j,tt,gen)*2)
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))

	** Tchannel bremsstrahlung 12

	radT = radT + NCCF(2gs2(ER(i,tt,gen)(2muu1muu3sqrtS*2ELc(j,tt,gen)(FL(i,tt,gen)
	- ((-2p1k2 + sqrtS2 + muu32(-4p1k1 + 4p1k2 + 2sqrtS2) - sqrtS2x1 + 2p1k1x2 - sqrtS*2x2)*
	- FL(j,tt,gen) + 2muu1muu3sqrtS2(4muu32 - x1 - x2)FR(j,tt,gen)) + FR(i,tt,gen)*
	- (2muu1muu3sqrtS2(4muu32 - x1 - x2)FL(j,tt,gen) + (-2p1k2 + sqrtS2 + muu32
	- (-4p1k1 + 4p1k2 + 2sqrtS2) - sqrtS2x1 + 2p1k1x2 - sqrtS*2x2)*FR(j,tt,gen))) +
	- ER(j,tt,gen)(2FR(i,tt,gen)(muu1muu3sqrtS2((2 - 4muu32)p1k1 + sqrtS*2 - 2p1k2*x1 -
	- 2sqrtS2x1 + 2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - sqrtS*2x2 + sqrtS*2x1x2)FL(j,tt,gen) +
	- (p1k1((4 - 8muu3*2)p1k2 + 2muu32sqrtS*2(-1 + x2)) + sqrtS*2(-2p1k2x1 + muu3*2
	- (2p1k2(1 + x1) - sqrtS*2(-1 + x2)) + muu1*2sqrtS*2(-1 + x1)(-1 + x2)))FR(j,tt,gen)) +
	- FL(i,tt,gen)((2p1k2sqrtS2 - 4p1k2*2x1 - 6p1k2sqrtS*2x1 + sqrtS*4x1 + 2p1k2sqrtS*2x1**2 -
	- sqrtS*4x1*2 + 2muu3*2(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) -
	- 2p1k2sqrtS*2(-1 + x2)) - 4p1k12x2 - 2p1k2sqrtS*2x2 + 4p1k2sqrtS*2x1x2 - sqrtS4x1*x2 +
	- 2p1k1(-2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + x2 + x1x2)))FL(j,tt,gen) +
	- 2muu1muu3sqrtS2((2 - 4muu32)p1k1 + sqrtS*2 - 2p1k2x1 - 2sqrtS*2x1 +
	- 2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - sqrtS*2x2 + sqrtS*2x1x2)FR(j,tt,gen)))) +
	c
	- ELc(i,tt,gen)(2muu1muu3sqrtS*2ER(j,tt,gen)(FL(i,tt,gen)
	- ((-2p1k2 + sqrtS2 + muu32(-4p1k1 + 4p1k2 + 2sqrtS2) - sqrtS2x1 + 2p1k1x2 - sqrtS*2x2)*
	- FL(j,tt,gen) + 2muu1muu3sqrtS2(4muu32 - x1 - x2)FR(j,tt,gen)) + FR(i,tt,gen)*
	- (2muu1muu3sqrtS2(4muu32 - x1 - x2)FL(j,tt,gen) + (-2p1k2 + sqrtS2 + muu32
	- (-4p1k1 + 4p1k2 + 2sqrtS2) - sqrtS2x1 + 2p1k1x2 - sqrtS*2x2)*FR(j,tt,gen))) +
	- ELc(j,tt,gen)(-2FL(i,tt,gen)((2p1k1((-2 + 4muu3*2)p1k2 - muu3*2sqrtS*2(-1 + x2)) +
	- sqrtS*2(2p1k2x1 - muu3*2(2p1k2(1 + x1) - sqrtS*2(-1 + x2)) + muu1*2sqrtS*2
	- (-1 + x1 + x2 - x1x2)))FL(j,tt,gen) - muu1muu3sqrtS*2((2 - 4muu32)p1k1 + sqrtS*2 - 2p1k2*x1 -
	- 2sqrtS2x1 + 2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - sqrtS*2x2 + sqrtS*2x1x2)
	- FR(j,tt,gen)) + FR(i,tt,gen)(2muu1muu3sqrtS*2((2 - 4muu32)p1k1 + sqrtS*2 - 2p1k2*x1 -
	- 2sqrtS2x1 + 2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - sqrtS*2x2 + sqrtS*2x1x2)
	- FL(j,tt,gen) + (2p1k2sqrtS*2 - 4p1k2*2x1 - 6p1k2sqrtS*2x1 + sqrtS*4x1 + 2p1k2sqrtS*2x1**2 -
	- sqrtS*4x1*2 + 2muu3*2(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) -
	- 2p1k2sqrtS*2(-1 + x2)) - 4p1k12x2 - 2p1k2sqrtS*2x2 + 4p1k2sqrtS*2x1x2 - sqrtS4x1*x2 +
	- 2p1k1(-2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + x2 + x1x2)))FR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(-1 + x2)(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (-2p1k1 - sqrtS2(muu12 + muu32 - x1) + MSf(j,tt,gen)**2))


	** Tchannel bremsstrahlung 13

	radT = radT + NCCF(4gs2(ELc(i,tt,gen)((2p1k1(4p1k2 - sqrtS*2(-1 + 2muu3*2 + x2)) +
	- sqrtS*2(2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - 2muu12sqrtS*2(-1 + x2) -
	- x1(4p1k2 + sqrtS2 - sqrtS2x2)))ELc(j,tt,gen) + 2muu1muu3sqrtS2(2p1k1 - 2p1k2 +
	- sqrtS*2(1 + 4muu32 - x1 - x2))ER(j,tt,gen)) + ER(i,tt,gen)(2muu1muu3sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))ELc(j,tt,gen) + (2p1k1(4*p1k2 -
	- sqrtS*2(-1 + 2muu3*2 + x2)) +
	- sqrtS*2(2muu32(2p1k2 + sqrtS2(1 + x1 - x2)) - 2muu12sqrtS*2(-1 + x2) - x1(4p1k2 +
	- sqrtS2 - sqrtS2x2)))ER(j,tt,gen)))*
	- (FR(i,tt,gen)(muu1muu3sqrtS2FL(j,tt,gen) + p1k2FR(j,tt,gen)) + FL(i,tt,gen)(p1k2*FL(j,tt,gen) +
	- muu1muu3sqrtS*2FR(j,tt,gen))))/
	- (sqrtS*2(-1 + x2)(2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(j,tt,gen)*2)
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))

	** Tchannel bremsstrahlung 23

	radT = radT + NCCF(4gs2(ELc(i,tt,gen)((-2p1k1 + sqrtS*2x1)*ELc(j,tt,gen) +
	- 2muu1muu3sqrtS2ER(j,tt,gen)) + ER(i,tt,gen)(2muu1muu3sqrtS*2ELc(j,tt,gen) +
	- (-2p1k1 + sqrtS2x1)ER(j,tt,gen)))(FL(i,tt,gen)((p1k1(-4p1k2 + 2muu3*2sqrtS2) + sqrtS2*
	- (sqrtS*2(-muu32 + muu12(-1 + x1)) + p1k2(1 - 2muu32 + x1)))FL(j,tt,gen) -
	- muu1muu3sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - 2x1))*FR(j,tt,gen)) +
	- FR(i,tt,gen)(-(muu1muu3sqrtS2(2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - 2x1))*FL(j,tt,gen)) +
	- (p1k1(-4p1k2 + 2muu32sqrtS2) + sqrtS2(sqrtS2(-muu32 + muu12*(-1 + x1)) +
	- p1k2(1 - 2muu3*2 + x1)))FR(j,tt,gen))))/
	- (sqrtS*2(-1 + x1)(2p1k1 + sqrtS*2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(j,tt,gen)*2)
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))

	** Tchannel dipole 312

	radT = radT - NCCF(-32GS2(-4muu32 + x22)((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x22))(4muu34 + muu32(4 - 6x2) + (-1 + x2)(-1 + x1 + 2x2)))/
	- ((-1 + 2muu32)(-1 + x2)))(muu1muu3sqrtS2(ELc(j,tt,gen)ER(i,tt,gen) + ELc(i,tt,gen)ER(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) + dsqrt(-4muu32 + x22)))
	- (ELc(i,tt,gen)ELc(j,tt,gen) + ER(i,tt,gen)ER(j,tt,gen)))/(4.dsqrt(-4muu32 + x22)))*
	- (muu1muu3sqrtS*2(FL(j,tt,gen)FR(i,tt,gen) + FL(i,tt,gen)FR(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) + dsqrt(-4muu32 + x22)))
	- (FL(i,tt,gen)FL(j,tt,gen) + FR(i,tt,gen)FR(j,tt,gen)))/(4.dsqrt(-4muu32 + x22))))/
	- (sqrtS*2(-1 + x2)(-4dsqrt(1 - 4muu32)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 +
	- (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x2*2)) -
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)2)(-4dsqrt(1 - 4muu3*2)p1k2 +
	- sqrtS*2(dsqrt(1 - 4muu32)x2 + (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x2*2)) -
	- 2dsqrt(-4muu32 + x22)MSf(j,tt,gen)*2))

	** Tchannel dipole 321

	radT = radT - NCCF(-32GS2(-4muu32 + x12)((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x12))(4muu34 + muu32(4 - 6x1) + (-1 + x1)(-1 + 2x1 + x2)))/
	- ((-1 + 2muu32)(-1 + x1)))(muu1muu3sqrtS2(ELc(j,tt,gen)ER(i,tt,gen) + ELc(i,tt,gen)ER(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12)))
	- (ELc(i,tt,gen)ELc(j,tt,gen) + ER(i,tt,gen)ER(j,tt,gen)))/(4.dsqrt(-4muu32 + x12)))*
	- (muu1muu3sqrtS*2(FL(j,tt,gen)FR(i,tt,gen) + FL(i,tt,gen)FR(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12)))
	- (FL(i,tt,gen)FL(j,tt,gen) + FR(i,tt,gen)FR(j,tt,gen)))/(4.dsqrt(-4muu32 + x12))))/
	- (sqrtS*2(-1 + x1)(-4dsqrt(1 - 4muu32)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 +
	- (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x1*2)) +
	- 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)2)(-4dsqrt(1 - 4muu3*2)p1k1 +
	- sqrtS*2(dsqrt(1 - 4muu32)x1 + (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x1*2)) +
	- 2dsqrt(-4muu32 + x12)MSf(j,tt,gen)*2))


	enddo
	enddo

	end

	************************************************

	double precision function radU()

	#include "ChiChi2QQ_Kinematics.h"
	#include "ChiChi2QQ_Model.h"
	#include "ChiChi2QQ_GenCouplings.h"

	integer i,j

	radU = 0d0

	do j = 1, 2
	do i = 1, 2

	** Uchannel bremsstrahlung 11

	radU = radU + NCCF(-2gs2(GL(i,tt,gen)((2p1k2 - sqrtS*2x2)*GL(j,tt,gen) -
	- 2muu1muu3sqrtS2GR(j,tt,gen)) - GR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen) +
	- (-2p1k2 + sqrtS2x2)GR(j,tt,gen)))
	- (HL(i,tt,gen)((2p1k1(-1 + x2) + (2p1k2 - sqrtS*2)(-1 + 2muu32 + x2))HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(-1 - 2muu32 + x2)HR(j,tt,gen)) +
	- HR(i,tt,gen)(-2muu1muu3sqrtS*2(1 + 2muu32 - x2)HL(j,tt,gen) + (2p1k1(-1 + x2) +
	- (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))HR(j,tt,gen))))/
	- (sqrtS*2(-1 + x2)*2(2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)*2)
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)**2))


	** Uchannel bremsstrahlung 22

	radU = radU + NCCF(-4gs2(GL(i,tt,gen)((2p1k1(-1 + 2muu3*2 + x1) + (-1 + x1)
	- (2p1k2 - sqrtS2(-1 + 2muu32 + x1 + x2)))GL(j,tt,gen) +
	- 2muu1muu3sqrtS2(1 + 2muu32 - x1)GR(j,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3sqrtS*2(1 + 2muu32 - x1)GL(j,tt,gen) +
	- (2p1k1(-1 + 2muu32 + x1) + (-1 + x1)(2p1k2 - sqrtS2(-1 + 2muu32 + x1 + x2)))GR(j,tt,gen)))*
	- (HR(i,tt,gen)(muu1muu3sqrtS2HL(j,tt,gen) + p1k1HR(j,tt,gen)) + HL(i,tt,gen)
	- (p1k1HL(j,tt,gen) + muu1muu3sqrtS2HR(j,tt,gen))))/
	- (sqrtS*2(-1 + x1)*2(2p1k1 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)*
	- (2p1k1 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))


	** Uchannel bremsstrahlung 33

	radU = radU + (-4CFgs*2NC(2p1k1 - 2p1k2 + sqrtS2(-2muu12 - 2muu3*2 + x2))
	- (GL(i,tt,gen)((2p1k2 - sqrtS*2x2)GL(j,tt,gen) - 2muu1muu3sqrtS*2GR(j,tt,gen)) -
	- GR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen) + (-2p1k2 + sqrtS2x2)GR(j,tt,gen)))
	- (HR(i,tt,gen)(muu1muu3sqrtS2HL(j,tt,gen) + p1k1HR(j,tt,gen)) + HL(i,tt,gen)
	- (p1k1HL(j,tt,gen) + muu1muu3sqrtS2HR(j,tt,gen))))/
	- ((2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)*2)
	- (2p1k1 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)*2)
	- (2p1k1 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))


	** Uchannel bremsstrahlung 12

	radU = radU + NCCF(2gs2(GL(i,tt,gen)(2muu1muu3sqrtS*2GR(j,tt,gen)*
	- (HL(i,tt,gen)(((-2 + 4muu3*2)p1k1 + sqrtS2 + muu32(-4p1k2 + 2sqrtS*2) +
	- 2p1k2x1 - sqrtS*2x1 - sqrtS*2x2)*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(4muu32 - x1 - x2)HR(j,tt,gen)) +
	- HR(i,tt,gen)(2muu1muu3sqrtS*2(4muu32 - x1 - x2)HL(j,tt,gen) +
	- ((-2 + 4muu32)p1k1 + sqrtS2 + muu32(-4p1k2 + 2sqrtS2) + 2p1k2*x1 -
	- sqrtS*2x1 - sqrtS*2x2)*HR(j,tt,gen))) +
	- GL(j,tt,gen)(-2HL(i,tt,gen)((sqrtS2(-1 + x1)(muu32(-2p1k2 + sqrtS*2) -
	- muu1*2sqrtS*2(-1 + x2)) +
	- 2p1k1((-2 + 4muu32)p1k2 - sqrtS*2(-x2 + muu3*2(1 + x2))))*HL(j,tt,gen) +
	- muu1muu3sqrtS*2(-2p1k2 - sqrtS2 + sqrtS2x1 + 2p1k1x2 + 2sqrtS2x2 - sqrtS*2x1*x2 -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HR(j,tt,gen)) +
	- HR(i,tt,gen)(2muu1muu3sqrtS*2(2p1k2 + sqrtS2 - sqrtS2x1 - 2p1k1x2 - 2sqrtS2x2 +
	- sqrtS*2x1x2 + 2muu3*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HL(j,tt,gen) +
	- (-4p1k12x2 - (2p1k2 - sqrtS2x2)(2p1k2x1 - sqrtS2(-1 + x1 + x2)) +
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) +
	- 2p1k1(-2p1k2(-2 + x1 + x2) + sqrtS*2(1 - x1 - 3x2 + 2x1x2 + x22)))HR(j,tt,gen)))) +
	- GR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen)(HL(i,tt,gen)
	- (((-2 + 4muu32)p1k1 + sqrtS2 + muu32(-4p1k2 + 2sqrtS2) + 2p1k2*x1 -
	- sqrtS*2x1 - sqrtS*2x2)*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(4muu32 - x1 - x2)HR(j,tt,gen)) +
	- HR(i,tt,gen)(2muu1muu3sqrtS*2(4muu32 - x1 - x2)HL(j,tt,gen) +
	- ((-2 + 4muu32)p1k1 + sqrtS2 + muu32(-4p1k2 + 2sqrtS2) + 2p1k2x1 - sqrtS2x1 -
	- sqrtS*2x2)*HR(j,tt,gen))) +
	- GR(j,tt,gen)(HL(i,tt,gen)((-4p1k12x2 - (2p1k2 - sqrtS2x2)(2p1k2x1 - sqrtS2(-1 + x1 + x2)) +
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) +
	- 2p1k1(-2p1k2(-2 + x1 + x2) + sqrtS*2(1 - x1 - 3x2 + 2x1x2 + x22)))HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(2p1k2 + sqrtS2 - sqrtS2x1 - 2p1k1x2 - 2sqrtS2x2 + sqrtS*2x1*x2 +
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HR(j,tt,gen)) +
	- 2HR(i,tt,gen)(muu1muu3sqrtS*2(2p1k2 + sqrtS2 - sqrtS2x1 - 2p1k1x2 - 2sqrtS2x2 +
	- sqrtS*2x1x2 + 2muu3*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HL(j,tt,gen) +
	- (sqrtS*2(-1 + x1)(muu32(2p1k2 - sqrtS2) + muu12sqrtS*2(-1 + x2)) +
	- p1k1((4 - 8muu3*2)p1k2 + 2sqrtS2(-x2 + muu3*2(1 + x2))))*HR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(-1 + x2)(-2p1k2 - sqrtS2(muu12 + muu32 - x2) + MSf(i,tt,gen)*2)
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2))



	** Uchannel bremsstrahlung 13

	radU = radU + NCCF(4gs2(GL(i,tt,gen)((-2p1k2 + sqrtS*2x2)*GL(j,tt,gen) +
	- 2muu1muu3sqrtS2GR(j,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen) + (-2p1k2 + sqrtS2x2)GR(j,tt,gen)))
	- (HL(i,tt,gen)((muu32(2p1k2sqrtS2 - sqrtS4) + muu1*2sqrtS*4(-1 + x2) + p1k1*
	- (-4p1k2 + sqrtS2(1 - 2muu32 + x2)))HL(j,tt,gen) -
	- muu1muu3sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - 2x2))*HR(j,tt,gen)) +
	- HR(i,tt,gen)(-(muu1muu3sqrtS2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - 2x2))*HL(j,tt,gen)) +
	- (muu3*2(2p1k2sqrtS2 - sqrtS4) + muu1*2sqrtS*4(-1 + x2) + p1k1(-4p1k2 + sqrtS*2
	- (1 - 2muu32 + x2)))HR(j,tt,gen))))/
	- (sqrtS*2(-1 + x2)(2p1k2 + sqrtS*2(muu12 + muu32 - x2) - MSf(i,tt,gen)*2)
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)*2)
	- (2p1k1 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))



	** Uchannel bremsstrahlung 23

	radU = radU + NCCF(4gs2(GL(i,tt,gen)((4p1k1(2p1k2 + sqrtS*2(muu3**2 - x2)) +
	- sqrtS*2(-2p1k2(-1 + 2muu32 + x1) + sqrtS2(-2muu12(-1 + x1) + (-1 + x1)x2 + muu32
	- (2 - 2x1 + 2x2))))*GL(j,tt,gen) +
	- 2muu1muu3sqrtS2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))GR(j,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))GL(j,tt,gen) +
	- (4p1k1(2p1k2 + sqrtS2(muu32 - x2)) + sqrtS2*
	- (-2p1k2(-1 + 2muu32 + x1) + sqrtS2(-2muu12(-1 + x1) + (-1 + x1)x2 + muu32
	- (2 - 2x1 + 2x2))))GR(j,tt,gen)))
	- (HR(i,tt,gen)(muu1muu3sqrtS2HL(j,tt,gen) + p1k1HR(j,tt,gen)) + HL(i,tt,gen)
	- (p1k1HL(j,tt,gen) + muu1muu3sqrtS2HR(j,tt,gen))))/
	- (sqrtS*2(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)(-1 + x1)
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)**2))


	** Uchannel dipole 312

	radU = radU - NCCF(-32GS2(-4muu32 + x22)((-2 + 4muu3*2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x22))(4muu34 + muu32(4 - 6x2) + (-1 + x2)(-1 + x1 + 2x2)))/
	- ((-1 + 2muu32)(-1 + x2)))*
	- (muu1muu3sqrtS*2(GL(j,tt,gen)GR(i,tt,gen) + GL(i,tt,gen)GR(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22)))
	- (GL(i,tt,gen)GL(j,tt,gen) + GR(i,tt,gen)GR(j,tt,gen)))/
	- (4d0dsqrt(-4muu32 + x22)))(muu1muu3sqrtS2(HL(j,tt,gen)HR(i,tt,gen) + HL(i,tt,gen)HR(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22)))
	- (HL(i,tt,gen)HL(j,tt,gen) + HR(i,tt,gen)HR(j,tt,gen)))/
	- (4d0dsqrt(-4muu32 + x22))))/
	- (sqrtS*2(-1 + x2)(-4dsqrt(1 - 4muu32)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 +
	- (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x2*2)) +
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)2)(-4dsqrt(1 - 4muu3*2)p1k2 +
	- sqrtS*2(dsqrt(1 - 4muu32)x2 + (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x2*2)) +
	- 2dsqrt(-4muu32 + x22)MSf(j,tt,gen)*2))

	** Uchannel dipole 321

	radU = radU - NCCF(8GS2((-2 + 4muu32)/(-2 + x1 + x2) + (dsqrt((-1 + 4muu3**2)/
	- (4muu32 - x12))(4muu34 + muu32(4 - 6x1) + (-1 + x1)(-1 + 2*x1 + x2)))/
	- ((-1 + 2muu32)(-1 + x1)))*
	- (muu1muu3sqrtS*2(GL(j,tt,gen)GR(i,tt,gen) + GL(i,tt,gen)GR(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + dsqrt(-4muu32 + x12)))
	- (GL(i,tt,gen)GL(j,tt,gen) + GR(i,tt,gen)GR(j,tt,gen)))/(4.dsqrt(-4muu32 + x12)))*
	- (muu1muu3sqrtS*2(HL(j,tt,gen)HR(i,tt,gen) + HL(i,tt,gen)HR(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + dsqrt(-4muu32 + x12)))
	- (HL(i,tt,gen)HL(j,tt,gen) + HR(i,tt,gen)HR(j,tt,gen)))/(4.dsqrt(-4muu32 + x12))))/
	- (sqrtS*2(1 - x1)((2dsqrt(1 - 4muu32)p1k1)/dsqrt(-4muu32 + x1*2) -
	- (sqrtS*2(-1 + 2muu12 + 2muu3**2 +
	- (dsqrt(1 - 4muu32)x1)/dsqrt(-4muu32 + x12)))/2d0 + MSf(i,tt,gen)2)
	- ((2dsqrt(1 - 4muu3*2)p1k1)/dsqrt(-4muu32 + x12) - (sqrtS2(-1 + 2muu12 + 2muu3**2 +
	- (dsqrt(1 - 4muu32)x1)/dsqrt(-4muu32 + x12)))/2d0 + MSf(j,tt,gen)*2))

	enddo
	enddo

	end


	************************************************

	double precision function radTS()

	#include "ChiChi2QQ_Kinematics.h"
	#include "ChiChi2QQ_Model.h"
	#include "ChiChi2QQ_GenCouplings.h"

	integer i,j
	complex*16 DenSmZ, DenSmH(4)
	double precision DenSmZ1, DenSmH1(4)

	DenSmZ = (sqrtS*2 - MZ2) + cIMZ*wZ
	DenSmZ1 = DenSmZdconjg(DenSmZ)/(sqrtS*2 - MZ2)

	do i = 1, 4
	DenSmH(i) = (sqrtS2 - MHk0(i)2) + cIMHk0(i)wHk0(i)
	DenSmH1(i) = DenSmH(i)dconjg(DenSmH(i))/(sqrtS2 - MHk0(i)*2)
	enddo

	radTS = 0d0

	do i = 1, 2

	** TSchannel bremsstrahlung 11

	radTS = radTS + NCCF(-8ALgs*2(BL(-(muu3sqrtS*3(1 + 2muu32 - x2)ELc(i,tt,gen)*
	- (muu1(-6p1k2sqrtS + sqrtS3x2)FL(i,tt,gen) + (1 - 6muu1*2)muu3sqrtS3FR(i,tt,gen))) +
	- ER(i,tt,gen)((muu12sqrtS*6(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2) +
	- 2p1k2sqrtS*2(2p1k1(-1 + x2) + 2p1k2(-1 + 2muu3*2 + x2) -
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2)))FL(i,tt,gen) +
	- muu1muu3sqrtS*4(2muu32(6p1k2 + sqrtS2(1 - 2x2)) + (6p1k1 + 6*p1k2 +
	- sqrtS*2(1 - 2x1 - 2x2))(-1 + x2))FR(i,tt,gen))) +
	- BR(muu3sqrtS*3(1 + 2muu32 - x2)ER(i,tt,gen)((1 - 6muu1*2)muu3sqrtS3FL(i,tt,gen) +
	- muu1sqrtS(-6p1k2 + sqrtS2x2)*FR(i,tt,gen)) +
	- ELc(i,tt,gen)(muu1muu3sqrtS4(-2muu32(6p1k2 + sqrtS2(1 - 2*x2)) -
	- (6p1k1 + 6p1k2 + sqrtS*2(1 - 2x1 - 2x2))(-1 + x2))
	- FL(i,tt,gen) + (muu1*2sqrtS*6(-1 - 4muu34 + x1 + x2 + 2muu3*2x2 - x1*x2) +
	- 2p1k2sqrtS*2(-2p1k1(-1 + x2) - 2p1k2(-1 + 2muu32 + x2) + sqrtS2(-1 + x2)*
	- (-1 + 2muu32 + x1 + x2)))FR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x2)2(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2))


	** TSchannel bremsstrahlung 12

	radTS = radTS + NCCF(4ALgs*2(BL(ER(i,tt,gen)(sqrtS*2(4p1k22x1 -
	- 2muu32(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) - 2p1k2sqrtS*2(-1 + x2)) +
	- 4p1k12x2 + sqrtS*4x1(-1 + x1 + x2) + 2muu1*2sqrtS*4
	- (1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) -
	- 2p1k2sqrtS*2(1 + x1*2 - x2 + x1(-3 + 2x2)) + 2p1k1(2p1k2(-2 + x1 + x2) - sqrtS2
	- (-1 + x1 + x2 + x1x2)))FL(i,tt,gen) +
	- muu1muu3sqrtS*4(6p1k2x1 + 2p1k1(-4 + x2) + sqrtS*2(-1 + x1(3 - 2x2) + x2) + 4muu32
	- (2p1k1 - 4p1k2 + sqrtS*2(-x1 + x2)))*
	- FR(i,tt,gen)) + muu3sqrtSELc(i,tt,gen)(muu1sqrtS*3
	- (-2p1k2(-4 + 8muu32 + x1) + p1k1(4 + 8muu32 - 4x1 - 6x2) + sqrtS2(4muu32(-1 + x2) +
	- (-1 + 2x1)(-1 + x1 + x2)))FL(i,tt,gen) - 2muu3sqrtS(4p1k12 + 4p1k2**2 -
	- 2p1k2sqrtS*2(-1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + x1 - 3muu1*2x1 + x2 - 3muu12x2) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)))*FR(i,tt,gen))) +
	- BR(muu3sqrtSER(i,tt,gen)(2muu3sqrtS(4p1k1*2 + 4p1k2*2 - 2p1k2sqrtS2(-1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + x1 - 3muu1*2x1 + x2 - 3muu12x2) + 2p1k1
	- (4p1k2 - sqrtS2(-1 + x1 + x2)))*FL(i,tt,gen) +
	- muu1sqrtS3(2p1k2(-4 + x1) - 4muu32(2p1k1 - 4p1k2 + sqrtS*2(-1 + x2)) - sqrtS*2(-1 + 2x1)
	- (-1 + x1 + x2) + p1k1(-4 + 4x1 + 6x2))FR(i,tt,gen)) + ELc(i,tt,gen)*
	- (muu1muu3sqrtS*4(sqrtS*2 - 6p1k2x1 - 3sqrtS*2x1 + 4muu32(4p1k2 + sqrtS2(x1 - x2)) -
	- sqrtS*2x2 + 2sqrtS2x1x2 - 2p1k1(-4 + 4muu3*2 + x2))FL(i,tt,gen) +
	- sqrtS*2(-4p1k22x1 + 2muu32(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) -
	- 2p1k2sqrtS*2(-1 + x2)) - 4p1k12x2 - sqrtS*4x1*(-1 + x1 + x2) +
	- 2muu12sqrtS*4(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)) +
	- 2p1k2sqrtS*2(1 + x1*2 - x2 + x1(-3 + 2x2)) - 2p1k1(2p1k2(-2 + x1 + x2) - sqrtS2
	- (-1 + x1 + x2 + x1x2)))FR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x1)(-1 + x2)(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)*2))


	** TSchannel bremsstrahlung 21

	radTS = radTS + NCCF(4ALgs*2(BL(ER(i,tt,gen)(sqrtS*2(4p1k22x1 -
	- 2muu32(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) - 2p1k2sqrtS*2(-1 + x2)) +
	- 4p1k12x2 +
	- sqrtS*4x1(-1 + x1 + x2) + 2muu1*2sqrtS*4(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) -
	- 2p1k2sqrtS*2(1 + x1*2 - x2 + x1(-3 + 2x2)) + 2p1k1(2p1k2*(-2 + x1 + x2) -
	- sqrtS*2(-1 + x1 + x2 + x1x2)))FL(i,tt,gen) +
	- muu1muu3sqrtS*4(2p1k1(-4 + 8muu32 + x2) + p1k2(-4 - 8muu32 + 6x1 + 4*x2) -
	- sqrtS*2(-1 - 3x1 + 4muu3*2(1 + x1 - x2) + x2 + 2x1x2))*FR(i,tt,gen)) +
	- muu3sqrtSELc(i,tt,gen)(muu1sqrtS*3(-2p1k2(-4 + x1) + 4muu32(4p1k1 - 2p1k2 - sqrtS*2x1) -
	- 6p1k1x2 +
	- sqrtS*2(-1 + x1 - x2 + 2x1x2))*FL(i,tt,gen) +
	- 2muu3sqrtS(-4p1k1*2 - 4p1k2*2 + 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(1 + (2 - 12muu12)muu3*2 - 2x1 + 3muu12x1 - 2x2 + 3muu1*2x2) + 2p1k1
	- (-4p1k2 + sqrtS2(1 + x1 + x2)))*FR(i,tt,gen))) +
	- BR(muu3sqrtSER(i,tt,gen)(2muu3sqrtS(4p1k1*2 + 4p1k2*2 - 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + 2x1 - 3muu12x1 + 2x2 - 3muu1*2x2) +
	- p1k1(8p1k2 - 2sqrtS2(1 + x1 + x2)))*FL(i,tt,gen)
	- + muu1sqrtS3(2p1k2(-4 + x1) + muu3*2(-16p1k1 + 8p1k2 + 4sqrtS2x1) + 6p1k1x2 +
	- sqrtS*2(1 + x2 - x1(1 + 2x2)))*FR(i,tt,gen)) +
	- ELc(i,tt,gen)(muu1muu3sqrtS4(p1k2(4 + 8muu3*2 - 6x1 - 4x2) - 2p1k1(-4 + 8muu3**2 + x2) +
	- sqrtS*2(-1 - 3x1 + 4muu3*2(1 + x1 - x2) + x2 + 2x1x2))*FL(i,tt,gen) +
	- sqrtS*2(-4p1k22x1 + 2muu32(4p1k12 + 4p1k22 + sqrtS4x1 - 2p1k1sqrtS2(1 + x1) -
	- 2p1k2sqrtS*2(-1 + x2)) - 4p1k12x2 -
	- sqrtS*4x1(-1 + x1 + x2) + 2muu1*2sqrtS*4(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)) +
	- 2p1k2sqrtS*2(1 + x1*2 - x2 + x1(-3 + 2x2)) - 2p1k1(2p1k2*(-2 + x1 + x2) -
	- sqrtS*2(-1 + x1 + x2 + x1x2)))FR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(2p1k1 + sqrtS*2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)(-1 + x1)*(-1 + x2))



	** TSchannel bremsstrahlung 22

	radTS = radTS + NCCF(8ALgs*2(BL(muu3sqrtSELc(i,tt,gen)(muu1sqrtS3
	- (2muu32(6p1k1 - sqrtS2(1 + x1)) - (-1 + x1)*
	- (-6p1k1 - 6p1k2 + sqrtS*2(1 + x1 + x2)))*FL(i,tt,gen) +
	- (-1 + 6muu12)muu3sqrtS5(-1 - 2muu32 + x1)FR(i,tt,gen)) +
	- ER(i,tt,gen)((-(sqrtS2(muu3*2(4p1k1 - 2sqrtS*2) + (2p1k1 + 2p1k2 - sqrtS2)(-1 + x1))*
	- (2p1k1 - sqrtS2x1)) +
	- muu1*2sqrtS*6(-1 - 4muu34 + x1 + 2muu3*2x1 + x2 - x1x2))FL(i,tt,gen) -
	- 2muu1muu3sqrtS4(1 + 2muu32 - x1)(-3p1k1 + sqrtS2x1)*FR(i,tt,gen))) +
	- BR(ELc(i,tt,gen)(2muu1muu3sqrtS4(1 + 2muu32 - x1)(-3p1k1 + sqrtS2x1)*FL(i,tt,gen) +
	- (sqrtS*2(muu3*2(4p1k1 - 2sqrtS*2) + (2p1k1 + 2p1k2 - sqrtS2)(-1 + x1))*
	- (2p1k1 - sqrtS2x1) +
	- muu1*2sqrtS*6(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2)))FR(i,tt,gen)) +
	- muu3sqrtSER(i,tt,gen)((-1 + 6muu1*2)muu3sqrtS5(1 + 2muu32 - x1)FL(i,tt,gen) +
	- muu1sqrtS3(2muu32(-6p1k1 + sqrtS2(1 + x1)) + (-1 + x1)(-6p1k1 - 6*p1k2 +
	- sqrtS*2(1 + x1 + x2)))*FR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x1)2(2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)**2))



	** TSchannel bremsstrahlung 31

	radTS = radTS + NCCF(-4ALgs*2(-(BL(muu3sqrtSELc(i,tt,gen)(-(muu1sqrtS
	- (8p1k22 + 2p1k1(-8p1k2 + sqrtS*2(2muu32 + x2)) + 2p1k2sqrtS2
	- (-10muu32 + 3(-1 + x1 + x2)) +
	- sqrtS*4(1 + 2muu12(-1 + x1) - x1 - x22 + muu32(2 - 2x1 + 2x2)))FL(i,tt,gen)) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + 2muu1*2(6p1k1 - 6p1k2 +
	- sqrtS*2(5 + 12muu32 - 4x1 - 4x2)) - 2p1k1*x2 +
	- sqrtS*2(-2 + x1 + 2x2))FR(i,tt,gen)) +
	- ER(i,tt,gen)((2p1k2(2muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + x1 - x2)) +
	- (2p1k1 - sqrtS2x1)(4p1k2 - sqrtS*2(-1 + x2))) +
	- muu1*2sqrtS*4(8muu34sqrtS*2 + 2muu3*2(2p1k1 - 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) +
	- 2p1k2(2 + x1 - 2*x2) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2))))*FL(i,tt,gen) +
	- muu1muu3sqrtS*2(2p1k1(12p1k2 + sqrtS2(1 - 2*x2)) +
	- 2muu32sqrtS*2(-6p1k1 + 6p1k2 + sqrtS*2(1 + 2x1 - 2x2)) +
	- sqrtS*2(-10p1k2x1 + sqrtS*2(-1 - x1 - 6muu12(-1 + x2) + x2 + 2x1x2)))*FR(i,tt,gen)))) +
	- BR(muu3sqrtSER(i,tt,gen)(muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + 2muu1*2
	- (6p1k1 - 6p1k2 + sqrtS*2(5 + 12muu32 - 4x1 - 4*x2)) -
	- 2p1k1x2 + sqrtS*2(-2 + x1 + 2x2))FL(i,tt,gen) -
	- muu1sqrtS(8p1k22 + 2p1k1(-8p1k2 + sqrtS*2(2muu32 + x2)) + 2p1k2sqrtS2
	- (-10muu32 + 3(-1 + x1 + x2)) +
	- sqrtS*4(1 + 2muu12(-1 + x1) - x1 - x22 + muu32(2 - 2x1 + 2x2)))FR(i,tt,gen)) +
	- ELc(i,tt,gen)(muu1muu3sqrtS2(2p1k1(12p1k2 + sqrtS2(1 - 2*x2)) +
	- 2muu32sqrtS*2(-6p1k1 + 6p1k2 + sqrtS*2(1 + 2x1 - 2x2)) +
	- sqrtS*2(-10p1k2x1 + sqrtS*2(-1 - x1 - 6muu12(-1 + x2) + x2 + 2x1x2)))*FL(i,tt,gen) +
	- (2p1k2(2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + x1 - x2)) + (2p1k1 - sqrtS2x1)*
	- (4p1k2 - sqrtS2(-1 + x2))) +
	- muu1*2sqrtS*4(8muu34sqrtS*2 + 2muu3*2(2p1k1 - 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) +
	- 2p1k2(2 + x1 - 2*x2) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2))))*FR(i,tt,gen)))))/
	- (sqrtS*2(-DenSmZ1)(-1 + x2)(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)**2))



	** TSchannel bremsstrahlung 32

	radTS = radTS + NCCF(-4ALgs*2(BL(ER(i,tt,gen)((-2(2p1k1 - sqrtS*2x1)*
	- (-(muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2)) + p1k2(-4p1k1 + sqrtS2(1 + x1))) -
	- muu1*2sqrtS*4(8muu34sqrtS*2 + p1k1(-4 + 4muu32 + 4x1 - 2x2) + x1(2*p1k2 +
	- sqrtS*2(1 - x1 + x2)) -
	- 2muu32(2p1k2 + sqrtS2(-1 + 3x1 + x2))))FL(i,tt,gen) +
	- muu1muu3sqrtS*2(8p1k12 + 2muu3*2sqrtS*2(10p1k1 - 2p1k2 + sqrtS*2(1 - 4*x1)) -
	- 2p1k1(8p1k2 + sqrtS2(-3 + 7*x1 - x2)) +
	- sqrtS*2(6p1k2x1 + sqrtS*2(1 - 3x1 + 4x1*2 + 2muu1*2(-1 + x2) - x2)))*FR(i,tt,gen)) -
	- muu3sqrtSELc(i,tt,gen)(muu1sqrtS(p1k1(24p1k2 - 2sqrtS*2x2) -
	- 2muu32sqrtS*2(6p1k1 - 6p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2(1 + 4x1) + sqrtS2(-1 - 6muu12(-1 + x1) + x1 + x1x2)))FL(i,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + 2muu1*2(6p1k1 - 6p1k2 +
	- sqrtS*2(5 + 12muu32 - 7x1 - x2)) - 2p1k1x2 +
	- sqrtS*2(-2 + 2x1 + x2))FR(i,tt,gen))) + BR*
	- (muu3sqrtSER(i,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 + 2p1k2*x1 +
	- 2muu12(6p1k1 - 6p1k2 + sqrtS*2(5 + 12muu32 - 7x1 - x2)) -
	- 2p1k1x2 + sqrtS*2(-2 + 2x1 + x2))FL(i,tt,gen) +
	- muu1sqrtS(p1k1(24p1k2 - 2sqrtS2x2) - 2muu32sqrtS*2(6p1k1 - 6p1k2 +
	- sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2(1 + 4x1) + sqrtS2(-1 - 6muu12(-1 + x1) + x1 + x1x2)))FR(i,tt,gen)) +
	- ELc(i,tt,gen)(-(muu1muu3sqrtS2(8p1k12 + 2muu3*2sqrtS*2
	- (10p1k1 - 2p1k2 + sqrtS*2(1 - 4*x1)) -
	- 2p1k1(8p1k2 + sqrtS2(-3 + 7x1 - x2)) + sqrtS2(6p1k2x1 +
	- sqrtS*2(1 - 3x1 + 4x1*2 + 2muu1*2(-1 + x2) - x2)))*FL(i,tt,gen)) +
	- (2(2p1k1 - sqrtS*2x1)(-(muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2)) +
	- p1k2(-4p1k1 + sqrtS*2(1 + x1))) +
	- muu1*2sqrtS*4(8muu34sqrtS*2 + p1k1(-4 + 4muu32 + 4x1 - 2*x2) +
	- x1(2p1k2 + sqrtS*2(1 - x1 + x2)) -
	- 2muu32(2p1k2 + sqrtS2(-1 + 3x1 + x2))))FR(i,tt,gen)))))/
	- (sqrtS*2(-DenSmZ1)(-1 + x1)(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)**2))




	** TSchannel dipole with Z 312

	radTS = radTS - NCCF(-8ALgs*2((-1 + x2)(2 + 2dsqrt((1 - 4muu32)/(-4muu32 + x22)) + x1*2
	- dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- 3x1(-1 + x2)dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 5x2dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 2x22dsqrt((1 - 4muu32)/(-4muu32 + x22))) +
	- 4muu34(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) + x1dsqrt((1 - 4muu3**2)/
	- (-4muu32 + x2*2)) +
	- x2(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x22)))) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 2x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 3x22dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- x2(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x22)) + 3x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x22)))))
	- (-8(-1 + 4muu3*2)p1k2*2(BLER(i,tt,gen)FL(i,tt,gen) - BRELc(i,tt,gen)FR(i,tt,gen)) +
	- 4p1k2sqrtS*2(-(BR(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)ER(i,tt,gen)*FR(i,tt,gen) +
	- ELc(i,tt,gen)(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)FL(i,tt,gen) -
	- x2FR(i,tt,gen) + 4muu3*2x2*FR(i,tt,gen) +
	- dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)FR(i,tt,gen)))) +
	- BL(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)ELc(i,tt,gen)*FL(i,tt,gen) +
	- ER(i,tt,gen)(-(x2FL(i,tt,gen)) + 4muu32x2FL(i,tt,gen) + dsqrt(1 - 4muu3*2)
	- dsqrt(-4muu32 + x22)FL(i,tt,gen) +
	- 3muu1muu3dsqrt(1 - 4muu3*2)dsqrt(-4muu32 + x22)FR(i,tt,gen)))) +
	- sqrtS*4(BR(ELc(i,tt,gen)(muu1muu3(4muu32 - x22 + 3dsqrt(1 - 4muu32)x2*
	- dsqrt(-4muu32 + x22))FL(i,tt,gen) +
	- (2muu32(1 + muu1*2(-4 + 8muu32)) - (-1 + 2muu1*2)(-1 + 2muu32)x2**2 +
	- dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*
	- FR(i,tt,gen)) + muu3ER(i,tt,gen)(-2(-1 + 6muu1*2)muu3(-4muu32 + x22)*FL(i,tt,gen) +
	- muu1(4muu32 - x22 + 3dsqrt(1 - 4muu3*2)x2dsqrt(-4muu32 + x22))*FR(i,tt,gen))) -
	- BL(muu3ELc(i,tt,gen)(muu1(4muu32 - x22 + 3dsqrt(1 - 4muu32)x2*
	- dsqrt(-4muu32 + x22))FL(i,tt,gen) -
	- 2(-1 + 6muu1*2)muu3(-4muu32 + x22)*FR(i,tt,gen)) +
	- ER(i,tt,gen)((2muu3*2(1 + muu1*2(-4 + 8muu32)) - (-1 + 2muu1*2)(-1 + 2muu32)x2**2 +
	- dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*FL(i,tt,gen) +
	- muu1muu3(4muu32 - x22 + 3dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*FR(i,tt,gen))))))/
	- ((-1 + 2muu32)sqrtS*2(-DenSmZ1)(-1 + x2)2(-2 + x1 + x2)dsqrt(-4muu32 + x22)*
	- (-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + (-1 + 2muu12 + 2muu3*2)
	- dsqrt(-4muu32 + x22)) - 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)**2))


	** TSchannel dipole with Z 321

	radTS = radTS - NCCF(8ALgs*2(4muu34(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) +
	- x1(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x12))) +
	- dsqrt((1 - 4muu32)/(-4muu32 + x12))x2) + (-1 + x1)
	- (2 + 2dsqrt((-1 + 4muu3*2)/(4muu32 - x12)) + 2x12dsqrt((1 - 4muu32)/(-4muu32 + x12)) -
	- 3dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2 + dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2*2 +
	- x1dsqrt((1 - 4muu3*2)/(-4muu32 + x12))(-5 + 3x2)) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3x12dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x1*2)) -
	- 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12))*x2 +
	- x1(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x12)) + 3dsqrt((1 - 4muu3*2)/(-4muu32 + x12))x2)))
	- (8(-1 + 4muu3*2)p1k1*2(BLER(i,tt,gen)FL(i,tt,gen) - BRELc(i,tt,gen)FR(i,tt,gen)) +
	- 4p1k1sqrtS*2(BL(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)ELc(i,tt,gen)*FL(i,tt,gen) +
	- ER(i,tt,gen)((x1 - 4muu3*2x1 + dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- 3muu1muu3dsqrt(1 - 4muu3*2)dsqrt(-4muu32 + x12)FR(i,tt,gen))) -
	- BR(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)ER(i,tt,gen)*FR(i,tt,gen) +
	- ELc(i,tt,gen)(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)FL(i,tt,gen) +
	- (x1 - 4muu32x1 + dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12))FR(i,tt,gen)))) +
	- sqrtS*4(BR(ELc(i,tt,gen)(muu1muu3(-4muu32 + x12 + 3dsqrt(1 - 4muu32)x1*
	- dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- (-2muu32(1 + muu1*2(-4 + 8muu32)) + (-1 + 2muu1*2)(-1 + 2muu32)x1**2 +
	- dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*
	- FR(i,tt,gen)) + muu3ER(i,tt,gen)(2(-1 + 6muu1*2)muu3(-4muu32 + x12)*FL(i,tt,gen) +
	- muu1(-4muu32 + x12 + 3dsqrt(1 - 4muu3*2)x1dsqrt(-4muu32 + x12))*FR(i,tt,gen))) -
	- BL(muu3ELc(i,tt,gen)(muu1(-4muu32 + x12 + 3dsqrt(1 - 4muu32)x1*
	- dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- 2(-1 + 6muu1*2)muu3(-4muu32 + x12)*FR(i,tt,gen)) +
	- ER(i,tt,gen)((-2muu3*2(1 + muu1*2(-4 + 8muu32)) + (-1 + 2muu1*2)(-1 + 2muu32)x1**2 +
	- dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*FL(i,tt,gen) +
	- muu1muu3(-4muu32 + x12 + 3dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*FR(i,tt,gen))))))/
	- ((-1 + 2muu32)sqrtS*2(-DenSmZ1)(-1 + x1)2dsqrt(-4muu32 + x12)(-2 + x1 + x2)*
	- (4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + (-1 + 2muu12 + 2muu3*2)
	- dsqrt(-4muu32 + x12)) - 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)**2))


	do j = 1, 4

	** TSchannel bremsstrahlung 13

	radTS = radTS + NCCF(-2gs2(CL(j)(2DR(j)(muu1sqrtSER(i,tt,gen)(muu3sqrtS3
	- (2p1k1(-1 + x2) + (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))FL(i,tt,gen) +
	- muu1sqrtS5(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2)*FR(i,tt,gen)) -
	- muu3sqrtS3(1 + 2muu32 - x2)ELc(i,tt,gen)((1 - 2muu1*2)muu3sqrtS3FL(i,tt,gen) +
	- muu1sqrtS(-2p1k2 + sqrtS2x2)*FR(i,tt,gen))) +
	- sqrtS*4DL(j)(4muu1muu3(1 + 2muu32 - x2)ER(i,tt,gen)*
	- (p1k2FL(i,tt,gen) + muu1muu3sqrtS2FR(i,tt,gen)) +
	- ELc(i,tt,gen)((8muu3*2p1k2 - 4muu34sqrtS*2 + 2muu1*2sqrtS*2
	- (4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2) -
	- (-1 + x2)(2p1k2(-2 + x1) - 2p1k1x2 + sqrtS2(-1 + x1 + x2)))*FL(i,tt,gen) +
	- 2muu1muu3(2p1k1(-1 + x2) + 2p1k2(-1 + 2muu3**2 + x2) -
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))FR(i,tt,gen)))) +
	- CR(j)(2DL(j)(-(muu3sqrtS*3(1 + 2muu32 - x2)ER(i,tt,gen)*
	- (muu1(-2p1k2sqrtS + sqrtS3x2)FL(i,tt,gen) + (1 - 2muu1*2)muu3sqrtS3FR(i,tt,gen))) +
	- muu1sqrtSELc(i,tt,gen)(muu1sqrtS*5(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2)*
	- FL(i,tt,gen) +
	- muu3sqrtS3(2p1k1(-1 + x2) + (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))FR(i,tt,gen))) +
	- sqrtS*4DR(j)(4muu1muu3(1 + 2muu32 - x2)ELc(i,tt,gen)*
	- (muu1muu3sqrtS*2FL(i,tt,gen) + p1k2*FR(i,tt,gen)) +
	- ER(i,tt,gen)(2muu1muu3(2p1k1(-1 + x2) + 2p1k2(-1 + 2muu3*2 + x2) -
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))FL(i,tt,gen) +
	- (8muu32p1k2 - 4muu34sqrtS*2 + 2muu1*2sqrtS*2(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2) -
	- (-1 + x2)(2p1k2(-2 + x1) - 2p1k1x2 + sqrtS2(-1 + x1 + x2)))*FR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x2)*2(-DenSmH1(j))(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)*2))


	** TSchannel bremsstrahlung 14

	radTS = radTS + NCCF(-2gs2(CL(j)(DR(j)(muu1sqrtSER(i,tt,gen)(muu3sqrtS*3
	- (sqrtS*2(1 - 4muu32 + x1 - x2) + 2p1k1x2 + p1k2(-4 + 8muu32 - 2x1 + 4x2))FL(i,tt,gen) +
	- 2muu1sqrtS*5(4muu34 + x1x2 - muu3*2(2 + x1 + x2))*FR(i,tt,gen)) +
	- muu3sqrtS4ELc(i,tt,gen)(2muu3(2p1k1 + 2p1k2 + sqrtS2((-2 + 4muu12)muu3**2 -
	- muu1*2(x1 + x2)))*FL(i,tt,gen) +
	- muu1(-2p1k2x1 + 2p1k1x2 + sqrtS2(-1 + x1 + x2) + muu3*2(8p1k2 - 4sqrtS*2x2))*FR(i,tt,gen))) +
	- DL(j)(muu1muu3sqrtS2ER(i,tt,gen)(sqrtS2(8muu32p1k2 - 2p1k2x1 + 2p1k1x2 -
	- sqrtS*2(-1 + x1 + x2))*FL(i,tt,gen) +
	- 2muu1muu3sqrtS4(4muu32 - x1 - x2)FR(i,tt,gen)) +
	- sqrtS*4ELc(i,tt,gen)(2(-2muu34sqrtS*2 - 2p1k2x1 + muu32(2p1k1 + 6p1k2 - sqrtS*2(-1 + x2)) +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)))*FL(i,tt,gen) +
	- muu1muu3(2p1k1x2 + p1k2(-4 + 8muu3*2 - 2x1 + 4*x2) -
	- sqrtS*2(1 + x1 + 4muu32(-1 + x2) - 3x2 + 2x2*2))FR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu1sqrtSELc(i,tt,gen)(2muu1sqrtS5(4muu34 + x1x2 - muu3*2(2 + x1 + x2))*
	- FL(i,tt,gen) + muu3sqrtS3(sqrtS*2(1 - 4muu3*2 + x1 - x2) +
	- 2p1k1x2 + p1k2(-4 + 8muu3*2 - 2x1 + 4x2))FR(i,tt,gen)) +
	- muu3sqrtS4ER(i,tt,gen)(muu1(-2p1k2x1 + 2p1k1x2 + sqrtS*2(-1 + x1 + x2) +
	- muu3*2(8p1k2 - 4sqrtS*2x2))*FL(i,tt,gen) +
	- 2muu3(2p1k1 + 2p1k2 + sqrtS*2((-2 + 4muu12)muu32 - muu12(x1 + x2)))FR(i,tt,gen))) +
	- DR(j)(muu1muu3sqrtS2ELc(i,tt,gen)(2muu1muu3sqrtS*4(4muu32 - x1 - x2)FL(i,tt,gen) +
	- sqrtS*2(8muu32p1k2 - 2p1k2x1 + 2p1k1x2 - sqrtS*2(-1 + x1 + x2))*FR(i,tt,gen)) +
	- sqrtS*4ER(i,tt,gen)(muu1muu3(2p1k1x2 + p1k2(-4 + 8muu32 - 2x1 + 4*x2) -
	- sqrtS*2(1 + x1 + 4muu32(-1 + x2) - 3x2 + 2x2*2))
	- FL(i,tt,gen) + 2(-2muu3*4sqrtS*2 - 2p1k2x1 + muu32(2p1k1 + 6p1k2 - sqrtS*2(-1 + x2)) +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)))*FR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x1)(-1 + x2)(-DenSmH1(j))*
	- (-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2))



	** TSchannel bremsstrahlung 23

	radTS = radTS + NCCF(2gs2(CL(j)(DL(j)(muu1muu3sqrtS*2ER(i,tt,gen)*
	- (sqrtS*2(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2*x2) +
	- sqrtS*2(1 - 4muu32 - x1 + x2))FL(i,tt,gen) +
	- 2muu1muu3sqrtS4(-4muu32 + x1 + x2)FR(i,tt,gen)) +
	- sqrtS*4ELc(i,tt,gen)(-2(-2muu34sqrtS*2 + (2p1k1 - sqrtS*2x1)*x2 +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)) +
	- muu3*2(-6p1k1 - 2p1k2 + sqrtS*2(1 + 2x1 + x2)))FL(i,tt,gen) +
	- muu1muu3(2p1k2x1 + muu3*2(8p1k1 - 4sqrtS*2x1) - 2p1k1x2 +
	- sqrtS*2(-1 + x1 + x2))*FR(i,tt,gen))) +
	- DR(j)(muu1sqrtSER(i,tt,gen)(muu3sqrtS3(8muu32p1k1 + 2p1k2x1 - 2p1k1x2 -
	- sqrtS*2(-1 + x1 + x2))*FL(i,tt,gen) +
	- 2muu1sqrtS*5(-4muu34 - x1x2 + muu3*2(2 + x1 + x2))*FR(i,tt,gen)) +
	- muu3sqrtS4ELc(i,tt,gen)(2muu3(2p1k1 + 2p1k2 - sqrtS2((-2 + 4muu12)muu3**2 -
	- (-1 + muu1*2)(x1 + x2)))*FL(i,tt,gen) +
	- muu1(2p1k2x1 + p1k1(-4 + 8muu32 + 4x1 - 2*x2) -
	- sqrtS*2(1 + 4muu32(-1 + x1) - 3x1 + 2x1*2 + x2))FR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu1sqrtSELc(i,tt,gen)(2muu1sqrtS5(-4muu34 - x1x2 +
	- muu3*2(2 + x1 + x2))*FL(i,tt,gen) +
	- muu3sqrtS3(8muu32p1k1 + 2p1k2x1 - 2p1k1x2 - sqrtS*2(-1 + x1 + x2))*FR(i,tt,gen)) +
	- muu3sqrtS4ER(i,tt,gen)(muu1(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2*x2) -
	- sqrtS*2(1 + 4muu32(-1 + x1) - 3x1 + 2x1*2 + x2))FL(i,tt,gen) +
	- 2muu3(2p1k1 + 2p1k2 - sqrtS*2((-2 + 4muu12)muu32 - (-1 + muu12)(x1 + x2)))FR(i,tt,gen))) +
	- DR(j)(muu1muu3sqrtS2ELc(i,tt,gen)(2muu1muu3sqrtS*4(-4muu32 + x1 + x2)FL(i,tt,gen) +
	- sqrtS*2(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2*x2) +
	- sqrtS*2(1 - 4muu32 - x1 + x2))FR(i,tt,gen)) +
	- sqrtS*4ER(i,tt,gen)(muu1muu3(2p1k2x1 + muu32(8p1k1 - 4sqrtS*2x1) - 2p1k1x2 +
	- sqrtS*2(-1 + x1 + x2))*FL(i,tt,gen) -
	- 2(-2muu3*4sqrtS*2 + (2p1k1 - sqrtS*2x1)x2 + muu12sqrtS*2(4muu34 + x1x2 -
	- muu3*2(2 + x1 + x2)) +
	- muu3*2(-6p1k1 - 2p1k2 + sqrtS*2(1 + 2x1 + x2)))FR(i,tt,gen))))))/
	- (sqrtS*4(MSf(i,tt,gen)*2 - 2p1k1 - sqrtS*2(muu12 + muu32 - x1))*
	- (-1 + x1)(-1 + x2)(DenSmH1(j)))


	** TSchannel bremsstrahlung 24

	radTS = radTS + NCCF(-2gs2(CR(j)(2DL(j)(muu1sqrtSELc(i,tt,gen)
	- (muu1sqrtS5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))FL(i,tt,gen) -
	- 2muu3p1k1sqrtS3(1 + 2muu32 - x1)FR(i,tt,gen)) +
	- muu3sqrtSER(i,tt,gen)(muu1sqrtS*3(-2p1k1(-1 + 2muu3*2 + x1) +
	- (-1 + x1)(-2p1k2 + sqrtS*2(-1 + 2muu32 + x1 + x2)))FL(i,tt,gen) +
	- (-1 + 2muu12)muu3sqrtS5(1 + 2muu32 - x1)FR(i,tt,gen))) +
	- DR(j)(2muu1muu3sqrtS*2ELc(i,tt,gen)(2muu1muu3sqrtS*4(1 + 2muu32 - x1)FL(i,tt,gen) +
	- sqrtS*2(-2p1k2(-1 + x1) - 2p1k1(-1 + 2muu3*2 + x1) +
	- sqrtS*2(-1 + 2muu32 + x1))FR(i,tt,gen)) +
	- sqrtS*4ER(i,tt,gen)(2muu1muu3(1 + 2muu32 - x1)(-2p1k1 + sqrtS2x1)*FL(i,tt,gen) +
	- (-4muu34sqrtS*2 + 4muu3*2(-2p1k1 + sqrtS2x1) + (-1 + x1)(-2p1k2*x1 +
	- sqrtS*2(1 + x1 - x2) + 2p1k1(-2 + x2)) +
	- 2muu12sqrtS*2(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2)))FR(i,tt,gen)))) +
	- CL(j)(DL(j)(2muu1muu3sqrtS2ER(i,tt,gen)(sqrtS2(-2p1k2(-1 + x1) -
	- 2p1k1(-1 + 2muu32 + x1) + sqrtS2(-1 + 2muu32 + x1))
	- FL(i,tt,gen) + 2muu1muu3sqrtS4(1 + 2muu32 - x1)FR(i,tt,gen)) +
	- sqrtS*4ELc(i,tt,gen)((-4muu3*4sqrtS*2 + 4muu3*2(-2p1k1 + sqrtS2x1) +
	- (-1 + x1)(-2p1k2x1 + sqrtS2(1 + x1 - x2) + 2p1k1(-2 + x2)) + 2muu12sqrtS*2
	- (4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2)))
	- FL(i,tt,gen) + 2muu1muu3(1 + 2muu3*2 - x1)(-2p1k1 + sqrtS2x1)*FR(i,tt,gen))) +
	- 2DR(j)(muu1sqrtSER(i,tt,gen)(-2muu3p1k1sqrtS*3(1 + 2muu32 - x1)FL(i,tt,gen) +
	- muu1sqrtS5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))FR(i,tt,gen)) +
	- muu3sqrtSELc(i,tt,gen)((-1 + 2muu1*2)muu3sqrtS5(1 + 2muu32 - x1)FL(i,tt,gen) +
	- muu1sqrtS3(-2p1k1(-1 + 2muu32 + x1) + (-1 + x1)(-2*p1k2 +
	- sqrtS*2(-1 + 2muu32 + x1 + x2)))FR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x1)*2(-DenSmH1(j))*
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)**2))


	** TSchannel bremsstrahlung 33

	radTS = radTS + NCCF(gs*2(CL(j)(sqrtS2DL(j)(ELc(i,tt,gen)((-8muu34sqrtS**4 +
	- 2p1k1(8p1k2 + sqrtS2(-1 + x2)*x2) -
	- sqrtS*2(2p1k2x1(3 + x2) + sqrtS2(-1 + x2)*(-1 + x1 + x2)) -
	- muu3*2sqrtS*2(4p1k1(2 + x2) - 2(2p1k2(2 + x1) + sqrtS2(-2 + 3x1 + 2x2))) +
	- 2muu12sqrtS*2(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1*x2 +
	- muu3*2(4p1k1 - 4p1k2 - 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3x2))))FL(i,tt,gen) +
	- 2muu1muu3(2muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + x1 - x2)) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x2)) +
	- sqrtS*2(-(x1(4p1k2 - sqrtS*2(-1 + x2))) - 2muu12sqrtS*2(-1 + x2)))*FR(i,tt,gen)) -
	- 4muu1muu3ER(i,tt,gen)((-4p1k1p1k2 - muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2) +
	- sqrtS*2(muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x2)))
	- FL(i,tt,gen) + muu1muu3sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(-1 - 4muu32 + x1 + x2))FR(i,tt,gen))) +
	- 2DR(j)(muu1sqrtSER(i,tt,gen)(muu3sqrtS(p1k1(8p1k2 - 2sqrtS**2) -
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*FL(i,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 + 2p1k2x1 + 2muu3*2(2p1k1 - 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2)))*FR(i,tt,gen)) +
	- muu3sqrtSELc(i,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2*x2)) -
	- 2p1k1x2 + sqrtS*2(-2 + x1 + 2x2))FL(i,tt,gen) +
	- muu1sqrtS(p1k1(8p1k2 - 2sqrtS2x2) - 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(-2p1k2(-1 + x1 + x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2*2)))FR(i,tt,gen)))) +
	- CR(j)(2DL(j)(muu3sqrtSER(i,tt,gen)(muu1sqrtS(p1k1(8p1k2 - 2sqrtS2x2) - 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(-2p1k2(-1 + x1 + x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2*2)))FL(i,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) - 2p1k1*x2 +
	- sqrtS*2(-2 + x1 + 2x2))FR(i,tt,gen)) +
	- muu1sqrtSELc(i,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 + 2p1k2*x1 +
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2)))*FL(i,tt,gen) +
	- muu3sqrtS(p1k1(8p1k2 - 2sqrtS2) - 2muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*FR(i,tt,gen))) +
	- sqrtS*2DR(j)(-4muu1muu3ELc(i,tt,gen)(muu1muu3sqrtS2(-2p1k1 + 2p1k2 +
	- sqrtS*2(-1 - 4muu32 + x1 + x2))FL(i,tt,gen) +
	- (-4p1k1p1k2 - muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS2) + sqrtS2*
	- (muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x2)))FR(i,tt,gen)) +
	- ER(i,tt,gen)(2muu1muu3(2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + x1 - x2)) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x2)) +
	- sqrtS*2(-(x1(4p1k2 - sqrtS*2(-1 + x2))) - 2muu12sqrtS*2(-1 + x2)))*FL(i,tt,gen) +
	- (-8muu34sqrtS*4 + 2p1k1(8p1k2 + sqrtS*2(-1 + x2)x2) - sqrtS2(2p1k2x1*(3 + x2) +
	- sqrtS*2(-1 + x2)*(-1 + x1 + x2)) -
	- muu3*2sqrtS*2(4p1k1(2 + x2) - 2(2p1k2(2 + x1) + sqrtS2(-2 + 3x1 + 2x2))) +
	- 2muu12sqrtS*2(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1*x2 +
	- muu3*2(4p1k1 - 4p1k2 - 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3x2))))FR(i,tt,gen))))))/
	- (sqrtS*2(-1 + x2)(-DenSmH1(j))(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)**2))



	** TSchannel bremsstrahlung 34

	radTS = radTS + NCCF(gs*2(CL(j)(sqrtS2DL(j)(-4muu1muu3ER(i,tt,gen)*
	- ((-4p1k1p1k2 - muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS2) + sqrtS2*
	- (muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x1)))FL(i,tt,gen) -
	- muu1muu3sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - 2x1))*FR(i,tt,gen)) +
	- ELc(i,tt,gen)((-8muu3*4sqrtS*4 + 2p1k1(8p1k2 + sqrtS*2(-1 + x1)*x2) -
	- sqrtS*2(2p1k2x1(3 + x1) + sqrtS2(-1 + x1)*(-1 + x1 + x2)) -
	- muu3*2sqrtS*2(4p1k1(2 + x2) - 2(2p1k2(2 + x1) + sqrtS2(-2 + 4*x1 + x2))) +
	- 2muu12sqrtS*2(8muu34sqrtS*2 + 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) -
	- 2p1k1x2 + muu3*2(4p1k1 - 4p1k2 - 2sqrtS2(-3 + 4x1 + 2x2))))*FL(i,tt,gen) -
	- 2muu1muu3(-2p1k1(4p1k2 + sqrtS*2(1 + x1 - 2x2)) + sqrtS2(x1(4p1k2 +
	- sqrtS*2(1 + x1 - 2x2)) + 2muu1*2sqrtS*2(-1 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)))*FR(i,tt,gen))) +
	- 2DR(j)(muu1sqrtSER(i,tt,gen)(muu3sqrtS(-2muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2) +
	- 2p1k1(4p1k2 + sqrtS2(-1 + x1 - x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*FL(i,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1x2 + sqrtS2x1*(-1 + x1 + x2) +
	- 2muu32(2p1k1 - 2p1k2 - sqrtS*2(-1 + 3x1 + x2)))FR(i,tt,gen)) +
	- muu3sqrtSELc(i,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 3x1 - x2)) -
	- 2p1k1x2 + sqrtS*2(-2 + 2x1 + x2))FL(i,tt,gen) +
	- muu1sqrtS(p1k1(8p1k2 - 2sqrtS2x2) - 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(p1k2(2 - 4x1) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x1x2)))FR(i,tt,gen)))) +
	- CR(j)(2DL(j)(muu3sqrtSER(i,tt,gen)(muu1sqrtS(p1k1(8p1k2 - 2sqrtS2x2) - 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(p1k2(2 - 4x1) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x1x2)))FL(i,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 3x1 - x2)) - 2p1k1x2 +
	- sqrtS*2(-2 + 2x1 + x2))FR(i,tt,gen)) +
	- muu1sqrtSELc(i,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1*x2 +
	- sqrtS*2x1*(-1 + x1 + x2) +
	- 2muu32(2p1k1 - 2p1k2 - sqrtS*2(-1 + 3x1 + x2)))FL(i,tt,gen) +
	- muu3sqrtS(-2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2) + 2p1k1(4p1k2 + sqrtS*2(-1 + x1 - x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*FR(i,tt,gen))) +
	- sqrtS*2DR(j)(4muu1muu3ELc(i,tt,gen)(muu1muu3sqrtS2(2p1k1 - 2p1k2 +
	- sqrtS*2(1 + 4muu32 - 2x1))*FL(i,tt,gen) +
	- (4p1k1p1k2 + muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2) -
	- sqrtS*2(muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x1)))FR(i,tt,gen)) +
	- ER(i,tt,gen)(-2muu1muu3(-2p1k1(4p1k2 + sqrtS2(1 + x1 - 2*x2)) +
	- sqrtS*2(x1(4p1k2 + sqrtS*2(1 + x1 - 2x2)) + 2muu1*2sqrtS*2(-1 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)))*FL(i,tt,gen) +
	- (-8muu34sqrtS*4 + 2p1k1(8p1k2 + sqrtS*2(-1 + x1)*x2) -
	- sqrtS*2(2p1k2x1(3 + x1) + sqrtS2(-1 + x1)*(-1 + x1 + x2)) -
	- muu3*2sqrtS*2(4p1k1(2 + x2) - 2(2p1k2(2 + x1) + sqrtS2(-2 + 4*x1 + x2))) +
	- 2muu12sqrtS*2(8muu34sqrtS*2 + 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) -
	- 2p1k1x2 + muu3*2(4p1k1 - 4p1k2 - 2sqrtS2(-3 + 4x1 + 2x2))))*FR(i,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(DenSmH1(j))(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (-2p1k1 - sqrtS2(muu12 + muu32 - x1) + MSf(i,tt,gen)**2))

	** TSchannel dipole from Higgs 312

	radTS = radTS - (-4CFGS*2NC((-1 + x2)(2 + 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) +
	- x1*2dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- 3x1(-1 + x2)dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 5x2dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 2x22dsqrt((1 - 4muu32)/(-4muu32 + x22))) +
	- 4muu34(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) + x1dsqrt((1 - 4muu3**2)/
	- (-4muu32 + x2*2)) +
	- x2(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x22)))) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 2x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 3x22dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- x2(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x22)) + 3x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x22)))))
	- (-4dsqrt(1 - 4muu3*2)p1k2(CL(j)(muu1muu3DR(j)(ER(i,tt,gen)FL(i,tt,gen) +
	- ELc(i,tt,gen)*FR(i,tt,gen)) +
	- DL(j)(ELc(i,tt,gen)FL(i,tt,gen) + muu1muu3ER(i,tt,gen)*FL(i,tt,gen) +
	- muu1muu3ELc(i,tt,gen)*FR(i,tt,gen))) +
	- CR(j)(muu1muu3DL(j)(ER(i,tt,gen)FL(i,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)) +
	- DR(j)(muu1muu3ER(i,tt,gen)FL(i,tt,gen) + muu1muu3ELc(i,tt,gen)*FR(i,tt,gen) +
	- ER(i,tt,gen)*FR(i,tt,gen)))) +
	- sqrtS*2(CL(j)(DL(j)(muu1muu3ER(i,tt,gen)(dsqrt(1 - 4muu3*2)x2*FL(i,tt,gen) -
	- dsqrt(-4muu32 + x22)FL(i,tt,gen) -
	- 4muu1muu3dsqrt(-4muu32 + x22)*FR(i,tt,gen)) +
	- ELc(i,tt,gen)((dsqrt(1 - 4muu3*2)x2 - (-1 + 2muu12)(-1 + 2muu32)
	- dsqrt(-4muu32 + x22))FL(i,tt,gen) +
	- muu1muu3(dsqrt(1 - 4muu32)x2 - dsqrt(-4muu32 + x22))FR(i,tt,gen))) +
	- DR(j)(muu1ER(i,tt,gen)(muu3(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22))FL(i,tt,gen) +
	- 2muu1(1 - 2muu32)dsqrt(-4muu32 + x22)FR(i,tt,gen)) +
	- muu3ELc(i,tt,gen)(2(1 - 2muu1*2)muu3dsqrt(-4muu32 + x22)*FL(i,tt,gen) +
	- muu1(dsqrt(1 - 4muu3*2)x2 + dsqrt(-4muu32 + x22))FR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu3ER(i,tt,gen)(muu1(dsqrt(1 - 4muu3*2)x2 +
	- dsqrt(-4muu32 + x22))FL(i,tt,gen) +
	- 2(1 - 2muu1*2)muu3dsqrt(-4muu32 + x22)*FR(i,tt,gen)) +
	- muu1ELc(i,tt,gen)(2muu1(1 - 2muu32)dsqrt(-4muu32 + x22)FL(i,tt,gen) +
	- muu3(dsqrt(1 - 4muu3*2)x2 + dsqrt(-4muu32 + x22))FR(i,tt,gen))) +
	- DR(j)(muu1muu3ELc(i,tt,gen)(-4muu1muu3dsqrt(-4muu32 + x22)*FL(i,tt,gen) +
	- dsqrt(1 - 4muu32)x2*FR(i,tt,gen) -
	- dsqrt(-4muu32 + x22)FR(i,tt,gen)) +
	- ER(i,tt,gen)(muu1muu3(dsqrt(1 - 4muu3*2)x2 - dsqrt(-4muu32 + x22))FL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x2 - (-1 + 2muu12)(-1 + 2muu32)
	- dsqrt(-4muu32 + x22))FR(i,tt,gen)))))))/
	- ((-1 + 2muu32)(-1 + x2)*2(-2 + x1 + x2)(DenSmH1(j))
	- (-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 +
	- (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x2*2)) -
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)*2))


	** TSchannel dipole from Higgs 321

	radTS = radTS - (-4CFGS*2NC(4muu3*4(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) +
	- x1(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x12))) +
	- dsqrt((1 - 4muu32)/(-4muu32 + x12))x2) + (-1 + x1)
	- (2 + 2dsqrt((-1 + 4muu3*2)/(4muu32 - x12)) + 2x12dsqrt((1 - 4muu32)/(-4muu32 + x12)) -
	- 3dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2 + dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2*2 +
	- x1dsqrt((1 - 4muu3*2)/(-4muu32 + x12))(-5 + 3x2)) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3x12dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x1*2)) -
	- 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12))*x2 +
	- x1(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x12)) + 3dsqrt((1 - 4muu3**2)/
	- (-4muu32 + x12))x2)))*
	- (-4dsqrt(1 - 4muu3*2)p1k1(CL(j)(muu1muu3DR(j)(ER(i,tt,gen)FL(i,tt,gen) +
	- ELc(i,tt,gen)*FR(i,tt,gen)) +
	- DL(j)(ELc(i,tt,gen)FL(i,tt,gen) + muu1muu3ER(i,tt,gen)*FL(i,tt,gen) +
	- muu1muu3ELc(i,tt,gen)*FR(i,tt,gen))) +
	- CR(j)(muu1muu3DL(j)(ER(i,tt,gen)FL(i,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)) +
	- DR(j)(muu1muu3ER(i,tt,gen)FL(i,tt,gen) + muu1muu3ELc(i,tt,gen)*FR(i,tt,gen) +
	- ER(i,tt,gen)*FR(i,tt,gen)))) +
	- sqrtS*2(CL(j)(DL(j)(muu1muu3ER(i,tt,gen)((dsqrt(1 - 4muu3*2)x1 +
	- dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- 4muu1muu3dsqrt(-4muu32 + x12)*FR(i,tt,gen)) +
	- ELc(i,tt,gen)((dsqrt(1 - 4muu3*2)x1 + (-1 + 2muu12)(-1 + 2muu32)
	- dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- muu1muu3(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12))FR(i,tt,gen))) +
	- DR(j)(muu1ER(i,tt,gen)((muu3dsqrt(1 - 4muu32)x1 -
	- muu3dsqrt(-4muu32 + x12))*FL(i,tt,gen) +
	- 2muu1(-1 + 2muu32)dsqrt(-4muu32 + x12)FR(i,tt,gen)) +
	- muu3ELc(i,tt,gen)(2(-1 + 2muu1*2)muu3dsqrt(-4muu32 + x12)*FL(i,tt,gen) +
	- (muu1dsqrt(1 - 4muu3*2)x1 - muu1dsqrt(-4muu32 + x12))*FR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu3ER(i,tt,gen)((muu1dsqrt(1 - 4muu3*2)x1 -
	- muu1dsqrt(-4muu32 + x12))*FL(i,tt,gen) +
	- 2(-1 + 2muu1*2)muu3dsqrt(-4muu32 + x12)*FR(i,tt,gen)) +
	- muu1ELc(i,tt,gen)(2muu1(-1 + 2muu32)dsqrt(-4muu32 + x12)FL(i,tt,gen) +
	- (muu3dsqrt(1 - 4muu3*2)x1 - muu3dsqrt(-4muu32 + x12))*FR(i,tt,gen))) +
	- DR(j)(muu1muu3ELc(i,tt,gen)(4muu1muu3dsqrt(-4muu32 + x12)*FL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12))FR(i,tt,gen)) +
	- ER(i,tt,gen)(muu1muu3(dsqrt(1 - 4muu3*2)x1 + dsqrt(-4muu32 + x12))FL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x1 + (-1 + 2muu12)(-1 + 2muu32)
	- dsqrt(-4muu32 + x12))FR(i,tt,gen)))))))/
	- ((-1 + 2muu32)(-1 + x1)*2(-2 + x1 + x2)(DenSmH1(j))
	- (-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + (1 - 2muu12 - 2muu3*2)
	- dsqrt(-4muu32 + x1*2)) +
	- 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)*2))

	enddo
	enddo

	end


	************************************************

	double precision function radUS()

	#include "mykinematics.h"
	#include "model.h"
	#include "gencoupl.h"

	integer i,j
	complex*16 DenSmZ, DenSmH(4)
	double precision DenSmZ1, DenSmH1(4)

	DenSmZ = (sqrtS*2 - MZ2) + cIMZ*wZ
	DenSmZ1 = DenSmZdconjg(DenSmZ)/(sqrtS*2 - MZ2)

	do i = 1, 4
	DenSmH(i) = (sqrtS2 - MHk0(i)2) + cIMHk0(i)wHk0(i)
	DenSmH1(i) = DenSmH(i)dconjg(DenSmH(i))/(sqrtS2 - MHk0(i)*2)
	enddo

	radUS = 0d0

	do i = 1, 2

	** USchannel bremsstrahlung 11

	radUS = radUS + NCCF(8ALgs*2(BR(muu3sqrtSGL(i,tt,gen)(muu1sqrtS3
	- (2muu32(6p1k2 - sqrtS2(1 + x2)) - (-1 + x2)(-6p1k1 - 6p1k2 + sqrtS2(1 + x1 + x2)))*HL(i,tt,gen) +
	- (-1 + 6muu12)muu3sqrtS5(-1 - 2muu32 + x2)HR(i,tt,gen)) +
	- GR(i,tt,gen)((muu12sqrtS*6(-1 - 4muu34 + x1 + x2 + 2muu3*2x2 - x1*x2) -
	- sqrtS*2(2p1k2 - sqrtS2x2)(2p1k1(-1 + x2) + (2p1k2 - sqrtS*2)(-1 + 2muu32 + x2)))HL(i,tt,gen) -
	- 2muu1muu3sqrtS4(1 + 2muu32 - x2)(-3p1k2 + sqrtS2x2)*HR(i,tt,gen))) +
	- BL(GL(i,tt,gen)(2muu1muu3sqrtS4(1 + 2muu32 - x2)(-3p1k2 + sqrtS2x2)*HL(i,tt,gen) +
	- (muu1*2sqrtS*6(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2) +
	- sqrtS*2(2p1k2 - sqrtS2x2)(2p1k1(-1 + x2) + (2p1k2 - sqrtS*2)(-1 + 2muu32 + x2)))HR(i,tt,gen)) +
	- muu3sqrtSGR(i,tt,gen)((-1 + 6muu1*2)muu3sqrtS5(1 + 2muu32 - x2)HL(i,tt,gen) +
	- muu1sqrtS3(2muu32(-6p1k2 + sqrtS2(1 + x2)) + (-1 + x2)(-6p1k1 - 6p1k2 + sqrtS2(1 + x1 + x2)))*
	- HR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x2)2(2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))

	** USchannel bremsstrahlung 12

	radUS = radUS + NCCF(4ALgs*2(BR(GR(i,tt,gen)(sqrtS*2(4p1k12x2 + (2p1k2 - sqrtS2x2)*
	- (2p1k2x1 - sqrtS*2(-1 + x1 + x2)) +
	- 2muu12sqrtS*4(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) -
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) +
	- 2p1k1(2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + 3x2 - 2x1x2 - x22)))HL(i,tt,gen) +
	- muu1muu3sqrtS*4(2p1k2(-4 + 8muu32 + x1) + p1k1(-4 - 8muu32 + 4x1 + 6*x2) +
	- sqrtS*2(1 - x1 + 4muu32(-1 + x1 - x2) + 3x2 - 2x1x2))HR(i,tt,gen)) +
	- muu3sqrtSGL(i,tt,gen)(-(muu1sqrtS*3(6p1k2x1 + 2p1k1(-4 + x2) + 4muu32(2p1k1 - 4p1k2 + sqrtS*2x2) +
	- sqrtS*2(1 + x1 - x2 - 2x1x2))*HL(i,tt,gen)) +
	- 2muu3sqrtS(-4p1k1*2 - 4p1k2*2 + 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(1 + (2 - 12muu12)muu3*2 - 2x1 + 3muu12x1 - 2x2 + 3muu1*2x2) +
	- 2p1k1(-4p1k2 + sqrtS2(1 + x1 + x2)))*HR(i,tt,gen))) +
	- BL(muu3sqrtSGR(i,tt,gen)(2muu3sqrtS(4p1k1*2 + 4p1k2*2 - 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + 2x1 - 3muu12x1 + 2x2 - 3muu1*2x2) +
	- p1k1(8p1k2 - 2sqrtS2(1 + x1 + x2)))*
	- HL(i,tt,gen) + muu1sqrtS3(6p1k2x1 + 2p1k1(-4 + x2) + 4muu32(2p1k1 - 4p1k2 + sqrtS*2x2) +
	- sqrtS*2(1 + x1 - x2 - 2x1x2))*
	- HR(i,tt,gen)) + GL(i,tt,gen)(muu1muu3sqrtS4
	- (-2p1k2(-4 + 8muu32 + x1) + p1k1(4 + 8muu32 - 4x1 - 6*x2) +
	- sqrtS*2(-1 + x1 - 3x2 + 2x1x2 + muu32(4 - 4x1 + 4x2)))*
	- HL(i,tt,gen) + sqrtS*2(-4p1k12x2 - (2p1k2 - sqrtS2x2)(2p1k2x1 - sqrtS2(-1 + x1 + x2)) +
	- 2muu12sqrtS*4(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)) +
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) -
	- 2p1k1(2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + 3x2 - 2x1x2 - x22)))HR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x1)(-1 + x2)(2p1k2 + sqrtS*2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))


	** USchannel bremsstrahlung 21

	radUS = radUS + NCCF(4ALgs*2(BR(GR(i,tt,gen)(sqrtS*2(4p1k12x2 + (2p1k2 - sqrtS2x2)*
	- (2p1k2x1 - sqrtS*2(-1 + x1 + x2)) +
	- 2muu12sqrtS*4(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) -
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) +
	- 2p1k1(2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + 3x2 - 2x1x2 - x22)))HL(i,tt,gen) +
	- muu1muu3sqrtS*4(2p1k2(-4 + x1) + 6p1k1x2 + sqrtS*2(-1 + x1 + 3x2 - 2x1*x2) -
	- 4muu32(4p1k1 - 2p1k2 + sqrtS*2(-x1 + x2)))*
	- HR(i,tt,gen)) + muu3sqrtSGL(i,tt,gen)(muu1sqrtS*3
	- (p1k2(4 + 8muu3*2 - 6x1 - 4x2) - 2p1k1(-4 + 8muu3**2 + x2) +
	- sqrtS*2(4muu32(-1 + x1) + (-1 + x1 + x2)(-1 + 2x2)))*HL(i,tt,gen) -
	- 2muu3sqrtS(4p1k1*2 + 4p1k2*2 - 2p1k2sqrtS2(-1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + x1 - 3muu1*2x1 + x2 - 3muu12x2) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)))*HR(i,tt,gen))) +
	- BL(muu3sqrtSGR(i,tt,gen)(2muu3sqrtS(4p1k1*2 + 4p1k2*2 - 2p1k2sqrtS2(-1 + x1 + x2) +
	- sqrtS*4(-1 + 2(-1 + 6muu1*2)muu3*2 + x1 - 3muu1*2x1 + x2 - 3muu12x2) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)))*HL(i,tt,gen) +
	- muu1sqrtS3(2p1k1(-4 + 8muu32 + x2) + p1k2(-4 - 8muu32 + 6x1 + 4*x2) +
	- sqrtS*2(-4muu32(-1 + x1) - (-1 + x1 + x2)(-1 + 2x2)))*HR(i,tt,gen)) +
	- GL(i,tt,gen)(muu1muu3sqrtS4(sqrtS*2 - 2p1k2(-4 + x1) - sqrtS2x1 - 6p1k1x2 -
	- 3sqrtS2x2 + 2sqrtS2x1*x2 +
	- 4muu32(4p1k1 - 2p1k2 + sqrtS*2(-x1 + x2)))*HL(i,tt,gen) +
	- sqrtS*2(-4p1k12x2 - (2p1k2 - sqrtS2x2)(2p1k2x1 - sqrtS2(-1 + x1 + x2)) +
	- 2muu12sqrtS*4(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)) +
	- 2muu32(4p1k12 - 2p1k1sqrtS2(-1 + x1) + (2p1k2 - sqrtS2)(2p1k2 - sqrtS2x2)) -
	- 2p1k1(2p1k2(-2 + x1 + x2) + sqrtS*2(-1 + x1 + 3x2 - 2x1x2 - x22)))HR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-MSf(i,tt,gen)2 - 2p1k1 + (muu12 + muu32)sqrtS2)(-1 + x1)*(-1 + x2))


	** USchannel bremsstrahlung 22

	radUS = radUS + NCCF(8ALgs*2(BR(muu3sqrtS*3(1 + 2muu32 - x1)GL(i,tt,gen)*
	- (muu1(-6p1k1sqrtS + sqrtS3x1)HL(i,tt,gen) + (1 - 6muu1*2)muu3sqrtS3HR(i,tt,gen)) +
	- GR(i,tt,gen)((muu12sqrtS*6(-1 - 4muu34 + x1 + 2muu3*2x1 + x2 - x1*x2) +
	- 2p1k1sqrtS*2(-2p1k1(-1 + 2muu32 + x1) + (-1 + x1)(-2*p1k2 +
	- sqrtS*2(-1 + 2muu32 + x1 + x2))))HL(i,tt,gen) +
	- muu1muu3sqrtS*4(-2muu32(6p1k1 + sqrtS2(1 - 2*x1)) -
	- (-1 + x1)(6p1k1 + 6p1k2 + sqrtS2(1 - 2x1 - 2x2)))*HR(i,tt,gen))) +
	- BL(-(muu3sqrtS*3(1 + 2muu32 - x1)GR(i,tt,gen)*
	- ((1 - 6muu12)muu3sqrtS3HL(i,tt,gen) + muu1sqrtS(-6p1k1 + sqrtS2x1)*HR(i,tt,gen))) +
	- GL(i,tt,gen)(muu1muu3sqrtS4(2muu32(6p1k1 + sqrtS2(1 - 2*x1)) +
	- (-1 + x1)(6p1k1 + 6p1k2 + sqrtS2(1 - 2x1 - 2x2)))*HL(i,tt,gen) +
	- (muu1*2sqrtS*6(4muu34 - 2muu3*2x1 + (-1 + x1)*(-1 + x2)) +
	- 2p1k1sqrtS*2(2p1k1(-1 + 2muu3*2 + x1) +
	- (-1 + x1)(2p1k2 - sqrtS*2(-1 + 2muu32 + x1 + x2))))HR(i,tt,gen)))))/
	- (sqrtS*4(-DenSmZ1)(-1 + x1)2(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2))


	** USchannel bremsstrahlung 31

	radUS = radUS + (4ALCFgs2NC(BR(GR(i,tt,gen)((2(2p1k2 - sqrtS2x2)(-(muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2)) + p1k1(-4p1k2 + sqrtS2(1 + x2))) +
	- muu1*2sqrtS*4(8muu34sqrtS*2 - 2p1k2(2 + x1 - 2x2) + (2p1k1 + sqrtS2(1 + x1 - x2))*x2 -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2))))HL(i,tt,gen) -
	- muu1muu3sqrtS*2(8p1k22 - 2p1k1(8p1k2 + sqrtS*2(2muu32 - 3x2)) +
	- 2p1k2sqrtS*2(3 + 10muu32 + x1 - 7x2) +
	- sqrtS*4(1 + 2muu12(-1 + x1) - x1 + muu3*2(2 - 8x2) - 3x2 + 4x22))HR(i,tt,gen)) +
	- muu3sqrtSGL(i,tt,gen)(muu1sqrtS(2muu3*2sqrtS*2(6p1k1 - 6p1k2 + sqrtS*2(1 - x1 + x2)) +
	- p1k1(24p1k2 - 2sqrtS2(1 + 4*x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 6muu12(-1 + x2) + x2 + x1x2)))HL(i,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + 2p1k1x2 + sqrtS2(-2 + x1 + 2*x2) -
	- 2muu12(6p1k1 - 6p1k2 + sqrtS*2(-5 - 12muu32 + x1 + 7x2)))*HR(i,tt,gen))) +
	- BL(-(muu3sqrtSGR(i,tt,gen)(muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + 2p1k1*x2 +
	- sqrtS*2(-2 + x1 + 2*x2) -
	- 2muu12(6p1k1 - 6p1k2 + sqrtS*2(-5 - 12muu32 + x1 + 7x2)))*HL(i,tt,gen) +
	- muu1sqrtS(2muu32sqrtS*2(6p1k1 - 6p1k2 + sqrtS*2(1 - x1 + x2)) + p1k1(24p1k2 -
	- 2sqrtS2(1 + 4*x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 6muu12(-1 + x2) + x2 + x1x2)))HR(i,tt,gen))) +
	- GL(i,tt,gen)(muu1muu3sqrtS2(8p1k22 - 2p1k1(8p1k2 + sqrtS*2(2muu32 - 3x2)) +
	- 2p1k2sqrtS*2(3 + 10muu32 + x1 - 7x2) +
	- sqrtS*4(1 + 2muu12(-1 + x1) - x1 + muu3*2(2 - 8x2) - 3x2 + 4x22))HL(i,tt,gen) +
	- (-2(2p1k2 - sqrtS*2x2)(-(muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2)) +
	- p1k1(-4p1k2 + sqrtS*2(1 + x2))) -
	- muu1*2sqrtS*4(8muu34sqrtS*2 - 2p1k2(2 + x1 - 2x2) + (2p1k1 + sqrtS2(1 + x1 - x2))*x2 -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2))))HR(i,tt,gen)))))/
	- (sqrtS*2(-DenSmZ1)(-1 + x2)(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))


	** USchannel bremsstrahlung 32

	radUS = radUS + (4ALCFgs2NC(BR(muu3sqrtSGL(i,tt,gen)(-(muu1sqrtS*
	- (8p1k12 + 2muu3*2sqrtS*2(-10p1k1 + 2p1k2 + sqrtS*2(1 + x1 - x2)) +
	- p1k1(-16p1k2 + 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HL(i,tt,gen)) +
	- muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + 2muu1*2(-6p1k1 + 6p1k2 +
	- sqrtS*2(5 + 12muu32 - 4x1 - 4x2)) + 2p1k1*x2 +
	- sqrtS*2(-2 + 2x1 + x2))HR(i,tt,gen)) + GR(i,tt,gen)*
	- ((2p1k1((4p1k1 - sqrtS2(-1 + x1))(2p1k2 - sqrtS*2x2) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2))) +
	- muu1*2sqrtS*4(8muu34sqrtS**2 +
	- 2muu32(-2p1k1 + 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) + p1k1(4 - 4x1 + 2*x2) +
	- x1(-2p1k2 + sqrtS*2(-1 + x1 + x2))))*HL(i,tt,gen) +
	- muu1muu3sqrtS*2(2p1k1(12p1k2 - 5sqrtS*2x2) +
	- 2muu32sqrtS*2(6p1k1 - 6p1k2 + sqrtS*2(1 - 2x1 + 2x2)) +
	- sqrtS*2(p1k2(2 - 4x1) + sqrtS*2(-1 - 6muu12
	- (-1 + x1) + x1 - x2 + 2x1x2)))*HR(i,tt,gen))) -
	- BL(muu3sqrtSGR(i,tt,gen)(muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2*x1 +
	- 2muu12(-6p1k1 + 6p1k2 + sqrtS*2(5 + 12muu32 - 4x1 - 4*x2)) +
	- 2p1k1x2 + sqrtS*2(-2 + 2x1 + x2))HL(i,tt,gen) -
	- muu1sqrtS(8p1k12 + 2muu3*2sqrtS*2(-10p1k1 + 2p1k2 +
	- sqrtS*2(1 + x1 - x2)) + p1k1(-16p1k2 + 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HR(i,tt,gen)) +
	- GL(i,tt,gen)(muu1muu3sqrtS2(2p1k1(12p1k2 - 5sqrtS*2x2) +
	- 2muu32sqrtS*2(6p1k1 - 6p1k2 + sqrtS*2(1 - 2x1 + 2x2)) +
	- sqrtS*2(p1k2(2 - 4x1) + sqrtS*2(-1 - 6muu12(-1 + x1) + x1 - x2 + 2x1x2)))*HL(i,tt,gen) +
	- (2p1k1((4p1k1 - sqrtS2(-1 + x1))(2p1k2 - sqrtS*2x2) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2))) +
	- muu1*2sqrtS*4(8muu34sqrtS*2 + 2muu3*2(-2p1k1 + 2p1k2 +
	- sqrtS*2(1 - 2x1 - 2x2)) + p1k1(4 - 4x1 + 2*x2) +
	- x1(-2p1k2 + sqrtS*2(-1 + x1 + x2))))*HR(i,tt,gen)))))/
	- (sqrtS*2(-DenSmZ1)(-1 + x1)(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))


	** USchannel dipole Z 312

	radUS = radUS - (8ALCFGS2NC((-1 + x2)(2 + 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) +
	- x1*2dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- 3x1(-1 + x2)dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 5x2dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 2x22dsqrt((1 - 4muu32)/(-4muu32 + x22))) +
	- 4muu34(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) + x1dsqrt((1 - 4muu3**2)/
	- (-4muu32 + x2*2)) +
	- x2(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x22)))) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 2x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 3x22dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- x2(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x22)) + 3x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x22)))))
	- (8(-1 + 4muu3*2)p1k2*2(BRGR(i,tt,gen)HL(i,tt,gen) - BLGL(i,tt,gen)HR(i,tt,gen)) +
	- 4p1k2sqrtS*2(BR(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)GL(i,tt,gen)*HL(i,tt,gen) +
	- GR(i,tt,gen)((x2 - 4muu3*2x2 + dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22))HL(i,tt,gen) +
	- 3muu1muu3dsqrt(1 - 4muu3*2)dsqrt(-4muu32 + x22)HR(i,tt,gen))) -
	- BL(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)GR(i,tt,gen)*HR(i,tt,gen) +
	- GL(i,tt,gen)(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22)HL(i,tt,gen) +
	- (x2 - 4muu32x2 + dsqrt(1 - 4muu32)dsqrt(-4muu32 + x22))HR(i,tt,gen)))) +
	- sqrtS*4(BL(GL(i,tt,gen)(muu1muu3(-4muu32 + x2*2 +
	- 3dsqrt(1 - 4muu3*2)x2dsqrt(-4muu32 + x22))*HL(i,tt,gen) +
	- (-2muu32(1 + muu1*2(-4 + 8muu32)) + (-1 + 2muu1*2)(-1 + 2muu32)x2**2 +
	- dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*
	- HR(i,tt,gen)) + muu3GR(i,tt,gen)(2(-1 + 6muu1*2)muu3(-4muu32 + x22)*HL(i,tt,gen) +
	- muu1(-4muu32 + x22 + 3dsqrt(1 - 4muu3*2)x2dsqrt(-4muu32 + x22))*HR(i,tt,gen))) -
	- BR(muu3GL(i,tt,gen)(muu1(-4muu32 + x22 + 3dsqrt(1 - 4muu32)x2*
	- dsqrt(-4muu32 + x22))HL(i,tt,gen) +
	- 2(-1 + 6muu1*2)muu3(-4muu32 + x22)*HR(i,tt,gen)) +
	- GR(i,tt,gen)((-2muu3*2(1 + muu1*2(-4 + 8muu32)) + (-1 + 2muu1*2)(-1 + 2muu32)x2**2 +
	- dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*HL(i,tt,gen) +
	- muu1muu3(-4muu32 + x22 + 3dsqrt(1 - 4muu32)x2dsqrt(-4muu32 + x22))*HR(i,tt,gen))))))/
	- ((-1 + 2muu32)sqrtS*2(-DenSmZ1)(-1 + x2)2(-2 + x1 + x2)dsqrt(-4muu32 + x22)*
	- (4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) +
	- (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x2*2)) -
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)*2))


	** USchannel dipole Z 321

	radUS = radUS - (8ALCFGS2NC(4muu3*4(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) +
	- x1(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x12))) +
	- dsqrt((1 - 4muu32)/(-4muu32 + x12))x2) + (-1 + x1)
	- (2 + 2dsqrt((-1 + 4muu3*2)/(4muu32 - x12)) + 2x12dsqrt((1 - 4muu32)/(-4muu32 + x12)) -
	- 3dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2 + dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2*2 +
	- x1dsqrt((1 - 4muu3*2)/(-4muu32 + x12))(-5 + 3x2)) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3x12dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x1*2)) -
	- 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12))x2 + x1
	- (4 - 8dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3dsqrt((1 - 4muu3*2)/(-4muu32 + x12))x2)))
	- (8(-1 + 4muu3*2)p1k1*2(BRGR(i,tt,gen)HL(i,tt,gen) - BLGL(i,tt,gen)HR(i,tt,gen)) -
	- 4p1k1sqrtS*2(-(BL(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)GR(i,tt,gen)*HR(i,tt,gen) +
	- GL(i,tt,gen)(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)HL(i,tt,gen) -
	- x1HR(i,tt,gen) + 4muu3*2x1*HR(i,tt,gen) +
	- dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)HR(i,tt,gen)))) +
	- BR(3muu1muu3dsqrt(1 - 4muu32)dsqrt(-4muu32 + x12)GL(i,tt,gen)*HL(i,tt,gen) +
	- GR(i,tt,gen)(-(x1HL(i,tt,gen)) + 4muu32x1HL(i,tt,gen) + dsqrt(1 - 4muu3*2)
	- dsqrt(-4muu32 + x12)HL(i,tt,gen) +
	- 3muu1muu3dsqrt(1 - 4muu3*2)dsqrt(-4muu32 + x12)HR(i,tt,gen)))) +
	- sqrtS*4(-(BL(GL(i,tt,gen)(muu1muu3(4muu32 - x12 + 3dsqrt(1 - 4muu32)x1*
	- dsqrt(-4muu32 + x12))HL(i,tt,gen) +
	- (2muu32(1 + muu1*2(-4 + 8muu32)) - (-1 + 2muu1*2)(-1 + 2muu32)x1**2 +
	- dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*
	- HR(i,tt,gen)) + muu3GR(i,tt,gen)(-2(-1 + 6muu1*2)muu3(-4muu32 + x12)*HL(i,tt,gen) +
	- muu1(4muu32 - x12 + 3dsqrt(1 - 4muu3*2)x1dsqrt(-4muu32 + x12))*HR(i,tt,gen)))) +
	- BR(muu3GL(i,tt,gen)(muu1(4muu32 - x1*2 +
	- 3dsqrt(1 - 4muu3*2)x1dsqrt(-4muu32 + x12))*HL(i,tt,gen) -
	- 2(-1 + 6muu1*2)muu3(-4muu32 + x12)*HR(i,tt,gen)) +
	- GR(i,tt,gen)((2muu3*2(1 + muu1*2(-4 + 8muu32)) - (-1 + 2muu1*2)(-1 + 2muu32)x1**2 +
	- dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*HL(i,tt,gen) +
	- muu1muu3(4muu32 - x12 + 3dsqrt(1 - 4muu32)x1dsqrt(-4muu32 + x12))*HR(i,tt,gen))))))/
	- ((-1 + 2muu32)sqrtS*2(-DenSmZ1)(-1 + x1)2dsqrt(-4muu32 + x12)(-2 + x1 + x2)*
	- (-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 +
	- (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x1*2)) -
	- 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)*2))


	do j = 1, 4

	** USchannel bremsstrahlung 13

	radUS = radUS + (2CFgs*2NC(CL(j)(DL(j)(2muu1muu3sqrtS*2GR(i,tt,gen)*
	- (sqrtS*2(-2p1k1(-1 + x2) - (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))HL(i,tt,gen) +
	- 2muu1muu3sqrtS4(1 + 2muu32 - x2)HR(i,tt,gen))
	- + sqrtS*4GL(i,tt,gen)((-4muu3*4sqrtS*2 + 2muu1*2sqrtS*2(4muu3*4 +
	- (-1 + x1)(-1 + x2) - 2muu3*2x2) +
	- 4muu32(-2p1k2 + sqrtS2x2) + (-1 + x2)(2p1k2(-2 + x1) - 2p1k1*x2 +
	- sqrtS*2(1 - x1 + x2)))*HL(i,tt,gen) +
	- 2muu1muu3(1 + 2muu3*2 - x2)(-2p1k2 + sqrtS2x2)*HR(i,tt,gen))) +
	- 2DR(j)(muu1sqrtSGR(i,tt,gen)(-2muu3p1k2sqrtS*3(1 + 2muu32 - x2)HL(i,tt,gen) +
	- muu1sqrtS5(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2)*HR(i,tt,gen)) +
	- muu3sqrtSGL(i,tt,gen)((-1 + 2muu1*2)muu3sqrtS5(1 + 2muu32 - x2)HL(i,tt,gen) +
	- muu1sqrtS3(-2p1k1(-1 + x2) - 2p1k2(-1 + 2muu3*2 + x2) +
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))HR(i,tt,gen)))) +
	- CR(j)(2DL(j)(muu1sqrtSGL(i,tt,gen)(muu1sqrtS5(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2)*HL(i,tt,gen) -
	- 2muu3p1k2sqrtS3(1 + 2muu32 - x2)HR(i,tt,gen)) +
	- muu3sqrtSGR(i,tt,gen)(muu1sqrtS*3(-2p1k1(-1 + x2) - 2p1k2(-1 + 2muu3*2 + x2) +
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))
	- HL(i,tt,gen) + (-1 + 2muu12)muu3sqrtS5(1 + 2muu32 - x2)HR(i,tt,gen))) +
	- DR(j)(2muu1muu3sqrtS*2GL(i,tt,gen)(2muu1muu3sqrtS*4(1 + 2muu32 - x2)HL(i,tt,gen) +
	- sqrtS*2(-2p1k1(-1 + x2) - (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))HR(i,tt,gen)) +
	- sqrtS*4GR(i,tt,gen)(2muu1muu3(1 + 2muu32 - x2)(-2p1k2 + sqrtS2x2)*HL(i,tt,gen) +
	- (-4muu34sqrtS*2 + 2muu1*2sqrtS*2(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2) +
	- 4muu32(-2p1k2 + sqrtS2x2) +
	- (-1 + x2)(2p1k2(-2 + x1) - 2p1k1x2 + sqrtS2(1 - x1 + x2)))*HR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x2)*2(-DenSmH1(j))(2p1k2 + sqrtS*2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))



	** USchannel bremsstrahlung 14

	radUS = radUS + (-2CFgs*2NC(CL(j)(DL(j)(muu1muu3sqrtS2GR(i,tt,gen)*
	- (sqrtS*2(sqrtS*2(1 - 4muu32 + x1 - x2) + 2p1k1*x2 +
	- p1k2(-4 + 8muu3*2 - 2x1 + 4x2))HL(i,tt,gen) +
	- 2muu1muu3sqrtS4(-4muu32 + x1 + x2)HR(i,tt,gen)) +
	- sqrtS*4GL(i,tt,gen)(-2(-2muu34sqrtS*2 + x1(2p1k2 - sqrtS2x2) +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)) +
	- muu3*2(-2p1k1 - 6p1k2 + sqrtS*2(1 + x1 + 2x2)))HL(i,tt,gen) +
	- muu1muu3(-2p1k2x1 + 2p1k1x2 + sqrtS*2(-1 + x1 + x2) + muu3*2
	- (8p1k2 - 4sqrtS*2x2))*HR(i,tt,gen))) +
	- DR(j)(muu1sqrtSGR(i,tt,gen)(muu3sqrtS3(8muu32p1k2 - 2p1k2x1 + 2p1k1x2 -
	- sqrtS*2(-1 + x1 + x2))*HL(i,tt,gen) +
	- 2muu1sqrtS*5(-4muu34 - x1x2 + muu3*2(2 + x1 + x2))*HR(i,tt,gen)) +
	- muu3sqrtS4GL(i,tt,gen)(2muu3(2p1k1 + 2p1k2 - sqrtS2((-2 + 4muu12)muu3**2 -
	- (-1 + muu1*2)(x1 + x2)))*HL(i,tt,gen) +
	- muu1(2p1k1x2 + p1k2(-4 + 8muu32 - 2x1 + 4x2) - sqrtS2(1 + x1 +
	- 4muu32(-1 + x2) - 3x2 + 2x2*2))HR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu1sqrtSGL(i,tt,gen)(2muu1sqrtS5(-4muu34 - x1x2 +
	- muu3*2(2 + x1 + x2))*HL(i,tt,gen) +
	- muu3sqrtS3(8muu32p1k2 - 2p1k2x1 + 2p1k1x2 - sqrtS*2(-1 + x1 + x2))*HR(i,tt,gen)) +
	- muu3sqrtS4GR(i,tt,gen)(muu1(2p1k1x2 + p1k2(-4 + 8muu3*2 - 2x1 + 4*x2) -
	- sqrtS*2(1 + x1 + 4muu32(-1 + x2) - 3x2 + 2x2*2))
	- HL(i,tt,gen) + 2muu3(2p1k1 + 2p1k2 - sqrtS*2((-2 + 4muu12)muu3**2 -
	- (-1 + muu1*2)(x1 + x2)))*HR(i,tt,gen))) +
	- DR(j)(muu1muu3sqrtS2GL(i,tt,gen)(2muu1muu3sqrtS*4(-4muu32 + x1 + x2)HL(i,tt,gen) +
	- sqrtS*2(sqrtS*2(1 - 4muu32 + x1 - x2) + 2p1k1*x2 +
	- p1k2(-4 + 8muu3*2 - 2x1 + 4x2))HR(i,tt,gen)) +
	- sqrtS*4GR(i,tt,gen)(muu1muu3(-2p1k2x1 + 2p1k1x2 + sqrtS2(-1 + x1 + x2) +
	- muu3*2(8p1k2 - 4sqrtS*2x2))*HL(i,tt,gen) -
	- 2(-2muu3*4sqrtS*2 + x1(2p1k2 - sqrtS2x2) + muu1*2sqrtS*2(4muu3*4 +
	- x1x2 - muu32(2 + x1 + x2)) +
	- muu3*2(-2p1k1 - 6p1k2 + sqrtS*2(1 + x1 + 2x2)))HR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x1)(-1 + x2)(DenSmH1(j))(-2p1k2 - sqrtS*2(muu12 + muu32 - x2) + MSf(i,tt,gen)**2))


	** USchannel bremsstrahlung 23

	radUS = radUS + (2CFgs*2NC(CL(j)(DR(j)(muu3sqrtS*4GL(i,tt,gen)(2muu3(2p1k1 + 2*p1k2 +
	- sqrtS*2((-2 + 4muu12)muu32 - muu12(x1 + x2)))
	- HL(i,tt,gen) + muu1(2p1k2x1 + muu32(8p1k1 - 4sqrtS*2x1) - 2p1k1x2 +
	- sqrtS*2(-1 + x1 + x2))*HR(i,tt,gen)) +
	- muu1sqrtSGR(i,tt,gen)(muu3sqrtS*3(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2*x2) +
	- sqrtS*2(1 - 4muu32 - x1 + x2))HL(i,tt,gen) +
	- 2muu1sqrtS*5(4muu34 + x1x2 - muu3*2(2 + x1 + x2))*HR(i,tt,gen))) +
	- DL(j)(muu1muu3sqrtS2GR(i,tt,gen)(sqrtS2(8muu32p1k1 + 2p1k2x1 - 2p1k1x2 -
	- sqrtS*2(-1 + x1 + x2))*HL(i,tt,gen) +
	- 2muu1muu3sqrtS4(4muu32 - x1 - x2)HR(i,tt,gen)) +
	- sqrtS*4GL(i,tt,gen)(2(-2muu34sqrtS2 + muu32(6p1k1 + 2p1k2 - sqrtS2(-1 + x1)) -
	- 2p1k1x2 +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)))*HL(i,tt,gen) +
	- muu1muu3(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2x2) - sqrtS2(1 + 4muu32(-1 + x1) -
	- 3x1 + 2x1*2 + x2))HR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu1sqrtSGL(i,tt,gen)(2muu1sqrtS5(4muu34 + x1x2 - muu3*2(2 + x1 + x2))*HL(i,tt,gen) +
	- muu3sqrtS3(2p1k2x1 + p1k1(-4 + 8muu3*2 + 4x1 - 2*x2) +
	- sqrtS*2(1 - 4muu32 - x1 + x2))HR(i,tt,gen)) +
	- muu3sqrtS4GR(i,tt,gen)(muu1(2p1k2x1 + muu3*2(8p1k1 - 4sqrtS*2x1) -
	- 2p1k1x2 + sqrtS*2(-1 + x1 + x2))*HL(i,tt,gen) +
	- 2muu3(2p1k1 + 2p1k2 + sqrtS*2((-2 + 4muu12)muu32 - muu12(x1 + x2)))HR(i,tt,gen))) +
	- DR(j)(muu1muu3sqrtS2GL(i,tt,gen)(2muu1muu3sqrtS*4(4muu32 - x1 - x2)HL(i,tt,gen) +
	- sqrtS*2(8muu32p1k1 + 2p1k2x1 - 2p1k1x2 - sqrtS*2(-1 + x1 + x2))*HR(i,tt,gen)) +
	- sqrtS*4GR(i,tt,gen)(muu1muu3(2p1k2x1 + p1k1(-4 + 8muu32 + 4x1 - 2*x2) -
	- sqrtS*2(1 + 4muu32(-1 + x1) - 3x1 + 2x1*2 + x2))
	- HL(i,tt,gen) + 2(-2muu3*4sqrtS2 + muu32(6p1k1 + 2p1k2 - sqrtS2(-1 + x1)) - 2p1k1x2 +
	- muu1*2sqrtS*2(4muu34 + x1x2 - muu3*2(2 + x1 + x2)))*HR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x1)(-1 + x2)(-DenSmH1(j))(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)*2))

	** USchannel bremsstrahlung 24

	radUS = radUS + (2CFgs*2NC(CR(j)(2DL(j)(-(muu3sqrtS3(1 + 2muu32 - x1)GR(i,tt,gen)*
	- (muu1(-2p1k1sqrtS + sqrtS3x1)HL(i,tt,gen) + (1 - 2muu1*2)muu3sqrtS3HR(i,tt,gen))) +
	- muu1sqrtSGL(i,tt,gen)(muu1sqrtS*5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))HL(i,tt,gen) +
	- muu3sqrtS3(muu3*2(4p1k1 - 2sqrtS*2) + (2p1k1 + 2p1k2 - sqrtS2)(-1 + x1))*HR(i,tt,gen))) +
	- sqrtS*4DR(j)(4muu1muu3(1 + 2muu32 - x1)GL(i,tt,gen)*
	- (muu1muu3sqrtS*2HL(i,tt,gen) + p1k1*HR(i,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3(2p1k1(-1 + 2muu32 + x1) + (-1 + x1)(2*p1k2 -
	- sqrtS*2(-1 + 2muu32 + x1 + x2)))HL(i,tt,gen) +
	- (8muu32p1k1 - 4muu34sqrtS*2 + 2muu1*2sqrtS*2
	- (4muu34 - 2muu3*2x1 + (-1 + x1)*(-1 + x2)) -
	- (-1 + x1)(-2p1k2x1 + 2p1k1(-2 + x2) + sqrtS2(-1 + x1 + x2)))*HR(i,tt,gen)))) +
	- CL(j)(2DR(j)(-(muu3sqrtS*3(1 + 2muu32 - x1)GL(i,tt,gen)*
	- ((1 - 2muu12)muu3sqrtS3HL(i,tt,gen) + muu1sqrtS(-2p1k1 + sqrtS2x1)*HR(i,tt,gen))) +
	- muu1sqrtSGR(i,tt,gen)(muu3sqrtS*3(muu3*2(4p1k1 - 2sqrtS**2) +
	- (2p1k1 + 2p1k2 - sqrtS*2)(-1 + x1))*HL(i,tt,gen) +
	- muu1sqrtS5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))HR(i,tt,gen))) +
	- sqrtS*4DL(j)(4muu1muu3(1 + 2muu32 - x1)GR(i,tt,gen)*
	- (p1k1HL(i,tt,gen) + muu1muu3sqrtS2HR(i,tt,gen)) +
	- GL(i,tt,gen)((8muu3*2p1k1 - 4muu34sqrtS*2 + 2muu1*2sqrtS*2
	- (4muu34 - 2muu3*2x1 + (-1 + x1)*(-1 + x2)) -
	- (-1 + x1)(-2p1k2x1 + 2p1k1(-2 + x2) + sqrtS2(-1 + x1 + x2)))*HL(i,tt,gen) +
	- 2muu1muu3(2p1k1(-1 + 2muu3*2 + x1) + (-1 + x1)(2p1k2 - sqrtS2
	- (-1 + 2muu32 + x1 + x2)))HR(i,tt,gen))))))/
	- (sqrtS*4(-1 + x1)*2(-DenSmH1(j))(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)*2))

	** USchannel bremsstrahlung 33

	radUS = radUS + (CFgs2NC(CL(j)(sqrtS*2DL(j)(4muu1muu3GR(i,tt,gen)*
	- ((-(muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS2)) + muu12sqrtS4(-1 + x2) + p1k1*
	- (-4p1k2 + sqrtS2(1 + x2)))*HL(i,tt,gen) -
	- muu1muu3sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - 2x2))*HR(i,tt,gen)) +
	- GL(i,tt,gen)((8muu3*4sqrtS*4 + p1k1(-16p1k2 + 2sqrtS*2x2*(3 + x2)) +
	- sqrtS*2(-1 + x2)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32sqrtS*2(-2p1k2(2 + x1) + 2p1k1(2 + x2) + sqrtS*2(-2 + x1 + 4*x2)) -
	- 2muu12sqrtS*2(8muu34sqrtS*2 - 2p1k2x1 + 2p1k1*x2 +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3*x2)) +
	- muu3*2(-4p1k1 + 4p1k2 - 2sqrtS2(-3 + 2x1 + 4x2))))*HL(i,tt,gen) +
	- 2muu1muu3(p1k1(-8p1k2 + 4sqrtS*2x2) - 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(2muu12sqrtS*2(-1 + x1) + (-1 + 2x1 - x2)(2p1k2 - sqrtS2x2)))*HR(i,tt,gen))) -
	- 2DR(j)(muu1sqrtSGR(i,tt,gen)(muu3sqrtS(-2muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2) +
	- p1k1(8p1k2 - 2sqrtS2x2) +
	- sqrtS*2(-2p1k2(1 + x1 - x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HL(i,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 - 2p1k2x1 + x2(2p1k1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2)))HR(i,tt,gen)) +
	- muu3sqrtSGL(i,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 - 2p1k2x1 + muu12
	- (-4p1k1 + 4p1k2 + 2sqrtS2(3 + 4muu32 - x1 - 3x2)) +
	- 2p1k1x2 + sqrtS*2(-2 + x1 + 2x2))HL(i,tt,gen) +
	- muu1sqrtS(2p1k1(4p1k2 + sqrtS2(1 - 2x2)) + 2muu3*2sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 2muu12(-1 + x2) + x2 + x1x2)))HR(i,tt,gen)))) +
	- CR(j)(-2DL(j)(muu3sqrtSGR(i,tt,gen)(muu1sqrtS(2p1k1(4p1k2 + sqrtS2(1 - 2*x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 2muu12(-1 + x2) + x2 + x1x2)))
	- HL(i,tt,gen) + muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2*x1 +
	- muu1*2(-4p1k1 + 4p1k2 + 2sqrtS2(3 + 4muu32 - x1 - 3x2)) +
	- 2p1k1x2 + sqrtS*2(-2 + x1 + 2x2))HR(i,tt,gen)) +
	- muu1sqrtSGL(i,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 - 2p1k2x1 + x2(2p1k1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2)))HL(i,tt,gen) +
	- muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2(1 + x1 - x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HR(i,tt,gen))) +
	- sqrtS*2DR(j)(4muu1muu3GL(i,tt,gen)(-(muu1muu3sqrtS2
	- (-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - 2x2))*HL(i,tt,gen)) +
	- (-(muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS2)) + muu12sqrtS4(-1 + x2) +
	- p1k1(-4p1k2 + sqrtS*2(1 + x2)))*HR(i,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3(p1k1(-8p1k2 + 4sqrtS2x2) - 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(2muu12sqrtS*2(-1 + x1) + (-1 + 2x1 - x2)(2p1k2 - sqrtS2x2)))*HL(i,tt,gen) +
	- (8muu34sqrtS*4 + p1k1(-16p1k2 + 2sqrtS*2x2(3 + x2)) + sqrtS2
	- (-1 + x2)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32sqrtS*2(-2p1k2(2 + x1) + 2p1k1(2 + x2) + sqrtS*2(-2 + x1 + 4*x2)) -
	- 2muu12sqrtS*2(8muu34sqrtS*2 - 2p1k2x1 + 2p1k1*x2 +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3*x2)) +
	- muu3*2(-4p1k1 + 4p1k2 - 2sqrtS2(-3 + 2x1 + 4x2))))*HR(i,tt,gen))))))/
	- (sqrtS*2(-1 + x2)(DenSmH1(j))(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (-2p1k2 - sqrtS2(muu12 + muu32 - x2) + MSf(i,tt,gen)**2))


	** USchannel bremsstrahlung 34

	radUS = radUS + (CFgs2NC(CL(j)(sqrtS*2DL(j)(4muu1muu3GR(i,tt,gen)*
	- ((-(muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2)) + p1k1(-4p1k2 + sqrtS2(1 + x1)) +
	- muu1*2sqrtS*4(-1 + x2))*HL(i,tt,gen) +
	- muu1muu3sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - 4muu32 + x1 + x2))HR(i,tt,gen)) +
	- GL(i,tt,gen)((8muu3*4sqrtS*4 + p1k1(-16p1k2 + 2sqrtS*2(3 + x1)*x2) +
	- sqrtS*2(-1 + x1)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) +
	- muu3*2sqrtS*2(4p1k2(2 + x1) - 4p1k1(2 + x2) - 2sqrtS2(-2 + 2x1 + 3x2)) -
	- 2muu12sqrtS*2(8muu34sqrtS*2 - 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) + 2p1k1x2 +
	- muu3*2(-4p1k1 + 4p1k2 - 6sqrtS2(-1 + x1 + x2))))*HL(i,tt,gen) +
	- 2muu1muu3(sqrtS2(-1 + x1)(2p1k2 + sqrtS*2(2muu12 - x2)) + p1k1(-8p1k2 + 4sqrtS*2x2) -
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HR(i,tt,gen))) -
	- 2DR(j)(muu1sqrtSGR(i,tt,gen)(muu3sqrtS(-2muu3*2sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2 + sqrtS2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HL(i,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 - 2p1k2x1 + 2muu3*2(-2p1k1 + 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) + 2p1k1x2 +
	- sqrtS*2x1(-1 + x1 + x2))HR(i,tt,gen)) +
	- muu3sqrtSGL(i,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 - 2p1k2*x1 +
	- muu1*2(-4p1k1 + 4p1k2 + 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) + 2p1k1x2 + sqrtS2(-2 + 2x1 + x2))HL(i,tt,gen) +
	- muu1sqrtS(2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) + 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HR(i,tt,gen)))) +
	- CR(j)(-2DL(j)(muu3sqrtSGR(i,tt,gen)(muu1sqrtS(2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)) + sqrtS*2(-2p1k2x1 +
	- sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HL(i,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + muu12(-4p1k1 + 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) + 2p1k1*x2 +
	- sqrtS*2(-2 + 2x1 + x2))HR(i,tt,gen)) +
	- muu1sqrtSGL(i,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 - 2p1k2*x1 +
	- 2muu32(-2p1k1 + 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) + 2p1k1x2 +
	- sqrtS*2x1(-1 + x1 + x2))HL(i,tt,gen) +
	- muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2 + sqrtS2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HR(i,tt,gen))) +
	- sqrtS*2DR(j)(4muu1muu3GL(i,tt,gen)(muu1muu3sqrtS2(2p1k1 - 2p1k2 +
	- sqrtS*2(-1 - 4muu32 + x1 + x2))HL(i,tt,gen) +
	- (-(muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2)) + p1k1(-4p1k2 + sqrtS2(1 + x1)) +
	- muu1*2sqrtS*4(-1 + x2))*HR(i,tt,gen)) +
	- GR(i,tt,gen)(2muu1muu3(sqrtS*2(-1 + x1)(2p1k2 + sqrtS*2(2muu12 - x2)) + p1k1(-8p1k2 + 4sqrtS*2x2) -
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HL(i,tt,gen) +
	- (8muu34sqrtS*4 + p1k1(-16p1k2 + 2sqrtS*2(3 + x1)*x2) +
	- sqrtS*2(-1 + x1)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) +
	- muu3*2sqrtS*2(4p1k2(2 + x1) - 4p1k1(2 + x2) - 2sqrtS2(-2 + 2x1 + 3x2)) -
	- 2muu12sqrtS*2(8muu34sqrtS*2 - 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) + 2p1k1x2 +
	- muu3*2(-4p1k1 + 4p1k2 - 6sqrtS2(-1 + x1 + x2))))*HR(i,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(-DenSmH1(j))(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(i,tt,gen)**2))

	** USchannel dipole with H 312

	radUS = radUS - (4CFGS*2NC((-1 + x2)(2 + 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) +
	- x1*2dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- 3x1(-1 + x2)dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 5x2dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 2x22dsqrt((1 - 4muu32)/(-4muu32 + x22))) +
	- 4muu34(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) + x1dsqrt((1 - 4muu3**2)/
	- (-4muu32 + x2*2)) +
	- x2(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x22)))) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x22)) - 2x1dsqrt((1 - 4muu3*2)/
	- (-4muu32 + x2*2)) +
	- 3x22dsqrt((1 - 4muu32)/(-4muu32 + x22)) +
	- x2(4 - 8dsqrt((1 - 4muu32)/(-4muu32 + x22)) + 3x1dsqrt((1 - 4muu32)/(-4muu32 + x22)))))*
	- (-4dsqrt(1 - 4muu3*2)p1k2(CL(j)(muu1muu3DR(j)(GR(i,tt,gen)HL(i,tt,gen) + GL(i,tt,gen)*HR(i,tt,gen)) +
	- DL(j)(GL(i,tt,gen)HL(i,tt,gen) + muu1muu3GR(i,tt,gen)HL(i,tt,gen) + muu1muu3GL(i,tt,gen)HR(i,tt,gen))) +
	- CR(j)(muu1muu3DL(j)(GR(i,tt,gen)HL(i,tt,gen) + GL(i,tt,gen)HR(i,tt,gen)) +
	- DR(j)(muu1muu3GR(i,tt,gen)HL(i,tt,gen) + muu1muu3GL(i,tt,gen)HR(i,tt,gen) + GR(i,tt,gen)HR(i,tt,gen)))) +
	- sqrtS*2(CL(j)(DL(j)(muu1muu3GR(i,tt,gen)((dsqrt(1 - 4muu3*2)x2 + dsqrt(-4muu32 + x22))HL(i,tt,gen) +
	- 4muu1muu3dsqrt(-4muu32 + x22)*HR(i,tt,gen)) +
	- GL(i,tt,gen)((dsqrt(1 - 4muu3*2)x2 + (-1 + 2muu12)(-1 + 2muu32)dsqrt(-4muu32 + x22))HL(i,tt,gen) +
	- muu1muu3(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22))HR(i,tt,gen))) +
	- DR(j)(muu1GR(i,tt,gen)((muu3dsqrt(1 - 4muu32)x2 - muu3dsqrt(-4muu32 + x22))*HL(i,tt,gen) +
	- 2muu1(-1 + 2muu32)dsqrt(-4muu32 + x22)HR(i,tt,gen)) +
	- muu3GL(i,tt,gen)(2(-1 + 2muu1*2)muu3dsqrt(-4muu32 + x22)*HL(i,tt,gen) +
	- (muu1dsqrt(1 - 4muu3*2)x2 - muu1dsqrt(-4muu32 + x22))*HR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu3GR(i,tt,gen)((muu1dsqrt(1 - 4muu3*2)x2 - muu1dsqrt(-4muu32 + x22))*HL(i,tt,gen) +
	- 2(-1 + 2muu1*2)muu3dsqrt(-4muu32 + x22)*HR(i,tt,gen)) +
	- muu1GL(i,tt,gen)(2muu1(-1 + 2muu32)dsqrt(-4muu32 + x22)HL(i,tt,gen) +
	- (muu3dsqrt(1 - 4muu3*2)x2 - muu3dsqrt(-4muu32 + x22))*HR(i,tt,gen))) +
	- DR(j)(muu1muu3GL(i,tt,gen)(4muu1muu3dsqrt(-4muu32 + x22)*HL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22))HR(i,tt,gen)) +
	- GR(i,tt,gen)(muu1muu3(dsqrt(1 - 4muu3*2)x2 + dsqrt(-4muu32 + x22))HL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x2 + (-1 + 2muu12)(-1 + 2muu32)dsqrt(-4muu32 + x22))HR(i,tt,gen)))))))/
	- ((-1 + 2muu32)(-1 + x2)*2(-2 + x1 + x2)(DenSmH1(j))
	- (-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x2*2)) +
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)*2))


	** USchannel dipole with H 321

	radUS = radUS - (4CFGS*2NC(4muu3*4(-2 - 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) +
	- x1(2 + dsqrt((1 - 4muu3*2)/(-4muu32 + x12))) +
	- dsqrt((1 - 4muu32)/(-4muu32 + x12))x2) + (-1 + x1)
	- (2 + 2dsqrt((-1 + 4muu3*2)/(4muu32 - x12)) + 2x12dsqrt((1 - 4muu32)/(-4muu32 + x12)) -
	- 3dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2 + dsqrt((-1 + 4muu3*2)/(4muu32 - x12))x2*2 +
	- x1dsqrt((1 - 4muu3*2)/(-4muu32 + x12))(-5 + 3x2)) -
	- 2muu32(-4 + 4dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3x12dsqrt((1 - 4muu32)/(-4muu32 + x12)) -
	- 2dsqrt((1 - 4muu3*2)/(-4muu32 + x12))x2 + x1
	- (4 - 8dsqrt((1 - 4muu3*2)/(-4muu32 + x12)) + 3dsqrt((1 - 4muu3*2)/(-4muu32 + x12))x2)))
	- (-4dsqrt(1 - 4muu3*2)p1k1(CL(j)(muu1muu3DR(j)(GR(i,tt,gen)HL(i,tt,gen) + GL(i,tt,gen)*HR(i,tt,gen)) +
	- DL(j)(GL(i,tt,gen)HL(i,tt,gen) + muu1muu3GR(i,tt,gen)HL(i,tt,gen) + muu1muu3GL(i,tt,gen)HR(i,tt,gen))) +
	- CR(j)(muu1muu3DL(j)(GR(i,tt,gen)HL(i,tt,gen) + GL(i,tt,gen)HR(i,tt,gen)) +
	- DR(j)(muu1muu3GR(i,tt,gen)HL(i,tt,gen) + muu1muu3GL(i,tt,gen)HR(i,tt,gen) + GR(i,tt,gen)HR(i,tt,gen)))) +
	- sqrtS*2(CL(j)(DL(j)(muu1muu3GR(i,tt,gen)(dsqrt(1 - 4muu3*2)x1HL(i,tt,gen) - dsqrt(-4muu32 + x12)*HL(i,tt,gen) -
	- 4muu1muu3dsqrt(-4muu32 + x12)*HR(i,tt,gen)) +
	- GL(i,tt,gen)((dsqrt(1 - 4muu3*2)x1 - (-1 + 2muu12)(-1 + 2muu32)dsqrt(-4muu32 + x12))HL(i,tt,gen) +
	- muu1muu3(dsqrt(1 - 4muu32)x1 - dsqrt(-4muu32 + x12))HR(i,tt,gen))) +
	- DR(j)(muu1GR(i,tt,gen)(muu3(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12))HL(i,tt,gen) +
	- 2muu1(1 - 2muu32)dsqrt(-4muu32 + x12)HR(i,tt,gen)) +
	- muu3GL(i,tt,gen)(2(1 - 2muu1*2)muu3dsqrt(-4muu32 + x12)*HL(i,tt,gen) +
	- muu1(dsqrt(1 - 4muu3*2)x1 + dsqrt(-4muu32 + x12))HR(i,tt,gen)))) +
	- CR(j)(DL(j)(muu3GR(i,tt,gen)(muu1(dsqrt(1 - 4muu3*2)x1 + dsqrt(-4muu32 + x12))HL(i,tt,gen) +
	- 2(1 - 2muu1*2)muu3dsqrt(-4muu32 + x12)*HR(i,tt,gen)) +
	- muu1GL(i,tt,gen)(2muu1(1 - 2muu32)dsqrt(-4muu32 + x12)HL(i,tt,gen) +
	- muu3(dsqrt(1 - 4muu3*2)x1 + dsqrt(-4muu32 + x12))HR(i,tt,gen))) +
	- DR(j)(muu1muu3GL(i,tt,gen)(-4muu1muu3dsqrt(-4muu32 + x12)HL(i,tt,gen) + dsqrt(1 - 4muu3*2)x1*HR(i,tt,gen) -
	- dsqrt(-4muu32 + x12)HR(i,tt,gen)) + GR(i,tt,gen)*
	- (muu1muu3(dsqrt(1 - 4muu32)x1 - dsqrt(-4muu32 + x12))HL(i,tt,gen) +
	- (dsqrt(1 - 4muu32)x1 - (-1 + 2muu12)(-1 + 2muu32)dsqrt(-4muu32 + x12))HR(i,tt,gen)))))))/
	- ((-1 + 2muu32)(-1 + x1)*2(-2 + x1 + x2)(DenSmH1(j))
	- (-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x1*2)) -
	- 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)*2))


	enddo
	enddo

	end


	************************************************

	double precision function radTU()

	#include "ChiChi2QQ_Kinematics.h"
	#include "ChiChi2QQ_Model.h"
	#include "ChiChi2QQ_GenCouplings.h"

	integer i,j

	radTU = 0d0

	do i = 1, 2
	do j = 1, 2

	** TUchannel bremsstrahlung 11

	radTU = radTU + (2CFgs*2NC(ELc(i,tt,gen)(2FR(i,tt,gen)(muu1sqrtSGR(j,tt,gen)*
	- (muu1sqrtS5(4muu34 + (-1 + x1)(-1 + x2) - 2muu32x2)*HL(j,tt,gen) -
	- 2muu3p1k2sqrtS3(1 + 2muu32 - x2)HR(j,tt,gen)) +
	- muu3sqrtSGL(j,tt,gen)(muu1sqrtS*3(-2p1k1(-1 + x2) - 2p1k2(-1 + 2muu3*2 + x2) +
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))
	- HL(j,tt,gen) + (-1 + 2muu12)muu3sqrtS5(1 + 2muu32 - x2)HR(j,tt,gen))) +
	- FL(i,tt,gen)(-2muu1muu3sqrtS*2GR(j,tt,gen)(sqrtS2(-2p1k1(-1 + x2) -
	- (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))HL(j,tt,gen) +
	- 2muu1muu3sqrtS4(1 + 2muu32 - x2)HR(j,tt,gen)) +
	- GL(j,tt,gen)(sqrtS2(-4muu34sqrtS*4 + 2muu1*2sqrtS*4(4muu3*4 +
	- (-1 + x1)(-1 + x2) - 2muu3*2x2) +
	- 8muu32p1k2(-2p1k2 + sqrtS*2x2) - (-1 + x2)*
	- (8p1k22 + sqrtS4(-1 + x1 + x2) - 2p1k2sqrtS*2(x1 + 2*x2) +
	- p1k1(8p1k2 - 2sqrtS2x2)))*HL(j,tt,gen) -
	- 2muu1muu3sqrtS4(1 + 2muu32 - x2)(-2p1k2 + sqrtS2x2)*HR(j,tt,gen)))) -
	- ER(i,tt,gen)(-2FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)
	- (-2muu3p1k2sqrtS3(1 + 2muu32 - x2)HL(j,tt,gen) + muu1sqrtS5(4muu3*4 +
	- (-1 + x1)(-1 + x2) - 2muu3*2x2)*HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)((-1 + 2muu1*2)muu3sqrtS5(1 + 2muu32 - x2)HL(j,tt,gen) +
	- muu1sqrtS3(-2p1k1(-1 + x2) - 2p1k2(-1 + 2muu3*2 + x2) +
	- sqrtS*2(-1 + x2)(-1 + 2muu3*2 + x1 + x2))HR(j,tt,gen))) +
	- FR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen)(2muu1muu3sqrtS*4(1 + 2muu32 - x2)HL(j,tt,gen) +
	- sqrtS*2(-2p1k1(-1 + x2) - (2p1k2 - sqrtS2)(-1 + 2muu32 + x2))HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*4(1 + 2muu32 - x2)(-2p1k2 + sqrtS2x2)*HL(j,tt,gen) +
	- sqrtS*2(4muu34sqrtS*4 - 2muu1*2sqrtS*4(4muu3*4 +
	- (-1 + x1)(-1 + x2) - 2muu3*2x2) + 8muu32p1k2(2p1k2 - sqrtS*2x2) +
	- (-1 + x2)(8p1k22 + sqrtS4(-1 + x1 + x2) - 2p1k2sqrtS2(x1 + 2x2) + p1k1
	- (8p1k2 - 2sqrtS*2x2)))*HR(j,tt,gen))))))/
	- (sqrtS*4(-1 + x2)*2(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)**2))



	** TUchannel bremsstrahlung 12

	radTU = radTU + (2CFgs*2NC(ELc(i,tt,gen)(FL(i,tt,gen)(muu1muu3sqrtS2GR(j,tt,gen)*
	- (sqrtS*2(muu3*2(8p1k2 - 4sqrtS*2) - 2p1k2x1 - 2p1k1*(-2 + x2) +
	- sqrtS*2(-1 + x1 + x2))*HL(j,tt,gen) +
	- 2muu1muu3sqrtS4(-4muu32 + x1 + x2)HR(j,tt,gen)) +
	- GL(j,tt,gen)(2(-2sqrtS2(2p1k1p1k2 + muu3*4sqrtS**4 +
	- muu3*2(p1k2sqrtS2(-2 + x1) + p1k1(-4p1k2 + sqrtS*2(-2 + x2)))) +
	- muu1*2sqrtS*6(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)))*HL(j,tt,gen) +
	- muu1muu3sqrtS*4(2p1k2(-2 + x1) - sqrtS*2(-3 + 2x1)(-1 + x1 + x2) +
	- 2p1k1(-2 + 2x1 + x2) + muu32(8p1k2 - 4sqrtS*2x2))*
	- HR(j,tt,gen))) + FR(i,tt,gen)(muu1sqrtSGR(j,tt,gen)
	- (2muu1sqrtS*5(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2))*HL(j,tt,gen) +
	- muu3sqrtS3(muu3*2(8p1k1 - 4sqrtS*2) - 2p1k2(-2 + x1) - 2p1k1*x2 +
	- sqrtS*2(-1 + x1 + x2))*HR(j,tt,gen)) +
	- muu3sqrtSGL(j,tt,gen)(muu1sqrtS*3(muu3*2(8p1k1 - 4sqrtS*2x1) + 2p1k1(-2 + x2) -
	- sqrtS*2(-1 + x1 + x2)(-3 + 2x2) +
	- 2p1k2(-2 + x1 + 2x2))HL(j,tt,gen) + 2muu3sqrtS*
	- (4p1k12 + 4p1k2*2 - 2p1k2sqrtS2(-1 + x1 + x2) +
	- sqrtS*4(-1 + (-2 + 4muu12)muu32 + x1 - muu12x1 + x2 - muu12x2) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)))*HR(j,tt,gen)))) +
	- ER(i,tt,gen)(FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)(muu3sqrtS*3
	- (muu3*2(8p1k1 - 4sqrtS*2) - 2p1k2(-2 + x1) - 2p1k1x2 + sqrtS2(-1 + x1 + x2))*HL(j,tt,gen) +
	- 2muu1sqrtS*5(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2))*HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)(2muu3sqrtS(4p1k12 + 4p1k2*2 - 2p1k2sqrtS2(-1 + x1 + x2) +
	- sqrtS*4(-1 + (-2 + 4muu12)muu32 + x1 - muu12x1 + x2 - muu12x2) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)))*HL(j,tt,gen) +
	- muu1sqrtS3(muu3*2(8p1k1 - 4sqrtS*2x1) + 2p1k1(-2 + x2) -
	- sqrtS*2(-1 + x1 + x2)(-3 + 2x2) + 2p1k2(-2 + x1 + 2x2))
	- HR(j,tt,gen))) + FR(i,tt,gen)(muu1muu3sqrtS2GL(j,tt,gen)*
	- (2muu1muu3sqrtS4(-4muu32 + x1 + x2)HL(j,tt,gen) +
	- sqrtS*2(muu3*2(8p1k2 - 4sqrtS*2) - 2p1k2x1 - 2p1k1*(-2 + x2) +
	- sqrtS*2(-1 + x1 + x2))*HR(j,tt,gen)) +
	- GR(j,tt,gen)(muu1muu3sqrtS4(2p1k2(-2 + x1) - sqrtS*2(-3 + 2x1)(-1 + x1 + x2) +
	- 2p1k1(-2 + 2*x1 + x2) +
	- muu3*2(8p1k2 - 4sqrtS*2x2))*HL(j,tt,gen) +
	- 2(-2sqrtS*2(2p1k1p1k2 + muu3*4sqrtS4 + muu32(p1k2sqrtS*2(-2 + x1) +
	- p1k1(-4p1k2 + sqrtS*2(-2 + x2)))) +
	- muu1*2sqrtS*6(-1 + 4muu34 + x1 + x2 - muu32(2 + x1 + x2)))*HR(j,tt,gen))))))/
	- (sqrtS*4(-1 + x1)(-1 + x2)(-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)*
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2))



	** TUchannel bremsstrahlung 13

	radTU = radTU + (CFgs2NC(ELc(i,tt,gen)(-2FR(i,tt,gen)(muu3sqrtSGL(j,tt,gen)*
	- (muu1sqrtS(2p1k1(4p1k2 + sqrtS2(1 - 2x2)) + 2muu3*2sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 2muu12(-1 + x2) + x2 + x1x2)))HL(j,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + muu12(-4p1k1 + 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - x1 - 3x2)) + 2p1k1x2 +
	- sqrtS*2(-2 + x1 + 2x2))HR(j,tt,gen)) +
	- muu1sqrtSGR(j,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 - 2p1k2x1 + x2(2p1k1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2)))HL(j,tt,gen) +
	- muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2(1 + x1 - x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HR(j,tt,gen))) +
	- FL(i,tt,gen)(-4muu1muu3sqrtS*2GR(j,tt,gen)((-(muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2)) +
	- muu1*2sqrtS*4(-1 + x2) +
	- p1k1(-4p1k2 + sqrtS*2(1 + x2)))HL(j,tt,gen) - muu1muu3sqrtS2
	- (-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - 2x2))*HR(j,tt,gen)) +
	- GL(j,tt,gen)((8(1 - 2muu12)muu3*4sqrtS**6 -
	- 2p1k1(16p1k22 + 4p1k2sqrtS2(1 + 2muu32 - 3x2) +
	- sqrtS*4(-2muu32(-2 + 2muu12 + x2) + x2(-1 + 2muu1*2 + x2))) +
	- sqrtS*4((-1 + x2)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) + 2muu12
	- (2p1k2(-2 + x1 + 2x2) - sqrtS2(-1 + x1 + x2)(-2 + 3x2))) +
	- 2muu32(8p1k22sqrtS*2 + 2p1k2sqrtS4(-2muu12 + x1 - 2x2) +
	- sqrtS*6(2 - x1 - 2x2 + muu12(-6 + 4x1 + 8x2))))*
	- HL(j,tt,gen) + 2muu1muu3sqrtS2(-8p1k22 + 2p1k2sqrtS2(-3 + 2x1 + 5x2) +
	- sqrtS*4(2muu12(-1 + x1) + (3 - 2x1 - 3x2)x2) - 2muu3*2sqrtS*2(2p1k1 + 6p1k2 -
	- sqrtS*2(-1 + x1 + 3x2)))HR(j,tt,gen)))) +
	- ER(i,tt,gen)(-2FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)
	- (muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2(1 + x1 - x2) + sqrtS*2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HL(j,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 - 2p1k2x1 + x2(2p1k1 + sqrtS2(-1 + x1 + x2)) -
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(-1 + x1 + 3x2)))HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 - 2p1k2x1 + muu12
	- (-4p1k1 + 4p1k2 + 2sqrtS2(3 + 4muu32 - x1 - 3x2)) +
	- 2p1k1x2 + sqrtS*2(-2 + x1 + 2x2))HL(j,tt,gen) +
	- muu1sqrtS(2p1k1(4p1k2 + sqrtS2(1 - 2*x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 - 2muu12(-1 + x2) + x2 + x1x2)))HR(j,tt,gen))) +
	- FR(i,tt,gen)(4muu1muu3sqrtS*2GL(j,tt,gen)(muu1muu3sqrtS2(-2p1k1 + 2p1k2 +
	- sqrtS*2(1 + 4muu32 - 2x2))*HL(j,tt,gen) +
	- (muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS2) - muu12sqrtS4(-1 + x2) +
	- p1k1(4p1k2 - sqrtS*2(1 + x2)))*HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*2(-8p1k22 + 2p1k2sqrtS2(-3 + 2x1 + 5x2) +
	- sqrtS*4(2muu12(-1 + x1) + (3 - 2x1 - 3x2)*x2) -
	- 2muu32sqrtS*2(2p1k1 + 6p1k2 - sqrtS*2(-1 + x1 + 3x2)))HL(j,tt,gen) +
	- (8(1 - 2muu1*2)muu3*4sqrtS*6 - 2p1k1*
	- (16p1k22 + 4p1k2sqrtS2(1 + 2muu32 - 3x2) +
	- sqrtS*4(-2muu32(-2 + 2muu12 + x2) + x2(-1 + 2muu1*2 + x2))) +
	- sqrtS*4((-1 + x2)(-2p1k2x1 + sqrtS2(-1 + x1 + x2)) + 2muu12(2p1k2(-2 + x1 + 2*x2) -
	- sqrtS*2(-1 + x1 + x2)(-2 + 3x2))) +
	- 2muu32(8p1k22sqrtS*2 + 2p1k2sqrtS4(-2muu12 + x1 - 2x2) +
	- sqrtS*6(2 - x1 - 2x2 + muu12(-6 + 4x1 + 8x2))))*
	- HR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x2)(2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)
	- (-2p1k2 - sqrtS2(muu12 + muu32 - x2) + MSf(j,tt,gen)*2)
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2))


	** TUchannel bremsstrahlung 21

	radTU = radTU + (-2CFgs*2NC(ELc(i,tt,gen)(FL(i,tt,gen)(muu1muu3sqrtS2GR(j,tt,gen)*
	- (sqrtS*2(2p1k1(-2 + 4muu32 + x2) + 2p1k2(-2 + x1 + 2x2) -
	- sqrtS*2(-3 + x1 + 3x2))HL(j,tt,gen) +
	- 2muu1muu3sqrtS4(4muu32 - x1 - x2)HR(j,tt,gen)) +
	- GL(j,tt,gen)(2sqrtS*2(2muu34sqrtS*4 + (-2p1k1 + sqrtS*2x1)(-2p1k2 + sqrtS*2x2) +
	- muu1*2sqrtS*4(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) +
	- muu3*2(p1k1(-8p1k2 + 2sqrtS2(2 + x2)) - 2sqrtS2(-(p1k2*(2 + x1)) +
	- sqrtS*2(x1 + x2))))*HL(j,tt,gen) +
	- muu1muu3sqrtS*4(sqrtS*2 + 4muu3*2(2p1k1 - sqrtS2(-1 + x1)) - 2p1k2(-2 + x1) -
	- sqrtS*2x1 - 2p1k1x2 - 3sqrtS2x2 +
	- 2sqrtS2x1x2)HR(j,tt,gen))) + FR(i,tt,gen)*
	- (muu1sqrtSGR(j,tt,gen)(2muu1sqrtS5(1 - 4muu3*4 - x1 - x2 +
	- muu3*2(2 + x1 + x2))*HL(j,tt,gen) +
	- muu3sqrtS3(2p1k2(-2 + 4muu32 + x1) + 2p1k1(-2 + 2x1 + x2) -
	- sqrtS*2(-3 + 3x1 + x2))HR(j,tt,gen)) +
	- muu3sqrtSGL(j,tt,gen)(muu1sqrtS*3(sqrtS*2 - 2p1k2x1 - 3sqrtS*2x1 - 2p1k1(-2 + x2) -
	- sqrtS*2x2 + 2sqrtS2x1*x2 +
	- 4muu32(2p1k2 + sqrtS2 - sqrtS2x2))*HL(j,tt,gen) +
	- 2muu3sqrtS(-4p1k1*2 - 4p1k2*2 + 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(1 + (2 - 4muu12)muu32 + (-2 + muu12)x1 - 2x2 + muu1*2x2) +
	- 2p1k1(-4p1k2 + sqrtS2(1 + x1 + x2)))*HR(j,tt,gen)))) +
	- ER(i,tt,gen)(FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)*
	- (muu3sqrtS3(2p1k2(-2 + 4muu32 + x1) + 2p1k1(-2 + 2x1 + x2) -
	- sqrtS*2(-3 + 3x1 + x2))HL(j,tt,gen) +
	- 2muu1sqrtS*5(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2))*HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)(2muu3sqrtS(-4p1k12 - 4p1k2*2 + 2p1k2sqrtS2(1 + x1 + x2) +
	- sqrtS*4(1 + (2 - 4muu12)muu32 + (-2 + muu12)x1 - 2x2 + muu1*2x2) +
	- 2p1k1(-4p1k2 + sqrtS2(1 + x1 + x2)))*HL(j,tt,gen) +
	- muu1sqrtS3(sqrtS*2 - 2p1k2x1 - 3sqrtS*2x1 - 2p1k1(-2 + x2) -
	- sqrtS*2x2 + 2sqrtS2x1*x2 +
	- 4muu32(2p1k2 + sqrtS2 - sqrtS2x2))*HR(j,tt,gen))) +
	- FR(i,tt,gen)(muu1muu3sqrtS2GL(j,tt,gen)(2muu1muu3sqrtS*4(4muu32 - x1 - x2)HL(j,tt,gen) +
	- sqrtS*2(2p1k1(-2 + 4muu32 + x2) + 2p1k2(-2 + x1 + 2x2) -
	- sqrtS*2(-3 + x1 + 3x2))HR(j,tt,gen)) +
	- GR(j,tt,gen)(muu1muu3sqrtS4(sqrtS*2 + 4muu3*2(2p1k1 - sqrtS2(-1 + x1)) -
	- 2p1k2(-2 + x1) - sqrtS*2x1 - 2p1k1x2 -
	- 3sqrtS2x2 + 2sqrtS2x1x2)HL(j,tt,gen) +
	- 2sqrtS2(2muu34sqrtS*4 + (-2p1k1 + sqrtS*2x1)(-2p1k2 + sqrtS*2x2) +
	- muu1*2sqrtS*4(1 - 4muu34 - x1 - x2 + muu32(2 + x1 + x2)) +
	- muu3*2(p1k1(-8p1k2 + 2sqrtS2(2 + x2)) - 2sqrtS2(-(p1k2*(2 + x1)) +
	- sqrtS*2(x1 + x2))))*HR(j,tt,gen))))))/
	- (sqrtS*4(-1 + x1)(-1 + x2)(-2p1k1 - sqrtS2(muu12 + muu32 - x1) + MSf(i,tt,gen)*2)
	- (-2p1k2 - sqrtS2(muu12 + muu32 - x2) + MSf(j,tt,gen)**2))


	** TUchannel bremsstrahlung 22

	radTU = radTU + (2CFgs*2NC(ELc(i,tt,gen)(2FR(i,tt,gen)(-(muu3sqrtS3(1 + 2muu32 - x1)GL(j,tt,gen)*
	- (muu1(-2p1k1sqrtS + sqrtS3x1)HL(j,tt,gen) + (1 - 2muu1*2)muu3sqrtS3HR(j,tt,gen))) +
	- muu1sqrtSGR(j,tt,gen)(muu1sqrtS*5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))HL(j,tt,gen) +
	- muu3sqrtS3(muu3*2(4p1k1 - 2sqrtS*2) + (2p1k1 + 2p1k2 - sqrtS2)(-1 + x1))*HR(j,tt,gen))) +
	- FL(i,tt,gen)(-4muu1muu3sqrtS*4(1 + 2muu32 - x1)GR(j,tt,gen)(p1k1HL(j,tt,gen) +
	- muu1muu3sqrtS*2HR(j,tt,gen)) +
	- GL(j,tt,gen)(sqrtS2(-4muu34sqrtS*4 - 8muu3*2p1k1(2p1k1 - sqrtS*2x1) +
	- 2muu12sqrtS*4(4muu34 - 2muu3*2x1 + (-1 + x1)*(-1 + x2)) -
	- (-1 + x1)(8p1k1*2 - 2p1k2sqrtS2x1 + sqrtS*4(-1 + x1 + x2) +
	- p1k1(8p1k2 - 2sqrtS2(2x1 + x2))))HL(j,tt,gen) +
	- 2muu1muu3sqrtS4(-2p1k1(-1 + 2muu32 + x1) + (-1 + x1)(-2*p1k2 +
	- sqrtS*2(-1 + 2muu32 + x1 + x2)))HR(j,tt,gen)))) -
	- ER(i,tt,gen)(-2FL(i,tt,gen)(-(muu3sqrtS*3(1 + 2muu32 - x1)GR(j,tt,gen)*
	- ((1 - 2muu12)muu3sqrtS3HL(j,tt,gen) + muu1sqrtS(-2p1k1 + sqrtS2x1)*HR(j,tt,gen))) +
	- muu1sqrtSGL(j,tt,gen)(muu3sqrtS*3(muu3*2(4p1k1 - 2sqrtS**2) +
	- (2p1k1 + 2p1k2 - sqrtS*2)(-1 + x1))*HL(j,tt,gen) +
	- muu1sqrtS5(4muu34 - 2muu3*2x1 + (-1 + x1)(-1 + x2))HR(j,tt,gen))) +
	- FR(i,tt,gen)(4muu1muu3sqrtS*4(1 + 2muu32 - x1)GL(j,tt,gen)*
	- (muu1muu3sqrtS*2HL(j,tt,gen) + p1k1*HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*4(2p1k1(-1 + 2muu3*2 + x1) +
	- (-1 + x1)(2p1k2 - sqrtS*2(-1 + 2muu32 + x1 + x2)))HL(j,tt,gen) +
	- sqrtS*2(4muu34sqrtS*4 + 8muu3*2p1k1(2p1k1 - sqrtS*2x1) -
	- 2muu12sqrtS*4(4muu34 - 2muu3*2x1 + (-1 + x1)*(-1 + x2)) +
	- (-1 + x1)(8p1k1*2 - 2p1k2sqrtS2x1 +
	- sqrtS*4(-1 + x1 + x2) + p1k1(8p1k2 - 2sqrtS2(2x1 + x2))))HR(j,tt,gen))))))/
	- (sqrtS*4(-1 + x1)*2(2p1k1 + sqrtS2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2))

	** TUchannel bremsstrahlung 23

	radTU = radTU + (CFgs2NC(ELc(i,tt,gen)(FL(i,tt,gen)(4muu1muu3sqrtS*2GR(j,tt,gen)*
	- ((-4p1k12 + muu32sqrtS*2(6p1k1 + 2p1k2 - sqrtS2) + muu12sqrtS4(-1 + x2) -
	- p1k1sqrtS2(-3 + x1 + 2x2))HL(j,tt,gen) +
	- muu1muu3sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))HR(j,tt,gen)) +
	- GL(j,tt,gen)((32p1k1*2p1k2 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu3*4sqrtS**6 -
	- 8p1k1p1k2sqrtS2x1 + 2p1k2sqrtS*4x1 -
	- 2sqrtS6x1 - 2p1k2sqrtS*4x12 + sqrtS6x12 - 16p1k1*2sqrtS*2x2 +
	- 2p1k1sqrtS*4x2 - sqrtS*6x2 +
	- 6p1k1sqrtS*4x1x2 + sqrtS6x1x2 + 2muu3*2sqrtS*2
	- (8p1k12 + sqrtS2(2p1k2(2 + x1) + sqrtS*2(2 - 2x1 - 3x2)) + p1k1(-8p1k2 +
	- 2sqrtS2(-2*x1 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 - 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) +
	- 2p1k1(-2 + 2*x1 + x2) +
	- muu3*2(-4p1k1 + 4p1k2 - 6sqrtS2(-1 + x1 + x2))))*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(p1k1(8p1k2 - 4sqrtS2x2) + sqrtS*2(-1 + x1)(-2p1k2 +
	- sqrtS*2(-2muu1*2 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HR(j,tt,gen))) -
	- 2FR(i,tt,gen)(muu3sqrtSGL(j,tt,gen)(muu1sqrtS*
	- (2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HL(j,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 - 2p1k2x1 + muu12(-4p1k1 + 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) + 2p1k1*x2 +
	- sqrtS*2(-2 + 2x1 + x2))HR(j,tt,gen)) +
	- muu1sqrtSGR(j,tt,gen)(muu1sqrtS*3(8muu34sqrtS*2 - 2p1k2*x1 +
	- 2muu32(-2p1k1 + 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) + 2p1k1x2 +
	- sqrtS*2x1(-1 + x1 + x2))HL(j,tt,gen) +
	- muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2 + sqrtS2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HR(j,tt,gen)))) +
	- ER(i,tt,gen)(-2FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)
	- (muu3sqrtS(-2muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS*2) + p1k1(8p1k2 - 2sqrtS*2x2) +
	- sqrtS*2(-2p1k2 + sqrtS2(-1 - 2muu12(-1 + x1) + x1 + x2)))*HL(j,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 - 2p1k2x1 + 2muu3*2(-2p1k1 + 2p1k2 +
	- sqrtS*2(1 - 2x1 - 2x2)) + 2p1k1x2 +
	- sqrtS*2x1(-1 + x1 + x2))HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 - 2p1k2*x1 +
	- muu1*2(-4p1k1 + 4p1k2 + 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) + 2p1k1*x2 +
	- sqrtS*2(-2 + 2x1 + x2))HL(j,tt,gen) +
	- muu1sqrtS(2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)) +
	- 2p1k1(4p1k2 - sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(-2p1k2x1 + sqrtS*2(-1 + x1*2 - 2muu1*2(-1 + x2) + x2)))*HR(j,tt,gen))) +
	- FR(i,tt,gen)(4muu1muu3sqrtS*2GL(j,tt,gen)(muu1muu3sqrtS2
	- (-2p1k1 + 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))HL(j,tt,gen) +
	- (-4p1k12 + muu32sqrtS*2(6p1k1 + 2p1k2 - sqrtS2) + muu12sqrtS4(-1 + x2) -
	- p1k1sqrtS2(-3 + x1 + 2x2))HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*2(p1k1(8p1k2 - 4sqrtS2x2) +
	- sqrtS*2(-1 + x1)(-2p1k2 + sqrtS*2(-2muu1*2 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)))*HL(j,tt,gen) +
	- (32p1k12p1k2 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu3*4sqrtS*6 - 8p1k1p1k2sqrtS*2x1 +
	- 2p1k2sqrtS*4x1 - 2sqrtS6x1 -
	- 2p1k2sqrtS*4x12 + sqrtS6x12 - 16p1k1*2sqrtS*2x2 + 2p1k1sqrtS*4x2 - sqrtS*6x2 +
	- 6p1k1sqrtS*4x1*x2 +
	- sqrtS*6x1x2 + 2muu3*2sqrtS*2(8p1k12 + sqrtS2(2p1k2(2 + x1) + sqrtS*2(2 - 2x1 - 3x2)) +
	- p1k1(-8p1k2 + 2sqrtS2(-2*x1 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 - 2p1k2x1 + sqrtS2(2 + x1*2 + 3x1(-1 + x2) - 2x2) +
	- 2p1k1(-2 + 2*x1 + x2) +
	- muu3*2(-4p1k1 + 4p1k2 - 6sqrtS2(-1 + x1 + x2))))*HR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(2p1k2 + sqrtS*2(muu12 + muu32 - x2) - MSf(j,tt,gen)*2)
	- (-2p1k1 - sqrtS2(muu12 + muu32 - x1) + MSf(i,tt,gen)*2)
	- (2p1k1 - (muu12 + muu32)sqrtS2 + MSf(j,tt,gen)2))


	** TUchannel bremsstrahlung 31

	radTU = radTU + (CFgs2NC(ELc(i,tt,gen)(-2FR(i,tt,gen)(muu1sqrtSGR(j,tt,gen)*
	- (muu1sqrtS3(8muu34sqrtS*2 + 2p1k2x1 + 2muu3*2(2p1k1 - 2p1k2 +
	- sqrtS*2(1 - 2x1 - 2x2)) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2)))*HL(j,tt,gen) +
	- 2muu3sqrtS(4p1k22 + muu32sqrtS2(-2p1k1 - 6p1k2 + sqrtS**2) -
	- muu1*2sqrtS*4(-1 + x1) + p1k2sqrtS2(-3 + 2x1 + x2))
	- HR(j,tt,gen)) + muu3sqrtSGL(j,tt,gen)(muu1sqrtS*
	- (-2p1k1(4p1k2 - sqrtS2(-1 + x2)) + sqrtS*2(x1(4p1k2 - sqrtS*2(-1 + x2)) +
	- 2muu12sqrtS*2(-1 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)))*HL(j,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) - 2p1k1*x2 +
	- sqrtS*2(-2 + x1 + 2x2))HR(j,tt,gen))) +
	- FL(i,tt,gen)(2muu1muu3sqrtS*2GR(j,tt,gen)((2muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2) +
	- p1k1(-8p1k2 + 2sqrtS*2) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - x1 - x2))HR(j,tt,gen)) +
	- GL(j,tt,gen)((32p1k1p1k22 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu34sqrtS**6 -
	- 16p1k22sqrtS*2x1 + 2p1k2sqrtS*4x1 -
	- sqrtS*6x1 - 8p1k1p1k2sqrtS2x2 + 2p1k1sqrtS*4x2 - 2sqrtS6x2 +
	- 6p1k2sqrtS*4x1x2 + sqrtS6x1*x2 -
	- 2p1k1sqrtS*4x22 + sqrtS6x22 + 2muu3*2sqrtS*2
	- (8p1k22 + sqrtS4(2 - 3x1 - 2x2) + 2p1k2sqrtS*2(x1 - 2*x2) +
	- p1k1(-8p1k2 + 2sqrtS2(2 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 - 2p1k1x2 + 2p1k2(-2 + x1 + 2x2) +
	- muu3*2(4p1k1 - 4p1k2 - 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3x2))))HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(p1k1(-8p1k2 + 2sqrtS2x2) + 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(2p1k2(-1 + x1 + x2) + sqrtS*2(1 + 2muu12(-1 + x1) - x1 - x2*2)))HR(j,tt,gen)))) +
	- ER(i,tt,gen)(-2FL(i,tt,gen)(muu3sqrtSGR(j,tt,gen)
	- (muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 2x1 - 2x2)) - 2p1k1*x2 +
	- sqrtS*2(-2 + x1 + 2x2))HL(j,tt,gen) +
	- muu1sqrtS(-2p1k1(4p1k2 - sqrtS2(-1 + x2)) + sqrtS*2(x1(4p1k2 - sqrtS*2(-1 + x2)) +
	- 2muu12sqrtS*2(-1 + x2)) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(-1 - x1 + x2)))*HR(j,tt,gen)) +
	- muu1sqrtSGL(j,tt,gen)(2muu3sqrtS(4p1k22 + muu32sqrtS*2
	- (-2p1k1 - 6p1k2 + sqrtS2) - muu12sqrtS4(-1 + x1) +
	- p1k2sqrtS2(-3 + 2x1 + x2))HL(j,tt,gen) +
	- muu1sqrtS3(8muu34sqrtS*2 + 2p1k2*x1 +
	- 2muu32(2p1k1 - 2p1k2 + sqrtS*2(1 - 2x1 - 2x2)) +
	- x2(-2p1k1 + sqrtS*2(-1 + x1 + x2)))*HR(j,tt,gen))) +
	- FR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen)(2muu1muu3sqrtS*2(2p1k1 - 2p1k2 +
	- sqrtS*2(1 + 4muu32 - x1 - x2))HL(j,tt,gen) +
	- (2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2) + p1k1(-8p1k2 + 2sqrtS**2) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*2(p1k1(-8p1k2 + 2sqrtS2x2) +
	- 2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(2p1k2(-1 + x1 + x2) + sqrtS*2(1 + 2muu12(-1 + x1) - x1 - x2*2)))HL(j,tt,gen) +
	- (32p1k1p1k2*2 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu34sqrtS**6 -
	- 16p1k22sqrtS*2x1 + 2p1k2sqrtS*4x1 - sqrtS*6x1 -
	- 8p1k1p1k2sqrtS2x2 + 2p1k1sqrtS*4x2 - 2sqrtS6x2 + 6p1k2sqrtS*4x1*x2 +
	- sqrtS*6x1x2 - 2p1k1sqrtS4x2**2 +
	- sqrtS*6x2*2 + 2muu3*2sqrtS*2(8p1k22 + sqrtS4(2 - 3x1 - 2x2) + 2p1k2sqrtS*2(x1 - 2*x2) +
	- p1k1(-8p1k2 + 2sqrtS2(2 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 - 2p1k1x2 + 2p1k2(-2 + x1 + 2x2) +
	- muu3*2(4p1k1 - 4p1k2 - 6sqrtS2(-1 + x1 + x2)) +
	- sqrtS*2(2 - 3x2 + x22 + x1(-2 + 3x2))))HR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x2)(2p1k1 + sqrtS*2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)*
	- (-2p1k2 - sqrtS2(muu12 + muu32 - x2) + MSf(j,tt,gen)**2))

	** TUchannel bremsstrahlung 32

	radTU = radTU + (CFgs2NC(ELc(i,tt,gen)(-2FR(i,tt,gen)(muu1sqrtSGR(j,tt,gen)*
	- (muu1sqrtS3(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1x2 + sqrtS2x1*(-1 + x1 + x2) +
	- 2muu32(2p1k1 - 2p1k2 - sqrtS*2(-1 + 3x1 + x2)))HL(j,tt,gen) +
	- 2muu3sqrtS(-4p1k1p1k2 - muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2) +
	- sqrtS*2(muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x1)))
	- HR(j,tt,gen)) + muu3sqrtSGL(j,tt,gen)(muu1sqrtS*
	- (8p1k12 - 2p1k1sqrtS2(-3 + 5x1 + 2x2) + 2muu32sqrtS*2(6p1k1 + 2p1k2 -
	- sqrtS*2(-1 + 3*x1 + x2)) +
	- sqrtS*4(-2muu12(-1 + x2) + x1(-3 + 3x1 + 2x2)))HL(j,tt,gen) +
	- muu3sqrtS3(-4muu32sqrtS*2 + 2p1k2x1 + muu12(4p1k1 - 4p1k2 +
	- 2sqrtS2(3 + 4muu32 - 3x1 - x2)) - 2p1k1x2 +
	- sqrtS*2(-2 + 2x1 + x2))HR(j,tt,gen))) +
	- FL(i,tt,gen)(2muu1muu3sqrtS*2GR(j,tt,gen)((2muu3*2sqrtS*2(2p1k1 - 2p1k2 + sqrtS**2) +
	- 2p1k1(-4p1k2 + sqrtS2(1 - x1 + x2)) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - 2x1))*HR(j,tt,gen)) +
	- GL(j,tt,gen)((-32p1k1*2p1k2 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu3*4sqrtS**6 +
	- 24p1k1p1k2sqrtS2x1 + 2p1k2sqrtS*4x1 -
	- 2sqrtS6x1 - 2p1k2sqrtS*4x12 + sqrtS6x12 + 2p1k1sqrtS4x2 - sqrtS*6x2 -
	- 2p1k1sqrtS*4x1x2 + sqrtS6x1*x2 +
	- 2muu32sqrtS*2(8p1k12 + p1k1(-8p1k2 + 2sqrtS*2(-2*x1 + x2)) -
	- sqrtS*2(-2p1k2(-2 + x1) + sqrtS*2(-2 + 2*x1 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 + 2p1k2x1 + sqrtS2(-2 + 3x1)(-1 + x1 + x2) -
	- 2p1k1(-2 + 2*x1 + x2) +
	- muu3*2(4p1k1 - 4p1k2 - 2sqrtS2(-3 + 4x1 + 2x2))))*HL(j,tt,gen) +
	- 2muu1muu3sqrtS2(p1k1(-8p1k2 + 2sqrtS2x2) + 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(p1k2(-2 + 4x1) + sqrtS*2(1 + 2muu12(-1 + x1) - x1 - x1x2)))HR(j,tt,gen)))) +
	- ER(i,tt,gen)(-2FL(i,tt,gen)(muu1sqrtSGL(j,tt,gen)
	- (2muu3sqrtS(-4p1k1p1k2 - muu32sqrtS*2(-2p1k1 + 2p1k2 + sqrtS**2) +
	- sqrtS*2(muu1*2sqrtS*2(-1 + x1) + p1k2(1 + x1)))
	- HL(j,tt,gen) + muu1sqrtS3(8muu34sqrtS*2 + 2p1k2x1 - 2p1k1x2 + sqrtS2x1*(-1 + x1 + x2) +
	- 2muu32(2p1k1 - 2p1k2 - sqrtS*2(-1 + 3x1 + x2)))HR(j,tt,gen)) +
	- muu3sqrtSGR(j,tt,gen)(muu3sqrtS*3(-4muu32sqrtS*2 + 2p1k2x1 + muu12
	- (4p1k1 - 4p1k2 + 2sqrtS2(3 + 4muu32 - 3x1 - x2)) -
	- 2p1k1x2 + sqrtS*2(-2 + 2x1 + x2))HL(j,tt,gen) +
	- muu1sqrtS(8p1k12 - 2p1k1sqrtS2(-3 + 5x1 + 2x2) + 2muu32sqrtS*2
	- (6p1k1 + 2p1k2 - sqrtS*2(-1 + 3*x1 + x2)) +
	- sqrtS*4(-2muu12(-1 + x2) + x1(-3 + 3x1 + 2x2)))HR(j,tt,gen))) +
	- FR(i,tt,gen)(2muu1muu3sqrtS*2GL(j,tt,gen)(2muu1muu3sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 + 4muu32 - 2x1))*HL(j,tt,gen) +
	- (2muu32sqrtS*2(2p1k1 - 2p1k2 + sqrtS*2) + 2p1k1(-4p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(2p1k2x1 - sqrtS*2(-1 + x1 - 2muu12(-1 + x2) + x2)))*HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*2(p1k1(-8p1k2 + 2sqrtS2x2) + 2muu32sqrtS*2
	- (2p1k1 - 2p1k2 + sqrtS*2(1 - x1 + x2)) +
	- sqrtS*2(p1k2(-2 + 4x1) + sqrtS*2(1 + 2muu12(-1 + x1) - x1 - x1x2)))HL(j,tt,gen) +
	- (-32p1k12p1k2 - 8p1k1p1k2sqrtS2 + sqrtS6 + 8muu3*4sqrtS**6 +
	- 24p1k1p1k2sqrtS2x1 + 2p1k2sqrtS*4x1 - 2sqrtS6x1 -
	- 2p1k2sqrtS*4x12 + sqrtS6x12 + 2p1k1sqrtS4x2 - sqrtS*6x2 -
	- 2p1k1sqrtS*4x1x2 + sqrtS6x1*x2 +
	- 2muu32sqrtS*2(8p1k12 + p1k1(-8p1k2 + 2sqrtS*2(-2*x1 + x2)) -
	- sqrtS*2(-2p1k2(-2 + x1) + sqrtS*2(-2 + 2*x1 + x2))) -
	- 2muu12sqrtS*4(8muu34sqrtS*2 + 2p1k2x1 + sqrtS2(-2 + 3x1)(-1 + x1 + x2) -
	- 2p1k1(-2 + 2*x1 + x2) +
	- muu3*2(4p1k1 - 4p1k2 - 2sqrtS2(-3 + 4x1 + 2x2))))*HR(j,tt,gen))))))/
	- (sqrtS*2(-1 + x1)(2p1k2 - (muu12 + muu32)sqrtS2 + MSf(i,tt,gen)2)
	- (-2p1k1 - sqrtS2(muu12 + muu32 - x1) + MSf(i,tt,gen)*2)
	- (-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2))


	** TUchannel bremsstrahlung 33

	radTU = radTU - ((CFgs2NCsqrtS2(-2 + 4muu12 - 4muu3*2 + x1 + x2)
	- (ELc(i,tt,gen)(2FR(i,tt,gen)(muu1sqrtSGR(j,tt,gen)(muu1sqrtS3
	- (1 + 2muu32 - x1 - x2)HL(j,tt,gen) - 2muu3p1k2sqrtSHR(j,tt,gen)) +
	- muu3sqrtSGL(j,tt,gen)(muu1sqrtS(2p1k1 - sqrtS*2x1)*HL(j,tt,gen) +
	- (-1 + 2muu12)muu3sqrtS3HR(j,tt,gen))) +
	- FL(i,tt,gen)(-4muu1muu3sqrtS*2GR(j,tt,gen)(p1k1HL(j,tt,gen) + muu1muu3sqrtS*2HR(j,tt,gen)) +
	- GL(j,tt,gen)((-2p1k2sqrtS2x1 + (-1 + 2muu12)sqrtS*4(1 + 2muu3*2 - x1 - x2) +
	- p1k1(8p1k2 - 2sqrtS2x2))*HL(j,tt,gen) -
	- 2muu1muu3sqrtS2(-2p1k2 + sqrtS2x2)*HR(j,tt,gen)))) -
	- ER(i,tt,gen)(2FL(i,tt,gen)(muu3sqrtSGR(j,tt,gen)
	- ((1 - 2muu12)muu3sqrtS3HL(j,tt,gen) + muu1sqrtS(-2p1k1 + sqrtS2x1)*HR(j,tt,gen)) +
	- muu1sqrtSGL(j,tt,gen)(2muu3p1k2sqrtSHL(j,tt,gen) + muu1sqrtS*3(-1 - 2muu32 + x1 + x2)HR(j,tt,gen))) +
	- FR(i,tt,gen)(4muu1muu3sqrtS*2GL(j,tt,gen)(muu1muu3sqrtS2HL(j,tt,gen) + p1k1*HR(j,tt,gen)) +
	- GR(j,tt,gen)(2muu1muu3sqrtS*2(-2p1k2 + sqrtS2x2)*HL(j,tt,gen) +
	- (2p1k2sqrtS*2x1 - (-1 + 2muu12)sqrtS*4(1 + 2muu3*2 - x1 - x2) +
	- p1k1(-8p1k2 + 2sqrtS2x2))*HR(j,tt,gen))))))/
	- ((-2p1k2 + (muu12 + muu32)sqrtS2 - MSf(i,tt,gen)2)(2p1k1 + sqrtS*2(muu12 + muu32 - x1) - MSf(i,tt,gen)*2)
	- (2p1k2 + sqrtS2(muu12 + muu32 - x2) - MSf(j,tt,gen)*2)(-2p1k1 + (muu12 + muu32)sqrtS2 - MSf(j,tt,gen)2)))


	** TUchannel dipole 312

	radTU = radTU - (16CFGS*2NC(-4muu32 + x22)((-2 + 4muu3**2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x22))(4muu34 + muu32(4 - 6x2) + (-1 + x2)(-1 + x1 + 2x2)))/
	- ((-1 + 2muu32)(-1 + x2)))*
	- (-((-1 + 2muu32)sqrtS*4(ELc(i,tt,gen)((-1 + 2muu1*2)FL(i,tt,gen)*GL(j,tt,gen) +
	- 2muu12FR(i,tt,gen)GR(j,tt,gen))HL(j,tt,gen) +
	- ER(i,tt,gen)(2muu1*2FL(i,tt,gen)GL(j,tt,gen) + (-1 + 2muu1*2)FR(i,tt,gen)GR(j,tt,gen))HR(j,tt,gen)))/2. +
	- muu3*2sqrtS*4(ER(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen) +
	- 2muu12(FR(i,tt,gen)GL(j,tt,gen) - FL(i,tt,gen)GR(j,tt,gen))(ER(i,tt,gen)HL(j,tt,gen) - ELc(i,tt,gen)*HR(j,tt,gen))) +
	- (muu1muu3sqrtS*2((-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22)))
	- (ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen)) +
	- (4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) + dsqrt(-4muu32 + x22)))
	- (ER(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen))) +
	- (4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) + dsqrt(-4muu32 + x22)))
	- (muu1muu3sqrtS*2(ELc(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(-(dsqrt(1 - 4muu32)x2) + dsqrt(-4muu32 + x22)))
	- (ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)))/
	- (4.dsqrt(-4muu32 + x22))) +
	- (-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22)))
	- (muu1muu3sqrtS*2(ELc(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + dsqrt(-4muu32 + x22)))
	- (ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)))/
	- (4.dsqrt(-4muu32 + x22))))
	- /(2.dsqrt(-4muu32 + x22))))/
	- (sqrtS*2(-1 + x2)(-4dsqrt(1 - 4muu32)p1k2 + sqrtS*2(dsqrt(1 - 4muu32)x2 + (-1 + 2muu12 + 2muu3*2)
	- dsqrt(-4muu32 + x2*2)) -
	- 2dsqrt(-4muu32 + x22)MSf(i,tt,gen)2)(-4dsqrt(1 - 4muu3*2)p1k2 +
	- sqrtS*2(dsqrt(1 - 4muu32)x2 + (1 - 2muu12 - 2muu3*2)dsqrt(-4muu32 + x2*2)) +
	- 2dsqrt(-4muu32 + x22)MSf(j,tt,gen)*2))


	** TUchannel dipole 321

	radTU = radTU - (16CFGS*2NC(-4muu32 + x12)((-2 + 4muu3**2)/(-2 + x1 + x2) +
	- (dsqrt((-1 + 4muu32)/(4muu32 - x12))(4muu34 + muu32(4 - 6x1) + (-1 + x1)(-1 + 2x1 + x2)))/
	- ((-1 + 2muu32)(-1 + x1)))*
	- (-((-1 + 2muu32)sqrtS*4(ELc(i,tt,gen)((-1 + 2muu1*2)FL(i,tt,gen)*GL(j,tt,gen) +
	- 2muu12FR(i,tt,gen)GR(j,tt,gen))HL(j,tt,gen) +
	- ER(i,tt,gen)(2muu1*2FL(i,tt,gen)GL(j,tt,gen) + (-1 + 2muu1*2)FR(i,tt,gen)GR(j,tt,gen))HR(j,tt,gen)))/2. +
	- muu3*2sqrtS*4(ER(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen) +
	- 2muu12(FR(i,tt,gen)GL(j,tt,gen) - FL(i,tt,gen)GR(j,tt,gen))(ER(i,tt,gen)HL(j,tt,gen) - ELc(i,tt,gen)*HR(j,tt,gen))) +
	- (muu1muu3sqrtS*2((4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + dsqrt(-4muu32 + x12)))
	- (ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen)) +
	- (-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12)))
	- (ER(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ELc(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen))) +
	- (4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + dsqrt(-4muu32 + x12)))
	- (muu1muu3sqrtS*2(ELc(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HR(j,tt,gen)) +
	- ((4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(-(dsqrt(1 - 4muu32)x1) + dsqrt(-4muu32 + x12)))
	- (ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)))/
	- (4.dsqrt(-4muu32 + x12))) +
	- (-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12)))
	- (muu1muu3sqrtS*2(ELc(i,tt,gen)FR(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FL(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)) +
	- ((-4dsqrt(1 - 4muu3*2)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + dsqrt(-4muu32 + x12)))
	- (ELc(i,tt,gen)FL(i,tt,gen)GL(j,tt,gen)HL(j,tt,gen) + ER(i,tt,gen)FR(i,tt,gen)GR(j,tt,gen)HR(j,tt,gen)))/
	- (4.dsqrt(-4muu32 + x12))))
	- /(2.dsqrt(-4muu32 + x12))))/
	- (sqrtS*2(-1 + x1)(-4dsqrt(1 - 4muu32)p1k1 + sqrtS*2(dsqrt(1 - 4muu32)x1 + (1 - 2muu12 - 2muu3*2)
	- dsqrt(-4muu32 + x1*2)) +
	- 2dsqrt(-4muu32 + x12)MSf(i,tt,gen)2)(-4dsqrt(1 - 4muu3*2)p1k1 +
	- sqrtS*2(dsqrt(1 - 4muu32)x1 + (-1 + 2muu12 + 2muu3*2)dsqrt(-4muu32 + x1*2)) -
	- 2dsqrt(-4muu32 + x12)MSf(j,tt,gen)*2))

	enddo
	enddo

	end

File Metadata

Mime Type: text/x-c
Expires: Wed, May 14, 11:13 AM (20 h, 31 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 5111377
Default Alt Text: ChiChi2QQ_M2Radiation.F (162 KB)

ChiChi2QQ_M2Radiation.FNo OneTemporaryActions

ChiChi2QQ_M2Radiation.FView Options

File Metadata

Event Timeline

ChiChi2QQ_M2Radiation.F
No OneTemporary
Actions

ChiChi2QQ_M2Radiation.F
View Options