Page MenuHomeHEPForge
Files F10881134

tsums_omega.f
No OneTemporary
Actions

tsums_omega.f
View Options

!-----------------------
! madgraph - a Feynman Diagram package by Tim Stelzer and Bill Long
! (c) 1993
!
! Filename: tsums.f
!-----------------------
Subroutine ttsum(color,iterm,ifact1,element1,ifact2,element2)
!**************************************************************************
! Does the sum over two gluons in T matrices
!**************************************************************************
implicit none
! Constants
integer nc
parameter(nc=3)
integer maxlines, maxfactors, maxterms
parameter(maxlines=8,maxfactors=9,maxterms=250)
! Arguments
integer color(0:maxlines,0:maxfactors,0:maxterms)
integer ifact1,ifact2,element1,element2,iterm
! Local Variables
integer nterms,nfactors
integer ielement,nadded
!-----------
! Begin Code
!-----------
! Add First Term
call copyterm(color,iterm,ifact1,ifact2) !first make a copy
nterms=color(0,0,0)
! First check to see if this is a quark line or a gluon line
if (element1 .le. 2) then !this is a quark line
if (element2 .gt. 2) then
print*,'Error one quark and one glue in ttsum',
& element1,element2
endif
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
if (element1 .ne. 1 .or. element2 .ne. 2) then
print*,'Element 1 is wrong',element1,element2
endif
color(2,nfactors,nterms)=color(2,ifact1,iterm)
color(1,nfactors,nterms)=color(1,ifact2,iterm)
nadded=2
do ielement=3,-color(0,ifact2,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
do ielement=3,-color(0,ifact1,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
else
if (element2 .lt. 2) then
print*,'Error one quark and one glue in ttsum',
& element1,element2
stop
endif
! Add First Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact1,iterm)
color(2,nfactors,nterms)=color(2,ifact2,iterm)
nadded=2
do ielement=3,element1-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
do ielement=element2+1,-color(0,ifact2,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
! Add Second Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact2,iterm)
color(2,nfactors,nterms)=color(2,ifact1,iterm)
nadded=2
do ielement=3,element2-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
do ielement=element1+1,-color(0,ifact1,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
color(2,1,nterms) = color(2,1,nterms)*2
! Add Second Term
call copyterm(color,iterm,ifact1,ifact2) !first make a copy
nterms=color(0,0,0)
! Add First Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact1,iterm)
color(2,nfactors,nterms)=color(2,ifact1,iterm)
nadded=2
do ielement=3,element1-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
do ielement=element1+1,-color(0,ifact1,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
! Add Second Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact2,iterm)
color(2,nfactors,nterms)=color(2,ifact2,iterm)
nadded=2
do ielement=3,element2-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
do ielement=element2+1,-color(0,ifact2,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
color(2,1,nterms) = -color(2,1,nterms)*2*Nc
endif
color(1,1,iterm) = 0 !This means clear it
call ClearZeros(color)
end
Subroutine tsum(color,iterm,ifact1,element1,ifact2,element2)
!**************************************************************************
! Does the sum over two gluons in T matrices
!**************************************************************************
implicit none
! Constants
integer nc
parameter(nc=3)
integer maxlines, maxfactors, maxterms
parameter(maxlines=8,maxfactors=9,maxterms=250)
! Arguments
integer color(0:maxlines,0:maxfactors,0:maxterms)
integer ifact1,ifact2,element1,element2,iterm
! Local Variables
integer nterms,nfactors
integer ielement,nadded
!-----------
! Begin Code
!-----------
if (ifact1 .ne. ifact2) then
print*,'Tsum called when factors not right',ifact1,ifact2
stop
endif
if (element1 .le. 2 .and. element2 .le. 2) then
call SumQuarks(color)
else
! Add First Term
call copyterm(color,iterm,ifact1,ifact2) !first make a copy
nterms=color(0,0,0)
! Add First Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact1,iterm)
color(2,nfactors,nterms)=color(2,ifact2,iterm)
nadded=2
do ielement=3,element1-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
enddo
do ielement=element2+1,-color(0,ifact2,iterm) !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
color(0,nfactors,nterms) = -nadded !number of elements
! Add Second Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
nadded=0
do ielement=element1+1,element2-1 !add remaining glue's
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact2,iterm)
enddo
color(0,nfactors,nterms) = nadded !number of elements
color(2,1,nterms) = color(2,1,nterms)*2
! Add Second Term
call copyterm(color,iterm,ifact1,ifact2) !first make a copy
nterms=color(0,0,0)
! Add First Factor in Term
color(0,0,nterms)=color(0,0,nterms)+1 !add one factor
nfactors=color(0,0,nterms)
color(1,nfactors,nterms)=color(1,ifact1,iterm)
color(2,nfactors,nterms)=color(2,ifact1,iterm)
nadded=2
do ielement=3,abs(color(0,ifact1,iterm)) !add remaining glue's
if (ielement .ne. element1 .and.
& ielement .ne. element2) then
nadded=nadded+1
color(nadded,nfactors,nterms)=color(ielement,ifact1,iterm)
endif
enddo
color(0,nfactors,nterms) = -nadded !number of elements
! Add Second Factor in Term
color(2,1,nterms) = -color(2,1,nterms)*2*Nc
color(1,1,iterm) = 0 !this means clear it
endif
call ClearZeros(color)
end

File Metadata

Mime Type
text/plain
Expires
Sat, May 3, 5:55 AM (17 h, 56 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4982867
Default Alt Text
tsums_omega.f (8 KB)

Event Timeline