Merge lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping into lp:~madteam/mg5amcnlo/trunk

[Olivier Mattelaer]

Call the C routine if pair of majorana have the wrong input-output order according to the interactions.

Note this also require to modify the routine defining the sign of the coupling for the octet majorana.

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

191. By Johan Alwall on 2011-11-08

New implementation of check for hermition conjugate HELAS function based on Majorana order, including test

192. By Johan Alwall on 2011-11-08

Merged with latest version of the trunk

193. By Johan Alwall on 2011-11-08

Updated UpdateNotes and VERSION

[Johan Alwall (johan-alwall)]

review: Approve

[Johan Alwall (johan-alwall)]

Olivier and Tim,

I have redone the implementation in a minimal manner. This way, none of our previous functionality is affected, which I believe is the most safe way to do it (especially considering the *enormous* amount of time I have spent on getting the fermion flow correct for all cases). I certainly don't dare touch that part of the code unless it's absolutely necessary.

My new implementation is minimal and should be self-explanatory. None of the old tests failed (this was the most worrying part of Olivier's implementation - many of the old tests failed, although they should have no relation to this issue), and I added a test - please feel free to go through the test in detail and make sure this is exactly the desired behavior.

Olivier, can you also please run the full set of MSSM parallel tests (in particular, of course, those involving neutralinos), comparing with the old version of MG5 (or MG4, which I think should be equivalent). This is a considerable modification, which new helamp calls whenever there are two different neutralinos and more than 2 final state particles. So it's truly crucial to make sure that the results for regular MSSM have not changed.

I have tested e+ e- > n1 n2 z z, where all relevant conjugated vertices appear, and got identical result with the old MSSM. But I'm not sure enough that this guarantees correctness for all processes.

Many thanks,
Johan

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

194. By Johan Alwall on 2011-11-08

Removed erraneous is_majorana function (checking for identity of name and antiname instead of self_antipart)

195. By Johan Alwall on 2011-11-08

Sorry doc line on majorana_conjugates was not up-to-date

196. By Johan Alwall on 2011-11-08

Fixed one more tiny thing I missed

197. By omatt on 2011-11-08

simple improvment (nicer code and faster as well)

[Olivier Mattelaer (olivier-mattelaer)]

Hi Johan,

>I have redone the implementation in a minimal manner. This way, none of our previous >functionality is affected, which I believe is the most safe way to do it (especially >considering
>the *enormous* amount of time I have spent on getting the fermion flow correct for all cases).I
>certainly don't dare touch that part of the code unless it's absolutely necessary.

At the end this is quite similar to what I did.
- one modification to need_hermitian (both).
- one modification to need_charge_conjugate (both).
So I don't understand your comment above.

The main differences between those two versions are
1) that I didn't veto the Conjugate on identical majorana (which creates the trouble with the tests that you mention)
2) that you factorize everything in one function (see below)
3) that you treat 4 fermion case, (Excellent that you did)

But I really consider my implementation much more simpler to understand. (Off course I'm biased)

I would strongly suggest to split majorana_conjugates
in two separate functions. One for the wavefunctions and one for the amplitudes.
the number of if in this case checking in which case we are (and some of them quite implicit)
makes the code un-readable.

Please also split the long if. This will make the code more readable.
some if test in fact two completely different things
1) if both particle are majorana
2) if the flow is correct
make two if is maybe longer (and with additional else) but much more clear.

Another problem, but with wich we can life with, is that they the structure is now quite un-efficient.

in get_call_key:
1) it check if it need the hermitian
--> one call the majorana_conjugates
2) it ask for the conjugate index
--> call to get_conjugate_index
    This one check first if if need the hermitian conjugate
    --> second call to majorana_conjugates
    Then call itself majorana_conjugates
    --> third call to majorana_conjugates

Then at a point we will call get_used_lorentz
which will make a call to get_conjugate_index
So two additional call to majorana_conjugates

So five calls to this functions this is a bit pointless.
We can:
1) store the data (forbids the static method but fine if you split the routine in two as suggested above) such that the data is compute only once.
2) Refactorize the code to remove the need_conjugate routine.
(We will pass from 5 to 2 call). -- I prefer this one since the modif are very small at the end--
3) keep like that (As I said this is not critical for the moment at least)

Finally, your test is about MG4 Model, which are in fact insensitive to this problem in mg5.
(Since the call for such routine are not implemented) Are they a good reason for not using UFO syntax in your tests?

Cheers,

Olivier

review: Needs Fixing

[Johan Alwall (johan-alwall)]

Hello Olivier,

Many thanks for your comments.

> >I have redone the implementation in a minimal manner. This way, none of our
> previous >functionality is affected, which I believe is the most safe way to
> do it (especially >considering
> >the *enormous* amount of time I have spent on getting the fermion flow
> correct for all cases).I
> >certainly don't dare touch that part of the code unless it's absolutely
> necessary.
>
> At the end this is quite similar to what I did.
> - one modification to need_hermitian (both).
> - one modification to need_charge_conjugate (both).
> So I don't understand your comment above.

You also touched check_and_fix_fermion_flow, which is very dangerous (this I know from experience - as I said, this part of the code is among the most complicated and I have been working this over many times already without managing to simplify it, fixing many bugs). This, and the host of failing tests, and the obvious problem that you had to touch the sign for gluinos made me prefer to redo the implementation in a way that made more sense to me. Of course, the honor for solving the problem (using conjugate wfs) all goes to you. :-)

> The main differences between those two versions are
> 1) that I didn't veto the Conjugate on identical majorana (which creates the
> trouble with the tests that you mention)
> 2) that you factorize everything in one function (see below)
> 3) that you treat 4 fermion case, (Excellent that you did)
>
> But I really consider my implementation much more simpler to understand. (Off
> course I'm biased)

My solution exactly mirrors what is done in sorted_mothers, so basically, if you understand one you understand the other. So clearly also sorted_mothers is not clear (see below).

> I would strongly suggest to split majorana_conjugates
> in two separate functions. One for the wavefunctions and one for the
> amplitudes.
> the number of if in this case checking in which case we are (and some of them
> quite implicit)
> makes the code un-readable.

Do you want me to do this also for sorted_mothers? The idea was to make the code as DRY as possible, but indeed the repetition is quite minor at this point (since I already at some point extracted sort_by_pdg_codes from this function).

> Please also split the long if. This will make the code more readable.
> some if test in fact two completely different things
> 1) if both particle are majorana
> 2) if the flow is correct
> make two if is maybe longer (and with additional else) but much more clear.

Ok, sure.

> Another problem, but with wich we can life with, is that they the structure is
> now quite un-efficient.
>
> in get_call_key:
> 1) it check if it need the hermitian
> --> one call the majorana_conjugates
> 2) it ask for the conjugate index
> --> call to get_conjugate_index
> This one check first if if need the hermitian conjugate
> --> second call to majorana_conjugates
> Then call itself majorana_conjugates
> --> third call to majorana_conjugates
>
> Then at a point we will call get_used_lorentz
> which will make a call to get_conjugate_index
> So two additional call to majorana_conjugates
>
> So five calls to this functions this ...

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

198. By Johan Alwall on 2011-11-09

Simplified structure for Majorana conjugates, and included pretty complete test for this. Changed to UFO format for test_majorana_conjugate_process.

[Johan Alwall (johan-alwall)]

Hello again,

> > I would strongly suggest to split majorana_conjugates
> > in two separate functions. One for the wavefunctions and one for the
> > amplitudes.
> > the number of if in this case checking in which case we are (and some of
> them
> > quite implicit)
> > makes the code un-readable.

I have taken care of this, I hope you like my new implementation better. It should be considerably more intuitive and readable.

> > We can:
> > 1) store the data (forbids the static method but fine if you split the
> routine
> > in two as suggested above) such that the data is compute only once.

Done.

Thanks!
Johan

[Olivier Mattelaer (olivier-mattelaer)]

Excellent!!!!!!!
This is very very good.
Good work.

I like it a lot!!!!

Let me run the parralel test on it and I'll approve it when this is done.

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

199. By mattelaer-olivier on 2011-11-10

Implemetation of a new parralel tests

200. By mattelaer-olivier on 2011-11-10

Implement all the test in the new way (comparaison with 1.3.3)
Add a small parralel test in the create release (comparaison with 1.3.28 via piclke file)

[Olivier Mattelaer (olivier-mattelaer)]

I've include in this branch the modifications with the parralel tests.

This include:
1) most of the full suite is runned against 1.3.3 (with model version of 1.3.3)
2) a short suite is runned against pickle files (computed with the current version 1.3.28 with current model)
3) The release script launch automatically the short suite.

Otherwise the parralel test is still running (I launched it yesterday night).

[Johan Alwall (johan-alwall)]

Hello Olivier,

Awesome job for the parallel tests. I haven't tried it yet though. Just one question:
586 + filename = "short_sm.log",
isn't it supposed to be "short_mssm.log"? Same thing for heft.

Cheers,
Johan

[Olivier Mattelaer (olivier-mattelaer)]

Hi I've simplify/change the parralel test for the 2 to 3 mssm.
This should be faster.
and at least, you will have faster information if it works or not.
So I relaunched that parralel test. And so far, he succesfully passed on thousands of them.

So I'll approve the review. (Off course if a surprise appear in the tests, I will shout).

review: Approve

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

201. By mattelaer-olivier on 2011-11-11

simplify mssm parralel tests
support the case where no processes are valid in the parralel tests

202. By Johan Alwall on 2011-11-12

Fix color for multiparticle vertices (and regular vertices) in most general way in addmothers.f; fix bug for processes with multiple multiparticle vertices

[Johan Alwall (johan-alwall)]

Hello Olivier,

I have completely redone the color treatment in addmothers.f in what I believe is the most general way possible. Now, color should work for absolutely any process (with color 1,3,6,8 particles). Please do a set of MadEvent runs with different color combinations to make sure there's no trivial bug left.

Thanks,
Johan

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

203. By Johan Alwall on 2011-11-12

Added error-induced stop if color assignment fails

204. By Johan Alwall on 2011-11-13

Take into account min_events from run_config.inc in gridpack refine

[Olivier Mattelaer (olivier-mattelaer)]

Hi Johan,

All my tests passes without any trouble,
but I don't fully understand how the code works.
Could you add some comments in the code to explain
the philosophy of this new routines and the meaning
of the different variable.

I'm pretty sure that with such global comment, I'll be able to fully
understand your code. Maybe it's better to wait that I understand it before
release it. But If you want you can already release it.

On my side, I'll merge this version with the 1.4.0 such that I'll be able to run parralel tests during the night.

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/fix_majorana_flipping]

205. By Johan Alwall on 2011-11-14

Added more comments to appease Olivier :-)

=== modified file 'Template/Source/gen_ximprove.f'
--- Template/Source/gen_ximprove.f	2011-11-04 15:04:11 +0000
+++ Template/Source/gen_ximprove.f	2011-11-14 05:57:24 +0000
@@ -827,8 +827,8 @@
             write(26,20) 'if [[ ! -e ftn25 ]]; then'
 
 
-            write(26,'(9x,a,2i8,a)') 'echo "',npoints,max_iter,
-     $           '" >& input_sg.txt' 
+            write(26,'(9x,a,2i8,a)') 'echo "',max(npoints,min_events),
+     $           max_iter,'" >& input_sg.txt' 
 c
 c     tjs 8/7/2007  Allow stop when have enough events
 c
@@ -850,8 +850,8 @@
 
             write(26,25) 'rm -f $k'
             
-            write(26,'(9x,a,2i8,a)') 'echo "',npoints,max_iter,
-     $           '" >& input_sg.txt' 
+            write(26,'(9x,a,2i8,a)') 'echo "',max(npoints,min_events),
+     $           max_iter,'" >& input_sg.txt' 
 c
 c tjs 8/7/2007    Change to request events not accuracy
 c

=== modified file 'Template/SubProcesses/addmothers.f'
--- Template/SubProcesses/addmothers.f	2011-10-28 10:20:10 +0000
+++ Template/SubProcesses/addmothers.f	2011-11-14 05:57:24 +0000
@@ -25,6 +25,11 @@
       double precision xtarget
       data iseed/0/
 
+c     Variables for combination of color indices (including multipart. vert)
+      integer maxcolmp
+      parameter(maxcolmp=20)
+      integer ncolmp,icolmp(2,maxcolmp)
+
       double precision ZERO
       parameter (ZERO=0d0)
       double precision pmass(-nexternal:0,lmaxconfigs)
@@ -67,9 +72,9 @@
 c      common/to_colstats/ncols,ncolflow,ncolalt,icorg
 
       double precision pt
-      integer get_color
+      integer get_color,elim_indices
       real ran1
-      external pt,ran1,get_color
+      external pt,ran1,get_color,elim_indices
 
       if (first_time) then
          include 'props.inc'
@@ -188,8 +193,9 @@
           mo_color = get_color(jpart(1,i))
           da_color(1) = get_color(jpart(1,ida(1)))
           da_color(2) = get_color(jpart(1,ida(2)))
-          if(da_color(2).lt.da_color(1))then
-c            Order daughters according to color
+          if(da_color(1).ne.2.and.da_color(2).lt.da_color(1).or.
+     $         da_color(2).eq.2)then
+c            Order daughters according to color, but always color 2 first
              itmp=ida(1)
              ida(1)=ida(2)
              ida(2)=itmp
@@ -197,149 +203,48 @@
              da_color(1)=da_color(2)
              da_color(2)=itmp
           endif
-c          print *,'graph: ',iconfig
-c          print *,'Resonance: ',i,' daughters ',ida(1),ida(2),
-c     $         ' ids ',jpart(1,i),jpart(1,ida(1)),jpart(1,ida(2)),
-c     $         ' colors ',mo_color,da_color(1),da_color(2)
-          
+c     Reset list of color indices if not inside multipart. vertex
+c     (indicated by color 2)
+          if(da_color(1).ne.2)then
+             ncolmp=0
+          endif
+c     Add new color indices to list of color indices
+c     Note that color=2 means continued multiparticle index
+          do j=1,2
+             if(da_color(j).eq.2.or.da_color(j).eq.1) cycle
+             ncolmp=ncolmp+1
+             icolmp(1,ncolmp)=icolalt(1,ida(j))
+             icolmp(2,ncolmp)=icolalt(2,ida(j))
+c            Avoid color sextet-type negative indices
+             if(icolmp(1,ncolmp).lt.0)then
+                ncolmp=ncolmp+1
+                icolmp(2,ncolmp)=-icolmp(1,ncolmp-1)
+                icolmp(1,ncolmp-1)=0
+                icolmp(1,ncolmp)=0
+             elseif(icolmp(2,ncolmp).lt.0)then
+                ncolmp=ncolmp+1
+                icolmp(1,ncolmp)=-icolmp(2,ncolmp-1)
+                icolmp(2,ncolmp-1)=0
+                icolmp(2,ncolmp)=0
+             endif
+             if(ncolmp.gt.maxcolmp)
+     $            call write_error(1000,ncolmp,maxcolmp)
+          enddo
+
           if(mo_color.eq.1) then ! color singlet
              icolalt(1,i) = 0
              icolalt(2,i) = 0
-          elseif(mo_color.eq.2) then ! used as dummy for multipart. vert
-             if(icolalt(1,ida(1)).eq.icolalt(2,ida(2)).and.
-     $            icolalt(2,ida(1)).eq.icolalt(1,ida(2)).or.
-     $            icolalt(1,ida(1)).eq.-icolalt(1,ida(2)).and.
-     $            icolalt(2,ida(1)).eq.-icolalt(2,ida(2))) then ! color singlet
-                icolalt(1,i) = 0
-                icolalt(2,i) = 0            
-             elseif(icolalt(1,ida(1)).eq.icolalt(2,ida(2))) then ! 3bar 3 -> 8 or 8 8 -> 8
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(1))
-             else if(icolalt(1,ida(2)).eq.icolalt(2,ida(1))) then ! 3 3bar -> 8 or 8 8 -> 8
-                icolalt(1,i) = icolalt(1,ida(1))
-                icolalt(2,i) = icolalt(2,ida(2))
-             else if(icolalt(1,ida(1)).eq.0.and.icolalt(2,ida(1)).eq.0) then ! 1 3/6/8 -> 3/6/8
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(2))
-             else if(icolalt(1,ida(2)).eq.0.and.icolalt(2,ida(2)).eq.0) then ! 3/6/8 1 -> 3/6/8
-                icolalt(1,i) = icolalt(1,ida(1))
-                icolalt(2,i) = icolalt(2,ida(1))
-             elseif(icolalt(1,ida(1)).eq.-icolalt(1,ida(2))) then ! antisextet
-                icolalt(1,i) = -icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(1))
-             elseif(icolalt(2,ida(1)).eq.-icolalt(2,ida(2))) then ! sextet
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = -icolalt(1,ida(1))
-             else
-c               This might be due to multiparticle vertex with color combinations
-c               that can't be captured with just two sets of color indices
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
           elseif(mo_color.eq.-3) then ! color anti-triplet
-                icolalt(1,i) = 0
-             if(da_color(1).eq.-3.and.da_color(2).eq.1)then
-                icolalt(2,i) = icolalt(2,ida(1))
-             elseif(da_color(1).eq.-3.and.da_color(2).eq.8)then
-                icolalt(2,i) = icolalt(2,ida(2))
-             elseif(da_color(1).eq.-6.and.da_color(2).eq.3)then
-                if(-icolalt(1,ida(1)).eq.icolalt(1,ida(2)))then
-                   icolalt(2,i) = icolalt(2,ida(1))
-                else
-                   icolalt(2,i) = -icolalt(1,ida(1))
-                endif
-             elseif(da_color(1).eq.3.and.da_color(2).eq.3)then
-                maxcolor=maxcolor+1
-                icolalt(2,i) = maxcolor
-             else
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
+             maxcolor=elim_indices(0,1,ncolmp,icolmp,icolalt(1,i),maxcolor)
           elseif(mo_color.eq.3) then ! color triplet
-                icolalt(2,i) = 0
-             if(da_color(1).eq.1.and.da_color(2).eq.3)then
-                icolalt(1,i) = icolalt(1,ida(2))
-             elseif(da_color(1).eq.3.and.da_color(2).eq.8)then
-                icolalt(1,i) = icolalt(1,ida(2))
-             elseif(da_color(1).eq.-3.and.da_color(2).eq.6)then
-                if(icolalt(2,ida(1)).eq.icolalt(1,ida(2)))then
-                   icolalt(1,i) = -icolalt(2,ida(2))
-                else
-                   icolalt(1,i) = icolalt(1,ida(2))
-                endif
-             elseif(da_color(1).eq.-3.and.da_color(2).eq.-3)then
-                maxcolor=maxcolor+1
-                icolalt(1,i) = maxcolor
-             else
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
+             maxcolor=elim_indices(1,0,ncolmp,icolmp,icolalt(1,i),maxcolor)
           elseif(mo_color.eq.-6) then ! color anti-sextet
-             if(da_color(1).eq.-6.and.da_color(2).eq.1)then
-                icolalt(1,i) = icolalt(1,ida(1))
-                icolalt(2,i) = icolalt(2,ida(1))
-             elseif(da_color(1).eq.-6.and.da_color(2).eq.8)then
-                if(icolalt(2,ida(1)).eq.icolalt(1,ida(2)))then
-                   icolalt(1,i) = icolalt(1,ida(1))
-                   icolalt(2,i) = icolalt(2,ida(2))
-                else
-                   icolalt(1,i) = -icolalt(2,ida(2))
-                   icolalt(2,i) = icolalt(2,ida(1))
-                endif
-             elseif(da_color(1).eq.-3.and.da_color(2).eq.-3)then
-                icolalt(1,i) = -icolalt(2,ida(1))
-                icolalt(2,i) = icolalt(2,ida(2))
-             else
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
+             maxcolor=elim_indices(0,2,ncolmp,icolmp,icolalt(1,i),maxcolor)
           elseif(mo_color.eq.6) then ! color sextet
-             if(da_color(1).eq.1.and.da_color(2).eq.6)then
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(2))
-             elseif(da_color(1).eq.6.and.da_color(2).eq.8)then
-                if(icolalt(1,ida(1)).eq.icolalt(2,ida(2)))then
-                   icolalt(1,i) = icolalt(1,ida(2))
-                   icolalt(2,i) = icolalt(2,ida(1))
-                else
-                   icolalt(1,i) = icolalt(1,ida(1))
-                   icolalt(2,i) = -icolalt(1,ida(2))
-                endif
-             elseif(da_color(1).eq.3.and.da_color(2).eq.3)then
-                icolalt(1,i) = icolalt(1,ida(1))
-                icolalt(2,i) = -icolalt(1,ida(2))
-             else
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
+             maxcolor=elim_indices(2,0,ncolmp,icolmp,icolalt(1,i),maxcolor)
           elseif(mo_color.eq.8) then ! color octet
-             if(da_color(1).eq.-3.and.da_color(2).eq.3)then
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(1))             
-             elseif(da_color(1).eq.1.and.da_color(2).eq.8)then
-                icolalt(1,i) = icolalt(1,ida(2))
-                icolalt(2,i) = icolalt(2,ida(2))             
-             elseif(da_color(1).eq.8.and.da_color(2).eq.8)then
-                if(icolalt(1,ida(1)).eq.icolalt(2,ida(2)))then
-                   icolalt(1,i) = icolalt(1,ida(2))
-                   icolalt(2,i) = icolalt(2,ida(1))
-                else
-                   icolalt(1,i) = icolalt(1,ida(1))
-                   icolalt(2,i) = icolalt(2,ida(2))
-                endif
-             elseif(da_color(1).eq.-6.and.da_color(2).eq.6)then
-                if(-icolalt(1,ida(1)).eq.icolalt(1,ida(2)))then
-                   icolalt(1,i) = -icolalt(1,ida(2))
-                   icolalt(2,i) = icolalt(2,ida(1))
-                elseif(icolalt(1,ida(1)).eq.icolalt(2,ida(2)))then
-                   icolalt(1,i) = icolalt(1,ida(2))
-                   icolalt(2,i) = icolalt(2,ida(1))
-                elseif(icolalt(2,ida(1)).eq.icolalt(1,ida(2)))then
-                   icolalt(1,i) = -icolalt(2,ida(2))
-                   icolalt(2,i) = -icolalt(1,ida(1))
-                else
-                   icolalt(1,i) = icolalt(1,ida(2))
-                   icolalt(2,i) = -icolalt(1,ida(1))
-                endif
-             else
-                call write_error(da_color(1), da_color(2), mo_color)
-             endif
-          else
+             maxcolor=elim_indices(1,1,ncolmp,icolmp,icolalt(1,i),maxcolor)
+          elseif(mo_color.ne.2) then ! 2 indicates multipart. vertex
              call write_error(da_color(1), da_color(2), mo_color)
           endif
 c         Set tentative mothers
@@ -406,6 +311,18 @@
       integer ida1,ida2,imo
 
       open(unit=26,file='../../../error',status='unknown',err=999)
+      if (ida1.eq.1000)then
+         write(26,*) 'Error: too many particles in multipart. vertex,',
+     $        ' please increase maxcolmp in addmothers.f'
+         write(*,*) 'Error: too many particles in multipart. vertex,',
+     $        ' please increase maxcolmp in addmothers.f'
+         stop
+      endif
+      if (ida1.eq.1001)then
+         write(26,*) 'Error: failed to reduce to color indices: ',ida2,imo
+         write(*,*) 'Error: failed to reduce to color indices: ',ida2,imo
+         stop
+      endif
       write(26,*) 'Error: Color combination ',ida1,ida2,
      $     '->',imo,' not implemented in addmothers.f'
       write(*,*) 'Error: Color combination ',ida1,ida2,
@@ -414,3 +331,79 @@
 
  999  write(*,*) 'error'
       end
+
+c*******************************************************************
+      function elim_indices(n3,n3bar,ncolmp,icolmp,icolres,maxcolor)
+c*******************************************************************
+c     Successively eliminate identical pairwise color indices from the
+c     icolmp list, until only the wanted indices remain
+c     n3 gives the number of triplet indices, n3bar number of antitriplets
+c     n3=1 for triplet, n3bar=1 for antitriplet, 
+c     (n3,n3bar)=(1,1) for octet,
+c     n3=2 for sextet, n3bar=2 for antisextet 
+c     If there are epsilon^{ijk} or epsilonbar color couplings, we
+c     need to introduce new index based on maxcolor.
+c
+
+      implicit none
+      integer elim_indices
+      integer n3,n3bar,ncolmp,icolmp(2,*),icolres(2),maxcolor
+      integer i,j,i3,i3bar
+
+c     Successively eliminate color indices in pairs until only the wanted
+c     indices remain
+      do i=1,ncolmp
+         do j=1,ncolmp
+            if(icolmp(1,i).ne.0.and.icolmp(1,i).eq.icolmp(2,j)) then
+               icolmp(1,i)=0
+               icolmp(2,j)=0
+            endif
+         enddo
+      enddo
+      
+      i3=0
+      i3bar=0
+      icolres(1)=0
+      icolres(2)=0
+      do i=1,ncolmp
+         if(icolmp(1,i).gt.0)then
+            i3=i3+1
+            if(i3.eq.1) icolres(1)=icolmp(1,i)
+            if(i3.eq.2) icolres(2)=-icolmp(1,i)
+         endif
+         if(icolmp(2,i).gt.0)then
+            i3bar=i3bar+1
+            if(i3bar.eq.1) icolres(2)=icolmp(2,i)
+            if(i3bar.eq.2) icolres(1)=-icolmp(2,i)
+         endif
+      enddo
+
+      if(i3.ne.n3.or.i3bar.ne.n3bar) then
+         if(n3.gt.0.and.n3bar.eq.0)then
+c        This is an epsilon index interaction
+            maxcolor=maxcolor+1
+            icolres(1)=maxcolor
+            icolres(2)=0
+            if(n3.eq.2)then
+               maxcolor=maxcolor+1
+               icolres(2)=-maxcolor
+            endif
+         elseif(n3bar.gt.0.and.n3.eq.0)then
+c        This is an epsilonbar index interaction
+            maxcolor=maxcolor+1
+            icolres(1)=0
+            icolres(2)=maxcolor
+            if(n3.eq.2)then
+               maxcolor=maxcolor+1
+               icolres(1)=-maxcolor
+            endif
+         else
+c           Don't know how to deal with this
+            call write_error(1001,n3,n3bar)
+         endif
+      endif
+
+      elim_indices=maxcolor
+      
+      return
+      end

=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2011-11-05 03:47:52 +0000
+++ UpdateNotes.txt	2011-11-14 05:57:24 +0000
@@ -1,5 +1,18 @@
 Update notes for MadGraph 5 (in reverse time order)
 
+1.3.28 (12/11/11) OM+JA: Fix special case when Lorentz structure combining
+                      two different Majorana particles depends on the
+                      incoming/outgoing status of the Majorana particles
+                      (needed for MSSM with Goldstino).
+                  OM: Update the parralel test (now testing against MG5 1.3.3)
+                  OM: Include some parralel test in the release script.
+                  JA: Fixed problem with colors in addmothers.f for complicated
+                      multiparticle vertices and simplified color treatment 
+                      (thanks to Gauthier Durieux for pointing this out).
+                  JA: Fixed problem with colors in addmothers.f for complicated
+                      multiparticle vertices and simplified color treatment 
+                  JA: Further improved gridpack parameters
+
 1.3.27 (05/11/11) JA: Fix bug in mirrored amplitudes (sometimes
                       amplitudes that should not be flagged as
                       mirrored were flagged as mirrored). Thanks
@@ -21,10 +34,10 @@
 1.3.25 (27/10/11) JA: Ensure that the correct intermediate resonance
                       is always written in the event file, even when we
                       have resonances with identical properties.
- 		  OM: Fix the bug forcing to quit the web browse in order to
-		      have MG5 continuing to run.
-		  OM: Change the tutorial in order to allow open index.html
-		      after the output command. 
+                  OM: Fix the bug forcing to quit the web browse in order to
+                      have MG5 continuing to run.
+                  OM: Change the tutorial in order to allow open index.html
+                      after the output command. 
  
 1.3.24 (22/10/11) JA: Fix problem with getting enough events in gridpack
                       mode (this was broken in v. 1.3.11 when we moved
@@ -80,7 +93,7 @@
                       new scale treatment (mu_F for pp>hbb~ is (pT_b^max*m_Th)),
                       no clustering of gluons to final-state massive particles
                       in MadEvent.
-		  JA: Ensure that factorization scale settings and matching works
+                  JA: Ensure that factorization scale settings and matching works
                       also in singlet t-channel exchange processes like
                       single top and VBF. The dynamic factorization
                       scale is given by the pT of the scattered quark
@@ -129,7 +142,7 @@
 1.3.12 (01/09/11) JA: Fixed problem with decay chains when different decays
                       result in identical final states, such as
                       p p > go go, (go > b1/b1~ b/b~, b1/b1~ > b/b~ n1)
-		      (only one of the decay chains was chosen, instead of
+                      (only one of the decay chains was chosen, instead of
                       all 3 combinations (b1,b1), (b1,b1~), (b1~,b1~))
                   JA: Allow for overall orders also with grouped subprocesses
                   JA: Ensure that only leading color flows are included in event
@@ -185,9 +198,9 @@
                       pythia8 process
 
 1.3.2 (14/06/11): OM: Fix fortran output when a model is case sensitive 
-      		      (Bug if a coupling was depending of a case sensitive parameter)
+                        (Bug if a coupling was depending of a case sensitive parameter)
                   SdV: Remove a annoying print in the new cuts (added in 1.3.0)
-		  OM: Fix a compilation problem in the standalone cpp output
+                  OM: Fix a compilation problem in the standalone cpp output
 
 1.3.1 (02/06/11): JA: Fixed missing file bug with the introduction of
                       inclusive HT cut
@@ -218,13 +231,13 @@
                       QED/QCD orders, e.g. p p > t t~ QED=2. Note that this fix
                       doesn't work for p p > t t~ j j QED=4 which should still
                       be avoided.
-		  JA: Added the ptllmin/max options in the default run_card.dat
+                  JA: Added the ptllmin/max options in the default run_card.dat
 
 1.2.2 (09/05/11): OM: fix ALOHA symmetries creating not gauge invariant result 
                       for scalar octet
 
 1.2.1 (08/05/11): OM: reduce the quantity of RAM use by matrix.f
-      	          OM: support speed of psyco if this python module is installed
+                  OM: support speed of psyco if this python module is installed
                   OM: fix a minor bug in the model parsing
                   OM: add the check of valid model.pkl also for v4 model
                   OM: add a check that UpdatesNotes is up-to-date when

=== modified file 'aloha/create_aloha.py'
--- aloha/create_aloha.py	2011-09-30 13:55:43 +0000
+++ aloha/create_aloha.py	2011-11-14 05:57:24 +0000
@@ -151,7 +151,9 @@
         if self.routine_kernel is None:
             self.kernel_tag = set()
             self.routine_kernel = eval(self.lorentz_expr)
-           
+        
+        # We need to compute C Gamma^T C^-1 = C_ab G_cb (-1) C_cd 
+        #                  = C_ac G_bc (-1) C_bd = C_ac G_bc C_db
         self.routine_kernel = \
              C(new_id, old_id + 1) * self.routine_kernel * C(new_id + 1, old_id)
         self.name += 'C'

=== modified file 'bin/create_release.py'
--- bin/create_release.py	2011-05-08 07:27:55 +0000
+++ bin/create_release.py	2011-11-14 05:57:24 +0000
@@ -199,6 +199,24 @@
     os.remove(filename)
     exit()
 
+p_test_results = test_manager.run(['test_short_.*'],
+                                  re_opt=0,
+                                  package=os.path.join('tests',
+                                                       'parallel_tests')
+                                  )
+
+if not p_test_results.wasSuccessful():
+    logging.error("Failed %d parallel tests, please check!" % \
+                  (len(p_test_results.errors) + len(p_test_results.failures)))
+
+if p_test_results.errors or p_test_results.failures:
+    logging.error("Removing %s and quitting..." % filename)
+    os.remove(filename)
+    exit()
+
+
+
+
 try:
     os.remove("%s.asc" % filename)
 except:

=== modified file 'madgraph/VERSION'
--- madgraph/VERSION	2011-11-05 01:40:58 +0000
+++ madgraph/VERSION	2011-11-14 05:57:24 +0000
@@ -1,2 +1,2 @@
-version = 1.3.27
-date = 2011-11-05
+version = 1.3.28
+date = 2011-11-12

=== modified file 'madgraph/core/helas_objects.py'
--- madgraph/core/helas_objects.py	2011-10-28 10:20:10 +0000
+++ madgraph/core/helas_objects.py	2011-11-14 05:57:24 +0000
@@ -222,6 +222,9 @@
         # should be onshell (True), as well as for forbidden s-channels (False).
         # Default is None
         self['onshell'] = None
+        # conjugate_indices is a list [1,2,...] with fermion lines
+        # that need conjugates. Default is "None"
+        self['conjugate_indices'] = None
 
     # Customized constructor
     def __init__(self, *arguments):
@@ -398,6 +401,12 @@
                         "%s is not a valid bool" % str(value) + \
                         " for onshell"
 
+        if name == 'conjugate_indices':
+            if not isinstance(value, tuple) and value != None:
+                raise self.PhysicsObjectError, \
+                        "%s is not a valid tuple" % str(value) + \
+                        " for conjugate_indices"
+
         return True
 
     # Enhanced get function, where we can directly call the properties of the particle
@@ -405,6 +414,10 @@
         """When calling any property related to the particle,
         automatically call the corresponding property of the particle."""
 
+        # Set conjugate_indices if it's not already set
+        if name == 'conjugate_indices' and self[name] == None:
+            self['conjugate_indices'] = self.get_conjugate_index()
+
         if name in ['spin', 'mass', 'width', 'self_antipart']:
             return self['particle'].get(name)
         elif name == 'pdg_code':
@@ -488,6 +501,9 @@
     def is_boson(self):
         return not self.is_fermion()
 
+    def is_majorana(self):
+        return self.is_fermion() and self.get('self_antipart')
+
     def to_array(self):
         """Generate an array with the information needed to uniquely
         determine if a wavefunction has been used before: interaction
@@ -841,9 +857,7 @@
         """Returns true if any of the mothers have negative
         fermionflow"""
 
-        return any([wf.get('fermionflow') < 0 for wf in \
-                    self.get('mothers')]) or \
-                    (self.get('interaction_id') and self.get('fermionflow') < 0)
+        return self.get('conjugate_indices') != ()
 
     def get_with_flow(self, name):
         """Generate the is_part and state needed for writing out
@@ -920,7 +934,7 @@
 
         # Check if we need to append a charge conjugation flag
         if self.needs_hermitian_conjugate():
-            res.append(self.get_conjugate_index())
+            res.append(self.get('conjugate_indices'))
 
         return (tuple(res), tuple(self.get('lorentz')))
 
@@ -1153,20 +1167,37 @@
     def get_conjugate_index(self):
         """Return the index of the particle that should be conjugated."""
 
-        if self.needs_hermitian_conjugate():
-            fermions = [wf for wf in self.get('mothers') if \
-                        wf.is_fermion()]
-            indices = []
-            self_index = self.find_outgoing_number() - 1
-            if self.is_fermion():
-                fermions.insert(self_index, self)
-            for i in range(0,len(fermions), 2):
-                if fermions[i].get('fermionflow') < 0 or \
-                   fermions[i+1].get('fermionflow') < 0:
-                    indices.append(i/2 + 1)
-            return tuple(indices)
-        else:
+        if not any([(wf.get('fermionflow') < 0 or wf.is_majorana()) for wf in \
+                    self.get('mothers')]) and \
+                    (not self.get('interaction_id') or \
+                    self.get('fermionflow') >= 0):
             return ()
+        
+        # Pick out first sorted mothers, then fermions
+        mothers, self_index = \
+                      self.get('mothers').sort_by_pdg_codes(self.get('pdg_codes'),
+                                                            self.get_anti_pdg_code())
+        fermions = HelasWavefunctionList([wf for wf in mothers if wf.is_fermion()])
+
+        # Insert this wavefunction in list (in the right place)
+        if self.is_fermion():
+            me = copy.copy(self)
+            # Flip incoming/outgoing to make me equivalent to mother
+            # as needed by majorana_conjugates
+            me.set('state', [state for state in ['incoming', 'outgoing'] \
+                             if state != me.get('state')][0])
+            fermions.insert(self_index, me)
+
+        # Initialize indices with indices due to Majoranas with wrong order
+        indices = fermions.majorana_conjugates()
+
+        # Check for fermions with negative fermion flow
+        for i in range(0,len(fermions), 2):
+            if fermions[i].get('fermionflow') < 0 or \
+               fermions[i+1].get('fermionflow') < 0:
+                indices.append(i/2 + 1)
+
+        return tuple(sorted(indices))
 
     def get_vertex_leg_numbers(self):
         """Get a list of the number of legs in vertices in this diagram"""
@@ -1405,6 +1436,40 @@
 
         return HelasWavefunctionList(sorted_mothers), my_index
 
+    def majorana_conjugates(self):
+        """Returns a list [1,2,...] of fermion lines that need
+         conjugate wfs due to wrong order of I/O Majorana particles
+         compared to interaction order (or empty list if no Majorana
+         particles).  This is crucial if the Lorentz structure depends
+         on the direction of the Majorana particles, as in MSSM with
+         goldstinos."""
+
+        if len([m for m in self if m.is_majorana()]) < 2:
+            return []
+
+        conjugates = []
+        
+        # Check if the order for Majorana fermions is correct
+        for i in range(0, len(self), 2):
+            if self[i].is_majorana() and self[i+1].is_majorana() \
+                   and self[i].get_pdg_code() != \
+                   self[i+1].get_pdg_code():
+                # Check if mother I/O order is correct (IO)
+                if self[i].get_spin_state_number() > 0 and \
+                   self[i + 1].get_spin_state_number() < 0:
+                    # Order is wrong, we need a conjugate here
+                    conjugates.append(True)
+                else:
+                    conjugates.append(False)
+            elif self[i].is_fermion():
+                # For non-Majorana case, always False
+                conjugates.append(False)
+
+        # Return list 1,2,... for which indices are needed
+        conjugates = [i+1 for (i,c) in enumerate(conjugates) if c]
+
+        return conjugates
+
     @staticmethod
     def extract_wavefunctions(mothers):
         """Recursively extract the wavefunctions from mothers of mothers"""
@@ -1442,6 +1507,9 @@
         self['fermionfactor'] = 0
         self['color_indices'] = []
         self['mothers'] = HelasWavefunctionList()
+        # conjugate_indices is a list [1,2,...] with fermion lines
+        # that need conjugates. Default is "None"
+        self['conjugate_indices'] = None
 
     # Customized constructor
     def __init__(self, *arguments):
@@ -1559,6 +1627,12 @@
                       "%s is not a valid list of mothers for amplitude" % \
                       str(value)
 
+        if name == 'conjugate_indices':
+            if not isinstance(value, tuple) and value != None:
+                raise self.PhysicsObjectError, \
+                        "%s is not a valid tuple" % str(value) + \
+                        " for conjugate_indices"
+
         return True
 
     def __str__(self):
@@ -1592,6 +1666,10 @@
         if name == 'fermionfactor' and not self[name]:
             self.calculate_fermionfactor()
 
+        # Set conjugate_indices if it's not already set
+        if name == 'conjugate_indices' and self[name] == None:
+            self['conjugate_indices'] = self.get_conjugate_index()
+
         return super(HelasAmplitude, self).get(name)
 
     # Enhanced set function, where we can append a model
@@ -1667,8 +1745,7 @@
         """Returns true if any of the mothers have negative
         fermionflow"""
 
-        return any([wf.get('fermionflow') < 0 for wf in \
-                    self.get('mothers')])
+        return self.get('conjugate_indices') != ()
 
     def get_call_key(self):
         """Generate the (spin, state) tuples used as key for the helas call
@@ -1683,7 +1760,7 @@
 
         # Check if we need to append a charge conjugation flag
         if self.needs_hermitian_conjugate():
-            res.append(self.get_conjugate_index())
+            res.append(self.get('conjugate_indices'))
 
         return (tuple(res), tuple(self.get('lorentz')))
 
@@ -1728,9 +1805,7 @@
         for i in range(len(fermions) - 1):
             for j in range(i + 1, len(fermions)):
                 if fermions[j] < fermions[i]:
-                    tmp = fermions[i]
-                    fermions[i] = fermions[j]
-                    fermions[j] = tmp
+                    fermions[i], fermions[j] = fermions[j], fermions[i]
                     nflips = nflips + 1
 
         return (-1) ** nflips
@@ -1870,12 +1945,17 @@
             schannels.extend(mother_s)
             tchannels.extend(mother_t)
 
+        # Sort s-channels according to number
+        schannels.sort(lambda x1,x2: x2.get('legs')[-1].get('number') - \
+                       x1.get('legs')[-1].get('number'))
+
         # Split up multiparticle vertices using fake s-channel propagators
         multischannels = [(i, v) for (i, v) in enumerate(schannels) \
                           if len(v.get('legs')) > 3]
         multitchannels = [(i, v) for (i, v) in enumerate(tchannels) \
                           if len(v.get('legs')) > 3]
-
+        
+        increase = 0
         for channel in multischannels + multitchannels:
             newschannels = []
             vertex = channel[1]
@@ -1897,7 +1977,9 @@
 
                 # Insert the new s-channel before this vertex
                 if channel in multischannels:
-                    schannels.insert(channel[0], new_vertex)
+                    schannels.insert(channel[0]+increase, new_vertex)
+                    # Account for previous insertions
+                    increase += 1
                 else:
                     schannels.append(new_vertex)
                 legs = vertex.get('legs')
@@ -1906,10 +1988,6 @@
                 # Renumber resulting leg according to minimum leg number
                 legs[-1].set('number', min([l.get('number') for l in legs[:-1]]))
 
-        # Sort s-channels according to number
-        schannels.sort(lambda x1,x2: x2.get('legs')[-1].get('number') - \
-                       x1.get('legs')[-1].get('number'))
-
         # Finally go through all vertices, sort the legs and replace
         # leg number with propagator number -1, -2, ...
         number_dict = {}
@@ -1966,17 +2044,25 @@
     def get_conjugate_index(self):
         """Return the index of the particle that should be conjugated."""
 
-        if self.needs_hermitian_conjugate():
-            fermions = [wf for wf in self.get('mothers') if \
-                        wf.is_fermion()]
-            indices = []
-            for i in range(0,len(fermions), 2):
-                if fermions[i].get('fermionflow') < 0 or \
-                   fermions[i+1].get('fermionflow') < 0:
-                    indices.append(i/2 + 1)
-            return tuple(indices)
-        else:
+        if not any([(wf.get('fermionflow') < 0 or wf.is_majorana()) for wf in \
+                    self.get('mothers')]):
             return ()
+        
+        # Pick out first sorted mothers, then fermions
+        mothers, self_index = \
+                      self.get('mothers').sort_by_pdg_codes(self.get('pdg_codes'))
+        fermions = HelasWavefunctionList([wf for wf in mothers if wf.is_fermion()])
+
+        # Initialize indices with indices due to Majoranas with wrong order
+        indices = fermions.majorana_conjugates()
+
+        # Check for fermions with negative fermion flow
+        for i in range(0,len(fermions), 2):
+            if fermions[i].get('fermionflow') < 0 or \
+               fermions[i+1].get('fermionflow') < 0:
+                indices.append(i/2 + 1)
+                
+        return tuple(sorted(indices))
 
     def get_vertex_leg_numbers(self):
         """Get a list of the number of legs in vertices in this diagram"""

=== modified file 'models/__init__.py'
--- models/__init__.py	2010-09-23 12:41:20 +0000
+++ models/__init__.py	2011-11-14 05:57:24 +0000
@@ -29,7 +29,7 @@
         __import__(model_pos)
         return sys.modules[model_pos]
     else:
-        sys.path.append(os.sep.join(path_split[:-1]))
+        sys.path.insert(0, os.sep.join(path_split[:-1]))
         __import__(path_split[-1])
+        sys.path.pop(0)
         return sys.modules[path_split[-1]]
-        sys.path.pop()

=== modified file 'models/import_ufo.py'
--- models/import_ufo.py	2011-07-22 13:55:20 +0000
+++ models/import_ufo.py	2011-11-14 05:57:24 +0000
@@ -84,7 +84,7 @@
             restrict_file = os.path.join(model_path,'restrict_default.dat')
         else:
             restrict_file = None
-            
+    
     #import the FULL model
     model = import_full_model(model_path) 
     # restore the model name

=== added file 'tests/parallel_tests/compare_with_old_mg5_version.py'
--- tests/parallel_tests/compare_with_old_mg5_version.py	1970-01-01 00:00:00 +0000
+++ tests/parallel_tests/compare_with_old_mg5_version.py	2011-11-14 05:57:24 +0000
@@ -0,0 +1,451 @@
+################################################################################
+#
+# Copyright (c) 2009 The MadGraph Development team and Contributors
+#
+# This file is a part of the MadGraph 5 project, an application which 
+# automatically generates Feynman diagrams and matrix elements for arbitrary
+# high-energy processes in the Standard Model and beyond.
+#
+# It is subject to the MadGraph license which should accompany this 
+# distribution.
+#
+# For more information, please visit: http://madgraph.phys.ucl.ac.be
+#
+################################################################################
+""" A test suite to compare the current version of MG5 with a version of 
+reference. The model is taken from the reference version. So the only 
+constraints to this test is that the current version can still read the old 
+model format. 
+The reference version is given here as a argument which can be changed by hand.
+"""
+import logging
+import os
+import shutil
+import me_comparator
+import unittest
+import subprocess
+from madgraph import MG5DIR
+
+pjoin = os.path.join
+_file_path = os.path.dirname(os.path.realpath(__file__))
+_pickle_path = pjoin(_file_path, 'input_files')
+
+
+class OLDMG5Comparator(unittest.TestCase):
+    """A class to compare the value of a old MG5 version and the current one"""
+    
+    old_mg5 = None
+    reference_number = 149 #149 corresponds to 1.3.3 
+    nb_test = 0
+    
+    
+    @classmethod
+    def build_old_mg5(cls):
+        """build the directory for mg5 run in the old template"""
+        if cls.old_mg5:
+            return cls.old_mg5
+        
+        init_dir = os.getcwd()
+        os.chdir(MG5DIR)    
+        # 1. bzr branch the present directory to a new directory
+        #    parralel_tests_mg5
+
+        filepath = 'parralel_tests_mg5'
+        if os.path.exists(filepath):
+            logging.info("Removing existing directory " + filepath)
+            shutil.rmtree(filepath)
+        logging.info("Branching " + MG5DIR + " to directory " + filepath)
+        devnull = open(os.devnull,'w')   
+        status = subprocess.call(['bzr', 'branch', MG5DIR, filepath], stdout=devnull, stderr=devnull)   
+        if status:
+            raise Exception, 'Impossible to configure old run'
+        
+        #2. reverse the directory to its old status
+        os.chdir(filepath)
+        status = subprocess.call(['bzr', 'revert', '-r', str(cls.reference_number)], stdout=devnull, stderr=devnull)
+        if status:
+            raise Exception, 'Impossible to configure old run'
+        
+        cls.old_mg5 = pjoin(MG5DIR,filepath)
+        os.chdir(init_dir)
+        return cls.old_mg5 
+
+    def compare_processes(self, my_proc_list = [], orders = {}, model = 'sm',
+                        energy = 500, filename = "", pickle_file = "",
+                        tolerance = 1e-06):
+        """ """
+        mg5_path = self.build_old_mg5()
+        
+        if 'v4' in model:
+            old_mg5 = me_comparator.MG5OldRunner()
+            old_mg5.setup(mg5_path)
+            current_mg5 = me_comparator.MG5Runner()
+            current_mg5.setup(MG5DIR, MG5DIR)
+            current_mg5.store_proc_card = True
+        else:
+            old_mg5 = me_comparator.MG5_UFO_OldRunner()
+            old_mg5.setup(mg5_path)
+            current_mg5 = me_comparator.MG5_UFO_Runner()
+            current_mg5.setup(MG5DIR,MG5DIR)
+            current_mg5.store_proc_card = True
+        
+        self.nb_test +=1      
+        while os.path.exists(pjoin(MG5DIR,'models','paralel_test_model_%s' % self.nb_test )):
+            self.nb_test += 1
+        os.system('cp -rf %s %s' % (pjoin(mg5_path,'models',model) ,
+                                    pjoin(MG5DIR,'models','paralel_test_model_%s' % self.nb_test)))
+        
+        # Create and setup a comparator
+        my_comp = me_comparator.MEComparator()
+        my_comp.set_me_runners(current_mg5, old_mg5)
+
+        # Run the actual comparison
+        my_comp.run_comparison(my_proc_list,
+                               ['paralel_test_model_%s' %self.nb_test, model], orders, energy)
+
+        # Print the output
+        if filename:
+            my_comp.output_result(filename=filename)
+        
+                # Store output to a pickle file in the input_files directory
+        if pickle_file:
+            me_comparator.PickleRunner.store_comparison(\
+                os.path.join(_pickle_path, pickle_file),
+                my_comp.get_non_zero_processes(),
+                my_comp.me_runners[0].model,
+                my_comp.me_runners[0].name,
+                my_comp.me_runners[0].orders,
+                my_comp.me_runners[0].energy)
+
+        # Assert that all process comparisons passed the tolerance cut
+        my_comp.assert_processes(self, tolerance)
+            
+        # Do some cleanup
+        my_comp.cleanup()
+       
+    ############################################################################    
+    #  ROUTINE FOR CREATING THE SHORT TEST (USE by the release script)
+    ############################################################################    
+    def test_create_all_pickle(self):
+        """re-create all the pickle for the short test (those call in the release).
+           Note that if you need to redo this, this is potentially due to a change
+           in the model. In consequence, you need to change the old MG5 comparison
+           point. (Since the default use another model)."""
+        
+        return # By default no need this
+        
+        self.create_short_paralel_sm()
+        self.create_short_paralel_mssm()
+        self.create_short_paralel_heft()
+        
+        
+
+    def create_short_paralel_sm(self):
+        """Test a short list of sm processes"""
+        # Create a list of processes to check automatically
+        my_proc_list = ['e+ e- > e+ e-', 
+                        'h h > h h', 
+                        'g g > t t~',
+                        'w+ w- > w+ w-',
+                        'b b~ > t t~',
+                        'u u~ > z u u~',
+                        'g g > t t~ h',
+                        'u u~ > d d~ w+ w-']
+        
+        # Store list of non-zero processes and results in file
+        pickle_file = os.path.join(_pickle_path, "mg5_short_parraleltest_sm.pkl")
+        self.compare_processes(my_proc_list,
+                             orders = {'QED':99, 'QCD':99},
+                             filename = "short_sm.log",
+                             pickle_file = pickle_file)
+
+    def create_short_paralel_mssm(self):
+        """Test a short list of mssm processes"""
+        # Create a list of processes to check automatically
+        my_proc_list = [' g g > go go', 
+                        'u u~ > go go', 
+                        'e+ e- > n1 n2',
+                        'd d~ > el+ el-',
+                        'b b~ > h1 h2',
+                        'u u~ > z u u~',
+                        'e+ e- > n1 n2 z',
+                        'd d~ > x1+ x1- g',
+                        'h1 h1 > x1+ x1- h2']
+        
+        # Store list of non-zero processes and results in file
+        pickle_file = os.path.join(_pickle_path, "mg5_short_parraleltest_mssm.pkl")
+        self.compare_processes(my_proc_list,
+                             model='mssm',
+                             orders = {'QED':99, 'QCD':99},
+                             filename = "short_mssm.log",
+                             pickle_file = pickle_file)
+
+    def create_short_paralel_heft(self):
+        """Test a short list of mssm processes"""
+        # Create a list of processes to check automatically
+        my_proc_list = ['h h > w+ w-', 
+                        'g g > h g', 
+                        'g g > h g g',
+                        'g g > h g g g']
+        
+        # Store list of non-zero processes and results in file
+        pickle_file = os.path.join(_pickle_path, "mg5_short_parraleltest_heft.pkl")
+        self.compare_processes(my_proc_list,
+                             model='heft',
+                             orders = {'QED':99, 'QCD':99},
+                             filename = "short_heft.log",
+                             pickle_file = pickle_file)
+
+    ############################################################################    
+    #  ROUTINE FOR THE SHORT TEST (USE by the release script)
+    ############################################################################
+    def test_short_sm(self):
+        """Test a minimal list of sm 2->2 processes, mainly to test the test"""
+        # Create a list of processes to check automatically
+        comparisons = me_comparator.PickleRunner.find_comparisons(\
+            os.path.join(_pickle_path, "mg5_short_parraleltest_sm.pkl"))
+        for stored_runner in comparisons:
+            # Create a MERunner object for MG5
+            my_mg5 = me_comparator.MG5_UFO_Runner()
+            my_mg5.setup(MG5DIR, MG5DIR)
+
+            # Create and setup a comparator
+            my_comp = me_comparator.MEComparator()
+            my_comp.set_me_runners(stored_runner, my_mg5)
+
+            # Run the actual comparison
+            my_comp.run_comparison(stored_runner.proc_list,
+                                   'sm',
+                                   stored_runner.orders,
+                                   stored_runner.energy)
+
+            my_comp.assert_processes(self)
+
+            # Do some cleanup
+            my_comp.cleanup()
+            
+    def test_short_mssm(self):
+        """Test a minimal list of sm 2->2 processes, mainly to test the test"""
+        # Create a list of processes to check automatically
+        comparisons = me_comparator.PickleRunner.find_comparisons(\
+            os.path.join(_pickle_path, "mg5_short_parraleltest_mssm.pkl"))
+        for stored_runner in comparisons:
+            # Create a MERunner object for MG5
+            my_mg5 = me_comparator.MG5_UFO_Runner()
+            my_mg5.setup(MG5DIR, MG5DIR)
+
+            # Create and setup a comparator
+            my_comp = me_comparator.MEComparator()
+            my_comp.set_me_runners(stored_runner, my_mg5)
+
+            # Run the actual comparison
+            my_comp.run_comparison(stored_runner.proc_list,
+                                   'mssm',
+                                   stored_runner.orders,
+                                   stored_runner.energy)
+
+            my_comp.assert_processes(self)
+
+            # Do some cleanup
+            my_comp.cleanup()           
+
+    def test_short_heft(self):
+        """Test a minimal list of sm 2->2 processes, mainly to test the test"""
+        # Create a list of processes to check automatically
+        comparisons = me_comparator.PickleRunner.find_comparisons(\
+            os.path.join(_pickle_path, "mg5_short_parraleltest_heft.pkl"))
+        for stored_runner in comparisons:
+            # Create a MERunner object for MG5
+            my_mg5 = me_comparator.MG5_UFO_Runner()
+            my_mg5.setup(MG5DIR, MG5DIR)
+
+            # Create and setup a comparator
+            my_comp = me_comparator.MEComparator()
+            my_comp.set_me_runners(stored_runner, my_mg5)
+
+            # Run the actual comparison
+            my_comp.run_comparison(stored_runner.proc_list,
+                                   'heft',
+                                   stored_runner.orders,
+                                   stored_runner.energy)
+
+            my_comp.assert_processes(self)
+
+            # Do some cleanup
+            my_comp.cleanup()   
+        
+    ############################################################################    
+    #  ROUTINE FOR CHECKING THE PARRALEL TEST
+    ############################################################################           
+    def test_mg5_minitest_sm(self):
+        """Test a minimal list of sm 2->2 processes, mainly to test the test"""
+        # Create a list of processes to check automatically
+        my_proc_list = me_comparator.create_proc_list(\
+            ['u'],
+            initial=2, final=2)
+        my_proc_list = ['e+ e- > a > e+ e-', 'h h > h h', 'e+ e+ > e- e-']
+        # Store list of non-zero processes and results in file
+        #pickle_file = "mg4_sm_%sminitest.pkl" % self.suffix_name
+        self.compare_processes(my_proc_list, model='sm',
+                             orders = {'QED':2, 'QCD':2},
+                             filename = "sm_mini.log",
+                             energy = 1000)
+        
+    def test_mg5_minitest_mssm(self):
+        """Test a minimal list of sm 2->2 processes, mainly to test the test"""
+        # Create a list of processes to check automatically
+        my_proc_list = me_comparator.create_proc_list(\
+            ['u'],
+            initial=2, final=2)
+        my_proc_list = ['g g > go go', 'e+ e-  > n1 n2']
+        # Store list of non-zero processes and results in file
+        #pickle_file = "mg4_sm_%sminitest.pkl" % self.suffix_name
+        self.compare_processes(my_proc_list, model='mssm',
+                             orders = {'QED':2, 'QCD':2},
+                             filename = "mssm_mini.log",
+                             energy = 1000)
+
+    ############################################################################    
+    #  EXTENSIVE TEST FOR THE SM
+    ############################################################################ 
+    def test_mg5_sm_22(self):
+        """Test a semi-complete list of sm 2->2 processes"""
+        # Create a list of processes to check automatically
+        my_proc_list = me_comparator.create_proc_list(\
+            ['w+', 'w-', 'a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
+            'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
+            initial=2, final=2)
+
+        # Store list of non-zero processes and results in file
+        self.compare_processes(my_proc_list,
+                             orders = {'QED':2, 'QCD':2},
+                             filename = "sm_22.log")
+
+
+    def test_mg5_sm_23(self):
+        """Test a semi-complete list of sm 2->3 processes"""
+        # Create a list of processes to check automatically
+        my_proc_list = me_comparator.create_proc_list(\
+            ['w+', 'w-','a', 'z', 'h', 'g', 'u', 'u~', 'd', 'd~',
+             'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~'],
+            initial=2, final=3)
+
+        for i in range(len(my_proc_list)//500):
+            print 'step %s/%s' %(i+1,len(my_proc_list)//500 )
+        # Store list of non-zero processes and results in file
+            self.compare_processes(my_proc_list[500*i:500*(i+1)],
+                             orders = {'QED':3, 'QCD':3},
+                             filename = "sm_23.log")
+
+
+
+    ############################################################################    
+    #  EXTENSIVE TEST FOR THE MSSM
+    ############################################################################ 
+    def test_mg5_mssm_22(self):
+        """Test a semi-complete list of mssm 2->2 processes"""
+        # Create a list of processes to check automatically
+        sm_parts = ['w+', 'w-', 'a', 'z', 'h1', 'h+', 'h-', 'g', 'u', 'u~',
+            'd', 'd~', 'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~']
+        mssm_parts = ['dl', 'dl~', 'dr', 'dr~', 'ul', 'ul~', 'ur', 'ur~', 'b1',
+                      'b1~', 'b2', 'b2~', 't1', 't1~', 'ta1-', 'ta1+', 'ta2-',
+                      'ta2+', 'svt', 'svt~', 'x1-', 'x1+', 'x2-', 'x2+',
+                      'go', 'n1']
+        # Generate 2 -> 2 processes, with MSSM particles in pairs in
+        # final state
+        my_proc_list = me_comparator.create_proc_list_enhanced(\
+            sm_parts, mssm_parts)
+
+        # Store list of non-zero processes and results in file
+        self.compare_processes(my_proc_list,
+                             orders = {'QED':2, 'QCD':2},
+                             model = "mssm",
+                             energy = 2000,
+                             filename = "mssm_22.log")
+
+    def test_mg5_mssm_23_p1(self):
+        """Test a semi-complete list of mssm 2->3 processes"""
+        # Create a list of processes to check automatically
+        
+        init_part_list1 = ['w+', 'a', 'z', 'h1', 'h+', 'g', 'u~',
+            'd~', 'b~', 't~', 'ta+', 'vt~']
+        init_part_list2 = ['w-', 'a', 'z', 'h1', 'h-', 'g', 'u',
+            'd', 'b', 't', 'ta+', 'vt']        
+        final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g'] 
+        final_part_list2 = ['dl~', 'dr~', 'ul~', 'ur~']
+        final_part_list3 = ['dl', 'dr', 'ul', 'ur']
+        
+        my_proc_list = me_comparator.create_proc_list_2_3(
+                    init_part_list1, init_part_list2, final_part_list1,
+                      final_part_list2,final_part_list3)                                   
+                                                                  
+        sm_parts = ['w+', 'w-', 'a', 'z', 'h1', 'h+', 'h-', 'g', 'u', 'u~',
+            'd', 'd~', 'b', 'b~', 't', 't~', 'ta+', 'ta-', 'vt', 'vt~']
+        mssm_parts = ['dl', 'dl~', 'dr', 'dr~', 'ul', 'ul~', 'ur', 'ur~', 'b1',
+                      'b1~', 'b2', 'b2~', 't1', 't1~', 'ta1-', 'ta1+', 'ta2-',
+                      'ta2+', 'svt', 'svt~', 'x1-', 'x1+', 'x2-', 'x2+',
+                      'go', 'n1']
+        # Generate 2 -> 2+1 processes, with MSSM particles in pairs in
+        # final state
+        #my_proc_list = me_comparator.create_proc_list_enhanced(\
+        #    sm_parts, mssm_parts, sm_parts)
+        for i in range(len(my_proc_list)//500):
+            print 'step %s/%s' %(i+1,len(my_proc_list)//500 )
+        # Store list of non-zero processes and results in file
+            self.compare_processes(my_proc_list[500*i:500*(i+1)],
+                             orders = {'QED':3, 'QCD':3},
+                             model = "mssm",
+                             energy = 2000,
+                             filename = "mssm_23_p1.log")    
+
+    def test_mg5_mssm_23_p2(self):
+        """Test a semi-complete list of mssm 2->3 processes"""
+        # Create a list of processes to check automatically
+        
+        init_part_list1 = ['w+', 'a', 'z', 'h1', 'h+', 'g', 'u~',
+            'd~', 'b~', 't~', 'ta+', 'vt~']
+        init_part_list2 = ['w-', 'a', 'z', 'h1', 'h-', 'g', 'u',
+            'd', 'b', 't', 'ta+', 'vt']        
+        final_part_list1 = ['w-', 'w+', 'a', 'z', 'h1', 'h-', 'h+', 'g'] 
+        final_part_list2 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']
+        final_part_list3 = ['x1-', 'x1+', 'x2-', 'x2+', 'go', 'n1']
+        
+        my_proc_list = me_comparator.create_proc_list_2_3(
+                    init_part_list1, init_part_list2, final_part_list1,
+                      final_part_list2,final_part_list3)
+
+        for i in range(len(my_proc_list)//500):
+            print 'step %s/%s' %(i+1,len(my_proc_list)//500 )
+        # Store list of non-zero processes and results in file
+            self.compare_processes(my_proc_list[500*i:500*(i+1)],
+                             orders = {'QED':3, 'QCD':3},
+                             model = "mssm",
+                             energy = 2000,
+                             filename = "mssm_23_p2.log")   
+
+        
+    ############################################################################    
+    #  EXTENSIVE TEST FOR THE HEFT
+    ############################################################################     
+    def test_mg5_heft_23(self):
+        """Test a heft 2->3 processes"""
+        # Create a list of processes to check automatically
+        sm_parts = ['g', 'a', 'w+', 'w-']
+        heft_parts = ['h']
+        # Generate 2 -> 1+2 processes, with one Higgs in final state
+        # final state
+        my_proc_list = me_comparator.create_proc_list_enhanced(\
+            sm_parts, sm_parts, heft_parts)
+
+        # Store list of non-zero processes and results in file
+        self.compare_MG4_MG5(my_proc_list,
+                             orders = {'QED':2, 'QCD':0, 'HIG':1, 'HIW':1},
+                             model = "heft",
+                             energy = 500,
+                             filename = "heft_23.log")
+
+
+
+
+
+        
\ No newline at end of file

=== removed file 'tests/parallel_tests/input_files/mg4_heft_23.pkl'
--- tests/parallel_tests/input_files/mg4_heft_23.pkl	2010-12-18 22:39:25 +0000
+++ tests/parallel_tests/input_files/mg4_heft_23.pkl	1970-01-01 00:00:00 +0000
@@ -1,493 +0,0 @@
-ccopy_reg
-_reconstructor
-p0
-(ctests.parallel_tests.me_comparator
-PickleRunner
-p1
-c__builtin__
-object
-p2
-Ntp3
-Rp4
-(dp5
-S'res_list'
-p6
-(lp7
-((F4.552249565111725e-16
-I-2
-tp8
-(lp9
-(lp10
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp11
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp12
-F199.2265
-aF73.61784
-aF164.938
-aF-84.06588
-aa(lp13
-F158.1664
-aF-7.963291
-aF-151.0579
-aF46.20289
-aa(lp14
-F142.6071
-aF-65.65455
-aF-13.88019
-aF37.86299
-aatp15
-a((F5.157051974955813e-10
-I-2
-tp16
-(lp17
-(lp18
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp19
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp20
-F222.571
-aF69.26654
-aF155.1891
-aF-79.09703
-aa(lp21
-F169.1565
-aF-7.492608
-aF-142.1293
-aF43.47199
-aa(lp22
-F108.2725
-aF-61.77393
-aF-13.05978
-aF35.62504
-aatp23
-a((F3.183580321665399e-15
-I-2
-tp24
-(lp25
-(lp26
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp27
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp28
-F199.2265
-aF73.61784
-aF164.938
-aF-84.06588
-aa(lp29
-F158.1664
-aF-7.963291
-aF-151.0579
-aF46.20289
-aa(lp30
-F142.6071
-aF-65.65455
-aF-13.88019
-aF37.86299
-aatp31
-a((F5.839526791487125e-10
-I-2
-tp32
-(lp33
-(lp34
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp35
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp36
-F194.4067
-aF71.83682
-aF160.9477
-aF-82.0321
-aa(lp37
-F195.5014
-aF-7.770638
-aF-147.4034
-aF45.08512
-aa(lp38
-F110.0919
-aF-64.06619
-aF-13.54439
-aF36.94698
-aatp39
-a((F5.839526791487125e-10
-I-2
-tp40
-(lp41
-(lp42
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp43
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp44
-F194.4067
-aF71.83682
-aF160.9477
-aF-82.0321
-aa(lp45
-F195.5014
-aF-7.770638
-aF-147.4034
-aF45.08512
-aa(lp46
-F110.0919
-aF-64.06619
-aF-13.54439
-aF36.94698
-aatp47
-a((F2.913439721671504e-14
-I-2
-tp48
-(lp49
-(lp50
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp51
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp52
-F199.2265
-aF73.61784
-aF164.938
-aF-84.06588
-aa(lp53
-F158.1664
-aF-7.963291
-aF-151.0579
-aF46.20289
-aa(lp54
-F142.6071
-aF-65.65455
-aF-13.88019
-aF37.86299
-aatp55
-a((F1.2205648499991544e-08
-I-2
-tp56
-(lp57
-(lp58
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp59
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp60
-F222.571
-aF69.26654
-aF155.1891
-aF-79.09703
-aa(lp61
-F169.1565
-aF-7.492608
-aF-142.1293
-aF43.47199
-aa(lp62
-F108.2725
-aF-61.77393
-aF-13.05978
-aF35.62504
-aatp63
-a((F8.181054310911128e-10
-I-2
-tp64
-(lp65
-(lp66
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp67
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp68
-F194.4067
-aF71.83682
-aF160.9477
-aF-82.0321
-aa(lp69
-F195.5014
-aF-7.770638
-aF-147.4034
-aF45.08512
-aa(lp70
-F110.0919
-aF-64.06619
-aF-13.54439
-aF36.94698
-aatp71
-a((F8.181054310911127e-10
-I-2
-tp72
-(lp73
-(lp74
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp75
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp76
-F194.4067
-aF71.83682
-aF160.9477
-aF-82.0321
-aa(lp77
-F195.5014
-aF-7.770638
-aF-147.4034
-aF45.08512
-aa(lp78
-F110.0919
-aF-64.06619
-aF-13.54439
-aF36.94698
-aatp79
-a((F7.521239868390367e-09
-I-2
-tp80
-(lp81
-(lp82
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp83
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp84
-F222.571
-aF69.26654
-aF155.1891
-aF-79.09703
-aa(lp85
-F169.1565
-aF-7.492608
-aF-142.1293
-aF43.47199
-aa(lp86
-F108.2725
-aF-61.77393
-aF-13.05978
-aF35.62504
-aatp87
-a((F1.2890895736403965e-08
-I-2
-tp88
-(lp89
-(lp90
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp91
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp92
-F199.2265
-aF73.61784
-aF164.938
-aF-84.06588
-aa(lp93
-F158.1664
-aF-7.963291
-aF-151.0579
-aF46.20289
-aa(lp94
-F142.6071
-aF-65.65455
-aF-13.88019
-aF37.86299
-aatp95
-a((F9.44746991279531e-10
-I-2
-tp96
-(lp97
-(lp98
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp99
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp100
-F199.2265
-aF73.61784
-aF164.938
-aF-84.06588
-aa(lp101
-F158.1664
-aF-7.963291
-aF-151.0579
-aF46.20289
-aa(lp102
-F142.6071
-aF-65.65455
-aF-13.88019
-aF37.86299
-aatp103
-a((F8.117582898346285e-09
-I-2
-tp104
-(lp105
-(lp106
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp107
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp108
-F222.571
-aF69.26654
-aF155.1891
-aF-79.09703
-aa(lp109
-F169.1565
-aF-7.492608
-aF-142.1293
-aF43.47199
-aa(lp110
-F108.2725
-aF-61.77393
-aF-13.05978
-aF35.62504
-aatp111
-a((F7.52123986839037e-09
-I-2
-tp112
-(lp113
-(lp114
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp115
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp116
-F222.571
-aF69.26654
-aF155.1891
-aF-79.09703
-aa(lp117
-F169.1565
-aF-7.492608
-aF-142.1293
-aF43.47199
-aa(lp118
-F108.2725
-aF-61.77393
-aF-13.05978
-aF35.62504
-aatp119
-asS'name'
-p120
-S'Stored MadGraph v4'
-p121
-sS'energy'
-p122
-I500
-sS'proc_list'
-p123
-(lp124
-S'g g > a a h'
-p125
-aS'g g > h w+ w-'
-p126
-aS'a g > a g h'
-p127
-aS'g w+ > g h w+'
-p128
-aS'g w- > g h w-'
-p129
-aS'a a > g g h'
-p130
-aS'a a > h w+ w-'
-p131
-aS'a w+ > a h w+'
-p132
-aS'a w- > a h w-'
-p133
-aS'w+ w+ > h w+ w+'
-p134
-aS'w+ w- > g g h'
-p135
-aS'w+ w- > a a h'
-p136
-aS'w+ w- > h w+ w-'
-p137
-aS'w- w- > h w- w-'
-p138
-asS'model'
-p139
-S'heft'
-p140
-sS'orders'
-p141
-(dp142
-S'HIW'
-p143
-I1
-sS'HIG'
-p144
-I1
-sS'QED'
-p145
-I2
-sS'QCD'
-p146
-I0
-ssb.
\ No newline at end of file

=== removed file 'tests/parallel_tests/input_files/mg4_sm_22.pkl'
--- tests/parallel_tests/input_files/mg4_sm_22.pkl	2010-12-18 22:39:25 +0000
+++ tests/parallel_tests/input_files/mg4_sm_22.pkl	1970-01-01 00:00:00 +0000
@@ -1,12108 +0,0 @@
-ccopy_reg
-_reconstructor
-p0
-(ctests.parallel_tests.me_comparator
-PickleRunner
-p1
-c__builtin__
-object
-p2
-Ntp3
-Rp4
-(dp5
-S'res_list'
-p6
-(lp7
-((F2.607967194937274
-I0
-tp8
-(lp9
-(lp10
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp11
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp12
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp13
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp14
-a((F0.44485271984270963
-I0
-tp15
-(lp16
-(lp17
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp18
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp19
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp20
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp21
-a((F0.10891452740443075
-I0
-tp22
-(lp23
-(lp24
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp25
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp26
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp27
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp28
-a((F0.8251799073331111
-I0
-tp29
-(lp30
-(lp31
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp32
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp33
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp34
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp35
-a((F0.1162498764206988
-I0
-tp36
-(lp37
-(lp38
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp39
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp40
-F235.6
-aF52.26752
-aF209.6043
-aF-94.02934
-aa(lp41
-F264.4
-aF-52.26752
-aF-209.6043
-aF94.02934
-aatp42
-a((F1.5136546830891073
-I0
-tp43
-(lp44
-(lp45
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp46
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp47
-F250.0
-aF51.64106
-aF207.092
-aF-92.90234
-aa(lp48
-F250.0
-aF-51.64106
-aF-207.092
-aF92.90234
-aatp49
-a((F0.3060094780480738
-I0
-tp50
-(lp51
-(lp52
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp53
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp54
-F256.0847
-aF50.1885
-aF201.2669
-aF-90.28918
-aa(lp55
-F243.9153
-aF-50.1885
-aF-201.2669
-aF90.28918
-aatp56
-a((F0.025237530351830174
-I0
-tp57
-(lp58
-(lp59
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp60
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp61
-F250.0
-aF48.65519
-aF195.118
-aF-87.53075
-aa(lp62
-F250.0
-aF-48.65519
-aF-195.118
-aF87.53075
-aatp63
-a((F0.023894097966856313
-I0
-tp64
-(lp65
-(lp66
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp67
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp68
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp69
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp70
-a((F0.10226683540928873
-I0
-tp71
-(lp72
-(lp73
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp74
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp75
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp76
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp77
-a((F0.5452853056572509
-I0
-tp78
-(lp79
-(lp80
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp81
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp82
-F250.0
-aF55.45234
-aF222.3761
-aF-99.75883
-aa(lp83
-F250.0
-aF-55.45234
-aF-222.3761
-aF99.75883
-aatp84
-a((F0.518498879904332
-I0
-tp85
-(lp86
-(lp87
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp88
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp89
-F250.0
-aF39.75952
-aF159.4444
-aF-71.52742
-aa(lp90
-F250.0
-aF-39.75952
-aF-159.4444
-aF71.52742
-aatp91
-a((F0.008877132600434398
-I0
-tp92
-(lp93
-(lp94
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp95
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp96
-F250.0
-aF55.46074
-aF222.4098
-aF-99.77394
-aa(lp97
-F250.0
-aF-55.46074
-aF-222.4098
-aF99.77394
-aatp98
-a((F0.006766402984646188
-I0
-tp99
-(lp100
-(lp101
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp102
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp103
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp104
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp105
-a((F0.15302159061039164
-I0
-tp106
-(lp107
-(lp108
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp109
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp110
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp111
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp112
-a((F0.11390477519377805
-I0
-tp113
-(lp114
-(lp115
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp116
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp117
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp118
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp119
-a((F0.08155081107088635
-I0
-tp120
-(lp121
-(lp122
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp123
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp124
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp125
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp126
-a((F0.0036753222317665395
-I0
-tp127
-(lp128
-(lp129
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp130
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp131
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp132
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp133
-a((F0.07593651679585227
-I0
-tp134
-(lp135
-(lp136
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp137
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp138
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp139
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp140
-a((F1.790807005072219
-I0
-tp141
-(lp142
-(lp143
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp144
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp145
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp146
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp147
-a((F0.395475132436254
-I0
-tp148
-(lp149
-(lp150
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp151
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp152
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp153
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp154
-a((F0.003624860653933186
-I0
-tp155
-(lp156
-(lp157
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp158
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp159
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp160
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp161
-a((F0.1631016221417725
-I0
-tp162
-(lp163
-(lp164
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp165
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp166
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp167
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp168
-a((F1.186425397308747
-I0
-tp169
-(lp170
-(lp171
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp172
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp173
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp174
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp175
-a((F0.24736033503494492
-I0
-tp176
-(lp177
-(lp178
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp179
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp180
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp181
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp182
-a((F0.0072598782536768035
-I0
-tp183
-(lp184
-(lp185
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp186
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp187
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp188
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp189
-a((F0.2250348696812327
-I0
-tp190
-(lp191
-(lp192
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp193
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp194
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp195
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp196
-a((F0.026684962843085858
-I0
-tp197
-(lp198
-(lp199
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp200
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp201
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp202
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp203
-a((F0.023354820365466676
-I0
-tp204
-(lp205
-(lp206
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp207
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp208
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp209
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp210
-a((F0.22097276866415483
-I0
-tp211
-(lp212
-(lp213
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp214
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp215
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp216
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp217
-a((F0.061776299282655774
-I0
-tp218
-(lp219
-(lp220
-F243.5549
-aF0.0
-aF0.0
-aF243.5095
-aa(lp221
-F256.4451
-aF0.0
-aF0.0
-aF-243.5095
-aa(lp222
-F243.5549
-aF54.02224
-aF216.6411
-aF-97.18608
-aa(lp223
-F256.4451
-aF-54.02224
-aF-216.6411
-aF97.18608
-aatp224
-a((F0.5851054196065437
-I0
-tp225
-(lp226
-(lp227
-F243.5549
-aF0.0
-aF0.0
-aF243.5095
-aa(lp228
-F256.4451
-aF0.0
-aF0.0
-aF-243.5095
-aa(lp229
-F243.5549
-aF54.02224
-aF216.6411
-aF-97.18608
-aa(lp230
-F256.4451
-aF-54.02224
-aF-216.6411
-aF97.18608
-aatp231
-a((F1.7262624364246424
-I0
-tp232
-(lp233
-(lp234
-F273.9133
-aF0.0
-aF0.0
-aF211.3007
-aa(lp235
-F226.0867
-aF0.0
-aF0.0
-aF-211.3007
-aa(lp236
-F273.9133
-aF46.87676
-aF187.9861
-aF-84.33135
-aa(lp237
-F226.0867
-aF-46.87676
-aF-187.9861
-aF84.33135
-aatp238
-a((F0.20725559701898785
-I0
-tp239
-(lp240
-(lp241
-F273.9133
-aF0.0
-aF0.0
-aF211.3007
-aa(lp242
-F226.0867
-aF0.0
-aF0.0
-aF-211.3007
-aa(lp243
-F273.9133
-aF46.87676
-aF187.9861
-aF-84.33135
-aa(lp244
-F226.0867
-aF-46.87676
-aF-187.9861
-aF84.33135
-aatp245
-a((F0.23009623586761052
-I0
-tp246
-(lp247
-(lp248
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp249
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp250
-F243.5359
-aF54.02666
-aF216.6588
-aF-97.19403
-aa(lp251
-F256.4641
-aF-54.02666
-aF-216.6588
-aF97.19403
-aatp252
-a((F0.030961673336769788
-I0
-tp253
-(lp254
-(lp255
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp256
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp257
-F243.5359
-aF54.02666
-aF216.6588
-aF-97.19403
-aa(lp258
-F256.4641
-aF-54.02666
-aF-216.6588
-aF97.19403
-aatp259
-a((F0.2827267193870959
-I0
-tp260
-(lp261
-(lp262
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp263
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp264
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp265
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp266
-a((F0.11173628232582335
-I0
-tp267
-(lp268
-(lp269
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp270
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp271
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp272
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp273
-a((F0.10947501935636043
-I0
-tp274
-(lp275
-(lp276
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp277
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp278
-F264.4
-aF52.26752
-aF209.6043
-aF-94.02934
-aa(lp279
-F235.6
-aF-52.26752
-aF-209.6043
-aF94.02934
-aatp280
-a((F0.21784650821434426
-I0
-tp281
-(lp282
-(lp283
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp284
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp285
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp286
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp287
-a((F0.020961233306794114
-I0
-tp288
-(lp289
-(lp290
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp291
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp292
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp293
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp294
-a((F0.13129073989589382
-I0
-tp295
-(lp296
-(lp297
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp298
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp299
-F249.9968
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp300
-F250.0032
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp301
-a((F0.19957331011891252
-I0
-tp302
-(lp303
-(lp304
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp305
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp306
-F241.6879
-aF53.61667
-aF215.0146
-aF-96.45645
-aa(lp307
-F258.3121
-aF-53.61667
-aF-215.0146
-aF96.45645
-aatp308
-a((F0.10943521099365543
-I0
-tp309
-(lp310
-(lp311
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp312
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp313
-F264.3968
-aF52.26674
-aF209.6011
-aF-94.02793
-aa(lp314
-F235.6032
-aF-52.26674
-aF-209.6011
-aF94.02793
-aatp315
-a((F2.607967194937269
-I0
-tp316
-(lp317
-(lp318
-F250.0
-aF0.0
-aF0.0
-aF236.7125
-aa(lp319
-F250.0
-aF0.0
-aF0.0
-aF-236.7125
-aa(lp320
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp321
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp322
-a((F0.15302159061039167
-I0
-tp323
-(lp324
-(lp325
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp326
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp327
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp328
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp329
-a((F0.11390477519377806
-I0
-tp330
-(lp331
-(lp332
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp333
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp334
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp335
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp336
-a((F0.08155081107088634
-I0
-tp337
-(lp338
-(lp339
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp340
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp341
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp342
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp343
-a((F0.0036753222317665395
-I0
-tp344
-(lp345
-(lp346
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp347
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp348
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp349
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp350
-a((F0.07593651679585228
-I0
-tp351
-(lp352
-(lp353
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp354
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp355
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp356
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp357
-a((F1.790807005072219
-I0
-tp358
-(lp359
-(lp360
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp361
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp362
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp363
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp364
-a((F0.3954751324362539
-I0
-tp365
-(lp366
-(lp367
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp368
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp369
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp370
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp371
-a((F0.003624860653933185
-I0
-tp372
-(lp373
-(lp374
-F248.152
-aF0.0
-aF0.0
-aF234.7598
-aa(lp375
-F251.848
-aF0.0
-aF0.0
-aF-234.7598
-aa(lp376
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp377
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp378
-a((F0.16310162214177257
-I0
-tp379
-(lp380
-(lp381
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp382
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp383
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp384
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp385
-a((F1.1864253973087473
-I0
-tp386
-(lp387
-(lp388
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp389
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp390
-F248.152
-aF52.08113
-aF208.8568
-aF-93.69403
-aa(lp391
-F251.848
-aF-52.08113
-aF-208.8568
-aF93.69403
-aatp392
-a((F0.24736033503494492
-I0
-tp393
-(lp394
-(lp395
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp396
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp397
-F257.9328
-aF50.65212
-aF203.1262
-aF-91.12324
-aa(lp398
-F242.0672
-aF-50.65212
-aF-203.1262
-aF91.12324
-aatp399
-a((F0.0072598782536768035
-I0
-tp400
-(lp401
-(lp402
-F257.9328
-aF0.0
-aF0.0
-aF228.3185
-aa(lp403
-F242.0672
-aF0.0
-aF0.0
-aF-228.3185
-aa(lp404
-F250.0032
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp405
-F249.9968
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp406
-a((F0.02668496284308585
-I0
-tp407
-(lp408
-(lp409
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp410
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp411
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp412
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp413
-a((F0.22503486968123265
-I0
-tp414
-(lp415
-(lp416
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp417
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp418
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp419
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp420
-a((F0.22097276866415477
-I0
-tp421
-(lp422
-(lp423
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp424
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp425
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp426
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp427
-a((F0.023354820365466687
-I0
-tp428
-(lp429
-(lp430
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp431
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp432
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp433
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp434
-a((F0.5851054196065434
-I0
-tp435
-(lp436
-(lp437
-F243.5549
-aF0.0
-aF0.0
-aF243.5095
-aa(lp438
-F256.4451
-aF0.0
-aF0.0
-aF-243.5095
-aa(lp439
-F243.5549
-aF54.02224
-aF216.6411
-aF-97.18608
-aa(lp440
-F256.4451
-aF-54.02224
-aF-216.6411
-aF97.18608
-aatp441
-a((F0.0617762992826556
-I0
-tp442
-(lp443
-(lp444
-F243.5549
-aF0.0
-aF0.0
-aF243.5095
-aa(lp445
-F256.4451
-aF0.0
-aF0.0
-aF-243.5095
-aa(lp446
-F243.5549
-aF54.02224
-aF216.6411
-aF-97.18608
-aa(lp447
-F256.4451
-aF-54.02224
-aF-216.6411
-aF97.18608
-aatp448
-a((F0.20725559701898788
-I0
-tp449
-(lp450
-(lp451
-F273.9133
-aF0.0
-aF0.0
-aF211.3007
-aa(lp452
-F226.0867
-aF0.0
-aF0.0
-aF-211.3007
-aa(lp453
-F273.9133
-aF46.87676
-aF187.9861
-aF-84.33135
-aa(lp454
-F226.0867
-aF-46.87676
-aF-187.9861
-aF84.33135
-aatp455
-a((F1.7262624364246415
-I0
-tp456
-(lp457
-(lp458
-F273.9133
-aF0.0
-aF0.0
-aF211.3007
-aa(lp459
-F226.0867
-aF0.0
-aF0.0
-aF-211.3007
-aa(lp460
-F273.9133
-aF46.87676
-aF187.9861
-aF-84.33135
-aa(lp461
-F226.0867
-aF-46.87676
-aF-187.9861
-aF84.33135
-aatp462
-a((F0.030961673336769805
-I0
-tp463
-(lp464
-(lp465
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp466
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp467
-F243.5359
-aF54.02666
-aF216.6588
-aF-97.19403
-aa(lp468
-F256.4641
-aF-54.02666
-aF-216.6588
-aF97.19403
-aatp469
-a((F0.282726719387096
-I0
-tp470
-(lp471
-(lp472
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp473
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp474
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp475
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp476
-a((F0.1117362823258234
-I0
-tp477
-(lp478
-(lp479
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp480
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp481
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp482
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp483
-a((F0.10947501935636035
-I0
-tp484
-(lp485
-(lp486
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp487
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp488
-F264.4
-aF52.26752
-aF209.6043
-aF-94.02934
-aa(lp489
-F235.6
-aF-52.26752
-aF-209.6043
-aF94.02934
-aatp490
-a((F0.23009623586761052
-I0
-tp491
-(lp492
-(lp493
-F243.5359
-aF0.0
-aF0.0
-aF243.5295
-aa(lp494
-F256.4641
-aF0.0
-aF0.0
-aF-243.5295
-aa(lp495
-F243.5359
-aF54.02666
-aF216.6588
-aF-97.19403
-aa(lp496
-F256.4641
-aF-54.02666
-aF-216.6588
-aF97.19403
-aatp497
-a((F0.02096123330679411
-I0
-tp498
-(lp499
-(lp500
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp501
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp502
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp503
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp504
-a((F0.13129073989589382
-I0
-tp505
-(lp506
-(lp507
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp508
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp509
-F249.9968
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp510
-F250.0032
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp511
-a((F0.1995733101189124
-I0
-tp512
-(lp513
-(lp514
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp515
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp516
-F241.6879
-aF53.61667
-aF215.0146
-aF-96.45645
-aa(lp517
-F258.3121
-aF-53.61667
-aF-215.0146
-aF96.45645
-aatp518
-a((F0.10943521099365548
-I0
-tp519
-(lp520
-(lp521
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp522
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp523
-F264.3968
-aF52.26674
-aF209.6011
-aF-94.02793
-aa(lp524
-F235.6032
-aF-52.26674
-aF-209.6011
-aF94.02793
-aatp525
-a((F0.21784650821434415
-I0
-tp526
-(lp527
-(lp528
-F243.5328
-aF0.0
-aF0.0
-aF243.5328
-aa(lp529
-F256.4672
-aF0.0
-aF0.0
-aF-243.5328
-aa(lp530
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp531
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp532
-a((F0.4901153733199385
-I0
-tp533
-(lp534
-(lp535
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp536
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp537
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp538
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp539
-a((F0.029396815093657935
-I0
-tp540
-(lp541
-(lp542
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp543
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp544
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp545
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp546
-a((F0.001837300943353621
-I0
-tp547
-(lp548
-(lp549
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp550
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp551
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp552
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp553
-a((F0.001838208457272684
-I0
-tp554
-(lp555
-(lp556
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp557
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp558
-F250.0
-aF55.45234
-aF222.3761
-aF-99.75883
-aa(lp559
-F250.0
-aF-55.45234
-aF-222.3761
-aF99.75883
-aatp560
-a((F0.03574211792602636
-I0
-tp561
-(lp562
-(lp563
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp564
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp565
-F250.0
-aF39.75952
-aF159.4444
-aF-71.52742
-aa(lp566
-F250.0
-aF-39.75952
-aF-159.4444
-aF71.52742
-aatp567
-a((F0.0496106298749928
-I0
-tp568
-(lp569
-(lp570
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp571
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp572
-F250.0
-aF55.46074
-aF222.4098
-aF-99.77394
-aa(lp573
-F250.0
-aF-55.46074
-aF-222.4098
-aF99.77394
-aatp574
-a((F1.237769860999663
-I0
-tp575
-(lp576
-(lp577
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp578
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp579
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp580
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp581
-a((F0.020872161925933056
-I0
-tp582
-(lp583
-(lp584
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp585
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp586
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp587
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp588
-a((F0.006691468526472817
-I0
-tp589
-(lp590
-(lp591
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp592
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp593
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp594
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp595
-a((F0.00672304572338048
-I0
-tp596
-(lp597
-(lp598
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp599
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp600
-F250.0
-aF55.45234
-aF222.3761
-aF-99.75883
-aa(lp601
-F250.0
-aF-55.45234
-aF-222.3761
-aF99.75883
-aatp602
-a((F0.1896056386238811
-I0
-tp603
-(lp604
-(lp605
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp606
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp607
-F250.0
-aF39.75952
-aF159.4444
-aF-71.52742
-aa(lp608
-F250.0
-aF-39.75952
-aF-159.4444
-aF71.52742
-aatp609
-a((F0.013604934640539718
-I0
-tp610
-(lp611
-(lp612
-F241.6847
-aF0.0
-aF0.0
-aF241.6847
-aa(lp613
-F258.3153
-aF0.0
-aF0.0
-aF-241.6847
-aa(lp614
-F250.0
-aF55.46074
-aF222.4098
-aF-99.77394
-aa(lp615
-F250.0
-aF-55.46074
-aF-222.4098
-aF99.77394
-aatp616
-a((F0.5231244438931428
-I0
-tp617
-(lp618
-(lp619
-F235.6
-aF0.0
-aF0.0
-aF235.6
-aa(lp620
-F264.4
-aF0.0
-aF0.0
-aF-235.6
-aa(lp621
-F250.0
-aF52.51432
-aF210.594
-aF-94.47333
-aa(lp622
-F250.0
-aF-52.51432
-aF-210.594
-aF94.47333
-aatp623
-a((F9.880082938079737e-05
-I0
-tp624
-(lp625
-(lp626
-F235.6
-aF0.0
-aF0.0
-aF235.6
-aa(lp627
-F264.4
-aF0.0
-aF0.0
-aF-235.6
-aa(lp628
-F250.0
-aF55.45234
-aF222.3761
-aF-99.75883
-aa(lp629
-F250.0
-aF-55.45234
-aF-222.3761
-aF99.75883
-aatp630
-a((F0.2684799994270568
-I0
-tp631
-(lp632
-(lp633
-F235.6
-aF0.0
-aF0.0
-aF235.6
-aa(lp634
-F264.4
-aF0.0
-aF0.0
-aF-235.6
-aa(lp635
-F250.0
-aF39.75952
-aF159.4444
-aF-71.52742
-aa(lp636
-F250.0
-aF-39.75952
-aF-159.4444
-aF71.52742
-aatp637
-a((F5.3096175922586565e-05
-I0
-tp638
-(lp639
-(lp640
-F235.6
-aF0.0
-aF0.0
-aF235.6
-aa(lp641
-F264.4
-aF0.0
-aF0.0
-aF-235.6
-aa(lp642
-F250.0
-aF55.46074
-aF222.4098
-aF-99.77394
-aa(lp643
-F250.0
-aF-55.46074
-aF-222.4098
-aF99.77394
-aatp644
-a((F0.17236628114329705
-I0
-tp645
-(lp646
-(lp647
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp648
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp649
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp650
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp651
-a((F0.04309157028582426
-I0
-tp652
-(lp653
-(lp654
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp655
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp656
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp657
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp658
-a((F0.043112854877208236
-I0
-tp659
-(lp660
-(lp661
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp662
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp663
-F250.0
-aF55.45234
-aF222.3761
-aF-99.75883
-aa(lp664
-F250.0
-aF-55.45234
-aF-222.3761
-aF99.75883
-aatp665
-a((F0.20957154465428646
-I0
-tp666
-(lp667
-(lp668
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp669
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp670
-F250.0
-aF39.75952
-aF159.4444
-aF-71.52742
-aa(lp671
-F250.0
-aF-39.75952
-aF-159.4444
-aF71.52742
-aatp672
-a((F0.01289355401139826
-I0
-tp673
-(lp674
-(lp675
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp676
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp677
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp678
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp679
-a((F0.007166475491463663
-I0
-tp680
-(lp681
-(lp682
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp683
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp684
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp685
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp686
-a((F0.6048040098451131
-I0
-tp687
-(lp688
-(lp689
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp690
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp691
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp692
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp693
-a((F0.01289355401139826
-I0
-tp694
-(lp695
-(lp696
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp697
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp698
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp699
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp700
-a((F0.0071664754914636634
-I0
-tp701
-(lp702
-(lp703
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp704
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp705
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp706
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp707
-a((F0.6048040098451131
-I0
-tp708
-(lp709
-(lp710
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp711
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp712
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp713
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp714
-a((F0.0008058471257123913
-I0
-tp715
-(lp716
-(lp717
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp718
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp719
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp720
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp721
-a((F0.002297521711791924
-I0
-tp722
-(lp723
-(lp724
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp725
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp726
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp727
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp728
-a((F0.0887099794045066
-I0
-tp729
-(lp730
-(lp731
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp732
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp733
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp734
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp735
-a((F0.0008058471257123913
-I0
-tp736
-(lp737
-(lp738
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp739
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp740
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp741
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp742
-a((F0.002297521711791923
-I0
-tp743
-(lp744
-(lp745
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp746
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp747
-F241.6847
-aF53.61742
-aF215.0176
-aF-96.4578
-aa(lp748
-F258.3153
-aF-53.61742
-aF-215.0176
-aF96.4578
-aatp749
-a((F0.0887099794045066
-I0
-tp750
-(lp751
-(lp752
-F250.0
-aF0.0
-aF0.0
-aF250.0
-aa(lp753
-F250.0
-aF0.0
-aF0.0
-aF-250.0
-aa(lp754
-F250.0
-aF55.46214
-aF222.4154
-aF-99.77646
-aa(lp755
-F250.0
-aF-55.46214
-aF-222.4154
-aF99.77646
-aatp756
-a((F0.0008055261451803681
-I0
-tp757
-(lp758
-(lp759
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp760
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp761
-F249.9779
-aF55.45724
-aF222.3957
-aF-99.76765
-aa(lp762
-F250.0221
-aF-55.45724
-aF-222.3957
-aF99.76765
-aatp763
-a((F0.0022982520201042457
-I0
-tp764
-(lp765
-(lp766
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp767
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp768
-F241.7068
-aF53.61218
-aF214.9966
-aF-96.44838
-aa(lp769
-F258.2932
-aF-53.61218
-aF-214.9966
-aF96.44838
-aatp770
-a((F3.891322460982593e-07
-I0
-tp771
-(lp772
-(lp773
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp774
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp775
-F235.6221
-aF52.26202
-aF209.5822
-aF-94.01945
-aa(lp776
-F264.3779
-aF-52.26202
-aF-209.5822
-aF94.01945
-aatp777
-a((F0.08870250090799825
-I0
-tp778
-(lp779
-(lp780
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp781
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp782
-F250.0221
-aF55.45724
-aF222.3957
-aF-99.76765
-aa(lp783
-F249.9779
-aF-55.45724
-aF-222.3957
-aF99.76765
-aatp784
-a((F0.0008055261451803681
-I0
-tp785
-(lp786
-(lp787
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp788
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp789
-F249.9779
-aF55.45724
-aF222.3957
-aF-99.76765
-aa(lp790
-F250.0221
-aF-55.45724
-aF-222.3957
-aF99.76765
-aatp791
-a((F0.0022982520201042457
-I0
-tp792
-(lp793
-(lp794
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp795
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp796
-F241.7068
-aF53.61218
-aF214.9966
-aF-96.44838
-aa(lp797
-F258.2932
-aF-53.61218
-aF-214.9966
-aF96.44838
-aatp798
-a((F3.891322460982593e-07
-I0
-tp799
-(lp800
-(lp801
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp802
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp803
-F235.6221
-aF52.26202
-aF209.5822
-aF-94.01945
-aa(lp804
-F264.3779
-aF-52.26202
-aF-209.5822
-aF94.01945
-aatp805
-a((F0.08870250090799825
-I0
-tp806
-(lp807
-(lp808
-F249.9779
-aF0.0
-aF0.0
-aF249.9779
-aa(lp809
-F250.0221
-aF0.0
-aF0.0
-aF-249.9779
-aa(lp810
-F250.0221
-aF55.45724
-aF222.3957
-aF-99.76765
-aa(lp811
-F249.9779
-aF-55.45724
-aF-222.3957
-aF99.76765
-aatp812
-a((F0.008144164219845314
-I0
-tp813
-(lp814
-(lp815
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp816
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp817
-F219.6195
-aF48.72227
-aF195.387
-aF-87.65143
-aa(lp818
-F280.3805
-aF-48.72227
-aF-195.387
-aF87.65143
-aatp819
-a((F0.010506069043620324
-I0
-tp820
-(lp821
-(lp822
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp823
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp824
-F272.0652
-aF46.34406
-aF185.8499
-aF-83.37302
-aa(lp825
-F227.9348
-aF-46.34406
-aF-185.8499
-aF83.37302
-aatp826
-a((F0.04053214453020507
-I0
-tp827
-(lp828
-(lp829
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp830
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp831
-F234.0195
-aF44.57176
-aF178.7426
-aF-80.18466
-aa(lp832
-F265.9805
-aF-44.57176
-aF-178.7426
-aF80.18466
-aatp833
-a((F0.38202214631009856
-I0
-tp834
-(lp835
-(lp836
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp837
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp838
-F219.6195
-aF48.72227
-aF195.387
-aF-87.65143
-aa(lp839
-F280.3805
-aF-48.72227
-aF-195.387
-aF87.65143
-aatp840
-a((F0.008144164219845314
-I0
-tp841
-(lp842
-(lp843
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp844
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp845
-F219.6195
-aF48.72227
-aF195.387
-aF-87.65143
-aa(lp846
-F280.3805
-aF-48.72227
-aF-195.387
-aF87.65143
-aatp847
-a((F0.010506069043620328
-I0
-tp848
-(lp849
-(lp850
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp851
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp852
-F272.0652
-aF46.34406
-aF185.8499
-aF-83.37302
-aa(lp853
-F227.9348
-aF-46.34406
-aF-185.8499
-aF83.37302
-aatp854
-a((F0.04053214453020507
-I0
-tp855
-(lp856
-(lp857
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp858
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp859
-F234.0195
-aF44.57176
-aF178.7426
-aF-80.18466
-aa(lp860
-F265.9805
-aF-44.57176
-aF-178.7426
-aF80.18466
-aatp861
-a((F0.38202214631009856
-I0
-tp862
-(lp863
-(lp864
-F219.6195
-aF0.0
-aF0.0
-aF219.6195
-aa(lp865
-F280.3805
-aF0.0
-aF0.0
-aF-219.6195
-aa(lp866
-F219.6195
-aF48.72227
-aF195.387
-aF-87.65143
-aa(lp867
-F280.3805
-aF-48.72227
-aF-195.387
-aF87.65143
-aatp868
-a((F0.19693610984384063
-I0
-tp869
-(lp870
-(lp871
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp872
-F250.0032
-aF0.0
-aF0.0
-aF-249.9968
-aa(lp873
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp874
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp875
-a((F0.06526989955179051
-I0
-tp876
-(lp877
-(lp878
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp879
-F250.0032
-aF0.0
-aF0.0
-aF-249.9968
-aa(lp880
-F249.9968
-aF55.46144
-aF222.4126
-aF-99.7752
-aa(lp881
-F250.0032
-aF-55.46144
-aF-222.4126
-aF99.7752
-aatp882
-a((F0.014001276293925691
-I0
-tp883
-(lp884
-(lp885
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp886
-F250.0032
-aF0.0
-aF0.0
-aF-249.9968
-aa(lp887
-F241.6879
-aF53.61667
-aF215.0146
-aF-96.45645
-aa(lp888
-F258.3121
-aF-53.61667
-aF-215.0146
-aF96.45645
-aatp889
-a((F7.641132099338047e-06
-I0
-tp890
-(lp891
-(lp892
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp893
-F250.0032
-aF0.0
-aF0.0
-aF-249.9968
-aa(lp894
-F264.3968
-aF52.26674
-aF209.6011
-aF-94.02793
-aa(lp895
-F235.6032
-aF-52.26674
-aF-209.6011
-aF94.02793
-aatp896
-a((F0.19693610984384072
-I0
-tp897
-(lp898
-(lp899
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp900
-F250.0032
-aF0.0
-aF0.0
-aF-249.9968
-aa(lp901
-F243.5328
-aF54.0274
-aF216.6618
-aF-97.19536
-aa(lp902
-F256.4672
-aF-54.0274
-aF-216.6618
-aF97.19536
-aatp903
-a((F0.06526989955179051
-I0
-tp904
-(lp905
-(lp906
-F249.9968
-aF0.0
-aF0.0
-aF249.9968
-aa(lp907
-F250.0032