anomGamwt.f

this file is to print, not to load !
        Subroutine dipolgamma(iqed, E, theta, A, B, R)
        implicit none
        real*8 amz0,gamz0,sin2w0,alphaqed,gswr(7),gswi(7)
        integer iqed,channel,k,j,k0,j0
        real*8 e,theta,a,b,xnor,thet, buf
        real*8 r(1:4, 1:4),r0(1:4, 1:4),sign(4)
!       ReA0=0
!       ImA0=0
!       ReB=0
!       ImB=0
!       ReX=0
!       ImX=0
!       ReY=0
!       ImY=0
        sign(1)=-1.0
        sign(2)= 1.0
        sign(3)=-1.0
        sign(4)= 1.0
        thet=acos(cos(theta))
        
        call dipolgammarij(iqed, e, theta, a, b, r0)
 
 
        
! rotation from q ~q --> tau- tau+  to  ~q q --> tau+ tau- frames.
!               WARNING: we normalize R0, but  x-section is set into R0(4,4).        
        xnor=r0(4,4)
        do k=1,4
          do j=1,4
              r0(k,j)=r0(k,j)/xnor*sign(k)*sign(j)
           enddo
        enddo
        r0(4,4)=xnor
 
 ! rotation in second index by angle pi/2 around z-axis       
        do k=1,4
           do j=1,4
              if(j.eq.1) then
                 buf=r0(k,2)
                 r0(k,2)=r0(k,j)
                 r0(k,1)=-buf
              endif
              
           enddo
        enddo
 
! rotation in first index by angle pi/2  around z-axis 
        do j=1,4
           do k=1,4
             if(k.eq.1) then
                 buf=r0(2,j)
                 r0(2,j)=r0(k,j)
                 r0(1,j)=-buf
              endif
 
           enddo
        enddo
 
 
! shift order of entries  from x,y,z,t  to t,x,y,z  (necessary for C++ conventions). 
        do k0=1,4
           do j0=1,4
               k=k0+1
               j=j0+1
               if(j.eq.5) j=1
               if(k.eq.5) k=1
               r(k,j)=r0(k0,j0)
           enddo
        enddo
 !       if(R0(3,3)/R0(4,4).gt.1d0) then
  !         write(*,*) "z fortranu thet=", thet," energy= ",energy, " Rzz=", R0(3,3)/R0(4,4)," ", R0(3,3)," ",R0(4,4)
           
 !          write(*,*) "iqed,Energy, thet, channel ",iqed,Energy, thet, channel
        !   write(*,*)  "Amz0,Gamz0,sin2W0,alphaQED ",   Amz0,Gamz0,sin2W0,alphaQED
      !      write(*,*)  "ReA0, ImA0, ReB, ImB ",   ReA0, ImA0, ReB, ImB
      !      write(*,*)  "ReX, ImX, ReY, ImY ",    ReX, ImX, ReY, ImY
   !      write(*,*)  "GSWr ",   GSWr
   !      write(*,*)  "GSWi ",   GSWi
   !     endif
    !  IF(R(4,4)/R(1,1).gt.1.0) then
    !    write(*,*) "==============",R(4,4)/R(1,1)
    !    write(*,*) 'energuy,theta,channel= ',Energy,' ',cos(theta),' ',channel
    !    do k=1,4
    !       write(*,7) k, R(k,1)/R(1,1),R(k,2)/R(1,1),R(k,3)/R(1,1),R(k,4)/R(1,1)
    !   enddo
    !  endif
      
 
 7      format(1x,i4, 4(2x,f12.6))
        
      end