!--------------------------------------------------------------- ! This subroutine checks the number of external dipoles---- !---absorbing the correct number into the fragmenation functions ! it then returns the finite (msq_qcd*dip) --------------------- !--------------------------------------------------------------- !--- Author C. Williams Nov 2012 !----------------------------------------------------------------- ! modified under frag re-write March 2013 subroutine qqb_dm_monophot_fragdips(p,p_phys,qcd_tree,msq_out) implicit none include 'constants.f' include 'ewcouple.f' include 'ewcharge.f' include 'frag.f' double precision p(mxpart,4),p_phys(mxpart,4) double precision msq_qcd(-nf:nf,-nf:nf),msq_out(-nf:nf,-nf:nf) integer j,k double precision virt_dips,xl,dot,fsq double precision aewo2pi,fi_gaq external qcd_tree aewo2pi=esq/(fourpi*twopi) fsq=frag_scale**2 xl=dlog(-two*dot(p_phys,2,5)/fsq) virt_dips=+aewo2pi*(fi_gaq(z_frag,p_phys,xl,5,2,2)) do j=-nf,nf do k=-nf,nf msq_qcd(j,k)=0d0 msq_out(j,k)=0d0 enddo enddo call qcd_tree(p,msq_qcd) do j=-nf,nf do k=-nf,nf if((j.eq.0).and.(k.ne.0)) then msq_out(j,k)=Q(k)**2*msq_qcd(j,k)*virt_dips elseif((j.ne.0).and.(k.eq.0)) then msq_out(j,k)=Q(j)**2*msq_qcd(j,k)*virt_dips endif enddo enddo return end subroutine