Hi Rik,

Negate all the cut automatically is easy. I can implement that if you want to.
I do that already at LO.

Cheers,

Olivier
On 08 Oct 2015, at 19:24, Rikkert Frederix <email address hidden> wrote:

> Rikkert Frederix has proposed merging lp:~maddevelopers/mg5amcnlo/NLO_decay_process into lp:~maddevelopers/mg5amcnlo/2.3.3.
> 
> Requested reviews:
>  Hua-Sheng Shao (erdissshaw)
>  marco zaro (marco-zaro)
> 
> For more details, see:
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/NLO_decay_process/+merge/273880
> 
> Allow for the computation of (partial) decay widths at NLO accuracy, e.g.,
> 
> generate Z > b b~ [QCD]
> 
> or 
> 
> import model loop_sm-no_b_mass
> generate t > e+ ve j [QCD]
> 
> etc. 
> 
> The renormalisation (and Ellis-Sexton) scales are automatically set to the mass of the decaying particle. Cuts are not automatically set to zero, so, in general you have to negate them in the run_card by hand. 
> -- 
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/mg5amcnlo/2.3.3.
> === modified file 'Template/NLO/SubProcesses/driver_mintFO.f'
> --- Template/NLO/SubProcesses/driver_mintFO.f	2015-10-02 07:00:09 +0000
> +++ Template/NLO/SubProcesses/driver_mintFO.f	2015-10-08 17:23:54 +0000
> @@ -170,7 +170,7 @@
>       else
>         flat_grid=.false.
>       endif
> -      ndim = 3*(nexternal-2)-4
> +      ndim = 3*(nexternal-nincoming)-4
>       if (abs(lpp(1)) .ge. 1) ndim=ndim+1
>       if (abs(lpp(2)) .ge. 1) ndim=ndim+1
> c Don't proceed if muF1#muF2 (we need to work out the relevant formulae
> @@ -709,9 +709,14 @@
>             endif
>          elseif(buffer(1:12).eq.'SUM_HELICITY') then
>             read(buffer(15:),*) i
> -            if (i.eq.0) then
> -               mc_hel=0
> -               write (*,*) 'Explicitly summing over helicities'
> +            if (nincoming.eq.1) then
> +               write (*,*) 'Sum over helicities in the virtuals'/
> +     $              /' for decay process'
> +               mc_hel=0
> +            elseif (i.eq.0) then
> +               mc_hel=0
> +               write (*,*) 'Explicitly summing over helicities'/
> +     $              /' for the virtuals'
>             else
>                mc_hel=1
>                write(*,*) 'Do MC over helicities for the virtuals'
> 
> === modified file 'Template/NLO/SubProcesses/driver_mintMC.f'
> --- Template/NLO/SubProcesses/driver_mintMC.f	2015-09-25 14:53:36 +0000
> +++ Template/NLO/SubProcesses/driver_mintMC.f	2015-10-08 17:23:54 +0000
> @@ -108,6 +108,11 @@
>       call cpu_time(tBefore)
>       fixed_order=.false.
>       nlo_ps=.true.
> +      if (nincoming.ne.2) then
> +         write (*,*) 'Decay processes not supported for'/
> +     &        /' event generation'
> +         stop 1
> +      endif
> 
> c     Read general MadFKS parameters
> c
> @@ -158,7 +163,7 @@
>       else
>         flat_grid=.false.
>       endif
> -      ndim = 3*(nexternal-2)-4
> +      ndim = 3*(nexternal-nincoming)-4
>       if (abs(lpp(1)) .ge. 1) ndim=ndim+1
>       if (abs(lpp(2)) .ge. 1) ndim=ndim+1
> c Don''t proceed if muF1#muF2 (we need to work out the relevant formulae
> @@ -638,12 +643,17 @@
> 
>       write(*,10) 'Exact helicity sum (0 yes, n = number/event)? '
>       read(*,*) i
> -      if (i .eq. 0) then
> -         mc_hel= 0
> -         write(*,*) 'Explicitly summing over helicities for virt'
> +      if (nincoming.eq.1) then
> +         write (*,*) 'Sum over helicities in the virtuals'/
> +     $        /' for decay process'
> +         mc_hel=0
> +      elseif (i.eq.0) then
> +         mc_hel=0
> +         write (*,*) 'Explicitly summing over helicities'/
> +     $        /' for the virtuals'
>       else
> -         mc_hel= i
> -         write(*,*) 'Summing over',i,' helicities/event for virt'
> +         mc_hel=1
> +         write(*,*) 'Do MC over helicities for the virtuals'
>       endif
>       isum_hel = 0
> 
> 
> === modified file 'Template/NLO/SubProcesses/fks_Sij.f'
> --- Template/NLO/SubProcesses/fks_Sij.f	2012-08-28 21:06:34 +0000
> +++ Template/NLO/SubProcesses/fks_Sij.f	2015-10-08 17:23:54 +0000
> @@ -90,9 +90,13 @@
> 
> c Consistency check -- call to set_cms_stuff() must be done prior to
> c entering this function
> -      shattmp=2d0*dot(p(0,1),p(0,2))
> +      if (nincoming.eq.2) then
> +         shattmp=2d0*dot(p(0,1),p(0,2))
> +      else
> +         shattmp=p(0,1)**2
> +      endif
>       if(abs(shattmp/shat-1.d0).gt.1.d-5)then
> -        write(*,*)'Error in fks_Sij: inconsistent shat'
> +        write(*,*)'Error in fks_Sij: inconsistent shat #1'
>         write(*,*)shattmp,shat
>         stop
>       endif
> @@ -148,6 +152,7 @@
>            do j=1,fks_j_from_i(i,0)
>             kk = i
>             ll = fks_j_from_i(i,j)
> +            if (nincoming.ne.2 .and. ll.le.nincoming) cycle
>             if     ( ijskip(kk,ll).eq.0 .and. ijskip(ll,kk).eq.0 ) then
>                ijskip(kk,ll) = 1
>             elseif ( ijskip(kk,ll).eq.0 .and. ijskip(ll,kk).eq.1 ) then
> @@ -178,6 +183,7 @@
>         do j=1,fks_j_from_i(i,0)
>          kk = i
>          ll = fks_j_from_i(i,j)
> +         if (nincoming.ne.2 .and. ll.le.nincoming) cycle
>          if(ijskip(kk,ll).ne.1)goto 222
>          if(particle_type(ll).eq.8.and.particle_type(kk).ne.8.and.
>      #      ll.gt.nincoming)then
> @@ -567,7 +573,11 @@
> 
> c Consistency check -- call to set_cms_stuff() must be done prior to
> c entering this function
> -      shattmp=2d0*dot(p(0,1),p(0,2))
> +      if (nincoming.eq.2) then
> +         shattmp=2d0*dot(p(0,1),p(0,2))
> +      else
> +         shattmp=p(0,1)**2
> +      endif
>       if(abs(shattmp/shat-1.d0).gt.1.d-5)then
>         write(*,*)'Error in fks_Hij: inconsistent shat'
>         write(*,*)shattmp,shat
> 
> === modified file 'Template/NLO/SubProcesses/fks_singular.f'
> --- Template/NLO/SubProcesses/fks_singular.f	2015-09-14 07:31:51 +0000
> +++ Template/NLO/SubProcesses/fks_singular.f	2015-10-08 17:23:54 +0000
> @@ -1230,7 +1230,12 @@
> c PDG codes
>       niproc(ict)=iproc
>       do j=1,iproc
> -         parton_iproc(j,ict)=pd(j)*conv
> +         if (nincoming.eq.2) then
> +            parton_iproc(j,ict)=pd(j)*conv
> +         else
> +c           Keep GeV's for decay processes (no conv. factor needed)
> +            parton_iproc(j,ict)=pd(j)
> +         endif
>          do k=1,nexternal
>             parton_pdg(k,j,ict)=idup_d(iFKS,k,j)
>             if (k.lt.fks_j_d(iFKS)) then
> @@ -2199,9 +2204,15 @@
>             do ii=1,iproc_save(nFKS(ict))
>                if (eto(ii,nFKS(ict)).ne.ipr) cycle
>                n_ctr_found=n_ctr_found+1
> -               write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> -     &              (wgt(j,ict)*conv,j=1,3),
> -     &              nexternal
> +               if (nincoming.eq.2) then
> +                  write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> +     &                 (wgt(j,ict)*conv,j=1,3),
> +     &                 nexternal
> +               else
> +                  write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> +     &                 (wgt(j,ict),j=1,3),
> +     &                 nexternal
> +               endif
>                procid=''
>                do j=1,nexternal
>                   write (str_temp,*) parton_pdg(j,ii,ict)
> @@ -2227,9 +2238,15 @@
> c H-event
>             ipr=iproc_picked
>             n_ctr_found=n_ctr_found+1
> -            write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> -     &           (wgt(j,ict)*conv,j=1,3),
> -     &           nexternal
> +            if (nincoming.eq.2) then
> +               write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> +     &              (wgt(j,ict)*conv,j=1,3),
> +     &              nexternal
> +            else
> +               write (n_ctr_str(n_ctr_found),'(3(1x,d18.12),1x,i2)')
> +     &              (wgt(j,ict),j=1,3),
> +     &              nexternal
> +            endif
>             procid=''
>             do j=1,nexternal
>                write (str_temp,*) parton_pdg(j,ipr,ict)
> @@ -2342,7 +2359,7 @@
>       do i=0,3
>         xsum(i)=0.d0
>         xsuma(i)=0.d0
> -        do j=3,npart
> +        do j=nincoming+1,npart
>           xsum(i)=xsum(i)+xmom(i,j)
>           xsuma(i)=xsuma(i)+abs(xmom(i,j))
>         enddo
> @@ -2409,9 +2426,14 @@
>       pass=.true.
>       jflag=0
>       do i=0,3
> -        xsum(i)=-xmom(i,1)-xmom(i,2)
> -        xsuma(i)=abs(xmom(i,1))+abs(xmom(i,2))
> -        do j=3,npart
> +        if (nincoming.eq.2) then
> +          xsum(i)=-xmom(i,1)-xmom(i,2)
> +          xsuma(i)=abs(xmom(i,1))+abs(xmom(i,2))
> +        elseif(nincoming.eq.1) then
> +          xsum(i)=-xmom(i,1)
> +          xsuma(i)=abs(xmom(i,1))
> +        endif
> +        do j=nincoming+1,npart
>           xsum(i)=xsum(i)+xmom(i,j)
>           xsuma(i)=xsuma(i)+abs(xmom(i,j))
>         enddo
> @@ -2424,14 +2446,14 @@
>           write(*,*)'Momentum is not conserved [nocms]'
>           write(*,*)'i=',i
>           do j=1,npart
> -            write(*,'(4(d14.8,1x))') (xmom(jj,j),jj=0,3)
> +            write(*,'(i2,1x,4(d14.8,1x))') j,(xmom(jj,j),jj=0,3)
>           enddo
>           jflag=1
>         endif
>       enddo
>       if(jflag.eq.1)then
> -        write(*,'(4(d14.8,1x))') (xsum(jj),jj=0,3)
> -        write(*,'(4(d14.8,1x))') (xrat(jj),jj=0,3)
> +        write(*,'(a3,1x,4(d14.8,1x))') 'sum',(xsum(jj),jj=0,3)
> +        write(*,'(a3,1x,4(d14.8,1x))') 'rat',(xrat(jj),jj=0,3)
>         pass=.false.
>         return
>       endif
> @@ -2458,7 +2480,11 @@
>         endif
>       enddo
> c
> -      ecmtmp=sqrt(2d0*dot(xmom(0,1),xmom(0,2)))
> +      if (nincoming.eq.2) then
> +         ecmtmp=sqrt(2d0*dot(xmom(0,1),xmom(0,2)))
> +      elseif (nincoming.eq.1) then
> +         ecmtmp=xmom(0,1)
> +      endif
>       if(abs(ecm-ecmtmp).gt.vtiny)then
>         write(*,*)'Inconsistent shat [nocms]'
>         write(*,*)'ecm given=   ',ecm
> @@ -2716,7 +2742,11 @@
> 
> c Consistency check -- call to set_cms_stuff() must be done prior to
> c entering this function
> -      shattmp=2d0*dot(pp(0,1),pp(0,2))
> +      if (nincoming.eq.2) then
> +         shattmp=2d0*dot(pp(0,1),pp(0,2))
> +      else
> +         shattmp=pp(0,1)**2
> +      endif
>       if(abs(shattmp/shat-1.d0).gt.1.d-5)then
>         write(*,*)'Error in sreal: inconsistent shat'
>         write(*,*)shattmp,shat
> @@ -2724,9 +2754,9 @@
>       endif
> 
>       if (1d0-y_ij_fks.lt.tiny)then
> -         if (pmass(j_fks).eq.zero.and.j_fks.le.2)then
> +         if (pmass(j_fks).eq.zero.and.j_fks.le.nincoming)then
>             call sborncol_isr(pp,xi_i_fks,y_ij_fks,wgt)
> -         elseif (pmass(j_fks).eq.zero.and.j_fks.ge.3)then
> +         elseif (pmass(j_fks).eq.zero.and.j_fks.ge.nincoming+1)then
>             call sborncol_fsr(pp,xi_i_fks,y_ij_fks,wgt)
>          else
>             wgt=0d0
> @@ -2830,7 +2860,7 @@
>       E_i_fks = p(0,i_fks)
>       z = 1d0 - E_i_fks/(E_i_fks+E_j_fks)
>       t = z * shat/4d0
> -      if(rotategranny .and. nexternal-1.ne.3)then
> +      if(rotategranny .and. nexternal-1.ne.3 .and. nincoming.eq.2)then
> c Exclude 2->1 (at the Born level) processes: matrix elements are
> c independent of the PS point, but non-zero helicity configurations
> c might flip when rotating the momenta.
> @@ -2950,7 +2980,7 @@
> c Thus, an extra factor z (implicit in the flux of the reduced Born
> c in FKS) has to be inserted here
>       t = z*shat/4d0
> -      if(j_fks.eq.2 .and. nexternal-1.ne.3)then
> +      if(j_fks.eq.2 .and. nexternal-1.ne.3 .and. nincoming.eq.2)then
> c Rotation according to innerpin.m. Use rotate_invar() if a more 
> c general rotation is needed.
> c Exclude 2->1 (at the Born level) processes: matrix elements are
> @@ -2981,7 +3011,7 @@
>                pi(i)=p_i_fks_ev(i)
>                pj(i)=p(i,j_fks)
>             enddo
> -            if(j_fks.eq.2)then
> +            if(j_fks.eq.2 .and. nincoming.eq.2)then
> c Rotation according to innerpin.m. Use rotate_invar() if a more 
> c general rotation is needed
>                pi(1)=-pi(1)
> @@ -3002,7 +3032,7 @@
>             azifact=Wij_angle/Wij_recta
>          endif
> c Insert the extra factor due to Madgraph convention for polarization vectors
> -         if(j_fks.eq.2)then
> +         if(j_fks.eq.2 .and. nincoming.eq.2)then
>            cphi_mother=-1.d0
>            sphi_mother=0.d0
>          else
> @@ -3477,7 +3507,11 @@
> 
> c Consistency check -- call to set_cms_stuff() must be done prior to
> c entering this function
> -      shattmp=2d0*dot(p(0,1),p(0,2))
> +      if (nincoming.eq.2) then
> +         shattmp=2d0*dot(p(0,1),p(0,2))
> +      else
> +         shattmp=p(0,1)**2
> +      endif
>       if(abs(shattmp/shat-1.d0).gt.1.d-5)then
>         write(*,*)'Error in sreal: inconsistent shat'
>         write(*,*)shattmp,shat
> @@ -3786,7 +3820,7 @@
>         enddo
>       enddo
>       do i=0,3
> -        if(j_fks.gt.2)then
> +        if(j_fks.gt.nincoming)then
>           xnum=p1_cnt(i,i_fks,inum)+p1_cnt(i,j_fks,inum)
>           xden=p1_cnt(i,i_fks,iden)+p1_cnt(i,j_fks,iden)
>         else
> @@ -4386,7 +4420,11 @@
> 
> c Consistency check -- call to set_cms_stuff() must be done prior to
> c entering this function
> -         shattmp=2d0*dot(p(0,1),p(0,2))
> +         if (nincoming.eq.2) then
> +            shattmp=2d0*dot(p(0,1),p(0,2))
> +         else
> +            shattmp=p(0,1)**2
> +         endif
>          if(abs(shattmp/shat-1.d0).gt.1.d-5)then
>            write(*,*)'Error in sreal: inconsistent shat'
>            write(*,*)shattmp,shat
> @@ -4463,7 +4501,7 @@
>          enddo
> c
>          do i=1,nincoming
> -            if (particle_type(i).ne.1)then
> +            if (particle_type(i).ne.1 .and. pmass(i).eq.ZERO) then
>                if (particle_type(i).eq.8) then
>                   aj=0
>                elseif(abs(particle_type(i)).eq.3) then
> @@ -4988,8 +5026,13 @@
>           elseif(abrv.ne.'virt' .and. abrv.ne.'viSC' .and.
>      #           abrv.ne.'viLC')then
> c 1+2+3+4
> -            tmp=dlog(xicut_used**2*shat/QES2)-
> -     #          1/betai*dlog((1+betai)/(1-betai))
> +            if (betai.gt.1d-6) then
> +               tmp=dlog(xicut_used**2*shat/QES2)-
> +     &              1/betai*dlog((1+betai)/(1-betai))
> +            else
> +               tmp=dlog(xicut_used**2*shat/QES2)-
> +     &              2d0*(1d0+betai**2/3d0+betai**4/5d0)
> +            endif
>           else
>              write(*,*)'Error #14 in eikonal_Ireg',abrv
>              stop
> @@ -5216,7 +5259,11 @@
>               elseif(pmass(m).ne.zero.and.pmass(n).ne.zero)then
>                 kikj=dot(p(0,n),p(0,m))
>                 vij=sqrt(1-(pmass(n)*pmass(m)/kikj)**2)
> -                single=single+0.5d0*1/vij*log((1+vij)/(1-vij))*wgt
> +                if (vij .gt. 1d-6) then
> +                   single=single+0.5d0*1/vij*log((1+vij)/(1-vij))*wgt
> +                else
> +                   single=single+(1d0+vij**2/3d0+vij**4/5d0)*wgt
> +                endif
>               else
>                 write(*,*)'Error in getpoles',i,j,n,m,pmass(n),pmass(m)
>                 stop
> 
> === modified file 'Template/NLO/SubProcesses/genps_fks.f'
> --- Template/NLO/SubProcesses/genps_fks.f	2015-02-04 13:26:19 +0000
> +++ Template/NLO/SubProcesses/genps_fks.f	2015-10-08 17:23:54 +0000
> @@ -182,15 +182,19 @@
>          endif
>       enddo
>       if( firsttime .or. iconfig0.ne.iconfigsave ) then
> -         stot = 4d0*ebeam(1)*ebeam(2)
> +         if (nincoming.eq.2) then
> +            stot = 4d0*ebeam(1)*ebeam(2)
> +         else
> +            stot=pmass(1)**2
> +         endif
> c Make sure have enough mass for external particles
>          totmassin=0d0
> -         do i=3-nincoming,2
> +         do i=1,nincoming
>             totmassin=totmassin+m(i)
>          enddo
>          totmass=0d0
>          nbranch = nexternal-3 ! nexternal is for n+1-body, while itree uses n-body
> -         do i=3,nbranch+2
> +         do i=nincoming+1,nexternal-1
>             totmass=totmass+m(i)
>          enddo
>          fksmass=totmass
> @@ -294,7 +298,7 @@
>       else
>          pb(0,1)=sqrtshat_born
>          do i=1,2
> -            pb(0,1)=0d0
> +            pb(i,1)=0d0
>          enddo
>          p(3,1)=1e-14           ! For HELAS routine ixxxxx for neg. mass
>       endif
> @@ -453,22 +457,22 @@
>             isolsign=1
>             call generate_momenta_massless_final(icountevts,i_fks,j_fks
>      &           ,p_born(0,imother),shat,sqrtshat ,x(ixEi),xmrec2,xp
> -     &           ,phi_i_fks,   xiimax,xinorm,xi_i_fks,y_ij_fks,p_i_fks
> +     &           ,phi_i_fks,xiimax,xinorm,xi_i_fks,y_ij_fks,p_i_fks
>      &           ,xjac,xpswgt,pass)
>             if (.not.pass) goto 112
>          elseif(m_j_fks.gt.0d0) then
>             call generate_momenta_massive_final(icountevts,isolsign
> -     &           ,i_fks,j_fks,p_born(0,imother),shat,sqrtshat ,m_j_fks
> -     &           ,x(ixEi),xmrec2,xp,phi_i_fks,    xiimax,xinorm
> +     &           ,i_fks,j_fks,p_born(0,imother),shat,sqrtshat,m_j_fks
> +     &           ,x(ixEi),xmrec2,xp,phi_i_fks,xiimax,xinorm
>      &           ,xi_i_fks,y_ij_fks,p_i_fks,xjac,xpswgt,pass)
>             if (.not.pass) goto 112
>          endif
>       elseif(j_fks.le.nincoming) then
>          isolsign=1
>          call generate_momenta_initial(icountevts,i_fks,j_fks,xbjrk_born
> -     &        ,tau_born,ycm_born,ycmhat,shat_born,phi_i_fks ,xp ,x(ixEi)
> -     &        ,shat,stot ,sqrtshat,tau,ycm,xbjrk ,p_i_fks,xiimax,xinorm
> -     &        ,xi_i_fks,y_ij_fks ,xpswgt,xjac ,pass)
> +     &        ,tau_born,ycm_born,ycmhat,shat_born,phi_i_fks,xp ,x(ixEi)
> +     &        ,shat,stot,sqrtshat,tau,ycm,xbjrk,p_i_fks,xiimax,xinorm
> +     &        ,xi_i_fks,y_ij_fks ,xpswgt,xjac,pass)
>          if (.not.pass) goto 112
>       else
>          write (*,*) 'Error #2 in genps_fks.f',j_fks
> @@ -750,6 +754,7 @@
> c remove the following if no importance sampling towards soft
> c singularity is performed when integrating over xi_i_hat
>       xjac=xjac*2d0*x(1)
> +
> c Check that xii is in the allowed range
>       if( icountevts.eq.-100 .or. abs(icountevts).eq.1 )then
>          if(xi_i_fks.gt.(1-xmrec2/shat))then
> @@ -829,7 +834,11 @@
> c mother four momenta
>       do i=0,3
>          xp_mother(i)=xp(i,i_fks)+xp(i,j_fks)
> -         recoil(i)=xp(i,1)+xp(i,2)-xp_mother(i)
> +         if (nincoming.eq.2) then
> +            recoil(i)=xp(i,1)+xp(i,2)-xp_mother(i)
> +         else
> +            recoil(i)=xp(i,1)-xp_mother(i)
> +         endif
>       enddo
>       sumrec=recoil(0)+rho(recoil)
>       sumrec2=sumrec**2
> @@ -843,7 +852,7 @@
>          xdir(j)=xp_mother(j)/x3len_fks_mother
>       enddo
> c Perform the boost here
> -      do i=3,nexternal
> +      do i=nincoming+1,nexternal
>          if(i.ne.i_fks.and.i.ne.j_fks.and.shybst.ne.0.d0)
>      &      call boostwdir2(chybst,shybst,chybstmo,xdir,xp(0,i),xp(0,i))
>       enddo
> @@ -1145,7 +1154,11 @@
> c mother four momenta
>       do i=0,3
>          xp_mother(i)=xp(i,i_fks)+xp(i,j_fks)
> -         recoil(i)=xp(i,1)+xp(i,2)-xp_mother(i)
> +         if (nincoming.eq.2) then
> +            recoil(i)=xp(i,1)+xp(i,2)-xp_mother(i)
> +         else
> +            recoil(i)=xp(i,1)-xp_mother(i)
> +         endif
>       enddo
> c
>       sumrec=recoil(0)+rho(recoil)
> @@ -1168,7 +1181,7 @@
>          xdir(j)=xp_mother(j)/x3len_fks_mother
>       enddo
> c Boost the momenta
> -      do i=3,nexternal
> +      do i=nincoming+1,nexternal
>          if(i.ne.i_fks.and.i.ne.j_fks.and.shybst.ne.0.d0)
>      &      call boostwdir2(chybst,shybst,chybstmo,xdir,xp(0,i),xp(0,i))
>       enddo
> @@ -2357,7 +2370,11 @@
>       external dot
>       pass=.true.
>       do i=0,3
> -         recoilbar(i)=p_born(i,1)+p_born(i,2)-p_born(i,imother)
> +         if (nincoming.eq.2) then
> +            recoilbar(i)=p_born(i,1)+p_born(i,2)-p_born(i,imother)
> +         else
> +            recoilbar(i)=p_born(i,1)-p_born(i,imother)
> +         endif
>       enddo
>       xmrec2=dot(recoilbar,recoilbar)
>       if(xmrec2.lt.0.d0)then
> 
> === modified file 'Template/NLO/SubProcesses/setcuts.f'
> --- Template/NLO/SubProcesses/setcuts.f	2015-07-24 18:59:18 +0000
> +++ Template/NLO/SubProcesses/setcuts.f	2015-10-08 17:23:54 +0000
> @@ -91,8 +91,6 @@
>             lpp(2)=0
>             ebeam(1)=pmass(1)/2d0
>             ebeam(2)=pmass(1)/2d0
> -            scale=pmass(1)
> -            fixed_ren_scale=.true.
>          endif
> c-check consistency of maxjetflavor in the run_card and with Nf
> c specified in coupl.inc
> 
> === modified file 'Template/NLO/SubProcesses/setscales.f'
> --- Template/NLO/SubProcesses/setscales.f	2015-09-11 13:19:28 +0000
> +++ Template/NLO/SubProcesses/setscales.f	2015-10-08 17:23:54 +0000
> @@ -237,7 +237,10 @@
>      $     ,FxFx_fac_scale
> c
>       tmp=0
> -      if(ickkw.eq.3)then
> +      if (nincoming.eq.1) then
> +         tmp=pp(0,1) ! mass of the decaying particle
> +         temp_scale_id='Mass of decaying particle'
> +      elseif(ickkw.eq.3)then
> c FxFx merging scale:
> c     Note that nFxFx_ren_scales includes the one scale that corresponds
> c     to the real-emission one (and is zero for the n-body conf.). Skip
> @@ -490,7 +493,10 @@
>       parameter (iQEStype=1)
> c
>       tmp=0
> -      if(iQEStype.eq.1)then
> +      if (nincoming.eq.1) then
> +         tmp=pp(0,1) ! mass of the decaying particle
> +         temp_scale_id='Mass of decaying particle'
> +      elseif(iQEStype.eq.1)then
>         tmp=scale_global_reference(pp)
>       elseif(iQEStype.eq.2)then
>         do i=nincoming+1,nexternal
> 
> === modified file 'UpdateNotes.txt'
> --- UpdateNotes.txt	2015-10-05 14:14:43 +0000
> +++ UpdateNotes.txt	2015-10-08 17:23:54 +0000
> @@ -9,6 +9,7 @@
> 	    processes which have jets at Born level, but do not require generation cut (like t-channel single-top). 
> 	RF: Added an option to the run_card to allow for easier variation of the shower starting scale (NLO only).
> 	RF: Fixed a problem in the setting of the flavour map used for runs with iAPPL >= 1. 
> +	RF: Allow for decay processes to compute (partial) decay widths at NLO accuracy (fixed order only).
> 
> 2.3.2(20/08/15)   
>         OM: Improve reweighting module. (https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Reweight)
> 
> === modified file 'madgraph/fks/fks_base.py'
> --- madgraph/fks/fks_base.py	2015-10-01 16:00:08 +0000
> +++ madgraph/fks/fks_base.py	2015-10-08 17:23:54 +0000
> @@ -550,10 +550,19 @@
>         """finds the FKS real configurations for a given process"""
>         if range(len(self.leglist)) != [l['number']-1 for l in self.leglist]:
>             raise fks_common.FKSProcessError('Disordered numbers of leglist')
> +
> +        if [ i['state'] for i in self.leglist].count(False) == 1:
> +            decay_process=True
> +        else:
> +            decay_process=False
> +
>         for i in self.leglist:
>             i_i = i['number'] - 1
>             self.reals.append([])
> -            self.splittings[i_i] = fks_common.find_splittings(i, self.born_proc['model'], {}, pert_order)
> +            if decay_process and not i['state']:
> +                self.splittings[i_i]=[]
> +            else:
> +                self.splittings[i_i] = fks_common.find_splittings(i, self.born_proc['model'], {}, pert_order)
>             for split in self.splittings[i_i]:
>                 self.reals[i_i].append(
>                             fks_common.insert_legs(self.leglist, i, split,pert=pert_order))
> 
> === modified file 'madgraph/interface/amcatnlo_interface.py'
> --- madgraph/interface/amcatnlo_interface.py	2015-10-01 16:00:08 +0000
> +++ madgraph/interface/amcatnlo_interface.py	2015-10-08 17:23:54 +0000
> @@ -542,8 +542,13 @@
>                         me.get('processes')[0].set('uid', uid)
>                         for fksreal in me.real_processes:
>                         # Pick out all initial state particles for the two beams
> -                            initial_states.append(sorted(list(set((p.get_initial_pdg(1),p.get_initial_pdg(2)) for \
> -                                                             p in fksreal.matrix_element.get('processes')))))
> +                            try:
> +                                initial_states.append(sorted(list(set((p.get_initial_pdg(1),p.get_initial_pdg(2)) for \
> +                                                             p in fksreal.matrix_element.get('processes')))))
> +                            except IndexError:
> +                                initial_states.append(sorted(list(set((p.get_initial_pdg(1)) for \
> +                                                             p in fksreal.matrix_element.get('processes')))))
> +                                
> 
>                     # remove doubles from the list
>                     checked = []
> 
> === modified file 'madgraph/interface/amcatnlo_run_interface.py'
> --- madgraph/interface/amcatnlo_run_interface.py	2015-10-08 11:31:22 +0000
> +++ madgraph/interface/amcatnlo_run_interface.py	2015-10-08 17:23:54 +0000
> @@ -1332,6 +1332,8 @@
>             return
> 
>         elif mode in ['aMC@NLO','aMC@LO','noshower','noshowerLO']:
> +            if self.ninitial == 1:
> +                raise aMCatNLOError('Decay processes can only be run at fixed order.')
>             mode_dict = {'aMC@NLO': 'all', 'aMC@LO': 'born',\
>                          'noshower': 'all', 'noshowerLO': 'born'}
>             shower = self.run_card['parton_shower'].upper()
> @@ -2109,7 +2111,15 @@
>         proc_info = '\n      Process %s\n      Run at %s-%s collider (%s + %s GeV)' % \
>         (process[:-3], lpp[self.run_card['lpp1']], lpp[self.run_card['lpp2']], 
>                 self.run_card['ebeam1'], self.run_card['ebeam2'])
> -        
> +
> +        if self.ninitial == 1:
> +            self.cross_sect_dict['unit']='GeV'
> +            self.cross_sect_dict['xsec_string']='(Partial) decay width'
> +            self.cross_sect_dict['axsec_string']='(Partial) abs(decay width)'
> +        else:
> +            self.cross_sect_dict['unit']='pb'
> +            self.cross_sect_dict['xsec_string']='Total cross-section'
> +            self.cross_sect_dict['axsec_string']='Total abs(cross-section)'
>         # Gather some basic statistics for the run and extracted from the log files.
>         if mode in ['aMC@NLO', 'aMC@LO', 'noshower', 'noshowerLO']: 
>             log_GV_files =  glob.glob(pjoin(self.me_dir, \
> @@ -2135,13 +2145,13 @@
>             if step != 2:
>                 message = status[step] + '\n\n      Intermediate results:' + \
>                     ('\n      Random seed: %(randinit)d' + \
> -                     '\n      Total cross-section:      %(xsect)8.3e +- %(errt)6.1e pb' + \
> -                     '\n      Total abs(cross-section): %(xseca)8.3e +- %(erra)6.1e pb \n') \
> +                     '\n      %(xsec_string)s:      %(xsect)8.3e +- %(errt)6.1e %(unit)s' + \
> +                     '\n      %(axsec_string)s: %(xseca)8.3e +- %(erra)6.1e %(unit)s \n') \
>                      % self.cross_sect_dict
>             else:
> 
>                 message = '\n      ' + status[step] + proc_info + \
> -                          '\n      Total cross-section: %(xsect)8.3e +- %(errt)6.1e pb' % \
> +                          '\n      %(xsec_string)s: %(xsect)8.3e +- %(errt)6.1e %(unit)s' % \
>                         self.cross_sect_dict
> 
>                 if self.run_card['nevents']>=10000 and self.run_card['reweight_scale']:
> @@ -2170,15 +2180,15 @@
>                       'Final results and run summary:']
>             if (not done) and (step == 0):
>                 message = '\n      ' + status[0] + \
> -                     '\n      Total cross-section:      %(xsect)8.3e +- %(errt)6.1e pb' % \
> +                     '\n      %(xsec_string)s:      %(xsect)8.3e +- %(errt)6.1e %(unit)s' % \
>                              self.cross_sect_dict
>             elif not done:
>                 message = '\n      ' + status[1] + \
> -                     '\n      Total cross-section:      %(xsect)8.3e +- %(errt)6.1e pb' % \
> +                     '\n      %(xsec_string)s:      %(xsect)8.3e +- %(errt)6.1e %(unit)s' % \
>                              self.cross_sect_dict
>             elif done:
>                 message = '\n      ' + status[2] + proc_info + \
> -                     '\n      Total cross-section:      %(xsect)8.3e +- %(errt)6.1e pb' % \
> +                     '\n      %(xsec_string)s:      %(xsect)8.3e +- %(errt)6.1e %(unit)s' % \
>                              self.cross_sect_dict
>                 if self.run_card['reweight_scale']:
>                     if self.run_card['ickkw'] != -1:
> 
> === modified file 'madgraph/iolibs/export_fks.py'
> --- madgraph/iolibs/export_fks.py	2015-10-01 16:00:08 +0000
> +++ madgraph/iolibs/export_fks.py	2015-10-08 17:23:54 +0000
> @@ -383,8 +383,11 @@
>             text=text+str(i+1)+' '+str(len(e))
>             for t in e:
>                 text=text+'   '
> -                for p in t:
> -                    text=text+' '+str(p)
> +                try:
> +                    for p in t:
> +                        text=text+' '+str(p)
> +                except TypeError:
> +                        text=text+' '+str(t)
>             text=text+'\n'
> 
>         ff = open(file_pos, 'w')
> 
> === modified file 'madgraph/iolibs/template_files/parton_lum_n_fks.inc'
> --- madgraph/iolibs/template_files/parton_lum_n_fks.inc	2013-05-31 14:58:30 +0000
> +++ madgraph/iolibs/template_files/parton_lum_n_fks.inc	2015-10-08 17:23:54 +0000
> @@ -54,7 +54,11 @@
> %(pdf_lines)s
>       endif
>       do i=1,IPROC
> -        lum = lum + pd(i) * conv 
> +        if (nincoming.eq.2) then
> +           lum = lum + pd(i) * conv 
> +        else
> +           lum = lum + pd(i)
> +        endif
>       enddo
>       return
>       END                                                             
> 
> === modified file 'models/loop_qcd_qed_sm/object_library.py'
> --- models/loop_qcd_qed_sm/object_library.py	2015-10-01 16:00:08 +0000
> +++ models/loop_qcd_qed_sm/object_library.py	2015-10-08 17:23:54 +0000
> @@ -288,7 +288,7 @@
>                     #else:
>                 def substitution(matchedObj):
>                     return matchedObj.group('first')+"("+ctpar.pole(x)+")"+matchedObj.group('second')
> -                       pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\(|\s)(?P<name>"+ctpar.name+r")(?P<second>\Z|\*|\+|\-|\)|/|\\|\s)")
> +                pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\(|\s)(?P<name>"+ctpar.name+r")(?P<second>\Z|\*|\+|\-|\)|/|\\|\s)")
>                 tempvalue2 = tempvalue
>                 tempvalue = pattern.sub(substitution,tempvalue)
>                 num = numbermatch[i]-1		
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%parton_lum_0.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%parton_lum_0.f	2015-03-10 08:54:03 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%parton_lum_0.f	2015-10-08 17:23:54 +0000
> @@ -92,7 +92,11 @@
>         PD(IPROC) = G1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%parton_lum_0.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%parton_lum_0.f	2015-03-10 08:54:03 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%parton_lum_0.f	2015-10-08 17:23:54 +0000
> @@ -119,7 +119,11 @@
>         PD(IPROC) = S1*SX2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%parton_lum_0.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%parton_lum_0.f	2015-03-10 08:54:03 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%parton_lum_0.f	2015-10-08 17:23:54 +0000
> @@ -119,7 +119,11 @@
>         PD(IPROC) = SX1*S2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_1.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_1.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_1.f	2015-10-08 17:23:54 +0000
> @@ -92,7 +92,11 @@
>         PD(IPROC) = G1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_2.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_2.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_2.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = CX1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_3.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_3.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_3.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = C1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_4.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_4.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_4.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = G1*CX2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_5.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_5.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_5.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = G1*C2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_1.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_1.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_1.f	2015-10-08 17:23:54 +0000
> @@ -119,7 +119,11 @@
>         PD(IPROC) = S1*SX2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_2.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_2.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_2.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = G1*SX2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_3.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_3.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_3.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = S1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_1.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_1.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_1.f	2015-10-08 17:23:54 +0000
> @@ -119,7 +119,11 @@
>         PD(IPROC) = SX1*S2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_2.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_2.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_2.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = G1*S2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_3.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_3.f	2014-10-20 01:07:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%parton_lum_3.f	2015-10-08 17:23:54 +0000
> @@ -113,7 +113,11 @@
>         PD(IPROC) = SX1*G2
>       ENDIF
>       DO I=1,IPROC
> -        LUM = LUM + PD(I) * CONV
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
>       ENDDO
>       RETURN
>       END
> 
> === added directory 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal'
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_001.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_001.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_001.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,153 @@
> +      SUBROUTINE SB_SF_001(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 1 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_001
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_001(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_001(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 1 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/4D0*(-3D0*AMP(1)+1D0/3D0*AMP(1))
> +      B_SF_001 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_001 =B_SF_001+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_002.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_002.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_002.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,153 @@
> +      SUBROUTINE SB_SF_002(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 2 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_002
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_002(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_002(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 2 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(+3D0*AMP(1)-1D0/3D0*AMP(1))
> +      B_SF_002 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_002 =B_SF_002+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_003.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_003.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_003.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,154 @@
> +      SUBROUTINE SB_SF_003(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 3 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_003
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_003(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_003(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 3 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=2)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  2) /    9,    3/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(-1D0/3D0*AMP(1))
> +      JAMP2(2)=+1D0/2D0*(+AMP(1))
> +      B_SF_003 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_003 =B_SF_003+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_004.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_004.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_004.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,154 @@
> +      SUBROUTINE SB_SF_004(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_004
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_004(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_004(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 1 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=2)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  2) /    9,    3/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(+1D0/3D0*AMP(1))
> +      JAMP2(2)=+1D0/2D0*(-AMP(1))
> +      B_SF_004 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_004 =B_SF_004+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_005.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_005.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_005.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,153 @@
> +      SUBROUTINE SB_SF_005(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 2 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_005
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_005(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_005(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 2 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/4D0*(-3D0*AMP(1)+1D0/3D0*AMP(1))
> +      B_SF_005 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_005 =B_SF_005+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_006.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_006.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_006.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,154 @@
> +      SUBROUTINE SB_SF_006(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 3 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_006
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_006(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_006(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 3 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=2)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  2) /    9,    3/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(+1D0/3D0*AMP(1))
> +      JAMP2(2)=+1D0/2D0*(-AMP(1))
> +      B_SF_006 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_006 =B_SF_006+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_007.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_007.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_007.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,154 @@
> +      SUBROUTINE SB_SF_007(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_007
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_007(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_007(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 2 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=2)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  2) /    9,    3/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(-1D0/3D0*AMP(1))
> +      JAMP2(2)=+1D0/2D0*(+AMP(1))
> +      B_SF_007 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_007 =B_SF_007+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_008.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_008.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%b_sf_008.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,153 @@
> +      SUBROUTINE SB_SF_008(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 3 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 B_SF_008
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE(*,*) 'Error in sb_sf: momenta not the same in th'
> +     $       //'e born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in sb_sf: color_linked borns should b'
> +     $   //'e called only with calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          ANS=ANS+B_SF_008(P1,IHEL)
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION B_SF_008(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     spectators: 3 4 
> +
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER NCOLOR1, NCOLOR2
> +      PARAMETER (NCOLOR1=1, NCOLOR2=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in b_sf: color_linked borns should be calle'
> +     $   //'d only with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP1(1)=-AMP(1)
> +      JAMP2(1)=+1D0/2D0*(+3D0*AMP(1)-1D0/3D0*AMP(1))
> +      B_SF_008 = 0.D0
> +      DO I = 1, NCOLOR1
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR2
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
> +        ENDDO
> +        B_SF_008 =B_SF_008+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,298 @@
> +      SUBROUTINE SBORN(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P1(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +      INCLUDE 'genps.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS=   1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1)
> +      COMPLEX*16 ANS(2)
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,I,J,JJ,GLU_IJ
> +      REAL*8 BORN,BORNS(2)
> +      COMPLEX*16 BORNTILDE
> +      INTEGER NTRY(3)
> +      DATA NTRY /3*0/
> +      INTEGER NHEL(NEXTERNAL-1,NCOMB)
> +      DATA (NHEL(I,   1),I=1,4) / 1,-1,-1, 1/
> +      DATA (NHEL(I,   2),I=1,4) / 1,-1,-1,-1/
> +      DATA (NHEL(I,   3),I=1,4) / 1,-1, 1, 1/
> +      DATA (NHEL(I,   4),I=1,4) / 1,-1, 1,-1/
> +      DATA (NHEL(I,   5),I=1,4) / 1, 1,-1, 1/
> +      DATA (NHEL(I,   6),I=1,4) / 1, 1,-1,-1/
> +      DATA (NHEL(I,   7),I=1,4) / 1, 1, 1, 1/
> +      DATA (NHEL(I,   8),I=1,4) / 1, 1, 1,-1/
> +      DATA (NHEL(I,   9),I=1,4) /-1,-1,-1, 1/
> +      DATA (NHEL(I,  10),I=1,4) /-1,-1,-1,-1/
> +      DATA (NHEL(I,  11),I=1,4) /-1,-1, 1, 1/
> +      DATA (NHEL(I,  12),I=1,4) /-1,-1, 1,-1/
> +      DATA (NHEL(I,  13),I=1,4) /-1, 1,-1, 1/
> +      DATA (NHEL(I,  14),I=1,4) /-1, 1,-1,-1/
> +      DATA (NHEL(I,  15),I=1,4) /-1, 1, 1, 1/
> +      DATA (NHEL(I,  16),I=1,4) /-1, 1, 1,-1/
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +      INTEGER IJ_VALUES(3)
> +      DATA IJ_VALUES /2, 3, 4/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXAMPS)
> +      COMMON/TO_AMPS/  AMP2,       JAMP2
> +      DATA JAMP2(0) /   1/
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/GOODHEL
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      DOUBLE PRECISION HEL_FAC
> +      INTEGER GET_HEL,SKIP(3)
> +      COMMON/CBORN/HEL_FAC,GET_HEL,SKIP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      GLU_IJ = IJ_VALUES(NFKSPROCESS)
> +      NTRY(NFKSPROCESS)=NTRY(NFKSPROCESS)+1
> +      IF (NTRY(NFKSPROCESS).LT.2) THEN
> +        SKIP(NFKSPROCESS)=1
> +        DO WHILE(NHEL(GLU_IJ ,SKIP(NFKSPROCESS)).NE.-NHEL(GLU_IJ ,1))
> +          SKIP(NFKSPROCESS)=SKIP(NFKSPROCESS)+1
> +        ENDDO
> +        SKIP(NFKSPROCESS)=SKIP(NFKSPROCESS)-1
> +      ENDIF
> +      DO JJ=1,NGRAPHS
> +        AMP2(JJ)=0D0
> +      ENDDO
> +      DO JJ=1,INT(JAMP2(0))
> +        JAMP2(JJ)=0D0
> +      ENDDO
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE (*,*) 'momenta not the same in Born'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ENDIF
> +      IF (.NOT.CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          SAVEMOM(J,1)=P1(0,J)
> +          SAVEMOM(J,2)=P1(3,J)
> +        ENDDO
> +        DO J=1,MAX_BHEL
> +          DO JJ=1,NGRAPHS
> +            SAVEAMP(JJ,J)=(0D0,0D0)
> +          ENDDO
> +        ENDDO
> +      ENDIF
> +      ANS(1) = 0D0
> +      ANS(2) = 0D0
> +      HEL_FAC=1D0
> +      DO IHEL=1,NCOMB
> +        IF (NHEL(GLU_IJ,IHEL).EQ.NHEL(GLU_IJ,1)) THEN
> +          IF ((GOODHEL(IHEL,NFKSPROCESS) .OR. GOODHEL(IHEL+SKIP(NFKSPRO
> +     $     CESS),NFKSPROCESS) .OR. NTRY(NFKSPROCESS) .LT. 2) ) THEN
> +            ANS(1)=ANS(1)+BORN(P1,NHEL(1,IHEL),IHEL,BORNTILDE,BORNS)
> +            ANS(2)=ANS(2)+BORNTILDE
> +            IF ( BORNS(1).NE.0D0 .AND. .NOT. GOODHEL(IHEL,NFKSPROCESS
> +     $       ) ) THEN
> +              GOODHEL(IHEL,NFKSPROCESS)=.TRUE.
> +            ENDIF
> +            IF ( BORNS(2).NE.0D0 .AND. .NOT. GOODHEL(IHEL+SKIP(NFKSPROC
> +     $       ESS),NFKSPROCESS) ) THEN
> +              GOODHEL(IHEL+SKIP(NFKSPROCESS),NFKSPROCESS)=.TRUE.
> +            ENDIF
> +          ENDIF
> +        ENDIF
> +      ENDDO
> +      ANS(1)=ANS(1)/DBLE(IDEN)
> +      ANS(2)=ANS(2)/DBLE(IDEN)
> +      CALCULATEDBORN=.TRUE.
> +      END
> +
> +
> +      REAL*8 FUNCTION BORN(P,NHEL,HELL,BORNTILDE,BORNS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER    NGRAPHS,    NEIGEN
> +      PARAMETER (NGRAPHS=   1,NEIGEN=  1)
> +      INTEGER    NWAVEFUNCS, NCOLOR
> +      PARAMETER (NWAVEFUNCS=5, NCOLOR=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +      INCLUDE 'coupl.inc'
> +      INCLUDE 'genps.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1),BORNS(2)
> +      INTEGER NHEL(NEXTERNAL-1), HELL
> +      COMPLEX*16 BORNTILDE
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J,IHEL,BACK_HEL,GLU_IJ
> +      INTEGER IC(NEXTERNAL-1),NMO
> +      PARAMETER (NMO=NEXTERNAL-1)
> +      DATA IC /NMO*1/
> +      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS)
> +     $ , JAMPH(2, NCOLOR)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE PRECISION AMP2(MAXAMPS), JAMP2(0:MAXAMPS)
> +      COMMON/TO_AMPS/  AMP2,       JAMP2
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      DOUBLE PRECISION HEL_FAC
> +      INTEGER GET_HEL,SKIP(3)
> +      COMMON/CBORN/HEL_FAC,GET_HEL,SKIP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +      INTEGER STEP_HEL
> +      INTEGER IJ_VALUES(3)
> +      DATA IJ_VALUES /2, 3, 4/
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     1 T(2,1) T(3,4)
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      GLU_IJ = IJ_VALUES(NFKSPROCESS)
> +      BORN = 0D0
> +      BORNTILDE = (0D0,0D0)
> +      BACK_HEL = NHEL(GLU_IJ)
> +      BORNS(1) = 0D0
> +      BORNS(2) = 0D0
> +      IF (BACK_HEL.NE.0) THEN
> +        STEP_HEL=-2*BACK_HEL
> +      ELSE
> +        STEP_HEL=1
> +      ENDIF
> +      DO IHEL=BACK_HEL,-BACK_HEL,STEP_HEL
> +        IF (IHEL.EQ.BACK_HEL.OR.NHEL(GLU_IJ).NE.0) THEN
> +          IF (NHEL(GLU_IJ).NE.0) NHEL(GLU_IJ) = IHEL
> +          IF (.NOT. CALCULATEDBORN) THEN
> +            CALL IXXXXX(P(0,1),MDL_MT,NHEL(1),+1*IC(1),W(1,1))
> +            CALL OXXXXX(P(0,2),MDL_MB,NHEL(2),+1*IC(2),W(1,2))
> +            CALL OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
> +            CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
> +            CALL FFV2_3(W(1,4),W(1,3),GC_100,MDL_MW,MDL_WW,W(1,5))
> +C           Amplitude(s) for diagram number 1
> +            CALL FFV2_0(W(1,1),W(1,2),W(1,5),GC_100,AMP(1))
> +            DO I=1,NGRAPHS
> +              IF(IHEL.EQ.BACK_HEL)THEN
> +                SAVEAMP(I,HELL)=AMP(I)
> +              ELSEIF(IHEL.EQ.-BACK_HEL)THEN
> +                SAVEAMP(I,HELL+SKIP(NFKSPROCESS))=AMP(I)
> +              ELSE
> +                WRITE(*,*) 'ERROR #1 in born.f'
> +                STOP
> +              ENDIF
> +            ENDDO
> +          ELSEIF (CALCULATEDBORN) THEN
> +            DO I=1,NGRAPHS
> +              IF(IHEL.EQ.BACK_HEL)THEN
> +                AMP(I)=SAVEAMP(I,HELL)
> +              ELSEIF(IHEL.EQ.-BACK_HEL)THEN
> +                AMP(I)=SAVEAMP(I,HELL+SKIP(NFKSPROCESS))
> +              ELSE
> +                WRITE(*,*) 'ERROR #1 in born.f'
> +                STOP
> +              ENDIF
> +            ENDDO
> +          ENDIF
> +          JAMP(1)=-AMP(1)
> +          DO I = 1, NCOLOR
> +            ZTEMP = (0.D0,0.D0)
> +            DO J = 1, NCOLOR
> +              ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
> +            ENDDO
> +            BORNS(2-(1+BACK_HEL*IHEL)/2)=BORNS(2-(1+BACK_HEL*IHEL)/2)
> +     $       +ZTEMP*DCONJG(JAMP(I))/DENOM(I)
> +          ENDDO
> +          DO I = 1, NGRAPHS
> +            AMP2(I)=AMP2(I)+AMP(I)*DCONJG(AMP(I))
> +          ENDDO
> +          DO I = 1, NCOLOR
> +            JAMP2(I)=JAMP2(I)+JAMP(I)*DCONJG(JAMP(I))
> +            JAMPH(2-(1+BACK_HEL*IHEL)/2,I)=JAMP(I)
> +          ENDDO
> +        ENDIF
> +      ENDDO
> +      BORN=BORNS(1)+BORNS(2)
> +      DO I = 1, NCOLOR
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR
> +          ZTEMP = ZTEMP + CF(J,I)*JAMPH(2,J)
> +        ENDDO
> +        BORNTILDE = BORNTILDE + ZTEMP*DCONJG(JAMPH(1,I))/DENOM(I)
> +      ENDDO
> +      NHEL(GLU_IJ) = BACK_HEL
> +      END
> +
> +
> +      BLOCK DATA GOODHELS
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/GOODHEL
> +      DATA GOODHEL/THEL*.FALSE./
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_conf.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_conf.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_conf.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,8 @@
> +C     Diagram 1, Amplitude 1
> +      DATA MAPCONFIG(   1)/   1/
> +      DATA (IFOREST(I, -1,   1),I=1,2)/  4,  3/
> +      DATA SPROP(  -1,   1)/      24/
> +      DATA (IFOREST(I, -2,   1),I=1,2)/ -1,  2/
> +      DATA SPROP(  -2,   1)/       6/
> +C     Number of configs
> +      DATA MAPCONFIG(0)/   1/
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_decayBW.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_decayBW.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_decayBW.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,2 @@
> +      DATA GFORCEBW(-1,1)/.FALSE./
> +      DATA GFORCEBW(-2,1)/.FALSE./
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_hel.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_hel.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_hel.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,151 @@
> +      SUBROUTINE SBORN_HEL(P1,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     AND HELICITIES
> +C     FOR THE POINT IN PHASE SPACE P1(0:3,NEXTERNAL-1)
> +C     
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=  16 )
> +      INTEGER    THEL
> +      PARAMETER (THEL=NCOMB*3)
> +      INTEGER NGRAPHS
> +      PARAMETER (NGRAPHS = 1)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P1(0:3,NEXTERNAL-1),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,J
> +      REAL*8 BORN_HEL
> +      INTEGER IDEN_VALUES(3)
> +      DATA IDEN_VALUES /6, 6, 6/
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      LOGICAL GOODHEL(NCOMB,3)
> +      COMMON /C_GOODHEL/ GOODHEL
> +      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
> +      COMMON/TO_SAVEMOM/SAVEMOM
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +      DOUBLE PRECISION WGT_HEL(MAX_BHEL)
> +      COMMON/C_BORN_HEL/WGT_HEL
> +
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IDEN=IDEN_VALUES(NFKSPROCESS)
> +      IF (CALCULATEDBORN) THEN
> +        DO J=1,NEXTERNAL-1
> +          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
> +     $     ,J)) THEN
> +            CALCULATEDBORN=.FALSE.
> +            WRITE (*,*) 'momenta not the same in Born_hel'
> +            STOP
> +          ENDIF
> +        ENDDO
> +      ELSE
> +        WRITE(*,*) 'Error in born_hel: should be called only wit'
> +     $   //'h calculatedborn = true'
> +        STOP
> +      ENDIF
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        WGT_HEL(IHEL)=0D0
> +        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
> +          WGT_HEL(IHEL)=BORN_HEL(P1,IHEL)/DBLE(IDEN)
> +          ANS=ANS+WGT_HEL(IHEL)
> +        ENDIF
> +      ENDDO
> +      END
> +
> +
> +      REAL*8 FUNCTION BORN_HEL(P,HELL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
> +C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
> +
> +C     Process: t > b u d~ [ real = QCD ]
> +C     Process: t > b c s~ [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER     NGRAPHS
> +      PARAMETER ( NGRAPHS = 1 )
> +      INTEGER    NCOLOR
> +      PARAMETER (NCOLOR=1)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1 = (0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'born_nhel.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL-1)
> +      INTEGER HELL
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR)
> +C     
> +C     GLOBAL VARIABLES
> +C     
> +      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
> +      COMMON/TO_SAVEAMP/SAVEAMP
> +      LOGICAL CALCULATEDBORN
> +      COMMON/CCALCULATEDBORN/CALCULATEDBORN
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  1) /    9/
> +C     1 T(2,1) T(3,4)
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      IF (.NOT. CALCULATEDBORN) THEN
> +        WRITE(*,*) 'Error in born_hel.f: this should be called onl'
> +     $   //'y with calculatedborn = true'
> +        STOP
> +      ELSEIF (CALCULATEDBORN) THEN
> +        DO I=1,NGRAPHS
> +          AMP(I)=SAVEAMP(I,HELL)
> +        ENDDO
> +      ENDIF
> +      JAMP(1)=-AMP(1)
> +      BORN_HEL = 0.D0
> +      DO I = 1, NCOLOR
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
> +        ENDDO
> +        BORN_HEL =BORN_HEL+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_leshouche.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_leshouche.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_leshouche.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,8 @@
> +      DATA (IDUP(I,1),I=1,4)/6,5,2,-1/
> +      DATA (MOTHUP(1,I,  1),I=1, 4)/  0,  1,  1,  1/
> +      DATA (MOTHUP(2,I,  1),I=1, 4)/  0,  2,  2,  2/
> +      DATA (ICOLUP(1,I,  1),I=1, 4)/501,501,502,  0/
> +      DATA (ICOLUP(2,I,  1),I=1, 4)/  0,  0,  0,502/
> +      DATA (IDUP(I,2),I=1,4)/6,5,4,-3/
> +      DATA (MOTHUP(1,I,  2),I=1, 4)/  0,  1,  1,  1/
> +      DATA (MOTHUP(2,I,  2),I=1, 4)/  0,  2,  2,  2/
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_maxamps.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_maxamps.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_maxamps.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,3 @@
> +      INTEGER    MAXAMPS, MAXFLOW, MAXPROC, MAXSPROC
> +      PARAMETER (MAXAMPS=1, MAXFLOW=1)
> +      PARAMETER (MAXPROC=2, MAXSPROC=1)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_ngraphs.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_ngraphs.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_ngraphs.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,2 @@
> +      INTEGER    N_MAX_CG
> +      PARAMETER (N_MAX_CG=1)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_nhel.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_nhel.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_nhel.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,3 @@
> +      INTEGER    MAX_BHEL, MAX_BCOL
> +      PARAMETER (MAX_BHEL=16)
> +      PARAMETER(MAX_BCOL=1)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_props.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_props.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%born_props.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,6 @@
> +      PMASS( -1,   1)  = ABS(MDL_MW)
> +      PWIDTH( -1,   1) = ABS(MDL_WW)
> +      POW( -1,   1) = 2
> +      PMASS( -2,   1)  = ABS(MDL_MT)
> +      PWIDTH( -2,   1) = ABS(MDL_WT)
> +      POW( -2,   1) = 1
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%coloramps.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%coloramps.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%coloramps.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,2 @@
> +      LOGICAL ICOLAMP(1,1,1)
> +      DATA(ICOLAMP(I,1,1),I=1,1)/.TRUE./
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%config_subproc_map.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%config_subproc_map.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%config_subproc_map.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,1 @@
> +      DATA CONFSUB(1,1)/1/
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%configs_and_props_decl.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%configs_and_props_decl.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%configs_and_props_decl.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,12 @@
> +      INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
> +      PARAMETER (LMAXCONFIGS_USED=   4)
> +      PARAMETER (MAX_BRANCH_USED =   3)
> +      INTEGER MAPCONFIG_D(  3,0:LMAXCONFIGS_USED)
> +      INTEGER IFOREST_D(  3,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
> +      INTEGER SPROP_D(  3,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
> +      INTEGER TPRID_D(  3,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
> +      DOUBLE PRECISION PMASS_D(  3,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED
> +     $ )
> +      DOUBLE PRECISION PWIDTH_D(  3,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USE
> +     $ D)
> +      INTEGER POW_D(  3,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%fks_info.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%fks_info.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%fks_info.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,46 @@
> +      INTEGER IPOS, JPOS
> +      INTEGER FKS_I_D(3), FKS_J_D(3)
> +      INTEGER FKS_J_FROM_I_D(3, NEXTERNAL, 0:NEXTERNAL)
> +      INTEGER PARTICLE_TYPE_D(3, NEXTERNAL), PDG_TYPE_D(3, NEXTERNAL)
> +      REAL*8 PARTICLE_CHARGE_D(3, NEXTERNAL)
> +
> +      DATA FKS_I_D / 5, 5, 5 /
> +      DATA FKS_J_D / 2, 3, 4 /
> +
> +      DATA (FKS_J_FROM_I_D(1, 5, JPOS), JPOS = 0, 4)  / 4, 1, 2, 3, 4 /
> +
> +      DATA (FKS_J_FROM_I_D(2, 5, JPOS), JPOS = 0, 4)  / 4, 1, 2, 3, 4 /
> +
> +      DATA (FKS_J_FROM_I_D(3, 5, JPOS), JPOS = 0, 4)  / 4, 1, 2, 3, 4 /
> +
> +
> +C     
> +C     Particle type:
> +C     octet = 8, triplet = 3, singlet = 1
> +      DATA (PARTICLE_TYPE_D(1, IPOS), IPOS=1, NEXTERNAL) / 3, 3, 3, 
> +     $ -3, 8 /
> +      DATA (PARTICLE_TYPE_D(2, IPOS), IPOS=1, NEXTERNAL) / 3, 3, 3, 
> +     $ -3, 8 /
> +      DATA (PARTICLE_TYPE_D(3, IPOS), IPOS=1, NEXTERNAL) / 3, 3, 3, 
> +     $ -3, 8 /
> +
> +C     
> +C     Particle type according to PDG:
> +C     
> +      DATA (PDG_TYPE_D(1, IPOS), IPOS=1, NEXTERNAL) / 6, 5, 2, -1, 21 /
> +      DATA (PDG_TYPE_D(2, IPOS), IPOS=1, NEXTERNAL) / 6, 5, 2, -1, 21 /
> +      DATA (PDG_TYPE_D(3, IPOS), IPOS=1, NEXTERNAL) / 6, 5, 2, -1, 21 /
> +
> +C     
> +C     Particle charge:
> +C     charge is set 0. with QCD corrections, which is irrelevant
> +      DATA (PARTICLE_CHARGE_D(1, IPOS), IPOS=1, NEXTERNAL) /  
> +     $  0.000000000000000D0,   0.000000000000000D0,   0.000000000000000
> +     $ D0,   0.000000000000000D0,   0.000000000000000D0 /
> +      DATA (PARTICLE_CHARGE_D(2, IPOS), IPOS=1, NEXTERNAL) /  
> +     $  0.000000000000000D0,   0.000000000000000D0,   0.000000000000000
> +     $ D0,   0.000000000000000D0,   0.000000000000000D0 /
> +      DATA (PARTICLE_CHARGE_D(3, IPOS), IPOS=1, NEXTERNAL) /  
> +     $  0.000000000000000D0,   0.000000000000000D0,   0.000000000000000
> +     $ D0,   0.000000000000000D0,   0.000000000000000D0 /
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%get_color.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%get_color.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%get_color.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,54 @@
> +      FUNCTION GET_COLOR(IPDG)
> +      IMPLICIT NONE
> +      INTEGER GET_COLOR, IPDG
> +
> +      IF(IPDG.EQ.-24)THEN
> +        GET_COLOR=1
> +        RETURN
> +      ELSE IF(IPDG.EQ.-6)THEN
> +        GET_COLOR=-3
> +        RETURN
> +      ELSE IF(IPDG.EQ.-5)THEN
> +        GET_COLOR=-3
> +        RETURN
> +      ELSE IF(IPDG.EQ.-3)THEN
> +        GET_COLOR=-3
> +        RETURN
> +      ELSE IF(IPDG.EQ.-2)THEN
> +        GET_COLOR=-3
> +        RETURN
> +      ELSE IF(IPDG.EQ.-1)THEN
> +        GET_COLOR=-3
> +        RETURN
> +      ELSE IF(IPDG.EQ.1)THEN
> +        GET_COLOR=3
> +        RETURN
> +      ELSE IF(IPDG.EQ.2)THEN
> +        GET_COLOR=3
> +        RETURN
> +      ELSE IF(IPDG.EQ.4)THEN
> +        GET_COLOR=3
> +        RETURN
> +      ELSE IF(IPDG.EQ.5)THEN
> +        GET_COLOR=3
> +        RETURN
> +      ELSE IF(IPDG.EQ.6)THEN
> +        GET_COLOR=3
> +        RETURN
> +      ELSE IF(IPDG.EQ.21)THEN
> +        GET_COLOR=8
> +        RETURN
> +      ELSE IF(IPDG.EQ.24)THEN
> +        GET_COLOR=1
> +        RETURN
> +      ELSE IF(IPDG.EQ.7)THEN
> +C       This is dummy particle used in multiparticle vertices
> +        GET_COLOR=2
> +        RETURN
> +      ELSE
> +        WRITE(*,*)'Error: No color given for pdg ',IPDG
> +        GET_COLOR=0
> +        RETURN
> +      ENDIF
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%leshouche_decl.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%leshouche_decl.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%leshouche_decl.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,6 @@
> +      INTEGER MAXPROC_USED, MAXFLOW_USED
> +      PARAMETER (MAXPROC_USED = 2)
> +      PARAMETER (MAXFLOW_USED = 2)
> +      INTEGER IDUP_D(3,5,MAXPROC_USED)
> +      INTEGER MOTHUP_D(3,2,5,MAXPROC_USED)
> +      INTEGER ICOLUP_D(3,2,5,MAXFLOW_USED)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%matrix_1.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%matrix_1.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%matrix_1.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,194 @@
> +      SUBROUTINE SMATRIX_1(P,ANS)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     Returns amplitude squared summed/avg over colors
> +C     and helicities
> +C     for the point in phase space P(0:3,NEXTERNAL)
> +C     
> +C     Process: t > b u d~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     Process: t > b c s~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INCLUDE 'nexternal.inc'
> +      INTEGER     NCOMB
> +      PARAMETER ( NCOMB=32)
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL),ANS
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER IHEL,IDEN,I,T_IDENT(NCOMB)
> +      REAL*8 MATRIX_1
> +      REAL*8 T,T_SAVE(NCOMB)
> +      SAVE T_SAVE,T_IDENT
> +      INTEGER NHEL(NEXTERNAL,NCOMB)
> +      DATA (NHEL(I,   1),I=1,5) / 1,-1,-1, 1,-1/
> +      DATA (NHEL(I,   2),I=1,5) / 1,-1,-1, 1, 1/
> +      DATA (NHEL(I,   3),I=1,5) / 1,-1,-1,-1,-1/
> +      DATA (NHEL(I,   4),I=1,5) / 1,-1,-1,-1, 1/
> +      DATA (NHEL(I,   5),I=1,5) / 1,-1, 1, 1,-1/
> +      DATA (NHEL(I,   6),I=1,5) / 1,-1, 1, 1, 1/
> +      DATA (NHEL(I,   7),I=1,5) / 1,-1, 1,-1,-1/
> +      DATA (NHEL(I,   8),I=1,5) / 1,-1, 1,-1, 1/
> +      DATA (NHEL(I,   9),I=1,5) / 1, 1,-1, 1,-1/
> +      DATA (NHEL(I,  10),I=1,5) / 1, 1,-1, 1, 1/
> +      DATA (NHEL(I,  11),I=1,5) / 1, 1,-1,-1,-1/
> +      DATA (NHEL(I,  12),I=1,5) / 1, 1,-1,-1, 1/
> +      DATA (NHEL(I,  13),I=1,5) / 1, 1, 1, 1,-1/
> +      DATA (NHEL(I,  14),I=1,5) / 1, 1, 1, 1, 1/
> +      DATA (NHEL(I,  15),I=1,5) / 1, 1, 1,-1,-1/
> +      DATA (NHEL(I,  16),I=1,5) / 1, 1, 1,-1, 1/
> +      DATA (NHEL(I,  17),I=1,5) /-1,-1,-1, 1,-1/
> +      DATA (NHEL(I,  18),I=1,5) /-1,-1,-1, 1, 1/
> +      DATA (NHEL(I,  19),I=1,5) /-1,-1,-1,-1,-1/
> +      DATA (NHEL(I,  20),I=1,5) /-1,-1,-1,-1, 1/
> +      DATA (NHEL(I,  21),I=1,5) /-1,-1, 1, 1,-1/
> +      DATA (NHEL(I,  22),I=1,5) /-1,-1, 1, 1, 1/
> +      DATA (NHEL(I,  23),I=1,5) /-1,-1, 1,-1,-1/
> +      DATA (NHEL(I,  24),I=1,5) /-1,-1, 1,-1, 1/
> +      DATA (NHEL(I,  25),I=1,5) /-1, 1,-1, 1,-1/
> +      DATA (NHEL(I,  26),I=1,5) /-1, 1,-1, 1, 1/
> +      DATA (NHEL(I,  27),I=1,5) /-1, 1,-1,-1,-1/
> +      DATA (NHEL(I,  28),I=1,5) /-1, 1,-1,-1, 1/
> +      DATA (NHEL(I,  29),I=1,5) /-1, 1, 1, 1,-1/
> +      DATA (NHEL(I,  30),I=1,5) /-1, 1, 1, 1, 1/
> +      DATA (NHEL(I,  31),I=1,5) /-1, 1, 1,-1,-1/
> +      DATA (NHEL(I,  32),I=1,5) /-1, 1, 1,-1, 1/
> +      LOGICAL GOODHEL(NCOMB)
> +      DATA GOODHEL/NCOMB*.FALSE./
> +      INTEGER NTRY
> +      DATA NTRY/0/
> +      DATA IDEN/ 6/
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      NTRY=NTRY+1
> +      ANS = 0D0
> +      DO IHEL=1,NCOMB
> +        IF (GOODHEL(IHEL) .OR. NTRY .LT. 2) THEN
> +          IF (NTRY.LT.2) THEN
> +C           for the first ps-point, check for helicities that give
> +C           identical matrix elements
> +            T=MATRIX_1(P ,NHEL(1,IHEL))
> +            T_SAVE(IHEL)=T
> +            T_IDENT(IHEL)=-1
> +            DO I=1,IHEL-1
> +              IF (T.EQ.0D0) EXIT
> +              IF (T_SAVE(I).EQ.0D0) CYCLE
> +              IF (ABS(T/T_SAVE(I)-1D0) .LT. 1D-12) THEN
> +C               WRITE (*,*) 'FOUND IDENTICAL',T,IHEL,T_SAVE(I),I
> +                T_IDENT(IHEL) = I
> +              ENDIF
> +            ENDDO
> +          ELSE
> +            IF (T_IDENT(IHEL).GT.0) THEN
> +C             if two helicity states are identical, dont recompute
> +              T=T_SAVE(T_IDENT(IHEL))
> +              T_SAVE(IHEL)=T
> +            ELSE
> +              T=MATRIX_1(P ,NHEL(1,IHEL))
> +              T_SAVE(IHEL)=T
> +            ENDIF
> +          ENDIF
> +C         add to the sum of helicities
> +          ANS=ANS+T
> +          IF (T .NE. 0D0 .AND. .NOT. GOODHEL(IHEL)) THEN
> +            GOODHEL(IHEL)=.TRUE.
> +          ENDIF
> +        ENDIF
> +      ENDDO
> +      ANS=ANS/DBLE(IDEN)
> +      END
> +
> +
> +      REAL*8 FUNCTION MATRIX_1(P,NHEL)
> +C     
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     
> +C     Returns amplitude squared summed/avg over colors
> +C     for the point with external lines W(0:6,NEXTERNAL)
> +C     
> +C     Process: t > b u d~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     Process: t > b c s~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS
> +C     
> +      INTEGER    NGRAPHS
> +      PARAMETER (NGRAPHS=4)
> +      INTEGER    NWAVEFUNCS, NCOLOR
> +      PARAMETER (NWAVEFUNCS=7, NCOLOR=2)
> +      REAL*8     ZERO
> +      PARAMETER (ZERO=0D0)
> +      COMPLEX*16 IMAG1
> +      PARAMETER (IMAG1=(0D0,1D0))
> +      INCLUDE 'nexternal.inc'
> +      INCLUDE 'coupl.inc'
> +C     
> +C     ARGUMENTS 
> +C     
> +      REAL*8 P(0:3,NEXTERNAL)
> +      INTEGER NHEL(NEXTERNAL)
> +C     
> +C     LOCAL VARIABLES 
> +C     
> +      INTEGER I,J
> +      INTEGER IC(NEXTERNAL)
> +      DATA IC /NEXTERNAL*1/
> +      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
> +      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS)
> +C     
> +C     COLOR DATA
> +C     
> +      DATA DENOM(1)/1/
> +      DATA (CF(I,  1),I=  1,  2) /   12,    0/
> +C     1 T(2,1) T(5,3,4)
> +      DATA DENOM(2)/1/
> +      DATA (CF(I,  2),I=  1,  2) /    0,   12/
> +C     1 T(3,4) T(5,2,1)
> +C     ----------
> +C     BEGIN CODE
> +C     ----------
> +      CALL IXXXXX(P(0,1),MDL_MT,NHEL(1),+1*IC(1),W(1,1))
> +      CALL OXXXXX(P(0,2),MDL_MB,NHEL(2),+1*IC(2),W(1,2))
> +      CALL OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
> +      CALL IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
> +      CALL VXXXXX(P(0,5),ZERO,NHEL(5),+1*IC(5),W(1,5))
> +      CALL FFV1_1(W(1,2),W(1,5),GC_11,MDL_MB,ZERO,W(1,6))
> +      CALL FFV2_3(W(1,4),W(1,3),GC_100,MDL_MW,MDL_WW,W(1,7))
> +C     Amplitude(s) for diagram number 1
> +      CALL FFV2_0(W(1,1),W(1,6),W(1,7),GC_100,AMP(1))
> +      CALL FFV2_1(W(1,2),W(1,7),GC_100,MDL_MT,MDL_WT,W(1,6))
> +C     Amplitude(s) for diagram number 2
> +      CALL FFV1_0(W(1,1),W(1,6),W(1,5),GC_11,AMP(2))
> +      CALL FFV1_1(W(1,3),W(1,5),GC_11,ZERO,ZERO,W(1,6))
> +      CALL FFV2_3(W(1,4),W(1,6),GC_100,MDL_MW,MDL_WW,W(1,7))
> +C     Amplitude(s) for diagram number 3
> +      CALL FFV2_0(W(1,1),W(1,2),W(1,7),GC_100,AMP(3))
> +      CALL FFV1_2(W(1,4),W(1,5),GC_11,ZERO,ZERO,W(1,7))
> +      CALL FFV2_3(W(1,7),W(1,3),GC_100,MDL_MW,MDL_WW,W(1,5))
> +C     Amplitude(s) for diagram number 4
> +      CALL FFV2_0(W(1,1),W(1,2),W(1,5),GC_100,AMP(4))
> +      JAMP(1)=-AMP(3)-AMP(4)
> +      JAMP(2)=-AMP(1)-AMP(2)
> +      MATRIX_1 = 0.D0
> +      DO I = 1, NCOLOR
> +        ZTEMP = (0.D0,0.D0)
> +        DO J = 1, NCOLOR
> +          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
> +        ENDDO
> +        MATRIX_1 = MATRIX_1+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
> +      ENDDO
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nFKSconfigs.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nFKSconfigs.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nFKSconfigs.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,4 @@
> +      INTEGER FKS_CONFIGS
> +      PARAMETER (FKS_CONFIGS=3)
> +
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ncombs.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ncombs.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ncombs.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,2 @@
> +      INTEGER    N_MAX_CL
> +      PARAMETER (N_MAX_CL=32)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nexternal.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nexternal.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%nexternal.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,4 @@
> +      INTEGER    NEXTERNAL
> +      PARAMETER (NEXTERNAL=5)
> +      INTEGER    NINCOMING
> +      PARAMETER (NINCOMING=1)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ngraphs.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ngraphs.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%ngraphs.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,2 @@
> +      INTEGER    N_MAX_CG
> +      PARAMETER (N_MAX_CG=4)
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_1.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_1.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_1.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,98 @@
> +      SUBROUTINE DLUM_1(LUM)
> +C     ****************************************************            
> +C         
> +C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
> +C     By the MadGraph5_aMC@NLO Development Team
> +C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
> +C     RETURNS PARTON LUMINOSITIES FOR MADFKS                          
> +C        
> +C     
> +C     Process: t > b u d~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     Process: t > b c s~ g WEIGHTED=5 QED=2 [ real = QCD ]
> +C     
> +C     ****************************************************            
> +C         
> +      IMPLICIT NONE
> +C     
> +C     CONSTANTS                                                       
> +C         
> +C     
> +      INCLUDE 'genps.inc'
> +      INCLUDE 'nexternal.inc'
> +      DOUBLE PRECISION       CONV
> +      PARAMETER (CONV=389379660D0)  !CONV TO PICOBARNS             
> +C     
> +C     ARGUMENTS                                                       
> +C         
> +C     
> +      DOUBLE PRECISION PP(0:3,NEXTERNAL), LUM
> +C     
> +C     LOCAL VARIABLES                                                 
> +C         
> +C     
> +      INTEGER I, ICROSS,ITYPE,LP
> +      DOUBLE PRECISION P1(0:3,NEXTERNAL)
> +
> +      DOUBLE PRECISION XPQ(-7:7)
> +C     
> +C     EXTERNAL FUNCTIONS                                              
> +C         
> +C     
> +      DOUBLE PRECISION ALPHAS2,REWGT,PDG2PDF
> +C     
> +C     GLOBAL VARIABLES                                                
> +C         
> +C     
> +      INTEGER              IPROC
> +      DOUBLE PRECISION PD(0:MAXPROC)
> +      COMMON /SUBPROC/ PD, IPROC
> +      INCLUDE 'coupl.inc'
> +      INCLUDE 'run.inc'
> +      INTEGER IMIRROR
> +      COMMON/CMIRROR/IMIRROR
> +C     
> +C     DATA                                                            
> +C         
> +C     
> +
> +      DATA ICROSS/1/
> +C     ----------                                                      
> +C         
> +C     BEGIN CODE                                                      
> +C         
> +C     ----------                                                      
> +C         
> +      LUM = 0D0
> +      IF (IMIRROR.EQ.2) THEN
> +        PD(0) = 0D0
> +        IPROC = 0
> +        IPROC=IPROC+1  ! t > b u d~ g
> +        PD(IPROC) = 1D0
> +
> +        PD(0)=PD(0)+PD(IPROC)
> +        IPROC=IPROC+1  ! t > b c s~ g
> +        PD(IPROC) = 1D0
> +
> +        PD(0)=PD(0)+PD(IPROC)
> +      ELSE
> +        PD(0) = 0D0
> +        IPROC = 0
> +        IPROC=IPROC+1  ! t > b u d~ g
> +        PD(IPROC) = 1D0
> +
> +        PD(0)=PD(0)+PD(IPROC)
> +        IPROC=IPROC+1  ! t > b c s~ g
> +        PD(IPROC) = 1D0
> +
> +        PD(0)=PD(0)+PD(IPROC)
> +      ENDIF
> +      DO I=1,IPROC
> +        IF (NINCOMING.EQ.2) THEN
> +          LUM = LUM + PD(I) * CONV
> +        ELSE
> +          LUM = LUM + PD(I)
> +        ENDIF
> +      ENDDO
> +      RETURN
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_chooser.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_chooser.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%parton_lum_chooser.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,21 @@
> +      DOUBLE PRECISION FUNCTION DLUM()
> +      IMPLICIT NONE
> +      INCLUDE 'timing_variables.inc'
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +      CALL CPU_TIME(TBEFORE)
> +      IF (NFKSPROCESS.EQ.1) THEN
> +        CALL DLUM_1(DLUM)
> +      ELSEIF (NFKSPROCESS.EQ.2) THEN
> +        CALL DLUM_1(DLUM)
> +      ELSEIF (NFKSPROCESS.EQ.3) THEN
> +        CALL DLUM_1(DLUM)
> +      ELSE
> +        WRITE(*,*) 'ERROR: invalid n in dlum :', NFKSPROCESS
> +        STOP
> +      ENDIF
> +      CALL CPU_TIME(TAFTER)
> +      TPDF = TPDF + (TAFTER-TBEFORE)
> +      RETURN
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%pmass.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%pmass.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%pmass.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,5 @@
> +      PMASS(1)=ABS(MDL_MT)
> +      PMASS(2)=ABS(MDL_MB)
> +      PMASS(3)=ZERO
> +      PMASS(4)=ZERO
> +      PMASS(5)=ZERO
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_from_born_configs.inc'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_from_born_configs.inc	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_from_born_configs.inc	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,5 @@
> +      INTEGER IRFBC
> +      INTEGER REAL_FROM_BORN_CONF(1,3)
> +      DATA (REAL_FROM_BORN_CONF(IRFBC,1),IRFBC=1,1) /1/
> +      DATA (REAL_FROM_BORN_CONF(IRFBC,2),IRFBC=1,1) /3/
> +      DATA (REAL_FROM_BORN_CONF(IRFBC,3),IRFBC=1,1) /4/
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_me_chooser.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_me_chooser.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%real_me_chooser.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,20 @@
> +      SUBROUTINE SMATRIX_REAL(P, WGT)
> +      IMPLICIT NONE
> +      INCLUDE 'nexternal.inc'
> +      DOUBLE PRECISION P(0:3, NEXTERNAL)
> +      DOUBLE PRECISION WGT
> +      INTEGER NFKSPROCESS
> +      COMMON/C_NFKSPROCESS/NFKSPROCESS
> +      IF (NFKSPROCESS.EQ.1) THEN
> +        CALL SMATRIX_1(P, WGT)
> +      ELSEIF (NFKSPROCESS.EQ.2) THEN
> +        CALL SMATRIX_1(P, WGT)
> +      ELSEIF (NFKSPROCESS.EQ.3) THEN
> +        CALL SMATRIX_1(P, WGT)
> +      ELSE
> +        WRITE(*,*) 'ERROR: invalid n in real_matrix :', NFKSPROCESS
> +        STOP
> +      ENDIF
> +      RETURN
> +      END
> +
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,45 @@
> +      SUBROUTINE SBORN_SF(P_BORN,M,N,WGT)
> +      IMPLICIT NONE
> +      INCLUDE 'nexternal.inc'
> +      DOUBLE PRECISION P_BORN(0:3,NEXTERNAL-1),WGT
> +      DOUBLE COMPLEX WGT1(2)
> +      INTEGER M,N
> +
> +C     b_sf_001 links partons 1 and 1 
> +      IF (M.EQ.1 .AND. N.EQ.1) THEN
> +        CALL SB_SF_001(P_BORN,WGT)
> +
> +C       b_sf_002 links partons 1 and 2 
> +      ELSEIF ((M.EQ.1 .AND. N.EQ.2).OR.(M.EQ.2 .AND. N.EQ.1)) THEN
> +        CALL SB_SF_002(P_BORN,WGT)
> +
> +C       b_sf_003 links partons 1 and 3 
> +      ELSEIF ((M.EQ.1 .AND. N.EQ.3).OR.(M.EQ.3 .AND. N.EQ.1)) THEN
> +        CALL SB_SF_003(P_BORN,WGT)
> +
> +C       b_sf_004 links partons 1 and 4 
> +      ELSEIF ((M.EQ.1 .AND. N.EQ.4).OR.(M.EQ.4 .AND. N.EQ.1)) THEN
> +        CALL SB_SF_004(P_BORN,WGT)
> +
> +C       b_sf_005 links partons 2 and 2 
> +      ELSEIF (M.EQ.2 .AND. N.EQ.2) THEN
> +        CALL SB_SF_005(P_BORN,WGT)
> +
> +C       b_sf_006 links partons 2 and 3 
> +      ELSEIF ((M.EQ.2 .AND. N.EQ.3).OR.(M.EQ.3 .AND. N.EQ.2)) THEN
> +        CALL SB_SF_006(P_BORN,WGT)
> +
> +C       b_sf_007 links partons 2 and 4 
> +      ELSEIF ((M.EQ.2 .AND. N.EQ.4).OR.(M.EQ.4 .AND. N.EQ.2)) THEN
> +        CALL SB_SF_007(P_BORN,WGT)
> +
> +C       b_sf_008 links partons 3 and 4 
> +      ELSEIF ((M.EQ.3 .AND. N.EQ.4).OR.(M.EQ.4 .AND. N.EQ.3)) THEN
> +        CALL SB_SF_008(P_BORN,WGT)
> +
> +      ELSE
> +        WGT = 0D0
> +      ENDIF
> +
> +      RETURN
> +      END
> 
> === added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf_dum.f'
> --- tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf_dum.f	1970-01-01 00:00:00 +0000
> +++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_tdecay_fksreal/%SubProcesses%P0_t_budx%sborn_sf_dum.f	2015-10-08 17:23:54 +0000
> @@ -0,0 +1,13 @@
> +      SUBROUTINE SBORN_SF(P_BORN,M,N,WGT)
> +      IMPLICIT NONE
> +      INCLUDE 'nexternal.inc'
> +      DOUBLE PRECISION P_BORN(0:3,NEXTERNAL-1),WGT
> +      DOUBLE COMPLEX WGT1(2)
> +      INTEGER M,N
> +
> +C     This is a dummy function because
> +C     this subdir has no soft singularities
> +      WGT = 0D0
> +
> +      RETURN
> +      END
> 
> === modified file 'tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2014-03-05 09:04:59 +0000
> +++ tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/%..%..%Source%MODEL%model_functions.inc'
> --- tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2014-03-19 22:06:11 +0000
> +++ tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2015-10-08 17:23:54 +0000
> @@ -6,7 +6,11 @@
>       DOUBLE COMPLEX CONDIF
>       DOUBLE COMPLEX REGLOG
>       DOUBLE COMPLEX ARG
> +
> +
>       COMPLEX*32 MP_COND
>       COMPLEX*32 MP_CONDIF
>       COMPLEX*32 MP_REGLOG
>       COMPLEX*32 MP_ARG
> +
> +
> 
> === modified file 'tests/input_files/LoopSMEWTest/object_library.py'
> --- tests/input_files/LoopSMEWTest/object_library.py	2015-10-01 16:00:08 +0000
> +++ tests/input_files/LoopSMEWTest/object_library.py	2015-10-08 17:23:54 +0000
> @@ -294,7 +294,7 @@
>                     #else:
>                 def substitution(matchedObj):
>                     return matchedObj.group('first')+"("+ctpar.pole(x)+")"+matchedObj.group('second')
> -                       pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\(|\s)(?P<name>"+ctpar.name+r")(?P<second>\Z|\*|\+|\-|\)|/|\\|\s)")
> +                pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\(|\s)(?P<name>"+ctpar.name+r")(?P<second>\Z|\*|\+|\-|\)|/|\\|\s)")
>                 tempvalue2 = tempvalue
>                 tempvalue = pattern.sub(substitution,tempvalue)
>                 num = numbermatch[i]-1		
> 
> === modified file 'tests/unit_tests/iolibs/test_export_fks.py'
> --- tests/unit_tests/iolibs/test_export_fks.py	2015-10-01 16:00:08 +0000
> +++ tests/unit_tests/iolibs/test_export_fks.py	2015-10-08 17:23:54 +0000
> @@ -68,6 +68,10 @@
>         """ target: SubProcesses/[P0.*\/.+\.(inc|f)]"""
>         self.generate(['p p > t t~ [real=QCD]'], 'sm')
> 
> +    @IOTests.createIOTest()
> +    def testIO_test_tdecay_fksreal(self):
> +        """ target: SubProcesses/[P0.*\/.+\.(inc|f)]"""
> +        self.generate(['t > j j b [real=QCD]'], 'sm')
> 
>     @IOTests.createIOTest()
>     def testIO_test_pptt_fks_loonly(self):
>