Page MenuHomeHEPForge
Files F10275373

d217h0l1.f90
No OneTemporary
Actions

d217h0l1.f90
View Options

module uussbb_d217h0l1
! file: /home/thomasr/Work/samurai-processes/uussbb/helicity0/d217h0l1.f
! 90
! generator: haggies (1.1)
use precision, only: ki
use uussbb_config
use uussbb_model
use uussbb_kinematics
use uussbb_util, only: cond
use uussbb_color
use uussbb_abbrevh0l1
implicit none
private :: brack_1
private :: brack_2
private :: brack_3
private :: brack_4
private :: brack_5
complex(ki), parameter :: i_ = (0.0_ki, 1.0_ki)
private :: i_
contains
pure function brack_1(Q, mu2) result(brack)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
complex(ki) :: brack
brack = (i_/((es12+es34-(es56+es345))*es34))
end function brack_1
pure function brack_2(Q, mu2) result(brack)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
complex(ki), dimension(numcs) :: brack
real(ki) :: t1
t1 = NC*NC
brack = (TR*TR*TR*(((1.0_ki)-(2.0_ki)/t1)*c4+((2.0_ki)/NC-NC)*c3+c1/(t1*NC&
&)-c2/t1))
end function brack_2
pure function brack_3(Q, mu2) result(brack)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
complex(ki) :: brack
complex(ki) :: t1
complex(ki) :: t2
complex(ki) :: t3
t1 = dotproduct(Q, spvak6k5)
t2 = dotproduct(Q, spvak4k5)
t3 = dotproduct(Q, spvak1k2)
brack = ((4.0_ki)*(mu2*abb217n4+abb217n4*es234-(t2*t3*abb217n5+t1*t3*abb21&
&7n6+abb217n4*es61+abb217n4*es56+t3*abb217n3+t2*abb217n2+t1*abb217n1)))
end function brack_3
pure function brack_4(Q, mu2) result(brack)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
complex(ki) :: brack
complex(ki) :: t1
t1 = dotproduct(Q, spvak1k2)
brack = ((2.0_ki)*(t1*abb217n5*dotproduct(Q, spvak4k5)+t1*abb217n6*dotprod&
&uct(Q, spvak6k5)))
end function brack_4
pure function brack_5(Q, mu2) result(brack)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
integer :: brack
brack = 0
end function brack_5
pure subroutine diagram217(Q, mu2, epspow, res)
implicit none
complex(ki), dimension(4), intent(in) :: Q
real(ki), intent(in) :: mu2
integer, intent(in) :: epspow
complex(ki), dimension(1:numcs), intent(inout) :: res
complex(ki), dimension(1:numcs) :: acc
complex(ki) :: prefactor
complex(ki), dimension(numcs) :: cf1
! res is set to zero in the calling routines,
! therefore we sum to whatever is there already.
acc = 0.0_ki
prefactor = brack_1(Q, mu2)
cf1 = brack_2(Q, mu2)
acc(:) = acc(:) + cf1(:) * ((cond(epspow.eq.0,brack_3,Q,mu2)+cond(epspow.e&
&q.1,brack_4,Q,mu2)))
! d217: if non-zero, there is only one color structure
acc(1) = acc(1) + (brack_5(Q, mu2))
res = res + prefactor * acc
end subroutine diagram217
end module uussbb_d217h0l1

File Metadata

Mime Type
text/plain
Expires
Fri, Apr 4, 9:29 PM (13 m, 20 s)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4737386
Default Alt Text
d217h0l1.f90 (3 KB)

Event Timeline