Merge lp:~maddevelopers/mg5amcnlo/NLO_decay_process into lp:~maddevelopers/mg5amcnlo/2.3.3

[Rikkert Frederix]

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.

[Olivier Mattelaer (olivier-mattelaer)]

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 su...

[Rikkert Frederix (frederix)]

Hi Olivier,

That would be great. It's just the 'ptj', 'mll_sf' and 'ptgmin' that should be set to zero.

Thanks,
Rikkert

[Olivier Mattelaer (olivier-mattelaer)]

Hi Rik,

I put to “no cut” all the run_card parameter flagged as “cut” in the python object. (The same will apply to all future cut)
This is done at the default run_card level, so all the cut are still allowed to be modified but are not present by default.
This is the exact behaviour that we have @LO.

When trying to test, I face this error:
aMCatNLOError : Decay processes can only be run at fixed order.

So I have also change the default question to:

The following switches determine which operations are executed:
 1 Perturbative order of the calculation: order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=ON
 3 Shower the generated events: shower=Not available for decay
 4 Decay particles with the MadSpin module: madspin=Not available for decay
 5 Add weights to the events based on changing model parameters: reweight=Not available for decay
  Either type the switch number (1 to 5) to change its default setting,
  or set any switch explicitly (e.g. type 'order=LO' at the prompt)
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, 1, auto, done, order=LO, order=NLO, fixed_order=ON, LO, NLO][60s to answer]

Such that this is automatically a good setup which is used and all unsupported case are forbidden.

For both of them, I had to change the way aMC@NLO supports the proc_characteristic object/file.
So Marco please take a look.

Cheers,

Olivier

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

> Hi Olivier,
>
> That would be great. It's just the 'ptj', 'mll_sf' and 'ptgmin' that should be set to zero.
>
> Thanks,
> Rikkert
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/NLO_decay_process/+merge/273880
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/mg5amcnlo/2.3.3.

[lp:~maddevelopers/mg5amcnlo/NLO_decay_process]

317. By Olivier Mattelaer on 2015-10-09

improve the interface

[Rikkert Frederix (frederix)]

Hi Olivier,

Thanks.

One question though (maybe this was already there in the default version; I haven't checked). What is this:

"fixed_order=ON, LO, NLO"

Shouldn't it not only be "ON" (or "OFF")?

Cheers,
Rik

[Hua-Sheng Shao (erdissshaw)]

HI Rik,

It is really cool that decay process works at NLO. I quickly tested one process in order to compare with SMWidth. If one wants to calculate the total width for Z (which can be used by MadWidth for CMS@NLO for example), I

define f = j l+ l- vl vl~ ta+ ta-
generate z > f f [QCD]

It is apparent that there is no QCD corrections to some subprocesses like z > l+ l-. The program crashes since no virtual and real diagrams will be generated.

Is it possible to automatically realize if there is no QCD corrections and only calculates tree-level? Or you think it is better to do it in MadWidth to decide which subprocesses have QCD corrections.

A tiny thing: in summary.txt and also the output, it is

      Final results and run summary:
      Process z > j j [QCD]
      Run at l-l collider (6500.0 + 6500.0 GeV)
      (Partial) decay width: 1.761e+00 +- 2.0e-03 GeV
      Ren. and fac. scale uncertainty: +0.4% -0.3%

Can you remove "Run at l-l collider (6500.0 + 6500.0 GeV)". Moreover, if one wants to change the scale, how can one set in decay process?

@Marco: do you think how easy/difficult to implement Rik's implementation in EW branches ? It might be quite useful and have some complicated comparison between MG5 and SMWidth.

I will give a serious report later and compare the numbers with SMWidth.

Great work and awesome !

Cheers,

Hua-Sheng

On Oct 9, 2015, at 11:33 AM, Olivier Mattelaer <email address hidden> wrote:

> Hi Rik,
>
> I put to “no cut” all the run_card parameter flagged as “cut” in the python object. (The same will apply to all future cut)
> This is done at the default run_card level, so all the cut are still allowed to be modified but are not present by default.
> This is the exact behaviour that we have @LO.
>
> When trying to test, I face this error:
> aMCatNLOError : Decay processes can only be run at fixed order.
>
> So I have also change the default question to:
>
> The following switches determine which operations are executed:
> 1 Perturbative order of the calculation: order=NLO
> 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=ON
> 3 Shower the generated events: shower=Not available for decay
> 4 Decay particles with the MadSpin module: madspin=Not available for decay
> 5 Add weights to the events based on changing model parameters: reweight=Not available for decay
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'order=LO' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, auto, done, order=LO, order=NLO, fixed_order=ON, LO, NLO][60s to answer]
>
> Such that this is automatically a good setup which is used and all unsupported case are forbidden.
>
> For both of them, I had to change the way aMC@NLO supports the proc_characteristic object/file.
> So Marco please take a look.
>
> Cheers,
>
> Olivier
>
>
> On 08 Oct 2015, at 19:42, Rikkert Frederix <email address hidden> wrote:
>
>> Hi Olivier,
>>
>> That would be great. It's just the 'ptj', 'mll_sf' and 'ptgm...

[marco zaro (marco-zaro)]

Hi,
is fixed-order supposed to always be on?
I cannot toggle it off, neither entering 2, nor with “fixed_order=OFF”
@Hua-sheng, it should be rather easy to merge with the EW branch

Cheers,

Marco

The following switches determine which operations are executed:
 1 Perturbative order of the calculation: order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=ON
 3 Shower the generated events: shower=Not available for decay
 4 Decay particles with the MadSpin module: madspin=Not available for decay
 5 Add weights to the events based on changing model parameters: reweight=Not available for decay
  Either type the switch number (1 to 5) to change its default setting,
  or set any switch explicitly (e.g. type 'order=LO' at the prompt)
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, 1, auto, done, order=LO, order=NLO, fixed_order=ON, LO, NLO][60s to answer]
>2
2 not valid argument. Valid argument are in (0,1,auto,done,order=LO,order=NLO,fixed_order=ON,LO,NLO,fixed_order=on,order=nlo,order=lo).
please retry
>[timer stopped]
fixed_order=ON
INFO: For coherence 'shower' is set to 'OFF'
INFO: For coherence 'madspin' is set to 'OFF'
INFO: For coherence 'reweight' is set to 'OFF'
The following switches determine which operations are executed:
 1 Perturbative order of the calculation: order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=ON
 3 Shower the generated events: shower=OFF
 4 Decay particles with the MadSpin module: madspin=OFF
 5 Add weights to the events based on changing model parameters: reweight=OFF
  Either type the switch number (1 to 5) to change its default setting,
  or set any switch explicitly (e.g. type 'order=LO' at the prompt)
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, 1, auto, done, order=LO, order=NLO, fixed_order=ON, LO, NLO, ... ][60s to answer]
>fixed_order=OFF
fixed_order=OFF not valid argument. Valid argument are in (0,1,auto,done,order=LO,order=NLO,fixed_order=ON,LO,NLO,fixed_order=on,order=nlo,order=lo,fixed_order=on,order=nlo,order=lo).
please retry
>
use 0

On 09 Oct 2015, at 11:48, Hua-Sheng Shao <email address hidden> wrote:

> HI Rik,
>
> It is really cool that decay process works at NLO. I quickly tested one process in order to compare with SMWidth. If one wants to calculate the total width for Z (which can be used by MadWidth for CMS@NLO for example), I
>
> define f = j l+ l- vl vl~ ta+ ta-
> generate z > f f [QCD]
>
> It is apparent that there is no QCD corrections to some subprocesses like z > l+ l-. The program crashes since no virtual and real diagrams will be generated.
>
> Is it possible to automatically realize if there is no QCD corrections and only calculates tree-level? Or you think it is better to do it in MadWidth to decide which subprocesses have QCD corrections.
>
> A tiny thing: in summary.txt and also the output, it is
>
> Final results and run summary:
> Process z...

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

This is the list of all available valid input
you have three valid entry that you can type:
> fixed_order=ON
LO
NLO
This is NOT suggesting “fixed_order=LO”

But ok having fixed_order=ON is irrelevant, So I have remove it:

The following switches determine which operations are executed:
 1 Perturbative order of the calculation: order=NLO
 2 Fixed order (no event generation and no MC@[N]LO matching): fixed_order=ON
 3 Shower the generated events: shower=Not available for decay
 4 Decay particles with the MadSpin module: madspin=Not available for decay
 5 Add weights to the events based on changing model parameters: reweight=Not available for decay
  Either type the switch number (1 to 5) to change its default setting,
  or set any switch explicitly (e.g. type 'order=LO' at the prompt)
  Type '0', 'auto', 'done' or just press enter when you are done.
 [0, 1, auto, done, order=LO, order=NLO, LO, NLO][60s to answer]

Cheers,

Olivier
On 09 Oct 2015, at 11:42, Rikkert Frederix <email address hidden> wrote:

> Hi Olivier,
>
> Thanks.
>
> One question though (maybe this was already there in the default version; I haven't checked). What is this:
>
> "fixed_order=ON, LO, NLO"
>
> Shouldn't it not only be "ON" (or "OFF")?
>
> Cheers,
> Rik
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/NLO_decay_process/+merge/273880
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/mg5amcnlo/2.3.3.

[lp:~maddevelopers/mg5amcnlo/NLO_decay_process]

318. By Olivier Mattelaer on 2015-10-09

slight improvement

319. By Rikkert Frederix on 2015-10-09

removed the printing of the collider (energy) at the end of the run for decay processes

[Rikkert Frederix (frederix)]

Hi Hua-Sheng,

I don't know how to combine LO and NLO processes like that. I guess it would be some work to implement it. Right now, I think it's okay to have it only for QCD processes, just like collider processes are only for QCD processes (things like "e+ e- > mu+ mu- [ all = QCD ]") don't work).

I've removed the "Run at l-l collider (6500.0 + 6500.0 GeV)" stuff from the run summary.

Marco: Yeah, it is only supposed to work for fixed order. I don't think it makes so much sense to do event generation for decays.

Cheers,
Rik

[Olivier Mattelaer (olivier-mattelaer)]

Hi guys,

Could we finish this review soon such that we can release 2.3.3?

I take a look and it works nicely. Great work.
Here is a couple of superficial comment/question:
(I did not look at the code since I'm not a reviewer)

1) I tried:
generate t > b all all [QCD]

and I have the following warning:
WARNING: Process: t > b g/gh/gh~/d/u/s/c/d~/u~/s~/c~/a/ve/vm/vt/e-/mu-/ve~/vm~/vt~/e+/mu+/b/t/b~/t~/z/w+/h/w-/ta-/ta+ g/gh/gh~/d/u/s/c/d~/u~/s~/c~/a/ve/vm/vt/e-/mu-/ve~/vm~/vt~/e+/mu+/b/t/b~/t~/z/w+/h/w-/ta-/ta+ [ all = QCD ] can have real emission processes which are not finite.
To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
 g, c, s, u, d, d~, u~, s~, c~

This should not be raised right?

Finally it crashed with:

Command "launch auto " interrupted with error:
IOError : [Errno 2] No such file or directory: '/Users/omatt/Documents/eclipse/NLO_decay_process/topwidth/SubProcesses/P0_t_tbbx/channels.txt'
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/Users/omatt/Documents/eclipse/NLO_decay_process/topwidth/run_01_tag_1_debug.log'.
Please attach this file to your report.

This should be very simple to prevent.

2) I tried
generate w+ > all all [QCD]
and it also crashes:

Command "generate w+ > all all [QCD]" interrupted with error:
InvalidCmd : All loop diagrams discarded by user selection.
 Consider using a tree-level generation or relaxing the coupling order constraints.
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in 'MG5_debug'.
Please attach this file to your report.

3) I tried
generate t > w+ b [QCD]
then I see as output:

(Partial) decay width: 1.347e+00 +- 8.9e-03 GeV
Ren. and fac. scale uncertainty: +1.2% -1.5%

One small remark here, "factorization variation" is irrelevant here and can be remove from the log [but ok we can keep it do not hurt]

then I re-run with fix (renormalisation) scale computation (at 90GeV) and get
(Partial) decay width: 1.343e+00 +- 9.2e-03 GeV
Clearly my scale choice is not consider.

This is confusion (at minima) and should either works or at least raise a warning that we do not support that (but should be easy to handle).
I would propose that the default "dynamical" scale is the mass of the decaying particle.
and that at least static scale works.

4) [This should not be implemented now but I want to start the discussion].
For the moment the "only" output is a partial width/plot. It will be useful to have a param_card as output with updated information (that's typically the point of computing the width).
Concerning the BR, this is obviously less clear @NLO, but I would propose that we use the "proton"
pdg code for the branching ratio to avoid problem with the massless quark flavor. (this can be discussed with pythia/herwig author)

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/NLO_decay_process]

320. By Rikkert Frederix on 2015-10-16

address Olivier's suggestions and made sure that the width of initial
state particles is also set to zero for consistency.

[Rikkert Frederix (frederix)]

Hi Olivier,

Thanks for trying it out!

Let me address your questions.

1). About the multi-particle warning, it should indeed not be raised. However, this is not specific to this branch, and I don't know how these multi-particles labels are defined. Maybe Marco can have a quick look to this?
I "fixed" the crash with a WARNING written to screen. Note that computing "t > b all all" does not make sense at all from a physics point of view because you are double counting many contributions (for example "t > b e+ ve" with "t > b wp g").

2). Remember that QCD NLO corrections can only be computed to processes that have at least one QCD particle in them. Therefore "w+ > all all" cannot be computed with the current code. This is not a limitation of this branch, but much more general. I don't think this branch is the right place to address this problem.

3). This was already exactly the way you want it to be. Note that the renormalisation scale dependence is simply rather small (your integration uncertainty is about 1% and the scale uncertainty is also about 1%).

4). I agree that this would be useful at some point. But I think first the above points 1) and 2) should be addressed before we could even start thinking about this.

Cheers,
Rik

[Olivier Mattelaer (olivier-mattelaer)]

Thanks Rik,

2). Remember that QCD NLO corrections can only be computed to processes that have at least one QCD particle in them. Therefore "w+ > all all" cannot be computed with the current code. This is not a limitation of this branch, but much more general. I don't think this branch is the right place to address this problem.

Yes we can do that in a second stage.

3). This was already exactly the way you want it to be. Note that the renormalisation scale dependence is simply rather small (your integration uncertainty is about 1% and the scale uncertainty is also about 1%).

Ah ok did not look at the statistical error. Thanks.
Just for curiosity. Did you use MC technique for the integration or are you using simpson? (so much faster in this case)

Cheers and thanks,

Olivier

review: Approve

[Rikkert Frederix (frederix)]

Hi Olivier,

> Ah ok did not look at the statistical error. Thanks.
> Just for curiosity. Did you use MC technique for the integration or are you
> using simpson? (so much faster in this case)

Normal MC integration. Note that for NLO processes, you always have 3 more integration dimensions compared to LO, so I doubt that Simpson is much faster. In particular for anything beyond 1->2 decays.

Cheers,
Rikkert

[Hua-Sheng Shao (erdissshaw)]

Hi Rik,

I have looked at the changes. It looks fine. Moreover, I have checked the total width (LO and QCD corrections) in loop_sm, i.e. w, z, top. It looks same as SMWidth. Moroever, for a 1 > 3 process t > j j j [QCD]. It also works. Hence, I think it is good to merge into trunk.

Cheers,

Hua-Sheng

review: Approve

[lp:~maddevelopers/mg5amcnlo/NLO_decay_process]

321. By Rikkert Frederix on 2015-10-16

merge with latest 2.3.3

322. By marco zaro on 2015-10-20

the warning for multiparticle definition is not raised for decay processes

[marco zaro (marco-zaro)]

Hi Rik
sorry for the delay, somehow i forgot about this.
So i had a look at the changes in the fortran code and in the fks python part, everything looks great there.
Green light from me.
For the warning on multiparticle labels, it is not there any more for decay processes.

Cheers,

Marco

review: Approve

=== 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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-15 20:22:39 +0000
+++ UpdateNotes.txt	2015-10-20 10:56:39 +0000
@@ -11,6 +11,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).
         OM: (SysCalc interface) Allow to bypass the pdf reweighting/alpsfact reweighting
         MZ: fixed bug related to slurm clusters
 	OM: remove the addmasses.py script of running by default on gridpack mode. 

=== 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-20 10:56:39 +0000
@@ -132,6 +132,9 @@
         perturbation = []
         for procdef in self['process_definitions']:
             soft_particles = []
+            # do not warn for decay processes
+            if [ i['state'] for i in procdef['legs']].count(False) == 1:
+                continue
             for pert in procdef['perturbation_couplings']:
                 if pert not in perturbation:
                     perturbation.append(pert)
@@ -550,10 +553,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-11 19:18:37 +0000
+++ madgraph/interface/amcatnlo_interface.py	2015-10-20 10:56:39 +0000
@@ -543,8 +543,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 = []
@@ -570,10 +575,10 @@
 
             #_curr_matrix_element is a FKSHelasMultiProcess Object 
             self._fks_directories = []
-            proc_charac = banner_mod.ProcCharacteristic()
+            proc_charac = self._curr_exporter.proc_characteristic 
             for charac in ['has_isr', 'has_fsr', 'has_loops']:
                 proc_charac[charac] = self._curr_matrix_elements[charac]
-            proc_charac.write(pjoin(path, 'proc_characteristics'))
+
 
             for ime, me in \
                 enumerate(self._curr_matrix_elements.get('matrix_elements')):

=== modified file 'madgraph/interface/amcatnlo_run_interface.py'
--- madgraph/interface/amcatnlo_run_interface.py	2015-10-16 14:52:37 +0000
+++ madgraph/interface/amcatnlo_run_interface.py	2015-10-20 10:56:39 +0000
@@ -1364,6 +1364,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()
@@ -1448,8 +1450,12 @@
             npoints = self.run_card['npoints_FO_grid']
             niters = self.run_card['niters_FO_grid']
             for p_dir in p_dirs:
-                with open(pjoin(self.me_dir,'SubProcesses',p_dir,'channels.txt')) as chan_file:
-                    channels=chan_file.readline().split()
+                try:
+                    with open(pjoin(self.me_dir,'SubProcesses',p_dir,'channels.txt')) as chan_file:
+                        channels=chan_file.readline().split()
+                except IOError:
+                    logger.warning('No integration channels found for contribution %s' % p_dir)
+                    continue
                 for channel in channels:
                     job={}
                     job['p_dir']=p_dir
@@ -2152,10 +2158,21 @@
             if line.startswith('generate') or line.startswith('add process'):
                 process = process+(line.replace('generate ', '')).replace('add process ','')+' ; '
         lpp = {0:'l', 1:'p', -1:'pbar'}
-        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:
+            proc_info = '\n      Process %s' % process[:-3]
+        else:
+            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, \
@@ -2181,13 +2198,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']:
@@ -2216,15 +2233,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:
@@ -4280,11 +4297,17 @@
         else:
             switch.update(dict((k,value) for k,v in switch_default.items() if k not in switch))
         default_switch = ['ON', 'OFF']
+        
+
         allowed_switch_value = {'order': ['LO', 'NLO'],
                                 'fixed_order': default_switch,
                                 'shower': default_switch,
                                 'madspin': default_switch,
                                 'reweight': default_switch}
+
+            
+            
+        
         
         description = {'order':  'Perturbative order of the calculation:',
                        'fixed_order': 'Fixed order (no event generation and no MC@[N]LO matching):',
@@ -4300,22 +4323,35 @@
         special_values = ['LO', 'NLO', 'aMC@NLO', 'aMC@LO', 'noshower', 'noshowerLO']
 
         assign_switch = lambda key, value: switch.__setitem__(key, value if switch[key] != void else void )
-        
+
+        if self.proc_characteristics['ninitial'] == 1:
+            switch['fixed_order'] = 'ON'
+            switch['shower'] = 'Not available for decay'
+            switch['madspin'] = 'Not available for decay'
+            switch['reweight'] = 'Not available for decay'
+            allowed_switch_value['fixed_order'] = ['ON']
+            allowed_switch_value['shower'] = ['OFF']
+            allowed_switch_value['madspin'] = ['OFF']
+            allowed_switch_value['reweight'] = ['OFF']
+            available_mode = ['0','1']
+            special_values = ['LO', 'NLO']
+        else: 
+            # Init the switch value according to the current status
+            available_mode = ['0', '1', '2','3']
 
         if mode == 'auto': 
             mode = None
         if not mode and (options['parton'] or options['reweightonly']):
             mode = 'noshower'         
         
-        # Init the switch value according to the current status
-        available_mode = ['0', '1', '2']
-        available_mode.append('3')
-        if os.path.exists(pjoin(self.me_dir, 'Cards', 'shower_card.dat')):
-            switch['shower'] = 'ON'
-        else:
-            switch['shower'] = 'OFF'
+
+        if '3' in available_mode:
+            if os.path.exists(pjoin(self.me_dir, 'Cards', 'shower_card.dat')):
+                switch['shower'] = 'ON'
+            else:
+                switch['shower'] = 'OFF' 
                 
-        if not aMCatNLO or self.options['mg5_path']:
+        if (not aMCatNLO or self.options['mg5_path']) and '3' in available_mode:
             available_mode.append('4')
             if os.path.exists(pjoin(self.me_dir,'Cards','madspin_card.dat')):
                 switch['madspin'] = 'ON'
@@ -4330,8 +4366,7 @@
             else:
                 switch['reweight'] = 'Not available (requires NumPy)'
 
-        
-        if 'do_reweight' in options and options['do_reweight']:
+        if 'do_reweight' in options and options['do_reweight'] and '3' in available_mode:
             if switch['reweight'] == "OFF":
                 switch['reweight'] = "ON"
             elif switch['reweight'] != "ON":
@@ -4341,12 +4376,12 @@
                 switch['madspin'] = 'ON'
             elif switch['madspin'] != "ON":
                 logger.critical("Cannot run MadSpin module: %s" % switch['reweight'])
-                    
-                    
+                        
         answers = list(available_mode) + ['auto', 'done']
         alias = {}
         for id, key in enumerate(switch_order):
-            if switch[key] != void and switch[key] in allowed_switch_value[key]:
+            if switch[key] != void and switch[key] in allowed_switch_value[key] and \
+                len(allowed_switch_value[key]) >1:
                 answers += ['%s=%s' % (key, s) for s in allowed_switch_value[key]]
                 #allow lower case for on/off
                 alias.update(dict(('%s=%s' % (key, s.lower()), '%s=%s' % (key, s))
@@ -4408,7 +4443,6 @@
                     return
             return switch
 
-
         modify_switch(mode, self.last_mode, switch)
         if switch['madspin'] == 'OFF' and  os.path.exists(pjoin(self.me_dir,'Cards','madspin_card.dat')):
             assign_switch('madspin', 'ON')

=== modified file 'madgraph/interface/common_run_interface.py'
--- madgraph/interface/common_run_interface.py	2015-10-15 10:10:38 +0000
+++ madgraph/interface/common_run_interface.py	2015-10-20 10:56:39 +0000
@@ -658,7 +658,7 @@
 
             if amcatnlo and not keepwidth:
                 # force particle in final states to have zero width
-                pids = self.get_pid_final_states()
+                pids = self.get_pid_final_initial_states()
                 # check those which are charged under qcd
                 if not MADEVENT and pjoin(self.me_dir,'bin','internal') not in sys.path:
                         sys.path.insert(0,pjoin(self.me_dir,'bin','internal'))
@@ -1351,8 +1351,8 @@
         self.update_status('delphes done', level='delphes', makehtml=False)
 
     ############################################################################
-    def get_pid_final_states(self):
-        """Find the pid of all particles in the final states"""
+    def get_pid_final_initial_states(self):
+        """Find the pid of all particles in the final and initial states"""
         pids = set()
         subproc = [l.strip() for l in open(pjoin(self.me_dir,'SubProcesses',
                                                                  'subproc.mg'))]
@@ -1363,7 +1363,7 @@
             group = pat.findall(text)
             for particles in group:
                 particles = particles.split(',')
-                pids.update(set(particles[nb_init:]))
+                pids.update(set(particles))
 
         return pids
 

=== 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-20 10:56:39 +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')
@@ -536,6 +539,8 @@
         (nexternal, ninitial) = matrix_element.get_nexternal_ninitial()
         self.write_nexternal_file(writers.FortranWriter(filename),
                              nexternal, ninitial)
+        self.proc_characteristic['ninitial'] = ninitial
+        self.proc_characteristic['nexternal'] = max(self.proc_characteristic['nexternal'], nexternal)
     
         filename = 'pmass.inc'
         try:
@@ -701,7 +706,8 @@
         pages,proc_card_mg5.dat and madevent.tar.gz."""
         
         self.proc_characteristic['grouped_matrix'] = False
-        
+        self.create_proc_charac()
+
         self.create_run_card(matrix_elements, history)
 #        modelname = self.model.get('name')
 #        if modelname == 'mssm' or modelname.startswith('mssm-'):
@@ -1221,6 +1227,8 @@
     def write_real_matrix_elements(self, matrix_element, fortran_model):
         """writes the matrix_i.f files which contain the real matrix elements""" 
 
+
+
         for n, fksreal in enumerate(matrix_element.real_processes):
             filename = 'matrix_%d.f' % (n + 1)
             self.write_matrix_element_fks(writers.FortranWriter(filename),

=== 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-20 10:56:39 +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 'madgraph/madevent/sum_html.py'
--- madgraph/madevent/sum_html.py	2015-10-01 16:00:08 +0000
+++ madgraph/madevent/sum_html.py	2015-10-20 10:56:39 +0000
@@ -664,24 +664,27 @@
         P_path = pjoin(cmd.me_dir, 'SubProcesses', Pdir)
         G_dir = [G for G in os.listdir(P_path) if G.startswith('G') and 
                                                 os.path.isdir(pjoin(P_path,G))]
-        
-        for line in open(pjoin(P_path, 'symfact.dat')):
-            name, mfactor = line.split()
-            if float(mfactor) < 0:
-                continue
-            if os.path.exists(pjoin(P_path, 'ajob.no_ps.log')):
-                continue
-                                  
-            if not folder_names:
-                name = 'G' + name
-                P_comb.add_results(name, pjoin(P_path,name,'results.dat'), mfactor)
-            else:
-                for folder in folder_names:
-                    if 'G' in folder:
-                        dir = folder.replace('*', name)
-                    else:
-                        dir = folder.replace('*', '_G' + name)
-                    P_comb.add_results(dir, pjoin(P_path,dir,'results.dat'), mfactor)
+
+        try:
+            for line in open(pjoin(P_path, 'symfact.dat')):
+                name, mfactor = line.split()
+                if float(mfactor) < 0:
+                    continue
+                if os.path.exists(pjoin(P_path, 'ajob.no_ps.log')):
+                    continue
+                                      
+                if not folder_names:
+                    name = 'G' + name
+                    P_comb.add_results(name, pjoin(P_path,name,'results.dat'), mfactor)
+                else:
+                    for folder in folder_names:
+                        if 'G' in folder:
+                            dir = folder.replace('*', name)
+                        else:
+                            dir = folder.replace('*', '_G' + name)
+                        P_comb.add_results(dir, pjoin(P_path,dir,'results.dat'), mfactor)
+        except IOError:
+            continue
         P_comb.compute_values()
         all.append(P_comb)
     all.compute_values()

=== modified file 'madgraph/various/banner.py'
--- madgraph/various/banner.py	2015-10-12 12:35:23 +0000
+++ madgraph/various/banner.py	2015-10-20 10:56:39 +0000
@@ -1119,7 +1119,7 @@
          
     def write(self, outputpath):
         """write the file"""
-        
+
         template ="#    Information about the process      #\n"
         template +="#########################################\n"
         
@@ -1430,6 +1430,7 @@
             fsock.writelines(line)
         fsock.close()   
 
+
     def get_banner_init_information(self):
         """return a dictionary with the information needed to write
         the first line of the <init> block of the lhe file."""
@@ -1468,7 +1469,22 @@
         output["pdfsup1"] = get_pdf_id(self["pdlabel"])
         output["pdfsup2"] = get_pdf_id(self["pdlabel"])
         return output
-        
+
+    def remove_all_cut(self): 
+        """remove all the cut"""
+
+        for name in self.cuts_parameter:
+            targettype = type(self[name])
+            if targettype == bool:
+                self[name] = False
+            elif 'min' in name:
+                self[name] = 0
+            elif 'max' in name:
+                self[name] = -1
+            elif 'eta' in name:
+                self[name] = -1
+            else:
+                self[name] = 0       
 
 class RunCardLO(RunCard):
     """an object to handle in a nice way the run_card infomration"""
@@ -1798,22 +1814,6 @@
                 self['drjl'] = 0
                 self['sys_alpsfact'] = "0.5 1 2"
                 
-
-    def remove_all_cut(self): 
-        """remove all the cut"""
-
-        for name in self.cuts_parameter:
-            targettype = type(self[name])
-            if targettype == bool:
-                self[name] = False
-            elif 'min' in name:
-                self[name] = 0
-            elif 'max' in name:
-                self[name] = -1
-            elif 'eta' in name:
-                self[name] = -1
-            else:
-                self[name] = 0
             
     def write(self, output_file, template=None, python_template=False):
         """Write the run_card in output_file according to template 
@@ -1984,7 +1984,7 @@
         """Rules
           e+ e- beam -> lpp:0 ebeam:500  
           p p beam -> set maxjetflavor automatically
-         """
+        """
 
         # check for beam_id
         beam_id = set()
@@ -2004,6 +2004,10 @@
         else:
             self['lpp1'] = 0
             self['lpp2'] = 0  
+            
+        if proc_characteristic['ninitial'] == 1:
+            #remove all cut
+            self.remove_all_cut()
         
 class MadLoopParam(ConfigFile):
     """ a class for storing/dealing with the file MadLoopParam.dat

=== 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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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-20 10:56:39 +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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/long_ML_SMQCD_default/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/dux_mumvmxg/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       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	2015-08-14 03:43:29 +0000
+++ tests/input_files/IOTestsComparison/short_ML_SMQCD_optimized/gg_ttx/%..%..%Source%MODEL%model_functions.inc	2015-10-20 10:56:39 +0000
@@ -9,6 +9,8 @@
       DOUBLE COMPLEX REGLOGM
       DOUBLE COMPLEX RECMS
       DOUBLE COMPLEX ARG
+
+
       COMPLEX*32 MP_COND
       COMPLEX*32 MP_CONDIF
       COMPLEX*32 MP_REGLOG
@@ -16,3 +18,5 @@
       COMPLEX*32 MP_REGLOGM
       COMPLEX*32 MP_RECMS
       COMPLEX*32 MP_ARG
+
+

=== 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-20 10:56:39 +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):