Merge lp:~maddevelopers/mg5amcnlo/2.5.5_NLOandLOmerging into lp:~mg5core2/mg5amcnlo/2.6.1

[marco zaro]

The LOonly output mode has been modified in order to have the diagram and color informations needed to do FxFx merging including LO matrix elements after a given multiplicity.
Josh Bendavid suggested to make this functionality available in the next release even some more tests need to be done on the results.

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

That's cool, before starting the review:

So What would be the syntax in this case? (Do you have this documented somewhere?
Would it be:
generate p p > e+ e- [QCD]
add process p p > e+ e- j [LOonly]
add process p p > e+ e- j j [LOonly]
Or should I skip the one jet multiplicity?

Do I understand correctly that this is not working for UNLOPS?

Can you be more specific by "it required more test"? Should we release this functionality as "beta"?

Cheers,

Olivier

[Rikkert Frederix (frederix)]

Hi Olivier,

Yes, that's indeed the syntax.

All tests with kind of low statistics have been performed. However, we wanted to give the code to Josh Bendavid (who requested this feature) to perform some high-statistics tests and to confirm that's working within the CMS framework. He prefers that the code is part of the main branch, as he's loosing track of all separate contributions. That's why Marco wrote that some more testing is needed. Hence, we should keep this functionality somewhat hidden (or in beta) for the moment.

Indeed, it's not working with UNLOPS. However, with UNLOPS one can always generate separate samples and put them together at the end (at the cost of introducing large fractions of negative weights).

best,
Rikkert

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

I'm starting to look at this. My question are so far on the interface/physics limitation.
Actually a lot of the point below should also be consider for FxFx run in general since they are not specific to LOonly.

1) the syntax "add process p p > e+ e- j [LOonly]" does not work. Would be great to not be force to do "add process p p > e+ e- j [LOonly=QCD]" since the QCD here is meaningless (I guess)

2) It would be great that ickkw being on default on 3 in this case.

3) Would be great that the code crash if it is on 4.

4) What would be the meaning of having order=LO in this setup? (This is actually a valid question even without the LOonly option). I guess the answer is None. Should we prevent it?
(actually same question with fixed_order=ON)

5) Does it make sense to have the shower ON by default? should not we have turn it off?

6) If this is to be consider as beta, I would prefer that a warning appear when the matrix-element are generated. (Such that this is clear)

7) When running, I have the following:
   --------------------------------------------------------------
      Summary:
      Process p p > e+ e- [QCD] ; p p > e+ e- j [LOonly=QCD]
      Run at p-p collider (6500.0 + 6500.0 GeV)
      Number of events generated: 10000
      Total cross section: 2.103e+03 +- 8.4e+00 pb
   --------------------------------------------------------------
      Scale variation (computed from LHE events):
          Dynamical_scale_choice -1 (envelope of 9 values):
              2.130e+03 pb +6.3% -11.9%
   --------------------------------------------------------------

Can we trust the scale variation?

I did not test the physics but I checked the code and this sounds fine.
I would have prefer that you factorise the code to have less a huge function with a bunch of if.
But this is not critical and we will be able to split the function in smaller piece if needed in the future.

Thanks for this work,

Olivier

[Rikkert Frederix (frederix)]

Hi Olivier,

I see your points. However, most are very generally to FxFx and do not really apply specifically to this branch. Therefore, I'm not sure that this is the right time to implement them...

>
> 1) the syntax "add process p p > e+ e- j [LOonly]" does not work. Would be
> great to not be force to do "add process p p > e+ e- j [LOonly=QCD]" since the
> QCD here is meaningless (I guess)

I think this would be a good idea. Marco, do you know if this is easy? Will at then also work with add-ons like MadSpin?

> 2) It would be great that ickkw being on default on 3 in this case.
>

Boahh.. don't know. Is this relevant? And how do you know when to put this? Do you suggest to count the 'j' in each 'add process' and check if it's more than in the 'generate' command? Is this sufficient?

> 3) Would be great that the code crash if it is on 4.
>

Same as for previous comment.

> 4) What would be the meaning of having order=LO in this setup? (This is
> actually a valid question even without the LOonly option). I guess the answer
> is None. Should we prevent it?
> (actually same question with fixed_order=ON)
>

The results for order=LO are not necessarily wrong. Therefore I would keep it. Same for fixed_order=ON.

> 5) Does it make sense to have the shower ON by default? should not we have
> turn it off?
>

It should be "On" in my opinion. But also here: in practice this will only be used by experimentalists. They know how to deal with either On or Off.

>
> 6) If this is to be consider as beta, I would prefer that a warning appear
> when the matrix-element are generated. (Such that this is clear)
>

The [LOonly=QCD] option is not in beta. Rather the NLO+LO merging is. Hence at the level of the generation, one doesn't know yet if the user will use the merging or not.

>
> 7) When running, I have the following:
> --------------------------------------------------------------
> Summary:
> Process p p > e+ e- [QCD] ; p p > e+ e- j [LOonly=QCD]
> Run at p-p collider (6500.0 + 6500.0 GeV)
> Number of events generated: 10000
> Total cross section: 2.103e+03 +- 8.4e+00 pb
> --------------------------------------------------------------
> Scale variation (computed from LHE events):
> Dynamical_scale_choice -1 (envelope of 9 values):
> 2.130e+03 pb +6.3% -11.9%
> --------------------------------------------------------------
>
> Can we trust the scale variation?

No. Nor the cross section as the showering will/should reject a fraction of the events. Just as for NLO FxFx. However, I'm against removing this since it serves as a check that all required weights are written correctly in the event-file.

Best,
Rikkert

[marco zaro (marco-zaro)]

Hi,
thanks to Olivier for the review and to Rik for the answers
On 02 Oct 2017, at 15:38, Rikkert Frederix <email address hidden> wrote:

> Hi Olivier,
>
> I see your points. However, most are very generally to FxFx and do not really apply specifically to this branch. Therefore, I'm not sure that this is the right time to implement them...
>
>>
>> 1) the syntax "add process p p > e+ e- j [LOonly]" does not work. Would be
>> great to not be force to do "add process p p > e+ e- j [LOonly=QCD]" since the
>> QCD here is meaningless (I guess)
> I think this would be a good idea. Marco, do you know if this is easy? Will at then also work with add-ons like MadSpin?

hmm…
currently the code searches a regular expression like [<option> = <orders>]. If option is not there (like for [QCD]) the option is set to ‘all’. In our case, [LOonly] would be equivalent to ‘all=LOonly’, which is certainly something that does not make sense. So, i don’t have any smarter idea than checking that, when one has orders=LOonly and option =all, then one assigns orders=QCD, option=LOonly with an ‘if’ around line 177 inside master_interface.
It is rather ugly, but I have no smarter idea…
What would you say?

Cheers,

Marco

>
>> 2) It would be great that ickkw being on default on 3 in this case.
>>
>
> Boahh.. don't know. Is this relevant? And how do you know when to put this? Do you suggest to count the 'j' in each 'add process' and check if it's more than in the 'generate' command? Is this sufficient?
>
>
>> 3) Would be great that the code crash if it is on 4.
>>
>
> Same as for previous comment.
>
>> 4) What would be the meaning of having order=LO in this setup? (This is
>> actually a valid question even without the LOonly option). I guess the answer
>> is None. Should we prevent it?
>> (actually same question with fixed_order=ON)
>>
>
> The results for order=LO are not necessarily wrong. Therefore I would keep it. Same for fixed_order=ON.
>
>> 5) Does it make sense to have the shower ON by default? should not we have
>> turn it off?
>>
>
> It should be "On" in my opinion. But also here: in practice this will only be used by experimentalists. They know how to deal with either On or Off.
>
>>
>> 6) If this is to be consider as beta, I would prefer that a warning appear
>> when the matrix-element are generated. (Such that this is clear)
>>
>
> The [LOonly=QCD] option is not in beta. Rather the NLO+LO merging is. Hence at the level of the generation, one doesn't know yet if the user will use the merging or not.
>
>>
>> 7) When running, I have the following:
>> --------------------------------------------------------------
>> Summary:
>> Process p p > e+ e- [QCD] ; p p > e+ e- j [LOonly=QCD]
>> Run at p-p collider (6500.0 + 6500.0 GeV)
>> Number of events generated: 10000
>> Total cross section: 2.103e+03 +- 8.4e+00 pb
>> --------------------------------------------------------------
>> Scale variation (computed from LHE events):
>> Dynamical_scale_choice -1 (envelope of 9 values):
>> 2.130e+03 pb +6.3% -11.9%
>> --------------------------------------------------------------
>>
>> Can we trus...

[Stefano Frixione (stefano-frixione)]

Hi Olivier,
a few scattered comments. In general, this branch has been set up to
comply with CMS request of having a LO option on top of FxFx, in order
to extend somewhat the multiplicity reach. So I'm not sure a discussion
unrelated to FxFx does make sense (although technically it does).

> 1) the syntax "add process p p > e+ e- j [LOonly]" does not work. Would
> be great to not be force to do "add process p p > e+ e- j
> [LOonly=QCD]" since the QCD here is meaningless (I guess)
not necessarily: there could be processes mediated by either QCD
or QED, with the same final state.
Since this is going to be used by advanced users, I think the effort
spent into rendering this cleaner is not particularly worth.

> 2) It would be great that ickkw being on default on 3 in this case.
Same as before.

> 3) Would be great that the code crash if it is on 4.
Same as before.

> 5) Does it make sense to have the shower ON by default? should not we
> have turn it off?
I'm inclined to have it off, since this prevents the code from crashing
in absence of MC installations (or if they have been done incorrectly).

> 7) When running, I have the following:
> --------------------------------------------------------------
> Summary:
> Process p p > e+ e- [QCD] ; p p > e+ e- j [LOonly=QCD]
> Run at p-p collider (6500.0 + 6500.0 GeV)
> Number of events generated: 10000
> Total cross section: 2.103e+03 +- 8.4e+00 pb
> --------------------------------------------------------------
> Scale variation (computed from LHE events):
> Dynamical_scale_choice -1 (envelope of 9 values):
> 2.130e+03 pb +6.3% -11.9%
> --------------------------------------------------------------
>
> Can we trust the scale variation?
I agree with Rikkert: maybe not, but it's a double check.

Cheers, Stefano.

>
> I did not test the physics but I checked the code and this sounds fine.
> I would have prefer that you factorise the code to have less a huge function with a bunch of if.
> But this is not critical and we will be able to split the function in smaller piece if needed in the future.
>
> Thanks for this work,
>
> Olivier
>
> --
> https://code.launchpad.net/~maddevelopers/mg5amcnlo/2.5.5_NLOandLOmerging/+merge/331151
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/mg5amcnlo/2.5.5_NLOandLOmerging.
>

[marco zaro (marco-zaro)]

Hi all,
any news on this?

thanks!
cheers,

Marco

[Rikkert Frederix (frederix)]

review: Approve

[lp:~maddevelopers/mg5amcnlo/2.5.5_NLOandLOmerging]

303. By Olivier Mattelaer on 2017-11-26

merge with 2.6.1

304. By Olivier Mattelaer on 2017-11-26

merge with latest 2.6.1

305. By Olivier Mattelaer on 2017-11-30

allow [LOonly] syntax

[Olivier Mattelaer (olivier-mattelaer)]

I'm doing the merge

review: Approve

=== modified file 'Template/NLO/Source/setrun.f'
--- Template/NLO/Source/setrun.f	2017-07-25 14:13:14 +0000
+++ Template/NLO/Source/setrun.f	2017-11-30 23:06:43 +0000
@@ -125,7 +125,7 @@
          else
             idbmup(i)=lpp(i)
          endif
-        ebmup(i)=ebeam(i)
+         ebmup(i)=ebeam(i)
       enddo
       call get_pdfup(pdlabel,pdfgup,pdfsup,lhaid)
 c Fill the nmemPDF(i) array with the number of PDF error set. This we
 
=== modified file 'Template/NLO/SubProcesses/driver_mintFO.f'
--- Template/NLO/SubProcesses/driver_mintFO.f	2017-06-19 17:15:00 +0000
+++ Template/NLO/SubProcesses/driver_mintFO.f	2017-11-30 23:06:43 +0000
@@ -784,11 +784,11 @@
             write (*,*) 'Process generated with [LOonly=QCD]. '/
      $           /'Setting abrv to "born".'
             abrv='born'
-            if (ickkw.eq.3) then
-               write (*,*) 'FxFx merging not possible with'/
-     $              /' [LOonly=QCD] processes'
-               stop 1
-            endif
+c$$$            if (ickkw.eq.3) then
+c$$$               write (*,*) 'FxFx merging not possible with'/
+c$$$     $              /' [LOonly=QCD] processes'
+c$$$               stop 1
+c$$$            endif
          endif
       endif
 c
 
=== modified file 'Template/NLO/SubProcesses/driver_mintMC.f'
--- Template/NLO/SubProcesses/driver_mintMC.f	2017-07-29 08:03:34 +0000
+++ Template/NLO/SubProcesses/driver_mintMC.f	2017-11-30 23:06:43 +0000
@@ -711,11 +711,11 @@
             write (*,*) 'Process generated with [LOonly=QCD]. '/
      $           /'Setting abrv to "born".'
             abrv='born'
-            if (ickkw.eq.3) then
-               write (*,*) 'FxFx merging not possible with'/
-     $              /' [LOonly=QCD] processes'
-               stop 1
-            endif
+c$$$            if (ickkw.eq.3) then
+c$$$               write (*,*) 'FxFx merging not possible with'/
+c$$$     $              /' [LOonly=QCD] processes'
+c$$$               stop 1
+c$$$            endif
          endif
       endif
       if(nbody.and.abrv.ne.'born'.and.abrv(1:2).ne.'vi'
 
=== modified file 'Template/NLO/SubProcesses/fks_singular.f'
--- Template/NLO/SubProcesses/fks_singular.f	2017-07-26 11:45:39 +0000
+++ Template/NLO/SubProcesses/fks_singular.f	2017-11-30 23:06:43 +0000
@@ -54,6 +54,8 @@
       common /cxiScut_used/xiScut_used,xiBSVcut_used
       double precision fxfx_exp_rewgt
       common /c_fxfx_exp_regt/ fxfx_exp_rewgt
+      character*4      abrv
+      common /to_abrv/ abrv
       call cpu_time(tBefore)
       if (f_nb.eq.0d0) return
       if (xi_i_fks_ev .gt. xiBSVcut_used) return
@@ -62,7 +64,8 @@
       g22=g**(nint(2*wgtbpower+2))
       wgt1=wgtnstmp*f_nb/g22
       wgt4=wgtnstmp_avgvirt*f_nb/g22
-      if (ickkw.eq.3 .and. fxfx_exp_rewgt.ne.0d0) then
+      if (ickkw.eq.3 .and. fxfx_exp_rewgt.ne.0d0 .and. abrv.ne.'born')
+     $     then
          wgt1=wgt1 - fxfx_exp_rewgt*born_wgt*f_nb/g2/(4d0*pi)
       elseif (ickkw.eq.-1) then
          if (wgtbpower.ne.0) then
 
=== modified file 'Template/NLO/SubProcesses/reweight.f'
--- Template/NLO/SubProcesses/reweight.f	2017-06-19 14:42:53 +0000
+++ Template/NLO/SubProcesses/reweight.f	2017-11-30 23:06:43 +0000
@@ -48,7 +48,7 @@
          stop
       endif
 
-      if (iipdg.eq.21) then
+      if (abs(iipdg).eq.21) then
 c Gluon sudakov         
 c     g->gg contribution
          gamma=CA*alphasq0*log(Q1/q0)/pi                     ! A1
@@ -211,7 +211,9 @@
       integer ipdg, irfl
       integer get_color
 
-      isqcd=(iabs(get_color(ipdg)).gt.1)
+      irfl=ipdg
+      if (irfl.eq.-21)irfl=21
+      isqcd=(iabs(get_color(irfl)).gt.1)
 
       return
       end
@@ -276,6 +278,9 @@
       idmo=ipdg(imo)
       idda1=ipdg(ida1)
       idda2=ipdg(ida2)
+      if (idmo.eq.-21) idmo=21
+      if (idda1.eq.-21) idda1=21
+      if (idda2.eq.-21) idda2=21
 
       if (btest(mlevel,4)) then
         write(*,*) ' updating ipart for: ',ida1,ida2,' -> ',imo
@@ -429,6 +434,7 @@
       idmo=ipdg(imo)
       idda1=ipdg(ida1)
       idda2=ipdg(ida2)
+
 c     Check QCD vertex
       if(islast.or..not.isqcd(idmo).or..not.isqcd(idda1).or.
      &     .not.isqcd(idda2)) then
@@ -795,7 +801,7 @@
 c     line
                   if(goodjet(ida(3-i))) then
                      if(partonline(j).or.
-     $                    ipdgcl(ida(3-i),igraphs(1),nFKSprocess).eq.21)then
+     $                    abs(ipdgcl(ida(3-i),igraphs(1),nFKSprocess)).eq.21)then
 c     Need to include gluon to avoid soft singularity
                         iqjets(ipart(1,ida(3-i)))=1 ! 1 means for sure jet
                      else
 
=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2017-11-24 22:52:27 +0000
+++ UpdateNotes.txt	2017-11-30 23:06:43 +0000
@@ -1,4 +1,9 @@
 Update notes for MadGraph5_aMC@NLO (in reverse time order)
+2.6.1 (XX/XX/17)
+      New Functionalities:
+      --------------------
+      RF+MZ: It is now possible to add LO matrix elements (with [LOonly]) 
+             to Fx-Fx merged samples. Thanks to Josh Bendavid for testing.
 
 2.6.1
       OM: Re-factoring the code asking which program to use (both at LO and NLO)
 
=== modified file 'madgraph/fks/fks_helas_objects.py'
--- madgraph/fks/fks_helas_objects.py	2016-02-23 09:08:23 +0000
+++ madgraph/fks/fks_helas_objects.py	2017-11-30 23:06:43 +0000
@@ -125,6 +125,8 @@
                                 gen_color=False)
 
     processes = helasfull.born_matrix_element.get('processes')
+
+    max_configs = helasfull.born_matrix_element.get_num_configs()
     
     metag = helas_objects.IdentifyMETag.create_tag(helasfull.born_matrix_element.get('base_amplitude'))
     
@@ -134,7 +136,7 @@
     cPickle.dump(outdata,output,protocol=2)
     output.close()
     
-    return [output.name,metag,has_loops,processes]
+    return [output.name,metag,has_loops,processes,max_configs]
 
 
 def async_finalize_matrix_elements(args):
@@ -322,6 +324,8 @@
                 self['has_loops'] = self['has_loops'] or has_loops
                 processes = bornout[3]
                 self['processes'].extend(processes)
+                self['max_particles'] = max([self['max_configs']] + [len(p['legs']) + 1 for p in bornout[3]])
+                self['max_configs'] = max(self['max_configs'], bornout[4])
                 unique = True
                 for ime2,bornout2 in enumerate(unique_me_list):
                     mefile2 = bornout2[0]
@@ -394,12 +398,12 @@
                     break
             self['has_virtuals'] = has_virtuals
             
-            configs_list = []
+            configs_list = [self['max_configs']]
             for meout in realmapout:
                 configs_list.append(meout[1])
             self['max_configs'] = max(configs_list)
             
-            nparticles_list = []
+            nparticles_list = [self['max_particles']]
             for meout in realmapout:
                 nparticles_list.append(meout[2])
             self['max_particles'] = max(nparticles_list)        
@@ -453,14 +457,16 @@
 
     def get_max_configs(self):
         """Return max_configs"""
-
+            
         if self['max_configs'] < 0:
             try:
                 self['max_configs'] = max([me.get_num_configs() \
                                   for me in self['real_matrix_elements']])
             except ValueError:
-                self['max_configs'] = max([me.born_matrix_element.get_num_configs() \
-                                  for me in self['matrix_elements']])
+                pass
+            self['max_configs'] = max(self['max_configs'],\
+                                      max([me.born_matrix_element.get_num_configs() \
+                                           for me in self['matrix_elements']]))
 
         return self['max_configs']
 
 
=== modified file 'madgraph/interface/madgraph_interface.py'
--- madgraph/interface/madgraph_interface.py	2017-11-23 12:05:44 +0000
+++ madgraph/interface/madgraph_interface.py	2017-11-30 23:06:43 +0000
@@ -4535,7 +4535,7 @@
                 raise self.InvalidCmd, "No particle %s in model" % part_name
 
         # Apply the keyword 'all' for perturbed coupling orders.
-        if perturbation_couplings.lower()=='all':
+        if perturbation_couplings.lower() in ['all', 'loonly']:
             perturbation_couplings=' '.join(self._curr_model['perturbation_couplings'])
 
         if filter(lambda leg: leg.get('state') == True, myleglist):
 
=== modified file 'madgraph/interface/master_interface.py'
--- madgraph/interface/master_interface.py	2017-03-27 08:37:47 +0000
+++ madgraph/interface/master_interface.py	2017-11-30 23:06:43 +0000
@@ -184,8 +184,11 @@
             else:
                 # If not option is set the convention is that the mode is 'all'
                 # unless no perturbation orders is defined.
-                if len(orders)>0:
-                    return ('NLO','all',orders)
+                # if order is set to LOonly assume LOonly=QCD
+                if orders == ['LOonly']:
+                    return ('NLO', 'LOonly', ['QCD'])
+                elif len(orders)>0:
+                    return ('NLO','all',orders) 
                 else:
                     return ('tree',None,[])               
         else:
 
=== modified file 'madgraph/iolibs/export_fks.py'
--- madgraph/iolibs/export_fks.py	2017-06-16 15:25:01 +0000
+++ madgraph/iolibs/export_fks.py	2017-11-30 23:06:43 +0000
@@ -524,7 +524,7 @@
                          ncolor,maxflow,fortran_model)
 
         filename = 'configs_and_props_info.dat'
-        nconfigs,max_leg_number,nfksconfs=self.write_configs_and_props_info_file(
+        nconfigs,max_leg_number=self.write_configs_and_props_info_file(
                               filename, 
                               matrix_element)
 
@@ -670,12 +670,6 @@
         else:
             os.system("ln -s ../BinothLHA_user.f ./BinothLHA.f")
 
-
-        #import nexternal/leshouches in Source
-#        ln('nexternal.inc', '../../Source', log=False)
-#        ln('born_leshouche.inc', '../../Source', log=False)
-
-
         # Return to SubProcesses dir
         os.chdir(os.path.pardir)
         # Add subprocess to subproc.mg
@@ -683,7 +677,6 @@
         files.append_to_file(filename,
                              self.write_subproc,
                              borndir)
-
             
         os.chdir(cwd)
         # Generate info page
@@ -955,9 +948,28 @@
             lines.append("data (real_from_born_conf(irfbc,%d),irfbc=1,%d) /%s/" \
                              % (iFKS,len(links),real_configs))
 
+        # this is for 'LOonly' processes; in this case, a fake configuration 
+        # with all the born diagrams is written
+        if not matrix_element.get_fks_info_list():
+            # compute (again) the number of configurations at the born
+            base_diagrams = born_me.get('base_amplitude').get('diagrams')
+            minvert = min([max([len(vert.get('legs')) for vert in \
+                                    diag.get('vertices')]) for diag in base_diagrams])
+    
+            for idiag, diag in enumerate(base_diagrams):
+                if any([len(vert.get('legs')) > minvert for vert in
+                        diag.get('vertices')]):
+                # Only 3-vertices allowed in configs.inc
+                    continue
+                max_links = max_links + 1
+                
+            real_configs=', '.join(['%d' % i for i in range(1, max_links+1)])
+            lines.append("data (real_from_born_conf(irfbc,%d),irfbc=1,%d) /%s/" \
+                             % (1,max_links,real_configs))
+
         lines2.append("integer irfbc")
         lines2.append("integer real_from_born_conf(%d,%d)" \
-                         % (max_links,len(matrix_element.get_fks_info_list())))
+                         % (max_links, max(len(matrix_element.get_fks_info_list()),1)))
         # Write the file
         writer.writelines(lines2+lines)
 
@@ -1057,12 +1069,14 @@
         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.get_nexternal_ninitial()
 
         max_iconfig=0
         max_leg_number=0
 
+        ########################################################
+        # this is for standard processes with [(real=)XXX]
+        ########################################################
         for iFKS, conf in enumerate(matrix_element.get_fks_info_list()):
             iFKS=iFKS+1
             iconfig = 0
@@ -1147,11 +1161,141 @@
                                      (iFKS,leg.get('number'), iconf + 1, leg.get('id')))
                     lines2.append("P   %4d   %4d   %4d   %4d " % \
                                      (iFKS,leg.get('number'), iconf + 1, pow_part))
+
+        ########################################################
+        # this is for [LOonly=XXX]
+        ########################################################
+        if not matrix_element.get_fks_info_list():
+            born_me = matrix_element.born_matrix_element
+            # as usual, in this case we assume just one FKS configuration 
+            # exists with diagrams corresponding to born ones X the ij -> i,j
+            # splitting. Here j is chosen to be the last colored particle in
+            # the particle list
+            bornproc = born_me.get('processes')[0]
+            colors = [l.get('color') for l in bornproc.get('legs')] 
+
+            fks_i = len(colors)
+            # use the last colored particle if it exists, or 
+            # just the last
+            fks_j=1
+            for cpos, col in enumerate(colors):
+                if col != 1:
+                    fks_j = cpos+1
+                    fks_j_id = [l.get('id') for l in bornproc.get('legs')][cpos]
+
+            # for the moment, if j is initial-state, we do nothing
+            if fks_j > ninitial:
+                iFKS=1
+                iconfig = 0
+                s_and_t_channels = []
+                mapconfigs = []
+                base_diagrams = born_me.get('base_amplitude').get('diagrams')
+                model = born_me.get('base_amplitude').get('process').get('model')
+                minvert = min([max([len(vert.get('legs')) for vert in \
+                                        diag.get('vertices')]) for diag in base_diagrams])
+        
+                lines.append("# ")
+                lines.append("# nFKSprocess %d" % iFKS)
+                for idiag, diag in enumerate(base_diagrams):
+                    if any([len(vert.get('legs')) > minvert for vert in
+                            diag.get('vertices')]):
+                    # Only 3-vertices allowed in configs.inc
+                        continue
+                    iconfig = iconfig + 1
+                    helas_diag = born_me.get('diagrams')[idiag]
+                    mapconfigs.append(helas_diag.get('number'))
+                    lines.append("# Diagram %d for nFKSprocess %d" % \
+                                     (helas_diag.get('number'),iFKS))
+                    # Correspondance between the config and the amplitudes
+                    lines.append("C   %4d   %4d   %4d " % (iFKS,iconfig,
+                                                           helas_diag.get('number')))
+
+                    # Need to reorganize the topology so that we start with all
+                    # final state external particles and work our way inwards
+                    schannels, tchannels = helas_diag.get('amplitudes')[0].\
+                        get_s_and_t_channels(ninitial, model, 990)
+        
+                    s_and_t_channels.append([schannels, tchannels])
+
+                    #the first thing to write is the splitting ij -> i,j
+                    lines.append("F   %4d   %4d   %4d   %4d" % \
+                                     (iFKS,-1,iconfig,2))
+                                     #(iFKS,last_leg.get('number'), iconfig, len(daughters)))
+                    lines.append("D   %4d" % nexternal)
+                    lines.append("D   %4d" % fks_j)
+                    lines.append("S   %4d   %4d   %4d   %10d" % \
+                                         (iFKS,-1, iconfig,fks_j_id)) 
+                    # now we continue with all the other vertices of the diagrams;
+                    # we need to shift the 'last_leg' by 1 and replace leg fks_j with -1
+
+                    # Write out propagators for s-channel and t-channel vertices
+                    allchannels = schannels
+                    if len(tchannels) > 1:
+                        # Write out tchannels only if there are any non-trivial ones
+                        allchannels = schannels + tchannels
+
+                    for vert in allchannels:
+                        daughters = [leg.get('number') for leg in vert.get('legs')[:-1]]
+                        last_leg = vert.get('legs')[-1]
+                        lines.append("F   %4d   %4d   %4d   %4d" % \
+                                         (iFKS,last_leg.get('number')-1, iconfig, len(daughters)))
+
+                        # legs with negative number in daughters have to be shifted by -1
+                        for i_dau in range(len(daughters)):
+                            if daughters[i_dau] < 0:
+                                daughters[i_dau] += -1
+                        # finally relable fks with -1 if it appears in daughters
+                        if fks_j in daughters:
+                            daughters[daughters.index(fks_j)] = -1
+                        for d in daughters:
+                            lines.append("D   %4d" % d)
+                        if vert in schannels:
+                            lines.append("S   %4d   %4d   %4d   %10d" % \
+                                             (iFKS,last_leg.get('number')-1, iconfig,
+                                              last_leg.get('id')))
+                        elif vert in tchannels[:-1]:
+                            lines.append("T   %4d   %4d   %4d   %10d" % \
+                                             (iFKS,last_leg.get('number')-1, iconfig,
+                                              abs(last_leg.get('id'))))
+
+                        # update what the array sizes (mapconfig,iforest,etc) will be
+                        max_leg_number = min(max_leg_number,last_leg.get('number')-1)
+                    max_iconfig = max(max_iconfig,iconfig)
+        
+                # Write out number of configs
+                lines.append("# Number of configs for nFKSprocess %d" % iFKS)
+                lines.append("C   %4d   %4d   %4d" % (iFKS,0,iconfig))
+
+                # write the props.inc information
+                lines2.append("# ")
+                particle_dict = born_me.get('processes')[0].get('model').\
+                    get('particle_dict')
+        
+                for iconf, configs in enumerate(s_and_t_channels):
+                    lines2.append("M   %4d   %4d   %4d   %10d " % \
+                                    (iFKS,-1, iconf + 1, fks_j_id))
+                    pow_part = 1 + int(particle_dict[fks_j_id].is_boson())
+                    lines2.append("P   %4d   %4d   %4d   %4d " % \
+                                    (iFKS,-1, iconf + 1, pow_part))
+                    for vertex in configs[0] + configs[1][:-1]:
+                        leg = vertex.get('legs')[-1]
+                        if leg.get('id') not in particle_dict:
+                            # Fake propagator used in multiparticle vertices
+                            pow_part = 0
+                        else:
+                            particle = particle_dict[leg.get('id')]
+        
+                            pow_part = 1 + int(particle.is_boson())
+        
+                        lines2.append("M   %4d   %4d   %4d   %10d " % \
+                                         (iFKS,leg.get('number')-1, iconf + 1, leg.get('id')))
+                        lines2.append("P   %4d   %4d   %4d   %4d " % \
+                                         (iFKS,leg.get('number')-1, iconf + 1, pow_part))
     
         # Write the file
         open(filename,'w').write('\n'.join(lines+lines2))
 
-        return max_iconfig, max_leg_number, nconfs
+        return max_iconfig, max_leg_number
 
 
     def write_leshouche_info_declarations(self, writer, nfksconfs, 
@@ -1203,6 +1347,11 @@
             maxproc = max(maxproc, nprocs)
             maxflow = max(maxflow, nflows)
 
+        # this is for LOonly
+        if not matrix_element.get_fks_info_list():
+            (newlines, nprocs, nflows) = self.get_leshouche_lines_dummy(matrix_element.born_matrix_element, 1)
+            lines.extend(newlines)
+
         # Write the file
         open(filename,'w').write('\n'.join(lines))
 
@@ -2206,8 +2355,8 @@
             charges = [0.] * len(colors) 
 
             fks_i = len(colors)
-            # use the first colored particle if it exists, or 
-            # just the first
+            # use the last colored particle if it exists, or 
+            # just the last
             fks_j=1
             for cpos, col in enumerate(colors[:-1]):
                 if col != 1:
@@ -2709,6 +2858,109 @@
         return lines, nproc, nflow
 
 
+    def get_leshouche_lines_dummy(self, matrix_element, ime):
+        #test written
+        """As get_leshouche_lines, but for 'fake' real emission processes (LOonly
+        In this case, write born color structure times ij -> i,j splitting)
+        """
+
+        bornproc = matrix_element.get('processes')[0]
+        colors = [l.get('color') for l in bornproc.get('legs')] 
+
+        fks_i = len(colors)
+        # use the last colored particle if it exists, or 
+        # just the last
+        fks_j=1
+        for cpos, col in enumerate(colors):
+            if col != 1:
+                fks_j = cpos+1
+    
+        # Extract number of external particles
+        (nexternal, ninitial) = matrix_element.get_nexternal_ninitial()
+        nexternal+=1 # remember, in this case matrix_element is born
+    
+        lines = []
+        for iproc, proc in enumerate(matrix_element.get('processes')):
+            # add the fake extra leg
+            legs = proc.get_legs_with_decays() + \
+                   [fks_common.FKSLeg({'id': -21,
+                                       'number': nexternal,
+                                       'state': True,
+                                       'fks': 'i',
+                                       'color': 8,
+                                       'charge': 0.,
+                                       'massless': True,
+                                       'spin': 3,
+                                       'is_part': True,
+                                       'self_antipart': True})]
+                                        
+            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("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!
+            if iproc == 0:
+                # If no color basis, just output trivial color flow
+                if not matrix_element.get('color_basis'):
+                    for i in [1, 2]:
+                        lines.append("C   %4d   %4d   1      %s" % \
+                                 (ime, i, 
+                                  " ".join([ "%3d" % 0 ] * nexternal)))
+                    color_flow_list = []
+                    nflow = 1
+    
+                else:
+                    # in this case the last particle (-21) has two color indices
+                    # and it has to be emitted by j_fks
+                    # First build a color representation dictionnary
+                    repr_dict = {}
+                    for l in legs[:-1]:
+                        repr_dict[l.get('number')] = \
+                            proc.get('model').get_particle(l.get('id')).get_color()\
+                            * (-1)**(1+l.get('state'))
+                    # Get the list of color flows
+                    color_flow_list = \
+                        matrix_element.get('color_basis').color_flow_decomposition(repr_dict,
+                                                                                   ninitial)
+                    # And output them properly
+                    for cf_i, color_flow_dict in enumerate(color_flow_list):
+                        # we have to add the extra leg (-21), linked to the j_fks leg
+                        # first, find the maximum color label
+                        maxicol = max(sum(color_flow_dict.values(), []))
+                        #then, replace the color labels
+                        if color_flow_dict[fks_j][0] == 0:
+                            anti = True
+                            icol_j = color_flow_dict[fks_j][1]
+                        else:
+                            anti = False
+                            icol_j = color_flow_dict[fks_j][0]
+
+                        if anti:
+                            color_flow_dict[nexternal] = (maxicol + 1, color_flow_dict[fks_j][1])
+                            color_flow_dict[fks_j][1] = maxicol + 1
+                        else:
+                            color_flow_dict[nexternal] = (color_flow_dict[fks_j][0], maxicol + 1) 
+                            color_flow_dict[fks_j][0] = maxicol + 1
+
+                        for i in [0, 1]:
+                            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'))
+    
+        return lines, nproc, nflow
+
+
     #===============================================================================
     # get_den_factor_lines
     #===============================================================================
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc	2016-02-27 14:33:01 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_decl.inc	2017-11-30 23:06:43 +0000
@@ -1,12 +1,12 @@
       INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
-      PARAMETER (LMAXCONFIGS_USED=   0)
-      PARAMETER (MAX_BRANCH_USED =   0)
-      INTEGER MAPCONFIG_D(  0,0:LMAXCONFIGS_USED)
-      INTEGER IFOREST_D(  0,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER SPROP_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER TPRID_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PMASS_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PWIDTH_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      INTEGER POW_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      PARAMETER (LMAXCONFIGS_USED=   3)
+      PARAMETER (MAX_BRANCH_USED =   3)
+      INTEGER MAPCONFIG_D(  1,0:LMAXCONFIGS_USED)
+      INTEGER IFOREST_D(  1,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER SPROP_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER TPRID_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PMASS_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PWIDTH_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      INTEGER POW_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%configs_and_props_info.dat	2017-11-30 23:06:43 +0000
@@ -3,4 +3,57 @@
 # S -> SPROP_D
 # T -> TPRID_D
 # M -> PMASS_D/PWIDTH_D
-# P -> POW_D
 \ No newline at end of file
+# P -> POW_D
+# 
+# nFKSprocess 1
+# Diagram 1 for nFKSprocess 1
+C      1      1      1 
+F      1     -1      1      2
+D      5
+D      4
+S      1     -1      1           -6
+F      1     -2      1      2
+D     -1
+D      3
+S      1     -2      1           21
+# Diagram 2 for nFKSprocess 1
+C      1      2      2 
+F      1     -1      2      2
+D      5
+D      4
+S      1     -1      2           -6
+F      1     -2      2      2
+D      1
+D      3
+T      1     -2      2            6
+F      1     -3      2      2
+D     -2
+D     -1
+# Diagram 3 for nFKSprocess 1
+C      1      3      3 
+F      1     -1      3      2
+D      5
+D      4
+S      1     -1      3           -6
+F      1     -2      3      2
+D      1
+D     -1
+T      1     -2      3            6
+F      1     -3      3      2
+D     -2
+D      3
+# Number of configs for nFKSprocess 1
+C      1      0      3
+# 
+M      1     -1      1           -6 
+P      1     -1      1      1 
+M      1     -2      1           21 
+P      1     -2      1      2 
+M      1     -1      2           -6 
+P      1     -1      2      1 
+M      1     -2      2            6 
+P      1     -2      2      1 
+M      1     -1      3           -6 
+P      1     -1      3      1 
+M      1     -2      3           -6 
+P      1     -2      3      1 
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%leshouche_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%leshouche_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%leshouche_info.dat	2017-11-30 23:06:43 +0000
@@ -1,3 +1,10 @@
 # I -> IDUP_D
 # M -> MOTHUP_D
-# C -> ICOLUP_D
 \ No newline at end of file
+# C -> ICOLUP_D
+I      1      1       21 21 6 -6 -21
+M      1      1      1        0   0   1   1   1
+M      1      2      1        0   0   2   2   2
+C      1      1      1      501 502 501   0 504
+C      1      2      1      502 503   0 504 503
+C      1      1      2      503 501 501   0 504
+C      1      2      2      502 503   0 504 502
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%ngraphs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%ngraphs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%ngraphs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,2 @@
       INTEGER    N_MAX_CG
-      PARAMETER (N_MAX_CG=0)
+      PARAMETER (N_MAX_CG=3)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%real_from_born_configs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%real_from_born_configs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%real_from_born_configs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,3 @@
       INTEGER IRFBC
-      INTEGER REAL_FROM_BORN_CONF(0,0)
+      INTEGER REAL_FROM_BORN_CONF(3,1)
+      DATA (REAL_FROM_BORN_CONF(IRFBC,1),IRFBC=1,3) /1, 2, 3/
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc	2016-02-27 14:33:01 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_decl.inc	2017-11-30 23:06:43 +0000
@@ -1,12 +1,12 @@
       INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
-      PARAMETER (LMAXCONFIGS_USED=   0)
-      PARAMETER (MAX_BRANCH_USED =   0)
-      INTEGER MAPCONFIG_D(  0,0:LMAXCONFIGS_USED)
-      INTEGER IFOREST_D(  0,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER SPROP_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER TPRID_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PMASS_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PWIDTH_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      INTEGER POW_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      PARAMETER (LMAXCONFIGS_USED=   1)
+      PARAMETER (MAX_BRANCH_USED =   2)
+      INTEGER MAPCONFIG_D(  1,0:LMAXCONFIGS_USED)
+      INTEGER IFOREST_D(  1,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER SPROP_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER TPRID_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PMASS_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PWIDTH_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      INTEGER POW_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%configs_and_props_info.dat	2017-11-30 23:06:43 +0000
@@ -3,4 +3,23 @@
 # S -> SPROP_D
 # T -> TPRID_D
 # M -> PMASS_D/PWIDTH_D
-# P -> POW_D
 \ No newline at end of file
+# P -> POW_D
+# 
+# nFKSprocess 1
+# Diagram 1 for nFKSprocess 1
+C      1      1      1 
+F      1     -1      1      2
+D      5
+D      4
+S      1     -1      1           -6
+F      1     -2      1      2
+D     -1
+D      3
+S      1     -2      1           21
+# Number of configs for nFKSprocess 1
+C      1      0      1
+# 
+M      1     -1      1           -6 
+P      1     -1      1      1 
+M      1     -2      1           21 
+P      1     -2      1      2 
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%leshouche_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%leshouche_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%leshouche_info.dat	2017-11-30 23:06:43 +0000
@@ -1,3 +1,19 @@
 # I -> IDUP_D
 # M -> MOTHUP_D
-# C -> ICOLUP_D
 \ No newline at end of file
+# C -> ICOLUP_D
+I      1      1       2 -2 6 -6 -21
+M      1      1      1        0   0   1   1   1
+M      1      2      1        0   0   2   2   2
+C      1      1      1      501   0 502   0 503
+C      1      2      1        0 501   0 503 502
+C      1      1      2      502   0 502   0 503
+C      1      2      2        0 501   0 503 501
+I      1      2       4 -4 6 -6 -21
+M      1      1      2        0   0   1   1   1
+M      1      2      2        0   0   2   2   2
+I      1      3       1 -1 6 -6 -21
+M      1      1      3        0   0   1   1   1
+M      1      2      3        0   0   2   2   2
+I      1      4       3 -3 6 -6 -21
+M      1      1      4        0   0   1   1   1
+M      1      2      4        0   0   2   2   2
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%ngraphs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%ngraphs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%ngraphs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,2 @@
       INTEGER    N_MAX_CG
-      PARAMETER (N_MAX_CG=0)
+      PARAMETER (N_MAX_CG=1)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%real_from_born_configs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%real_from_born_configs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%real_from_born_configs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,3 @@
       INTEGER IRFBC
-      INTEGER REAL_FROM_BORN_CONF(0,0)
+      INTEGER REAL_FROM_BORN_CONF(1,1)
+      DATA (REAL_FROM_BORN_CONF(IRFBC,1),IRFBC=1,1) /1/
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_decl.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_decl.inc	2016-02-27 14:33:01 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_decl.inc	2017-11-30 23:06:43 +0000
@@ -1,12 +1,12 @@
       INTEGER IFR,LMAXCONFIGS_USED,MAX_BRANCH_USED
-      PARAMETER (LMAXCONFIGS_USED=   0)
-      PARAMETER (MAX_BRANCH_USED =   0)
-      INTEGER MAPCONFIG_D(  0,0:LMAXCONFIGS_USED)
-      INTEGER IFOREST_D(  0,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER SPROP_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      INTEGER TPRID_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PMASS_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      DOUBLE PRECISION PWIDTH_D(  0,-MAX_BRANCH_USED:-1
-     $ ,LMAXCONFIGS_USED)
-      INTEGER POW_D(  0,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      PARAMETER (LMAXCONFIGS_USED=   1)
+      PARAMETER (MAX_BRANCH_USED =   2)
+      INTEGER MAPCONFIG_D(  1,0:LMAXCONFIGS_USED)
+      INTEGER IFOREST_D(  1,2,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER SPROP_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      INTEGER TPRID_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PMASS_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      DOUBLE PRECISION PWIDTH_D(  1,-MAX_BRANCH_USED:-1
+     $ ,LMAXCONFIGS_USED)
+      INTEGER POW_D(  1,-MAX_BRANCH_USED:-1,LMAXCONFIGS_USED)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%configs_and_props_info.dat	2017-11-30 23:06:43 +0000
@@ -3,4 +3,23 @@
 # S -> SPROP_D
 # T -> TPRID_D
 # M -> PMASS_D/PWIDTH_D
-# P -> POW_D
 \ No newline at end of file
+# P -> POW_D
+# 
+# nFKSprocess 1
+# Diagram 1 for nFKSprocess 1
+C      1      1      1 
+F      1     -1      1      2
+D      5
+D      4
+S      1     -1      1           -6
+F      1     -2      1      2
+D     -1
+D      3
+S      1     -2      1           21
+# Number of configs for nFKSprocess 1
+C      1      0      1
+# 
+M      1     -1      1           -6 
+P      1     -1      1      1 
+M      1     -2      1           21 
+P      1     -2      1      2 
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%leshouche_info.dat'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%leshouche_info.dat	2017-05-15 10:39:31 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%leshouche_info.dat	2017-11-30 23:06:43 +0000
@@ -1,3 +1,19 @@
 # I -> IDUP_D
 # M -> MOTHUP_D
-# C -> ICOLUP_D
 \ No newline at end of file
+# C -> ICOLUP_D
+I      1      1       -2 2 6 -6 -21
+M      1      1      1        0   0   1   1   1
+M      1      2      1        0   0   2   2   2
+C      1      1      1        0 501 502   0 503
+C      1      2      1      501   0   0 503 502
+C      1      1      2        0 502 502   0 503
+C      1      2      2      501   0   0 503 501
+I      1      2       -4 4 6 -6 -21
+M      1      1      2        0   0   1   1   1
+M      1      2      2        0   0   2   2   2
+I      1      3       -1 1 6 -6 -21
+M      1      1      3        0   0   1   1   1
+M      1      2      3        0   0   2   2   2
+I      1      4       -3 3 6 -6 -21
+M      1      1      4        0   0   1   1   1
+M      1      2      4        0   0   2   2   2
 \ No newline at end of file
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%ngraphs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%ngraphs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%ngraphs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,2 @@
       INTEGER    N_MAX_CG
-      PARAMETER (N_MAX_CG=0)
+      PARAMETER (N_MAX_CG=1)
 
=== modified file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%real_from_born_configs.inc'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%real_from_born_configs.inc	2015-03-10 08:54:03 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%real_from_born_configs.inc	2017-11-30 23:06:43 +0000
@@ -1,2 +1,3 @@
       INTEGER IRFBC
-      INTEGER REAL_FROM_BORN_CONF(0,0)
+      INTEGER REAL_FROM_BORN_CONF(1,1)
+      DATA (REAL_FROM_BORN_CONF(IRFBC,1),IRFBC=1,1) /1/
 
=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%genps.inc.moved'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%genps.inc.moved	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%genps.inc.moved	2017-11-30 23:06:43 +0000
@@ -0,0 +1,4 @@
+      INCLUDE 'maxparticles.inc'
+      INCLUDE 'maxconfigs.inc'
+      INTEGER MAXPROC,NGRAPHS,NCOLOR,MAXFLOW
+      PARAMETER (MAXPROC=4,NGRAPHS=3,NCOLOR=2,MAXFLOW=6)
 
=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%genps.inc.moved'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%genps.inc.moved	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%genps.inc.moved	2017-11-30 23:06:43 +0000
@@ -0,0 +1,4 @@
+      INCLUDE 'maxparticles.inc'
+      INCLUDE 'maxconfigs.inc'
+      INTEGER MAXPROC,NGRAPHS,NCOLOR,MAXFLOW
+      PARAMETER (MAXPROC=4,NGRAPHS=1,NCOLOR=2,MAXFLOW=4)
 
=== added file 'tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%genps.inc.moved'
--- tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%genps.inc.moved	1970-01-01 00:00:00 +0000
+++ tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%genps.inc.moved	2017-11-30 23:06:43 +0000
@@ -0,0 +1,4 @@
+      INCLUDE 'maxparticles.inc'
+      INCLUDE 'maxconfigs.inc'
+      INTEGER MAXPROC,NGRAPHS,NCOLOR,MAXFLOW
+      PARAMETER (MAXPROC=4,NGRAPHS=1,NCOLOR=2,MAXFLOW=4)
 
=== modified file 'tests/unit_tests/iolibs/test_export_fks.py'
--- tests/unit_tests/iolibs/test_export_fks.py	2016-08-31 20:35:13 +0000
+++ tests/unit_tests/iolibs/test_export_fks.py	2017-11-30 23:06:43 +0000
@@ -92,6 +92,53 @@
                 self.assertEqual(old_l, new_l)
 
 
+    def test_wj_loonly_gen(self):
+        """check that the new (memory and cpu efficient) and old generation
+        mode at NLO give the same results for p p > w j [LOonly=QCD]
+        """
+        path = tempfile.mkdtemp('', 'TMPWTest', None)
+
+        def run_cmd(cmd):
+            interface.exec_cmd(cmd, errorhandling=False, printcmd=False, 
+                               precmd=True, postcmd=True)
+
+        interface = MGCmd.MasterCmd()
+        
+        run_cmd('generate p p > w+ j [LOonly=QCD]')
+        run_cmd('output %s' % os.path.join(path, 'W-oldway'))
+        run_cmd('set low_mem_multicore_nlo_generation True')
+        run_cmd('generate p p > w+ j [LOonly=QCD]')
+        run_cmd('output %s' % os.path.join(path, 'W-newway'))
+        run_cmd('set low_mem_multicore_nlo_generation False')
+        
+        # the P0 dirs
+        for oldf in \
+          (glob.glob(os.path.join(path, 'W-oldway', 'SubProcesses', 'P0*', '*.inc')) + \
+           glob.glob(os.path.join(path, 'W-oldway', 'SubProcesses', 'P0*', '*.f')) + \
+           [os.path.join(path, 'W-oldway', 'SubProcesses', 'proc_characteristics')]):
+            
+            if os.path.islink(oldf): 
+                continue
+
+            newf = oldf.replace('oldway', 'newway')
+
+            for old_l, new_l in zip(open(oldf), open(newf)):
+                self.assertEqual(old_l, new_l)
+
+        # the V0 dirs
+        for oldf in \
+          (glob.glob(os.path.join(path, 'W-oldway', 'SubProcesses', 'P0*', 'V0*', '*.inc')) + \
+           glob.glob(os.path.join(path, 'W-oldway', 'SubProcesses', 'P0*', 'V0*', '*.f'))):
+            
+            if os.path.islink(oldf): 
+                continue
+
+            newf = oldf.replace('oldway', 'newway')
+
+            for old_l, new_l in zip(open(oldf), open(newf)):
+                self.assertEqual(old_l, new_l)
+
+
     def test_w_nlo_gen_gosam(self):
         """check that the new generation mode works when gosam is set 
         for p p > w [QCD]