modhiggsj Module Reference

Functions/Subroutines
subroutine, public	evalamp_hj (p, res)
subroutine	me2_ggg_tree (j1, j2, j3, sprod, me2)
subroutine	me2_qbqg_tree (j1, j2, j3, sprod, me2)
subroutine	spinoru2 (n, p, za, zb, s)

Function/Subroutine Documentation

evalamp_hj()

subroutine, public modhiggsj::evalamp_hj	(	real(dp), dimension(4,4), intent(in)	p,
		real(dp), dimension(-5:5,-5:5), intent(out)	res
	)

Definition at line 12 of file mod_HiggsJ.F90.

    real(dp), intent(in) :: p(4,4)
    complex(dp) :: heftcoupl
    real(dp), intent(out) :: res(-5:5,-5:5)
    real(dp) :: sprod(4,4)
    complex(dp) :: za(4,4), zb(4,4)
    real(dp) :: restmp!, restmpid
    integer :: i,j
 
    heftcoupl = gs*alphas/(6.0_dp * pi * vev)
        res=zero
 
    call spinoru2(3,(/-p(:,1),-p(:,2),p(:,4)/),za,zb,sprod)
 
    call me2_ggg_tree(1,2,3,sprod,res(0,0))
    res(0,0) = res(0,0)*avegg
 
    do i=1,5
      call me2_qbqg_tree(2,3,1,sprod,res(0,i))
      res(0,i) = -res(0,i)*aveqg
      !print *,res(0,i)
    enddo
 
    do i=1,5
      call me2_qbqg_tree(1,3,2,sprod,res(i,0))
      res(i,0) = -res(i,0)*aveqg
      !print *,res(i,0)
    enddo
 
    do i=1,5
      call me2_qbqg_tree(3,2,1,sprod,res(0,-i))
      res(0,-i) = -res(0,-i)*aveqg
      !print *,res(0,-i)
    enddo
 
    do i=1,5
      call me2_qbqg_tree(1,3,2,sprod,res(-i,0))
      res(-i,0) = -res(-i,0)*aveqg
    !if(res(-i,0).le.0d0) then
    !print *,res(-i,0)
    !endif
    enddo
 
    do i=1,5
      call me2_qbqg_tree(1,2,3,sprod,res(-i,i))
      res(-i,i) = res(-i,i) * aveqq
!      res(i,-i) = res(i,-i) + restmp * aveqq
!      res(0, i) = res(0, i) + restmp * aveqg
!      res(0,-i) = res(0,-i) + restmp * aveqg
!      res(i, 0) = res(i, 0) + restmp * aveqg
!      res(-i,0) = res(-i,0) + restmp * aveqg
    enddo
 
    do i=1,5
      call me2_qbqg_tree(2,1,3,sprod,res(i,-i))
      res(i,-i) = res(i,-i)*aveqq
    enddo
 
    res = res * (heftcoupl**2)
!if(res(0,0).gt.10000d0)then
    !do i = -5,5
    !do j = -5,5
    !print *, res(i,j),i,j
    !enddo
    !enddo
!    print *, dsqrt(scr(p(:,1),p(:,1))), dsqrt(scr(p(:,2),p(:,2))), dsqrt(scr(p(:,3),p(:,3))), dsqrt(scr(p(:,4),p(:,4)))
!    print*, p
!endif
 

me2_ggg_tree()

subroutine modhiggsj::me2_ggg_tree ( integer, intent(in)  j1, integer, intent(in)  j2, integer, intent(in)  j3, real(dp), dimension(3,3), intent(in)  sprod, real(dp), intent(out)  me2  )

private

Definition at line 87 of file mod_HiggsJ.F90.

   integer, intent(in) :: j1,j2,j3
   real(dp), intent(in) :: sprod(3,3)
   real(dp), intent(out) :: me2
   real(dp) :: s12, s13, s23, qsq
   real(dp), parameter :: col_gg = 4.0d0 * ca**2 * cf
 
   me2 = zero
 
   s12 = sprod(j1,j2)
   s13 = sprod(j1,j3)
   s23 = sprod(j2,j3)
 
   qsq = s12 + s13 + s23
 
   me2 = (s12**4 + s13**4 + s23**4 + qsq**4)/s12/s13/s23
 
   me2 = me2 * two * col_gg
 
   return
 

me2_qbqg_tree()

subroutine modhiggsj::me2_qbqg_tree ( integer, intent(in)  j1, integer, intent(in)  j2, integer, intent(in)  j3, real(dp), dimension(3,3), intent(in)  sprod, real(dp), intent(out)  me2  )

private

Definition at line 112 of file mod_HiggsJ.F90.

   integer, intent(in) :: j1, j2, j3
   real(dp), intent(in) :: sprod(3,3)
   real(dp), intent(out) :: me2
   real(dp) :: s12, s13,s23
   real(dp), parameter :: col_qg = two * xn * cf
 
   me2 = zero
 
   s12 = sprod(j1,j2)
   s13 = sprod(j1,j3)
   s23 = sprod(j2,j3)
 
   me2 = two * (s13**2 + s23**2)/s12 * col_qg
!if(me2.le.0d0)then
!print *, s12,s13,s23
!endif
   return
 

spinoru2()

subroutine modhiggsj::spinoru2 ( integer, intent(in)  n, real(dp), dimension(4,n), intent(in)  p, complex(dp), dimension(n,n), intent(out)  za, complex(dp), dimension(n,n), intent(out)  zb, real(dp), dimension(n,n), intent(out)  s  )

private

Definition at line 137 of file mod_HiggsJ.F90.

    implicit none
    integer, intent(in) :: n
    real(dp), intent(in) :: p(4,n)
    complex(dp), intent(out) :: za(n,n), zb(n,n)
    real(dp), intent(out) :: s(n,n)
    integer :: i,j
    complex(dp) :: c23(n), f(n)
    real(dp) :: rt(n)
 
    !---if one of the vectors happens to be zero this routine fails.
    do j=1,n
       za(j,j)=czero
       zb(j,j)=za(j,j)
 
       !-----positive energy case
       if (p(1,j) .gt. zero) then
          rt(j)=sqrt(p(2,j)+p(1,j))
          c23(j)=cmplx(p(4,j),-p(3,j),kind=dp)
          f(j)=(one,zero)
       else
       !-----negative energy case
          rt(j)=sqrt(-p(1,j)-p(2,j))
          c23(j)=cmplx(-p(4,j),p(3,j),kind=dp)
          f(j)=ci
       endif
    enddo
 
    do i=2,n
 
     do j=1,i-1
          s(i,j)=two*scr(p(:,i),p(:,j))
          za(i,j)=f(i)*f(j) &
               *(c23(i)*cmplx(rt(j)/rt(i),kind=dp)-c23(j)*cmplx(rt(i)/rt(j),kind=dp))
 
          if (abs(s(i,j)).lt.1d-5) then
             zb(i,j)=-(f(i)*f(j))**2*conjg(za(i,j))
          else
             zb(i,j)=-cmplx(s(i,j),kind=dp)/za(i,j)
          endif
 
          za(j,i)=-za(i,j)
          zb(j,i)=-zb(i,j)
          s(j,i)=s(i,j)
 
       enddo
 
    enddo
 
    return