Merge lp:~maddevelopers/mg5amcnlo/2.2.2 into lp:mg5amcnlo/lts

[Olivier Mattelaer]

I think that this is time for this version to be released since it contains a couple of important fixes already. I'm also investigating a problem at LO which might have huge impact on the efficiency of the code. This is not yet pushed in this version but I suggest that you take a look already.

If nobody disagree, I plan to release this one on Friday.

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/2.2.2]

276. By Rikkert Frederix on 2014-10-30

Added the 'strip' command to the makefiles to reduce the sizes of the
executables.

277. By Rikkert Frederix on 2014-10-30

updated the UpdateNotes

278. By Rikkert Frederix on 2014-10-30

another update for the UpdateNotes

[Rikkert Frederix (frederix)]

Please, everybody who made changes in the branch, put a comment in the UpdateNotes.txt

Cheers,
Rik

review: Needs Fixing

[lp:~maddevelopers/mg5amcnlo/2.2.2]

279. By marco zaro on 2014-10-30

added possibility to install lhapdf sets in LHAPDF_DATA_PATH if it is set as an environmental variable.

[Paolo Torrielli (paolo-torrielli)]

Hi Olivier,
tomorrow morning I’ll merge the 2.2.2_py82 branch,
so please wait for this merging before releasing.
Cheers.
Paolo

On 29 Oct 2014, at 16:02, Olivier Mattelaer <email address hidden> wrote:

> Olivier Mattelaer has proposed merging lp:~maddevelopers/mg5amcnlo/2.2.2 into lp:mg5amcnlo.
>
> Requested reviews:
> MadTeam (madteam)
>
> For more details, see:
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
>
> I think that this is time for this version to be released since it contains a couple of important fixes already. I'm also investigating a problem at LO which might have huge impact on the efficiency of the code. This is not yet pushed in this version but I suggest that you take a look already.
>
> If nobody disagree, I plan to release this one on Friday.
>
> Cheers,
>
> Olivier
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/mg5amcnlo/2.2.2.
> === modified file 'MadSpin/decay.py'
> --- MadSpin/decay.py 2014-09-17 09:53:23 +0000
> +++ MadSpin/decay.py 2014-10-29 15:01:21 +0000
> @@ -1185,7 +1185,7 @@
> decaying_base = [id for id in final if abs(id) in self.decay_ids]
> decaying_base.sort()
> topo['base_order'] = (initial , final)
> - topo['matrix_element'] = me
> +# topo['matrix_element'] = me
> tags = []
> topo['tag2order'] = {}
>
> @@ -2662,6 +2662,8 @@
> decay_text = []
> for decays in self.mscmd.list_branches.values():
> for decay in decays:
> + if '=' not in decay:
> + decay += ' QCD=99'
> if ',' in decay:
> decay_text.append('(%s)' % decay)
> else:
> @@ -2791,11 +2793,11 @@
> 'processes': matrix_element.get('processes'),
> 'tag': me.shell_string(pdg_order=True)}
> #
> - if __debug__:
> - #check that all decay matrix element correspond to a decay only
> - for prod in self.all_ME.values():
> - for decay in prod['matrix_element']['base_amplitude']['process']['decay_chains']:
> - assert decay.shell_string() in self.all_decay
> +# if __debug__:
> +# #check that all decay matrix element correspond to a decay only
> +# for prod in self.all_ME.values():
> +# for decay in prod['matrix_element']['base_amplitude']['process']['decay_chains']:
> +# assert decay.shell_string() in self.all_decay
>
>
> def get_branching_ratio(self):
>
> === modified file 'Template/LO/Source/PDF/pdf.f'
> --- Template/LO/Source/PDF/pdf.f 2014-05-16 08:15:16 +0000
> +++ Template/LO/Source/PDF/pdf.f 2014-10-29 15:01:21 +0000
> @@ -45,6 +45,7 @@
> do i=-5,5
> fx(i)=nnfx(i)/x
> enddo
> + fx(7)=nnfx(7)/x
> elseif ((pdlabel(1:3) .eq. 'mrs')
> . .or. (pdlabel(2:4) .eq. 'mrs')) then
>
>
> === modified file 'Template/LO/Source/rw_events.f'
> --- Template/LO/S...

[lp:~maddevelopers/mg5amcnlo/2.2.2]

280. By marco zaro on 2014-10-31

the default accuracy for check_poles is taken from FKSparams.dat
Fix in ajob relative to HelConfigs.dat (see Question #256309)

[marco zaro (marco-zaro)]

Hi,
btw there are some unit tests failing:
======================================================================
ERROR: test_write_lh_order_EW (tests.unit_tests.iolibs.test_export_fks_EW.IOExportFKSEWTest)
tests the correct writing of the B-LH order file
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/tests/unit_tests/iolibs/test_export_fks_EW.py", line 145, in test_write_lh_order_EW
    self.myfks_me)
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/madgraph/iolibs/export_fks.py", line 1474, in write_lh_order
    replace_dict['corr'] = ' '.join(matrix_elements[0].get('processes')[0].\
TypeError: 'FKSHelasProcess' object does not support indexing

======================================================================
FAIL: test_export_matrix_element_v4_madevent_group (tests.unit_tests.iolibs.test_export_v4.IOExportV4Test)
Test the result of exporting a subprocess group matrix element
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/tests/unit_tests/iolibs/test_export_v4.py", line 2134, in test_export_matrix_element_v4_madevent_group
    self.assertFileContains('test', goal_matrix1)
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/tests/unit_tests/iolibs/test_file_writers.py", line 52, in assertFileContains
    self.assertEqual(a,b)
AssertionError: ' JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))' != ' JAMP2(I)=JAMP2(I)+DABS(DBLE(JAMP(I,M)*DCONJG(JAMP(I,N))))'

======================================================================
FAIL: test_export_matrix_element_v4_madevent_nogroup (tests.unit_tests.iolibs.test_export_v4.IOExportV4Test)
Test the result of exporting a subprocess group matrix element
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/tests/unit_tests/iolibs/test_export_v4.py", line 1139, in test_export_matrix_element_v4_madevent_nogroup
    self.assertFileContains('test', goal_matrix1)
  File "/Users/marcozaro/Physics/MadGraph/2.2.2/tests/unit_tests/iolibs/test_file_writers.py", line 52, in assertFileContains
    self.assertEqual(a,b)
AssertionError: ' JAMP2(I)=JAMP2(I)+JAMP(I,M)*DCONJG(JAMP(I,N))' != ' JAMP2(I)=JAMP2(I)+DABS(DBLE(JAMP(I,M)*DCONJG(JAMP(I,N))))'

----------------------------------------------------------------------

Marco

review: Needs Fixing

[lp:~maddevelopers/mg5amcnlo/2.2.2]

281. By Paolo Torrielli on 2014-10-31

merged with the py82 branch

282. By marco zaro on 2014-10-31

merged with 2.2.2_read_conf_props

[Rikkert Frederix (frederix)]

I've got the same 3 unit tests failing as Marco. I guess the DABS() problem is not serious, but the first one related to the order file might be. Although also that one is only used when not linking with MadLoop.

Valentin, do you have an idea what the problem is?

Rikkert

review: Approve

[Rikkert Frederix (frederix)]

review: Needs Fixing

[Valentin Hirschi (valentin-hirschi)]

So, these three changes are related to two small modifications I brought
directly to 2.2.2 a week ago or so.

The new 'DABS(DBLE(' is correct as we recently realized with Olivier, that
the former expression was not correct as it could yields negative JAMP**2
in presence of interferences.

I think the first crash is due to the '[0]' I appended to
'join(matrix_elements'
because it was necessary when writing the order file for the GoSam
interface. However, apparently, in the context of this test 'IOExportFKSEWTest'
the type of the object 'matrix_elements' passed is not the same as in an
actual run. So I think the test should be updated so as to pass a list of
matrix_elements instead.
Marco, do you confirm that this would be ok?

Cheers,

On Fri, Oct 31, 2014 at 6:20 AM, Rikkert Frederix <email address hidden>
wrote:

> Review: Needs Fixing
>
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
> Your team MadDevelopers is subscribed to branch
> lp:~maddevelopers/mg5amcnlo/2.2.2.
>

--
Valentin

[lp:~maddevelopers/mg5amcnlo/2.2.2]

283. By Olivier Mattelaer on 2014-11-01

1. saver use of gzip module (not used for file > 4Gb)
2. add security in MadWidth for scale close to Lambda_QCD
3. improve the security associate to the pickle of the model
4. improve cluster support for lsf/nlo run (thanks Josh)

284. By marco zaro on 2014-11-03

fixed failing test test_write_lh_order_EW

[marco zaro (marco-zaro)]

ciao,
i have fixed the
test_weite_lh_order_ew passing the correct ‘matrix_elements’, and pushed to 2.2.2
Cheers,

Marco
On 31 Oct 2014, at 19:45, Valentin Hirschi <email address hidden> wrote:

> So, these three changes are related to two small modifications I brought
> directly to 2.2.2 a week ago or so.
>
> The new 'DABS(DBLE(' is correct as we recently realized with Olivier, that
> the former expression was not correct as it could yields negative JAMP**2
> in presence of interferences.
>
> I think the first crash is due to the '[0]' I appended to
> 'join(matrix_elements'
> because it was necessary when writing the order file for the GoSam
> interface. However, apparently, in the context of this test 'IOExportFKSEWTest'
> the type of the object 'matrix_elements' passed is not the same as in an
> actual run. So I think the test should be updated so as to pass a list of
> matrix_elements instead.
> Marco, do you confirm that this would be ok?
>
> Cheers,
>
> On Fri, Oct 31, 2014 at 6:20 AM, Rikkert Frederix <email address hidden>
> wrote:
>
>> Review: Needs Fixing
>>
>>
>> --
>> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
>> Your team MadDevelopers is subscribed to branch
>> lp:~maddevelopers/mg5amcnlo/2.2.2.
>>
>
>
>
> --
> Valentin
>
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
> You are reviewing the proposed merge of lp:~maddevelopers/mg5amcnlo/2.2.2 into lp:mg5amcnlo.

[lp:~maddevelopers/mg5amcnlo/2.2.2]

285. By Marco Zaro <email address hidden> on 2014-11-03

fix in check_poles.f (case error) and FKSParamsReader (problem with declaration of filename variable)

286. By Marco Zaro <email address hidden> on 2014-11-03

added required output for fNLO runs

287. By Olivier Mattelaer on 2014-11-03

Fix a small crash in MadWidth

288. By Olivier Mattelaer on 2014-11-03

fix the remaining unittest

289. By marco zaro on 2014-11-04

randinit is also symlinked to the P* dirs

290. By Olivier Mattelaer on 2014-11-04

safer version of checking if MG5 is on a readonly disk

[Olivier Mattelaer (olivier-mattelaer)]

Hi Valentin,

What is the status with the IOTest which were crashing?
Please do it asap.
I found a real deep bug, (I'm evaluating how deep it is with Tim and Fabio) and need to have this version ready as soon as we have finish our investigation).

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/2.2.2]

291. By Olivier Mattelaer on 2014-11-05

- correct the bug of a a > mu+ mu-
- bypass the call to do_plots for pythia8/herwig++ shower
- fix a problem with the set command (in AskforEdit)for HC model
- fix a problem with the name of the particle for the compute_width command

292. By Valentin Hirschi on 2014-11-06

1. Updated IOTests

[Valentin Hirschi (valentin-hirschi)]

Hi Olivier,

IOTest pass on my computer, so what do you mean by fixing it? Make them
pass as well on your computer despite the fact that you don't have PJFry or
Golem? If so, I think it can be done for the next release no? Also I have
to update the parsing for the date and version appearing in the MadLoop
banner.

Cheers

On Wed, Nov 5, 2014 at 6:00 AM, Olivier Mattelaer <
<email address hidden>> wrote:

> Hi Valentin,
>
> What is the status with the IOTest which were crashing?
> Please do it asap.
> I found a real deep bug, (I'm evaluating how deep it is with Tim and
> Fabio) and need to have this version ready as soon as we have finish our
> investigation).
>
> Cheers,
>
> Olivier
>
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.2.2/+merge/240012
> Your team MadTeam is subscribed to branch lp:mg5amcnlo.
>

--
Valentin

[Rikkert Frederix (frederix)]

I think we can make this version public.

Rikkert

review: Approve

=== modified file 'MadSpin/decay.py'
--- MadSpin/decay.py	2014-09-17 09:53:23 +0000
+++ MadSpin/decay.py	2014-11-06 06:25:51 +0000
@@ -1185,7 +1185,7 @@
             decaying_base = [id for id in final if abs(id) in self.decay_ids]
             decaying_base.sort()
             topo['base_order'] = (initial , final)
-            topo['matrix_element'] = me
+#            topo['matrix_element'] = me
             tags = []
             topo['tag2order'] = {}
              
@@ -2662,6 +2662,8 @@
         decay_text = []
         for decays in self.mscmd.list_branches.values():
             for decay in  decays:
+                if '=' not in decay:
+                    decay += ' QCD=99'
                 if ',' in decay:
                     decay_text.append('(%s)' % decay)
                 else:
@@ -2791,11 +2793,11 @@
                                          'processes': matrix_element.get('processes'),
                                          'tag': me.shell_string(pdg_order=True)}
         #
-        if __debug__:
-            #check that all decay matrix element correspond to a decay only
-            for prod in self.all_ME.values():
-                for decay in prod['matrix_element']['base_amplitude']['process']['decay_chains']:
-                    assert decay.shell_string() in self.all_decay
+#        if __debug__:
+#            #check that all decay matrix element correspond to a decay only
+#            for prod in self.all_ME.values():
+#                for decay in prod['matrix_element']['base_amplitude']['process']['decay_chains']:
+#                    assert decay.shell_string() in self.all_decay
             
         
     def get_branching_ratio(self):

=== modified file 'Template/LO/Source/PDF/pdf.f'
--- Template/LO/Source/PDF/pdf.f	2014-05-16 08:15:16 +0000
+++ Template/LO/Source/PDF/pdf.f	2014-11-06 06:25:51 +0000
@@ -45,6 +45,7 @@
          do i=-5,5
             fx(i)=nnfx(i)/x
          enddo
+         fx(7)=nnfx(7)/x
       elseif     ((pdlabel(1:3) .eq. 'mrs')
      .   .or. (pdlabel(2:4) .eq. 'mrs')) then
 

=== modified file 'Template/LO/Source/dsample.f'
--- Template/LO/Source/dsample.f	2014-06-23 14:11:08 +0000
+++ Template/LO/Source/dsample.f	2014-11-06 06:25:51 +0000
@@ -974,7 +974,7 @@
          im = ng -1
          ddum(j) = ng
       endif
-      if (im.le.0) im = 1
+      if (im.lt.0) im = 0
       ip = im + 1
       ij = Minvar(j,ipole)
 c------

=== modified file 'Template/LO/Source/rw_events.f'
--- Template/LO/Source/rw_events.f	2014-06-27 13:31:53 +0000
+++ Template/LO/Source/rw_events.f	2014-11-06 06:25:51 +0000
@@ -142,7 +142,7 @@
       integer nexternal, ic(7,*)
       double precision P(0:4,*),wgt
       double precision aqcd, aqed, scale
-      character*300 buff
+      character*1000 buff
       logical u_syst
       character*(s_bufflen) s_buff(*)
       integer nclus

=== modified file 'Template/LO/SubProcesses/setcuts.f'
--- Template/LO/SubProcesses/setcuts.f	2014-05-05 15:00:26 +0000
+++ Template/LO/SubProcesses/setcuts.f	2014-11-06 06:25:51 +0000
@@ -95,13 +95,17 @@
      $     ,IS_A_ONIUM_SAVE(NEXTERNAL),IS_A_NU_SAVE(NEXTERNAL)
      $     ,IS_HEAVY_SAVE(NEXTERNAL),DO_CUTS_SAVE(NEXTERNAL)
       double precision r2min_save(nexternal,nexternal)
-     $     ,s_min_save(nexternal,nexternal),s_max_save(nexternal
-     $     ,nexternal),ptll_min_save(nexternal,nexternal)
-     $     ,ptll_max_save(nexternal,nexternal)
+     $     ,r2max_save(nexternal,nexternal),s_min_save(nexternal
+     $     ,nexternal),s_max_save(nexternal,nexternal)
+     $     ,ptll_min_save(nexternal,nexternal),ptll_max_save(nexternal
+     $     ,nexternal),etmin_save(nexternal),etmax_save(nexternal)
+     $     ,emin_save(nexternal) ,emax_save(nexternal)
+     $     ,etamin_save(nexternal),etamax_save(nexternal)
       save  IS_A_J_SAVE,IS_A_L_SAVE,IS_A_B_SAVE,IS_A_A_SAVE
      $     ,IS_A_ONIUM_SAVE,IS_A_NU_SAVE,IS_HEAVY_SAVE,DO_CUTS
-     $     ,r2min_save,s_min_save,s_max_save,ptll_min_save
-     $     ,ptll_max_save
+     $     ,r2min_save,r2max_save,s_min_save,s_max_save,ptll_min_save
+     $     ,ptll_max_save,etmin_save,etmax_save,emin_save,emax_save
+     $     ,etamin_save,etamax_save
 c
 c     setup masses for the final-state particles
 c
@@ -419,14 +423,14 @@
             ptll_min(j,i)=0.0d0**2
             ptll_max(j,i)=-1
             if(((is_a_l(i).and.is_a_l(j)).and.
-     &      (abs(idup(i,1,iproc)).eq.abs(idup(j,1,iproc))).and.
-     &      (idup(i,1,iproc)*idup(j,1,iproc).lt.0))! Leptons from same flavor but different charge
-     &       .or.(is_a_nu(i).and.is_a_l(j))
-     &       .or.(is_a_l(i).and.is_a_nu(j))   !a lepton and a neutrino
-     &       .or.(is_a_nu(i).and.is_a_nu(j))) then ! two neutrinos 
+     &           (abs(idup(i,1,iproc)).eq.abs(idup(j,1,iproc))).and.
+     &           (idup(i,1,iproc)*idup(j,1,iproc).lt.0)) ! Leptons from same flavor but different charge
+     &           .or.(is_a_nu(i).and.is_a_l(j))
+     &           .or.(is_a_l(i).and.is_a_nu(j)) !a lepton and a neutrino
+     &           .or.(is_a_nu(i).and.is_a_nu(j))) then ! two neutrinos 
                ptll_min(j,i)=ptllmin*dabs(ptllmin)
                ptll_max(j,i)=ptllmax*dabs(ptllmax)
-           endif
+            endif
          enddo
       enddo
 
@@ -506,16 +510,23 @@
       if (iproc.eq.1) then
          do i=nincoming+1,nexternal
             do_cuts_save(i)=do_cuts(i)
-            is_a_j_save(i) = is_a_j(i)
-            is_a_l_save(i) = is_a_l(i)
-            is_a_b_save(i) = is_a_b(i)
-            is_a_a_save(i) = is_a_a(i)
-            is_a_nu_save(i) = is_a_nu(i)
-            is_heavy_save(i) = is_heavy(i)
+            is_a_j_save(i)=is_a_j(i)
+            is_a_b_save(i)=is_a_b(i)
+            is_a_l_save(i)=is_a_l(i)
+            is_a_a_save(i)=is_a_a(i)
+            is_a_nu_save(i)=is_a_nu(i)
             is_a_onium_save(i)=is_a_onium(i)
+            is_heavy_save(i)=is_heavy(i) 
+            etmin_save(i)=etmin(i)
+            etmax_save(i)=etmax(i)
+            emin_save(i)=emin(i)
+            emax_save(i)=emax(i)
+            etamin_save(i)=etamin(i)
+            etamax_save(i)=etamax(i)
             if (i.eq.nexternal) cycle
             do j=i+1,nexternal
                r2min_save(j,i) = r2min(j,i)
+               r2max_save(j,i) = r2max(j,i)
                s_min_save(j,i) = s_min(j,i)
                s_max_save(j,i) = s_max(j,i)
                ptll_min_save(j,i) = ptll_min(j,i)
@@ -526,16 +537,55 @@
          equal=.true.
          do i=nincoming+1,nexternal
             if (do_cuts_save(i).neqv.do_cuts(i)) equal=.false.
-            if (is_a_j_save(i).neqv.is_a_j(i)) equal=.false.
-            if (is_a_l_save(i).neqv.is_a_l(i)) equal=.false.
-            if (is_a_b_save(i).neqv.is_a_b(i)) equal=.false.
-            if (is_a_a_save(i).neqv.is_a_a(i)) equal=.false.
-            if (is_a_nu_save(i).neqv.is_a_nu(i)) equal=.false.
-            if (is_heavy_save(i).neqv.is_heavy(i)) equal=.false.
-            if (is_a_onium_save(i).neqv.is_a_onium(i)) equal=.false.
+            if (is_a_j_save(i).neqv.is_a_j(i)) then
+               if (ptjmin4(1).gt.0d0 .or. ptjmax4(1).ge.0d0) equal=.false.
+               if (ptjmin4(2).gt.0d0 .or. ptjmax4(2).ge.0d0) equal=.false.
+               if (ptjmin4(3).gt.0d0 .or. ptjmax4(3).ge.0d0) equal=.false.
+               if (ptjmin4(4).gt.0d0 .or. ptjmax4(4).ge.0d0) equal=.false.
+               if (Htjmin4(2).gt.0d0 .or. Htjmax4(2).ge.0d0) equal=.false.
+               if (Htjmin4(3).gt.0d0 .or. Htjmax4(3).ge.0d0) equal=.false.
+               if (Htjmin4(4).gt.0d0 .or. Htjmax4(4).ge.0d0) equal=.false.
+               if (inclHtmin.gt.0d0 .or. inclHtmax.ge.0d0) equal=.false.
+               if (htjmin.gt.0d0 .or. htjmax.ge.0d0) equal=.false.
+               if (xptj.gt.0d0) equal=.false.
+               if (xetamin.gt.0d0 .or. deltaeta.gt.0d0) equal=.false.
+               if (ptgmin.ne.0d0) equal=.false.
+               if (kt_durham.gt.0d0) equal=.false.
+            endif
+            if (is_a_b_save(i).neqv.is_a_b(i)) then
+               if (inclHtmin.gt.0d0 .or. inclHtmax.ge.0d0) equal=.false.
+               if (xptb.gt.0d0) equal=.false.
+            endif
+            if (is_a_a_save(i).neqv.is_a_a(i)) then
+               if (xpta.gt.0d0) equal=.false.
+               if (ptgmin.ne.0d0) equal=.false.
+            endif
+            if (is_a_l_save(i).neqv.is_a_l(i)) then
+               if (ptlmin4(1).gt.0d0 .or. ptlmax4(1).ge.0d0) equal=.false.
+               if (ptlmin4(2).gt.0d0 .or. ptlmax4(2).ge.0d0) equal=.false.
+               if (ptlmin4(3).gt.0d0 .or. ptlmax4(3).ge.0d0) equal=.false.
+               if (ptlmin4(4).gt.0d0 .or. ptlmax4(4).ge.0d0) equal=.false.
+               if (mmnl.gt.0d0 .or. mmnlmax.ge.0d0) equal=.false.
+               if (xptl.gt.0d0) equal=.false.
+               if (ptgmin.ne.0d0 .and. isoEM) equal=.false.
+            endif
+            if (is_a_nu_save(i).neqv.is_a_nu(i)) then
+               if (misset.gt.0d0 .or. missetmax.ge.0d0) equal=.false.
+               if (mmnl.gt.0d0 .or. mmnlmax.ge.0d0) equal=.false.
+            endif
+            if (is_heavy_save(i).neqv.is_heavy(i)) then
+               if (ptheavy.gt.0d0) equal=.false.
+            endif
+            if (etmin_save(i).ne.etmin(i)) equal=.false.
+            if (etmax_save(i).ne.etmax(i)) equal=.false.
+            if (emin_save(i).ne.emin(i)) equal=.false.
+            if (emax_save(i).ne.emax(i)) equal=.false.
+            if (etamin_save(i).ne.etamin(i)) equal=.false.
+            if (etamax_save(i).ne.etamax(i)) equal=.false.
             if (i.eq.nexternal) cycle
             do j=i+1,nexternal
                if (r2min_save(j,i).ne.r2min(j,i)) equal=.false.
+               if (r2max_save(j,i).ne.r2max(j,i)) equal=.false.
                if (s_min_save(j,i).ne.s_min(j,i)) equal=.false.
                if (s_max_save(j,i).ne.s_max(j,i)) equal=.false.
                if (ptll_min_save(j,i).ne.ptll_min(j,i)) equal=.false.
@@ -626,7 +676,7 @@
 
 c     set possible xqcut combinations (for better grid preparation)
         if(xqcut.gt.0)
-     $       call setxqcuts()
+     $       call setxqcuts(do_cuts)
 
 c      call write_cuts()
       RETURN
@@ -634,7 +684,7 @@
       END
 
 
-      subroutine setxqcuts()
+      subroutine setxqcuts(do_cuts)
 c**************************************************
 c     Set xqcuti and xqcutij between all particles
 c     to allow for grid preparation based on xqcut
@@ -645,6 +695,8 @@
       include 'nexternal.inc'
       include 'cuts.inc'
 
+      logical  do_cuts(nexternal)
+
       double precision pmass(nexternal)
       common/to_mass/  pmass
       integer iforest(2,-max_branch:-1,lmaxconfigs)
@@ -682,13 +734,15 @@
      $                 idup(iforest(2,k,j),1,1).eq.21.or.
      $                 iabs(idup(iforest(2,k,j),1,1)).le.maxjetflavor.or.
      $                 idup(iforest(2,k,j),1,1).eq.21)then
-                     xqcutij(iforest(2,k,j),iforest(1,k,j))=xqcut
-                     xqcutij(iforest(1,k,j),iforest(2,k,j))=xqcut
+                     if (do_cuts(iforest(2,k,j)).and. do_cuts(iforest(1,k,j)))then
+                        xqcutij(iforest(2,k,j),iforest(1,k,j))=xqcut
+                        xqcutij(iforest(1,k,j),iforest(2,k,j))=xqcut
+                     endif
                   endif
                endif
             enddo
             if(.not.foundpartner.and.(iabs(idup(i,1,1)).le.maxjetflavor.or.
-     $           idup(i,1,1).eq.21)) xqcuti(i)=max(0d0,sqrt(xqcut**2-pmass(i)**2))
+     $           idup(i,1,1).eq.21).and.do_cuts(i)) xqcuti(i)=max(0d0,sqrt(xqcut**2-pmass(i)**2))
          enddo
       enddo
 

=== modified file 'Template/LO/SubProcesses/unwgt.f'
--- Template/LO/SubProcesses/unwgt.f	2014-02-19 01:00:53 +0000
+++ Template/LO/SubProcesses/unwgt.f	2014-11-06 06:25:51 +0000
@@ -241,7 +241,7 @@
       integer iseed, nover, nstore
       double precision scale,aqcd,aqed
       integer ievent
-      character*300 buff
+      character*1000 buff
       logical u_syst
       character*(s_bufflen) s_buff(7)
       integer nclus
@@ -408,7 +408,7 @@
       real ran1
       external ran1
 
-      character*300 buff
+      character*1000 buff
       character*(s_bufflen) s_buff(7)
       integer nclus
       character*(clus_bufflen) buffclus(nexternal)

=== modified file 'Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIG6.Script'
--- Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIG6.Script	2014-07-10 16:03:27 +0000
+++ Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIG6.Script	2014-11-06 06:25:51 +0000
@@ -92,7 +92,6 @@
 EXTRAPOLATE|extrapolate) PDFGROUP=LHAEXT ;;
 *) echo "no such option; failure in linklhapdf" ; exit 1 ;;
 esac
-source ../Source/fj_lhapdf_opts
 }    
 
 # creates the executable for NLO and MC; the entry of the function
@@ -295,7 +294,7 @@
   fi
 fi
 echo " "
-echo "*****   Now ready for showering" $NEVENTS "events with $MC_TYPE "
+echo "*****   Now ready for showering" $NEVENTS "events with $MCMODE "
 echo " "
 runMCMadFKS
 teststatus runMC

=== modified file 'Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIGPP.Script'
--- Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIGPP.Script	2014-05-02 11:09:41 +0000
+++ Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIGPP.Script	2014-11-06 06:25:51 +0000
@@ -92,7 +92,6 @@
 EXTRAPOLATE|extrapolate) PDFGROUP=LHAEXT ;;
 *) echo "no such option; failure in linklhapdf" ; exit 1 ;;
 esac
-source ../Source/fj_lhapdf_opts
 }    
 
 # converts pdg numbers to names, relevant for decays
@@ -337,7 +336,7 @@
 let "ERR_NUM_HWPP += 1"
 
 echo " "
-echo "*****   Now ready for showering" $NEVENTS "events with $MC_TYPE "
+echo "*****   Now ready for showering" $NEVENTS "events with $MCMODE "
 echo " "
 runMCMadFKS
 teststatus runMC

=== modified file 'Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6PT.Script'
--- Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6PT.Script	2014-05-02 11:09:41 +0000
+++ Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6PT.Script	2014-11-06 06:25:51 +0000
@@ -92,7 +92,6 @@
 EXTRAPOLATE|extrapolate) PDFGROUP=LHAEXT ;;
 *) echo "no such option; failure in linklhapdf" ; exit 1 ;;
 esac
-source ../Source/fj_lhapdf_opts
 }    
 
 # creates the executable for NLO and MC; the entry of the function
@@ -296,7 +295,7 @@
   fi
 fi
 echo " "
-echo "*****   Now ready for showering" $NEVENTS "events with $MC_TYPE "
+echo "*****   Now ready for showering" $NEVENTS "events with $MCMODE "
 echo " "
 runMCMadFKS
 teststatus runMC

=== modified file 'Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6Q.Script'
--- Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6Q.Script	2014-05-02 11:09:41 +0000
+++ Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA6Q.Script	2014-11-06 06:25:51 +0000
@@ -92,7 +92,6 @@
 EXTRAPOLATE|extrapolate) PDFGROUP=LHAEXT ;;
 *) echo "no such option; failure in linklhapdf" ; exit 1 ;;
 esac
-source ../Source/fj_lhapdf_opts
 }    
 
 # creates the executable for NLO and MC; the entry of the function
@@ -296,7 +295,7 @@
   fi
 fi
 echo " "
-echo "*****   Now ready for showering" $NEVENTS "events with $MC_TYPE "
+echo "*****   Now ready for showering" $NEVENTS "events with $MCMODE "
 echo " "
 runMCMadFKS
 teststatus runMC

=== modified file 'Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA8.Script'
--- Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA8.Script	2014-07-16 16:23:08 +0000
+++ Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_PYTHIA8.Script	2014-11-06 06:25:51 +0000
@@ -92,7 +92,6 @@
 EXTRAPOLATE|extrapolate) PDFGROUP=LHAEXT ;;
 *) echo "no such option; failure in linklhapdf" ; exit 1 ;;
 esac
-source ../Source/fj_lhapdf_opts
 }    
 
 # creates the executable for NLO and MC; the entry of the function
@@ -211,23 +210,43 @@
 # new compilation step for PY8 mode
 export PYTHIA8LOCATION=$PY8PATH
 
+# check the PYTHIA8 version by the position of certain files
+
+if [ ! -d $PYTHIA8LOCATION/xmldoc ]
+then
+ PY8VER="8.2"
+ echo ' USING PYTHIA 8, version >= 8.200'
+ chmod +x $PYTHIA8LOCATION/bin/pythia8-config
+ read HEPMCINCLIB <<< $($PYTHIA8LOCATION/bin/pythia8-config --hepmc2)
+else
+ PY8VER="8.1"
+ echo ' USING PYTHIA 8, version < 8.200'
+fi
+
 if [ "$PY8UTI" != "" ]
 then
- cp -rf $thisdir/srcPythia8/Pythia8.cc $thisdir/Pythia8.cc
+ if [ $PY8VER == "8.1" ]; then cp -rf $thisdir/srcPythia8/Pythia8.cc $thisdir/Pythia8.cc; fi
+ if [ $PY8VER == "8.2" ]; then cp -rf $thisdir/srcPythia8/Pythia82.cc $thisdir/Pythia82.cc; fi
  cp -rf $thisdir/srcPythia8/Makefile $thisdir/Makefile
 else
- cp -rf $thisdir/srcPythia8/Pythia8_hep.cc $thisdir/Pythia8.cc
+ if [ $PY8VER == "8.1" ]; then cp -rf $thisdir/srcPythia8/Pythia8_hep.cc $thisdir/Pythia8.cc; fi
+ if [ $PY8VER == "8.2" ]; then cp -rf $thisdir/srcPythia8/Pythia82_hep.cc $thisdir/Pythia82.cc; fi
  cp -rf $thisdir/srcPythia8/Makefile_hep $thisdir/Makefile
 fi
 if [ -f $thisdir/Pythia8.exe ]
 then
   rm -rf $thisdir/Pythia8.exe
 fi
-cp -rf $PYTHIA8LOCATION/xmldoc $thisdir
+
+if [ $PY8VER == "8.1" ]; then cp -rf $PYTHIA8LOCATION/xmldoc $thisdir; fi
+if [ $PY8VER == "8.2" ]; then cp -rf $PYTHIA8LOCATION/share/Pythia8/xmldoc $thisdir; fi
 PY8UTIMAKE="PY8UTI="$PY8UTI
 FORFILEMAKE="FORFILE="$PY8UTI
+HEPMCMAKE="HEPMCINCLIB="$HEPMCINCLIB
 LGFORTRANMAKE="LIBGFORTRAN=-lgfortran"
-make Pythia8 "$PY8UTIMAKE" "$FORFILEMAKE" "$LGFORTRANMAKE" "$LIBSMAKE" "$LIBSMAKEP" "$INCLMAKE"
+
+if [ $PY8VER == "8.1" ]; then make Pythia8 "$PY8UTIMAKE" "$FORFILEMAKE" "$LGFORTRANMAKE" "$LIBSMAKE" "$LIBSMAKEP" "$INCLMAKE"; fi
+if [ $PY8VER == "8.2" ]; then make Pythia82 "$PY8UTIMAKE" "$FORFILEMAKE" "$LGFORTRANMAKE" "$LIBSMAKE" "$LIBSMAKEP" "$INCLMAKE" "$HEPMCMAKE"; fi
 
 cd $thisdir
 if [ ! -f $thisdir/Pythia8.exe ]
@@ -296,7 +315,7 @@
 let "ERR_NUM_PY8 += 1"
 
 echo " "
-echo "*****   Now ready for showering" $NEVENTS "events with $MC_TYPE "
+echo "*****   Now ready for showering" $NEVENTS "events with $MCMODE "
 echo " "
 runMCMadFKS
 teststatus runMC
@@ -320,8 +339,6 @@
 Main:spareWord1 = $EVENT_NORM         ! Event weights are normalized to sum 
                                       ! or average to the cross section
 Main:timesAllowErrors = $ERR_NUM_PY8  ! Number of allowed errors
-Main:showChangedSettings = on         ! Shows all non-default settings
-Main:showChangedParticleData = off    ! Shows all non-default particle settings
 
 ! 2) Settings related to output in init(), next(), amd stat().
 Init:showChangedSettings = on         ! Shows all non-default settings
@@ -344,7 +361,9 @@
 EOF
 if [ $UsedPdfLib = "LHAPDF" ]
 then
-cat <<EOF >> $ifile
+ if [ $PY8VER == "8.1" ]
+ then
+  cat <<EOF >> $ifile
 PDF:useLHAPDF = on                 ! Use of LHAPDF
 PDF:LHAPDFSET = $PDFSETN           ! PDF set
 EOF
@@ -362,9 +381,41 @@
   echo "Unknown PDFGROUP " $PDFGROUP
   exit 1
   fi
+ else
+  if [ -f $PYTHIA8LOCATION/lib/libpythia8lhapdf5.so ]
+  then
+  cat <<EOF >> $ifile
+PDF:pSet = LHAPDF5:$PDFSETN
+EOF
+   if [ $PDFGROUP = "LHAPDF" ]
+   then
+   cat <<EOF >> $ifile
+PDF:extrapolate = off        ! extrapolate PDF set outside the boundaries
+EOF
+   elif [ $PDFGROUP = "LHAEXT" ]
+   then
+   cat <<EOF >> $ifile
+PDF:extrapolate = on         ! extrapolate PDF set outside the boundaries
+EOF
+   else
+   echo "Unknown PDFGROUP " $PDFGROUP
+   exit 1
+   fi
+  elif [ -f $PYTHIA8LOCATION/lib/libpythia8lhapdf6.so ]
+  then
+  cat <<EOF >> $ifile
+PDF:pSet = LHAPDF6:$PDFSETN
+EOF
+  fi
+ fi
 else
-cat <<EOF >> $ifile
+ if [ $PY8VER == "8.1" ]
+ then
+ cat <<EOF >> $ifile
 PDF:useLHAPDF = off                ! Use of LHAPDF
+EOF
+ fi
+ cat <<EOF >> $ifile
 PDF:pSet = 7                       ! CTEQ6L
 EOF
 fi

=== modified file 'Template/NLO/MCatNLO/include/LHEF.h'
--- Template/NLO/MCatNLO/include/LHEF.h	2014-01-23 17:04:40 +0000
+++ Template/NLO/MCatNLO/include/LHEF.h	2014-11-06 06:25:51 +0000
@@ -1698,9 +1698,9 @@
     version = 1;
     if ( currentLine.find("version=\"2" ) != std::string::npos )
       version = 2;
-    else if ( currentLine.find("version=\"1" ) == std::string::npos )
-      throw std::runtime_error
-	("Tried to read a LesHouchesEvents file which is not version 1.0.");
+    //else if ( currentLine.find("version=\"1" ) == std::string::npos )
+    //  throw std::runtime_error
+    //("Tried to read a LesHouchesEvents file which is not version 1.0.");
 
     // Loop over all lines until we hit the </init> tag.
     while ( getline() && currentLine.find("</init>") == std::string::npos ) {

=== modified file 'Template/NLO/MCatNLO/shower_template.sh'
--- Template/NLO/MCatNLO/shower_template.sh	2014-06-24 13:08:32 +0000
+++ Template/NLO/MCatNLO/shower_template.sh	2014-11-06 06:25:51 +0000
@@ -35,7 +35,7 @@
     ./MCATNLO_$SHOWER\_EXE < MCATNLO_$SHOWER\_input > mcatnlo_run.log 2>&1
 
 elif [ $SHOWER == "PYTHIA8" ] ; then
-    source config.sh
+    if [ -f config.sh ] ; then source config.sh ; fi
     ./Pythia8.exe Pythia8.cmd > mcatnlo_run.log 2>&1
 fi
 

=== modified file 'Template/NLO/MCatNLO/srcPythia8/Makefile'
--- Template/NLO/MCatNLO/srcPythia8/Makefile	2014-01-23 17:04:40 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Makefile	2014-11-06 06:25:51 +0000
@@ -66,11 +66,19 @@
 Pythia8: $(FORFILE) $(PYTHIA8LOCATION)/$(LIBDIRARCH)/libpythia8.a \
 	$(PYTHIA8LOCATION)/$(LIBDIRARCH)/libpythia8tohepmc.a
 	$(CXX) $(CXXFLAGS) -I$(PYTHIA8LOCATION)/$(INCDIR) \
-        -I$(HEPMCLOCATION)/$(INCDIR) -I$(TOPDIR)/$(INCDIR) \
-        $(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o $@.exe $(FORFILE) \
+	-I$(HEPMCLOCATION)/$(INCDIR) -I$(TOPDIR)/$(INCDIR) \
+	$(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o Pythia8.exe $(FORFILE) \
 	-L$(PYTHIA8LOCATION)/$(LIBDIRARCH) -lpythia8 $(LIBLHAPDF) $(LIBGZIP) \
 	-lpythia8tohepmc -L$(HEPMCLOCATION)/$(LIBDIR) -lHepMC -lHepMCfio $(LIBGFORTRAN) \
 
+Pythia82: $(FORFILE) $(PYTHIA8LOCATION)/$(LIBDIR)/libpythia8.a
+	  $(CXX) $(CXXFLAGS) -I$(PYTHIA8LOCATION)/$(INCDIR) \
+	  $(HEPMCINCLIB) -I$(TOPDIR)/$(INCDIR) \
+	  -I$(TOPDIR)/$(INCDIR) $(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o Pythia8.exe $(FORFILE) \
+	  -L$(PYTHIA8LOCATION)/$(LIBDIR) -lpythia8 $(LIBGZIP) \
+	  -lHepMCfio $(LIBGFORTRAN) \
+
+
 .PHONY: clean
 clean:
 	rm -rf *.exe

=== modified file 'Template/NLO/MCatNLO/srcPythia8/Makefile_hep'
--- Template/NLO/MCatNLO/srcPythia8/Makefile_hep	2014-01-23 17:04:40 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Makefile_hep	2014-11-06 06:25:51 +0000
@@ -64,12 +64,18 @@
 	$(CXX) -c $< -fPIC $(INCLOPTION) -o $@
 
 Pythia8: $(FORFILE) $(PYTHIA8LOCATION)/$(LIBDIRARCH)/libpythia8.a \
-	$(PYTHIA8LOCATION)/$(LIBDIRARCH)/libpythia8tohepmc.a
-	$(CXX) $(CXXFLAGS) -I$(PYTHIA8LOCATION)/$(INCDIR) \
-        -I$(HEPMCLOCATION)/$(INCDIR) -I$(TOPDIR)/$(INCDIR) \
-        $(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o $@.exe $(FORFILE) \
-	-L$(PYTHIA8LOCATION)/$(LIBDIRARCH) -lpythia8 $(LIBLHAPDF) $(LIBGZIP) \
-	-lpythia8tohepmc -L$(HEPMCLOCATION)/$(LIBDIR) -lHepMC
+	  $(PYTHIA8LOCATION)/$(LIBDIRARCH)/libpythia8tohepmc.a
+	  $(CXX) $(CXXFLAGS) -I$(PYTHIA8LOCATION)/$(INCDIR) \
+	  -I$(HEPMCLOCATION)/$(INCDIR) -I$(TOPDIR)/$(INCDIR) \
+	  $(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o Pythia8.exe $(FORFILE) \
+	  -L$(PYTHIA8LOCATION)/$(LIBDIRARCH) -lpythia8 $(LIBLHAPDF) $(LIBGZIP) \
+	  -lpythia8tohepmc -L$(HEPMCLOCATION)/$(LIBDIR) -lHepMC
+
+Pythia82: $(FORFILE) $(PYTHIA8LOCATION)/$(LIBDIR)/libpythia8.a
+	  $(CXX) $(CXXFLAGS) -I$(PYTHIA8LOCATION)/$(INCDIR) \
+	  $(HEPMCINCLIB) -I$(TOPDIR)/$(INCDIR) \
+	  $(EXTRAPATHS) $(EXTRALIBS) $(INCLOPTION) $@.cc -o Pythia8.exe $(FORFILE) \
+	  -L$(PYTHIA8LOCATION)/$(LIBDIR) -lpythia8 $(LIBGZIP)
 
 .PHONY: clean
 clean:

=== added file 'Template/NLO/MCatNLO/srcPythia8/Pythia8.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia8.cc	1970-01-01 00:00:00 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia8.cc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,152 @@
+// Driver for Pythia 8. Reads an input file dynamically created on
+// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
+#include "Pythia8/Pythia.h"
+#include "Pythia8/Pythia8ToHepMC.h"
+#include "HepMC/GenEvent.h"
+#include "HepMC/IO_GenEvent.h"
+#include "HepMC/IO_BaseClass.h"
+#include "HepMC/IO_HEPEVT.h"
+#include "HepMC/HEPEVT_Wrapper.h"
+#include "fstream"
+#include "LHEFRead.h"
+#include "../examples/CombineMatchingInput.h"
+
+using namespace Pythia8;
+
+extern "C" {
+  extern struct {
+    double EVWGT;
+  } cevwgt_;
+}
+#define cevwgt cevwgt_
+
+extern "C" { 
+  void pyabeg_(int&,char(*)[15]);
+  void pyaend_(double&);
+  void pyanal_(int&,double(*));
+}
+
+int main() {
+  Pythia pythia;
+
+  int cwgtinfo_nn;
+  char cwgtinfo_weights_info[250][15];
+  double cwgt_ww[250];
+
+  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
+
+  pythia.readFile(inputname.c_str());
+
+  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
+  CombineMatchingInput combined;
+  UserHooks* matching = combined.getHook(pythia);
+  if (!matching) return 1;
+  pythia.setUserHooksPtr(matching);
+
+  pythia.init();
+  string filename = pythia.word("Beams:LHEF");
+
+  MyReader read(filename);
+  read.lhef_read_wgtsinfo_(cwgtinfo_nn,cwgtinfo_weights_info);
+  pyabeg_(cwgtinfo_nn,cwgtinfo_weights_info);
+
+  int nAbort=10;
+  int nPrintLHA=1;
+  int iAbort=0;
+  int iPrintLHA=0;
+  int nstep=5000;
+  int iEventtot=pythia.mode("Main:numberOfEvents");
+  int iEventshower=pythia.mode("Main:spareMode1");
+  string evt_norm=pythia.word("Main:spareWord1");
+  int iEventtot_norm=iEventtot;
+  if (evt_norm == "average"){
+    iEventtot_norm=1;
+  }
+
+  //FxFx merging
+  bool isFxFx=pythia.flag("JetMatching:doFxFx");
+  if (isFxFx) {
+    int nJmax=pythia.mode("JetMatching:nJetMax");
+    double Qcut=pythia.parm("JetMatching:qCut");
+    double PTcut=pythia.parm("JetMatching:qCutME");
+    if (Qcut <= PTcut || Qcut <= 0.) {
+      std::cout << " \n";
+      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
+		<< Qcut << " " << PTcut << "\n";
+      return 0;
+    }
+  }
+
+  HepMC::IO_BaseClass *_hepevtio;
+  HepMC::Pythia8ToHepMC ToHepMC;
+  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
+  double nSelected;
+  double norm;
+
+  // Cross section
+  double sigmaTotal  = 0.;
+
+  for (int iEvent = 0; ; ++iEvent) {
+    if (!pythia.next()) {
+      if (++iAbort < nAbort) continue;
+      break;
+    }
+    // the number of events read by Pythia so far
+    nSelected=double(pythia.info.nSelected());
+    // normalisation factor for the default analyses defined in pyanal_
+    norm=iEventtot_norm*iEvent/nSelected;
+
+    if (nSelected >= iEventshower) break;
+    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
+      pythia.LHAeventList();
+      pythia.info.list();
+      pythia.process.list();
+      pythia.event.list();
+      ++iPrintLHA;
+    }
+
+    double evtweight = pythia.info.weight();
+    double normhepmc;
+    // Add the weight of the current event to the cross section.
+    normhepmc = 1. / double(iEventshower);
+    if (evt_norm == "average") {
+      sigmaTotal  += evtweight*normhepmc;
+    } else {
+      sigmaTotal  += evtweight*normhepmc*iEventtot;
+    }
+
+    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
+    ToHepMC.fill_next_event( pythia, hepmcevt );
+
+    //define the IO_HEPEVT
+    _hepevtio = new HepMC::IO_HEPEVT;
+    _hepevtio->write_event(hepmcevt);
+    
+    //event weight
+    cevwgt.EVWGT=hepmcevt->weights()[0];
+
+    //call the FORTRAN analysis for this event
+    read.lhef_read_wgts_(cwgt_ww);
+    pyanal_(cwgtinfo_nn,cwgt_ww);
+
+    if (iEvent % nstep == 0 && iEvent >= 100){
+      pyaend_(norm);
+    }
+    delete hepmcevt;
+  }
+  pyaend_(norm);
+
+  pythia.stat();
+  if (isFxFx){
+    std::cout << " \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "Cross section, including FxFx merging is: "
+	      << sigmaTotal << "\n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+  }
+  delete matching;
+
+  return 0;
+}

=== removed file 'Template/NLO/MCatNLO/srcPythia8/Pythia8.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia8.cc	2014-05-20 11:47:02 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia8.cc	1970-01-01 00:00:00 +0000
@@ -1,152 +0,0 @@
-// Driver for Pythia 8. Reads an input file dynamically created on
-// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
-#include "Pythia8/Pythia.h"
-#include "Pythia8/Pythia8ToHepMC.h"
-#include "HepMC/GenEvent.h"
-#include "HepMC/IO_GenEvent.h"
-#include "HepMC/IO_BaseClass.h"
-#include "HepMC/IO_HEPEVT.h"
-#include "HepMC/HEPEVT_Wrapper.h"
-#include "fstream"
-#include "LHEFRead.h"
-#include "../examples/CombineMatchingInput.h"
-
-using namespace Pythia8;
-
-extern "C" {
-  extern struct {
-    double EVWGT;
-  } cevwgt_;
-}
-#define cevwgt cevwgt_
-
-extern "C" { 
-  void pyabeg_(int&,char(*)[15]);
-  void pyaend_(double&);
-  void pyanal_(int&,double(*));
-}
-
-int main() {
-  Pythia pythia;
-
-  int cwgtinfo_nn;
-  char cwgtinfo_weights_info[250][15];
-  double cwgt_ww[250];
-
-  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
-
-  pythia.readFile(inputname.c_str());
-
-  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
-  CombineMatchingInput combined;
-  UserHooks* matching = combined.getHook(pythia);
-  if (!matching) return 1;
-  pythia.setUserHooksPtr(matching);
-
-  pythia.init();
-  string filename = pythia.word("Beams:LHEF");
-
-  MyReader read(filename);
-  read.lhef_read_wgtsinfo_(cwgtinfo_nn,cwgtinfo_weights_info);
-  pyabeg_(cwgtinfo_nn,cwgtinfo_weights_info);
-
-  int nAbort=10;
-  int nPrintLHA=1;
-  int iAbort=0;
-  int iPrintLHA=0;
-  int nstep=5000;
-  int iEventtot=pythia.mode("Main:numberOfEvents");
-  int iEventshower=pythia.mode("Main:spareMode1");
-  string evt_norm=pythia.word("Main:spareWord1");
-  int iEventtot_norm=iEventtot;
-  if (evt_norm == "average"){
-    iEventtot_norm=1;
-  }
-
-  //FxFx merging
-  bool isFxFx=pythia.flag("JetMatching:doFxFx");
-  if (isFxFx) {
-    int nJmax=pythia.mode("JetMatching:nJetMax");
-    double Qcut=pythia.parm("JetMatching:qCut");
-    double PTcut=pythia.parm("JetMatching:qCutME");
-    if (Qcut <= PTcut || Qcut <= 0.) {
-      std::cout << " \n";
-      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
-		<< Qcut << " " << PTcut << "\n";
-      return 0;
-    }
-  }
-
-  HepMC::IO_BaseClass *_hepevtio;
-  HepMC::Pythia8ToHepMC ToHepMC;
-  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
-  double nSelected;
-  double norm;
-
-  // Cross section
-  double sigmaTotal  = 0.;
-
-  for (int iEvent = 0; ; ++iEvent) {
-    if (!pythia.next()) {
-      if (++iAbort < nAbort) continue;
-      break;
-    }
-    // the number of events read by Pythia so far
-    nSelected=double(pythia.info.nSelected());
-    // normalisation factor for the default analyses defined in pyanal_
-    norm=iEventtot_norm*iEvent/nSelected;
-
-    if (nSelected >= iEventshower) break;
-    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
-      pythia.LHAeventList();
-      pythia.info.list();
-      pythia.process.list();
-      pythia.event.list();
-      ++iPrintLHA;
-    }
-
-    double evtweight = pythia.info.weight();
-    double normhepmc;
-    // Add the weight of the current event to the cross section.
-    normhepmc = 1. / double(iEventshower);
-    if (evt_norm == "average") {
-      sigmaTotal  += evtweight*normhepmc;
-    } else {
-      sigmaTotal  += evtweight*normhepmc*iEventtot;
-    }
-
-    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
-    ToHepMC.fill_next_event( pythia, hepmcevt );
-
-    //define the IO_HEPEVT
-    _hepevtio = new HepMC::IO_HEPEVT;
-    _hepevtio->write_event(hepmcevt);
-    
-    //event weight
-    cevwgt.EVWGT=hepmcevt->weights()[0];
-
-    //call the FORTRAN analysis for this event
-    read.lhef_read_wgts_(cwgt_ww);
-    pyanal_(cwgtinfo_nn,cwgt_ww);
-
-    if (iEvent % nstep == 0 && iEvent >= 100){
-      pyaend_(norm);
-    }
-    delete hepmcevt;
-  }
-  pyaend_(norm);
-
-  pythia.stat();
-  if (isFxFx){
-    std::cout << " \n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "Cross section, including FxFx merging is: "
-	      << sigmaTotal << "\n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "*********************************************************************** \n";
-  }
-  delete matching;
-
-  return 0;
-}

=== added file 'Template/NLO/MCatNLO/srcPythia8/Pythia82.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia82.cc	1970-01-01 00:00:00 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia82.cc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,153 @@
+// Driver for Pythia 8. Reads an input file dynamically created on
+// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
+#include "Pythia8/Pythia.h"
+#include "Pythia8Plugins/HepMC2.h"
+#include "Pythia8Plugins/aMCatNLOHooks.h"
+#include "Pythia8Plugins/CombineMatchingInput.h"
+#include "HepMC/GenEvent.h"
+#include "HepMC/IO_GenEvent.h"
+#include "HepMC/IO_BaseClass.h"
+#include "HepMC/IO_HEPEVT.h"
+#include "HepMC/HEPEVT_Wrapper.h"
+#include "fstream"
+#include "LHEFRead.h"
+
+using namespace Pythia8;
+
+extern "C" {
+  extern struct {
+    double EVWGT;
+  } cevwgt_;
+}
+#define cevwgt cevwgt_
+
+extern "C" { 
+  void pyabeg_(int&,char(*)[15]);
+  void pyaend_(double&);
+  void pyanal_(int&,double(*));
+}
+
+int main() {
+  Pythia pythia;
+
+  int cwgtinfo_nn;
+  char cwgtinfo_weights_info[250][15];
+  double cwgt_ww[250];
+
+  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
+
+  pythia.readFile(inputname.c_str());
+
+  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
+  CombineMatchingInput combined;
+  UserHooks* matching = combined.getHook(pythia);
+  if (!matching) return 1;
+  pythia.setUserHooksPtr(matching);
+
+  pythia.init();
+  string filename = pythia.word("Beams:LHEF");
+
+  MyReader read(filename);
+  read.lhef_read_wgtsinfo_(cwgtinfo_nn,cwgtinfo_weights_info);
+  pyabeg_(cwgtinfo_nn,cwgtinfo_weights_info);
+
+  int nAbort=10;
+  int nPrintLHA=1;
+  int iAbort=0;
+  int iPrintLHA=0;
+  int nstep=5000;
+  int iEventtot=pythia.mode("Main:numberOfEvents");
+  int iEventshower=pythia.mode("Main:spareMode1");
+  string evt_norm=pythia.word("Main:spareWord1");
+  int iEventtot_norm=iEventtot;
+  if (evt_norm == "average"){
+    iEventtot_norm=1;
+  }
+
+  //FxFx merging
+  bool isFxFx=pythia.flag("JetMatching:doFxFx");
+  if (isFxFx) {
+    int nJmax=pythia.mode("JetMatching:nJetMax");
+    double Qcut=pythia.parm("JetMatching:qCut");
+    double PTcut=pythia.parm("JetMatching:qCutME");
+    if (Qcut <= PTcut || Qcut <= 0.) {
+      std::cout << " \n";
+      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
+		<< Qcut << " " << PTcut << "\n";
+      return 0;
+    }
+  }
+
+  HepMC::IO_BaseClass *_hepevtio;
+  HepMC::Pythia8ToHepMC ToHepMC;
+  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
+  double nSelected;
+  double norm;
+
+  // Cross section
+  double sigmaTotal  = 0.;
+
+  for (int iEvent = 0; ; ++iEvent) {
+    if (!pythia.next()) {
+      if (++iAbort < nAbort) continue;
+      break;
+    }
+    // the number of events read by Pythia so far
+    nSelected=double(pythia.info.nSelected());
+    // normalisation factor for the default analyses defined in pyanal_
+    norm=iEventtot_norm*iEvent/nSelected;
+
+    if (nSelected >= iEventshower) break;
+    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
+      pythia.LHAeventList();
+      pythia.info.list();
+      pythia.process.list();
+      pythia.event.list();
+      ++iPrintLHA;
+    }
+
+    double evtweight = pythia.info.weight();
+    double normhepmc;
+    // Add the weight of the current event to the cross section.
+    normhepmc = 1. / double(iEventshower);
+    if (evt_norm == "average") {
+      sigmaTotal  += evtweight*normhepmc;
+    } else {
+      sigmaTotal  += evtweight*normhepmc*iEventtot;
+    }
+
+    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
+    ToHepMC.fill_next_event( pythia, hepmcevt );
+
+    //define the IO_HEPEVT
+    _hepevtio = new HepMC::IO_HEPEVT;
+    _hepevtio->write_event(hepmcevt);
+    
+    //event weight
+    cevwgt.EVWGT=hepmcevt->weights()[0];
+
+    //call the FORTRAN analysis for this event
+    read.lhef_read_wgts_(cwgt_ww);
+    pyanal_(cwgtinfo_nn,cwgt_ww);
+
+    if (iEvent % nstep == 0 && iEvent >= 100){
+      pyaend_(norm);
+    }
+    delete hepmcevt;
+  }
+  pyaend_(norm);
+
+  pythia.stat();
+  if (isFxFx){
+    std::cout << " \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "Cross section, including FxFx merging is: "
+	      << sigmaTotal << "\n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+  }
+  delete matching;
+
+  return 0;
+}

=== added file 'Template/NLO/MCatNLO/srcPythia8/Pythia82_hep.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia82_hep.cc	1970-01-01 00:00:00 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia82_hep.cc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,112 @@
+// Driver for Pythia 8. Reads an input file dynamically created on
+// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
+#include "Pythia8/Pythia.h"
+#include "Pythia8Plugins/HepMC2.h"
+#include "Pythia8Plugins/aMCatNLOHooks.h"
+#include "Pythia8Plugins/CombineMatchingInput.h"
+#include "HepMC/GenEvent.h"
+#include "HepMC/IO_GenEvent.h"
+
+using namespace Pythia8;
+
+int main() {
+  Pythia pythia;
+
+  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
+
+  pythia.readFile(inputname.c_str());
+
+  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
+  CombineMatchingInput combined;
+  UserHooks* matching = combined.getHook(pythia);
+  if (!matching) return 1;
+  pythia.setUserHooksPtr(matching);
+
+  pythia.init();
+
+  int nAbort=10;
+  int nPrintLHA=1;
+  int iAbort=0;
+  int iPrintLHA=0;
+  int iEventtot=pythia.mode("Main:numberOfEvents");
+  int iEventshower=pythia.mode("Main:spareMode1");
+  string evt_norm=pythia.word("Main:spareWord1");
+
+  //FxFx merging
+  bool isFxFx=pythia.flag("JetMatching:doFxFx");
+  if (isFxFx) {
+    int nJmax=pythia.mode("JetMatching:nJetMax");
+    double Qcut=pythia.parm("JetMatching:qCut");
+    double PTcut=pythia.parm("JetMatching:qCutME");
+    if (Qcut <= PTcut || Qcut <= 0.) {
+      std::cout << " \n";
+      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
+		<< Qcut << " " << PTcut << "\n";
+      return 0;
+    }
+  }
+
+  HepMC::Pythia8ToHepMC ToHepMC;
+  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
+  // Do not store cross section information, as this will be done manually.
+  ToHepMC.set_store_pdf(false);
+  ToHepMC.set_store_proc(false);
+  ToHepMC.set_store_xsec(false);
+
+  // Cross section an error.
+  double sigmaTotal  = 0.;
+  double errorTotal  = 0.;
+
+  for (int iEvent = 0; ; ++iEvent) {
+    if (!pythia.next()) {
+      if (++iAbort < nAbort) continue;
+      break;
+    }
+    // the number of events read by Pythia so far
+    int nSelected=pythia.info.nSelected();
+
+    if (nSelected >= iEventshower) break;
+    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
+      pythia.LHAeventList();
+      pythia.info.list();
+      pythia.process.list();
+      pythia.event.list();
+      ++iPrintLHA;
+    }
+
+    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
+    double evtweight = pythia.info.weight();
+    double normhepmc;
+    // ALWAYS NORMALISE HEPMC WEIGHTS TO SUM TO THE CROSS SECTION
+    if (evt_norm == "average") {
+      normhepmc = 1. / double(iEventshower);
+    } else {
+      normhepmc = double(iEventtot) / double(iEventshower);
+    }
+    hepmcevt->weights().push_back(evtweight*normhepmc);
+    ToHepMC.fill_next_event( pythia, hepmcevt );
+    // Add the weight of the current event to the cross section.
+    // Report cross section to hepmc
+    HepMC::GenCrossSection xsec;
+    xsec.set_cross_section( sigmaTotal, pythia.info.sigmaErr() );
+    hepmcevt->set_cross_section( xsec );
+    // Write the HepMC event to file. Done with it.
+    ascii_io << hepmcevt;    
+    delete hepmcevt;
+  }
+
+  pythia.stat();
+  if (isFxFx){
+    std::cout << " \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "Cross section, including FxFx merging is: "
+	      << sigmaTotal << "\n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+  }
+
+  delete matching;
+
+  return 0;
+}

=== added file 'Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc	1970-01-01 00:00:00 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,111 @@
+// Driver for Pythia 8. Reads an input file dynamically created on
+// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
+#include "Pythia8/Pythia.h"
+#include "Pythia8/Pythia8ToHepMC.h"
+#include "HepMC/GenEvent.h"
+#include "HepMC/IO_GenEvent.h"
+#include "../examples/CombineMatchingInput.h"
+
+using namespace Pythia8;
+
+int main() {
+  Pythia pythia;
+
+  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
+
+  pythia.readFile(inputname.c_str());
+
+  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
+  CombineMatchingInput combined;
+  UserHooks* matching = combined.getHook(pythia);
+  if (!matching) return 1;
+  pythia.setUserHooksPtr(matching);
+
+  pythia.init();
+
+  int nAbort=10;
+  int nPrintLHA=1;
+  int iAbort=0;
+  int iPrintLHA=0;
+  int iEventtot=pythia.mode("Main:numberOfEvents");
+  int iEventshower=pythia.mode("Main:spareMode1");
+  string evt_norm=pythia.word("Main:spareWord1");
+
+  //FxFx merging
+  bool isFxFx=pythia.flag("JetMatching:doFxFx");
+  if (isFxFx) {
+    int nJmax=pythia.mode("JetMatching:nJetMax");
+    double Qcut=pythia.parm("JetMatching:qCut");
+    double PTcut=pythia.parm("JetMatching:qCutME");
+    if (Qcut <= PTcut || Qcut <= 0.) {
+      std::cout << " \n";
+      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
+		<< Qcut << " " << PTcut << "\n";
+      return 0;
+    }
+  }
+
+  HepMC::Pythia8ToHepMC ToHepMC;
+  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
+  // Do not store cross section information, as this will be done manually.
+  ToHepMC.set_store_pdf(false);
+  ToHepMC.set_store_proc(false);
+  ToHepMC.set_store_xsec(false);
+
+  // Cross section an error.
+  double sigmaTotal  = 0.;
+  double errorTotal  = 0.;
+
+  for (int iEvent = 0; ; ++iEvent) {
+    if (!pythia.next()) {
+      if (++iAbort < nAbort) continue;
+      break;
+    }
+    // the number of events read by Pythia so far
+    int nSelected=pythia.info.nSelected();
+
+    if (nSelected >= iEventshower) break;
+    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
+      pythia.LHAeventList();
+      pythia.info.list();
+      pythia.process.list();
+      pythia.event.list();
+      ++iPrintLHA;
+    }
+
+    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
+    double evtweight = pythia.info.weight();
+    double normhepmc;
+    // ALWAYS NORMALISE HEPMC WEIGHTS TO SUM TO THE CROSS SECTION
+    if (evt_norm == "average") {
+      normhepmc = 1. / double(iEventshower);
+    } else {
+      normhepmc = double(iEventtot) / double(iEventshower);
+    }
+    hepmcevt->weights().push_back(evtweight*normhepmc);
+    ToHepMC.fill_next_event( pythia, hepmcevt );
+    // Add the weight of the current event to the cross section.
+    // Report cross section to hepmc
+    HepMC::GenCrossSection xsec;
+    xsec.set_cross_section( sigmaTotal, pythia.info.sigmaErr() );
+    hepmcevt->set_cross_section( xsec );
+    // Write the HepMC event to file. Done with it.
+    ascii_io << hepmcevt;    
+    delete hepmcevt;
+  }
+
+  pythia.stat();
+  if (isFxFx){
+    std::cout << " \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "Cross section, including FxFx merging is: "
+	      << sigmaTotal << "\n";
+    std::cout << "*********************************************************************** \n";
+    std::cout << "*********************************************************************** \n";
+  }
+
+  delete matching;
+
+  return 0;
+}

=== removed file 'Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc'
--- Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc	2014-05-20 11:47:02 +0000
+++ Template/NLO/MCatNLO/srcPythia8/Pythia8_hep.cc	1970-01-01 00:00:00 +0000
@@ -1,111 +0,0 @@
-// Driver for Pythia 8. Reads an input file dynamically created on
-// the basis of the inputs specified in MCatNLO_MadFKS_PY8.Script 
-#include "Pythia8/Pythia.h"
-#include "Pythia8/Pythia8ToHepMC.h"
-#include "HepMC/GenEvent.h"
-#include "HepMC/IO_GenEvent.h"
-#include "../examples/CombineMatchingInput.h"
-
-using namespace Pythia8;
-
-int main() {
-  Pythia pythia;
-
-  string inputname="Pythia8.cmd",outputname="Pythia8.hep";
-
-  pythia.readFile(inputname.c_str());
-
-  //Create UserHooks pointer for the FxFX matching. Stop if it failed. Pass pointer to Pythia.
-  CombineMatchingInput combined;
-  UserHooks* matching = combined.getHook(pythia);
-  if (!matching) return 1;
-  pythia.setUserHooksPtr(matching);
-
-  pythia.init();
-
-  int nAbort=10;
-  int nPrintLHA=1;
-  int iAbort=0;
-  int iPrintLHA=0;
-  int iEventtot=pythia.mode("Main:numberOfEvents");
-  int iEventshower=pythia.mode("Main:spareMode1");
-  string evt_norm=pythia.word("Main:spareWord1");
-
-  //FxFx merging
-  bool isFxFx=pythia.flag("JetMatching:doFxFx");
-  if (isFxFx) {
-    int nJmax=pythia.mode("JetMatching:nJetMax");
-    double Qcut=pythia.parm("JetMatching:qCut");
-    double PTcut=pythia.parm("JetMatching:qCutME");
-    if (Qcut <= PTcut || Qcut <= 0.) {
-      std::cout << " \n";
-      std::cout << "Merging scale (shower_card.dat) smaller than pTcut (run_card.dat)"
-		<< Qcut << " " << PTcut << "\n";
-      return 0;
-    }
-  }
-
-  HepMC::Pythia8ToHepMC ToHepMC;
-  HepMC::IO_GenEvent ascii_io(outputname.c_str(), std::ios::out);
-  // Do not store cross section information, as this will be done manually.
-  ToHepMC.set_store_pdf(false);
-  ToHepMC.set_store_proc(false);
-  ToHepMC.set_store_xsec(false);
-
-  // Cross section an error.
-  double sigmaTotal  = 0.;
-  double errorTotal  = 0.;
-
-  for (int iEvent = 0; ; ++iEvent) {
-    if (!pythia.next()) {
-      if (++iAbort < nAbort) continue;
-      break;
-    }
-    // the number of events read by Pythia so far
-    int nSelected=pythia.info.nSelected();
-
-    if (nSelected >= iEventshower) break;
-    if (pythia.info.isLHA() && iPrintLHA < nPrintLHA) {
-      pythia.LHAeventList();
-      pythia.info.list();
-      pythia.process.list();
-      pythia.event.list();
-      ++iPrintLHA;
-    }
-
-    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
-    double evtweight = pythia.info.weight();
-    double normhepmc;
-    // ALWAYS NORMALISE HEPMC WEIGHTS TO SUM TO THE CROSS SECTION
-    if (evt_norm == "average") {
-      normhepmc = 1. / double(iEventshower);
-    } else {
-      normhepmc = double(iEventtot) / double(iEventshower);
-    }
-    hepmcevt->weights().push_back(evtweight*normhepmc);
-    ToHepMC.fill_next_event( pythia, hepmcevt );
-    // Add the weight of the current event to the cross section.
-    // Report cross section to hepmc
-    HepMC::GenCrossSection xsec;
-    xsec.set_cross_section( sigmaTotal, pythia.info.sigmaErr() );
-    hepmcevt->set_cross_section( xsec );
-    // Write the HepMC event to file. Done with it.
-    ascii_io << hepmcevt;    
-    delete hepmcevt;
-  }
-
-  pythia.stat();
-  if (isFxFx){
-    std::cout << " \n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "Cross section, including FxFx merging is: "
-	      << sigmaTotal << "\n";
-    std::cout << "*********************************************************************** \n";
-    std::cout << "*********************************************************************** \n";
-  }
-
-  delete matching;
-
-  return 0;
-}

=== modified file 'Template/NLO/Source/PDF/pdf.f'
--- Template/NLO/Source/PDF/pdf.f	2014-05-16 08:09:51 +0000
+++ Template/NLO/Source/PDF/pdf.f	2014-11-06 06:25:51 +0000
@@ -45,6 +45,7 @@
          do i=-5,5
             fx(i)=nnfx(i)/x
          enddo
+         fx(7)=nnfx(7)/x
       elseif     ((pdlabel(1:3) .eq. 'mrs')
      .   .or. (pdlabel(2:4) .eq. 'mrs')) then
 

=== modified file 'Template/NLO/Source/banner_header.txt'
--- Template/NLO/Source/banner_header.txt	2013-11-30 09:22:32 +0000
+++ Template/NLO/Source/banner_header.txt	2014-11-06 06:25:51 +0000
@@ -1,4 +1,4 @@
-<LesHouchesEvents version="1.0">
+<LesHouchesEvents version="3.0">
 <header>
 <!--
 #*********************************************************************

=== modified file 'Template/NLO/Source/make_opts.inc'
--- Template/NLO/Source/make_opts.inc	2014-09-02 12:20:56 +0000
+++ Template/NLO/Source/make_opts.inc	2014-11-06 06:25:51 +0000
@@ -16,8 +16,7 @@
   CFLAGS = -O 
 endif
 
-  
-
+# Increase the number of allowed charcters in a Fortran line
 FFLAGS+= -ffixed-line-length-132
 
 # Set FC unless it's defined by an environment variable

=== modified file 'Template/NLO/Source/setrun.f'
--- Template/NLO/Source/setrun.f	2014-09-11 05:33:30 +0000
+++ Template/NLO/Source/setrun.f	2014-11-06 06:25:51 +0000
@@ -56,7 +56,7 @@
      &     xmaxup(maxpup),lprup(maxpup)
 c
       include 'nexternal.inc'
-      include 'leshouche_info.inc'
+      include 'leshouche_decl.inc'
       logical gridrun,gridpack
       integer          iseed
       common /to_seed/ iseed
@@ -73,7 +73,8 @@
       
 c MZ add the possibility to have shower_MC input lowercase
       call to_upper(shower_MC)
-
+C
+      call read_leshouche_info(idup_d,mothup_d,icolup_d)
 
 c merging cuts
       xqcut=0d0

=== modified file 'Template/NLO/SubProcesses/FKSParamReader.f'
--- Template/NLO/SubProcesses/FKSParamReader.f	2013-11-27 09:02:57 +0000
+++ Template/NLO/SubProcesses/FKSParamReader.f	2014-11-06 06:25:51 +0000
@@ -6,7 +6,8 @@
       data HasReadOnce/.FALSE./
 
 
-      CHARACTER*64 fileName, buff, buff2, mode
+      character(*) filename
+      CHARACTER*64 buff, buff2, mode
       include "FKSParams.inc"
 
       logical printParam, couldRead, paramPrinted
@@ -18,36 +19,36 @@
           goto 901
       endif
 
-      open(666, file=fileName, err=676, action='READ')
+      open(68, file=fileName, err=676, action='READ')
       do
-         read(666,*,end=999) buff
+         read(68,*,end=999) buff
          if(index(buff,'#').eq.1) then
 
            if (buff .eq. '#IRPoleCheckThreshold') then
-             read(666,*,end=999) IRPoleCheckThreshold
+             read(68,*,end=999) IRPoleCheckThreshold
              if (IRPoleCheckThreshold .lt. -1.01d0 ) then
                stop 'IRPoleCheckThreshold must be >= -1.0d0.'
              endif 
 
            elseif (buff .eq. '#PrecisionVirtualAtRunTime') then
-             read(666,*,end=999) PrecisionVirtualAtRunTime
+             read(68,*,end=999) PrecisionVirtualAtRunTime
              if (IRPoleCheckThreshold .lt. -1.01d0 ) then
                stop 'PrecisionVirtualAtRunTime must be >= -1.0d0.'
              endif 
 
            else if (buff .eq. '#NHelForMCoverHels') then
-             read(666,*,end=999) NHelForMCoverHels
+             read(68,*,end=999) NHelForMCoverHels
              if (NHelForMCoverHels .lt. -1) then
                stop 'NHelForMCoverHels must be >= -1.'
              endif 
            else if (buff .eq. '#VirtualFraction') then
-             read(666,*,end=999) Virt_fraction
+             read(68,*,end=999) Virt_fraction
              if (Virt_fraction .lt. 0 .or. virt_fraction .gt.1) then
                 stop 'VirtualFraction should be a fraction'/
      $               /' between 0 and 1'
              endif 
            else if (buff .eq. '#MinVirtualFraction') then
-             read(666,*,end=999) Min_Virt_fraction
+             read(68,*,end=999) Min_Virt_fraction
              if (min_virt_fraction .lt. 0 .or. min_virt_fraction .gt.1)
      $            then
                 stop 'VirtualFraction should be a fraction'/
@@ -100,7 +101,7 @@
       paramPrinted=.TRUE.
       endif
 
-      close(666)
+      close(68)
       HasReadOnce=.TRUE.
   901 continue
       end

=== modified file 'Template/NLO/SubProcesses/add_write_info.f'
--- Template/NLO/SubProcesses/add_write_info.f	2014-08-05 13:04:11 +0000
+++ Template/NLO/SubProcesses/add_write_info.f	2014-11-06 06:25:51 +0000
@@ -9,7 +9,7 @@
       include "coloramps.inc"
       include "reweight0.inc"
       include "nFKSconfigs.inc"
-      include "leshouche_info.inc"
+      include "leshouche_decl.inc"
       include "run.inc"
 
 c Arguments

=== modified file 'Template/NLO/SubProcesses/ajob_template'
--- Template/NLO/SubProcesses/ajob_template	2014-09-12 10:53:25 +0000
+++ Template/NLO/SubProcesses/ajob_template	2014-11-06 06:25:51 +0000
@@ -8,7 +8,7 @@
 }
 
 tarCounter=0
-while [[ (-e MadLoop5_resources.tar.gz) && (! -e MadLoop5_resources/HelConfigs.dat) && ($tarCounter < 10) ]]; do
+while [[ (-f MadLoop5_resources.tar.gz) && (! -f MadLoop5_resources/HelConfigs.dat) && ($tarCounter < 10) ]]; do
     if [[ $tarCounter > 0 ]]; then
 	sleep 2s
     fi
@@ -90,6 +90,8 @@
     fi
 
     link1up FKS_params.dat
+    link1up configs_and_props_info.dat
+    link1up leshouche_info.dat
 #   Not necessary anymore
 #    link1up MadLoop5_resources
     link1up OLE_order.olc

=== modified file 'Template/NLO/SubProcesses/check_poles.f'
--- Template/NLO/SubProcesses/check_poles.f	2014-09-24 01:01:34 +0000
+++ Template/NLO/SubProcesses/check_poles.f	2014-11-06 06:25:51 +0000
@@ -15,7 +15,6 @@
       double precision tolerance, tolerance_default
       double precision, allocatable :: accuracies(:)
       double precision accuracy
-      parameter (tolerance_default = 1d-5)
       double precision ren_scale, energy
       include 'genps.inc'
       include 'nexternal.inc'
@@ -54,6 +53,7 @@
       integer nfail
       logical first_time
       data first_time/.TRUE./
+      include 'FKSParams.inc'
       
 C-----
 C  BEGIN CODE
@@ -72,13 +72,15 @@
       call run_printout          !Prints out a summary of the run settings
       include 'pmass.inc'
      
+      call FKSParamReader('FKS_params.dat',.TRUE.,.FALSE.)
+      tolerance_default = IRPoleCheckThreshold
 
 c     Set the energy to be characteristic of the run
       totmass = 0.0d0
       do i=1,nexternal
         totmass = totmass + pmass(i)
       enddo
-      energy = max((ebeam(1)+ebeam(2))/2.0d0,2.0d0*totmass)
+      energy = max((ebeam(1)+ebeam(2))/20.0d0,2.0d0*totmass)
 c     Set the renormalization scale to be of the order of sqrt(s) but
 c     not equal to it so as to be sensitive to all logs in the check.
       ren_scale = energy/2.0d0

=== modified file 'Template/NLO/SubProcesses/configs_and_props_inc_chooser.f'
--- Template/NLO/SubProcesses/configs_and_props_inc_chooser.f	2013-01-11 08:13:45 +0000
+++ Template/NLO/SubProcesses/configs_and_props_inc_chooser.f	2014-11-06 06:25:51 +0000
@@ -22,7 +22,11 @@
       integer prow(-max_branch:-1,n_max_cg)
       common/c_props_inc/prmass,prwidth,prow
       double precision pmass(nexternal)
-      include 'configs_and_props_info.inc'
+      logical firsttime
+      data firsttime /.true./
+      include 'configs_and_props_decl.inc'
+      save mapconfig_d, iforest_d, sprop_d, tprid_d, pmass_d, pwidth_d
+     $ , pow_d
       include "pmass.inc"
 c     
       if (max_branch_used.gt.max_branch) then
@@ -35,6 +39,15 @@
      $        /' increase n_max_cg' ,n_max_cg,lmaxconfigs_used
          stop
       endif
+
+C the configurations and propagators infos are read at the first
+C evaluation
+      if (firsttime) then
+        call read_configs_and_props_info(mapconfig_d,iforest_d,sprop_d,
+     1                                   tprid_d,pmass_d,pwidth_d,pow_d)
+        firsttime = .false.
+      endif
+
 c
 c Fill the arrays of the c_configs_inc and c_props_inc common
 c blocks. Some of the information might not be available in the
@@ -65,3 +78,73 @@
       return
       end
 
+
+      subroutine read_configs_and_props_info(mapconfig_d,iforest_d,sprop_d,
+     1                                   tprid_d,pmass_d,pwidth_d,pow_d)
+C read the various information from the configs_and_props_info.dat file
+      implicit none
+      integer i,j,k
+      integer ndau, idau, dau, id
+      character *200 buff
+      double precision get_mass_from_id, get_width_from_id
+      include 'configs_and_props_decl.inc'
+
+      open(unit=78, file='configs_and_props_info.dat', status='old')
+      do while (.true.)
+        read(78,'(a)',end=999) buff
+        if (buff(:1).eq.'#') cycle
+        if (buff(:1).eq.'C') then
+        ! mapconfig
+        ! C  i   j   k -> MAPCONFIG_D(i,j)=k
+          read(buff(2:),*) i,j,k
+          mapconfig_d(i,j) = k
+        else if (buff(:1).eq.'F') then
+        ! iforest
+        ! after the first line there are as many lines
+        ! as the daughters
+        ! F  i   j   k  ndau
+        ! D dau_1
+        ! D ...
+        ! D dau_ndau        -> IFORREST_D(i,idau,i,k)=dau_idau
+          read(buff(2:),*) i,j,k,ndau
+          do idau=1,ndau
+            read(78,'(a)') buff
+            if (buff(:1).ne.'D') then
+              write(*,*) 'ERROR #1 in read_configs_and_props_info',
+     1                    i,j,k,ndau,buff
+              stop 
+            endif
+            read(buff(2:),*) dau
+            iforest_d(i,idau,j,k) = dau
+          enddo
+        else if (buff(:1).eq.'S') then
+        ! sprop
+        ! S  i   j   k  id -> SPROP_D(i,j,k)=id
+          read(buff(2:),*) i,j,k,id
+          sprop_d(i,j,k) = id
+        else if (buff(:1).eq.'T') then
+        ! tprid
+        ! T  i   j   k  id -> TPRID_D(i,j,k)=id
+          read(buff(2:),*) i,j,k,id
+          tprid_d(i,j,k) = id
+        else if (buff(:1).eq.'M') then
+        ! pmass and pwidth
+          read(buff(2:),*) i,j,k,id
+        ! M  i   j   k  id -> gives id of particle of which 
+        ! the mass/width is stored in PMASS/WIDTH_D(i,j,k)
+          pmass_d(i,j,k) = get_mass_from_id(id)
+          pwidth_d(i,j,k) = get_width_from_id(id)
+        else if (buff(:1).eq.'P') then
+        ! pow
+        ! P  i   j   k  id -> POW_D(i,j,k)=id
+          read(buff(2:),*) i,j,k,id
+          pow_d(i,j,k) = id
+        endif
+      enddo
+ 999  continue
+      close(78)
+
+      return 
+      end
+
+

=== modified file 'Template/NLO/SubProcesses/cuts.f'
--- Template/NLO/SubProcesses/cuts.f	2014-07-16 16:23:08 +0000
+++ Template/NLO/SubProcesses/cuts.f	2014-11-06 06:25:51 +0000
@@ -463,29 +463,6 @@
       END
 
 
-      subroutine unweight_function(p_born,unwgtfun)
-c This is a user-defined function to which to unweight the events
-c A non-flat distribution will generate events with a certain
-c weight. This is particularly useful to generate more events
-c (with smaller weight) in tails of distributions.
-c It computes the unwgt factor from the momenta and multiplies
-c the weight that goes into MINT (or vegas) with this factor.
-c Before writing out the events (or making the plots), this factor
-c is again divided out.
-c This function should be called with the Born momenta to be sure
-c that it stays the same for the events, counter-events, etc.
-c A value different from 1 makes that MINT (or vegas) does not list
-c the correct cross section.
-      implicit none
-      include 'nexternal.inc'
-      double precision unwgtfun,p_born(0:3,nexternal-1)
-
-      unwgtfun=1d0
-
-      return
-      end
-
-
       function chi_gamma_iso(dr,R0,xn,epsgamma,pTgamma)
 c Eq.(3.4) of Phys.Lett. B429 (1998) 369-374 [hep-ph/9801442]
       implicit none
@@ -909,3 +886,33 @@
       end
 
 
+      subroutine unweight_function(p_born,unwgtfun)
+c This is a user-defined function to which to unweight the events
+c A non-flat distribution will generate events with a certain
+c weight. This is particularly useful to generate more events
+c (with smaller weight) in tails of distributions.
+c It computes the unwgt factor from the momenta and multiplies
+c the weight that goes into MINT (or vegas) with this factor.
+c Before writing out the events (or making the plots), this factor
+c is again divided out.
+c This function should be called with the Born momenta to be sure
+c that it stays the same for the events, counter-events, etc.
+c A value different from 1 makes that MINT (or vegas) does not list
+c the correct cross section.
+      implicit none
+      include 'nexternal.inc'
+      double precision unwgtfun,p_born(0:3,nexternal-1),shat,sumdot
+      external sumdot
+
+      unwgtfun=1d0
+
+c How to enhance the tails is very process dependent. But, it is
+c probably easiest to enhance the tails using shat, e.g.:
+c      shat=sumdot(p_born(0,1),p_born(0,2),1d0)
+c      unwgtfun=max(100d0**2,shat)/100d0**2
+c      unwgtfun=unwgtfun**2
+
+      return
+      end
+
+

=== modified file 'Template/NLO/SubProcesses/fks_singular.f'
--- Template/NLO/SubProcesses/fks_singular.f	2014-09-11 05:33:30 +0000
+++ Template/NLO/SubProcesses/fks_singular.f	2014-11-06 06:25:51 +0000
@@ -955,7 +955,8 @@
      &              cnt_swgt*fkssymmetryfactor +
      &              (bsv_wgt-born_wgt)*fkssymmetryfactorBorn +
      &              deg_wgt*fkssymmetryfactorDeg +
-     &              deg_swgt*fkssymmetryfactorDeg )*vegaswgt + virt_wgt
+     &              deg_swgt*fkssymmetryfactorDeg )*vegaswgt +
+     &              virt_wgt/unwgtfun
          if(abs(plot_wgt).gt.1.d-20) then
             if(iplot.eq.-3)then
                write(*,*)'Error #1 in dsig'

=== modified file 'Template/NLO/SubProcesses/handling_lhe_events.f'
--- Template/NLO/SubProcesses/handling_lhe_events.f	2014-06-17 08:36:10 +0000
+++ Template/NLO/SubProcesses/handling_lhe_events.f	2014-11-06 06:25:51 +0000
@@ -32,7 +32,7 @@
      #                         muF2_id_str,QES_id_str
 c
       write(ifile,'(a)')
-     #     '<LesHouchesEvents version="1.0">'
+     #     '<LesHouchesEvents version="3.0">'
       write(ifile,'(a)')
      #     '  <!--'
       write(ifile,'(a)')'  <scalesfunctionalform>'
@@ -78,7 +78,7 @@
 c     statement.
       event_id=0
 c
-      write(ifile,'(a)') '<LesHouchesEvents version="1.0">'
+      write(ifile,'(a)') '<LesHouchesEvents version="3.0">'
       write(ifile,'(a)') '  <header>'
       write(ifile,'(a)') '  <MG5ProcCard>'
       open (unit=92,file=path(1:index(path," ")-1)//'proc_card_mg5.dat'

=== modified file 'Template/NLO/SubProcesses/iproc_map.f'
--- Template/NLO/SubProcesses/iproc_map.f	2014-06-26 08:45:41 +0000
+++ Template/NLO/SubProcesses/iproc_map.f	2014-11-06 06:25:51 +0000
@@ -206,7 +206,7 @@
       include 'nFKSconfigs.inc'
       include 'nexternal.inc'
       include 'genps.inc'
-      include "leshouche_info.inc"
+      include "leshouche_decl.inc"
       include "reweight_appl.inc"
       include "appl_common.inc"
 *

=== modified file 'Template/NLO/SubProcesses/leshouche_inc_chooser.f'
--- Template/NLO/SubProcesses/leshouche_inc_chooser.f	2012-08-28 21:06:34 +0000
+++ Template/NLO/SubProcesses/leshouche_inc_chooser.f	2014-11-06 06:25:51 +0000
@@ -4,7 +4,6 @@
       implicit none
       include 'nexternal.inc'
       include 'genps.inc'
-      include 'leshouche_info.inc'
       integer i,j,k
       INTEGER NFKSPROCESS
       COMMON/C_NFKSPROCESS/NFKSPROCESS
@@ -13,6 +12,11 @@
       integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc),
      &     icolup(2,nexternal,maxflow)
       common /c_leshouche_inc/idup,mothup,icolup
+      logical firsttime
+      data firsttime /.true./
+      include 'leshouche_decl.inc'
+      save idup_d, mothup_d, icolup_d
+      
 c
       if (maxproc_used.gt.maxproc) then
          write (*,*) 'ERROR in leshouche_inc_chooser: increase maxproc',
@@ -24,6 +28,12 @@
      &        maxflow,maxflow_used
          stop
       endif
+
+      if (firsttime) then
+        call read_leshouche_info(idup_d,mothup_d,icolup_d)
+        firsttime = .false.
+      endif
+
       do j=1,maxproc_used
          do i=1,nexternal
             IDUP(i,j)=IDUP_D(nFKSprocess,i,j)
@@ -42,3 +52,46 @@
       return
       end
 
+
+      subroutine read_leshouche_info(idup_d,mothup_d,icolup_d)
+C read the various information from the configs_and_props_info.dat file
+      implicit none
+      include "nexternal.inc"
+      integer itmp_array(nexternal)
+      integer i,j,k,l
+      character *200 buff
+      include 'leshouche_decl.inc'
+
+      open(unit=78, file='leshouche_info.dat', status='old')
+      do while (.true.)
+        read(78,'(a)',end=999) buff
+        if (buff(:1).eq.'#') cycle
+        if (buff(:1).eq.'I') then
+        ! idup
+        ! I  i   j   id1 ..idn -> IDUP_D(i,k,j)=idk
+          read(buff(2:),*) i,j,(itmp_array(k),k=1,nexternal)
+          do k=1,nexternal
+            idup_d(i,k,j)=itmp_array(k)
+          enddo
+        else if (buff(:1).eq.'M') then
+        ! idup
+        ! I  i   j   l   id1 ..idn -> MOTHUP_D(i,j,k,l)=idk
+          read(buff(2:),*) i,j,l,(itmp_array(k),k=1,nexternal)
+          do k=1,nexternal
+            mothup_d(i,j,k,l)=itmp_array(k)
+          enddo
+        else if (buff(:1).eq.'C') then
+        ! idup
+        ! I  i   j   l   id1 ..idn -> ICOLUP_D(i,j,k,l)=idk
+          read(buff(2:),*) i,j,l,(itmp_array(k),k=1,nexternal)
+          do k=1,nexternal
+            icolup_d(i,j,k,l)=itmp_array(k)
+          enddo
+        endif
+      enddo
+ 999  continue
+      close(78)
+
+      return 
+      end
+

=== modified file 'Template/NLO/SubProcesses/makefile_fks_dir'
--- Template/NLO/SubProcesses/makefile_fks_dir	2014-09-11 05:33:30 +0000
+++ Template/NLO/SubProcesses/makefile_fks_dir	2014-11-06 06:25:51 +0000
@@ -94,35 +94,44 @@
 gensym: $(SYM)
 	$(FC) $(LDFLAGS) -o gensym $(SYM) $(APPLLIBS) $(LINKLIBS) $(FJLIBS)
 	rm fks_singular.o
+	strip gensym
 
 test_Sij: symmetry_fks_test_Sij.o $(TEST)
 	$(FC) $(LDFLAGS) -o test_Sij symmetry_fks_test_Sij.o $(TEST) $(APPLLIBS) $(LINKLIBS) $(FJLIBS)
+	strip test_Sij
 
 test_ME: symmetry_fks_test_ME.o $(TEST)
 	$(FC) $(LDFLAGS) -o test_ME symmetry_fks_test_ME.o $(TEST) $(APPLLIBS) $(LINKLIBS) $(FJLIBS)
 	rm symmetry_fks_test_ME.o
+	strip test_ME
 
 test_MC: symmetry_fks_test_MC.o $(TEST) 
 	$(FC) $(LDFLAGS) -o test_MC symmetry_fks_test_MC.o $(TEST) $(APPLLIBS) $(LINKLIBS) $(FJLIBS)
 	rm symmetry_fks_test_MC.o
+	strip test_MC
 
 check_poles: $(POLES) $(libmadloop)
 	$(FC) $(LDFLAGS) -o check_poles $(POLES) $(NLOLIBS) $(APPLLIBS) $(LINKLIBS) $(FJLIBS)
+	strip check_poles
 
 madevent_vegas: $(VEGAS) $(libmadloop) makefile $(LIBS)
 	$(FC) $(LDFLAGS) -o madevent_vegas $(VEGAS) $(NLOLIBS) $(APPLLIBS) $(LINKLIBS) $(FJLIBS) $(FO_EXTRAPATHS) $(FO_EXTRALIBS)
+	strip madevent_vegas
 
 madevent_mintMC: $(MINTMC) $(libmadloop) makefile $(LIBS)
 	$(FC) $(LDFLAGS) -o madevent_mintMC $(MINTMC) $(NLOLIBS) $(APPLLIBS) $(LINKLIBS) $(FJLIBS) $(FO_EXTRAPATHS) $(FO_EXTRALIBS)
+	strip madevent_mintMC
 
 madevent_mintFO: $(MINTFO) $(libmadloop) makefile $(LIBS)
 	$(FC) $(LDFLAGS) -o madevent_mintFO $(MINTFO) $(NLOLIBS) $(APPLLIBS) $(LINKLIBS) $(FJLIBS) $(FO_EXTRAPATHS) $(FO_EXTRALIBS)
+	strip madevent_mintFO
 
 madevent_reweight: $(BORN_REWEIGHT) $(libmadloop) makefile $(LIBS)
 	$(FC) $(LDFLAGS) -o madevent_reweight $(BORN_REWEIGHT) $(APPLLIBS) $(LINKLIBS) $(FJLIBS) $(FO_EXTRAPATHS) $(FO_EXTRALIBS)
 
 reweight_xsec_events: $(RWGFILES) makefile $(LIBS)
 	$(FC) $(LDFLAGS) -o reweight_xsec_events $(RWGFILES) $(LINKLIBS) $(FJLIBS)
+	strip reweight_xsec_events
 
 libMadLoop.a:
 	for dir in `ls -d V*`; do cd $$dir; make; cd ../; done

=== modified file 'Template/NLO/SubProcesses/montecarlocounter.f'
--- Template/NLO/SubProcesses/montecarlocounter.f	2014-08-05 13:04:11 +0000
+++ Template/NLO/SubProcesses/montecarlocounter.f	2014-11-06 06:25:51 +0000
@@ -2517,6 +2517,7 @@
      &                       xm12,dum(1),dum(2),dum(3),dum(4),dum(5),qMC,.true.)
 
       emsca=2d0*sqrt(ebeam(1)*ebeam(2))
+      scalemax=sqrt((1-xi_i_fks)*shat)
       if(dampMCsubt)then
          call assign_scaleminmax(shat,xi_i_fks,scalemin,scalemax,ileg,xm12)
          emscasharp=(scalemax-scalemin).lt.(1d-3*scalemax)

=== modified file 'Template/NLO/SubProcesses/symmetry_fks_test_MC.f'
--- Template/NLO/SubProcesses/symmetry_fks_test_MC.f	2014-06-26 08:45:41 +0000
+++ Template/NLO/SubProcesses/symmetry_fks_test_MC.f	2014-11-06 06:25:51 +0000
@@ -144,6 +144,13 @@
       
       character*10 MonteCarlo
       common/cMonteCarloType/MonteCarlo
+
+      double precision shower_S_scale(fks_configs*2)
+     &     ,shower_H_scale(fks_configs*2),ref_H_scale(fks_configs*2)
+     &     ,pt_hardness
+      common /cshowerscale2/shower_S_scale,shower_H_scale,ref_H_scale
+     &     ,pt_hardness
+
 c      integer icomp
 c
 c     DATA
@@ -371,6 +378,10 @@
 
          call set_cms_stuff(0)
          calculatedBorn=.false.
+
+c Initialise shower_S_scale to a large value, not to get spurious dead zones
+         shower_S_scale=1d20
+
          if(ilim.eq.0)then
            call xmcsubt_wrap(p1_cnt(0,1,0),zero,y_ij_fks_ev,fxl)
          else

=== modified file 'Template/NLO/SubProcesses/write_event.f'
--- Template/NLO/SubProcesses/write_event.f	2014-06-17 08:37:03 +0000
+++ Template/NLO/SubProcesses/write_event.f	2014-11-06 06:25:51 +0000
@@ -176,7 +176,7 @@
       endif
       NPRUP=1
 
-      write(lunlhe,'(a)')'<LesHouchesEvents version="1.0">'
+      write(lunlhe,'(a)')'<LesHouchesEvents version="3.0">'
       write(lunlhe,'(a)')'  <!--'
       write(lunlhe,'(a)')'  <scalesfunctionalform>'
       write(lunlhe,'(2a)')'    muR  ',muR_id_str(1:len_trim(muR_id_str))

=== added symlink 'Template/NLO/Utilities/appl_interface_dummy.f'
=== target is u'../SubProcesses/appl_interface_dummy.f'
=== modified file 'Template/NLO/Utilities/makefile'
--- Template/NLO/Utilities/makefile	2014-03-07 11:01:10 +0000
+++ Template/NLO/Utilities/makefile	2014-11-06 06:25:51 +0000
@@ -5,7 +5,7 @@
 FF=gfortran -g -fno-automatic -ffixed-line-length-132 $(INCLUDE)
 FILES=handling_lhe_events.o fill_MC_mshell.o
 FILES2=madfks_plot.o dbook.o setcuts.o analysis_td_template.o \
-       open_output_files.o open_output_files_dummy.o
+       open_output_files.o open_output_files_dummy.o appl_interface_dummy.o
 
 %.o : %.f
 	$(FF) -c $<
@@ -13,7 +13,7 @@
 all: check_events combine_plots compare_events merge_events plot_events \
 	plot_events_lhef select_events split_events sum_plots uncertainties
 
-check_events : $(FILES) check_events.o dbook.o
+check_events : $(FILES) check_events.o dbook.o appl_interface_dummy.o
 	$(FF) $^ -o $@
 
 combine_plots : combine_plots.o

=== modified file 'Template/loop_material/StandAlone/SubProcesses/makefile.inc'
--- Template/loop_material/StandAlone/SubProcesses/makefile.inc	2014-07-24 07:10:22 +0000
+++ Template/loop_material/StandAlone/SubProcesses/makefile.inc	2014-11-06 06:25:51 +0000
@@ -4,11 +4,6 @@
     ROOT = ../..
 endif
 
-
-if [ -a ../../Source/make_opts]; then ROOT = ../..; fi;
-
-if [ -a ../Source/make_opts]; then ROOT = ../; fi;
-
 include $(ROOT)/Source/make_opts
 LIBDIR = $(ROOT)/lib/
 

=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2014-09-25 14:55:20 +0000
+++ UpdateNotes.txt	2014-11-06 06:25:51 +0000
@@ -1,5 +1,24 @@
 Update notes for MadGraph5_aMC@NLO (in reverse time order)
 
+2.2.2(XX/XX/XX) OM: Correct a bug in the integration grid (introduces in 2.1.2). This was biasing the cross-section of 
+                    processes like a a > mu+ mu- in the Effective Photon Approximation by three order of magnitude.
+                    For LHC processes no sizeable effect have been observe so far.
+                MZ: some informations for aMC@NLO runs which were before passed via include files are
+                    now read at runtime. The size of executables as well as compilation time / memory usage
+                    is reduced for complicated processes
+                RF: Fix crash #1377187 (check that cuts were consistent with the grouping was too restrictive) 
+		RF: For NLO running: added 'strip' to the makefiles to reduce executable sizes (removes symbol info)
+		Stefano Carrazza (by RF): fix for the photon PDF for the internal NNPDF sets
+		RF: Improved the check on the consistency of the cuts and the grouping of subprocesse (LO running)
+                PT: enabled PYTHIA8.2
+                OM: restore the usage of external gzip library for file larger than 4Gb which were crashing with
+                    the python gzip library
+                OM: Fixing the default card for Delphes
+                OM: Improve support of lsf cluster (thanks Josh) 
+                OM: Adding support for the UFO file functions.py (which was ignored before)
+                OM: Reduce the amount of RAM used by MadSpin in gridpack mode.
+                OM: discard in MadWidth partial width lower than \Lambda_QCD for colored particle.
+
 2.2.1(25/09/14) OM: Fix a bug preventing the generation of events at LO due to a wrong treatment of 
                       the color-flow.
 

=== modified file 'aloha/aloha_lib.py'
--- aloha/aloha_lib.py	2014-03-01 05:04:13 +0000
+++ aloha/aloha_lib.py	2014-11-06 06:25:51 +0000
@@ -209,7 +209,11 @@
                 continue
             tag = tuple(term.sort())
             if tag in items:
+                orig_prefac = items[tag].prefactor # to assume to zero 0.33333 -0.3333
                 items[tag].prefactor += term.prefactor
+                if items[tag].prefactor and \
+                    abs(items[tag].prefactor) / (abs(orig_prefac)+abs(term.prefactor)) < 1e-8:
+                    items[tag].prefactor = 0
                 del self[pos]
                 pos -=1
             else:

=== modified file 'input/.mg5_configuration_default.txt'
--- input/.mg5_configuration_default.txt	2014-09-10 11:15:57 +0000
+++ input/.mg5_configuration_default.txt	2014-11-06 06:25:51 +0000
@@ -54,12 +54,10 @@
 # timeout = 60
 
 #! Pythia8 path.
-#!  Defines the path to the pythia8 directory.
-#!  This is the directory that contains the subdirectory examples;
-#!  in particular, it is crucial that the file
-#   pythia8_path/examples/config.sh be present.
+#!  Defines the path to the main pythia8 directory (i.e. that containing
+#!  the pythia8 configure script) .
 #!  If using a relative path, that starts from the mg5 directory
-# pythia8_path = ./pythia8
+# pythia8_path =
 
 #! Herwig++ paths
 #!  specify here the paths also to HepMC ant ThePEG

=== modified file 'madgraph/__init__.py'
--- madgraph/__init__.py	2013-11-29 07:28:53 +0000
+++ madgraph/__init__.py	2014-11-06 06:25:51 +0000
@@ -35,10 +35,5 @@
     Please consider changing the path location of the code\033[0m''' % MG5DIR)
    time.sleep(1)
 MG4DIR = MG5DIR
-ReadWrite = True
+ReadWrite = os.access(MG5DIR, os.W_OK) # W_OK is for writing
 
-try:
-    open(os.path.join(MG5DIR,'.test'),'w').write('test')
-    os.remove(os.path.join(MG5DIR,'.test'))
-except IOError:
-    ReadWrite = False

=== modified file 'madgraph/interface/amcatnlo_run_interface.py'
--- madgraph/interface/amcatnlo_run_interface.py	2014-09-24 19:41:18 +0000
+++ madgraph/interface/amcatnlo_run_interface.py	2014-11-06 06:25:51 +0000
@@ -2322,7 +2322,7 @@
     def run_mcatnlo(self, evt_file):
         """runs mcatnlo on the generated event file, to produce showered-events
         """
-        logger.info('Prepairing MCatNLO run')
+        logger.info('Preparing MCatNLO run')
         try:     
             misc.gunzip(evt_file)
         except Exception:
@@ -2401,6 +2401,9 @@
         if shower == 'HERWIGPP':
             extrapaths.append(pjoin(self.options['hepmc_path'], 'lib'))
 
+        if shower == 'PYTHIA8' and not os.path.exists(pjoin(self.options['pythia8_path'], 'xmldoc')):
+            extrapaths.append(pjoin(self.options['pythia8_path'], 'lib'))
+
         if 'LD_LIBRARY_PATH' in os.environ.keys():
             ldlibrarypath = os.environ['LD_LIBRARY_PATH']
         else:
@@ -2464,8 +2467,11 @@
         # special treatment for pythia8
             files.mv(pjoin(self.me_dir, 'MCatNLO', 'Pythia8.cmd'), rundir)
             files.mv(pjoin(self.me_dir, 'MCatNLO', 'Pythia8.exe'), rundir)
-            files.ln(pjoin(self.options['pythia8_path'], 'examples', 'config.sh'), rundir)
-            files.ln(pjoin(self.options['pythia8_path'], 'xmldoc'), rundir)
+            if os.path.exists(pjoin(self.options['pythia8_path'], 'xmldoc')):
+                files.ln(pjoin(self.options['pythia8_path'], 'examples', 'config.sh'), rundir)
+                files.ln(pjoin(self.options['pythia8_path'], 'xmldoc'), rundir)
+            else:
+                files.ln(pjoin(self.options['pythia8_path'], 'share/Pythia8/xmldoc'), rundir)
         #link the hwpp exe in the rundir
         if shower == 'HERWIGPP':
             try:
@@ -2550,11 +2556,12 @@
              ' More information in %s' % pjoin(os.getcwd(), 'amcatnlo_run.log'))
 
             # run the plot creation in a secure way
-            try:
-                self.do_plot('%s -f' % self.run_name)
-            except Exception, error:
-                logger.info("Fail to make the plot. Continue...")
-                pass
+            if hep_format == 'StdHEP':
+                try:
+                    self.do_plot('%s -f' % self.run_name)
+                except Exception, error:
+                    logger.info("Fail to make the plot. Continue...")
+                    pass
 
         elif out_id == 'TOP':
             #copy the topdrawer file(s) back in events
@@ -3244,26 +3251,30 @@
             if os.path.exists(pdfinput):
                 input_files.append(pdfinput)
             input_files.append(pjoin(os.path.dirname(exe), os.path.pardir, 'reweight_xsec_events'))
+            input_files.append(pjoin(cwd, os.path.pardir, 'leshouche_info.dat'))
             input_files.append(args[0])
             output_files.append('%s.rwgt' % os.path.basename(args[0]))
             output_files.append('reweight_xsec_events.output')
             output_files.append('scale_pdf_dependence.dat')
 
             return self.cluster.submit2(exe, args, cwd=cwd, 
-                             input_files=input_files, output_files=output_files) 
+                             input_files=input_files, output_files=output_files,
+                             required_output=output_files) 
 
         elif 'ajob' in exe:
             # the 'standard' amcatnlo job
             # check if args is a list of string 
             if type(args[0]) == str:
-                input_files, output_files, args = self.getIO_ajob(exe,cwd, args)
+                input_files, output_files, required_output, args = self.getIO_ajob(exe,cwd, args)
                 #submitting
                 self.cluster.submit2(exe, args, cwd=cwd, 
-                             input_files=input_files, output_files=output_files)
+                             input_files=input_files, output_files=output_files,
+                             required_output=required_output)
 
                 # keep track of folders and arguments for splitted evt gen
-                if len(args) == 4 and '_' in output_files[-1]:
-                    self.split_folders[pjoin(cwd,output_files[-1])] = [exe] + args
+                subfolder=output_files[-1].split('/')[0]
+                if len(args) == 4 and '_' in subfolder:
+                    self.split_folders[pjoin(cwd,subfolder)] = [exe] + args
 
         elif 'shower' in exe:
             # a shower job
@@ -3276,7 +3287,10 @@
                 input_files.append(pjoin(cwd, 'Pythia8.exe'))
                 input_files.append(pjoin(cwd, 'Pythia8.cmd'))
                 input_files.append(pjoin(cwd, 'config.sh'))
-                input_files.append(pjoin(self.options['pythia8_path'], 'xmldoc'))
+                if os.path.exists(pjoin(self.options['pythia8_path'], 'xmldoc')):
+                    input_files.append(pjoin(self.options['pythia8_path'], 'xmldoc'))
+                else:
+                    input_files.append(pjoin(self.options['pythia8_path'], 'share/Pythia8/xmldoc'))
             else:
                 input_files.append(pjoin(cwd, 'MCATNLO_%s_EXE' % shower))
                 input_files.append(pjoin(cwd, 'MCATNLO_%s_input' % shower))
@@ -3328,11 +3342,14 @@
         
         keep_fourth_arg = False
         output_files = []
+        required_output = []
         input_files = [pjoin(self.me_dir, 'MGMEVersion.txt'),
                      pjoin(self.me_dir, 'SubProcesses', 'randinit'),
                      pjoin(cwd, 'symfact.dat'),
                      pjoin(cwd, 'iproc.dat'),
                      pjoin(cwd, 'initial_states_map.dat'),
+                     pjoin(cwd, 'configs_and_props_info.dat'),
+                     pjoin(cwd, 'leshouche_info.dat'),
                      pjoin(cwd, 'param_card.dat'),
                      pjoin(cwd, 'FKS_params.dat')]
 
@@ -3370,7 +3387,13 @@
                 if os.path.exists(pjoin(cwd,current)):
                     input_files.append(pjoin(cwd, current))
                 output_files.append(current)
+
+                required_output.append('%s/results.dat' % current)
+                required_output.append('%s/log.txt' % current)
+                required_output.append('%s/mint_grids' % current)
+                required_output.append('%s/grid.MC_integer' % current)
                 if len(args) == 4:
+                    required_output.append('%s/scale_pdf_dependence.dat' % current)
                     args[2] = '-1'
                     # use a grid train on another part
                     base = '%s_G%s' % (args[3],i)
@@ -3414,6 +3437,11 @@
                     keep_fourth_arg = True
                     # this is for the split event generation
                     output_files.append('G%s%s_%s' % (args[1], i, args[3]))
+                required_output.append('%s/log_MINT%s.txt' % (current,args[2]))
+                if args[2] in ['0','1']:
+                    required_output.append('%s/results.dat' % current)
+                if args[2] == '1':
+                    output_files.append('%s/results.dat' % current)
 
         else:
             raise aMCatNLOError, 'not valid arguments: %s' %(', '.join(args))
@@ -3426,7 +3454,7 @@
         if len(args) == 4 and not keep_fourth_arg:
             args = args[:3]
             
-        return input_files, output_files, args
+        return input_files, output_files, required_output,  args
             
     def write_madinMMC_file(self, path, run_mode, mint_mode):
         """writes the madinMMC_?.2 file"""

=== modified file 'madgraph/interface/common_run_interface.py'
--- madgraph/interface/common_run_interface.py	2014-09-24 19:41:18 +0000
+++ madgraph/interface/common_run_interface.py	2014-11-06 06:25:51 +0000
@@ -195,6 +195,7 @@
             raise self.ConfigurationError, '''Can\'t load MG5.
             The variable mg5_path should not be correctly configure.'''
         
+        
         ufo_path = pjoin(self.me_dir,'bin','internal', 'ufomodel')
         # Import model
         if not MADEVENT:
@@ -220,7 +221,7 @@
 #            raise self.InvalidCmd, 'The UFO model does not include partial widths information. Impossible to compute widths automatically'
             
         # check if the name are passed to default MG5
-        if '-modelname' in open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')).read():
+        if '-modelname' not in open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')).read():
             model.pass_particles_name_in_mg_default()        
         model = model_reader.ModelReader(model)
         particles_name = dict([(p.get('name'), p.get('pdg_code'))
@@ -949,7 +950,7 @@
             self.cluster.launch_and_wait(self.dirbin+'/run_hep2lhe',
                                          argument= [pydir],
                                         cwd=pjoin(self.me_dir,'Events'),
-                                        stdout=open(os.devnull,'w'))
+                                        stdout=os.devnull)
 
             logger.info('Warning! Never use this lhe file for detector studies!')
             # Creating ROOT file
@@ -2058,7 +2059,11 @@
     def get_lhapdf_pdfsetsdir(self):
         lhapdf_version = self.get_lhapdf_version()
 
-        if lhapdf_version.startswith('5.'):
+        # check if the LHAPDF_DATA_PATH variable is defined
+        if 'LHAPDF_DATA_PATH' in os.environ.keys() and os.environ['LHAPDF_DATA_PATH']:
+            datadir = os.environ['LHAPDF_DATA_PATH']
+
+        elif lhapdf_version.startswith('5.'):
             datadir = subprocess.Popen([self.options['lhapdf'], '--pdfsets-path'],
                          stdout = subprocess.PIPE).stdout.read().strip()
 
@@ -2141,7 +2146,12 @@
                     elif len(split) == 2:
                         if re.search(r'''\[[A-Z]\]eV\^''', split[1]):
                             all_var = [comment.strip().lower()]
+                    elif len(split) >=2 and split[1].startswith('('):
+                        all_var = [split[0].strip().lower()]
                     else:
+                        if not bname.startswith('qnumbers'):
+                            logger.debug("not recognize information for %s %s : %s",
+                                      bname, lha_id, comment)
                         # not recognized format
                         continue
 
@@ -2363,9 +2373,13 @@
         if args[0] in self.special_shortcut:
             if len(args) == 1:
                 values = {}
-            elif len(args) == 2:  
+            elif len(args) == 2: 
+                targettype = float
+                if args[1].strip().isdigit():
+                    targettype = int
+                 
                 try:  
-                    values = {'0': float(args[1])}
+                    values = {'0': targettype(args[1])}
                 except ValueError as e:
                     logger.warning("Wrong argument: The last entry should be a number.")
                     return
@@ -2885,4 +2899,4 @@
             except:
                 import internal.madweight.Cards as mwcards
             self.mw_card = mwcards.Card(pjoin(self.me_dir,'Cards','MadWeight_card.dat'))
- 
\ No newline at end of file
+ 

=== modified file 'madgraph/interface/madevent_interface.py'
--- madgraph/interface/madevent_interface.py	2014-09-24 19:41:18 +0000
+++ madgraph/interface/madevent_interface.py	2014-11-06 06:25:51 +0000
@@ -4064,11 +4064,10 @@
             
 
             logfile = pjoin(Pdir, 'gen_ximprove.log')
-            proc = misc.Popen([pjoin(bindir, 'gen_ximprove')],
+            misc.call([pjoin(bindir, 'gen_ximprove')],
                                     stdin=subprocess.PIPE,
                                     stdout=open(logfile,'w'),
                                     cwd=Pdir)
-            proc.communicate('%s 1 F\n' % (precision))
 
             if os.path.exists(pjoin(Pdir, 'ajob1')):
                 alljobs = glob.glob(pjoin(Pdir,'ajob*'))

=== modified file 'madgraph/interface/madgraph_interface.py'
--- madgraph/interface/madgraph_interface.py	2014-09-23 02:50:45 +0000
+++ madgraph/interface/madgraph_interface.py	2014-11-06 06:25:51 +0000
@@ -3541,7 +3541,7 @@
 
         # Now check for perturbation orders, specified in between squared brackets
         perturbation_couplings_pattern = \
-          re.compile("^(?P<proc>.+)\s*\[\s*((?P<option>\w+)\s*\=)?\s*"+\
+          re.compile("^(?P<proc>.+>.+)\s*\[\s*((?P<option>\w+)\s*\=)?\s*"+\
                                "(?P<pertOrders>(\w+\s*)*)\s*\]\s*(?P<rest>.*)$")
         perturbation_couplings_re = perturbation_couplings_pattern.match(line)
         perturbation_couplings = ""
@@ -3595,7 +3595,7 @@
             # Notice that one can have a negative value of the squared order to
             # indicate that one should take the N^{n}LO contribution into account.
             order_pattern = re.compile(\
-           "^(?P<before>.+)\s+(?P<name>(\w|(\^2))+)\s*(?P<type>"+\
+           "^(?P<before>.+>.+)\s+(?P<name>(\w|(\^2))+)\s*(?P<type>"+\
                     "(=|(<=)|(==)|(===)|(!=)|(>=)|<|>))\s*(?P<value>-?\d+)\s*$")
             order_re = order_pattern.match(line)
             while order_re:
@@ -4595,7 +4595,7 @@
             out = open(pjoin(MG5DIR, 'Template','Common', 'Cards', 'delphes_card_default.dat'), 'w')
             out.write(data)
         if args[0] == 'Delphes3':
-            files.cp(pjoin(MG5DIR, 'Delphes','examples','delphes_card_CMS_PileUp.tcl'),
+            files.cp(pjoin(MG5DIR, 'Delphes','examples','delphes_card_CMS.tcl'),
                      pjoin(MG5DIR,'Template', 'Common', 'Cards', 'delphes_card_default.dat'))
             files.cp(pjoin(MG5DIR, 'Delphes','examples','delphes_card_CMS.tcl'),
                      pjoin(MG5DIR,'Template', 'Common', 'Cards', 'delphes_card_CMS.dat'))
@@ -6288,6 +6288,13 @@
 
         for pid in particles:
             width = param['decay'].get((pid,)).value
+            particle = self._curr_model.get_particle(pid) 
+            #if particle['color'] !=1 and 0 < width.real < 0.1:
+            #    logger.warning("width of colored particle \"%s(%s)\" lower than QCD scale: %s. Set width to zero "
+            #                   % (particle.get('name'), pid, width.real))
+            #    width = 0
+                
+            
             if not pid in param['decay'].decay_table:
                 continue
             if pid not in decay_info:
@@ -6295,6 +6302,12 @@
             for BR in param['decay'].decay_table[pid]:
                 if len(BR.lhacode) == 3 and skip_2body:
                     continue
+                if BR.value * width <0.1 and particle['color'] !=1:
+                    logger.warning("partial width of particle %s lower than QCD scale:%s. Set it to zero. (%s)" \
+                                   % (particle.get('name'), BR.value * width, BR.lhacode[1:]))
+                                     
+                    continue
+                
                 decay_info[pid].append([BR.lhacode[1:], BR.value * width])
 
         madevent_interface.MadEventCmd.update_width_in_param_card(decay_info,

=== modified file 'madgraph/iolibs/export_fks.py'
--- madgraph/iolibs/export_fks.py	2014-09-17 07:43:49 +0000
+++ madgraph/iolibs/export_fks.py	2014-11-06 06:25:51 +0000
@@ -469,15 +469,25 @@
                                  matrix_element, 
                                  fortran_model)
 
-        filename = 'leshouche_info.inc'
-        self.write_leshouche_info_file(writers.FortranWriter(filename), 
-                                 matrix_element,
-                                 fortran_model)
-
-        filename = 'configs_and_props_info.inc'
-        nconfigs=self.write_configs_and_props_info_file(
-                              writers.FortranWriter(filename), 
-                              matrix_element,
+        filename = 'leshouche_info.dat'
+        nfksconfs,maxproc,maxflow,nexternal=\
+                self.write_leshouche_info_file(filename,matrix_element)
+
+        filename = 'leshouche_decl.inc'
+        self.write_leshouche_info_declarations(
+                              writers.FortranWriter(filename), 
+                              nfksconfs,maxproc,maxflow,nexternal,
+                              fortran_model)
+
+        filename = 'configs_and_props_info.dat'
+        nconfigs,max_leg_number,nfksconfs=self.write_configs_and_props_info_file(
+                              filename, 
+                              matrix_element)
+
+        filename = 'configs_and_props_decl.inc'
+        self.write_configs_and_props_info_declarations(
+                              writers.FortranWriter(filename), 
+                              nconfigs,max_leg_number,nfksconfs,
                               fortran_model)
         
         filename = 'real_from_born_configs.inc'
@@ -599,6 +609,7 @@
                      'cluster.inc',
                      'cluster.f',
                      'reweight.f',
+                     'randinit',
                      'sudakov.inc',
                      'maxconfigs.inc',
                      'timing_variables.inc']
@@ -619,7 +630,7 @@
 
         #import nexternal/leshouches in Source
         ln('nexternal.inc', '../../Source', log=False)
-        ln('leshouche_info.inc', '../../Source', log=False)
+        ln('leshouche_decl.inc', '../../Source', log=False)
 
 
         # Return to SubProcesses dir
@@ -653,6 +664,10 @@
 #            check_param_card.convert_to_mg5card(param_card, mg5_param)
 #            check_param_card.check_valid_param_card(mg5_param)
 
+#        # write the model functions get_mass/width_from_id
+        filename = os.path.join(self.dir_path,'Source','MODEL','get_mass_width_fcts.f')
+        makeinc = os.path.join(self.dir_path,'Source','MODEL','makeinc.inc')
+        self.write_get_mass_width_file(writers.FortranWriter(filename), makeinc, self.model)
 
 #        # Write maxconfigs.inc based on max of ME's/subprocess groups
         filename = os.path.join(self.dir_path,'Source','maxconfigs.inc')
@@ -826,25 +841,97 @@
         writer.writelines(lines2+lines)
 
 
-
-    def write_configs_and_props_info_file(self, writer, matrix_element, fortran_model):
+    #===============================================================================
+    # write_get_mass_width_file
+    #===============================================================================
+    #test written
+    def write_get_mass_width_file(self, writer, makeinc, model):
+        """Write the get_mass_width_file.f file for MG4.
+        Also update the makeinc.inc file
+        """
+        mass_particles = [p for p in model['particles'] if p['mass'].lower() != 'zero'] 
+        width_particles = [p for p in model['particles'] if p['width'].lower() != 'zero'] 
+        
+        iflines_mass = ''
+        iflines_width = ''
+
+        for i, part in enumerate(mass_particles):
+            if i == 0:
+                ifstring = 'if'
+            else:
+                ifstring = 'else if'
+            if part['self_antipart']:
+                iflines_mass += '%s (id.eq.%d) then\n' % \
+                        (ifstring, part.get_pdg_code())
+            else:
+                iflines_mass += '%s (id.eq.%d.or.id.eq.%d) then\n' % \
+                        (ifstring, part.get_pdg_code(), part.get_anti_pdg_code())
+            iflines_mass += 'get_mass_from_id=abs(%s)\n' % part.get('mass')
+
+        for i, part in enumerate(width_particles):
+            if i == 0:
+                ifstring = 'if'
+            else:
+                ifstring = 'else if'
+            if part['self_antipart']:
+                iflines_width += '%s (id.eq.%d) then\n' % \
+                        (ifstring, part.get_pdg_code())
+            else:
+                iflines_width += '%s (id.eq.%d.or.id.eq.%d) then\n' % \
+                        (ifstring, part.get_pdg_code(), part.get_anti_pdg_code())
+            iflines_width += 'get_width_from_id=abs(%s)\n' % part.get('width')
+
+        replace_dict = {'iflines_mass' : iflines_mass,
+                        'iflines_width' : iflines_width}
+
+        file = open(os.path.join(_file_path, \
+                          'iolibs/template_files/get_mass_width_fcts.inc')).read()
+        file = file % replace_dict
+        
+        # Write the file
+        writer.writelines(file)
+
+        # update the makeinc
+        makeinc_content = open(makeinc).read()
+        makeinc_content = makeinc_content.replace('MODEL = ', 'MODEL = get_mass_width_fcts.o ')
+        open(makeinc, 'w').write(makeinc_content)
+
+        return 
+
+
+    def write_configs_and_props_info_declarations(self, writer, max_iconfig, max_leg_number, nfksconfs, fortran_model):
+        """writes the declarations for the variables relevant for configs_and_props
+        """
+        lines = []
+        lines.append("integer ifr,lmaxconfigs_used,max_branch_used")
+        lines.append("parameter (lmaxconfigs_used=%4d)" % max_iconfig)
+        lines.append("parameter (max_branch_used =%4d)" % -max_leg_number)
+        lines.append("integer mapconfig_d(%3d,0:lmaxconfigs_used)" % nfksconfs)
+        lines.append("integer iforest_d(%3d,2,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+        lines.append("integer sprop_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+        lines.append("integer tprid_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+        lines.append("double precision pmass_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+        lines.append("double precision pwidth_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+        lines.append("integer pow_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nfksconfs)
+
+        writer.writelines(lines)
+
+
+    def write_configs_and_props_info_file(self, filename, matrix_element):
         """writes the configs_and_props_info.inc file that cointains
         all the (real-emission) configurations (IFOREST) as well as
         the masses and widths of intermediate particles"""
         lines = []
+        lines.append("# C -> MAPCONFIG_D") 
+        lines.append("# F/D -> IFOREST_D") 
+        lines.append("# S -> SPROP_D") 
+        lines.append("# T -> TPRID_D") 
+        lines.append("# M -> PMASS_D/PWIDTH_D") 
+        lines.append("# P -> POW_D") 
         lines2 = []
         nconfs = len(matrix_element.get_fks_info_list())
         (nexternal, ninitial) = matrix_element.real_processes[0].get_nexternal_ninitial()
 
-        lines.append("integer ifr,lmaxconfigs_used,max_branch_used")
-        lines.append("integer mapconfig_d(%3d,0:lmaxconfigs_used)" % nconfs)
-        lines.append("integer iforest_d(%3d,2,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-        lines.append("integer sprop_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-        lines.append("integer tprid_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-        lines.append("double precision pmass_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-        lines.append("double precision pwidth_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-        lines.append("integer pow_d(%3d,-max_branch_used:-1,lmaxconfigs_used)" % nconfs)
-
         max_iconfig=0
         max_leg_number=0
 
@@ -872,7 +959,7 @@
                 lines.append("# Diagram %d for nFKSprocess %d" % \
                                  (helas_diag.get('number'),iFKS))
                 # Correspondance between the config and the amplitudes
-                lines.append("data mapconfig_d(%3d,%4d)/%4d/" % (iFKS,iconfig,
+                lines.append("C   %4d   %4d   %4d " % (iFKS,iconfig,
                                                            helas_diag.get('number')))
     
                 # Need to reorganize the topology so that we start with all
@@ -891,15 +978,16 @@
                 for vert in allchannels:
                     daughters = [leg.get('number') for leg in vert.get('legs')[:-1]]
                     last_leg = vert.get('legs')[-1]
-                    lines.append("data (iforest_d(%3d, ifr,%3d,%4d),ifr=1,%d)/%s/" % \
-                                     (iFKS,last_leg.get('number'), iconfig, len(daughters),
-                                      ",".join(["%3d" % d for d in daughters])))
+                    lines.append("F   %4d   %4d   %4d   %4d" % \
+                                     (iFKS,last_leg.get('number'), iconfig, len(daughters)))
+                    for d in daughters:
+                        lines.append("D   %4d" % d)
                     if vert in schannels:
-                        lines.append("data sprop_d(%3d,%4d,%4d)/%8d/" % \
+                        lines.append("S   %4d   %4d   %4d   %10d" % \
                                          (iFKS,last_leg.get('number'), iconfig,
                                           last_leg.get('id')))
                     elif vert in tchannels[:-1]:
-                        lines.append("data tprid_d(%3d,%4d,%4d)/%8d/" % \
+                        lines.append("T   %4d   %4d   %4d   %10d" % \
                                          (iFKS,last_leg.get('number'), iconfig,
                                           abs(last_leg.get('id'))))
 
@@ -909,7 +997,7 @@
     
             # Write out number of configs
             lines.append("# Number of configs for nFKSprocess %d" % iFKS)
-            lines.append("data mapconfig_d(%3d,0)/%4d/" % (iFKS,iconfig))
+            lines.append("C   %4d   %4d   %4d" % (iFKS,0,iconfig))
             
             # write the props.inc information
             lines2.append("# ")
@@ -921,59 +1009,49 @@
                     leg = vertex.get('legs')[-1]
                     if leg.get('id') == 21 and 21 not in particle_dict:
                         # Fake propagator used in multiparticle vertices
-                        mass = 'zero'
-                        width = 'zero'
                         pow_part = 0
                     else:
                         particle = particle_dict[leg.get('id')]
-                    # Get mass
-                        if particle.get('mass').lower() == 'zero':
-                            mass = particle.get('mass')
-                        else:
-                            mass = "abs(%s)" % particle.get('mass')
-                    # Get width
-                        if particle.get('width').lower() == 'zero':
-                            width = particle.get('width')
-                        else:
-                            width = "abs(%s)" % particle.get('width')
     
                         pow_part = 1 + int(particle.is_boson())
     
-                    lines2.append("pmass_d (%3d,%3d,%4d) = %s " % \
-                                     (iFKS,leg.get('number'), iconf + 1, mass))
-                    lines2.append("pwidth_d(%3d,%3d,%4d) = %s " % \
-                                     (iFKS,leg.get('number'), iconf + 1, width))
-                    lines2.append("pow_d   (%3d,%3d,%4d) = %d " % \
+                    lines2.append("M   %4d   %4d   %4d   %10d " % \
+                                     (iFKS,leg.get('number'), iconf + 1, leg.get('id')))
+                    lines2.append("P   %4d   %4d   %4d   %4d " % \
                                      (iFKS,leg.get('number'), iconf + 1, pow_part))
-
-
-
     
-        lines.append("# ")
-        # insert the declaration of the sizes arrays at the beginning of the file
-        lines.insert(1,"parameter (lmaxconfigs_used=%4d)" % max_iconfig)
-        lines.insert(2,"parameter (max_branch_used =%4d)" % -max_leg_number)
-
         # Write the file
-        writer.writelines(lines+lines2)
-
-        return max_iconfig
-
-
-
-    def write_leshouche_info_file(self, writer, matrix_element, fortran_model):
-        """writes the leshouche_info.inc file which contains the LHA informations
-        for all the real emission processes"""
-        lines = []
-        nconfs = len(matrix_element.get_fks_info_list())
+        open(filename,'w').write('\n'.join(lines+lines2))
+
+        return max_iconfig, max_leg_number, nconfs
+
+
+    def write_leshouche_info_declarations(self, writer, nfksconfs, 
+                                  maxproc, maxflow, nexternal, fortran_model):
+        """writes the declarations for the variables relevant for leshouche_info
+        """
+        lines = []
+        lines.append('integer maxproc_used, maxflow_used')
+        lines.append('parameter (maxproc_used = %d)' % maxproc)
+        lines.append('parameter (maxflow_used = %d)' % maxflow)
+        lines.append('integer idup_d(%d,%d,maxproc_used)' % (nfksconfs, nexternal))
+        lines.append('integer mothup_d(%d,%d,%d,maxproc_used)' % (nfksconfs, 2, nexternal))
+        lines.append('integer icolup_d(%d,%d,%d,maxflow_used)' % (nfksconfs, 2, nexternal))
+
+        writer.writelines(lines)
+
+
+    def write_leshouche_info_file(self, filename, matrix_element):
+        """writes the leshouche_info.inc file which contains 
+        the LHA informations for all the real emission processes
+        """
+        lines = []
+        lines.append("# I -> IDUP_D")
+        lines.append("# M -> MOTHUP_D")
+        lines.append("# C -> ICOLUP_D")
+        nfksconfs = len(matrix_element.get_fks_info_list())
         (nexternal, ninitial) = matrix_element.real_processes[0].get_nexternal_ninitial()
 
-        lines.append('integer idup_d(%d,%d,maxproc_used)' % (nconfs, nexternal))
-        lines.append('integer mothup_d(%d,%d,%d,maxproc_used)' % (nconfs, 2, nexternal))
-        lines.append('integer icolup_d(%d,%d,%d,maxflow_used)' % (nconfs, 2, nexternal))
-        lines.append('integer ilh')
-        lines.append('')
-
         maxproc = 0
         maxflow = 0
         for i, conf in enumerate(matrix_element.get_fks_info_list()):
@@ -984,10 +1062,10 @@
             maxproc = max(maxproc, nprocs)
             maxflow = max(maxflow, nflows)
 
-        firstlines = ['integer maxproc_used, maxflow_used',
-                      'parameter (maxproc_used = %d)' % maxproc,
-                      'parameter (maxflow_used = %d)' % maxflow ]
-        writer.writelines(firstlines + lines)
+        # Write the file
+        open(filename,'w').write('\n'.join(lines))
+
+        return nfksconfs, maxproc, maxflow, nexternal
 
 
     def write_pdf_wrapper(self, writer, matrix_element, fortran_model):
@@ -1471,7 +1549,7 @@
 
         replace_dict = {}
         replace_dict['mesq'] = 'CHaveraged'
-        replace_dict['corr'] = ' '.join(matrix_elements.get('processes')[0].\
+        replace_dict['corr'] = ' '.join(matrix_elements[0].get('processes')[0].\
                                                   get('perturbation_couplings'))
         replace_dict['irreg'] = 'CDR'
         replace_dict['aspow'] = QCD
@@ -2364,14 +2442,14 @@
         lines = []
         for iproc, proc in enumerate(matrix_element.get('processes')):
             legs = proc.get_legs_with_decays()
-            lines.append("DATA (IDUP_D(%d,ilh,%d),ilh=1,%d)/%s/" % \
-                         (ime, iproc + 1, nexternal,
-                          ",".join([str(l.get('id')) for l in legs])))
+            lines.append("I   %4d   %4d       %s" % \
+                         (ime, iproc + 1,
+                          " ".join([str(l.get('id')) for l in legs])))
             for i in [1, 2]:
-                lines.append("DATA (MOTHUP_D(%d,%d,ilh,%3r),ilh=1,%2r)/%s/" % \
-                         (ime, i, iproc + 1, nexternal,
-                          ",".join([ "%3r" % 0 ] * ninitial + \
-                                   [ "%3r" % i ] * (nexternal - ninitial))))
+                lines.append("M   %4d   %4d   %4d      %s" % \
+                         (ime, i, iproc + 1,
+                          " ".join([ "%3d" % 0 ] * ninitial + \
+                                   [ "%3d" % i ] * (nexternal - ninitial))))
     
             # Here goes the color connections corresponding to the JAMPs
             # Only one output, for the first subproc!
@@ -2379,9 +2457,9 @@
                 # If no color basis, just output trivial color flow
                 if not matrix_element.get('color_basis'):
                     for i in [1, 2]:
-                        lines.append("DATA (ICOLUP_D(%d,%d,ilh,  1),ilh=1,%2r)/%s/" % \
-                                 (ime, i, nexternal,
-                                  ",".join([ "%3r" % 0 ] * nexternal)))
+                        lines.append("C   %4d   %4d   1      %s" % \
+                                 (ime, i, 
+                                  " ".join([ "%3d" % 0 ] * nexternal)))
                     color_flow_list = []
                     nflow = 1
     
@@ -2399,15 +2477,14 @@
                     # And output them properly
                     for cf_i, color_flow_dict in enumerate(color_flow_list):
                         for i in [0, 1]:
-                            lines.append("DATA (ICOLUP_D(%d,%d,ilh,%3r),ilh=1,%2r)/%s/" % \
-                                 (ime, i + 1, cf_i + 1, nexternal,
-                                  ",".join(["%3r" % color_flow_dict[l.get('number')][i] \
+                            lines.append("C   %4d   %4d   %4d      %s" % \
+                                 (ime, i + 1, cf_i + 1,
+                                  " ".join(["%3d" % color_flow_dict[l.get('number')][i] \
                                             for l in legs])))
 
                     nflow = len(color_flow_list)
 
         nproc = len(matrix_element.get('processes'))
-        lines.append('')
     
         return lines, nproc, nflow
 

=== modified file 'madgraph/iolibs/export_v4.py'
--- madgraph/iolibs/export_v4.py	2014-09-08 11:03:40 +0000
+++ madgraph/iolibs/export_v4.py	2014-11-06 06:25:51 +0000
@@ -4698,13 +4698,12 @@
                  
                 double complex gal(2)
                 common/weak/ gal
+                
+                double precision MU_R
+                common/rscale/ MU_R
 
                 """        
-        if self.model.get('expansion_order'):
-            header=header+"""double precision MU_R
-                common/rscale/ MU_R
 
-                """
         header = header+"""double precision Nf
                 parameter(Nf=%d)
                 """ % self.model.get_nflav()
@@ -4717,13 +4716,11 @@
                      
                     %(complex_mp_format)s %(mp_prefix)sgal(2)
                     common/MP_weak/ %(mp_prefix)sgal
-    
+                    
+                    %(complex_mp_format)s %(mp_prefix)sMU_R
+                    common/MP_rscale/ %(mp_prefix)sMU_R
                     """        
-            if self.model.get('expansion_order'):
-                header=header+"""%(complex_mp_format)s %(mp_prefix)sMU_R
-                    common/MP_rscale/ %(mp_prefix)sMU_R
-    
-                    """            
+                              
             mp_fsock.writelines(header%{'real_mp_format':self.mp_real_format,
                                   'complex_mp_format':self.mp_complex_format,
                                   'mp_prefix':self.mp_prefix})
@@ -5240,6 +5237,58 @@
             input = pjoin(model_path,'Fortran','functions.f')
             file.writelines(fsock, open(input).read())
             fsock.write_comment_line(' END USER DEFINE FUNCTIONS ')
+            
+        # check for functions define in the UFO model
+        ufo_fct = self.model.get('functions')
+        if ufo_fct:
+            fsock.write_comment_line(' START UFO DEFINE FUNCTIONS ')
+            for fct in ufo_fct:
+                # already handle by default
+                if fct.name not in ["complexconjugate", "re", "im", "sec", "csc", "asec", "acsc",
+                                    "theta_function", "cond", "reglog", "arg"]:
+                    ufo_fct_template = """
+          double complex function %(name)s(%(args)s)
+          implicit none
+          double complex %(args)s
+          %(name)s = %(fct)s
+
+          return
+          end
+          """
+                    text = ufo_fct_template % {
+                                'name': fct.name,
+                                'args': ", ".join(fct.arguments),                
+                                'fct': self.p_to_f.parse(fct.expr)
+                                 }
+                    fsock.writelines(text)
+            if self.opt['mp']:
+                fsock.write_comment_line(' START UFO DEFINE FUNCTIONS FOR MP')
+                for fct in ufo_fct:
+                    # already handle by default
+                    if fct.name not in ["complexconjugate", "re", "im", "sec", "csc", "asec", "acsc",
+                                        "theta_function", "cond", "reglog", "arg"]:
+                        ufo_fct_template = """
+          %(complex_mp_format)s function mp__%(name)s(mp__%(args)s)
+          implicit none
+          %(complex_mp_format)s mp__%(args)s
+          mp__%(name)s = %(fct)s
+
+          return
+          end
+          """
+                        text = ufo_fct_template % {
+                                'name': fct.name,
+                                'args': ", mp__".join(fct.arguments),                
+                                'fct': self.mp_p_to_f.parse(fct.expr),
+                                'complex_mp_format': self.mp_complex_format
+                                 }
+                        fsock.writelines(text)
+
+
+                    
+            fsock.write_comment_line(' STOP UFO DEFINE FUNCTIONS ')                    
+
+        
 
     def create_makeinc(self):
         """create makeinc.inc containing the file to compile """

=== modified file 'madgraph/iolibs/template_files/addmothers.f'
--- madgraph/iolibs/template_files/addmothers.f	2014-09-25 14:36:15 +0000
+++ madgraph/iolibs/template_files/addmothers.f	2014-11-06 06:25:51 +0000
@@ -14,7 +14,7 @@
       integer jpart(7,-nexternal+3:2*nexternal-3),npart,ip,numproc
       double precision pb(0:4,-nexternal+3:2*nexternal-3)
       double precision rscale,aqcd,aqed,targetamp(maxflow)
-      character*300 buff
+      character*1000 buff
       character*20 cform
       logical flip ! If .true., initial state is mirrored
 

=== added file 'madgraph/iolibs/template_files/get_mass_width_fcts.inc'
--- madgraph/iolibs/template_files/get_mass_width_fcts.inc	1970-01-01 00:00:00 +0000
+++ madgraph/iolibs/template_files/get_mass_width_fcts.inc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,23 @@
+DOUBLE PRECISION FUNCTION GET_MASS_FROM_ID(ID)
+IMPLICIT NONE
+INTEGER ID
+INCLUDE 'coupl.inc'
+
+%(iflines_mass)sELSE
+GET_MASS_FROM_ID=0d0
+ENDIF
+RETURN
+END
+
+
+DOUBLE PRECISION FUNCTION GET_WIDTH_FROM_ID(ID)
+IMPLICIT NONE
+INTEGER ID
+INCLUDE 'coupl.inc'
+
+%(iflines_width)sELSE
+GET_WIDTH_FROM_ID=0d0
+ENDIF
+RETURN
+END
+

=== modified file 'madgraph/iolibs/template_files/madevent_combine_events.f'
--- madgraph/iolibs/template_files/madevent_combine_events.f	2014-08-05 16:12:48 +0000
+++ madgraph/iolibs/template_files/madevent_combine_events.f	2014-11-06 06:25:51 +0000
@@ -57,7 +57,7 @@
       character*30 param_card_name
       common/to_param_card_name/param_card_name
 
-      character*300 buff
+      character*1000 buff
       logical u_syst, has_negative
       character*(s_bufflen) s_buff(7)
       integer nclus
@@ -119,7 +119,7 @@
       I4 = 4
       R8 = 8
       record_length = 4*I4+maxexternal*I4*7+maxexternal*5*R8+3*R8+
-     &   300+7*s_bufflen+max_particles*clus_bufflen
+     &   1000+7*s_bufflen+max_particles*clus_bufflen
 C $B$ scratch_name $B$ !this is tag for automatic modification by MW
       filename='scratch'
 C $E$ scratch_name $E$ !this is tag for automatic modification by MW
@@ -625,7 +625,7 @@
       double precision sscale,aqcd,aqed,tmpsum
       integer ievent,jseed
       logical done,found
-      character*300 buff
+      character*1000 buff
       logical u_syst
       character*(s_bufflen) s_buff(7)
       character*300 fullname

=== modified file 'madgraph/iolibs/template_files/matrix_madevent_group_v4.inc'
--- madgraph/iolibs/template_files/matrix_madevent_group_v4.inc	2014-05-23 02:30:58 +0000
+++ madgraph/iolibs/template_files/matrix_madevent_group_v4.inc	2014-11-06 06:25:51 +0000
@@ -258,7 +258,7 @@
 	  DO M = 1, NAMPSO	
 		DO N = 1, NAMPSO
           IF (CHOSEN_SO_CONFIGS(SQSOINDEX%(proc_id)s(M,N))) THEN
-            Jamp2(i)=Jamp2(i)+Jamp(i,m)*dconjg(Jamp(i,n))
+            Jamp2(i)=Jamp2(i)+DABS(DBLE(Jamp(i,m)*dconjg(Jamp(i,n))))
 		  ENDIF
 		enddo
 	  enddo

=== modified file 'madgraph/iolibs/template_files/matrix_madevent_v4.inc'
--- madgraph/iolibs/template_files/matrix_madevent_v4.inc	2014-05-23 02:30:58 +0000
+++ madgraph/iolibs/template_files/matrix_madevent_v4.inc	2014-11-06 06:25:51 +0000
@@ -249,7 +249,7 @@
 	  DO M = 1, NAMPSO	
 		DO N = 1, NAMPSO
           IF (CHOSEN_SO_CONFIGS(SQSOINDEX%(proc_id)s(M,N))) THEN		
-            Jamp2(i)=Jamp2(i)+Jamp(i,m)*dconjg(Jamp(i,n))
+            Jamp2(i)=Jamp2(i)+DABS(DBLE(Jamp(i,m)*dconjg(Jamp(i,n))))
 		  ENDIF
 		enddo
 	  enddo

=== modified file 'madgraph/iolibs/ufo_expression_parsers.py'
--- madgraph/iolibs/ufo_expression_parsers.py	2014-07-03 18:01:46 +0000
+++ madgraph/iolibs/ufo_expression_parsers.py	2014-11-06 06:25:51 +0000
@@ -61,7 +61,7 @@
     # List of tokens and literals
     tokens = (
         'LOGICAL','LOGICALCOMB','POWER', 'CSC', 'SEC', 'ACSC', 'ASEC',
-        'SQRT', 'CONJ', 'RE', 'IM', 'PI', 'COMPLEX', 'FUNCTION', 'IF','ELSE',
+        'SQRT', 'CONJ', 'RE', 'RE2', 'IM', 'PI', 'COMPLEX', 'FUNCTION', 'IF','ELSE',
         'VARIABLE', 'NUMBER','COND','REGLOG', 'ARG'
         )
     literals = "=+-*/(),"
@@ -115,6 +115,10 @@
     def t_RE(self, t):
         r'(?<!\w)re(?=\()'
         return t
+    def t_RE2(self, t):
+        r'\.real|\.imag'
+        return t
+    
     def t_COMPLEX(self, t):
         r'(?<!\w)complex(?=\()'
         return t
@@ -125,7 +129,7 @@
         r'[a-zA-Z_][0-9a-zA-Z_]*'
         return t
 
-    t_NUMBER = r'([0-9]+\.[0-9]*|\.[0-9]+|[0-9]+)([eE][+-]{0,1}[0-9]+){0,1}'
+    t_NUMBER = r'([0-9]+\.[0-9]*|\.[0-9]+|[0-9]+)([eE][+-]{0,1}[0-9]+){0,1}j{0,1}'
     t_POWER  = r'\*\*'
 
     t_ignore = " \t"
@@ -154,6 +158,7 @@
         ('left','='),
         ('left','+','-'),
         ('left','*','/'),
+        ('left', 'RE2'),
         ('right','UMINUS'),
         ('left','POWER'),
         ('right','REGLOG'),
@@ -255,7 +260,10 @@
 
     def p_expression_number(self, p):
         "expression : NUMBER"
-        p[0] = ('%e' % float(p[1])).replace('e', 'd')
+        if p[1].endswith('j'):
+            p[0] = ('DCOMPLX(0d0, %e)' % float(p[1][:-1])).replace('e', 'd')
+        else:
+            p[0] = ('%e' % float(p[1])).replace('e', 'd')
 
     def p_expression_variable(self, p):
         "expression : VARIABLE"
@@ -313,6 +321,21 @@
         elif p[1] == 'complexconjugate': p[0] = 'conjg(DCMPLX' + p[2]+')'
         elif p[1] == 'reglog': p[0] = 'reglog(DCMPLX' + p[2] +')'
 
+
+    def p_expression_real(self, p):
+        ''' expression : expression RE2 '''
+        
+        if p[2] == '.real':
+            if p[1].startswith('('):
+                p[0] = 'dble' +p[1]
+            else:
+                p[0] = 'dble(%s)' % p[1]
+        elif p[2] == '.imag':
+            if p[1].startswith('('):
+                p[0] = 'dimag' +p[1]
+            else:
+                p[0] = 'dimag(%s)' % p[1]            
+
     def p_expression_pi(self, p):
         '''expression : PI'''
         p[0] = 'pi'
@@ -328,7 +351,11 @@
 
     def p_expression_number(self, p):
         "expression : NUMBER"
-        p[0] = '%e_16' % float(p[1])
+        
+        if p[1].endswith('j'):
+            p[0] = 'CMPLX(0.000000e+00_16, %e_16 ,KIND=16)' % float(p[1][:-1])
+        else:
+            p[0] = '%e_16' % float(p[1])
 
     def p_expression_variable(self, p):
         "expression : VARIABLE"
@@ -384,6 +411,21 @@
         elif p[1] == 'complexconjugate': p[0] = 'conjg(CMPLX(' + p[2] + ',KIND=16))'
         elif p[1] == 'reglog': p[0] = 'mp_reglog(CMPLX(' + p[2] +',KIND=16))'
 
+    def p_expression_real(self, p):
+        ''' expression : expression RE2 '''
+        
+        if p[2] == '.real':
+            if p[1].startswith('('):
+                p[0] = 'real' +p[1]
+            else:
+                p[0] = 'real(%s)' % p[1]
+        elif p[2] == '.imag':
+            if p[1].startswith('('):
+                p[0] = 'imag' +p[1]
+            else:
+                p[0] = 'imag(%s)' % p[1]  
+
+
     def p_expression_pi(self, p):
         '''expression : PI'''
         p[0] = self.mp_prefix+'pi'
@@ -406,6 +448,13 @@
 
     def p_expression_number(self, p):
         'expression : NUMBER'
+        
+        if p[1].endswith('j'):
+            p[0] = 'std::complex<double>(0., %e)'  % float(p[1][:-1]) 
+        else:
+            p[0] = ('%e' % float(p[1])).replace('e', 'd')
+        
+        
         p[0] = p[1]
         # Check number is an integer, if so add "."
         if float(p[1]) == int(float(p[1])) and float(p[1]) < 1000:
@@ -462,6 +511,21 @@
         elif p[1] == 'cmath.sqrt' or p[1] == 'sqrt': p[0] = 'sqrt' + p[2]
         elif p[1] == 'complexconjugate': p[0] = 'conj' + p[2]
         elif p[1] == 'reglog': p[0] = 'reglog' + p[2]
+
+    def p_expression_real(self, p):
+        ''' expression : expression RE2 '''
+        
+        if p[2] == '.real':
+            if p[1].startswith('('):
+                p[0] = 'real' +p[1]
+            else:
+                p[0] = 'real(%s)' % p[1]
+        elif p[2] == '.imag':
+            if p[1].startswith('('):
+                p[0] = 'imag' +p[1]
+            else:
+                p[0] = 'imag(%s)' % p[1]    
+
     
     def p_expression_pi(self, p):
         '''expression : PI'''

=== modified file 'madgraph/loop/loop_diagram_generation.py'
--- madgraph/loop/loop_diagram_generation.py	2014-08-09 09:07:56 +0000
+++ madgraph/loop/loop_diagram_generation.py	2014-11-06 06:25:51 +0000
@@ -650,7 +650,6 @@
 
         # Now we can generate the loop diagrams.
         totloopsuccessful=self.generate_loop_diagrams()
-        self['process']['forbidden_particles']=[]
 
         # If there is no born neither loop diagrams, return now.
         if not self['process']['has_born'] and not self['loop_diagrams']:

=== modified file 'madgraph/various/banner.py'
--- madgraph/various/banner.py	2014-09-17 07:02:56 +0000
+++ madgraph/various/banner.py	2014-11-06 06:25:51 +0000
@@ -588,8 +588,15 @@
                 return '.false.'
             
         elif format == 'int':
-            return str(int(value))
-        
+            try:
+                return str(int(value))
+            except ValueError:
+                fl = float(value)
+                if int(fl) == fl:
+                    return str(int(fl))
+                else:
+                    raise
+                
         elif format == 'float':
             if isinstance(value, str):
                 value = value.replace('d','e')

=== modified file 'madgraph/various/cluster.py'
--- madgraph/various/cluster.py	2014-09-14 18:16:47 +0000
+++ madgraph/various/cluster.py	2014-11-06 06:25:51 +0000
@@ -1290,7 +1290,7 @@
             me_dir = 'a' + me_dir[1:]
         
         text = ""
-        command = ['bsub', '-J', me_dir]
+        command = ['bsub', '-C0', '-J', me_dir]
         if cwd is None:
             cwd = os.getcwd()
         else: 
@@ -1364,7 +1364,7 @@
         cmd = "bjobs " + ' '.join(self.submitted_ids) 
         status = misc.Popen([cmd], shell=True, stdout=subprocess.PIPE)
 
-        idle, run, fail = 0, 0, 0
+        jobstatus = {}
         for line in status.stdout:
             line = line.strip()
             if 'JOBID' in line:
@@ -1373,19 +1373,24 @@
             id = splitline[0]
             if id not in self.submitted_ids:
                 continue
-            status = splitline[2]
+            jobstatus[id] = splitline[2]
+
+        idle, run, fail = 0, 0, 0
+        for id in self.submitted_ids[:]:
+            if id in jobstatus:
+                status = jobstatus[id]
+            else:
+                status = 'MISSING'
             if status == 'RUN':
                 run += 1
             elif status == 'PEND':
                 idle += 1
-            elif status == 'DONE':
+            else:
                 status = self.check_termination(id)
                 if status == 'wait':
                     run += 1
                 elif status == 'resubmit':
-                    idle += 1
-            else:
-                fail += 1
+                    idle += 1                
 
         return idle, run, self.submitted - (idle+run+fail), fail
 

=== modified file 'madgraph/various/lhe_parser.py'
--- madgraph/various/lhe_parser.py	2014-07-07 12:59:46 +0000
+++ madgraph/various/lhe_parser.py	2014-11-06 06:25:51 +0000
@@ -1,6 +1,9 @@
 import collections
 import re
 import misc
+if '__main__' == __name__:
+    import sys
+    sys.path.append('../../')
 
 import logging
 logger = logging.getLogger("madgraph.lhe_parser")

=== modified file 'madgraph/various/misc.py'
--- madgraph/various/misc.py	2014-09-17 07:02:56 +0000
+++ madgraph/various/misc.py	2014-11-06 06:25:51 +0000
@@ -710,21 +710,44 @@
         else:
             raise Exception, "%(path)s does not finish by .gz and the file %(path)s.gz does not exists" %\
                               {"path": path}         
+
+    
+    #for large file (>1G) it is faster and safer to use a separate thread
+    if os.path.getsize(path) > 1e8:
+        if stdout:
+            os.system('gunzip -c %s > %s' % (path, stdout))
+        else:
+            os.system('gunzip  %s') 
+        return
+    
     if not stdout:
-        stdout = path[:-3]
+        stdout = path[:-3]        
     open(stdout,'w').write(ziplib.open(path, "r").read())
     if not keep:
         os.remove(path)
 
-def gzip(path, stdout=None, error=True):
+def gzip(path, stdout=None, error=True, forceexternal=False):
     """ a standard replacement for os.system('gzip %s ' % path)"""
+
+
+    
+    #for large file (>1G) it is faster and safer to use a separate thread
+    if os.path.getsize(path) > 1e9:
+        call(['gzip', '-f', path])
+        if stdout:
+            shutil.move('%s.gz' % path, stdout)
+        return
+    
     if not stdout:
         stdout = "%s.gz" % path
     elif not stdout.endswith(".gz"):
         stdout = "%s.gz" % stdout
+
     try:
         ziplib.open(stdout,"w").write(open(path).read())
-    except:
+    except OverflowError:
+        gzip(path, stdout, error=error, forceexternal=True)
+    except Exception:
         if error:
             raise
     else:

=== modified file 'mg5decay/decay_objects.py'
--- mg5decay/decay_objects.py	2014-07-10 23:42:17 +0000
+++ mg5decay/decay_objects.py	2014-11-06 06:25:51 +0000
@@ -943,13 +943,33 @@
                     vlist_a = inter_part.get_vertexlist(vlevel, True)
                     vlist_b = inter_part.get_vertexlist(vlevel, False)
 
-
+                    minv_max = eval(self['mass']) - \
+                                    sum([eval(model.get_particle(abs(l['id']))['mass'])
+                                         for l in sub_c.get_final_legs() if l!=leg])
+                    
+                    allow_qcd=True # if two colored particle are lower than the pion. 
+                                   # those computation are meaningless.
+                    if 0 < minv_max.real < 0.100:
+                        logger.warning("WARNING: Mass gap lower than pion mass for decay of %s "
+                                                             % self['pdg_code'])
+                        logger.warning("decay into colored particle will be remove.")
+                        allow_qcd=False
+                        if model.get_particle(abs(leg['id']))['color'] != 1:
+                            continue
+                        
+                        
                     # Find appropriate vertex
                     for vert in (vlist_a + vlist_b):
                         # Connect sub_channel to the vertex
                         # the connect_channel_vertex will
                         # inherit the 'has_idpart' from sub_c
 
+                        if not allow_qcd:
+                            nb_colored = sum([1 for l in vert['legs'] if \
+                                     model.get_particle(abs(l['id']))['color'] != 1])
+                            if nb_colored >=1:
+                                continue
+
                         temp_c = self.connect_channel_vertex(sub_c, index, 
                                                              vert, model)
                         temp_c_o = temp_c.get_onshell(model)
@@ -1198,7 +1218,10 @@
         for amp in list(self.get_amplitudes(clevel)):
             approx_width = amp.get('apx_decaywidth')
             if min_br:
-                br = approx_width / total_width
+                if approx_width or total_width:
+                    br = approx_width / total_width
+                else:
+                    br = 0
                 if br.real < min_br:
                     self.decay_amplitudes[clevel].remove(amp)
                 

=== modified file 'models/import_ufo.py'
--- models/import_ufo.py	2014-09-17 15:44:51 +0000
+++ models/import_ufo.py	2014-11-06 06:25:51 +0000
@@ -151,12 +151,18 @@
         
     # Check the validity of the model
     files_list_prov = ['couplings.py','lorentz.py','parameters.py',
-                       'particles.py', 'vertices.py']
+                       'particles.py', 'vertices.py', 'function_library.py',
+                       'propagators.py' ]
+    
+    if decay:
+        files_list_prov.append('decays.py')    
+    
     files_list = []
     for filename in files_list_prov:
         filepath = os.path.join(model_path, filename)
         if not os.path.isfile(filepath):
-            raise UFOImportError,  "%s directory is not a valid UFO model: \n %s is missing" % \
+            if filename not in ['propagators.py', 'decays.py']:
+                raise UFOImportError,  "%s directory is not a valid UFO model: \n %s is missing" % \
                                                          (model_path, filename)
         files_list.append(filepath)
     

=== modified file 'models/taudecay_UFO/param_card.dat'
--- models/taudecay_UFO/param_card.dat	2013-12-20 11:49:27 +0000
+++ models/taudecay_UFO/param_card.dat	2014-11-06 06:25:51 +0000
@@ -12,6 +12,8 @@
 ## INFORMATION FOR MASS
 ###################################
 Block MASS 
+   11 5.110000e-04 # Me 
+   13 1.056600e-01 # MMU 
    15 1.776820e+00 # MTA 
   111 1.349766e-01 # Mpi0 
   211 1.395702e-01 # Mpic 
@@ -19,6 +21,8 @@
 ## Those values should be edited following analytical the 
 ## analytical expression. Some generator could simply ignore 
 ## those values and use the analytical expression
+  12 0.000000 # ve : 0.0 
+  14 0.000000 # vm : 0.0 
   16 0.000000 # vt : 0.0 
 
 ###################################
@@ -29,7 +33,11 @@
 ## Those values should be edited following analytical the 
 ## analytical expression. Some generator could simply ignore 
 ## those values and use the analytical expression
+DECAY  12 0.000000 # ve : 0.0 
+DECAY  14 0.000000 # vm : 0.0 
 DECAY  16 0.000000 # vt : 0.0 
+DECAY  11 0.000000 # e- : 0.0 
+DECAY  13 0.000000 # mu- : 0.0 
 DECAY  211 0.000000 # pi+ : 0.0 
 DECAY  111 0.000000 # pi0 : 0.0 
 
@@ -40,16 +48,17 @@
     1 2.277360e-01 # cabi 
 
 ###################################
-## INFORMATION FOR FRBLOCK
+## INFORMATION FOR TAUDECAY
 ###################################
-Block FRBlock 
-    1 1.304100e-01 # F1 
-    2 1.000000e+00 # F2 
-    3 1.000000e+00 # F3 
-    4 1.000000e+00 # Fr1 
-    5 1.000000e+00 # Fr2 
-    6 1.000000e+00 # G1 
-    7 1.000000e+00 # G2 
+Block taudecay 
+    1 1.000000e+00 # F0 
+    2 1.304100e-01 # F1 
+    3 1.000000e+00 # F2 
+    4 1.000000e+00 # F3 
+    5 1.000000e+00 # Fr1 
+    6 1.000000e+00 # Fr2 
+    7 1.000000e+00 # Gr1 
+    8 1.000000e+00 # Gr2 
 #===========================================================
 # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
 #===========================================================

=== modified file 'models/taudecay_UFO/parameters.py'
--- models/taudecay_UFO/parameters.py	2014-05-09 08:09:59 +0000
+++ models/taudecay_UFO/parameters.py	2014-11-06 06:25:51 +0000
@@ -38,7 +38,7 @@
                type = 'real',
                value = 1.,
                texname = 'F_0',
-               lhablock = 'FRBlock',
+               lhablock = 'taudecay',
                lhacode = [ 1 ])
 
 F1 = Parameter(name = 'F1',
@@ -46,7 +46,7 @@
                type = 'real',
                value = 0.13041,
                texname = 'F_1',
-               lhablock = 'FRBlock',
+               lhablock = 'taudecay',
                lhacode = [ 2 ])
 
 F2 = Parameter(name = 'F2',
@@ -54,7 +54,7 @@
                type = 'real',
                value = 1.,
                texname = 'F_2',
-               lhablock = 'FRBlock',
+               lhablock = 'taudecay',
                lhacode = [ 3 ])
 
 F3 = Parameter(name = 'F3',
@@ -62,7 +62,7 @@
                type = 'real',
                value = 1.,
                texname = 'F_3',
-               lhablock = 'FRBlock',
+               lhablock = 'taudecay',
                lhacode = [ 4 ])
 
 Fr1 = Parameter(name = 'Fr1',
@@ -70,7 +70,7 @@
                 type = 'real',
                 value = 1.,
                 texname = 'F_{23}',
-                lhablock = 'FRBlock',
+                lhablock = 'taudecay',
                 lhacode = [ 5 ])
 
 Fr2 = Parameter(name = 'Fr2',
@@ -78,7 +78,7 @@
                 type = 'real',
                 value = 1.,
                 texname = 'F_{13}',
-                lhablock = 'FRBlock',
+                lhablock = 'taudecay',
                 lhacode = [ 6 ])
 
 Gr1 = Parameter(name = 'Gr1',
@@ -86,7 +86,7 @@
                 type = 'real',
                 value = 1.,
                 texname = 'G_{23}',
-                lhablock = 'FRBlock',
+                lhablock = 'taudecay',
                 lhacode = [ 7 ])
 
 Gr2 = Parameter(name = 'Gr2',
@@ -94,7 +94,7 @@
                 type = 'real',
                 value = 1.,
                 texname = 'G_{13}',
-                lhablock = 'FRBlock',
+                lhablock = 'taudecay',
                 lhacode = [ 8 ])
 
 Me = Parameter(name = 'Me',

=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,12 @@
+      INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
+      PARAMETER (LMAXCONFIGS_USED=  15)
+      PARAMETER (MAX_BRANCH_USED =   3)
+      INTEGER MAPCONFIG_D(  8,0:LMAXCONFIGS_USED)
+      INTEGER IFOREST_D(  8,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER SPROP_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER TPRID_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PMASS_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED
+     $ )
+      DOUBLE PRECISION PWIDTH_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USE
+     $ D)
+      INTEGER POW_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)

=== removed file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_info.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_info.inc	2014-03-19 08:01:45 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%configs_and_props_info.inc	1970-01-01 00:00:00 +0000
@@ -1,1077 +0,0 @@
-      INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
-      PARAMETER (LMAXCONFIGS_USED=  15)
-      PARAMETER (MAX_BRANCH_USED =   3)
-      INTEGER MAPCONFIG_D(  8,0:LMAXCONFIGS_USED)
-      INTEGER IFOREST_D(  8,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER SPROP_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER TPRID_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PMASS_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED
-     $ )
-      DOUBLE PRECISION PWIDTH_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USE
-     $ D)
-      INTEGER POW_D(  8,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-C     
-C     nFKSprocess 1
-C     Diagram 1 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   1)/   1/
-      DATA (IFOREST_D(  1, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,   1)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  1,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   2)/   2/
-      DATA (IFOREST_D(  1, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  1,  -1,   2)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  1,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   3)/   3/
-      DATA (IFOREST_D(  1, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  1,  -1,   3)/      -6/
-      DATA (IFOREST_D(  1, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  1,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   4)/   4/
-      DATA (IFOREST_D(  1, IFR, -1,   4),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  1,  -1,   4)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   4),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  1,  -2,   4)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,   4),IFR=1,2)/ -2,  4/
-C     Diagram 5 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   5)/   5/
-      DATA (IFOREST_D(  1, IFR, -1,   5),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  1,  -1,   5)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   5),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  1,  -2,   5)/      21/
-      DATA (IFOREST_D(  1, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Diagram 6 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   6)/   6/
-      DATA (IFOREST_D(  1, IFR, -1,   6),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  1,  -1,   6)/      -6/
-      DATA (IFOREST_D(  1, IFR, -2,   6),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  1,  -2,   6)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,   6),IFR=1,2)/ -2, -1/
-C     Diagram 7 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   7)/   7/
-      DATA (IFOREST_D(  1, IFR, -1,   7),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  1,  -1,   7)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   7),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  1,  -2,   7)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,   7),IFR=1,2)/ -2,  3/
-C     Diagram 8 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   8)/   8/
-      DATA (IFOREST_D(  1, IFR, -1,   8),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  1,  -1,   8)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   8),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  1,  -2,   8)/      21/
-      DATA (IFOREST_D(  1, IFR, -3,   8),IFR=1,2)/ -2,  5/
-C     Diagram 9 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   9)/   9/
-      DATA (IFOREST_D(  1, IFR, -1,   9),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  1,  -1,   9)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   9),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  1,  -2,   9)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,   9),IFR=1,2)/ -2, -1/
-C     Diagram 10 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  10)/  10/
-      DATA (IFOREST_D(  1, IFR, -1,  10),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  1,  -1,  10)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,  10),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  1,  -2,  10)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,  10),IFR=1,2)/ -2,  3/
-C     Diagram 11 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  11)/  11/
-      DATA (IFOREST_D(  1, IFR, -1,  11),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  1,  -1,  11)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,  11),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  1,  -2,  11)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,  11),IFR=1,2)/ -2,  4/
-C     Diagram 12 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  12)/  12/
-      DATA (IFOREST_D(  1, IFR, -1,  12),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,  12)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,  12),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  1,  -2,  12)/      21/
-      DATA (IFOREST_D(  1, IFR, -3,  12),IFR=1,2)/ -2, -1/
-C     Diagram 13 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  13)/  13/
-      DATA (IFOREST_D(  1, IFR, -1,  13),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  1,  -1,  13)/      -6/
-      DATA (IFOREST_D(  1, IFR, -2,  13),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  1,  -2,  13)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,  13),IFR=1,2)/ -2,  3/
-C     Diagram 14 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  14)/  14/
-      DATA (IFOREST_D(  1, IFR, -1,  14),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  1,  -1,  14)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,  14),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  1,  -2,  14)/       6/
-      DATA (IFOREST_D(  1, IFR, -3,  14),IFR=1,2)/ -2,  4/
-C     Diagram 15 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,  15)/  15/
-      DATA (IFOREST_D(  1, IFR, -1,  15),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,  15)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,  15),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  1,  -2,  15)/      21/
-      DATA (IFOREST_D(  1, IFR, -3,  15),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,0)/  15/
-C     
-C     nFKSprocess 2
-C     Diagram 1 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   1)/   1/
-      DATA (IFOREST_D(  2, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,   1)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  2,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   2)/   2/
-      DATA (IFOREST_D(  2, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  2,  -1,   2)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  2,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   3)/   3/
-      DATA (IFOREST_D(  2, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  2,  -1,   3)/      -6/
-      DATA (IFOREST_D(  2, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  2,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   4)/   4/
-      DATA (IFOREST_D(  2, IFR, -1,   4),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  2,  -1,   4)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   4),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  2,  -2,   4)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,   4),IFR=1,2)/ -2,  4/
-C     Diagram 5 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   5)/   5/
-      DATA (IFOREST_D(  2, IFR, -1,   5),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  2,  -1,   5)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   5),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  2,  -2,   5)/      21/
-      DATA (IFOREST_D(  2, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Diagram 6 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   6)/   6/
-      DATA (IFOREST_D(  2, IFR, -1,   6),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  2,  -1,   6)/      -6/
-      DATA (IFOREST_D(  2, IFR, -2,   6),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  2,  -2,   6)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,   6),IFR=1,2)/ -2, -1/
-C     Diagram 7 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   7)/   7/
-      DATA (IFOREST_D(  2, IFR, -1,   7),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  2,  -1,   7)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   7),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  2,  -2,   7)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,   7),IFR=1,2)/ -2,  3/
-C     Diagram 8 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   8)/   8/
-      DATA (IFOREST_D(  2, IFR, -1,   8),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  2,  -1,   8)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   8),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  2,  -2,   8)/      21/
-      DATA (IFOREST_D(  2, IFR, -3,   8),IFR=1,2)/ -2,  5/
-C     Diagram 9 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   9)/   9/
-      DATA (IFOREST_D(  2, IFR, -1,   9),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  2,  -1,   9)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   9),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  2,  -2,   9)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,   9),IFR=1,2)/ -2, -1/
-C     Diagram 10 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  10)/  10/
-      DATA (IFOREST_D(  2, IFR, -1,  10),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  2,  -1,  10)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,  10),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  2,  -2,  10)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,  10),IFR=1,2)/ -2,  3/
-C     Diagram 11 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  11)/  11/
-      DATA (IFOREST_D(  2, IFR, -1,  11),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  2,  -1,  11)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,  11),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  2,  -2,  11)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,  11),IFR=1,2)/ -2,  4/
-C     Diagram 12 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  12)/  12/
-      DATA (IFOREST_D(  2, IFR, -1,  12),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,  12)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,  12),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  2,  -2,  12)/      21/
-      DATA (IFOREST_D(  2, IFR, -3,  12),IFR=1,2)/ -2, -1/
-C     Diagram 13 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  13)/  13/
-      DATA (IFOREST_D(  2, IFR, -1,  13),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  2,  -1,  13)/      -6/
-      DATA (IFOREST_D(  2, IFR, -2,  13),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  2,  -2,  13)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,  13),IFR=1,2)/ -2,  3/
-C     Diagram 14 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  14)/  14/
-      DATA (IFOREST_D(  2, IFR, -1,  14),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  2,  -1,  14)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,  14),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  2,  -2,  14)/       6/
-      DATA (IFOREST_D(  2, IFR, -3,  14),IFR=1,2)/ -2,  4/
-C     Diagram 15 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,  15)/  15/
-      DATA (IFOREST_D(  2, IFR, -1,  15),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,  15)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,  15),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  2,  -2,  15)/      21/
-      DATA (IFOREST_D(  2, IFR, -3,  15),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,0)/  15/
-C     
-C     nFKSprocess 3
-C     Diagram 1 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   1)/   1/
-      DATA (IFOREST_D(  3, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,   1)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  3,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   2)/   2/
-      DATA (IFOREST_D(  3, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  3,  -1,   2)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  3,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   3)/   3/
-      DATA (IFOREST_D(  3, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  3,  -1,   3)/      -6/
-      DATA (IFOREST_D(  3, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  3,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   4)/   4/
-      DATA (IFOREST_D(  3, IFR, -1,   4),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  3,  -1,   4)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   4),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  3,  -2,   4)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,   4),IFR=1,2)/ -2,  4/
-C     Diagram 5 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   5)/   5/
-      DATA (IFOREST_D(  3, IFR, -1,   5),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  3,  -1,   5)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   5),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  3,  -2,   5)/      21/
-      DATA (IFOREST_D(  3, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Diagram 6 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   6)/   6/
-      DATA (IFOREST_D(  3, IFR, -1,   6),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  3,  -1,   6)/      -6/
-      DATA (IFOREST_D(  3, IFR, -2,   6),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  3,  -2,   6)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,   6),IFR=1,2)/ -2, -1/
-C     Diagram 7 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   7)/   7/
-      DATA (IFOREST_D(  3, IFR, -1,   7),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  3,  -1,   7)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   7),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  3,  -2,   7)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,   7),IFR=1,2)/ -2,  3/
-C     Diagram 8 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   8)/   8/
-      DATA (IFOREST_D(  3, IFR, -1,   8),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  3,  -1,   8)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   8),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  3,  -2,   8)/      21/
-      DATA (IFOREST_D(  3, IFR, -3,   8),IFR=1,2)/ -2,  5/
-C     Diagram 9 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   9)/   9/
-      DATA (IFOREST_D(  3, IFR, -1,   9),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  3,  -1,   9)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   9),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  3,  -2,   9)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,   9),IFR=1,2)/ -2, -1/
-C     Diagram 10 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  10)/  10/
-      DATA (IFOREST_D(  3, IFR, -1,  10),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  3,  -1,  10)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,  10),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  3,  -2,  10)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,  10),IFR=1,2)/ -2,  3/
-C     Diagram 11 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  11)/  11/
-      DATA (IFOREST_D(  3, IFR, -1,  11),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  3,  -1,  11)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,  11),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  3,  -2,  11)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,  11),IFR=1,2)/ -2,  4/
-C     Diagram 12 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  12)/  12/
-      DATA (IFOREST_D(  3, IFR, -1,  12),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,  12)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,  12),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  3,  -2,  12)/      21/
-      DATA (IFOREST_D(  3, IFR, -3,  12),IFR=1,2)/ -2, -1/
-C     Diagram 13 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  13)/  13/
-      DATA (IFOREST_D(  3, IFR, -1,  13),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  3,  -1,  13)/      -6/
-      DATA (IFOREST_D(  3, IFR, -2,  13),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  3,  -2,  13)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,  13),IFR=1,2)/ -2,  3/
-C     Diagram 14 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  14)/  14/
-      DATA (IFOREST_D(  3, IFR, -1,  14),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  3,  -1,  14)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,  14),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  3,  -2,  14)/       6/
-      DATA (IFOREST_D(  3, IFR, -3,  14),IFR=1,2)/ -2,  4/
-C     Diagram 15 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,  15)/  15/
-      DATA (IFOREST_D(  3, IFR, -1,  15),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,  15)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,  15),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  3,  -2,  15)/      21/
-      DATA (IFOREST_D(  3, IFR, -3,  15),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,0)/  15/
-C     
-C     nFKSprocess 4
-C     Diagram 1 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   1)/   1/
-      DATA (IFOREST_D(  4, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,   1)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  4,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   2)/   2/
-      DATA (IFOREST_D(  4, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  4,  -1,   2)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  4,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   3)/   3/
-      DATA (IFOREST_D(  4, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  4,  -1,   3)/      -6/
-      DATA (IFOREST_D(  4, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  4,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   4)/   4/
-      DATA (IFOREST_D(  4, IFR, -1,   4),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  4,  -1,   4)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   4),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  4,  -2,   4)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,   4),IFR=1,2)/ -2,  4/
-C     Diagram 5 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   5)/   5/
-      DATA (IFOREST_D(  4, IFR, -1,   5),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  4,  -1,   5)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   5),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  4,  -2,   5)/      21/
-      DATA (IFOREST_D(  4, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Diagram 6 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   6)/   6/
-      DATA (IFOREST_D(  4, IFR, -1,   6),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  4,  -1,   6)/      -6/
-      DATA (IFOREST_D(  4, IFR, -2,   6),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  4,  -2,   6)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,   6),IFR=1,2)/ -2, -1/
-C     Diagram 7 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   7)/   7/
-      DATA (IFOREST_D(  4, IFR, -1,   7),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  4,  -1,   7)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   7),IFR=1,2)/ -1,  5/
-      DATA TPRID_D(  4,  -2,   7)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,   7),IFR=1,2)/ -2,  3/
-C     Diagram 8 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   8)/   8/
-      DATA (IFOREST_D(  4, IFR, -1,   8),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  4,  -1,   8)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   8),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  4,  -2,   8)/      21/
-      DATA (IFOREST_D(  4, IFR, -3,   8),IFR=1,2)/ -2,  5/
-C     Diagram 9 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   9)/   9/
-      DATA (IFOREST_D(  4, IFR, -1,   9),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  4,  -1,   9)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   9),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  4,  -2,   9)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,   9),IFR=1,2)/ -2, -1/
-C     Diagram 10 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  10)/  10/
-      DATA (IFOREST_D(  4, IFR, -1,  10),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  4,  -1,  10)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,  10),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  4,  -2,  10)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,  10),IFR=1,2)/ -2,  3/
-C     Diagram 11 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  11)/  11/
-      DATA (IFOREST_D(  4, IFR, -1,  11),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  4,  -1,  11)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,  11),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  4,  -2,  11)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,  11),IFR=1,2)/ -2,  4/
-C     Diagram 12 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  12)/  12/
-      DATA (IFOREST_D(  4, IFR, -1,  12),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,  12)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,  12),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  4,  -2,  12)/      21/
-      DATA (IFOREST_D(  4, IFR, -3,  12),IFR=1,2)/ -2, -1/
-C     Diagram 13 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  13)/  13/
-      DATA (IFOREST_D(  4, IFR, -1,  13),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  4,  -1,  13)/      -6/
-      DATA (IFOREST_D(  4, IFR, -2,  13),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  4,  -2,  13)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,  13),IFR=1,2)/ -2,  3/
-C     Diagram 14 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  14)/  14/
-      DATA (IFOREST_D(  4, IFR, -1,  14),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  4,  -1,  14)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,  14),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  4,  -2,  14)/       6/
-      DATA (IFOREST_D(  4, IFR, -3,  14),IFR=1,2)/ -2,  4/
-C     Diagram 15 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,  15)/  15/
-      DATA (IFOREST_D(  4, IFR, -1,  15),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,  15)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,  15),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  4,  -2,  15)/      21/
-      DATA (IFOREST_D(  4, IFR, -3,  15),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,0)/  15/
-C     
-C     nFKSprocess 5
-C     Diagram 1 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   1)/   1/
-      DATA (IFOREST_D(  5, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  5,  -1,   1)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  5,  -2,   1)/      -1/
-C     Diagram 2 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   2)/   2/
-      DATA (IFOREST_D(  5, IFR, -1,   2),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  5,  -1,   2)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   2),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  5,  -2,   2)/       6/
-      DATA (IFOREST_D(  5, IFR, -3,   2),IFR=1,2)/ -2,  3/
-C     Diagram 3 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   3)/   3/
-      DATA (IFOREST_D(  5, IFR, -1,   3),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  5,  -1,   3)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  5,  -2,   3)/       6/
-      DATA (IFOREST_D(  5, IFR, -3,   3),IFR=1,2)/ -2,  4/
-C     Diagram 4 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   4)/   4/
-      DATA (IFOREST_D(  5, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  5,  -1,   4)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  5,  -2,   4)/      21/
-      DATA (IFOREST_D(  5, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   5)/   5/
-      DATA (IFOREST_D(  5, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  5,  -1,   5)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  5,  -2,   5)/       1/
-      DATA (IFOREST_D(  5, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,0)/   5/
-C     
-C     nFKSprocess 6
-C     Diagram 1 for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,   1)/   1/
-      DATA (IFOREST_D(  6, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  6,  -1,   1)/      21/
-      DATA (IFOREST_D(  6, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  6,  -2,   1)/       1/
-C     Diagram 2 for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,   2)/   2/
-      DATA (IFOREST_D(  6, IFR, -1,   2),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  6,  -1,   2)/      21/
-      DATA (IFOREST_D(  6, IFR, -2,   2),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  6,  -2,   2)/       6/
-      DATA (IFOREST_D(  6, IFR, -3,   2),IFR=1,2)/ -2,  3/
-C     Diagram 3 for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,   3)/   3/
-      DATA (IFOREST_D(  6, IFR, -1,   3),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  6,  -1,   3)/      21/
-      DATA (IFOREST_D(  6, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  6,  -2,   3)/       6/
-      DATA (IFOREST_D(  6, IFR, -3,   3),IFR=1,2)/ -2,  4/
-C     Diagram 4 for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,   4)/   4/
-      DATA (IFOREST_D(  6, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  6,  -1,   4)/      21/
-      DATA (IFOREST_D(  6, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  6,  -2,   4)/      21/
-      DATA (IFOREST_D(  6, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,   5)/   5/
-      DATA (IFOREST_D(  6, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  6,  -1,   5)/      21/
-      DATA (IFOREST_D(  6, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  6,  -2,   5)/       1/
-      DATA (IFOREST_D(  6, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 6
-      DATA MAPCONFIG_D(  6,0)/   5/
-C     
-C     nFKSprocess 7
-C     Diagram 1 for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,   1)/   1/
-      DATA (IFOREST_D(  7, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  7,  -1,   1)/      21/
-      DATA (IFOREST_D(  7, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  7,  -2,   1)/      -1/
-C     Diagram 2 for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,   2)/   2/
-      DATA (IFOREST_D(  7, IFR, -1,   2),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  7,  -1,   2)/       6/
-      DATA (IFOREST_D(  7, IFR, -2,   2),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  7,  -2,   2)/      21/
-      DATA (IFOREST_D(  7, IFR, -3,   2),IFR=1,2)/ -2,  5/
-C     Diagram 3 for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,   3)/   3/
-      DATA (IFOREST_D(  7, IFR, -1,   3),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  7,  -1,   3)/       6/
-      DATA (IFOREST_D(  7, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  7,  -2,   3)/      21/
-      DATA (IFOREST_D(  7, IFR, -3,   3),IFR=1,2)/ -2,  5/
-C     Diagram 4 for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,   4)/   4/
-      DATA (IFOREST_D(  7, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  7,  -1,   4)/      21/
-      DATA (IFOREST_D(  7, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  7,  -2,   4)/       1/
-      DATA (IFOREST_D(  7, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,   5)/   5/
-      DATA (IFOREST_D(  7, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  7,  -1,   5)/      21/
-      DATA (IFOREST_D(  7, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  7,  -2,   5)/      21/
-      DATA (IFOREST_D(  7, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 7
-      DATA MAPCONFIG_D(  7,0)/   5/
-C     
-C     nFKSprocess 8
-C     Diagram 1 for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,   1)/   1/
-      DATA (IFOREST_D(  8, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  8,  -1,   1)/      21/
-      DATA (IFOREST_D(  8, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  8,  -2,   1)/       1/
-C     Diagram 2 for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,   2)/   2/
-      DATA (IFOREST_D(  8, IFR, -1,   2),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  8,  -1,   2)/       6/
-      DATA (IFOREST_D(  8, IFR, -2,   2),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  8,  -2,   2)/      21/
-      DATA (IFOREST_D(  8, IFR, -3,   2),IFR=1,2)/ -2,  5/
-C     Diagram 3 for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,   3)/   3/
-      DATA (IFOREST_D(  8, IFR, -1,   3),IFR=1,2)/  1,  4/
-      DATA TPRID_D(  8,  -1,   3)/       6/
-      DATA (IFOREST_D(  8, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA TPRID_D(  8,  -2,   3)/      21/
-      DATA (IFOREST_D(  8, IFR, -3,   3),IFR=1,2)/ -2,  5/
-C     Diagram 4 for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,   4)/   4/
-      DATA (IFOREST_D(  8, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  8,  -1,   4)/      21/
-      DATA (IFOREST_D(  8, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  8,  -2,   4)/       1/
-      DATA (IFOREST_D(  8, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,   5)/   5/
-      DATA (IFOREST_D(  8, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  8,  -1,   5)/      21/
-      DATA (IFOREST_D(  8, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  8,  -2,   5)/      21/
-      DATA (IFOREST_D(  8, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 8
-      DATA MAPCONFIG_D(  8,0)/   5/
-C     
-C     
-      PMASS_D (  1, -1,   1) = ZERO
-      PWIDTH_D(  1, -1,   1) = ZERO
-      POW_D   (  1, -1,   1) = 2
-      PMASS_D (  1, -2,   1) = ZERO
-      PWIDTH_D(  1, -2,   1) = ZERO
-      POW_D   (  1, -2,   1) = 2
-      PMASS_D (  1, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   2) = ABS(MDL_WT)
-      POW_D   (  1, -1,   2) = 1
-      PMASS_D (  1, -2,   2) = ZERO
-      PWIDTH_D(  1, -2,   2) = ZERO
-      POW_D   (  1, -2,   2) = 2
-      PMASS_D (  1, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   3) = ABS(MDL_WT)
-      POW_D   (  1, -1,   3) = 1
-      PMASS_D (  1, -2,   3) = ZERO
-      PWIDTH_D(  1, -2,   3) = ZERO
-      POW_D   (  1, -2,   3) = 2
-      PMASS_D (  1, -1,   4) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   4) = ABS(MDL_WT)
-      POW_D   (  1, -1,   4) = 1
-      PMASS_D (  1, -2,   4) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,   4) = ABS(MDL_WT)
-      POW_D   (  1, -2,   4) = 1
-      PMASS_D (  1, -1,   5) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   5) = ABS(MDL_WT)
-      POW_D   (  1, -1,   5) = 1
-      PMASS_D (  1, -2,   5) = ZERO
-      PWIDTH_D(  1, -2,   5) = ZERO
-      POW_D   (  1, -2,   5) = 2
-      PMASS_D (  1, -1,   6) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   6) = ABS(MDL_WT)
-      POW_D   (  1, -1,   6) = 1
-      PMASS_D (  1, -2,   6) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,   6) = ABS(MDL_WT)
-      POW_D   (  1, -2,   6) = 1
-      PMASS_D (  1, -1,   7) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   7) = ABS(MDL_WT)
-      POW_D   (  1, -1,   7) = 1
-      PMASS_D (  1, -2,   7) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,   7) = ABS(MDL_WT)
-      POW_D   (  1, -2,   7) = 1
-      PMASS_D (  1, -1,   8) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   8) = ABS(MDL_WT)
-      POW_D   (  1, -1,   8) = 1
-      PMASS_D (  1, -2,   8) = ZERO
-      PWIDTH_D(  1, -2,   8) = ZERO
-      POW_D   (  1, -2,   8) = 2
-      PMASS_D (  1, -1,   9) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,   9) = ABS(MDL_WT)
-      POW_D   (  1, -1,   9) = 1
-      PMASS_D (  1, -2,   9) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,   9) = ABS(MDL_WT)
-      POW_D   (  1, -2,   9) = 1
-      PMASS_D (  1, -1,  10) = ZERO
-      PWIDTH_D(  1, -1,  10) = ZERO
-      POW_D   (  1, -1,  10) = 2
-      PMASS_D (  1, -2,  10) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,  10) = ABS(MDL_WT)
-      POW_D   (  1, -2,  10) = 1
-      PMASS_D (  1, -1,  11) = ZERO
-      PWIDTH_D(  1, -1,  11) = ZERO
-      POW_D   (  1, -1,  11) = 2
-      PMASS_D (  1, -2,  11) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,  11) = ABS(MDL_WT)
-      POW_D   (  1, -2,  11) = 1
-      PMASS_D (  1, -1,  12) = ZERO
-      PWIDTH_D(  1, -1,  12) = ZERO
-      POW_D   (  1, -1,  12) = 2
-      PMASS_D (  1, -2,  12) = ZERO
-      PWIDTH_D(  1, -2,  12) = ZERO
-      POW_D   (  1, -2,  12) = 2
-      PMASS_D (  1, -1,  13) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,  13) = ABS(MDL_WT)
-      POW_D   (  1, -1,  13) = 1
-      PMASS_D (  1, -2,  13) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,  13) = ABS(MDL_WT)
-      POW_D   (  1, -2,  13) = 1
-      PMASS_D (  1, -1,  14) = ABS(MDL_MT)
-      PWIDTH_D(  1, -1,  14) = ABS(MDL_WT)
-      POW_D   (  1, -1,  14) = 1
-      PMASS_D (  1, -2,  14) = ABS(MDL_MT)
-      PWIDTH_D(  1, -2,  14) = ABS(MDL_WT)
-      POW_D   (  1, -2,  14) = 1
-      PMASS_D (  1, -1,  15) = ZERO
-      PWIDTH_D(  1, -1,  15) = ZERO
-      POW_D   (  1, -1,  15) = 2
-      PMASS_D (  1, -2,  15) = ZERO
-      PWIDTH_D(  1, -2,  15) = ZERO
-      POW_D   (  1, -2,  15) = 2
-C     
-      PMASS_D (  2, -1,   1) = ZERO
-      PWIDTH_D(  2, -1,   1) = ZERO
-      POW_D   (  2, -1,   1) = 2
-      PMASS_D (  2, -2,   1) = ZERO
-      PWIDTH_D(  2, -2,   1) = ZERO
-      POW_D   (  2, -2,   1) = 2
-      PMASS_D (  2, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   2) = ABS(MDL_WT)
-      POW_D   (  2, -1,   2) = 1
-      PMASS_D (  2, -2,   2) = ZERO
-      PWIDTH_D(  2, -2,   2) = ZERO
-      POW_D   (  2, -2,   2) = 2
-      PMASS_D (  2, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   3) = ABS(MDL_WT)
-      POW_D   (  2, -1,   3) = 1
-      PMASS_D (  2, -2,   3) = ZERO
-      PWIDTH_D(  2, -2,   3) = ZERO
-      POW_D   (  2, -2,   3) = 2
-      PMASS_D (  2, -1,   4) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   4) = ABS(MDL_WT)
-      POW_D   (  2, -1,   4) = 1
-      PMASS_D (  2, -2,   4) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,   4) = ABS(MDL_WT)
-      POW_D   (  2, -2,   4) = 1
-      PMASS_D (  2, -1,   5) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   5) = ABS(MDL_WT)
-      POW_D   (  2, -1,   5) = 1
-      PMASS_D (  2, -2,   5) = ZERO
-      PWIDTH_D(  2, -2,   5) = ZERO
-      POW_D   (  2, -2,   5) = 2
-      PMASS_D (  2, -1,   6) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   6) = ABS(MDL_WT)
-      POW_D   (  2, -1,   6) = 1
-      PMASS_D (  2, -2,   6) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,   6) = ABS(MDL_WT)
-      POW_D   (  2, -2,   6) = 1
-      PMASS_D (  2, -1,   7) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   7) = ABS(MDL_WT)
-      POW_D   (  2, -1,   7) = 1
-      PMASS_D (  2, -2,   7) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,   7) = ABS(MDL_WT)
-      POW_D   (  2, -2,   7) = 1
-      PMASS_D (  2, -1,   8) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   8) = ABS(MDL_WT)
-      POW_D   (  2, -1,   8) = 1
-      PMASS_D (  2, -2,   8) = ZERO
-      PWIDTH_D(  2, -2,   8) = ZERO
-      POW_D   (  2, -2,   8) = 2
-      PMASS_D (  2, -1,   9) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,   9) = ABS(MDL_WT)
-      POW_D   (  2, -1,   9) = 1
-      PMASS_D (  2, -2,   9) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,   9) = ABS(MDL_WT)
-      POW_D   (  2, -2,   9) = 1
-      PMASS_D (  2, -1,  10) = ZERO
-      PWIDTH_D(  2, -1,  10) = ZERO
-      POW_D   (  2, -1,  10) = 2
-      PMASS_D (  2, -2,  10) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,  10) = ABS(MDL_WT)
-      POW_D   (  2, -2,  10) = 1
-      PMASS_D (  2, -1,  11) = ZERO
-      PWIDTH_D(  2, -1,  11) = ZERO
-      POW_D   (  2, -1,  11) = 2
-      PMASS_D (  2, -2,  11) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,  11) = ABS(MDL_WT)
-      POW_D   (  2, -2,  11) = 1
-      PMASS_D (  2, -1,  12) = ZERO
-      PWIDTH_D(  2, -1,  12) = ZERO
-      POW_D   (  2, -1,  12) = 2
-      PMASS_D (  2, -2,  12) = ZERO
-      PWIDTH_D(  2, -2,  12) = ZERO
-      POW_D   (  2, -2,  12) = 2
-      PMASS_D (  2, -1,  13) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,  13) = ABS(MDL_WT)
-      POW_D   (  2, -1,  13) = 1
-      PMASS_D (  2, -2,  13) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,  13) = ABS(MDL_WT)
-      POW_D   (  2, -2,  13) = 1
-      PMASS_D (  2, -1,  14) = ABS(MDL_MT)
-      PWIDTH_D(  2, -1,  14) = ABS(MDL_WT)
-      POW_D   (  2, -1,  14) = 1
-      PMASS_D (  2, -2,  14) = ABS(MDL_MT)
-      PWIDTH_D(  2, -2,  14) = ABS(MDL_WT)
-      POW_D   (  2, -2,  14) = 1
-      PMASS_D (  2, -1,  15) = ZERO
-      PWIDTH_D(  2, -1,  15) = ZERO
-      POW_D   (  2, -1,  15) = 2
-      PMASS_D (  2, -2,  15) = ZERO
-      PWIDTH_D(  2, -2,  15) = ZERO
-      POW_D   (  2, -2,  15) = 2
-C     
-      PMASS_D (  3, -1,   1) = ZERO
-      PWIDTH_D(  3, -1,   1) = ZERO
-      POW_D   (  3, -1,   1) = 2
-      PMASS_D (  3, -2,   1) = ZERO
-      PWIDTH_D(  3, -2,   1) = ZERO
-      POW_D   (  3, -2,   1) = 2
-      PMASS_D (  3, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   2) = ABS(MDL_WT)
-      POW_D   (  3, -1,   2) = 1
-      PMASS_D (  3, -2,   2) = ZERO
-      PWIDTH_D(  3, -2,   2) = ZERO
-      POW_D   (  3, -2,   2) = 2
-      PMASS_D (  3, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   3) = ABS(MDL_WT)
-      POW_D   (  3, -1,   3) = 1
-      PMASS_D (  3, -2,   3) = ZERO
-      PWIDTH_D(  3, -2,   3) = ZERO
-      POW_D   (  3, -2,   3) = 2
-      PMASS_D (  3, -1,   4) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   4) = ABS(MDL_WT)
-      POW_D   (  3, -1,   4) = 1
-      PMASS_D (  3, -2,   4) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,   4) = ABS(MDL_WT)
-      POW_D   (  3, -2,   4) = 1
-      PMASS_D (  3, -1,   5) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   5) = ABS(MDL_WT)
-      POW_D   (  3, -1,   5) = 1
-      PMASS_D (  3, -2,   5) = ZERO
-      PWIDTH_D(  3, -2,   5) = ZERO
-      POW_D   (  3, -2,   5) = 2
-      PMASS_D (  3, -1,   6) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   6) = ABS(MDL_WT)
-      POW_D   (  3, -1,   6) = 1
-      PMASS_D (  3, -2,   6) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,   6) = ABS(MDL_WT)
-      POW_D   (  3, -2,   6) = 1
-      PMASS_D (  3, -1,   7) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   7) = ABS(MDL_WT)
-      POW_D   (  3, -1,   7) = 1
-      PMASS_D (  3, -2,   7) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,   7) = ABS(MDL_WT)
-      POW_D   (  3, -2,   7) = 1
-      PMASS_D (  3, -1,   8) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   8) = ABS(MDL_WT)
-      POW_D   (  3, -1,   8) = 1
-      PMASS_D (  3, -2,   8) = ZERO
-      PWIDTH_D(  3, -2,   8) = ZERO
-      POW_D   (  3, -2,   8) = 2
-      PMASS_D (  3, -1,   9) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,   9) = ABS(MDL_WT)
-      POW_D   (  3, -1,   9) = 1
-      PMASS_D (  3, -2,   9) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,   9) = ABS(MDL_WT)
-      POW_D   (  3, -2,   9) = 1
-      PMASS_D (  3, -1,  10) = ZERO
-      PWIDTH_D(  3, -1,  10) = ZERO
-      POW_D   (  3, -1,  10) = 2
-      PMASS_D (  3, -2,  10) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,  10) = ABS(MDL_WT)
-      POW_D   (  3, -2,  10) = 1
-      PMASS_D (  3, -1,  11) = ZERO
-      PWIDTH_D(  3, -1,  11) = ZERO
-      POW_D   (  3, -1,  11) = 2
-      PMASS_D (  3, -2,  11) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,  11) = ABS(MDL_WT)
-      POW_D   (  3, -2,  11) = 1
-      PMASS_D (  3, -1,  12) = ZERO
-      PWIDTH_D(  3, -1,  12) = ZERO
-      POW_D   (  3, -1,  12) = 2
-      PMASS_D (  3, -2,  12) = ZERO
-      PWIDTH_D(  3, -2,  12) = ZERO
-      POW_D   (  3, -2,  12) = 2
-      PMASS_D (  3, -1,  13) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,  13) = ABS(MDL_WT)
-      POW_D   (  3, -1,  13) = 1
-      PMASS_D (  3, -2,  13) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,  13) = ABS(MDL_WT)
-      POW_D   (  3, -2,  13) = 1
-      PMASS_D (  3, -1,  14) = ABS(MDL_MT)
-      PWIDTH_D(  3, -1,  14) = ABS(MDL_WT)
-      POW_D   (  3, -1,  14) = 1
-      PMASS_D (  3, -2,  14) = ABS(MDL_MT)
-      PWIDTH_D(  3, -2,  14) = ABS(MDL_WT)
-      POW_D   (  3, -2,  14) = 1
-      PMASS_D (  3, -1,  15) = ZERO
-      PWIDTH_D(  3, -1,  15) = ZERO
-      POW_D   (  3, -1,  15) = 2
-      PMASS_D (  3, -2,  15) = ZERO
-      PWIDTH_D(  3, -2,  15) = ZERO
-      POW_D   (  3, -2,  15) = 2
-C     
-      PMASS_D (  4, -1,   1) = ZERO
-      PWIDTH_D(  4, -1,   1) = ZERO
-      POW_D   (  4, -1,   1) = 2
-      PMASS_D (  4, -2,   1) = ZERO
-      PWIDTH_D(  4, -2,   1) = ZERO
-      POW_D   (  4, -2,   1) = 2
-      PMASS_D (  4, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   2) = ABS(MDL_WT)
-      POW_D   (  4, -1,   2) = 1
-      PMASS_D (  4, -2,   2) = ZERO
-      PWIDTH_D(  4, -2,   2) = ZERO
-      POW_D   (  4, -2,   2) = 2
-      PMASS_D (  4, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   3) = ABS(MDL_WT)
-      POW_D   (  4, -1,   3) = 1
-      PMASS_D (  4, -2,   3) = ZERO
-      PWIDTH_D(  4, -2,   3) = ZERO
-      POW_D   (  4, -2,   3) = 2
-      PMASS_D (  4, -1,   4) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   4) = ABS(MDL_WT)
-      POW_D   (  4, -1,   4) = 1
-      PMASS_D (  4, -2,   4) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,   4) = ABS(MDL_WT)
-      POW_D   (  4, -2,   4) = 1
-      PMASS_D (  4, -1,   5) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   5) = ABS(MDL_WT)
-      POW_D   (  4, -1,   5) = 1
-      PMASS_D (  4, -2,   5) = ZERO
-      PWIDTH_D(  4, -2,   5) = ZERO
-      POW_D   (  4, -2,   5) = 2
-      PMASS_D (  4, -1,   6) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   6) = ABS(MDL_WT)
-      POW_D   (  4, -1,   6) = 1
-      PMASS_D (  4, -2,   6) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,   6) = ABS(MDL_WT)
-      POW_D   (  4, -2,   6) = 1
-      PMASS_D (  4, -1,   7) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   7) = ABS(MDL_WT)
-      POW_D   (  4, -1,   7) = 1
-      PMASS_D (  4, -2,   7) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,   7) = ABS(MDL_WT)
-      POW_D   (  4, -2,   7) = 1
-      PMASS_D (  4, -1,   8) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   8) = ABS(MDL_WT)
-      POW_D   (  4, -1,   8) = 1
-      PMASS_D (  4, -2,   8) = ZERO
-      PWIDTH_D(  4, -2,   8) = ZERO
-      POW_D   (  4, -2,   8) = 2
-      PMASS_D (  4, -1,   9) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,   9) = ABS(MDL_WT)
-      POW_D   (  4, -1,   9) = 1
-      PMASS_D (  4, -2,   9) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,   9) = ABS(MDL_WT)
-      POW_D   (  4, -2,   9) = 1
-      PMASS_D (  4, -1,  10) = ZERO
-      PWIDTH_D(  4, -1,  10) = ZERO
-      POW_D   (  4, -1,  10) = 2
-      PMASS_D (  4, -2,  10) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,  10) = ABS(MDL_WT)
-      POW_D   (  4, -2,  10) = 1
-      PMASS_D (  4, -1,  11) = ZERO
-      PWIDTH_D(  4, -1,  11) = ZERO
-      POW_D   (  4, -1,  11) = 2
-      PMASS_D (  4, -2,  11) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,  11) = ABS(MDL_WT)
-      POW_D   (  4, -2,  11) = 1
-      PMASS_D (  4, -1,  12) = ZERO
-      PWIDTH_D(  4, -1,  12) = ZERO
-      POW_D   (  4, -1,  12) = 2
-      PMASS_D (  4, -2,  12) = ZERO
-      PWIDTH_D(  4, -2,  12) = ZERO
-      POW_D   (  4, -2,  12) = 2
-      PMASS_D (  4, -1,  13) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,  13) = ABS(MDL_WT)
-      POW_D   (  4, -1,  13) = 1
-      PMASS_D (  4, -2,  13) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,  13) = ABS(MDL_WT)
-      POW_D   (  4, -2,  13) = 1
-      PMASS_D (  4, -1,  14) = ABS(MDL_MT)
-      PWIDTH_D(  4, -1,  14) = ABS(MDL_WT)
-      POW_D   (  4, -1,  14) = 1
-      PMASS_D (  4, -2,  14) = ABS(MDL_MT)
-      PWIDTH_D(  4, -2,  14) = ABS(MDL_WT)
-      POW_D   (  4, -2,  14) = 1
-      PMASS_D (  4, -1,  15) = ZERO
-      PWIDTH_D(  4, -1,  15) = ZERO
-      POW_D   (  4, -1,  15) = 2
-      PMASS_D (  4, -2,  15) = ZERO
-      PWIDTH_D(  4, -2,  15) = ZERO
-      POW_D   (  4, -2,  15) = 2
-C     
-      PMASS_D (  5, -1,   1) = ZERO
-      PWIDTH_D(  5, -1,   1) = ZERO
-      POW_D   (  5, -1,   1) = 2
-      PMASS_D (  5, -2,   1) = ZERO
-      PWIDTH_D(  5, -2,   1) = ZERO
-      POW_D   (  5, -2,   1) = 1
-      PMASS_D (  5, -1,   2) = ZERO
-      PWIDTH_D(  5, -1,   2) = ZERO
-      POW_D   (  5, -1,   2) = 2
-      PMASS_D (  5, -2,   2) = ABS(MDL_MT)
-      PWIDTH_D(  5, -2,   2) = ABS(MDL_WT)
-      POW_D   (  5, -2,   2) = 1
-      PMASS_D (  5, -1,   3) = ZERO
-      PWIDTH_D(  5, -1,   3) = ZERO
-      POW_D   (  5, -1,   3) = 2
-      PMASS_D (  5, -2,   3) = ABS(MDL_MT)
-      PWIDTH_D(  5, -2,   3) = ABS(MDL_WT)
-      POW_D   (  5, -2,   3) = 1
-      PMASS_D (  5, -1,   4) = ZERO
-      PWIDTH_D(  5, -1,   4) = ZERO
-      POW_D   (  5, -1,   4) = 2
-      PMASS_D (  5, -2,   4) = ZERO
-      PWIDTH_D(  5, -2,   4) = ZERO
-      POW_D   (  5, -2,   4) = 2
-      PMASS_D (  5, -1,   5) = ZERO
-      PWIDTH_D(  5, -1,   5) = ZERO
-      POW_D   (  5, -1,   5) = 2
-      PMASS_D (  5, -2,   5) = ZERO
-      PWIDTH_D(  5, -2,   5) = ZERO
-      POW_D   (  5, -2,   5) = 1
-C     
-      PMASS_D (  6, -1,   1) = ZERO
-      PWIDTH_D(  6, -1,   1) = ZERO
-      POW_D   (  6, -1,   1) = 2
-      PMASS_D (  6, -2,   1) = ZERO
-      PWIDTH_D(  6, -2,   1) = ZERO
-      POW_D   (  6, -2,   1) = 1
-      PMASS_D (  6, -1,   2) = ZERO
-      PWIDTH_D(  6, -1,   2) = ZERO
-      POW_D   (  6, -1,   2) = 2
-      PMASS_D (  6, -2,   2) = ABS(MDL_MT)
-      PWIDTH_D(  6, -2,   2) = ABS(MDL_WT)
-      POW_D   (  6, -2,   2) = 1
-      PMASS_D (  6, -1,   3) = ZERO
-      PWIDTH_D(  6, -1,   3) = ZERO
-      POW_D   (  6, -1,   3) = 2
-      PMASS_D (  6, -2,   3) = ABS(MDL_MT)
-      PWIDTH_D(  6, -2,   3) = ABS(MDL_WT)
-      POW_D   (  6, -2,   3) = 1
-      PMASS_D (  6, -1,   4) = ZERO
-      PWIDTH_D(  6, -1,   4) = ZERO
-      POW_D   (  6, -1,   4) = 2
-      PMASS_D (  6, -2,   4) = ZERO
-      PWIDTH_D(  6, -2,   4) = ZERO
-      POW_D   (  6, -2,   4) = 2
-      PMASS_D (  6, -1,   5) = ZERO
-      PWIDTH_D(  6, -1,   5) = ZERO
-      POW_D   (  6, -1,   5) = 2
-      PMASS_D (  6, -2,   5) = ZERO
-      PWIDTH_D(  6, -2,   5) = ZERO
-      POW_D   (  6, -2,   5) = 1
-C     
-      PMASS_D (  7, -1,   1) = ZERO
-      PWIDTH_D(  7, -1,   1) = ZERO
-      POW_D   (  7, -1,   1) = 2
-      PMASS_D (  7, -2,   1) = ZERO
-      PWIDTH_D(  7, -2,   1) = ZERO
-      POW_D   (  7, -2,   1) = 1
-      PMASS_D (  7, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  7, -1,   2) = ABS(MDL_WT)
-      POW_D   (  7, -1,   2) = 1
-      PMASS_D (  7, -2,   2) = ZERO
-      PWIDTH_D(  7, -2,   2) = ZERO
-      POW_D   (  7, -2,   2) = 2
-      PMASS_D (  7, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  7, -1,   3) = ABS(MDL_WT)
-      POW_D   (  7, -1,   3) = 1
-      PMASS_D (  7, -2,   3) = ZERO
-      PWIDTH_D(  7, -2,   3) = ZERO
-      POW_D   (  7, -2,   3) = 2
-      PMASS_D (  7, -1,   4) = ZERO
-      PWIDTH_D(  7, -1,   4) = ZERO
-      POW_D   (  7, -1,   4) = 2
-      PMASS_D (  7, -2,   4) = ZERO
-      PWIDTH_D(  7, -2,   4) = ZERO
-      POW_D   (  7, -2,   4) = 1
-      PMASS_D (  7, -1,   5) = ZERO
-      PWIDTH_D(  7, -1,   5) = ZERO
-      POW_D   (  7, -1,   5) = 2
-      PMASS_D (  7, -2,   5) = ZERO
-      PWIDTH_D(  7, -2,   5) = ZERO
-      POW_D   (  7, -2,   5) = 2
-C     
-      PMASS_D (  8, -1,   1) = ZERO
-      PWIDTH_D(  8, -1,   1) = ZERO
-      POW_D   (  8, -1,   1) = 2
-      PMASS_D (  8, -2,   1) = ZERO
-      PWIDTH_D(  8, -2,   1) = ZERO
-      POW_D   (  8, -2,   1) = 1
-      PMASS_D (  8, -1,   2) = ABS(MDL_MT)
-      PWIDTH_D(  8, -1,   2) = ABS(MDL_WT)
-      POW_D   (  8, -1,   2) = 1
-      PMASS_D (  8, -2,   2) = ZERO
-      PWIDTH_D(  8, -2,   2) = ZERO
-      POW_D   (  8, -2,   2) = 2
-      PMASS_D (  8, -1,   3) = ABS(MDL_MT)
-      PWIDTH_D(  8, -1,   3) = ABS(MDL_WT)
-      POW_D   (  8, -1,   3) = 1
-      PMASS_D (  8, -2,   3) = ZERO
-      PWIDTH_D(  8, -2,   3) = ZERO
-      POW_D   (  8, -2,   3) = 2
-      PMASS_D (  8, -1,   4) = ZERO
-      PWIDTH_D(  8, -1,   4) = ZERO
-      POW_D   (  8, -1,   4) = 2
-      PMASS_D (  8, -2,   4) = ZERO
-      PWIDTH_D(  8, -2,   4) = ZERO
-      POW_D   (  8, -2,   4) = 1
-      PMASS_D (  8, -1,   5) = ZERO
-      PWIDTH_D(  8, -1,   5) = ZERO
-      POW_D   (  8, -1,   5) = 2
-      PMASS_D (  8, -2,   5) = ZERO
-      PWIDTH_D(  8, -2,   5) = ZERO
-      POW_D   (  8, -2,   5) = 2

=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_decl.inc	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_decl.inc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,6 @@
+      INTEGER MAXPROC_USED, MAXFLOW_USED
+      PARAMETER (MAXPROC_USED = 4)
+      PARAMETER (MAXFLOW_USED = 6)
+      INTEGER IDUP_D(8,5,MAXPROC_USED)
+      INTEGER MOTHUP_D(8,2,5,MAXPROC_USED)
+      INTEGER ICOLUP_D(8,2,5,MAXFLOW_USED)

=== removed file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_info.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_info.inc	2014-03-19 08:01:45 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%leshouche_info.inc	1970-01-01 00:00:00 +0000
@@ -1,156 +0,0 @@
-      INTEGER MAXPROC_USED, MAXFLOW_USED
-      PARAMETER (MAXPROC_USED = 4)
-      PARAMETER (MAXFLOW_USED = 6)
-      INTEGER IDUP_D(8,5,MAXPROC_USED)
-      INTEGER MOTHUP_D(8,2,5,MAXPROC_USED)
-      INTEGER ICOLUP_D(8,2,5,MAXFLOW_USED)
-      INTEGER ILH
-
-      DATA (IDUP_D(1,ILH,1),ILH=1,5)/21,21,6,-6,21/
-      DATA (MOTHUP_D(1,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(1,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(1,1,ILH,  1),ILH=1, 5)/501,502,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  1),ILH=1, 5)/502,503,  0,504,503/
-      DATA (ICOLUP_D(1,1,ILH,  2),ILH=1, 5)/501,504,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  2),ILH=1, 5)/502,503,  0,503,502/
-      DATA (ICOLUP_D(1,1,ILH,  3),ILH=1, 5)/503,501,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  3),ILH=1, 5)/502,503,  0,504,502/
-      DATA (ICOLUP_D(1,1,ILH,  4),ILH=1, 5)/504,501,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  4),ILH=1, 5)/502,503,  0,502,503/
-      DATA (ICOLUP_D(1,1,ILH,  5),ILH=1, 5)/504,502,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  5),ILH=1, 5)/502,503,  0,503,501/
-      DATA (ICOLUP_D(1,1,ILH,  6),ILH=1, 5)/503,504,501,  0,504/
-      DATA (ICOLUP_D(1,2,ILH,  6),ILH=1, 5)/502,503,  0,502,501/
-
-      DATA (IDUP_D(2,ILH,1),ILH=1,5)/21,21,6,-6,21/
-      DATA (MOTHUP_D(2,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(2,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(2,1,ILH,  1),ILH=1, 5)/501,502,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  1),ILH=1, 5)/502,503,  0,504,503/
-      DATA (ICOLUP_D(2,1,ILH,  2),ILH=1, 5)/501,504,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  2),ILH=1, 5)/502,503,  0,503,502/
-      DATA (ICOLUP_D(2,1,ILH,  3),ILH=1, 5)/503,501,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  3),ILH=1, 5)/502,503,  0,504,502/
-      DATA (ICOLUP_D(2,1,ILH,  4),ILH=1, 5)/504,501,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  4),ILH=1, 5)/502,503,  0,502,503/
-      DATA (ICOLUP_D(2,1,ILH,  5),ILH=1, 5)/504,502,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  5),ILH=1, 5)/502,503,  0,503,501/
-      DATA (ICOLUP_D(2,1,ILH,  6),ILH=1, 5)/503,504,501,  0,504/
-      DATA (ICOLUP_D(2,2,ILH,  6),ILH=1, 5)/502,503,  0,502,501/
-
-      DATA (IDUP_D(3,ILH,1),ILH=1,5)/21,21,6,-6,21/
-      DATA (MOTHUP_D(3,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(3,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(3,1,ILH,  1),ILH=1, 5)/501,502,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  1),ILH=1, 5)/502,503,  0,504,503/
-      DATA (ICOLUP_D(3,1,ILH,  2),ILH=1, 5)/501,504,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  2),ILH=1, 5)/502,503,  0,503,502/
-      DATA (ICOLUP_D(3,1,ILH,  3),ILH=1, 5)/503,501,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  3),ILH=1, 5)/502,503,  0,504,502/
-      DATA (ICOLUP_D(3,1,ILH,  4),ILH=1, 5)/504,501,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  4),ILH=1, 5)/502,503,  0,502,503/
-      DATA (ICOLUP_D(3,1,ILH,  5),ILH=1, 5)/504,502,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  5),ILH=1, 5)/502,503,  0,503,501/
-      DATA (ICOLUP_D(3,1,ILH,  6),ILH=1, 5)/503,504,501,  0,504/
-      DATA (ICOLUP_D(3,2,ILH,  6),ILH=1, 5)/502,503,  0,502,501/
-
-      DATA (IDUP_D(4,ILH,1),ILH=1,5)/21,21,6,-6,21/
-      DATA (MOTHUP_D(4,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(4,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(4,1,ILH,  1),ILH=1, 5)/501,502,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  1),ILH=1, 5)/502,503,  0,504,503/
-      DATA (ICOLUP_D(4,1,ILH,  2),ILH=1, 5)/501,504,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  2),ILH=1, 5)/502,503,  0,503,502/
-      DATA (ICOLUP_D(4,1,ILH,  3),ILH=1, 5)/503,501,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  3),ILH=1, 5)/502,503,  0,504,502/
-      DATA (ICOLUP_D(4,1,ILH,  4),ILH=1, 5)/504,501,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  4),ILH=1, 5)/502,503,  0,502,503/
-      DATA (ICOLUP_D(4,1,ILH,  5),ILH=1, 5)/504,502,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  5),ILH=1, 5)/502,503,  0,503,501/
-      DATA (ICOLUP_D(4,1,ILH,  6),ILH=1, 5)/503,504,501,  0,504/
-      DATA (ICOLUP_D(4,2,ILH,  6),ILH=1, 5)/502,503,  0,502,501/
-
-      DATA (IDUP_D(5,ILH,1),ILH=1,5)/-1,21,6,-6,-1/
-      DATA (MOTHUP_D(5,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(5,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(5,1,ILH,  1),ILH=1, 5)/  0,502,502,  0,  0/
-      DATA (ICOLUP_D(5,2,ILH,  1),ILH=1, 5)/501,503,  0,501,503/
-      DATA (ICOLUP_D(5,1,ILH,  2),ILH=1, 5)/  0,502,502,  0,  0/
-      DATA (ICOLUP_D(5,2,ILH,  2),ILH=1, 5)/501,503,  0,503,501/
-      DATA (ICOLUP_D(5,1,ILH,  3),ILH=1, 5)/  0,501,502,  0,  0/
-      DATA (ICOLUP_D(5,2,ILH,  3),ILH=1, 5)/501,503,  0,503,502/
-      DATA (ICOLUP_D(5,1,ILH,  4),ILH=1, 5)/  0,501,502,  0,  0/
-      DATA (ICOLUP_D(5,2,ILH,  4),ILH=1, 5)/501,503,  0,502,503/
-      DATA (IDUP_D(5,ILH,2),ILH=1,5)/-3,21,6,-6,-3/
-      DATA (MOTHUP_D(5,1,ILH,  2),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(5,2,ILH,  2),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(5,ILH,3),ILH=1,5)/-2,21,6,-6,-2/
-      DATA (MOTHUP_D(5,1,ILH,  3),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(5,2,ILH,  3),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(5,ILH,4),ILH=1,5)/-4,21,6,-6,-4/
-      DATA (MOTHUP_D(5,1,ILH,  4),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(5,2,ILH,  4),ILH=1, 5)/  0,  0,  2,  2,  2/
-
-      DATA (IDUP_D(6,ILH,1),ILH=1,5)/1,21,6,-6,1/
-      DATA (MOTHUP_D(6,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(6,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(6,1,ILH,  1),ILH=1, 5)/503,501,501,  0,502/
-      DATA (ICOLUP_D(6,2,ILH,  1),ILH=1, 5)/  0,503,  0,502,  0/
-      DATA (ICOLUP_D(6,1,ILH,  2),ILH=1, 5)/502,501,501,  0,502/
-      DATA (ICOLUP_D(6,2,ILH,  2),ILH=1, 5)/  0,503,  0,503,  0/
-      DATA (ICOLUP_D(6,1,ILH,  3),ILH=1, 5)/503,502,501,  0,502/
-      DATA (ICOLUP_D(6,2,ILH,  3),ILH=1, 5)/  0,503,  0,501,  0/
-      DATA (ICOLUP_D(6,1,ILH,  4),ILH=1, 5)/501,502,501,  0,502/
-      DATA (ICOLUP_D(6,2,ILH,  4),ILH=1, 5)/  0,503,  0,503,  0/
-      DATA (IDUP_D(6,ILH,2),ILH=1,5)/3,21,6,-6,3/
-      DATA (MOTHUP_D(6,1,ILH,  2),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(6,2,ILH,  2),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(6,ILH,3),ILH=1,5)/2,21,6,-6,2/
-      DATA (MOTHUP_D(6,1,ILH,  3),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(6,2,ILH,  3),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(6,ILH,4),ILH=1,5)/4,21,6,-6,4/
-      DATA (MOTHUP_D(6,1,ILH,  4),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(6,2,ILH,  4),ILH=1, 5)/  0,  0,  2,  2,  2/
-
-      DATA (IDUP_D(7,ILH,1),ILH=1,5)/21,-1,6,-6,-1/
-      DATA (MOTHUP_D(7,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(7,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(7,1,ILH,  1),ILH=1, 5)/501,  0,502,  0,  0/
-      DATA (ICOLUP_D(7,2,ILH,  1),ILH=1, 5)/503,501,  0,503,502/
-      DATA (ICOLUP_D(7,1,ILH,  2),ILH=1, 5)/501,  0,502,  0,  0/
-      DATA (ICOLUP_D(7,2,ILH,  2),ILH=1, 5)/503,501,  0,502,503/
-      DATA (ICOLUP_D(7,1,ILH,  3),ILH=1, 5)/502,  0,502,  0,  0/
-      DATA (ICOLUP_D(7,2,ILH,  3),ILH=1, 5)/503,501,  0,503,501/
-      DATA (ICOLUP_D(7,1,ILH,  4),ILH=1, 5)/502,  0,502,  0,  0/
-      DATA (ICOLUP_D(7,2,ILH,  4),ILH=1, 5)/503,501,  0,501,503/
-      DATA (IDUP_D(7,ILH,2),ILH=1,5)/21,-3,6,-6,-3/
-      DATA (MOTHUP_D(7,1,ILH,  2),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(7,2,ILH,  2),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(7,ILH,3),ILH=1,5)/21,-2,6,-6,-2/
-      DATA (MOTHUP_D(7,1,ILH,  3),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(7,2,ILH,  3),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(7,ILH,4),ILH=1,5)/21,-4,6,-6,-4/
-      DATA (MOTHUP_D(7,1,ILH,  4),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(7,2,ILH,  4),ILH=1, 5)/  0,  0,  2,  2,  2/
-
-      DATA (IDUP_D(8,ILH,1),ILH=1,5)/21,1,6,-6,1/
-      DATA (MOTHUP_D(8,1,ILH,  1),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(8,2,ILH,  1),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (ICOLUP_D(8,1,ILH,  1),ILH=1, 5)/501,503,501,  0,502/
-      DATA (ICOLUP_D(8,2,ILH,  1),ILH=1, 5)/503,  0,  0,502,  0/
-      DATA (ICOLUP_D(8,1,ILH,  2),ILH=1, 5)/501,502,501,  0,502/
-      DATA (ICOLUP_D(8,2,ILH,  2),ILH=1, 5)/503,  0,  0,503,  0/
-      DATA (ICOLUP_D(8,1,ILH,  3),ILH=1, 5)/502,503,501,  0,502/
-      DATA (ICOLUP_D(8,2,ILH,  3),ILH=1, 5)/503,  0,  0,501,  0/
-      DATA (ICOLUP_D(8,1,ILH,  4),ILH=1, 5)/502,501,501,  0,502/
-      DATA (ICOLUP_D(8,2,ILH,  4),ILH=1, 5)/503,  0,  0,503,  0/
-      DATA (IDUP_D(8,ILH,2),ILH=1,5)/21,3,6,-6,3/
-      DATA (MOTHUP_D(8,1,ILH,  2),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(8,2,ILH,  2),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(8,ILH,3),ILH=1,5)/21,2,6,-6,2/
-      DATA (MOTHUP_D(8,1,ILH,  3),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(8,2,ILH,  3),ILH=1, 5)/  0,  0,  2,  2,  2/
-      DATA (IDUP_D(8,ILH,4),ILH=1,5)/21,4,6,-6,4/
-      DATA (MOTHUP_D(8,1,ILH,  4),ILH=1, 5)/  0,  0,  1,  1,  1/
-      DATA (MOTHUP_D(8,2,ILH,  4),ILH=1, 5)/  0,  0,  2,  2,  2/
-

=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc	2014-11-06 06:25:51 +0000
@@ -0,0 +1,12 @@
+      INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
+      PARAMETER (LMAXCONFIGS_USED=   5)
+      PARAMETER (MAX_BRANCH_USED =   3)
+      INTEGER MAPCONFIG_D(  6,0:LMAXCONFIGS_USED)
+      INTEGER IFOREST_D(  6,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER SPROP_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER TPRID_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PMASS_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED
+     $ )
+      DOUBLE PRECISION PWIDTH_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USE
+     $ D)
+      INTEGER POW_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)

=== removed file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_info.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_info.inc	2014-03-19 08:01:45 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%configs_and_props_info.inc	1970-01-01 00:00:00 +0000
@@ -1,419 +0,0 @@
-      INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
-      PARAMETER (LMAXCONFIGS_USED=   5)
-      PARAMETER (MAX_BRANCH_USED =   3)
-      INTEGER MAPCONFIG_D(  6,0:LMAXCONFIGS_USED)
-      INTEGER IFOREST_D(  6,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER SPROP_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER TPRID_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PMASS_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED
-     $ )
-      DOUBLE PRECISION PWIDTH_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USE
-     $ D)
-      INTEGER POW_D(  6,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-C     
-C     nFKSprocess 1
-C     Diagram 1 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   1)/   1/
-      DATA (IFOREST_D(  1, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,   1)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  1,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   2)/   2/
-      DATA (IFOREST_D(  1, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  1,  -1,   2)/       6/
-      DATA (IFOREST_D(  1, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  1,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   3)/   3/
-      DATA (IFOREST_D(  1, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  1,  -1,   3)/      -6/
-      DATA (IFOREST_D(  1, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  1,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   4)/   4/
-      DATA (IFOREST_D(  1, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,   4)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  1,  -2,   4)/       2/
-      DATA (IFOREST_D(  1, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,   5)/   5/
-      DATA (IFOREST_D(  1, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  1,  -1,   5)/      21/
-      DATA (IFOREST_D(  1, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  1,  -2,   5)/       2/
-      DATA (IFOREST_D(  1, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 1
-      DATA MAPCONFIG_D(  1,0)/   5/
-C     
-C     nFKSprocess 2
-C     Diagram 1 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   1)/   1/
-      DATA (IFOREST_D(  2, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,   1)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  2,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   2)/   2/
-      DATA (IFOREST_D(  2, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  2,  -1,   2)/       6/
-      DATA (IFOREST_D(  2, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  2,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   3)/   3/
-      DATA (IFOREST_D(  2, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  2,  -1,   3)/      -6/
-      DATA (IFOREST_D(  2, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  2,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   4)/   4/
-      DATA (IFOREST_D(  2, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,   4)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  2,  -2,   4)/       2/
-      DATA (IFOREST_D(  2, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,   5)/   5/
-      DATA (IFOREST_D(  2, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  2,  -1,   5)/      21/
-      DATA (IFOREST_D(  2, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  2,  -2,   5)/       2/
-      DATA (IFOREST_D(  2, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 2
-      DATA MAPCONFIG_D(  2,0)/   5/
-C     
-C     nFKSprocess 3
-C     Diagram 1 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   1)/   1/
-      DATA (IFOREST_D(  3, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,   1)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  3,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   2)/   2/
-      DATA (IFOREST_D(  3, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  3,  -1,   2)/       6/
-      DATA (IFOREST_D(  3, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  3,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   3)/   3/
-      DATA (IFOREST_D(  3, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  3,  -1,   3)/      -6/
-      DATA (IFOREST_D(  3, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  3,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   4)/   4/
-      DATA (IFOREST_D(  3, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,   4)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  3,  -2,   4)/       2/
-      DATA (IFOREST_D(  3, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,   5)/   5/
-      DATA (IFOREST_D(  3, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  3,  -1,   5)/      21/
-      DATA (IFOREST_D(  3, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  3,  -2,   5)/       2/
-      DATA (IFOREST_D(  3, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 3
-      DATA MAPCONFIG_D(  3,0)/   5/
-C     
-C     nFKSprocess 4
-C     Diagram 1 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   1)/   1/
-      DATA (IFOREST_D(  4, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,   1)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  4,  -2,   1)/      21/
-C     Diagram 2 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   2)/   2/
-      DATA (IFOREST_D(  4, IFR, -1,   2),IFR=1,2)/  5,  3/
-      DATA SPROP_D(  4,  -1,   2)/       6/
-      DATA (IFOREST_D(  4, IFR, -2,   2),IFR=1,2)/  4, -1/
-      DATA SPROP_D(  4,  -2,   2)/      21/
-C     Diagram 3 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   3)/   3/
-      DATA (IFOREST_D(  4, IFR, -1,   3),IFR=1,2)/  5,  4/
-      DATA SPROP_D(  4,  -1,   3)/      -6/
-      DATA (IFOREST_D(  4, IFR, -2,   3),IFR=1,2)/ -1,  3/
-      DATA SPROP_D(  4,  -2,   3)/      21/
-C     Diagram 4 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   4)/   4/
-      DATA (IFOREST_D(  4, IFR, -1,   4),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,   4)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,   4),IFR=1,2)/  1,  5/
-      DATA TPRID_D(  4,  -2,   4)/       2/
-      DATA (IFOREST_D(  4, IFR, -3,   4),IFR=1,2)/ -2, -1/
-C     Diagram 5 for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,   5)/   5/
-      DATA (IFOREST_D(  4, IFR, -1,   5),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  4,  -1,   5)/      21/
-      DATA (IFOREST_D(  4, IFR, -2,   5),IFR=1,2)/  1, -1/
-      DATA TPRID_D(  4,  -2,   5)/       2/
-      DATA (IFOREST_D(  4, IFR, -3,   5),IFR=1,2)/ -2,  5/
-C     Number of configs for nFKSprocess 4
-      DATA MAPCONFIG_D(  4,0)/   5/
-C     
-C     nFKSprocess 5
-C     Diagram 1 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   1)/   1/
-      DATA (IFOREST_D(  5, IFR, -1,   1),IFR=1,2)/  4,  3/
-      DATA SPROP_D(  5,  -1,   1)/      21/
-      DATA (IFOREST_D(  5, IFR, -2,   1),IFR=1,2)/  5, -1/
-      DATA SPROP_D(  5,  -2,   1)/      -2/
-C     Diagram 2 for nFKSprocess 5
-      DATA MAPCONFIG_D(  5,   2)/   2/
-      DATA (IFOREST_D(  5, IFR, -1,   2),IFR=1,2)/  1,  3/
-      DATA TPRID_D(  5,  -1,   2)/       6/
-      DATA (IFOREST_D(  5, IFR, -2,   2),IFR=1,2)/ -1,  4/
-      DATA TPRID_D(  5,  -2,   2)/      21/