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

[Olivier Mattelaer]

Add the install command

[Johan Alwall (johan-alwall)]

Hello Olivier,

This is really awesome! Very useful and helpful! I have found one thing right away: The packageinfo file is gotten from the UCL cluster (as far as I can tell below) but the package is being downloaded from the UIUC cluster. If there is a mismatch between the clusters, I expect that this doesn't work. Isn't it better to do everything from the UIUC cluster? Or is the file on the UCL cluster updated whenever the UIUC cluster is updated?

Cheers,
Johan

[Johan Alwall (johan-alwall)]

Hello again,

- I would recommend to remove the debug printout of the contents of the package_info file :-)
It would be useful however to have an INFO statement saying from where the presently downloaded file is downloaded.

- For Delphes, it's crucial that you get the version WITHOUT Frog, since this is not standard on any computer, and not used in the standard running: http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes_V_NoFROG_1.9.tar.gz

- Problems when running:
mg5>install MadAnalysis
gzip: stdin: not in gzip format
tar: Child returned status 1
tar: Error is not recoverable: exiting now
Command "install MadAnalysis" interrupted with error:
IndexError : list index out of range
MG5_debug:
Traceback (most recent call last):
  File "/Users/alwall/MadGraph/new_install_command/madgraph/interface/extended_cmd.py", line 151, in onecmd
    return cmd.Cmd.onecmd(self, line)
  File "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/cmd.py", line 219, in onecmd
    return func(arg)
  File "/Users/alwall/MadGraph/new_install_command/madgraph/interface/cmd_interface.py", line 2512, in do_install
    and not n.endswith('gz')][0]
IndexError: list index out of range

That's what I find so far :-)
Johan

review: Needs Fixing

[Olivier Mattelaer (olivier-mattelaer)]

Hi Johan, thanks for those very useful comment:

>If there is a mismatch between the clusters, I expect that this doesn't work.
Indeed that was the problem with MadAnalysis, where the last version of MadAnalysis was not present on uiuc. Stupid mistake when doing the copy past. Thanks a lot for this.

>- I would recommend to remove the debug printout of the contents of the package_info file :-)
>It would be useful however to have an INFO statement saying from where the presently downloaded >file is downloaded.

Sure Done

>- For Delphes, it's crucial that you get the version WITHOUT Frog, since this is not standard >on >any computer, and not used in the standard running: http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes>/files/Delphes_V_NoFROG_1.9.tar.gz

- Problems when running:
mg5>install MadAnalysis
See first point -> Done.
I also add more beautifull error message in case of no connection to the server/problem with the download.

Ok I'll modify this. But I'm confuse since the Downloads page of MG indicates the one with frog then... It also seems that it refer to version 1.8 and not to version 1.9.

[lp:~maddevelopers/mg5amcnlo/new_install_command]

173. By Olivier Mattelaer on 2011-10-04

include review modification

174. By Olivier Mattelaer on 2011-10-04

Adding TopEffTh (ask By Scott approved by Tim)

175. By Olivier Mattelaer on 2011-10-04

add propagating=False in TopEffTh model. Remove un-propagating particle in display/auto-completion

176. By Olivier Mattelaer on 2011-10-05

fix aloha bug 868645.

177. By Olivier Mattelaer on 2011-10-05

Update Update Note.

178. By Olivier Mattelaer on 2011-10-05

Forbid to install ExRootAnalysis if root not installed or ROOTSYS not configured

[Johan Alwall (johan-alwall)]

Hello Olivier,

Note that in the diff below, it seems that there are unresolved conflicts in the UpdateNotes file.

Issue: I downloaded MadAnalysis, and td was automatically downloaded. However, when I ran a production, no plots were generated, and the plots.log file says "td not found". It seems that you miss a ln or mv in the installation to copy the correct td_xxx file to td. If you can't do this automatically, there should at least be a clear warning message to the user.
Also, if MadAnalysis exists and plot fails, I think we should pipe the error message also to the run log, so the user can see what the problem is.

There's also an issue with Delphes now:
mg5>install Delphes
Downloading http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes_V_No_FROG_1.9.tar.gz
  % Total % Received % Xferd Average Speed Time Time Time Current
                                 Dload Upload Total Spent Left Speed
112 337 112 337 0 0 119 0 0:00:02 0:00:02 --:--:-- 492

gzip: stdin: not in gzip format
tar: Child returned status 1
tar: Error is not recoverable: exiting now
Command "install Delphes" interrupted with error:
MadGraph5Error : The file was not loaded correctly. Stop
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in 'MG5_debug'.
Please attach this file to your report.

The correct address is http://www.fynu.ucl.ac.be/users/s.ovyn/Delphes/files/Delphes_V_NoFROG_1.9.tar.gz

Besides this (and the issues with compiling Pythia/PGS on Mac, which should also be solved before the school), everything now seems to work very well. Thanks!

Johan

review: Needs Fixing

[Johan Alwall (johan-alwall)]

Hi again Olivier,

You might want a more informative error message if a download fails - this might happen also in the future, and

gzip: stdin: not in gzip format
tar: Child returned status 1
tar: Error is not recoverable: exiting now
Command "install Delphes" interrupted with error:
MadGraph5Error : The file was not loaded correctly. Stop

is not very informative :-)

[Johan Alwall (johan-alwall)]

Very nice job checking for root and ROOTSYS before installing ExRootAnalysis by the way :-)

[Olivier Mattelaer (olivier-mattelaer)]

Hi johan,

I didn't know about the -m32 -m64.
One thing is not clear in your message, is it only for Mac or also for linux as well?
(It should be both right?)

Sorry for the wrong path for Delphes,
And I will create a nicer error message.
For td I'm going to check what the problem is.

Thanks,

Olivier

[Johan Alwall (johan-alwall)]

Hello Olivier,

You see it in the stdhep_Arch file. It's only for Mac.

For td, the problem is that in the tar, there are two td exe's, one for intel and one for g4, called td_intel and td_g4. The correct one needs to be copied to "td".

Thanks,
Johan

[Olivier Mattelaer (olivier-mattelaer)]

Ok Thanks,

I assume that I can 'safely' always assume that all computer are intel
now.

On Oct 6, 2011, at 10:22 AM, Johan Alwall wrote:

> Hello Olivier,
>
> You see it in the stdhep_Arch file. It's only for Mac.
>
> For td, the problem is that in the tar, there are two td exe's, one
> for intel and one for g4, called td_intel and td_g4. The correct one
> needs to be copied to "td".
>
> Thanks,
> Johan
>
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/
> new_install_command/+merge/78015
> You proposed lp:~maddevelopers/madgraph5/new_install_command for
> merging.

[lp:~maddevelopers/mg5amcnlo/new_install_command]

179. By mattelaer-olivier on 2011-10-06

fix the td problem. change -m32 to -m64 in 64 bit computer
improve error message

180. By mattelaer-olivier on 2011-10-06

improve help for root

181. By mattelaer-olivier on 2011-10-06

merge with the trunk

182. By mattelaer-olivier on 2011-10-06

fix the problem between g77 and the HELAS routine (Thanks Tim)

183. By mattelaer-olivier on 2011-10-06

remove all include in HELAS routine (g77 compatibility)

[Olivier Mattelaer (olivier-mattelaer)]

Hi Johan,

by testing pythia-pgs on a mac, I've the following error:

WRONG LHAPDF set number = 0 given! STOP EXE!
 ==== PYTHIA WILL USE LHAPDF ====

I remember that you suggest to modify one line of pythia-pgs whne we was at fermilab and that resolve my problem. Do you remember what this is?

If you remember then I would be able to mofidy this line automaticaly on any Mac computer.

Cheers,

Olivier

[lp:~maddevelopers/mg5amcnlo/new_install_command]

184. By mattelaer-olivier on 2011-10-06

Fix problem name when download from linux

185. By mattelaer-olivier on 2011-10-07

merge with the last version

[Johan Alwall (johan-alwall)]

Hi Olivier,

I have fixed arch_mcfio.

However, I get the following error now when I download MadAnalysis:

Downloading TD for Mac
  % Total % Received % Xferd Average Speed Time Time Time Current
                                 Dload Upload Total Spent Left Speed
100 1949k 100 1949k 0 0 30652 0 0:01:05 0:01:05 --:--:-- 50303
libg2c.0.dylib
readme_intel.txt
td.1
td_g4
td_mac_intel
topdrawer.doc
Command "install MadAnalysis" interrupted with error:
IOError : [Errno 2] No such file or directory: 'td/td_mac_intel'
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in 'MG5_debug'.
Please attach this file to your report.

review: Needs Fixing

[lp:~maddevelopers/mg5amcnlo/new_install_command]

186. By Johan Alwall on 2011-10-08

Fix also arch_mcfio file for -m32/-m64

187. By Olivier Mattelaer on 2011-10-09

fix the path to be universal for the copy of td

[Johan Alwall (johan-alwall)]

Hello Olivier,

Thanks for fixing the copy of td. However, now td doesn't run, and the reason is that the td directory has to be in the DYLD_LIBRARY_PATH. So this needs to be set before td is run (on Mac).

Cheers,
Johan

review: Needs Fixing

[Olivier Mattelaer (olivier-mattelaer)]

> Thanks for fixing the copy of td. However, now td doesn't run, and
> the reason is that the td directory has to be in the
> DYLD_LIBRARY_PATH. So this needs to be set before td is run (on Mac).

This variable is not set on my computer...

[@olivier-mattelaers-macbook ~]$ echo $DYLD_LIBRARY_PATH
/Users/omatt/programme/root/lib:
[@olivier-mattelaers-macbook ~]$

Cheers,

Olivier

On 09-oct.-11, at 22:56, Johan Alwall wrote:

> Review: Needs Fixing
>
> Hello Olivier,
>
> Thanks for fixing the copy of td. However, now td doesn't run, and
> the reason is that the td directory has to be in the
> DYLD_LIBRARY_PATH. So this needs to be set before td is run (on Mac).
>
> Cheers,
> Johan
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/
> new_install_command/+merge/78015
> You proposed lp:~maddevelopers/madgraph5/new_install_command for
> merging.

[Johan Alwall (johan-alwall)]

Hello Olivier,

The missing library, which is included with td, is libg2c.0.dylib. It is possible that you have this library in your default path (e.g., /usr/local/lib/). For those who don't have it, it is found in the td directory, and the td directory must be added to the DYLD_LIBRARY_PATH before running td. This needs to be done in the plot script (or your equivalent script, where td is being run in the new scripts).

Cheers,
Johan

review: Needs Fixing

[lp:~maddevelopers/mg5amcnlo/new_install_command]

188. By Olivier Mattelaer on 2011-10-10

add the td path to DYLIB. Usefull for some computer

[Olivier Mattelaer (olivier-mattelaer)]

I've added it in the plot script.
Could you try if this works?

Thanks,

Olivier

[lp:~maddevelopers/mg5amcnlo/new_install_command]

189. By Johan Alwall on 2011-10-11

Fixed the LD_LIBRARY_PATH for td in plot (needed for Mac OS X)

[Johan Alwall (johan-alwall)]

Ok with this final point fixed, I approve the merge.

Johan

review: Approve

=== modified file 'HELAS/hstlxx.F'
--- HELAS/hstlxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/hstlxx.F	2011-10-11 00:36:19 +0000
@@ -15,7 +15,8 @@
 c       complex hst(3)         : scalar current              j(s':s1,s2)
 c     
       implicit none
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
 
       double complex s1(DIM),t2(DIM),hst(DIM),dg
       double precision q(0:3),smass,swidth,q2,gc

=== modified file 'HELAS/httaxx.F'
--- HELAS/httaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/httaxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex tc1(DIM),tc2(DIM),jsc(DIM)
       double complex vertex,dj,gc(2)
       double precision mass,width,q2,q(4)

=== modified file 'HELAS/httsxx.F'
--- HELAS/httsxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/httsxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include 'dimension.inc'
       double complex tc1(DIM),htts(DIM),sc(DIM),tc2(DIM)
       double complex dg,g1(2)
 

=== modified file 'HELAS/hvvhxx.F'
--- HELAS/hvvhxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/hvvhxx.F	2011-10-11 00:36:19 +0000
@@ -20,9 +20,9 @@
 c
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex v1(DIM),v2(DIM),jsvv(DIM),vertex1,vertex2,dj
       double complex v12,p2v1,p1v2	
       double complex v13,v14,v23,v24,v34

=== modified file 'HELAS/hvvshx.F'
--- HELAS/hvvshx.F	2010-10-30 03:26:37 +0000
+++ HELAS/hvvshx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
 c     
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex v1(DIM),v2(DIM),sc(DIM),hvvsh(DIM)
       double complex hvvsh1,hvvsh2,dg,g1(2)
       double complex v12,p2v1,p1v2,v13,v14,v23,v24,v34

=== modified file 'HELAS/hvvvxx.F'
--- HELAS/hvvvxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/hvvvxx.F	2011-10-11 00:36:19 +0000
@@ -22,9 +22,9 @@
       implicit none
 
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),gc(DIM),jhvvv(DIM)
 
       double complex dvertx, vertex, vertex1, vertex2,dj

=== modified file 'HELAS/iosxxx.F'
--- HELAS/iosxxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/iosxxx.F	2011-10-11 00:36:19 +0000
@@ -14,7 +14,10 @@
 c     
       implicit none
       double complex fi(6),fo(6),sc(3),gc(2),vertex
+      INTEGER DIM
+      PARAMETER(DIM=18)
 
+c      include 'dimension.inc'
 #ifdef HELAS_CHECK
       double precision p0,p1,p2,p3,q0,q1,q2,q3,r0,r1,r2,r3,pm
       double precision epsi

=== modified file 'HELAS/jioxxx.F'
--- HELAS/jioxxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jioxxx.F	2011-10-11 00:36:19 +0000
@@ -24,7 +24,6 @@
       parameter( rZero = 0.0d0, rOne = 1.0d0 )
       double complex cImag, cZero
       parameter( cImag = ( 0.0d0, 1.0d0 ), cZero = ( 0.0d0, 0.0d0 ) )
-
 #ifdef HELAS_CHECK
       integer stdo
       parameter( stdo = 6 )

=== modified file 'HELAS/jvshxx.F'
--- HELAS/jvshxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jvshxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
 c
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex vc(DIM),sc(DIM),jvs(DIM),jvs1(DIM),jvs2(DIM)
       double complex qvc,dg
       double precision qp1,p12,p13,p14,p23,p24,p34

=== modified file 'HELAS/jvsshx.F'
--- HELAS/jvsshx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jvsshx.F	2011-10-11 00:36:19 +0000
@@ -8,9 +8,9 @@
 c
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex vc(DIM),sc1(DIM),sc2(DIM),jvs(DIM)
       double complex jvs1(DIM),jvs2(DIM),dg,qvc
       double precision qp1,p12,p13,p14,p23,p24,p34

=== modified file 'HELAS/jvtaxx.F'
--- HELAS/jvtaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jvtaxx.F	2011-10-11 00:36:19 +0000
@@ -20,9 +20,9 @@
       implicit none
 
 c dimension of the current set to arbitrary length
-c      integer DIM
-c      parameter (DIM=18)
-      include "dimension.inc"
+      integer DIM
+      parameter (DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),jw(DIM),tc(DIM)
       double complex gt1(4), gt2(4)
       double precision q(0:3),g,q2,dv

=== modified file 'HELAS/jvvsxx.F'
--- HELAS/jvvsxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jvvsxx.F	2011-10-11 00:36:19 +0000
@@ -22,9 +22,9 @@
  
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),jggs(DIM),sc(DIM)
       double complex jggs1(DIM),jggs2(DIM)
       double complex sva,svb,vab,j12(0:3)

=== modified file 'HELAS/jvvtlx.F'
--- HELAS/jvvtlx.F	2010-10-30 03:26:37 +0000
+++ HELAS/jvvtlx.F	2011-10-11 00:36:19 +0000
@@ -21,9 +21,9 @@
  
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),jggs(DIM),sc(DIM),g2(2)
       double complex jggs1(DIM),jggs2(DIM)
       double complex sva,svb,vab,j12(0:3)

=== modified file 'HELAS/sstlxx.F'
--- HELAS/sstlxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/sstlxx.F	2011-10-11 00:36:19 +0000
@@ -13,8 +13,10 @@
 c       complex vertex         : amplitude               gamma(s1,s2,s3)
 c     
       implicit none
+      INTEGER DIM
+      PARAMETER(DIM=18)
 
-      include "dimension.inc"
+c      include "dimension.inc"
       double complex s1(DIM),s2(DIM),t3(DIM),vertex
       double precision gc
 

=== modified file 'HELAS/ttsaxx.F'
--- HELAS/ttsaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/ttsaxx.F	2011-10-11 00:36:19 +0000
@@ -17,9 +17,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex tc1(DIM),tc2(DIM),sc(DIM)
       double complex vertex, gc(2)
 

=== modified file 'HELAS/ttssxx.F'
--- HELAS/ttssxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/ttssxx.F	2011-10-11 00:36:19 +0000
@@ -18,9 +18,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex tc1(DIM),tc2(DIM),sc1(DIM),sc2(DIM)
 
       double complex vertex, g1(2)

=== modified file 'HELAS/usslxx.F'
--- HELAS/usslxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/usslxx.F	2011-10-11 00:36:19 +0000
@@ -17,7 +17,10 @@
 c
 
       implicit none
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+
+c      include "dimension.inc"
       double complex s1(DIM),s2(DIM),uss(DIM)
       double precision xm,xw,gc
 

=== modified file 'HELAS/utsaxx.F'
--- HELAS/utsaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/utsaxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex tc1(DIM),jts(DIM),sc(DIM), gt(2)
       double precision xm, xw
 

=== modified file 'HELAS/utssxx.F'
--- HELAS/utssxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/utssxx.F	2011-10-11 00:36:19 +0000
@@ -20,9 +20,9 @@
 
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex tc1(DIM),jts(DIM),sc1(DIM),sc2(DIM),g1(2)
       double precision g2, xm, xw
 

=== modified file 'HELAS/uvvaxx.F'
--- HELAS/uvvaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/uvvaxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
       implicit none
 
 c dimension of the current set to arbitrary length
-c      integer DIM
-c      parameter (DIM=18)
-      include "dimension.inc"
+      integer DIM
+      parameter (DIM=18)
+c      include "dimension.inc"
       double complex w1(DIM),w2(DIM),jt(DIM)
       double precision xm1,xm2,xw,g,s2g
       double precision sqrTwo

=== modified file 'HELAS/uvvvlx.F'
--- HELAS/uvvvlx.F	2010-10-30 03:26:37 +0000
+++ HELAS/uvvvlx.F	2011-10-11 00:36:19 +0000
@@ -22,9 +22,9 @@
       implicit none
 
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),gc(DIM),jhvvv(DIM),g2(2)
 
       double complex dvertx, vertex, vertex1, vertex2

=== modified file 'HELAS/vvshxx.F'
--- HELAS/vvshxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/vvshxx.F	2011-10-11 00:36:19 +0000
@@ -18,9 +18,9 @@
 c     
       implicit none
 c--   dimension of the current set to arbitrary length
-c     INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex v1(DIM),v2(DIM),sc(DIM),vertex,vertex1,vertex2
       double complex v12,p2v1,p1v2,v13,v14,v23,v24,v34
       double precision p12,p13,p14,p23,p24,p34

=== modified file 'HELAS/vvsshx.F'
--- HELAS/vvsshx.F	2010-10-30 03:26:37 +0000
+++ HELAS/vvsshx.F	2011-10-11 00:36:19 +0000
@@ -18,9 +18,9 @@
 c     
       implicit none
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex v1(DIM),v2(DIM),sc1(DIM),sc2(DIM),vertex
       double complex vertex1,vertex2,g1(2)
       double complex v12,p2v1,p1v2,v13,v14,v23,v24,v34

=== modified file 'HELAS/vvtaxx.F'
--- HELAS/vvtaxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/vvtaxx.F	2011-10-11 00:36:19 +0000
@@ -19,9 +19,9 @@
       implicit none
 
 c dimension of the current set to arbitrary length
-c      integer DIM
-c      parameter (DIM=18)
-      include "dimension.inc"
+      integer DIM
+      parameter (DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),tc(DIM)
       double precision xm,xw,g
       double precision sqrTwo

=== modified file 'HELAS/vvvsxx.F'
--- HELAS/vvvsxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/vvvsxx.F	2011-10-11 00:36:19 +0000
@@ -21,9 +21,9 @@
       implicit none
 
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),gc(DIM),sc(DIM)
 
       double complex dvertx, vertex, vertex1, vertex2

=== modified file 'HELAS/vvvtlx.F'
--- HELAS/vvvtlx.F	2010-10-30 03:26:37 +0000
+++ HELAS/vvvtlx.F	2011-10-11 00:36:19 +0000
@@ -20,9 +20,9 @@
       implicit none
 
 c--   dimension of the current set to arbitrary length
-c      INTEGER DIM
-c      PARAMETER(DIM=18)
-      include "dimension.inc"
+      INTEGER DIM
+      PARAMETER(DIM=18)
+c      include "dimension.inc"
       double complex ga(DIM),gb(DIM),gc(DIM),sc(DIM)
       double complex dvertx, vertex, vertex1, vertex2
       double complex vab, vbc, vca, v123, v124, v134, v234

=== modified file 'Template/bin/plot'
--- Template/bin/plot	2011-05-24 04:27:02 +0000
+++ Template/bin/plot	2011-10-11 00:36:19 +0000
@@ -20,6 +20,8 @@
 
 if [[ -x $td/td ]]; then
     echo 'topdrawer found'
+    # This is needed for Mac OS X
+    export DYLD_LIBRARY_PATH=$DYLD_LIBRARY_PATH:$td
 else
     echo 'topdrawer not found'
     exit

=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2011-10-10 07:34:37 +0000
+++ UpdateNotes.txt	2011-10-11 00:36:19 +0000
@@ -1,9 +1,21 @@
 Update notes for MadGraph 5 (in reverse time order)
 
+<<<<<<< TREE
 1.3.20 (09/10/11) JA: Fixed bug in myamp.f that created a hard stop
                       error for certain cases with many processes with
                       different propagators in the same subprocess dir.
 
+=======
+1.3.20 (07/10/11) OM: Add a new command: install. 
+                      This allow to install quite easily different package 
+                      devellop for Madgraph/MadEvent. The list of available 
+		      package are pythia-pgs/MadAnalysis/ExRootAnalysis/Delphes
+                  OM: Adding TopEffth Model
+                  OM: Improve display particles and autocompletion in presence
+  		      of nonpropagating particles
+                  OM: Fix Aloha bug linked to four fermion operator
+
+>>>>>>> MERGE-SOURCE
 1.3.19 (06/10/11) JA: Fixed problem with SubProcesses makefile on Linux,
                       introduced in 1.3.18.
 

=== modified file 'aloha/aloha_writers.py'
--- aloha/aloha_writers.py	2011-09-30 13:55:43 +0000
+++ aloha/aloha_writers.py	2011-10-11 00:36:19 +0000
@@ -517,11 +517,12 @@
             short_name, addon = name.split('C',1)
             if addon.split('_')[0].isdigit():
                 addon = 'C' +self.namestring.split('C',1)[1]
+            elif all([n.isdigit() for n in addon.split('_')[0].split('C')]):
+                addon = 'C' +self.namestring.split('C',1)[1]
             else:
                 addon = '_%s' % self.offshell
         else:
             addon = '_%s' % self.offshell
-
         # how to call the routine
         if not offshell:
             main = 'vertex'
@@ -907,6 +908,8 @@
             short_name, addon = name.split('C',1)
             if addon.split('_')[0].isdigit():
                 addon = 'C' +self.namestring.split('C',1)[1]
+            elif all([n.isdigit() for n in addon.split('_')[0].split('C')]):
+                addon = 'C' +self.namestring.split('C',1)[1]
             else:
                 addon = '_%s' % self.offshell
         else:
@@ -1239,6 +1242,8 @@
             short_name, addon = name.split('C',1)
             if addon.split('_')[0].isdigit():
                 addon = 'C' +self.namestring.split('C',1)[1]
+            elif all([n.isdigit() for n in addon.split('_')[0].split('C')]):
+                addon = 'C' +self.namestring.split('C',1)[1]
             else:
                 addon = '_%s' % self.offshell
         else:

=== modified file 'madgraph/VERSION'
--- madgraph/VERSION	2011-10-10 07:34:37 +0000
+++ madgraph/VERSION	2011-10-11 00:36:19 +0000
@@ -1,3 +1,8 @@
+<<<<<<< TREE
 version = 1.3.20
 date = 2011-10-09
+=======
+version = 1.3.20
+date = 2011-10-07
+>>>>>>> MERGE-SOURCE
 

=== modified file 'madgraph/interface/cmd_interface.py'
--- madgraph/interface/cmd_interface.py	2011-08-27 02:01:27 +0000
+++ madgraph/interface/cmd_interface.py	2011-10-11 00:36:19 +0000
@@ -275,6 +275,13 @@
         logger.info("")
         logger.info("   import command PATH :")
         logger.info("      Execute the list of command in the file at PATH")
+ 
+    def help_install(self):
+        logger.info("syntax: install " + "|".join(self._install_opts))
+        logger.info("-- Download the last version of the program and install it")
+        logger.info("   localy in the current Madgraph version. In order to have")
+        logger.info("   a sucessfull instalation, you will need to have up-to-date")
+        logger.info("   F77 and/or C and Root compiler.")
         
     def help_display(self):
         logger.info("syntax: display " + "|".join(self._display_opts))
@@ -634,6 +641,33 @@
             raise self.InvalidCmd('PATH is mandatory in the current context\n' + \
                                   'Did you forget to run the \"output\" command')
                         
+    def check_install(self, args):
+        """check that the install command is valid"""
+        
+        if len(args) != 1:
+            self.help_install()
+            raise self.InvalidCmd('install command require at least one argument')
+        
+        if args[0] not in self._install_opts:
+            if not args[0].startswith('td'):
+                self.help_install()
+                raise self.InvalidCmd('Not recognize program %s ' % args[0])
+            
+        if args[0] in ["ExRootAnalysis", "Delphes"]:
+            if not misc.which('root'):
+                raise self.InvalidCmd(
+'''In order to install ExRootAnalysis, you need to install Root on your computer first.
+please follow information on http://root.cern.ch/drupal/content/downloading-root''')
+            if 'ROOTSYS' not in os.environ:
+                raise self.InvalidCmd(
+'''The environment variable ROOTSYS is not configured.
+You can set it by adding the following lines in your .bashrc [.bash_profile for mac]:
+export ROOTSYS=%s
+export PATH=$PATH:$ROOTSYS/bin
+This will take effect only in a NEW terminal
+''' % os.path.realpath(os.path.join(misc.which('root'), \
+                                               os.path.pardir, os.path.pardir)))
+
                 
     def check_launch(self, args, options):
         """check the validity of the line"""
@@ -922,6 +956,10 @@
 
         CheckValidForCmd.check_import(self, args)
         
+    def check_install(self, args):
+        """ No possibility to install new software on the web """
+        raise self.WebRestriction('Impossible to install program on the cluster')
+        
     def check_load(self, args):
         """ check the validity of the line
         No Path authorize for the Web"""
@@ -988,6 +1026,7 @@
 
     def model_completion(self, text, process):
         """ complete the line with model information """
+
         while ',' in process:
             process = process[process.index(',')+1:]
         args = split_arg(process)
@@ -1421,6 +1460,14 @@
         # return
         return output
     
+    def complete_install(self, text, line, begidx, endidx):
+        "Complete the import command"
+
+        args = split_arg(line[0:begidx])
+        
+        # Format
+        if len(args) == 1:
+            return self.list_completion(text, self._install_opts)     
 #===============================================================================
 # MadGraphCmd
 #===============================================================================
@@ -1436,6 +1483,7 @@
     _tutorial_opts = ['start', 'stop']
     _check_opts = ['full', 'permutation', 'gauge', 'lorentz_invariance']
     _import_formats = ['model_v4', 'model', 'proc_v4', 'command']
+    _install_opts = ['pythia-pgs', 'Delphes', 'MadAnalysis', 'ExRootAnalysis']
     _v4_export_formats = ['madevent', 'standalone', 'matrix'] 
     _export_formats = _v4_export_formats + ['standalone_cpp', 'pythia8']
     _set_options = ['group_subprocesses',
@@ -1595,11 +1643,19 @@
             self.draw(' '.join(args[1:]))
 
         if args[0] == 'particles' and len(args) == 1:
+            propagating_particle = []
+            nb_unpropagating = 0
+            for particle in self._curr_model['particles']:
+                if particle.get('propagating'):
+                    propagating_particle.append(particle)
+                else:
+                    nb_unpropagating += 1
+                    
             print "Current model contains %i particles:" % \
-                    len(self._curr_model['particles'])
-            part_antipart = [part for part in self._curr_model['particles'] \
+                    len(propagating_particle)
+            part_antipart = [part for part in propagating_particle \
                              if not part['self_antipart']]
-            part_self = [part for part in self._curr_model['particles'] \
+            part_self = [part for part in propagating_particle \
                              if part['self_antipart']]
             for part in part_antipart:
                 print part['name'] + '/' + part['antiname'],
@@ -1607,6 +1663,9 @@
             for part in part_self:
                 print part['name'],
             print ''
+            if nb_unpropagating:
+                print 'In addition of %s un-physical particle mediating new interactions.' \
+                                     % nb_unpropagating
 
         elif args[0] == 'particles':
             for arg in args[1:]:
@@ -2316,6 +2375,7 @@
  
             #convert and excecute the card
             self.import_mg4_proc_card(proc_card)
+
     
     def import_ufo_model(self, model_name):
         """ import the UFO model """
@@ -2331,8 +2391,11 @@
         completion, define multiparticles"""
 
          # Set variables for autocomplete
-        self._particle_names = [p.get('name') for p in self._curr_model.get('particles')] + \
-             [p.get('antiname') for p in self._curr_model.get('particles')]
+        self._particle_names = [p.get('name') for p in self._curr_model.get('particles')\
+                                                    if p.get('propagating')] + \
+                 [p.get('antiname') for p in self._curr_model.get('particles') \
+                                                    if p.get('propagating')]
+        
         self._couplings = list(set(sum([i.get('orders').keys() for i in \
                                         self._curr_model.get('interactions')], [])))
         # Check if we can use case-independent particle names
@@ -2432,6 +2495,113 @@
         if removed_multiparticles:
             logger.info("Removed obsolete multiparticles %s" % \
                                          " / ".join(removed_multiparticles))
+
+    def do_install(self, line):
+        """Install optional package from the MG suite."""
+        
+        args = split_arg(line)
+        #check the validity of the arguments
+        self.check_install(args)
+
+        if sys.platform == "darwin":
+            program = "curl"
+        else:
+            program = "wget"
+            
+        # Load file with path of the different program:
+        import urllib
+        path = {}
+        try:
+            data = urllib.urlopen('http://madgraph.phys.ucl.ac.be/package_info.dat')
+        except:
+            raise MadGraph5Error, '''Impossible to connect the server. 
+            Please check your internet connection or retry later'''
+        for line in data: 
+            split = line.split()   
+            path[split[0]] = split[1]
+        
+        name = {'td_mac': 'td', 'td_linux':'td', 'Delphes':'Delphes', 
+                'pythia-pgs':'pythia-pgs', 'ExRootAnalysis': 'ExRootAnalysis',
+                'MadAnalysis':'MadAnalysis'}
+        name = name[args[0]]
+        
+        try:
+            os.system('rm -rf %s' % name)
+        except:
+            pass
+        
+        # Load that path
+        logger.info('Downloading %s' % path[args[0]])
+        if sys.platform == "darwin":
+            subprocess.call(['curl', path[args[0]], '-o%s.tgz' % name], cwd=MG5DIR)
+        else:
+            subprocess.call(['wget', path[args[0]], '--output-document=%s.tgz'% name], cwd=MG5DIR)
+        # Untar the file
+        returncode = subprocess.call(['tar', '-xzpvf', '%s.tgz' % name], cwd=MG5DIR)
+        if returncode:
+            raise MadGraph5Error, 'Fail to download correctly the File. Stop'
+        
+        # Check that the directory has the correct name
+        if not os.path.exists(os.path.join(MG5DIR, name)):
+            created_name = [n for n in os.listdir(MG5DIR) if n.startswith(name) 
+                                                  and not n.endswith('gz')]
+            if not created_name:
+                raise MadGraph5Error, 'The file was not loaded correctly. Stop'
+            else:
+                created_name = created_name[0]
+            files.mv(os.path.join(MG5DIR, created_name), os.path.join(MG5DIR, name))
+        logger.info('compile %s. This might takes a while.' % name)
+        
+        # Modify Makefile for pythia-pgs on Mac 64 bit
+        if args[0] == "pythia-pgs" and sys.maxsize > 2**32:
+            for path in [os.path.join(MG5DIR, 'pythia-pgs', 'libraries', \
+                         'PGS4', 'src', 'stdhep-dir', 'src', 'stdhep_Arch'),
+                         os.path.join(MG5DIR, 'pythia-pgs', 'libraries', \
+                         'PGS4', 'src', 'stdhep-dir', 'mcfio', 'arch_mcfio')]:
+                text = open(path).read()
+                text = text.replace('-m32','-m64')
+                open(path, 'w').writelines(text)
+            
+        # Compile the file
+        # Check for F77 compiler
+        if 'FC' not in os.environ or not os.environ['FC']:
+            if misc.which('g77'):
+                os.environ['FC'] = 'g77'
+            elif misc.which('gfortran'):
+                os.environ['FC'] = 'gfortran'
+            else:
+                raise self.InvalidCmd('Require g77 or Gfortran compiler')
+        subprocess.call(['make', 'clean'], cwd = os.path.join(MG5DIR, name))
+        subprocess.call(['make'], cwd = os.path.join(MG5DIR, name))
+
+
+        # Special treatment for TD program (require by MadAnalysis)
+        if args[0] == 'MadAnalysis':
+            try:
+                os.system('rm -rf td')
+                os.mkdir(os.path.join(MG5DIR, 'td'))
+            except Exception, error:
+                print error
+                pass
+            
+            if sys.platform == "darwin":
+                logger.info('Downloading TD for Mac')
+                target = 'http://theory.fnal.gov/people/parke/TD/td_mac_intel.tar.gz'
+                subprocess.call(['curl', target, '-otd.tgz'], 
+                                                  cwd=os.path.join(MG5DIR,'td'))      
+                subprocess.call(['tar', '-xzpvf', 'td.tgz'], 
+                                                  cwd=os.path.join(MG5DIR,'td'))
+                files.mv(MG5DIR + '/td/td_mac_intel',MG5DIR+'/td/td')
+            else:
+                logger.info('Downloading TD for Linux 32 bit')
+                target = 'http://cp3wks05.fynu.ucl.ac.be/twiki/pub/Software/TopDrawer/td'
+                subprocess.call(['wget', target], cwd=os.path.join(MG5DIR,'td'))      
+            
+                if sys.maxsize > 2**32:
+                    logger.warning('''td program (needed by MadAnalysis) is not compile for 64 bit computer
+                Please follow instruction in http://cp3wks05.fynu.ucl.ac.be/twiki/bin/view/Software/TopDrawer.''')
+
+
     
     def set_configuration(self, config_path=None):
         """ assign all configuration variable from file 

=== added directory 'models/TopEffTh'
=== added file 'models/TopEffTh/TopEffTh.log'
--- models/TopEffTh/TopEffTh.log	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/TopEffTh.log	2011-10-11 00:36:19 +0000
@@ -0,0 +1,79 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 11:59:51
+
+
+# 
+# This is the logfile for the model TopEffTh
+
+# Authors: C. Degrande
+# Model version: 1.1
+# Checking the Lagrangians
+   * All Lagrangians are ok.
+#
+# Particle definitions
+#
+
+   * No particles removed. All particles correspond to GenInt setup.
+
+# Automatically assigned PDG numbers
+   * Assigned PDG number 9000001 to particle ghA
+   * Assigned PDG number 9000002 to particle ghZ
+   * Assigned PDG number 9000003 to particle ghWp
+   * Assigned PDG number 9000004 to particle ghWm
+   * Assigned PDG number 9000005 to particle ghG
+   * Assigned PDG number 9000006 to particle Tri0
+   * Assigned PDG number 9000007 to particle Tri
+   * Assigned PDG number 9000008 to particle Tri80
+   * Assigned PDG number 9000009 to particle Tri8
+   * Assigned PDG number 9000010 to particle V8t
+   * Assigned PDG number 9000011 to particle V8Q
+
+
+# Compulsory PDG codes:
+   * Class SM leptons complete.
+   * Class SM neutrinos complete.
+   * Class SM quarks complete.
+   * Class SM gauge bosons complete.
+#
+# Parameter definitions
+#
+
+   * All parameters are ok.
+
+
+# Vertices
+   * Calling FeynmanRules for 1 Lagrangians.
+   * Number of classes vertices: 106
+   * Number of flavored vertices: 142
+   * Saved vertices in InterfaceRun[ 1 ].
+   * Checked QNumber conservation.
+      - Quantum number GhostNumber conserved in all vertices.
+      - Quantum number LeptonNumber conserved in all vertices.
+      - Quantum number Q conserved in all vertices.
+   * particles.py written.
+   * parameters.py written.
+#
+# Vertex definitions
+#
+
+   * 142 vertices written.
+   * vertices.py written.
+#
+# Lorentz structure definitions
+#
+
+   * 89 lorentz structures written.
+   * lorentz.py written.
+#
+# Coupling definitions
+#
+
+   * 131 couplings written.
+   * couplings.py written.
+#
+# Coupling order definitions
+#
+
+   * 0 couplings orders written.
+   * coupling_orders.py written.

=== added file 'models/TopEffTh/__init__.py'
--- models/TopEffTh/__init__.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/__init__.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,22 @@
+
+import particles
+import couplings
+import lorentz
+import parameters
+import vertices
+import coupling_orders
+import write_param_card
+
+
+all_particles = particles.all_particles
+all_vertices = vertices.all_vertices
+all_couplings = couplings.all_couplings
+all_lorentz = lorentz.all_lorentz
+all_parameters = parameters.all_parameters
+all_orders = coupling_orders.all_orders
+all_functions = function_library.all_functions
+
+
+__author__ = "C. Degrande"
+__version__ = "1.1"
+__email__ = "[e-mail]"

=== added file 'models/TopEffTh/coupling_orders.py'
--- models/TopEffTh/coupling_orders.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/coupling_orders.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,20 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:07
+
+
+from object_library import all_orders, CouplingOrder
+
+
+NP = CouplingOrder(name = 'NP',
+                   expansion_order = 2,
+                   hierarchy = 1)
+
+QCD = CouplingOrder(name = 'QCD',
+                    expansion_order = 99,
+                    hierarchy = 2)
+
+QED = CouplingOrder(name = 'QED',
+                    expansion_order = 99,
+                    hierarchy = 4)
+

=== added file 'models/TopEffTh/couplings.py'
--- models/TopEffTh/couplings.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/couplings.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,535 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:06
+
+
+from object_library import all_couplings, Coupling
+
+from function_library import complexconjugate, re, im, csc, sec, acsc, asec
+
+
+
+GC_1 = Coupling(name = 'GC_1',
+                value = '-(ee*complex(0,1))/3.',
+                order = {'QED':1})
+
+GC_2 = Coupling(name = 'GC_2',
+                value = '(2*ee*complex(0,1))/3.',
+                order = {'QED':1})
+
+GC_3 = Coupling(name = 'GC_3',
+                value = '-(ee*complex(0,1))',
+                order = {'QED':1})
+
+GC_4 = Coupling(name = 'GC_4',
+                value = '-G',
+                order = {'QCD':1})
+
+GC_5 = Coupling(name = 'GC_5',
+                value = 'complex(0,1)*G',
+                order = {'QCD':1})
+
+GC_6 = Coupling(name = 'GC_6',
+                value = 'complex(0,1)*G**2',
+                order = {'QCD':2})
+
+GC_7 = Coupling(name = 'GC_7',
+                value = 'complex(0,1)*g8Q',
+                order = {'NP':1})
+
+GC_8 = Coupling(name = 'GC_8',
+                value = 'complex(0,1)*g8Qd',
+                order = {'NP':1})
+
+GC_9 = Coupling(name = 'GC_9',
+                value = 'complex(0,1)*g8Qq',
+                order = {'NP':1})
+
+GC_10 = Coupling(name = 'GC_10',
+                 value = 'complex(0,1)*g8Qu',
+                 order = {'NP':1})
+
+GC_11 = Coupling(name = 'GC_11',
+                 value = 'complex(0,1)*g8t',
+                 order = {'NP':1})
+
+GC_12 = Coupling(name = 'GC_12',
+                 value = 'complex(0,1)*g8td',
+                 order = {'NP':1})
+
+GC_13 = Coupling(name = 'GC_13',
+                 value = 'complex(0,1)*g8tq',
+                 order = {'NP':1})
+
+GC_14 = Coupling(name = 'GC_14',
+                 value = 'complex(0,1)*g8tu',
+                 order = {'NP':1})
+
+GC_15 = Coupling(name = 'GC_15',
+                 value = '-(complex(0,1)*gT)',
+                 order = {'NP':1})
+
+GC_16 = Coupling(name = 'GC_16',
+                 value = 'complex(0,1)*gT',
+                 order = {'NP':1})
+
+GC_17 = Coupling(name = 'GC_17',
+                 value = 'complex(0,1)*gT*cmath.sqrt(2)',
+                 order = {'NP':1})
+
+GC_18 = Coupling(name = 'GC_18',
+                 value = '-(complex(0,1)*gT8)',
+                 order = {'NP':1})
+
+GC_19 = Coupling(name = 'GC_19',
+                 value = 'complex(0,1)*gT8',
+                 order = {'NP':1})
+
+GC_20 = Coupling(name = 'GC_20',
+                 value = 'complex(0,1)*gT8*cmath.sqrt(2)',
+                 order = {'NP':1})
+
+GC_21 = Coupling(name = 'GC_21',
+                 value = '-(complex(0,1)*gT8l)',
+                 order = {'NP':1})
+
+GC_22 = Coupling(name = 'GC_22',
+                 value = 'complex(0,1)*gT8l',
+                 order = {'NP':1})
+
+GC_23 = Coupling(name = 'GC_23',
+                 value = 'complex(0,1)*gT8l*cmath.sqrt(2)',
+                 order = {'NP':1})
+
+GC_24 = Coupling(name = 'GC_24',
+                 value = '-(complex(0,1)*gTl)',
+                 order = {'NP':1})
+
+GC_25 = Coupling(name = 'GC_25',
+                 value = 'complex(0,1)*gTl',
+                 order = {'NP':1})
+
+GC_26 = Coupling(name = 'GC_26',
+                 value = 'complex(0,1)*gTl*cmath.sqrt(2)',
+                 order = {'NP':1})
+
+GC_27 = Coupling(name = 'GC_27',
+                 value = 'cw*complex(0,1)*gw',
+                 order = {'QED':1})
+
+GC_28 = Coupling(name = 'GC_28',
+                 value = '-(complex(0,1)*gw**2)',
+                 order = {'QED':2})
+
+GC_29 = Coupling(name = 'GC_29',
+                 value = 'cw**2*complex(0,1)*gw**2',
+                 order = {'QED':2})
+
+GC_30 = Coupling(name = 'GC_30',
+                 value = '-6*complex(0,1)*lam',
+                 order = {'QED':2})
+
+GC_31 = Coupling(name = 'GC_31',
+                 value = '(-6*CG)/Lambda**2',
+                 order = {'NP':2})
+
+GC_32 = Coupling(name = 'GC_32',
+                 value = '(2*CphiG*complex(0,1))/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_33 = Coupling(name = 'GC_33',
+                 value = '(-2*CtG*complex(0,1)*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_34 = Coupling(name = 'GC_34',
+                 value = '(-2*CtW*complex(0,1))/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_35 = Coupling(name = 'GC_35',
+                 value = '-((CtW*cw*complex(0,1)*cmath.sqrt(2))/Lambda**2)',
+                 order = {'NP':2,'QED':1})
+
+GC_36 = Coupling(name = 'GC_36',
+                 value = '(6*CG*complex(0,1)*G)/Lambda**2',
+                 order = {'NP':2,'QCD':1})
+
+GC_37 = Coupling(name = 'GC_37',
+                 value = '(2*CphiG*G)/Lambda**2',
+                 order = {'NP':2,'QCD':1,'QED':1})
+
+GC_38 = Coupling(name = 'GC_38',
+                 value = '(-2*CtG*G*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2,'QCD':1,'QED':1})
+
+GC_39 = Coupling(name = 'GC_39',
+                 value = '(-3*CG*G**2)/Lambda**2',
+                 order = {'NP':2,'QCD':2})
+
+GC_40 = Coupling(name = 'GC_40',
+                 value = '(3*CG*G**2)/Lambda**2',
+                 order = {'NP':2,'QCD':2})
+
+GC_41 = Coupling(name = 'GC_41',
+                 value = '(-2*CphiG*complex(0,1)*G**2)/Lambda**2',
+                 order = {'NP':2,'QCD':2,'QED':1})
+
+GC_42 = Coupling(name = 'GC_42',
+                 value = '-((CG*complex(0,1)*G**3)/Lambda**2)',
+                 order = {'NP':2,'QCD':3})
+
+GC_43 = Coupling(name = 'GC_43',
+                 value = '(CG*complex(0,1)*G**3)/Lambda**2',
+                 order = {'NP':2,'QCD':3})
+
+GC_44 = Coupling(name = 'GC_44',
+                 value = '(CtW*complex(0,1)*gw*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2,'QED':2})
+
+GC_45 = Coupling(name = 'GC_45',
+                 value = '(2*CtW*cw*complex(0,1)*gw)/Lambda**2',
+                 order = {'NP':2,'QED':2})
+
+GC_46 = Coupling(name = 'GC_46',
+                 value = '(ee**2*complex(0,1))/(2.*sw**2)',
+                 order = {'QED':2})
+
+GC_47 = Coupling(name = 'GC_47',
+                 value = '(ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_48 = Coupling(name = 'GC_48',
+                 value = '(CKM11*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_49 = Coupling(name = 'GC_49',
+                 value = '(CKM12*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_50 = Coupling(name = 'GC_50',
+                 value = '(CKM13*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_51 = Coupling(name = 'GC_51',
+                 value = '(CKM21*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_52 = Coupling(name = 'GC_52',
+                 value = '(CKM22*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_53 = Coupling(name = 'GC_53',
+                 value = '(CKM23*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_54 = Coupling(name = 'GC_54',
+                 value = '(CKM31*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_55 = Coupling(name = 'GC_55',
+                 value = '(CKM32*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_56 = Coupling(name = 'GC_56',
+                 value = '(CKM33*ee*complex(0,1))/(sw*cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_57 = Coupling(name = 'GC_57',
+                 value = '-(cw*ee*complex(0,1))/(2.*sw)',
+                 order = {'QED':1})
+
+GC_58 = Coupling(name = 'GC_58',
+                 value = '(cw*ee*complex(0,1))/(2.*sw)',
+                 order = {'QED':1})
+
+GC_59 = Coupling(name = 'GC_59',
+                 value = '(C3phiq*ee*complex(0,1)*cmath.sqrt(2))/(Lambda**2*sw)',
+                 order = {'NP':2,'QED':2})
+
+GC_60 = Coupling(name = 'GC_60',
+                 value = '-(ee*complex(0,1)*sw)/(6.*cw)',
+                 order = {'QED':1})
+
+GC_61 = Coupling(name = 'GC_61',
+                 value = '(ee*complex(0,1)*sw)/(2.*cw)',
+                 order = {'QED':1})
+
+GC_62 = Coupling(name = 'GC_62',
+                 value = 'complex(0,1)*gw*sw',
+                 order = {'QED':1})
+
+GC_63 = Coupling(name = 'GC_63',
+                 value = '-2*cw*complex(0,1)*gw**2*sw',
+                 order = {'QED':2})
+
+GC_64 = Coupling(name = 'GC_64',
+                 value = '-((CtW*complex(0,1)*sw*cmath.sqrt(2))/Lambda**2)',
+                 order = {'NP':2,'QED':1})
+
+GC_65 = Coupling(name = 'GC_65',
+                 value = '(-2*CtW*complex(0,1)*gw*sw)/Lambda**2',
+                 order = {'NP':2,'QED':2})
+
+GC_66 = Coupling(name = 'GC_66',
+                 value = 'complex(0,1)*gw**2*sw**2',
+                 order = {'QED':2})
+
+GC_67 = Coupling(name = 'GC_67',
+                 value = '(cw*ee*complex(0,1))/(2.*sw) + (ee*complex(0,1)*sw)/(2.*cw)',
+                 order = {'QED':1})
+
+GC_68 = Coupling(name = 'GC_68',
+                 value = 'ee**2*complex(0,1) + (cw**2*ee**2*complex(0,1))/(2.*sw**2) + (ee**2*complex(0,1)*sw**2)/(2.*cw**2)',
+                 order = {'QED':2})
+
+GC_69 = Coupling(name = 'GC_69',
+                 value = '-6*complex(0,1)*lam*v',
+                 order = {'QED':1})
+
+GC_70 = Coupling(name = 'GC_70',
+                 value = '(2*CphiG*complex(0,1)*v)/Lambda**2',
+                 order = {'NP':2})
+
+GC_71 = Coupling(name = 'GC_71',
+                 value = '(-2*CtG*complex(0,1)*v*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2})
+
+GC_72 = Coupling(name = 'GC_72',
+                 value = '(-2*CtW*complex(0,1)*v)/Lambda**2',
+                 order = {'NP':2})
+
+GC_73 = Coupling(name = 'GC_73',
+                 value = '-((CtW*cw*complex(0,1)*v*cmath.sqrt(2))/Lambda**2)',
+                 order = {'NP':2})
+
+GC_74 = Coupling(name = 'GC_74',
+                 value = '(2*CphiG*G*v)/Lambda**2',
+                 order = {'NP':2,'QCD':1})
+
+GC_75 = Coupling(name = 'GC_75',
+                 value = '(-2*CtG*G*v*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2,'QCD':1})
+
+GC_76 = Coupling(name = 'GC_76',
+                 value = '(-2*CphiG*complex(0,1)*G**2*v)/Lambda**2',
+                 order = {'NP':2,'QCD':2})
+
+GC_77 = Coupling(name = 'GC_77',
+                 value = '(CtW*complex(0,1)*gw*v*cmath.sqrt(2))/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_78 = Coupling(name = 'GC_78',
+                 value = '(2*CtW*cw*complex(0,1)*gw*v)/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_79 = Coupling(name = 'GC_79',
+                 value = '(ee**2*complex(0,1)*v)/(2.*sw**2)',
+                 order = {'QED':1})
+
+GC_80 = Coupling(name = 'GC_80',
+                 value = '(C3phiq*ee*complex(0,1)*v*cmath.sqrt(2))/(Lambda**2*sw)',
+                 order = {'NP':2,'QED':1})
+
+GC_81 = Coupling(name = 'GC_81',
+                 value = '-((CtW*complex(0,1)*sw*v*cmath.sqrt(2))/Lambda**2)',
+                 order = {'NP':2})
+
+GC_82 = Coupling(name = 'GC_82',
+                 value = '(-2*CtW*complex(0,1)*gw*sw*v)/Lambda**2',
+                 order = {'NP':2,'QED':1})
+
+GC_83 = Coupling(name = 'GC_83',
+                 value = '(C3phiq*ee*complex(0,1)*v**2)/(Lambda**2*sw*cmath.sqrt(2))',
+                 order = {'NP':2})
+
+GC_84 = Coupling(name = 'GC_84',
+                 value = 'ee**2*complex(0,1)*v + (cw**2*ee**2*complex(0,1)*v)/(2.*sw**2) + (ee**2*complex(0,1)*sw**2*v)/(2.*cw**2)',
+                 order = {'QED':1})
+
+GC_85 = Coupling(name = 'GC_85',
+                 value = '-((complex(0,1)*yb)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_86 = Coupling(name = 'GC_86',
+                 value = '-((complex(0,1)*yc)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_87 = Coupling(name = 'GC_87',
+                 value = '-((complex(0,1)*ydo)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_88 = Coupling(name = 'GC_88',
+                 value = '-((complex(0,1)*ye)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_89 = Coupling(name = 'GC_89',
+                 value = '-((complex(0,1)*ym)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_90 = Coupling(name = 'GC_90',
+                 value = '-((complex(0,1)*ys)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_91 = Coupling(name = 'GC_91',
+                 value = '-((complex(0,1)*yt)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_92 = Coupling(name = 'GC_92',
+                 value = '-((complex(0,1)*ytau)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_93 = Coupling(name = 'GC_93',
+                 value = '-((complex(0,1)*yup)/cmath.sqrt(2))',
+                 order = {'QED':1})
+
+GC_94 = Coupling(name = 'GC_94',
+                 value = '(ee*complex(0,1)*complexconjugate(C3phiq)*cmath.sqrt(2))/(Lambda**2*sw)',
+                 order = {'NP':2,'QED':2})
+
+GC_95 = Coupling(name = 'GC_95',
+                 value = '(ee*complex(0,1)*v*complexconjugate(C3phiq)*cmath.sqrt(2))/(Lambda**2*sw)',
+                 order = {'NP':2,'QED':1})
+
+GC_96 = Coupling(name = 'GC_96',
+                 value = '(ee*complex(0,1)*v**2*complexconjugate(C3phiq))/(Lambda**2*sw*cmath.sqrt(2))',
+                 order = {'NP':2})
+
+GC_97 = Coupling(name = 'GC_97',
+                 value = '(C3phiq*complex(0,1))/(2.*Lambda**2) - (complex(0,1)*complexconjugate(C3phiq))/(2.*Lambda**2)',
+                 order = {'NP':2,'QED':1})
+
+GC_98 = Coupling(name = 'GC_98',
+                 value = '-(C3phiq*complex(0,1))/(2.*Lambda**2) + (complex(0,1)*complexconjugate(C3phiq))/(2.*Lambda**2)',
+                 order = {'NP':2,'QED':1})
+
+GC_99 = Coupling(name = 'GC_99',
+                 value = '-(C3phiq*cw*ee*complex(0,1))/(2.*Lambda**2*sw) - (C3phiq*ee*complex(0,1)*sw)/(2.*cw*Lambda**2) - (cw*ee*complex(0,1)*complexconjugate(C3phiq))/(2.*Lambda**2*sw) - (ee*complex(0,1)*sw*complexconjugate(C3phiq))/(2.*cw*Lambda**2)',
+                 order = {'NP':2,'QED':2})
+
+GC_100 = Coupling(name = 'GC_100',
+                  value = '(C3phiq*cw*ee*complex(0,1))/(2.*Lambda**2*sw) + (C3phiq*ee*complex(0,1)*sw)/(2.*cw*Lambda**2) + (cw*ee*complex(0,1)*complexconjugate(C3phiq))/(2.*Lambda**2*sw) + (ee*complex(0,1)*sw*complexconjugate(C3phiq))/(2.*cw*Lambda**2)',
+                  order = {'NP':2,'QED':2})
+
+GC_101 = Coupling(name = 'GC_101',
+                  value = '(C3phiq*complex(0,1)*v)/(2.*Lambda**2) - (complex(0,1)*v*complexconjugate(C3phiq))/(2.*Lambda**2)',
+                  order = {'NP':2})
+
+GC_102 = Coupling(name = 'GC_102',
+                  value = '-(C3phiq*complex(0,1)*v)/(2.*Lambda**2) + (complex(0,1)*v*complexconjugate(C3phiq))/(2.*Lambda**2)',
+                  order = {'NP':2})
+
+GC_103 = Coupling(name = 'GC_103',
+                  value = '-(C3phiq*cw*ee*complex(0,1)*v)/(2.*Lambda**2*sw) - (C3phiq*ee*complex(0,1)*sw*v)/(2.*cw*Lambda**2) - (cw*ee*complex(0,1)*v*complexconjugate(C3phiq))/(2.*Lambda**2*sw) - (ee*complex(0,1)*sw*v*complexconjugate(C3phiq))/(2.*cw*Lambda**2)',
+                  order = {'NP':2,'QED':1})
+
+GC_104 = Coupling(name = 'GC_104',
+                  value = '(C3phiq*cw*ee*complex(0,1)*v)/(2.*Lambda**2*sw) + (C3phiq*ee*complex(0,1)*sw*v)/(2.*cw*Lambda**2) + (cw*ee*complex(0,1)*v*complexconjugate(C3phiq))/(2.*Lambda**2*sw) + (ee*complex(0,1)*sw*v*complexconjugate(C3phiq))/(2.*cw*Lambda**2)',
+                  order = {'NP':2,'QED':1})
+
+GC_105 = Coupling(name = 'GC_105',
+                  value = '-(C3phiq*cw*ee*complex(0,1)*v**2)/(4.*Lambda**2*sw) - (C3phiq*ee*complex(0,1)*sw*v**2)/(4.*cw*Lambda**2) - (cw*ee*complex(0,1)*v**2*complexconjugate(C3phiq))/(4.*Lambda**2*sw) - (ee*complex(0,1)*sw*v**2*complexconjugate(C3phiq))/(4.*cw*Lambda**2)',
+                  order = {'NP':2})
+
+GC_106 = Coupling(name = 'GC_106',
+                  value = '(C3phiq*cw*ee*complex(0,1)*v**2)/(4.*Lambda**2*sw) + (C3phiq*ee*complex(0,1)*sw*v**2)/(4.*cw*Lambda**2) + (cw*ee*complex(0,1)*v**2*complexconjugate(C3phiq))/(4.*Lambda**2*sw) + (ee*complex(0,1)*sw*v**2*complexconjugate(C3phiq))/(4.*cw*Lambda**2)',
+                  order = {'NP':2})
+
+GC_107 = Coupling(name = 'GC_107',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM11))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_108 = Coupling(name = 'GC_108',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM12))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_109 = Coupling(name = 'GC_109',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM13))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_110 = Coupling(name = 'GC_110',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM21))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_111 = Coupling(name = 'GC_111',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM22))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_112 = Coupling(name = 'GC_112',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM23))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_113 = Coupling(name = 'GC_113',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM31))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_114 = Coupling(name = 'GC_114',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM32))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_115 = Coupling(name = 'GC_115',
+                  value = '(ee*complex(0,1)*complexconjugate(CKM33))/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
+
+GC_116 = Coupling(name = 'GC_116',
+                  value = '(-2*complex(0,1)*complexconjugate(CtG)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2,'QED':1})
+
+GC_117 = Coupling(name = 'GC_117',
+                  value = '(-2*G*complexconjugate(CtG)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2,'QCD':1,'QED':1})
+
+GC_118 = Coupling(name = 'GC_118',
+                  value = '(-2*complex(0,1)*v*complexconjugate(CtG)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2})
+
+GC_119 = Coupling(name = 'GC_119',
+                  value = '(-2*G*v*complexconjugate(CtG)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2,'QCD':1})
+
+GC_120 = Coupling(name = 'GC_120',
+                  value = '(-2*complex(0,1)*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2,'QED':1})
+
+GC_121 = Coupling(name = 'GC_121',
+                  value = '-((cw*complex(0,1)*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2)',
+                  order = {'NP':2,'QED':1})
+
+GC_122 = Coupling(name = 'GC_122',
+                  value = '(complex(0,1)*gw*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2,'QED':2})
+
+GC_123 = Coupling(name = 'GC_123',
+                  value = '(-2*cw*complex(0,1)*gw*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2,'QED':2})
+
+GC_124 = Coupling(name = 'GC_124',
+                  value = '-((complex(0,1)*sw*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2)',
+                  order = {'NP':2,'QED':1})
+
+GC_125 = Coupling(name = 'GC_125',
+                  value = '(2*complex(0,1)*gw*sw*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2,'QED':2})
+
+GC_126 = Coupling(name = 'GC_126',
+                  value = '(-2*complex(0,1)*v*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2})
+
+GC_127 = Coupling(name = 'GC_127',
+                  value = '-((cw*complex(0,1)*v*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2)',
+                  order = {'NP':2})
+
+GC_128 = Coupling(name = 'GC_128',
+                  value = '(complex(0,1)*gw*v*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2',
+                  order = {'NP':2,'QED':1})
+
+GC_129 = Coupling(name = 'GC_129',
+                  value = '(-2*cw*complex(0,1)*gw*v*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2,'QED':1})
+
+GC_130 = Coupling(name = 'GC_130',
+                  value = '-((complex(0,1)*sw*v*complexconjugate(CtW)*cmath.sqrt(2))/Lambda**2)',
+                  order = {'NP':2})
+
+GC_131 = Coupling(name = 'GC_131',
+                  value = '(2*complex(0,1)*gw*sw*v*complexconjugate(CtW))/Lambda**2',
+                  order = {'NP':2,'QED':1})
+

=== added file 'models/TopEffTh/function_library.py'
--- models/TopEffTh/function_library.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/function_library.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,54 @@
+# This file is part of the UFO.
+#
+# This file contains definitions for functions that
+# are extensions of the cmath library, and correspond
+# either to functions that are in cmath, but inconvenient
+# to access from there (e.g. z.conjugate()),
+# or functions that are simply not defined.
+#
+#
+
+__date__ = "22 July 2010"
+__author__ = "[e-mail]"
+
+import cmath
+from object_library import all_functions, Function
+
+#
+# shortcuts for functions from cmath
+#
+
+complexconjugate = Function(name = 'complexconjugate',
+                            arguments = ('z',),
+                            expression = 'z.conjugate()')
+
+
+re = Function(name = 're',
+              arguments = ('z',),
+              expression = 'z.real')
+
+im = Function(name = 'im',
+              arguments = ('z',),
+              expression = 'z.imag')
+
+# New functions (trigonometric)
+
+sec = Function(name = 'sec',
+             arguments = ('z',),
+             expression = '1./cmath.cos(z)')
+
+asec = Function(name = 'asec',
+             arguments = ('z',),
+             expression = 'cmath.acos(1./z)')
+
+csc = Function(name = 'csc',
+             arguments = ('z',),
+             expression = '1./cmath.sin(z)')
+
+acsc = Function(name = 'acsc',
+             arguments = ('z',),
+             expression = 'cmath.asin(1./z)')
+
+
+
+

=== added file 'models/TopEffTh/lorentz.py'
--- models/TopEffTh/lorentz.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/lorentz.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,367 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:06
+
+
+from object_library import all_lorentz, Lorentz
+
+from function_library import complexconjugate, re, im, csc, sec, acsc, asec
+
+
+
+UUV1 = Lorentz(name = 'UUV1',
+               spins = [ -1, -1, 3 ],
+               structure = 'P(3,2) + P(3,3)')
+
+SSS1 = Lorentz(name = 'SSS1',
+               spins = [ 1, 1, 1 ],
+               structure = '1')
+
+FFS1 = Lorentz(name = 'FFS1',
+               spins = [ 2, 2, 1 ],
+               structure = 'Identity(2,1)')
+
+FFS2 = Lorentz(name = 'FFS2',
+               spins = [ 2, 2, 1 ],
+               structure = 'P(-1,3)*Gamma(-1,2,-2)*ProjM(-2,1)')
+
+FFV1 = Lorentz(name = 'FFV1',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,1)')
+
+FFV2 = Lorentz(name = 'FFV2',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
+
+FFV3 = Lorentz(name = 'FFV3',
+               spins = [ 2, 2, 3 ],
+               structure = '-(P(-1,3)*Gamma(-1,-3,-2)*Gamma(3,2,-3)*ProjM(-2,1)) + P(3,3)*ProjM(2,1)')
+
+FFV4 = Lorentz(name = 'FFV4',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,-1)*ProjP(-1,1)')
+
+FFV5 = Lorentz(name = 'FFV5',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,-1)*ProjM(-1,1) - 2*Gamma(3,2,-1)*ProjP(-1,1)')
+
+FFV6 = Lorentz(name = 'FFV6',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 2*Gamma(3,2,-1)*ProjP(-1,1)')
+
+FFV7 = Lorentz(name = 'FFV7',
+               spins = [ 2, 2, 3 ],
+               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 4*Gamma(3,2,-1)*ProjP(-1,1)')
+
+FFV8 = Lorentz(name = 'FFV8',
+               spins = [ 2, 2, 3 ],
+               structure = '-(P(-1,3)*Gamma(-1,-3,-2)*Gamma(3,2,-3)*ProjP(-2,1)) + P(3,3)*ProjP(2,1)')
+
+VVS1 = Lorentz(name = 'VVS1',
+               spins = [ 3, 3, 1 ],
+               structure = 'Metric(1,2)')
+
+VVS2 = Lorentz(name = 'VVS2',
+               spins = [ 3, 3, 1 ],
+               structure = 'P(1,2)*P(2,1) - P(-1,1)*P(-1,2)*Metric(1,2)')
+
+VVV1 = Lorentz(name = 'VVV1',
+               spins = [ 3, 3, 3 ],
+               structure = 'P(3,1)*Metric(1,2) - P(3,2)*Metric(1,2) - P(2,1)*Metric(1,3) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3) - P(1,3)*Metric(2,3)')
+
+VVV2 = Lorentz(name = 'VVV2',
+               spins = [ 3, 3, 3 ],
+               structure = '-(P(1,2)*P(2,3)*P(3,1)) + P(1,3)*P(2,1)*P(3,2) + P(-1,2)*P(-1,3)*P(3,1)*Metric(1,2) - P(-1,1)*P(-1,3)*P(3,2)*Metric(1,2) - P(-1,2)*P(-1,3)*P(2,1)*Metric(1,3) + P(-1,1)*P(-1,2)*P(2,3)*Metric(1,3) + P(-1,1)*P(-1,3)*P(1,2)*Metric(2,3) - P(-1,1)*P(-1,2)*P(1,3)*Metric(2,3)')
+
+SSSS1 = Lorentz(name = 'SSSS1',
+                spins = [ 1, 1, 1, 1 ],
+                structure = '1')
+
+FFSS1 = Lorentz(name = 'FFSS1',
+                spins = [ 2, 2, 1, 1 ],
+                structure = 'P(-1,3)*Gamma(-1,2,-2)*ProjM(-2,1) + P(-1,4)*Gamma(-1,2,-2)*ProjM(-2,1)')
+
+FFVS1 = Lorentz(name = 'FFVS1',
+                spins = [ 2, 2, 3, 1 ],
+                structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
+
+FFVS2 = Lorentz(name = 'FFVS2',
+                spins = [ 2, 2, 3, 1 ],
+                structure = '-(P(-1,3)*Gamma(-1,-3,-2)*Gamma(3,2,-3)*ProjM(-2,1)) + P(3,3)*ProjM(2,1)')
+
+FFVS3 = Lorentz(name = 'FFVS3',
+                spins = [ 2, 2, 3, 1 ],
+                structure = '-(P(-1,3)*Gamma(-1,-3,-2)*Gamma(3,2,-3)*ProjP(-2,1)) + P(3,3)*ProjP(2,1)')
+
+FFVV1 = Lorentz(name = 'FFVV1',
+                spins = [ 2, 2, 3, 3 ],
+                structure = 'Gamma(3,2,-2)*Gamma(4,-2,-1)*ProjM(-1,1) - Metric(3,4)*ProjM(2,1)')
+
+FFVV2 = Lorentz(name = 'FFVV2',
+                spins = [ 2, 2, 3, 3 ],
+                structure = 'Gamma(3,2,-2)*Gamma(4,-2,-1)*ProjP(-1,1) - Metric(3,4)*ProjP(2,1)')
+
+VVSS1 = Lorentz(name = 'VVSS1',
+                spins = [ 3, 3, 1, 1 ],
+                structure = 'Metric(1,2)')
+
+VVSS2 = Lorentz(name = 'VVSS2',
+                spins = [ 3, 3, 1, 1 ],
+                structure = 'P(1,2)*P(2,1) - P(-1,1)*P(-1,2)*Metric(1,2)')
+
+VVVS1 = Lorentz(name = 'VVVS1',
+                spins = [ 3, 3, 3, 1 ],
+                structure = 'P(3,1)*Metric(1,2) - P(3,2)*Metric(1,2) - P(2,1)*Metric(1,3) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3) - P(1,3)*Metric(2,3)')
+
+VVVV1 = Lorentz(name = 'VVVV1',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
+
+VVVV2 = Lorentz(name = 'VVVV2',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'P(3,2)*P(4,1)*Metric(1,2) - P(3,1)*P(4,2)*Metric(1,2) + P(2,1)*P(4,2)*Metric(1,3) + P(2,4)*P(4,3)*Metric(1,3) - P(2,1)*P(3,2)*Metric(1,4) - P(2,3)*P(3,4)*Metric(1,4) - P(1,2)*P(4,1)*Metric(2,3) - P(1,4)*P(4,3)*Metric(2,3) + P(-1,1)*P(-1,2)*Metric(1,4)*Metric(2,3) + P(-1,3)*P(-1,4)*Metric(1,4)*Metric(2,3) + P(1,2)*P(3,1)*Metric(2,4) + P(1,3)*P(3,4)*Metric(2,4) - P(-1,1)*P(-1,2)*Metric(1,3)*Metric(2,4) - P(-1,3)*P(-1,4)*Metric(1,3)*Metric(2,4) + P(1,4)*P(2,3)*Metric(3,4) - P(1,3)*P(2,4)*Metric(3,4)')
+
+VVVV3 = Lorentz(name = 'VVVV3',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) - 2*Metric(1,2)*Metric(3,4)')
+
+VVVV4 = Lorentz(name = 'VVVV4',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
+
+VVVV5 = Lorentz(name = 'VVVV5',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
+
+VVVV6 = Lorentz(name = 'VVVV6',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'Metric(1,4)*Metric(2,3) - (Metric(1,3)*Metric(2,4))/2. - (Metric(1,2)*Metric(3,4))/2.')
+
+VVVV7 = Lorentz(name = 'VVVV7',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'P(3,4)*P(4,1)*Metric(1,2) + P(3,2)*P(4,3)*Metric(1,2) - P(2,4)*P(4,1)*Metric(1,3) - P(2,3)*P(4,2)*Metric(1,3) + P(2,4)*P(3,1)*Metric(1,4) - P(2,1)*P(3,4)*Metric(1,4) + P(1,3)*P(4,2)*Metric(2,3) - P(1,2)*P(4,3)*Metric(2,3) - P(1,4)*P(3,1)*Metric(2,4) - P(1,3)*P(3,2)*Metric(2,4) + P(-1,2)*P(-1,3)*Metric(1,3)*Metric(2,4) + P(-1,1)*P(-1,4)*Metric(1,3)*Metric(2,4) + P(1,4)*P(2,1)*Metric(3,4) + P(1,2)*P(2,3)*Metric(3,4) - P(-1,2)*P(-1,3)*Metric(1,2)*Metric(3,4) - P(-1,1)*P(-1,4)*Metric(1,2)*Metric(3,4)')
+
+VVVV8 = Lorentz(name = 'VVVV8',
+                spins = [ 3, 3, 3, 3 ],
+                structure = 'P(3,4)*P(4,2)*Metric(1,2) + P(3,1)*P(4,3)*Metric(1,2) + P(2,3)*P(4,1)*Metric(1,3) - P(2,1)*P(4,3)*Metric(1,3) - P(2,3)*P(3,1)*Metric(1,4) - P(2,4)*P(3,2)*Metric(1,4) - P(1,3)*P(4,1)*Metric(2,3) - P(1,4)*P(4,2)*Metric(2,3) + P(-1,1)*P(-1,3)*Metric(1,4)*Metric(2,3) + P(-1,2)*P(-1,4)*Metric(1,4)*Metric(2,3) + P(1,4)*P(3,2)*Metric(2,4) - P(1,2)*P(3,4)*Metric(2,4) + P(1,3)*P(2,1)*Metric(3,4) + P(1,2)*P(2,4)*Metric(3,4) - P(-1,1)*P(-1,3)*Metric(1,2)*Metric(3,4) - P(-1,2)*P(-1,4)*Metric(1,2)*Metric(3,4)')
+
+FFVSS1 = Lorentz(name = 'FFVSS1',
+                 spins = [ 2, 2, 3, 1, 1 ],
+                 structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
+
+FFVVS1 = Lorentz(name = 'FFVVS1',
+                 spins = [ 2, 2, 3, 3, 1 ],
+                 structure = 'Gamma(3,2,-2)*Gamma(4,-2,-1)*ProjM(-1,1) - Metric(3,4)*ProjM(2,1)')
+
+FFVVS2 = Lorentz(name = 'FFVVS2',
+                 spins = [ 2, 2, 3, 3, 1 ],
+                 structure = 'Gamma(3,2,-2)*Gamma(4,-2,-1)*ProjP(-1,1) - Metric(3,4)*ProjP(2,1)')
+
+VVVSS1 = Lorentz(name = 'VVVSS1',
+                 spins = [ 3, 3, 3, 1, 1 ],
+                 structure = 'P(3,1)*Metric(1,2) - P(3,2)*Metric(1,2) - P(2,1)*Metric(1,3) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3) - P(1,3)*Metric(2,3)')
+
+VVVVS1 = Lorentz(name = 'VVVVS1',
+                 spins = [ 3, 3, 3, 3, 1 ],
+                 structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
+
+VVVVS2 = Lorentz(name = 'VVVVS2',
+                 spins = [ 3, 3, 3, 3, 1 ],
+                 structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
+
+VVVVS3 = Lorentz(name = 'VVVVS3',
+                 spins = [ 3, 3, 3, 3, 1 ],
+                 structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
+
+VVVVV1 = Lorentz(name = 'VVVVV1',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(5,3)*Metric(1,4)*Metric(2,3) - P(4,3)*Metric(1,5)*Metric(2,3) - P(5,3)*Metric(1,3)*Metric(2,4) + P(4,3)*Metric(1,3)*Metric(2,5)')
+
+VVVVV2 = Lorentz(name = 'VVVVV2',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(5,4)*Metric(1,4)*Metric(2,3) - P(5,4)*Metric(1,3)*Metric(2,4) + P(3,4)*Metric(1,5)*Metric(2,4) - P(3,4)*Metric(1,4)*Metric(2,5)')
+
+VVVVV3 = Lorentz(name = 'VVVVV3',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(4,5)*Metric(1,5)*Metric(2,3) - P(3,5)*Metric(1,5)*Metric(2,4) - P(4,5)*Metric(1,3)*Metric(2,5) + P(3,5)*Metric(1,4)*Metric(2,5)')
+
+VVVVV4 = Lorentz(name = 'VVVVV4',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(5,1)*Metric(1,3)*Metric(2,4) - P(4,1)*Metric(1,3)*Metric(2,5) - P(5,1)*Metric(1,2)*Metric(3,4) + P(4,1)*Metric(1,2)*Metric(3,5)')
+
+VVVVV5 = Lorentz(name = 'VVVVV5',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(5,2)*Metric(1,4)*Metric(2,3) - P(4,2)*Metric(1,5)*Metric(2,3) - P(5,2)*Metric(1,2)*Metric(3,4) + P(4,2)*Metric(1,2)*Metric(3,5)')
+
+VVVVV6 = Lorentz(name = 'VVVVV6',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(3,1)*Metric(1,5)*Metric(2,4) - P(3,1)*Metric(1,4)*Metric(2,5) - P(2,1)*Metric(1,5)*Metric(3,4) + P(2,1)*Metric(1,4)*Metric(3,5)')
+
+VVVVV7 = Lorentz(name = 'VVVVV7',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(5,4)*Metric(1,4)*Metric(2,3) - P(5,4)*Metric(1,2)*Metric(3,4) + P(2,4)*Metric(1,5)*Metric(3,4) - P(2,4)*Metric(1,4)*Metric(3,5)')
+
+VVVVV8 = Lorentz(name = 'VVVVV8',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(4,5)*Metric(1,5)*Metric(2,3) - P(2,5)*Metric(1,5)*Metric(3,4) - P(4,5)*Metric(1,2)*Metric(3,5) + P(2,5)*Metric(1,4)*Metric(3,5)')
+
+VVVVV9 = Lorentz(name = 'VVVVV9',
+                 spins = [ 3, 3, 3, 3, 3 ],
+                 structure = 'P(3,2)*Metric(1,5)*Metric(2,4) - P(3,2)*Metric(1,4)*Metric(2,5) + P(1,2)*Metric(2,5)*Metric(3,4) - P(1,2)*Metric(2,4)*Metric(3,5)')
+
+VVVVV10 = Lorentz(name = 'VVVVV10',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(2,3)*Metric(1,5)*Metric(3,4) - P(1,3)*Metric(2,5)*Metric(3,4) - P(2,3)*Metric(1,4)*Metric(3,5) + P(1,3)*Metric(2,4)*Metric(3,5)')
+
+VVVVV11 = Lorentz(name = 'VVVVV11',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,4)*Metric(1,3)*Metric(2,4) - P(5,4)*Metric(1,2)*Metric(3,4) + P(1,4)*Metric(2,5)*Metric(3,4) - P(1,4)*Metric(2,4)*Metric(3,5)')
+
+VVVVV12 = Lorentz(name = 'VVVVV12',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,5)*Metric(1,3)*Metric(2,5) - P(1,5)*Metric(2,5)*Metric(3,4) - P(4,5)*Metric(1,2)*Metric(3,5) + P(1,5)*Metric(2,4)*Metric(3,5)')
+
+VVVVV13 = Lorentz(name = 'VVVVV13',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,1)*Metric(1,4)*Metric(2,3) - P(3,1)*Metric(1,4)*Metric(2,5) - P(5,1)*Metric(1,2)*Metric(3,4) + P(3,1)*Metric(1,2)*Metric(4,5)')
+
+VVVVV14 = Lorentz(name = 'VVVVV14',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,1)*Metric(1,5)*Metric(2,3) - P(3,1)*Metric(1,5)*Metric(2,4) - P(4,1)*Metric(1,2)*Metric(3,5) + P(3,1)*Metric(1,2)*Metric(4,5)')
+
+VVVVV15 = Lorentz(name = 'VVVVV15',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,2)*Metric(1,3)*Metric(2,4) - P(3,2)*Metric(1,5)*Metric(2,4) - P(5,2)*Metric(1,2)*Metric(3,4) + P(3,2)*Metric(1,2)*Metric(4,5)')
+
+VVVVV16 = Lorentz(name = 'VVVVV16',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,2)*Metric(1,3)*Metric(2,5) - P(3,2)*Metric(1,4)*Metric(2,5) - P(4,2)*Metric(1,2)*Metric(3,5) + P(3,2)*Metric(1,2)*Metric(4,5)')
+
+VVVVV17 = Lorentz(name = 'VVVVV17',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,1)*Metric(1,5)*Metric(2,3) - P(4,1)*Metric(1,3)*Metric(2,5) - P(2,1)*Metric(1,5)*Metric(3,4) + P(2,1)*Metric(1,3)*Metric(4,5)')
+
+VVVVV18 = Lorentz(name = 'VVVVV18',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,1)*Metric(1,4)*Metric(2,3) - P(5,1)*Metric(1,3)*Metric(2,4) - P(2,1)*Metric(1,4)*Metric(3,5) + P(2,1)*Metric(1,3)*Metric(4,5)')
+
+VVVVV19 = Lorentz(name = 'VVVVV19',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,3)*Metric(1,3)*Metric(2,4) - P(5,3)*Metric(1,2)*Metric(3,4) + P(2,3)*Metric(1,5)*Metric(3,4) - P(2,3)*Metric(1,3)*Metric(4,5)')
+
+VVVVV20 = Lorentz(name = 'VVVVV20',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,3)*Metric(1,3)*Metric(2,5) - P(4,3)*Metric(1,2)*Metric(3,5) + P(2,3)*Metric(1,4)*Metric(3,5) - P(2,3)*Metric(1,3)*Metric(4,5)')
+
+VVVVV21 = Lorentz(name = 'VVVVV21',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(3,4)*Metric(1,4)*Metric(2,5) - P(2,4)*Metric(1,4)*Metric(3,5) - P(3,4)*Metric(1,2)*Metric(4,5) + P(2,4)*Metric(1,3)*Metric(4,5)')
+
+VVVVV22 = Lorentz(name = 'VVVVV22',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(3,5)*Metric(1,5)*Metric(2,4) - P(2,5)*Metric(1,5)*Metric(3,4) - P(3,5)*Metric(1,2)*Metric(4,5) + P(2,5)*Metric(1,3)*Metric(4,5)')
+
+VVVVV23 = Lorentz(name = 'VVVVV23',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,2)*Metric(1,5)*Metric(2,3) - P(4,2)*Metric(1,3)*Metric(2,5) + P(1,2)*Metric(2,5)*Metric(3,4) - P(1,2)*Metric(2,3)*Metric(4,5)')
+
+VVVVV24 = Lorentz(name = 'VVVVV24',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,2)*Metric(1,4)*Metric(2,3) - P(5,2)*Metric(1,3)*Metric(2,4) + P(1,2)*Metric(2,4)*Metric(3,5) - P(1,2)*Metric(2,3)*Metric(4,5)')
+
+VVVVV25 = Lorentz(name = 'VVVVV25',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(5,3)*Metric(1,4)*Metric(2,3) - P(5,3)*Metric(1,2)*Metric(3,4) + P(1,3)*Metric(2,5)*Metric(3,4) - P(1,3)*Metric(2,3)*Metric(4,5)')
+
+VVVVV26 = Lorentz(name = 'VVVVV26',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(4,3)*Metric(1,5)*Metric(2,3) - P(4,3)*Metric(1,2)*Metric(3,5) + P(1,3)*Metric(2,4)*Metric(3,5) - P(1,3)*Metric(2,3)*Metric(4,5)')
+
+VVVVV27 = Lorentz(name = 'VVVVV27',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(3,4)*Metric(1,5)*Metric(2,4) - P(1,4)*Metric(2,4)*Metric(3,5) - P(3,4)*Metric(1,2)*Metric(4,5) + P(1,4)*Metric(2,3)*Metric(4,5)')
+
+VVVVV28 = Lorentz(name = 'VVVVV28',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(2,4)*Metric(1,5)*Metric(3,4) - P(1,4)*Metric(2,5)*Metric(3,4) - P(2,4)*Metric(1,3)*Metric(4,5) + P(1,4)*Metric(2,3)*Metric(4,5)')
+
+VVVVV29 = Lorentz(name = 'VVVVV29',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(3,5)*Metric(1,4)*Metric(2,5) - P(1,5)*Metric(2,5)*Metric(3,4) - P(3,5)*Metric(1,2)*Metric(4,5) + P(1,5)*Metric(2,3)*Metric(4,5)')
+
+VVVVV30 = Lorentz(name = 'VVVVV30',
+                  spins = [ 3, 3, 3, 3, 3 ],
+                  structure = 'P(2,5)*Metric(1,4)*Metric(3,5) - P(1,5)*Metric(2,4)*Metric(3,5) - P(2,5)*Metric(1,3)*Metric(4,5) + P(1,5)*Metric(2,3)*Metric(4,5)')
+
+VVVVSS1 = Lorentz(name = 'VVVVSS1',
+                  spins = [ 3, 3, 3, 3, 1, 1 ],
+                  structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
+
+VVVVSS2 = Lorentz(name = 'VVVVSS2',
+                  spins = [ 3, 3, 3, 3, 1, 1 ],
+                  structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
+
+VVVVSS3 = Lorentz(name = 'VVVVSS3',
+                  spins = [ 3, 3, 3, 3, 1, 1 ],
+                  structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
+
+VVVVVV1 = Lorentz(name = 'VVVVVV1',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,4)*Metric(3,5) - Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,5)*Metric(2,4)*Metric(3,6) + Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,3)*Metric(2,6)*Metric(4,5) + Metric(1,5)*Metric(2,3)*Metric(4,6) - Metric(1,3)*Metric(2,5)*Metric(4,6)')
+
+VVVVVV2 = Lorentz(name = 'VVVVVV2',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,5)*Metric(2,6)*Metric(3,4) + Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,5)*Metric(2,3)*Metric(4,6) - Metric(1,2)*Metric(3,5)*Metric(4,6)')
+
+VVVVVV3 = Lorentz(name = 'VVVVVV3',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,5)*Metric(2,6)*Metric(3,4) - Metric(1,6)*Metric(2,4)*Metric(3,5) + Metric(1,5)*Metric(2,4)*Metric(3,6) + Metric(1,3)*Metric(2,6)*Metric(4,5) - Metric(1,2)*Metric(3,6)*Metric(4,5) - Metric(1,3)*Metric(2,5)*Metric(4,6) + Metric(1,2)*Metric(3,5)*Metric(4,6)')
+
+VVVVVV4 = Lorentz(name = 'VVVVVV4',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,5)*Metric(2,6)*Metric(3,4) + Metric(1,5)*Metric(2,4)*Metric(3,6) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,3)*Metric(2,6)*Metric(4,5) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,3)*Metric(2,4)*Metric(5,6)')
+
+VVVVVV5 = Lorentz(name = 'VVVVVV5',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,4)*Metric(3,5) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,3)*Metric(2,5)*Metric(4,6) - Metric(1,2)*Metric(3,5)*Metric(4,6) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,3)*Metric(2,4)*Metric(5,6)')
+
+VVVVVV6 = Lorentz(name = 'VVVVVV6',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,5)*Metric(2,6)*Metric(3,4) - Metric(1,6)*Metric(2,4)*Metric(3,5) + Metric(1,4)*Metric(2,6)*Metric(3,5) + Metric(1,5)*Metric(2,3)*Metric(4,6) - Metric(1,3)*Metric(2,5)*Metric(4,6) - Metric(1,4)*Metric(2,3)*Metric(5,6) + Metric(1,3)*Metric(2,4)*Metric(5,6)')
+
+VVVVVV7 = Lorentz(name = 'VVVVVV7',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,5)*Metric(2,4)*Metric(3,6) - Metric(1,3)*Metric(2,6)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,5)*Metric(2,3)*Metric(4,6) - Metric(1,2)*Metric(3,5)*Metric(4,6) - Metric(1,4)*Metric(2,3)*Metric(5,6) + Metric(1,3)*Metric(2,4)*Metric(5,6)')
+
+VVVVVV8 = Lorentz(name = 'VVVVVV8',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,5)*Metric(2,6)*Metric(3,4) + Metric(1,6)*Metric(2,4)*Metric(3,5) - Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,5)*Metric(2,4)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV9 = Lorentz(name = 'VVVVVV9',
+                  spins = [ 3, 3, 3, 3, 3, 3 ],
+                  structure = 'Metric(1,5)*Metric(2,6)*Metric(3,4) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,3)*Metric(2,6)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) - Metric(1,5)*Metric(2,3)*Metric(4,6) + Metric(1,3)*Metric(2,5)*Metric(4,6) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV10 = Lorentz(name = 'VVVVVV10',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,6)*Metric(2,4)*Metric(3,5) + Metric(1,5)*Metric(2,4)*Metric(3,6) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,5)*Metric(2,3)*Metric(4,6) + Metric(1,2)*Metric(3,5)*Metric(4,6) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV11 = Lorentz(name = 'VVVVVV11',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,3)*Metric(2,6)*Metric(4,5) - Metric(1,3)*Metric(2,5)*Metric(4,6) + Metric(1,2)*Metric(3,5)*Metric(4,6) + Metric(1,4)*Metric(2,3)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV12 = Lorentz(name = 'VVVVVV12',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,6)*Metric(2,4)*Metric(3,5) + Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,3)*Metric(2,6)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,3)*Metric(2,4)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV13 = Lorentz(name = 'VVVVVV13',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,6)*Metric(2,5)*Metric(3,4) - Metric(1,5)*Metric(2,4)*Metric(3,6) - Metric(1,6)*Metric(2,3)*Metric(4,5) + Metric(1,2)*Metric(3,6)*Metric(4,5) + Metric(1,5)*Metric(2,3)*Metric(4,6) - Metric(1,3)*Metric(2,5)*Metric(4,6) + Metric(1,3)*Metric(2,4)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV14 = Lorentz(name = 'VVVVVV14',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,5)*Metric(2,6)*Metric(3,4) - Metric(1,6)*Metric(2,4)*Metric(3,5) + Metric(1,6)*Metric(2,3)*Metric(4,5) - Metric(1,3)*Metric(2,6)*Metric(4,5) - Metric(1,5)*Metric(2,3)*Metric(4,6) + Metric(1,2)*Metric(3,5)*Metric(4,6) + Metric(1,3)*Metric(2,4)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+
+VVVVVV15 = Lorentz(name = 'VVVVVV15',
+                   spins = [ 3, 3, 3, 3, 3, 3 ],
+                   structure = 'Metric(1,5)*Metric(2,6)*Metric(3,4) - Metric(1,4)*Metric(2,6)*Metric(3,5) - Metric(1,5)*Metric(2,4)*Metric(3,6) + Metric(1,4)*Metric(2,5)*Metric(3,6) - Metric(1,3)*Metric(2,5)*Metric(4,6) + Metric(1,2)*Metric(3,5)*Metric(4,6) + Metric(1,3)*Metric(2,4)*Metric(5,6) - Metric(1,2)*Metric(3,4)*Metric(5,6)')
+

=== added file 'models/TopEffTh/note.pdf'
Binary files models/TopEffTh/note.pdf	1970-01-01 00:00:00 +0000 and models/TopEffTh/note.pdf	2011-10-11 00:36:19 +0000 differ
=== added file 'models/TopEffTh/object_library.py'
--- models/TopEffTh/object_library.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/object_library.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,245 @@
+##
+##
+## Feynrules Header
+##
+##
+##
+##
+##
+
+import cmath
+
+
+class UFOBaseClass(object):
+    """The class from which all FeynRules classes are derived."""
+
+    require_args = []
+
+    def __init__(self, *args, **options):
+        assert(len(self.require_args) == len (args))
+    
+        for i, name in enumerate(self.require_args):
+            setattr(self, name, args[i])
+    
+        for (option, value) in options.items():
+            setattr(self, option, value)
+
+    def get(self, name):
+        return getattr(self, name)
+    
+    def set(self, name, value):
+        setattr(self, name, value)
+        
+    def get_all(self):
+        """Return a dictionary containing all the information of the object"""
+        return self.__dict__
+
+    def __str__(self):
+        return self.name
+
+    def nice_string(self):
+        """ return string with the full information """
+        return '\n'.join(['%s \t: %s' %(name, value) for name, value in self.__dict__.items()])
+
+    def __repr__(self):
+        replacements = [
+            ('+','__plus__'),
+            ('-','__minus__'),
+            ('@','__at__'),
+            ('!','__exclam__'),
+            ('?','__quest__'),
+            ('*','__star__'),
+            ('~','__tilde__')
+            ]
+        text = self.name
+        for orig,sub in replacements:
+            text = text.replace(orig,sub)
+        return text
+
+
+
+all_particles = []
+
+    
+
+class Particle(UFOBaseClass):
+    """A standard Particle"""
+
+    require_args=['pdg_code', 'name', 'antiname', 'spin', 'color', 'mass', 'width', 'texname', 'antitexname', 'charge']
+
+    require_args_all = ['pdg_code', 'name', 'antiname', 'spin', 'color', 'mass', 'width', 'texname', 'antitexname', 'charge', 'line', 'propagating', 'goldstoneboson']
+
+    def __init__(self, pdg_code, name, antiname, spin, color, mass, width, texname,
+                 antitexname, charge , line=None, propagating=True, goldstoneboson=False, **options):
+
+        args= (pdg_code, name, antiname, spin, color, mass, width, texname,
+                 antitexname, float(charge))
+
+        UFOBaseClass.__init__(self, *args,  **options)
+
+        global all_particles
+        all_particles.append(self)
+
+        self.propagating = propagating
+        self.goldstoneboson= goldstoneboson
+
+        self.selfconjugate = (name == antiname)
+        if 1: #not line:
+            self.line = self.find_line_type()
+        else:
+            self.line = line
+
+
+
+
+    def find_line_type(self):
+        """ find how we draw a line if not defined
+        valid output: dashed/straight/wavy/curly/double/swavy/scurly
+        """
+        
+        spin = self.spin
+        color = self.color
+        
+        #use default
+        if spin == 1:
+            return 'dashed'
+        elif spin == 2:
+            if not self.selfconjugate:
+                return 'straight'
+            elif color == 1:
+                return 'swavy'
+            else:
+                return 'scurly'
+        elif spin == 3:
+            if color == 1:
+                return 'wavy'
+            
+            else:
+                return 'curly'
+        elif spin == 5:
+            return 'double'
+        elif spin == -1:
+            return 'dotted'
+        else:
+            return 'dashed' # not supported yet
+
+    def anti(self):
+        if self.selfconjugate:
+            raise Exception('%s has no anti particle.' % self.name) 
+        outdic = {}
+        for k,v in self.__dict__.iteritems():
+            if k not in self.require_args_all:                
+                outdic[k] = -v
+        if self.color in [1,8]:
+            newcolor = self.color
+        else:
+            newcolor = -self.color
+                
+        return Particle(-self.pdg_code, self.antiname, self.name, self.spin, newcolor, self.mass, self.width,
+                        self.antitexname, self.texname, -self.charge, self.line, self.propagating, self.goldstoneboson, **outdic)
+
+
+
+all_parameters = []
+
+class Parameter(UFOBaseClass):
+
+    require_args=['name', 'nature', 'type', 'value', 'texname']
+
+    def __init__(self, name, nature, type, value, texname, lhablock=None, lhacode=None):
+
+        args = (name,nature,type,value,texname)
+
+        UFOBaseClass.__init__(self, *args)
+
+        args=(name,nature,type,value,texname)
+
+        global all_parameters
+        all_parameters.append(self)
+
+        if (lhablock is None or lhacode is None)  and nature == 'external':
+            raise Exception('Need LHA information for external parameter "%s".' % name)
+        self.lhablock = lhablock
+        self.lhacode = lhacode
+
+all_vertices = []
+
+class Vertex(UFOBaseClass):
+
+    require_args=['name', 'particles', 'color', 'lorentz', 'couplings']
+
+    def __init__(self, name, particles, color, lorentz, couplings, **opt):
+ 
+        args = (name, particles, color, lorentz, couplings)
+
+        UFOBaseClass.__init__(self, *args, **opt)
+
+        args=(particles,color,lorentz,couplings)
+
+        global all_vertices
+        all_vertices.append(self)
+
+all_couplings = []
+
+class Coupling(UFOBaseClass):
+
+    require_args=['name', 'value', 'order']
+
+    def __init__(self, name, value, order, **opt):
+
+        args =(name, value, order)	
+        UFOBaseClass.__init__(self, *args, **opt)
+        global all_couplings
+        all_couplings.append(self)
+  
+
+
+all_lorentz = []
+
+class Lorentz(UFOBaseClass):
+
+    require_args=['name','spins','structure']
+    
+    def __init__(self, name, spins, structure='external', **opt):
+        args = (name, spins, structure)
+        UFOBaseClass.__init__(self, *args, **opt)
+
+        global all_lorentz
+        all_lorentz.append(self)
+
+
+all_functions = []
+
+class Function(object):
+
+    def __init__(self, name, arguments, expression):
+
+        global all_functions
+        all_functions.append(self)
+
+        self.name = name
+        self.arguments = arguments
+        self.expr = expression
+    
+    def __call__(self, *opt):
+
+        for i, arg in enumerate(self.arguments):
+            exec('%s = %s' % (arg, opt[i] ))
+
+        return eval(self.expr)
+
+all_orders = []
+
+class CouplingOrder(object):
+
+    def __init__(self, name, expansion_order, hierarchy, perturbative_expansion = 0):
+        
+        global all_orders
+        all_orders.append(self)
+
+        self.name = name
+        self.expansion_order = expansion_order
+        self.hierarchy = hierarchy
+
+
+        

=== added file 'models/TopEffTh/parameters.py'
--- models/TopEffTh/parameters.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/parameters.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,705 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:06
+
+
+
+from object_library import all_parameters, Parameter
+
+
+from function_library import complexconjugate, re, im, csc, sec, acsc, asec
+
+# This is a default parameter object representing 0.
+ZERO = Parameter(name = 'ZERO',
+                 nature = 'internal',
+                 type = 'real',
+                 value = '0.0',
+                 texname = '0')
+
+# User-defined parameters.
+K = Parameter(name = 'K',
+              nature = 'external',
+              type = 'real',
+              value = 100,
+              texname = 'K',
+              lhablock = 'APPROX',
+              lhacode = [ 1 ])
+
+K1 = Parameter(name = 'K1',
+               nature = 'external',
+               type = 'real',
+               value = 100,
+               texname = '\\text{K1}',
+               lhablock = 'APPROX',
+               lhacode = [ 2 ])
+
+K2 = Parameter(name = 'K2',
+               nature = 'external',
+               type = 'real',
+               value = 100,
+               texname = '\\text{K2}',
+               lhablock = 'APPROX',
+               lhacode = [ 3 ])
+
+K3 = Parameter(name = 'K3',
+               nature = 'external',
+               type = 'real',
+               value = 100,
+               texname = '\\text{K3}',
+               lhablock = 'APPROX',
+               lhacode = [ 4 ])
+
+cabi = Parameter(name = 'cabi',
+                 nature = 'external',
+                 type = 'real',
+                 value = 0.227736,
+                 texname = '\\theta _c',
+                 lhablock = 'CKMBLOCK',
+                 lhacode = [ 1 ])
+
+Lambda = Parameter(name = 'Lambda',
+                   nature = 'external',
+                   type = 'real',
+                   value = 1000,
+                   texname = '\\Lambda ',
+                   lhablock = 'DIM6',
+                   lhacode = [ 1 ])
+
+RC3phiq = Parameter(name = 'RC3phiq',
+                    nature = 'external',
+                    type = 'real',
+                    value = 1,
+                    texname = '\\text{RC}_{\\text{$\\phi $q}}^{\\text{(3)}}',
+                    lhablock = 'DIM6',
+                    lhacode = [ 2 ])
+
+IC3phiq = Parameter(name = 'IC3phiq',
+                    nature = 'external',
+                    type = 'real',
+                    value = 1,
+                    texname = '\\text{IC}_{\\text{$\\phi $q}}^{\\text{(3)}}',
+                    lhablock = 'DIM6',
+                    lhacode = [ 3 ])
+
+RCtW = Parameter(name = 'RCtW',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = '\\text{RC}_{\\text{tW}}',
+                 lhablock = 'DIM6',
+                 lhacode = [ 4 ])
+
+ICtW = Parameter(name = 'ICtW',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = '\\text{IC}_{\\text{tW}}',
+                 lhablock = 'DIM6',
+                 lhacode = [ 5 ])
+
+RCtG = Parameter(name = 'RCtG',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = '\\text{RC}_{\\text{tG}}',
+                 lhablock = 'DIM6',
+                 lhacode = [ 6 ])
+
+ICtG = Parameter(name = 'ICtG',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = '\\text{IC}_{\\text{tG}}',
+                 lhablock = 'DIM6',
+                 lhacode = [ 7 ])
+
+CG = Parameter(name = 'CG',
+               nature = 'external',
+               type = 'real',
+               value = 1,
+               texname = 'C_G',
+               lhablock = 'DIM6',
+               lhacode = [ 8 ])
+
+CphiG = Parameter(name = 'CphiG',
+                  nature = 'external',
+                  type = 'real',
+                  value = 1,
+                  texname = 'C_{\\text{$\\phi $G}}',
+                  lhablock = 'DIM6',
+                  lhacode = [ 9 ])
+
+C13qq = Parameter(name = 'C13qq',
+                  nature = 'external',
+                  type = 'real',
+                  value = 1,
+                  texname = 'C_{\\text{qq}}^{\\text{(1,3)}}',
+                  lhablock = 'FourFermion',
+                  lhacode = [ 1 ])
+
+C81qq = Parameter(name = 'C81qq',
+                  nature = 'external',
+                  type = 'real',
+                  value = 1,
+                  texname = 'C_{\\text{qq}}^{\\text{(8,1)}}',
+                  lhablock = 'FourFermion',
+                  lhacode = [ 2 ])
+
+C83qq = Parameter(name = 'C83qq',
+                  nature = 'external',
+                  type = 'real',
+                  value = 1,
+                  texname = 'C_{\\text{qq}}^{\\text{(8,3)}}',
+                  lhablock = 'FourFermion',
+                  lhacode = [ 3 ])
+
+C8ut = Parameter(name = 'C8ut',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = 'C_{\\text{ut}}^{\\text{(8)}}',
+                 lhablock = 'FourFermion',
+                 lhacode = [ 4 ])
+
+C8dt = Parameter(name = 'C8dt',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = 'C_{\\text{dt}}^{\\text{(8)}}',
+                 lhablock = 'FourFermion',
+                 lhacode = [ 5 ])
+
+C1qu = Parameter(name = 'C1qu',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = 'C_{\\text{qu}}^{\\text{(1)}}',
+                 lhablock = 'FourFermion',
+                 lhacode = [ 6 ])
+
+C1qd = Parameter(name = 'C1qd',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = 'C_{\\text{qd}}^{\\text{(1)}}',
+                 lhablock = 'FourFermion',
+                 lhacode = [ 7 ])
+
+C1qt = Parameter(name = 'C1qt',
+                 nature = 'external',
+                 type = 'real',
+                 value = 1,
+                 texname = 'C_{\\text{qt}}^{\\text{(1)}}',
+                 lhablock = 'FourFermion',
+                 lhacode = [ 8 ])
+
+aEWM1 = Parameter(name = 'aEWM1',
+                  nature = 'external',
+                  type = 'real',
+                  value = 127.9,
+                  texname = '\\text{aEWM1}',
+                  lhablock = 'SMINPUTS',
+                  lhacode = [ 1 ])
+
+Gf = Parameter(name = 'Gf',
+               nature = 'external',
+               type = 'real',
+               value = 0.0000116637,
+               texname = 'G_f',
+               lhablock = 'SMINPUTS',
+               lhacode = [ 2 ])
+
+aS = Parameter(name = 'aS',
+               nature = 'external',
+               type = 'real',
+               value = 0.1184,
+               texname = '\\text{aS}',
+               lhablock = 'SMINPUTS',
+               lhacode = [ 3 ])
+
+ymdo = Parameter(name = 'ymdo',
+                 nature = 'external',
+                 type = 'real',
+                 value = 0.00504,
+                 texname = '\\text{ymdo}',
+                 lhablock = 'YUKAWA',
+                 lhacode = [ 1 ])
+
+ymup = Parameter(name = 'ymup',
+                 nature = 'external',
+                 type = 'real',
+                 value = 0.0025499999999999997,
+                 texname = '\\text{ymup}',
+                 lhablock = 'YUKAWA',
+                 lhacode = [ 2 ])
+
+yms = Parameter(name = 'yms',
+                nature = 'external',
+                type = 'real',
+                value = 0.101,
+                texname = '\\text{yms}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 3 ])
+
+ymc = Parameter(name = 'ymc',
+                nature = 'external',
+                type = 'real',
+                value = 1.27,
+                texname = '\\text{ymc}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 4 ])
+
+ymb = Parameter(name = 'ymb',
+                nature = 'external',
+                type = 'real',
+                value = 4.7,
+                texname = '\\text{ymb}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 5 ])
+
+ymt = Parameter(name = 'ymt',
+                nature = 'external',
+                type = 'real',
+                value = 172.,
+                texname = '\\text{ymt}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 6 ])
+
+yme = Parameter(name = 'yme',
+                nature = 'external',
+                type = 'real',
+                value = 0.0005110000000000001,
+                texname = '\\text{yme}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 11 ])
+
+ymm = Parameter(name = 'ymm',
+                nature = 'external',
+                type = 'real',
+                value = 0.10566,
+                texname = '\\text{ymm}',
+                lhablock = 'YUKAWA',
+                lhacode = [ 13 ])
+
+ymtau = Parameter(name = 'ymtau',
+                  nature = 'external',
+                  type = 'real',
+                  value = 1.777,
+                  texname = '\\text{ymtau}',
+                  lhablock = 'YUKAWA',
+                  lhacode = [ 15 ])
+
+Me = Parameter(name = 'Me',
+               nature = 'external',
+               type = 'real',
+               value = 0.0005110000000000001,
+               texname = '\\text{Me}',
+               lhablock = 'MASS',
+               lhacode = [ 11 ])
+
+MM = Parameter(name = 'MM',
+               nature = 'external',
+               type = 'real',
+               value = 0.10566,
+               texname = '\\text{MM}',
+               lhablock = 'MASS',
+               lhacode = [ 13 ])
+
+MTA = Parameter(name = 'MTA',
+                nature = 'external',
+                type = 'real',
+                value = 1.777,
+                texname = '\\text{MTA}',
+                lhablock = 'MASS',
+                lhacode = [ 15 ])
+
+MU = Parameter(name = 'MU',
+               nature = 'external',
+               type = 'real',
+               value = 0.0025499999999999997,
+               texname = 'M',
+               lhablock = 'MASS',
+               lhacode = [ 2 ])
+
+MC = Parameter(name = 'MC',
+               nature = 'external',
+               type = 'real',
+               value = 1.42,
+               texname = '\\text{MC}',
+               lhablock = 'MASS',
+               lhacode = [ 4 ])
+
+MT = Parameter(name = 'MT',
+               nature = 'external',
+               type = 'real',
+               value = 172,
+               texname = '\\text{MT}',
+               lhablock = 'MASS',
+               lhacode = [ 6 ])
+
+MD = Parameter(name = 'MD',
+               nature = 'external',
+               type = 'real',
+               value = 0.00504,
+               texname = '\\text{MD}',
+               lhablock = 'MASS',
+               lhacode = [ 1 ])
+
+MS = Parameter(name = 'MS',
+               nature = 'external',
+               type = 'real',
+               value = 0.101,
+               texname = '\\text{MS}',
+               lhablock = 'MASS',
+               lhacode = [ 3 ])
+
+MB = Parameter(name = 'MB',
+               nature = 'external',
+               type = 'real',
+               value = 4.7,
+               texname = '\\text{MB}',
+               lhablock = 'MASS',
+               lhacode = [ 5 ])
+
+MZ = Parameter(name = 'MZ',
+               nature = 'external',
+               type = 'real',
+               value = 91.1876,
+               texname = '\\text{MZ}',
+               lhablock = 'MASS',
+               lhacode = [ 23 ])
+
+MH = Parameter(name = 'MH',
+               nature = 'external',
+               type = 'real',
+               value = 120,
+               texname = '\\text{MH}',
+               lhablock = 'MASS',
+               lhacode = [ 25 ])
+
+WT = Parameter(name = 'WT',
+               nature = 'external',
+               type = 'real',
+               value = 1.50833649,
+               texname = '\\text{WT}',
+               lhablock = 'DECAY',
+               lhacode = [ 6 ])
+
+WZ = Parameter(name = 'WZ',
+               nature = 'external',
+               type = 'real',
+               value = 2.4952,
+               texname = '\\text{WZ}',
+               lhablock = 'DECAY',
+               lhacode = [ 23 ])
+
+WW = Parameter(name = 'WW',
+               nature = 'external',
+               type = 'real',
+               value = 2.085,
+               texname = '\\text{WW}',
+               lhablock = 'DECAY',
+               lhacode = [ 24 ])
+
+WH = Parameter(name = 'WH',
+               nature = 'external',
+               type = 'real',
+               value = 0.00575308848,
+               texname = '\\text{WH}',
+               lhablock = 'DECAY',
+               lhacode = [ 25 ])
+
+aEW = Parameter(name = 'aEW',
+                nature = 'internal',
+                type = 'real',
+                value = '1/aEWM1',
+                texname = '\\text{aEW}')
+
+G = Parameter(name = 'G',
+              nature = 'internal',
+              type = 'real',
+              value = '2*cmath.sqrt(aS)*cmath.sqrt(cmath.pi)',
+              texname = 'G')
+
+CKM11 = Parameter(name = 'CKM11',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = 'cmath.cos(cabi)',
+                  texname = '\\text{CKM11}')
+
+CKM12 = Parameter(name = 'CKM12',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = 'cmath.sin(cabi)',
+                  texname = '\\text{CKM12}')
+
+CKM13 = Parameter(name = 'CKM13',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '0',
+                  texname = '\\text{CKM13}')
+
+CKM21 = Parameter(name = 'CKM21',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '-cmath.sin(cabi)',
+                  texname = '\\text{CKM21}')
+
+CKM22 = Parameter(name = 'CKM22',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = 'cmath.cos(cabi)',
+                  texname = '\\text{CKM22}')
+
+CKM23 = Parameter(name = 'CKM23',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '0',
+                  texname = '\\text{CKM23}')
+
+CKM31 = Parameter(name = 'CKM31',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '0',
+                  texname = '\\text{CKM31}')
+
+CKM32 = Parameter(name = 'CKM32',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '0',
+                  texname = '\\text{CKM32}')
+
+CKM33 = Parameter(name = 'CKM33',
+                  nature = 'internal',
+                  type = 'complex',
+                  value = '1',
+                  texname = '\\text{CKM33}')
+
+C3phiq = Parameter(name = 'C3phiq',
+                   nature = 'internal',
+                   type = 'complex',
+                   value = 'complex(0,1)*IC3phiq + RC3phiq',
+                   texname = 'C_{\\text{$\\phi $q}}^{\\text{(3)}}')
+
+CtW = Parameter(name = 'CtW',
+                nature = 'internal',
+                type = 'complex',
+                value = 'complex(0,1)*ICtW + RCtW',
+                texname = 'C_{\\text{tW}}')
+
+CtG = Parameter(name = 'CtG',
+                nature = 'internal',
+                type = 'complex',
+                value = 'complex(0,1)*ICtG + RCtG',
+                texname = 'C_{\\text{tG}}')
+
+MTri = Parameter(name = 'MTri',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'K*Lambda',
+                 texname = 'M_{\\text{Tri}}')
+
+MTri8 = Parameter(name = 'MTri8',
+                  nature = 'internal',
+                  type = 'real',
+                  value = 'K1*Lambda',
+                  texname = 'M_{\\text{Tri8}}')
+
+M8t = Parameter(name = 'M8t',
+                nature = 'internal',
+                type = 'real',
+                value = 'K2*Lambda',
+                texname = 'M_{8 t}')
+
+M8Q = Parameter(name = 'M8Q',
+                nature = 'internal',
+                type = 'real',
+                value = 'K3*Lambda',
+                texname = 'M_{8 Q}')
+
+gT = Parameter(name = 'gT',
+               nature = 'internal',
+               type = 'real',
+               value = '-(K*cmath.sqrt(C13qq))',
+               texname = 'g_T')
+
+gTl = Parameter(name = 'gTl',
+                nature = 'internal',
+                type = 'real',
+                value = 'K*cmath.sqrt(C13qq)',
+                texname = 'g_{\\text{Tl}}')
+
+gT8 = Parameter(name = 'gT8',
+                nature = 'internal',
+                type = 'real',
+                value = '-(K1*cmath.sqrt(C83qq))',
+                texname = 'g_{\\text{T8}}')
+
+gT8l = Parameter(name = 'gT8l',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'K1*cmath.sqrt(C83qq)',
+                 texname = 'g_{\\text{T8l}}')
+
+g8t = Parameter(name = 'g8t',
+                nature = 'internal',
+                type = 'real',
+                value = '-K2',
+                texname = 'g_{8 t}')
+
+g8tu = Parameter(name = 'g8tu',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C8ut*K2',
+                 texname = 'g_{8 \\text{tu}}')
+
+g8td = Parameter(name = 'g8td',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C8dt*K2',
+                 texname = 'g_{8 \\text{td}}')
+
+g8tq = Parameter(name = 'g8tq',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C1qt*K2',
+                 texname = 'g_{8 \\text{tq}}')
+
+g8Q = Parameter(name = 'g8Q',
+                nature = 'internal',
+                type = 'real',
+                value = '-K3',
+                texname = 'g_{8 Q}')
+
+g8Qu = Parameter(name = 'g8Qu',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C1qu*K3',
+                 texname = 'g_{8 \\text{Qu}}')
+
+g8Qd = Parameter(name = 'g8Qd',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C1qd*K3',
+                 texname = 'g_{8 \\text{Qd}}')
+
+g8Qq = Parameter(name = 'g8Qq',
+                 nature = 'internal',
+                 type = 'real',
+                 value = 'C81qq*K3',
+                 texname = 'g_{8 \\text{Qq}}')
+
+MW = Parameter(name = 'MW',
+               nature = 'internal',
+               type = 'real',
+               value = 'cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2))))',
+               texname = 'M_W')
+
+ee = Parameter(name = 'ee',
+               nature = 'internal',
+               type = 'real',
+               value = '2*cmath.sqrt(aEW)*cmath.sqrt(cmath.pi)',
+               texname = 'e')
+
+sw2 = Parameter(name = 'sw2',
+                nature = 'internal',
+                type = 'real',
+                value = '1 - MW**2/MZ**2',
+                texname = '\\text{sw2}')
+
+cw = Parameter(name = 'cw',
+               nature = 'internal',
+               type = 'real',
+               value = 'cmath.sqrt(1 - sw2)',
+               texname = 'c_w')
+
+sw = Parameter(name = 'sw',
+               nature = 'internal',
+               type = 'real',
+               value = 'cmath.sqrt(sw2)',
+               texname = 's_w')
+
+g1 = Parameter(name = 'g1',
+               nature = 'internal',
+               type = 'real',
+               value = 'ee/cw',
+               texname = 'g_1')
+
+gw = Parameter(name = 'gw',
+               nature = 'internal',
+               type = 'real',
+               value = 'ee/sw',
+               texname = 'g_w')
+
+v = Parameter(name = 'v',
+              nature = 'internal',
+              type = 'real',
+              value = '(2*MW*sw)/ee',
+              texname = 'v')
+
+lam = Parameter(name = 'lam',
+                nature = 'internal',
+                type = 'real',
+                value = 'MH**2/(2.*v**2)',
+                texname = '\\text{lam}')
+
+yb = Parameter(name = 'yb',
+               nature = 'internal',
+               type = 'real',
+               value = '(ymb*cmath.sqrt(2))/v',
+               texname = '\\text{yb}')
+
+yc = Parameter(name = 'yc',
+               nature = 'internal',
+               type = 'real',
+               value = '(ymc*cmath.sqrt(2))/v',
+               texname = '\\text{yc}')
+
+ydo = Parameter(name = 'ydo',
+                nature = 'internal',
+                type = 'real',
+                value = '(ymdo*cmath.sqrt(2))/v',
+                texname = '\\text{ydo}')
+
+ye = Parameter(name = 'ye',
+               nature = 'internal',
+               type = 'real',
+               value = '(yme*cmath.sqrt(2))/v',
+               texname = '\\text{ye}')
+
+ym = Parameter(name = 'ym',
+               nature = 'internal',
+               type = 'real',
+               value = '(ymm*cmath.sqrt(2))/v',
+               texname = '\\text{ym}')
+
+ys = Parameter(name = 'ys',
+               nature = 'internal',
+               type = 'real',
+               value = '(yms*cmath.sqrt(2))/v',
+               texname = '\\text{ys}')
+
+yt = Parameter(name = 'yt',
+               nature = 'internal',
+               type = 'real',
+               value = '(ymt*cmath.sqrt(2))/v',
+               texname = '\\text{yt}')
+
+ytau = Parameter(name = 'ytau',
+                 nature = 'internal',
+                 type = 'real',
+                 value = '(ymtau*cmath.sqrt(2))/v',
+                 texname = '\\text{ytau}')
+
+yup = Parameter(name = 'yup',
+                nature = 'internal',
+                type = 'real',
+                value = '(ymup*cmath.sqrt(2))/v',
+                texname = '\\text{yup}')
+
+muH = Parameter(name = 'muH',
+                nature = 'internal',
+                type = 'real',
+                value = 'cmath.sqrt(lam*v**2)',
+                texname = '\\mu ')
+

=== added file 'models/TopEffTh/particles.py'
--- models/TopEffTh/particles.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/particles.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,449 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:06
+
+
+from __future__ import division
+from object_library import all_particles, Particle
+import parameters as Param
+
+ve = Particle(pdg_code = 12,
+              name = 've',
+              antiname = 've~',
+              spin = 2,
+              color = 1,
+              mass = Param.ZERO,
+              width = Param.ZERO,
+              texname = 've',
+              antitexname = 've',
+              charge = 0,
+              LeptonNumber = 1,
+              GhostNumber = 0)
+
+ve__tilde__ = ve.anti()
+
+vm = Particle(pdg_code = 14,
+              name = 'vm',
+              antiname = 'vm~',
+              spin = 2,
+              color = 1,
+              mass = Param.ZERO,
+              width = Param.ZERO,
+              texname = 'vm',
+              antitexname = 'vm',
+              charge = 0,
+              LeptonNumber = 1,
+              GhostNumber = 0)
+
+vm__tilde__ = vm.anti()
+
+vt = Particle(pdg_code = 16,
+              name = 'vt',
+              antiname = 'vt~',
+              spin = 2,
+              color = 1,
+              mass = Param.ZERO,
+              width = Param.ZERO,
+              texname = 'vt',
+              antitexname = 'vt',
+              charge = 0,
+              LeptonNumber = 1,
+              GhostNumber = 0)
+
+vt__tilde__ = vt.anti()
+
+e__minus__ = Particle(pdg_code = 11,
+                      name = 'e-',
+                      antiname = 'e+',
+                      spin = 2,
+                      color = 1,
+                      mass = Param.Me,
+                      width = Param.ZERO,
+                      texname = 'e-',
+                      antitexname = 'e-',
+                      charge = -1,
+                      LeptonNumber = 1,
+                      GhostNumber = 0)
+
+e__plus__ = e__minus__.anti()
+
+m__minus__ = Particle(pdg_code = 13,
+                      name = 'm-',
+                      antiname = 'm+',
+                      spin = 2,
+                      color = 1,
+                      mass = Param.MM,
+                      width = Param.ZERO,
+                      texname = 'm-',
+                      antitexname = 'm-',
+                      charge = -1,
+                      LeptonNumber = 1,
+                      GhostNumber = 0)
+
+m__plus__ = m__minus__.anti()
+
+tt__minus__ = Particle(pdg_code = 15,
+                       name = 'tt-',
+                       antiname = 'tt+',
+                       spin = 2,
+                       color = 1,
+                       mass = Param.MTA,
+                       width = Param.ZERO,
+                       texname = 'tt-',
+                       antitexname = 'tt-',
+                       charge = -1,
+                       LeptonNumber = 1,
+                       GhostNumber = 0)
+
+tt__plus__ = tt__minus__.anti()
+
+u = Particle(pdg_code = 2,
+             name = 'u',
+             antiname = 'u~',
+             spin = 2,
+             color = 3,
+             mass = Param.MU,
+             width = Param.ZERO,
+             texname = 'u',
+             antitexname = 'u',
+             charge = 2/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+u__tilde__ = u.anti()
+
+c = Particle(pdg_code = 4,
+             name = 'c',
+             antiname = 'c~',
+             spin = 2,
+             color = 3,
+             mass = Param.MC,
+             width = Param.ZERO,
+             texname = 'c',
+             antitexname = 'c',
+             charge = 2/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+c__tilde__ = c.anti()
+
+t = Particle(pdg_code = 6,
+             name = 't',
+             antiname = 't~',
+             spin = 2,
+             color = 3,
+             mass = Param.MT,
+             width = Param.WT,
+             texname = 't',
+             antitexname = 't',
+             charge = 2/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+t__tilde__ = t.anti()
+
+d = Particle(pdg_code = 1,
+             name = 'd',
+             antiname = 'd~',
+             spin = 2,
+             color = 3,
+             mass = Param.MD,
+             width = Param.ZERO,
+             texname = 'd',
+             antitexname = 'd',
+             charge = -1/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+d__tilde__ = d.anti()
+
+s = Particle(pdg_code = 3,
+             name = 's',
+             antiname = 's~',
+             spin = 2,
+             color = 3,
+             mass = Param.MS,
+             width = Param.ZERO,
+             texname = 's',
+             antitexname = 's',
+             charge = -1/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+s__tilde__ = s.anti()
+
+b = Particle(pdg_code = 5,
+             name = 'b',
+             antiname = 'b~',
+             spin = 2,
+             color = 3,
+             mass = Param.MB,
+             width = Param.ZERO,
+             texname = 'b',
+             antitexname = 'b',
+             charge = -1/3,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+b__tilde__ = b.anti()
+
+ghA = Particle(pdg_code = 9000001,
+               name = 'ghA',
+               antiname = 'ghA~',
+               spin = -1,
+               color = 1,
+               mass = Param.ZERO,
+               width = Param.ZERO,
+               texname = 'ghA',
+               antitexname = 'ghA',
+               charge = 0,
+               LeptonNumber = 0,
+               GhostNumber = 1)
+
+ghA__tilde__ = ghA.anti()
+
+ghZ = Particle(pdg_code = 9000002,
+               name = 'ghZ',
+               antiname = 'ghZ~',
+               spin = -1,
+               color = 1,
+               mass = Param.MZ,
+               width = Param.ZERO,
+               texname = 'ghZ',
+               antitexname = 'ghZ',
+               charge = 0,
+               LeptonNumber = 0,
+               GhostNumber = 1)
+
+ghZ__tilde__ = ghZ.anti()
+
+ghWp = Particle(pdg_code = 9000003,
+                name = 'ghWp',
+                antiname = 'ghWp~',
+                spin = -1,
+                color = 1,
+                mass = Param.MW,
+                width = Param.ZERO,
+                texname = 'ghWp',
+                antitexname = 'ghWp',
+                charge = 1,
+                LeptonNumber = 0,
+                GhostNumber = 1)
+
+ghWp__tilde__ = ghWp.anti()
+
+ghWm = Particle(pdg_code = 9000004,
+                name = 'ghWm',
+                antiname = 'ghWm~',
+                spin = -1,
+                color = 1,
+                mass = Param.MW,
+                width = Param.ZERO,
+                texname = 'ghWm',
+                antitexname = 'ghWm',
+                charge = -1,
+                LeptonNumber = 0,
+                GhostNumber = 1)
+
+ghWm__tilde__ = ghWm.anti()
+
+ghG = Particle(pdg_code = 9000005,
+               name = 'ghG',
+               antiname = 'ghG~',
+               spin = -1,
+               color = 8,
+               mass = Param.ZERO,
+               width = Param.ZERO,
+               texname = 'ghG',
+               antitexname = 'ghG',
+               charge = 0,
+               LeptonNumber = 0,
+               GhostNumber = 1)
+
+ghG__tilde__ = ghG.anti()
+
+A = Particle(pdg_code = 22,
+             name = 'A',
+             antiname = 'A',
+             spin = 3,
+             color = 1,
+             mass = Param.ZERO,
+             width = Param.ZERO,
+             texname = 'A',
+             antitexname = 'A',
+             charge = 0,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+Z = Particle(pdg_code = 23,
+             name = 'Z',
+             antiname = 'Z',
+             spin = 3,
+             color = 1,
+             mass = Param.MZ,
+             width = Param.WZ,
+             texname = 'Z',
+             antitexname = 'Z',
+             charge = 0,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+W__plus__ = Particle(pdg_code = 24,
+                     name = 'W+',
+                     antiname = 'W-',
+                     spin = 3,
+                     color = 1,
+                     mass = Param.MW,
+                     width = Param.WW,
+                     texname = 'W+',
+                     antitexname = 'W+',
+                     charge = 1,
+                     LeptonNumber = 0,
+                     GhostNumber = 0)
+
+W__minus__ = W__plus__.anti()
+
+G = Particle(pdg_code = 21,
+             name = 'G',
+             antiname = 'G',
+             spin = 3,
+             color = 8,
+             mass = Param.ZERO,
+             width = Param.ZERO,
+             texname = 'G',
+             antitexname = 'G',
+             charge = 0,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+H = Particle(pdg_code = 25,
+             name = 'H',
+             antiname = 'H',
+             spin = 1,
+             color = 1,
+             mass = Param.MH,
+             width = Param.WH,
+             texname = '\\phi',
+             antitexname = '\\phi',
+             charge = 0,
+             LeptonNumber = 0,
+             GhostNumber = 0)
+
+phi0 = Particle(pdg_code = 250,
+                name = 'phi0',
+                antiname = 'phi0',
+                spin = 1,
+                color = 1,
+                mass = Param.MZ,
+                width = Param.ZERO,
+                texname = 'phi0',
+                antitexname = 'phi0',
+                GoldstoneBoson = True,
+                charge = 0,
+                LeptonNumber = 0,
+                GhostNumber = 0)
+
+phi__plus__ = Particle(pdg_code = 251,
+                       name = 'phi+',
+                       antiname = 'phi-',
+                       spin = 1,
+                       color = 1,
+                       mass = Param.MW,
+                       width = Param.ZERO,
+                       texname = '\\phi^+',
+                       antitexname = '\\phi^+',
+                       GoldstoneBoson = True,
+                       charge = 1,
+                       LeptonNumber = 0,
+                       GhostNumber = 0)
+
+phi__minus__ = phi__plus__.anti()
+
+Tri0 = Particle(pdg_code = 9000006,
+                name = 'Tri0',
+                antiname = 'Tri0',
+                spin = 3,
+                color = 1,
+                mass = Param.MTri,
+                width = Param.ZERO,
+                texname = 'Tri0',
+                antitexname = 'Tri0',
+                charge = 0,
+                LeptonNumber = 0,
+                GhostNumber = 0,
+                propagating=False)
+
+Tri = Particle(pdg_code = 9000007,
+               name = 'Tri',
+               antiname = 'Tri~',
+               spin = 3,
+               color = 1,
+               mass = Param.MTri,
+               width = Param.ZERO,
+               texname = 'Tri',
+               antitexname = 'Tri',
+               charge = 1,
+               LeptonNumber = 0,
+               GhostNumber = 0,
+                propagating=False)
+
+Tri__tilde__ = Tri.anti()
+
+Tri80 = Particle(pdg_code = 9000008,
+                 name = 'Tri80',
+                 antiname = 'Tri80',
+                 spin = 3,
+                 color = 8,
+                 mass = Param.MTri8,
+                 width = Param.ZERO,
+                 texname = 'Tri80',
+                 antitexname = 'Tri80',
+                 charge = 0,
+                 LeptonNumber = 0,
+                 GhostNumber = 0,
+                propagating=False)
+
+Tri8 = Particle(pdg_code = 9000009,
+                name = 'Tri8',
+                antiname = 'Tri8~',
+                spin = 3,
+                color = 8,
+                mass = Param.MTri8,
+                width = Param.ZERO,
+                texname = 'Tri8',
+                antitexname = 'Tri8',
+                charge = 1,
+                LeptonNumber = 0,
+                GhostNumber = 0,
+                propagating=False)
+
+Tri8__tilde__ = Tri8.anti()
+
+V8t = Particle(pdg_code = 9000010,
+               name = 'V8t',
+               antiname = 'V8t',
+               spin = 3,
+               color = 8,
+               mass = Param.M8t,
+               width = Param.ZERO,
+               texname = 'V8t',
+               antitexname = 'V8t',
+               charge = 0,
+               LeptonNumber = 0,
+               GhostNumber = 0,
+               propagating=False)
+
+V8Q = Particle(pdg_code = 9000011,
+               name = 'V8Q',
+               antiname = 'V8Q',
+               spin = 3,
+               color = 8,
+               mass = Param.M8Q,
+               width = Param.ZERO,
+               texname = 'V8Q',
+               antitexname = 'V8Q',
+               charge = 0,
+               LeptonNumber = 0,
+               GhostNumber = 0,
+               propagating=False)
+

=== added file 'models/TopEffTh/vertices.py'
--- models/TopEffTh/vertices.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/vertices.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,863 @@
+# This file was automatically created by FeynRules $Revision: 821 $
+# Mathematica version: 7.0 for Microsoft Windows (32-bit) (February 18, 2009)
+# Date: Mon 3 Oct 2011 13:27:06
+
+
+from object_library import all_vertices, Vertex
+import particles as P
+import couplings as C
+import lorentz as L
+
+
+V_1 = Vertex(name = 'V_1',
+             particles = [ P.H, P.H, P.H, P.H ],
+             color = [ '1' ],
+             lorentz = [ L.SSSS1 ],
+             couplings = {(0,0):C.GC_30})
+
+V_2 = Vertex(name = 'V_2',
+             particles = [ P.H, P.H, P.H ],
+             color = [ '1' ],
+             lorentz = [ L.SSS1 ],
+             couplings = {(0,0):C.GC_69})
+
+V_3 = Vertex(name = 'V_3',
+             particles = [ P.G, P.G, P.H, P.H ],
+             color = [ 'Identity(1,2)' ],
+             lorentz = [ L.VVSS2 ],
+             couplings = {(0,0):C.GC_32})
+
+V_4 = Vertex(name = 'V_4',
+             particles = [ P.G, P.G, P.H ],
+             color = [ 'Identity(1,2)' ],
+             lorentz = [ L.VVS2 ],
+             couplings = {(0,0):C.GC_70})
+
+V_5 = Vertex(name = 'V_5',
+             particles = [ P.ghG, P.ghG__tilde__, P.G ],
+             color = [ 'f(3,1,2)' ],
+             lorentz = [ L.UUV1 ],
+             couplings = {(0,0):C.GC_4})
+
+V_6 = Vertex(name = 'V_6',
+             particles = [ P.G, P.G, P.G ],
+             color = [ 'f(1,2,3)' ],
+             lorentz = [ L.VVV1, L.VVV2 ],
+             couplings = {(0,1):C.GC_31,(0,0):C.GC_4})
+
+V_7 = Vertex(name = 'V_7',
+             particles = [ P.G, P.G, P.G, P.G ],
+             color = [ 'f(-1,1,2)*f(3,4,-1)', 'f(-1,1,3)*f(2,4,-1)', 'f(-1,1,4)*f(2,3,-1)' ],
+             lorentz = [ L.VVVV1, L.VVVV2, L.VVVV4, L.VVVV5, L.VVVV7, L.VVVV8 ],
+             couplings = {(0,1):C.GC_36,(1,5):C.GC_36,(2,4):C.GC_36,(1,2):C.GC_6,(0,0):C.GC_6,(2,3):C.GC_6})
+
+V_8 = Vertex(name = 'V_8',
+             particles = [ P.G, P.G, P.G, P.H, P.H ],
+             color = [ 'f(1,2,3)' ],
+             lorentz = [ L.VVVSS1 ],
+             couplings = {(0,0):C.GC_37})
+
+V_9 = Vertex(name = 'V_9',
+             particles = [ P.G, P.G, P.G, P.H ],
+             color = [ 'f(1,2,3)' ],
+             lorentz = [ L.VVVS1 ],
+             couplings = {(0,0):C.GC_74})
+
+V_10 = Vertex(name = 'V_10',
+              particles = [ P.G, P.G, P.G, P.G, P.H, P.H ],
+              color = [ 'f(-1,1,2)*f(3,4,-1)', 'f(-1,1,3)*f(2,4,-1)', 'f(-1,1,4)*f(2,3,-1)' ],
+              lorentz = [ L.VVVVSS1, L.VVVVSS2, L.VVVVSS3 ],
+              couplings = {(1,1):C.GC_41,(0,0):C.GC_41,(2,2):C.GC_41})
+
+V_11 = Vertex(name = 'V_11',
+              particles = [ P.G, P.G, P.G, P.G, P.H ],
+              color = [ 'f(-1,1,2)*f(3,4,-1)', 'f(-1,1,3)*f(2,4,-1)', 'f(-1,1,4)*f(2,3,-1)' ],
+              lorentz = [ L.VVVVS1, L.VVVVS2, L.VVVVS3 ],
+              couplings = {(1,1):C.GC_76,(0,0):C.GC_76,(2,2):C.GC_76})
+
+V_12 = Vertex(name = 'V_12',
+              particles = [ P.G, P.G, P.G, P.G, P.G ],
+              color = [ 'f(-1,1,-2)*f(2,3,-1)*f(4,5,-2)', 'f(-1,1,-2)*f(2,3,-2)*f(4,5,-1)', 'f(-1,1,-2)*f(2,4,-1)*f(3,5,-2)', 'f(-1,1,-2)*f(2,4,-2)*f(3,5,-1)', 'f(-1,1,-2)*f(2,5,-1)*f(3,4,-2)', 'f(-1,1,-2)*f(2,5,-2)*f(3,4,-1)', 'f(-1,1,2)*f(3,-2,-1)*f(4,5,-2)', 'f(-1,1,2)*f(3,4,-2)*f(5,-2,-1)', 'f(-1,1,2)*f(3,5,-2)*f(4,-2,-1)', 'f(-1,1,3)*f(2,-2,-1)*f(4,5,-2)', 'f(-1,1,3)*f(2,4,-2)*f(5,-2,-1)', 'f(-1,1,3)*f(2,5,-2)*f(4,-2,-1)', 'f(-1,1,4)*f(2,-2,-1)*f(3,5,-2)', 'f(-1,1,4)*f(2,3,-2)*f(5,-2,-1)', 'f(-1,1,4)*f(2,5,-2)*f(3,-2,-1)', 'f(-1,1,5)*f(2,-2,-1)*f(3,4,-2)', 'f(-1,1,5)*f(2,3,-2)*f(4,-2,-1)', 'f(-1,1,5)*f(2,4,-2)*f(3,-2,-1)', 'f(-1,2,-2)*f(1,3,-2)*f(4,5,-1)', 'f(-1,2,-2)*f(1,4,-2)*f(3,5,-1)', 'f(-1,2,-2)*f(1,5,-2)*f(3,4,-1)', 'f(-1,2,3)*f(1,4,-2)*f(5,-2,-1)', 'f(-1,2,3)*f(1,5,-2)*f(4,-2,-1)', 'f(-1,2,4)*f(1,3,-2)*f(5,-2,-1)', 'f(-1,2,4)*f(1,5,-2)*f(3,-2,-1)', 'f(-1,2,5)*f(1,3,-2)*f(4,-2,-1)', 'f(-1,2,5)*f(1,4,-2)*f(3,-2,-1)', 'f(-1,3,-2)*f(1,2,-2)*f(4,5,-1)', 'f(-1,3,4)*f(1,2,-2)*f(5,-2,-1)', 'f(-1,3,5)*f(1,2,-2)*f(4,-2,-1)' ],
+              lorentz = [ L.VVVVV1, L.VVVVV10, L.VVVVV11, L.VVVVV12, L.VVVVV13, L.VVVVV14, L.VVVVV15, L.VVVVV16, L.VVVVV17, L.VVVVV18, L.VVVVV19, L.VVVVV2, L.VVVVV20, L.VVVVV21, L.VVVVV22, L.VVVVV23, L.VVVVV24, L.VVVVV25, L.VVVVV26, L.VVVVV27, L.VVVVV28, L.VVVVV29, L.VVVVV3, L.VVVVV30, L.VVVVV4, L.VVVVV5, L.VVVVV6, L.VVVVV7, L.VVVVV8, L.VVVVV9 ],
+              couplings = {(2,4):C.GC_39,(4,9):C.GC_40,(5,9):C.GC_39,(3,4):C.GC_40,(15,16):C.GC_39,(9,25):C.GC_39,(20,16):C.GC_40,(18,25):C.GC_40,(17,17):C.GC_39,(6,0):C.GC_39,(24,17):C.GC_40,(27,0):C.GC_40,(25,27):C.GC_39,(29,11):C.GC_39,(8,11):C.GC_40,(11,27):C.GC_40,(2,8):C.GC_40,(4,5):C.GC_39,(5,5):C.GC_40,(3,8):C.GC_39,(12,15):C.GC_39,(19,15):C.GC_40,(14,18):C.GC_39,(26,18):C.GC_40,(23,28):C.GC_39,(28,22):C.GC_39,(7,22):C.GC_40,(10,28):C.GC_40,(0,24):C.GC_39,(1,24):C.GC_40,(12,6):C.GC_39,(19,6):C.GC_40,(26,10):C.GC_39,(14,10):C.GC_40,(16,2):C.GC_39,(22,2):C.GC_40,(0,26):C.GC_40,(1,26):C.GC_39,(9,29):C.GC_39,(18,29):C.GC_40,(11,19):C.GC_39,(25,19):C.GC_40,(21,14):C.GC_39,(13,14):C.GC_40,(15,7):C.GC_39,(20,7):C.GC_40,(24,12):C.GC_39,(17,12):C.GC_40,(13,3):C.GC_39,(21,3):C.GC_40,(22,13):C.GC_39,(16,13):C.GC_40,(10,21):C.GC_39,(23,21):C.GC_40,(6,1):C.GC_39,(27,1):C.GC_40,(8,20):C.GC_39,(29,20):C.GC_40,(7,23):C.GC_39,(28,23):C.GC_40})
+
+V_13 = Vertex(name = 'V_13',
+              particles = [ P.G, P.G, P.G, P.G, P.G, P.G ],
+              color = [ 'f(-2,-3,-1)*f(-1,1,2)*f(3,4,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,1,2)*f(3,4,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,1,2)*f(3,5,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,1,2)*f(3,5,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,1,2)*f(3,6,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,1,2)*f(3,6,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,4,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,4,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,5,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,5,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,6,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,1,3)*f(2,6,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,3,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,3,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,5,-2)*f(3,6,-3)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,5,-3)*f(3,6,-2)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,6,-2)*f(3,5,-3)', 'f(-2,-3,-1)*f(-1,1,4)*f(2,6,-3)*f(3,5,-2)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,3,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,3,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,4,-2)*f(3,6,-3)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,4,-3)*f(3,6,-2)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,6,-2)*f(3,4,-3)', 'f(-2,-3,-1)*f(-1,1,5)*f(2,6,-3)*f(3,4,-2)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,3,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,3,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,4,-2)*f(3,5,-3)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,4,-3)*f(3,5,-2)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,5,-2)*f(3,4,-3)', 'f(-2,-3,-1)*f(-1,1,6)*f(2,5,-3)*f(3,4,-2)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,4,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,4,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,5,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,5,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,6,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,2,3)*f(1,6,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,3,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,3,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,5,-2)*f(3,6,-3)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,5,-3)*f(3,6,-2)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,6,-2)*f(3,5,-3)', 'f(-2,-3,-1)*f(-1,2,4)*f(1,6,-3)*f(3,5,-2)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,3,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,3,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,4,-2)*f(3,6,-3)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,4,-3)*f(3,6,-2)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,6,-2)*f(3,4,-3)', 'f(-2,-3,-1)*f(-1,2,5)*f(1,6,-3)*f(3,4,-2)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,3,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,3,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,4,-2)*f(3,5,-3)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,4,-3)*f(3,5,-2)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,5,-2)*f(3,4,-3)', 'f(-2,-3,-1)*f(-1,2,6)*f(1,5,-3)*f(3,4,-2)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,2,-2)*f(5,6,-3)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,2,-3)*f(5,6,-2)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,5,-2)*f(2,6,-3)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,5,-3)*f(2,6,-2)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,6,-2)*f(2,5,-3)', 'f(-2,-3,-1)*f(-1,3,4)*f(1,6,-3)*f(2,5,-2)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,2,-2)*f(4,6,-3)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,2,-3)*f(4,6,-2)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,4,-2)*f(2,6,-3)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,4,-3)*f(2,6,-2)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,6,-2)*f(2,4,-3)', 'f(-2,-3,-1)*f(-1,3,5)*f(1,6,-3)*f(2,4,-2)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,2,-2)*f(4,5,-3)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,2,-3)*f(4,5,-2)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,4,-2)*f(2,5,-3)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,4,-3)*f(2,5,-2)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,5,-2)*f(2,4,-3)', 'f(-2,-3,-1)*f(-1,3,6)*f(1,5,-3)*f(2,4,-2)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,2,-2)*f(3,6,-3)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,2,-3)*f(3,6,-2)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,3,-2)*f(2,6,-3)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,3,-3)*f(2,6,-2)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,6,-2)*f(2,3,-3)', 'f(-2,-3,-1)*f(-1,4,5)*f(1,6,-3)*f(2,3,-2)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,2,-2)*f(3,5,-3)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,2,-3)*f(3,5,-2)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,3,-2)*f(2,5,-3)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,3,-3)*f(2,5,-2)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,5,-2)*f(2,3,-3)', 'f(-2,-3,-1)*f(-1,4,6)*f(1,5,-3)*f(2,3,-2)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,2,-2)*f(3,4,-3)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,2,-3)*f(3,4,-2)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,3,-2)*f(2,4,-3)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,3,-3)*f(2,4,-2)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,4,-2)*f(2,3,-3)', 'f(-2,-3,-1)*f(-1,5,6)*f(1,4,-3)*f(2,3,-2)' ],
+              lorentz = [ L.VVVVVV1, L.VVVVVV10, L.VVVVVV11, L.VVVVVV12, L.VVVVVV13, L.VVVVVV14, L.VVVVVV15, L.VVVVVV2, L.VVVVVV3, L.VVVVVV4, L.VVVVVV5, L.VVVVVV6, L.VVVVVV7, L.VVVVVV8, L.VVVVVV9 ],
+              couplings = {(5,11):C.GC_42,(4,11):C.GC_43,(3,9):C.GC_43,(2,9):C.GC_42,(11,1):C.GC_43,(10,1):C.GC_42,(7,7):C.GC_43,(6,7):C.GC_42,(17,4):C.GC_43,(16,4):C.GC_42,(13,8):C.GC_43,(12,8):C.GC_42,(21,2):C.GC_42,(20,2):C.GC_43,(19,3):C.GC_42,(18,3):C.GC_43,(33,3):C.GC_43,(32,3):C.GC_42,(31,8):C.GC_42,(30,8):C.GC_43,(39,2):C.GC_43,(38,2):C.GC_42,(37,7):C.GC_42,(36,7):C.GC_43,(51,4):C.GC_42,(50,4):C.GC_43,(49,1):C.GC_42,(48,1):C.GC_43,(63,4):C.GC_43,(62,4):C.GC_42,(61,9):C.GC_42,(60,9):C.GC_43,(71,2):C.GC_42,(70,2):C.GC_43,(67,11):C.GC_43,(66,11):C.GC_42,(75,1):C.GC_43,(74,1):C.GC_42,(73,11):C.GC_42,(72,11):C.GC_43,(83,3):C.GC_42,(82,3):C.GC_43,(79,9):C.GC_43,(78,9):C.GC_42,(89,8):C.GC_43,(88,8):C.GC_42,(87,7):C.GC_43,(86,7):C.GC_42,(9,13):C.GC_43,(8,13):C.GC_42,(15,5):C.GC_43,(14,5):C.GC_42,(27,14):C.GC_42,(26,14):C.GC_43,(25,6):C.GC_42,(24,6):C.GC_43,(35,6):C.GC_43,(34,6):C.GC_42,(41,14):C.GC_43,(40,14):C.GC_42,(45,5):C.GC_42,(44,5):C.GC_43,(43,13):C.GC_42,(42,13):C.GC_43,(65,14):C.GC_42,(64,14):C.GC_43,(69,5):C.GC_43,(68,5):C.GC_42,(77,6):C.GC_42,(76,6):C.GC_43,(81,13):C.GC_43,(80,13):C.GC_42,(1,0):C.GC_43,(0,0):C.GC_42,(23,10):C.GC_43,(22,10):C.GC_42,(53,10):C.GC_42,(52,10):C.GC_43,(57,10):C.GC_43,(56,10):C.GC_42,(55,0):C.GC_42,(54,0):C.GC_43,(85,0):C.GC_43,(84,0):C.GC_42,(29,12):C.GC_42,(28,12):C.GC_43,(47,12):C.GC_43,(46,12):C.GC_42,(59,12):C.GC_42,(58,12):C.GC_43})
+
+V_14 = Vertex(name = 'V_14',
+              particles = [ P.A, P.W__minus__, P.W__plus__ ],
+              color = [ '1' ],
+              lorentz = [ L.VVV1 ],
+              couplings = {(0,0):C.GC_62})
+
+V_15 = Vertex(name = 'V_15',
+              particles = [ P.W__minus__, P.W__plus__, P.H, P.H ],
+              color = [ '1' ],
+              lorentz = [ L.VVSS1 ],
+              couplings = {(0,0):C.GC_46})
+
+V_16 = Vertex(name = 'V_16',
+              particles = [ P.W__minus__, P.W__plus__, P.H ],
+              color = [ '1' ],
+              lorentz = [ L.VVS1 ],
+              couplings = {(0,0):C.GC_79})
+
+V_17 = Vertex(name = 'V_17',
+              particles = [ P.A, P.A, P.W__minus__, P.W__plus__ ],
+              color = [ '1' ],
+              lorentz = [ L.VVVV3 ],
+              couplings = {(0,0):C.GC_66})
+
+V_18 = Vertex(name = 'V_18',
+              particles = [ P.W__minus__, P.W__plus__, P.Z ],
+              color = [ '1' ],
+              lorentz = [ L.VVV1 ],
+              couplings = {(0,0):C.GC_27})
+
+V_19 = Vertex(name = 'V_19',
+              particles = [ P.W__minus__, P.W__minus__, P.W__plus__, P.W__plus__ ],
+              color = [ '1' ],
+              lorentz = [ L.VVVV3 ],
+              couplings = {(0,0):C.GC_28})
+
+V_20 = Vertex(name = 'V_20',
+              particles = [ P.A, P.W__minus__, P.W__plus__, P.Z ],
+              color = [ '1' ],
+              lorentz = [ L.VVVV6 ],
+              couplings = {(0,0):C.GC_63})
+
+V_21 = Vertex(name = 'V_21',
+              particles = [ P.Z, P.Z, P.H, P.H ],
+              color = [ '1' ],
+              lorentz = [ L.VVSS1 ],
+              couplings = {(0,0):C.GC_68})
+
+V_22 = Vertex(name = 'V_22',
+              particles = [ P.Z, P.Z, P.H ],
+              color = [ '1' ],
+              lorentz = [ L.VVS1 ],
+              couplings = {(0,0):C.GC_84})
+
+V_23 = Vertex(name = 'V_23',
+              particles = [ P.W__minus__, P.W__plus__, P.Z, P.Z ],
+              color = [ '1' ],
+              lorentz = [ L.VVVV3 ],
+              couplings = {(0,0):C.GC_29})
+
+V_24 = Vertex(name = 'V_24',
+              particles = [ P.t__tilde__, P.b, P.W__plus__, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVSS1 ],
+              couplings = {(0,0):C.GC_94})
+
+V_25 = Vertex(name = 'V_25',
+              particles = [ P.t__tilde__, P.b, P.W__plus__, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVS1, L.FFVS3 ],
+              couplings = {(0,0):C.GC_95,(0,1):C.GC_34})
+
+V_26 = Vertex(name = 'V_26',
+              particles = [ P.t__tilde__, P.b, P.W__plus__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2, L.FFV8 ],
+              couplings = {(0,0):[ C.GC_96, C.GC_115 ],(0,1):C.GC_72})
+
+V_27 = Vertex(name = 'V_27',
+              particles = [ P.b__tilde__, P.b, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFSS1 ],
+              couplings = {(0,0):C.GC_98})
+
+V_28 = Vertex(name = 'V_28',
+              particles = [ P.b__tilde__, P.b, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFS1, L.FFS2 ],
+              couplings = {(0,0):C.GC_85,(0,1):C.GC_102})
+
+V_29 = Vertex(name = 'V_29',
+              particles = [ P.b__tilde__, P.b, P.Z, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVSS1 ],
+              couplings = {(0,0):C.GC_99})
+
+V_30 = Vertex(name = 'V_30',
+              particles = [ P.b__tilde__, P.b, P.Z, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVS1 ],
+              couplings = {(0,0):C.GC_103})
+
+V_31 = Vertex(name = 'V_31',
+              particles = [ P.b__tilde__, P.b, P.Z ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2, L.FFV5 ],
+              couplings = {(0,0):[ C.GC_57, C.GC_105 ],(0,1):C.GC_60})
+
+V_32 = Vertex(name = 'V_32',
+              particles = [ P.t__tilde__, P.t, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFSS1 ],
+              couplings = {(0,0):C.GC_97})
+
+V_33 = Vertex(name = 'V_33',
+              particles = [ P.t__tilde__, P.t, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFS1, L.FFS2 ],
+              couplings = {(0,0):C.GC_91,(0,1):C.GC_101})
+
+V_34 = Vertex(name = 'V_34',
+              particles = [ P.t__tilde__, P.t, P.Z, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVSS1 ],
+              couplings = {(0,0):C.GC_100})
+
+V_35 = Vertex(name = 'V_35',
+              particles = [ P.t__tilde__, P.t, P.Z, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVS1, L.FFVS2, L.FFVS3 ],
+              couplings = {(0,0):C.GC_104,(0,1):C.GC_121,(0,2):C.GC_35})
+
+V_36 = Vertex(name = 'V_36',
+              particles = [ P.t__tilde__, P.t, P.Z ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2, L.FFV3, L.FFV7, L.FFV8 ],
+              couplings = {(0,0):[ C.GC_58, C.GC_106 ],(0,2):C.GC_60,(0,1):C.GC_127,(0,3):C.GC_73})
+
+V_37 = Vertex(name = 'V_37',
+              particles = [ P.b__tilde__, P.t, P.W__minus__, P.H, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVSS1 ],
+              couplings = {(0,0):C.GC_59})
+
+V_38 = Vertex(name = 'V_38',
+              particles = [ P.b__tilde__, P.t, P.W__minus__, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVS1, L.FFVS2 ],
+              couplings = {(0,0):C.GC_80,(0,1):C.GC_120})
+
+V_39 = Vertex(name = 'V_39',
+              particles = [ P.b__tilde__, P.t, P.W__minus__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2, L.FFV3 ],
+              couplings = {(0,0):[ C.GC_56, C.GC_83 ],(0,1):C.GC_126})
+
+V_40 = Vertex(name = 'V_40',
+              particles = [ P.s__tilde__, P.c, P.Tri__tilde__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_26})
+
+V_41 = Vertex(name = 'V_41',
+              particles = [ P.d__tilde__, P.u, P.Tri__tilde__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_26})
+
+V_42 = Vertex(name = 'V_42',
+              particles = [ P.c__tilde__, P.c, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_25})
+
+V_43 = Vertex(name = 'V_43',
+              particles = [ P.d__tilde__, P.d, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_24})
+
+V_44 = Vertex(name = 'V_44',
+              particles = [ P.s__tilde__, P.s, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_24})
+
+V_45 = Vertex(name = 'V_45',
+              particles = [ P.u__tilde__, P.u, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_25})
+
+V_46 = Vertex(name = 'V_46',
+              particles = [ P.u__tilde__, P.d, P.Tri ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_26})
+
+V_47 = Vertex(name = 'V_47',
+              particles = [ P.c__tilde__, P.s, P.Tri ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_26})
+
+V_48 = Vertex(name = 'V_48',
+              particles = [ P.b__tilde__, P.t, P.Tri__tilde__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_17})
+
+V_49 = Vertex(name = 'V_49',
+              particles = [ P.b__tilde__, P.b, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_15})
+
+V_50 = Vertex(name = 'V_50',
+              particles = [ P.t__tilde__, P.t, P.Tri0 ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_16})
+
+V_51 = Vertex(name = 'V_51',
+              particles = [ P.t__tilde__, P.b, P.Tri ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_17})
+
+V_52 = Vertex(name = 'V_52',
+              particles = [ P.c__tilde__, P.c, P.V8t ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_13,(0,1):C.GC_14})
+
+V_53 = Vertex(name = 'V_53',
+              particles = [ P.d__tilde__, P.d, P.V8t ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_13,(0,1):C.GC_12})
+
+V_54 = Vertex(name = 'V_54',
+              particles = [ P.s__tilde__, P.s, P.V8t ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_13,(0,1):C.GC_12})
+
+V_55 = Vertex(name = 'V_55',
+              particles = [ P.u__tilde__, P.u, P.V8t ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_13,(0,1):C.GC_14})
+
+V_56 = Vertex(name = 'V_56',
+              particles = [ P.t__tilde__, P.t, P.V8t ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV4 ],
+              couplings = {(0,0):C.GC_11})
+
+V_57 = Vertex(name = 'V_57',
+              particles = [ P.b__tilde__, P.b, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_7})
+
+V_58 = Vertex(name = 'V_58',
+              particles = [ P.c__tilde__, P.c, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_9,(0,1):C.GC_10})
+
+V_59 = Vertex(name = 'V_59',
+              particles = [ P.d__tilde__, P.d, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_9,(0,1):C.GC_8})
+
+V_60 = Vertex(name = 'V_60',
+              particles = [ P.s__tilde__, P.s, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_9,(0,1):C.GC_8})
+
+V_61 = Vertex(name = 'V_61',
+              particles = [ P.t__tilde__, P.t, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_7})
+
+V_62 = Vertex(name = 'V_62',
+              particles = [ P.u__tilde__, P.u, P.V8Q ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2, L.FFV4 ],
+              couplings = {(0,0):C.GC_9,(0,1):C.GC_10})
+
+V_63 = Vertex(name = 'V_63',
+              particles = [ P.s__tilde__, P.c, P.Tri8__tilde__ ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_23})
+
+V_64 = Vertex(name = 'V_64',
+              particles = [ P.d__tilde__, P.u, P.Tri8__tilde__ ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_23})
+
+V_65 = Vertex(name = 'V_65',
+              particles = [ P.c__tilde__, P.c, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_22})
+
+V_66 = Vertex(name = 'V_66',
+              particles = [ P.d__tilde__, P.d, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_21})
+
+V_67 = Vertex(name = 'V_67',
+              particles = [ P.s__tilde__, P.s, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_21})
+
+V_68 = Vertex(name = 'V_68',
+              particles = [ P.u__tilde__, P.u, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_22})
+
+V_69 = Vertex(name = 'V_69',
+              particles = [ P.u__tilde__, P.d, P.Tri8 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_23})
+
+V_70 = Vertex(name = 'V_70',
+              particles = [ P.c__tilde__, P.s, P.Tri8 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_23})
+
+V_71 = Vertex(name = 'V_71',
+              particles = [ P.b__tilde__, P.t, P.Tri8__tilde__ ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_20})
+
+V_72 = Vertex(name = 'V_72',
+              particles = [ P.b__tilde__, P.b, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_18})
+
+V_73 = Vertex(name = 'V_73',
+              particles = [ P.t__tilde__, P.t, P.Tri80 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_19})
+
+V_74 = Vertex(name = 'V_74',
+              particles = [ P.t__tilde__, P.b, P.Tri8 ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV2 ],
+              couplings = {(0,0):C.GC_20})
+
+V_75 = Vertex(name = 'V_75',
+              particles = [ P.t__tilde__, P.t, P.A, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVS2, L.FFVS3 ],
+              couplings = {(0,0):C.GC_124,(0,1):C.GC_64})
+
+V_76 = Vertex(name = 'V_76',
+              particles = [ P.t__tilde__, P.t, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1, L.FFV3, L.FFV8 ],
+              couplings = {(0,0):C.GC_2,(0,1):C.GC_130,(0,2):C.GC_81})
+
+V_77 = Vertex(name = 'V_77',
+              particles = [ P.t__tilde__, P.b, P.A, P.W__plus__, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVVS2 ],
+              couplings = {(0,0):C.GC_65})
+
+V_78 = Vertex(name = 'V_78',
+              particles = [ P.t__tilde__, P.b, P.A, P.W__plus__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVV2 ],
+              couplings = {(0,0):C.GC_82})
+
+V_79 = Vertex(name = 'V_79',
+              particles = [ P.t__tilde__, P.t, P.W__minus__, P.W__plus__, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVVS1, L.FFVVS2 ],
+              couplings = {(0,0):C.GC_122,(0,1):C.GC_44})
+
+V_80 = Vertex(name = 'V_80',
+              particles = [ P.t__tilde__, P.t, P.W__minus__, P.W__plus__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVV1, L.FFVV2 ],
+              couplings = {(0,0):C.GC_128,(0,1):C.GC_77})
+
+V_81 = Vertex(name = 'V_81',
+              particles = [ P.t__tilde__, P.b, P.W__plus__, P.Z, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVVS2 ],
+              couplings = {(0,0):C.GC_45})
+
+V_82 = Vertex(name = 'V_82',
+              particles = [ P.t__tilde__, P.b, P.W__plus__, P.Z ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVV2 ],
+              couplings = {(0,0):C.GC_78})
+
+V_83 = Vertex(name = 'V_83',
+              particles = [ P.t__tilde__, P.t, P.G, P.H ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFVS2, L.FFVS3 ],
+              couplings = {(0,0):C.GC_116,(0,1):C.GC_33})
+
+V_84 = Vertex(name = 'V_84',
+              particles = [ P.t__tilde__, P.t, P.G ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV1, L.FFV3, L.FFV8 ],
+              couplings = {(0,0):C.GC_5,(0,1):C.GC_118,(0,2):C.GC_71})
+
+V_85 = Vertex(name = 'V_85',
+              particles = [ P.t__tilde__, P.t, P.G, P.G, P.H ],
+              color = [ 'f(3,4,-1)*T(-1,2,1)' ],
+              lorentz = [ L.FFVVS1, L.FFVVS2 ],
+              couplings = {(0,0):C.GC_117,(0,1):C.GC_38})
+
+V_86 = Vertex(name = 'V_86',
+              particles = [ P.t__tilde__, P.t, P.G, P.G ],
+              color = [ 'f(3,4,-1)*T(-1,2,1)' ],
+              lorentz = [ L.FFVV1, L.FFVV2 ],
+              couplings = {(0,0):C.GC_119,(0,1):C.GC_75})
+
+V_87 = Vertex(name = 'V_87',
+              particles = [ P.b__tilde__, P.t, P.A, P.W__minus__, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVVS1 ],
+              couplings = {(0,0):C.GC_125})
+
+V_88 = Vertex(name = 'V_88',
+              particles = [ P.b__tilde__, P.t, P.A, P.W__minus__ ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVV1 ],
+              couplings = {(0,0):C.GC_131})
+
+V_89 = Vertex(name = 'V_89',
+              particles = [ P.b__tilde__, P.t, P.W__minus__, P.Z, P.H ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVVS1 ],
+              couplings = {(0,0):C.GC_123})
+
+V_90 = Vertex(name = 'V_90',
+              particles = [ P.b__tilde__, P.t, P.W__minus__, P.Z ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFVV1 ],
+              couplings = {(0,0):C.GC_129})
+
+V_91 = Vertex(name = 'V_91',
+              particles = [ P.d__tilde__, P.d, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_1})
+
+V_92 = Vertex(name = 'V_92',
+              particles = [ P.s__tilde__, P.s, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_1})
+
+V_93 = Vertex(name = 'V_93',
+              particles = [ P.b__tilde__, P.b, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_1})
+
+V_94 = Vertex(name = 'V_94',
+              particles = [ P.e__plus__, P.e__minus__, P.A ],
+              color = [ '1' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_3})
+
+V_95 = Vertex(name = 'V_95',
+              particles = [ P.m__plus__, P.m__minus__, P.A ],
+              color = [ '1' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_3})
+
+V_96 = Vertex(name = 'V_96',
+              particles = [ P.tt__plus__, P.tt__minus__, P.A ],
+              color = [ '1' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_3})
+
+V_97 = Vertex(name = 'V_97',
+              particles = [ P.u__tilde__, P.u, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_2})
+
+V_98 = Vertex(name = 'V_98',
+              particles = [ P.c__tilde__, P.c, P.A ],
+              color = [ 'Identity(1,2)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_2})
+
+V_99 = Vertex(name = 'V_99',
+              particles = [ P.d__tilde__, P.d, P.G ],
+              color = [ 'T(3,2,1)' ],
+              lorentz = [ L.FFV1 ],
+              couplings = {(0,0):C.GC_5})
+
+V_100 = Vertex(name = 'V_100',
+               particles = [ P.s__tilde__, P.s, P.G ],
+               color = [ 'T(3,2,1)' ],
+               lorentz = [ L.FFV1 ],
+               couplings = {(0,0):C.GC_5})
+
+V_101 = Vertex(name = 'V_101',
+               particles = [ P.b__tilde__, P.b, P.G ],
+               color = [ 'T(3,2,1)' ],
+               lorentz = [ L.FFV1 ],
+               couplings = {(0,0):C.GC_5})
+
+V_102 = Vertex(name = 'V_102',
+               particles = [ P.d__tilde__, P.d, P.H ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_87})
+
+V_103 = Vertex(name = 'V_103',
+               particles = [ P.s__tilde__, P.s, P.H ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_90})
+
+V_104 = Vertex(name = 'V_104',
+               particles = [ P.d__tilde__, P.d, P.Z ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2, L.FFV5 ],
+               couplings = {(0,0):C.GC_57,(0,1):C.GC_60})
+
+V_105 = Vertex(name = 'V_105',
+               particles = [ P.s__tilde__, P.s, P.Z ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2, L.FFV5 ],
+               couplings = {(0,0):C.GC_57,(0,1):C.GC_60})
+
+V_106 = Vertex(name = 'V_106',
+               particles = [ P.d__tilde__, P.u, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_48})
+
+V_107 = Vertex(name = 'V_107',
+               particles = [ P.d__tilde__, P.c, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_51})
+
+V_108 = Vertex(name = 'V_108',
+               particles = [ P.d__tilde__, P.t, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_54})
+
+V_109 = Vertex(name = 'V_109',
+               particles = [ P.s__tilde__, P.u, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_49})
+
+V_110 = Vertex(name = 'V_110',
+               particles = [ P.s__tilde__, P.c, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_52})
+
+V_111 = Vertex(name = 'V_111',
+               particles = [ P.s__tilde__, P.t, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_55})
+
+V_112 = Vertex(name = 'V_112',
+               particles = [ P.b__tilde__, P.u, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_50})
+
+V_113 = Vertex(name = 'V_113',
+               particles = [ P.b__tilde__, P.c, P.W__minus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_53})
+
+V_114 = Vertex(name = 'V_114',
+               particles = [ P.u__tilde__, P.d, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_107})
+
+V_115 = Vertex(name = 'V_115',
+               particles = [ P.c__tilde__, P.d, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_110})
+
+V_116 = Vertex(name = 'V_116',
+               particles = [ P.t__tilde__, P.d, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_113})
+
+V_117 = Vertex(name = 'V_117',
+               particles = [ P.u__tilde__, P.s, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_108})
+
+V_118 = Vertex(name = 'V_118',
+               particles = [ P.c__tilde__, P.s, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_111})
+
+V_119 = Vertex(name = 'V_119',
+               particles = [ P.t__tilde__, P.s, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_114})
+
+V_120 = Vertex(name = 'V_120',
+               particles = [ P.u__tilde__, P.b, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_109})
+
+V_121 = Vertex(name = 'V_121',
+               particles = [ P.c__tilde__, P.b, P.W__plus__ ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_112})
+
+V_122 = Vertex(name = 'V_122',
+               particles = [ P.u__tilde__, P.u, P.G ],
+               color = [ 'T(3,2,1)' ],
+               lorentz = [ L.FFV1 ],
+               couplings = {(0,0):C.GC_5})
+
+V_123 = Vertex(name = 'V_123',
+               particles = [ P.c__tilde__, P.c, P.G ],
+               color = [ 'T(3,2,1)' ],
+               lorentz = [ L.FFV1 ],
+               couplings = {(0,0):C.GC_5})
+
+V_124 = Vertex(name = 'V_124',
+               particles = [ P.e__plus__, P.e__minus__, P.H ],
+               color = [ '1' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_88})
+
+V_125 = Vertex(name = 'V_125',
+               particles = [ P.m__plus__, P.m__minus__, P.H ],
+               color = [ '1' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_89})
+
+V_126 = Vertex(name = 'V_126',
+               particles = [ P.tt__plus__, P.tt__minus__, P.H ],
+               color = [ '1' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_92})
+
+V_127 = Vertex(name = 'V_127',
+               particles = [ P.u__tilde__, P.u, P.H ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_93})
+
+V_128 = Vertex(name = 'V_128',
+               particles = [ P.c__tilde__, P.c, P.H ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFS1 ],
+               couplings = {(0,0):C.GC_86})
+
+V_129 = Vertex(name = 'V_129',
+               particles = [ P.e__plus__, P.e__minus__, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2, L.FFV6 ],
+               couplings = {(0,0):C.GC_57,(0,1):C.GC_61})
+
+V_130 = Vertex(name = 'V_130',
+               particles = [ P.m__plus__, P.m__minus__, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2, L.FFV6 ],
+               couplings = {(0,0):C.GC_57,(0,1):C.GC_61})
+
+V_131 = Vertex(name = 'V_131',
+               particles = [ P.tt__plus__, P.tt__minus__, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2, L.FFV6 ],
+               couplings = {(0,0):C.GC_57,(0,1):C.GC_61})
+
+V_132 = Vertex(name = 'V_132',
+               particles = [ P.e__plus__, P.ve, P.W__minus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_133 = Vertex(name = 'V_133',
+               particles = [ P.m__plus__, P.vm, P.W__minus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_134 = Vertex(name = 'V_134',
+               particles = [ P.tt__plus__, P.vt, P.W__minus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_135 = Vertex(name = 'V_135',
+               particles = [ P.ve__tilde__, P.e__minus__, P.W__plus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_136 = Vertex(name = 'V_136',
+               particles = [ P.vm__tilde__, P.m__minus__, P.W__plus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_137 = Vertex(name = 'V_137',
+               particles = [ P.vt__tilde__, P.tt__minus__, P.W__plus__ ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_47})
+
+V_138 = Vertex(name = 'V_138',
+               particles = [ P.u__tilde__, P.u, P.Z ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2, L.FFV7 ],
+               couplings = {(0,0):C.GC_58,(0,1):C.GC_60})
+
+V_139 = Vertex(name = 'V_139',
+               particles = [ P.c__tilde__, P.c, P.Z ],
+               color = [ 'Identity(1,2)' ],
+               lorentz = [ L.FFV2, L.FFV7 ],
+               couplings = {(0,0):C.GC_58,(0,1):C.GC_60})
+
+V_140 = Vertex(name = 'V_140',
+               particles = [ P.ve__tilde__, P.ve, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_67})
+
+V_141 = Vertex(name = 'V_141',
+               particles = [ P.vm__tilde__, P.vm, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_67})
+
+V_142 = Vertex(name = 'V_142',
+               particles = [ P.vt__tilde__, P.vt, P.Z ],
+               color = [ '1' ],
+               lorentz = [ L.FFV2 ],
+               couplings = {(0,0):C.GC_67})
+

=== added file 'models/TopEffTh/write_param_card.py'
--- models/TopEffTh/write_param_card.py	1970-01-01 00:00:00 +0000
+++ models/TopEffTh/write_param_card.py	2011-10-11 00:36:19 +0000
@@ -0,0 +1,181 @@
+
+__date__ = "3 june 2010"
+__author__ = '[e-mail]'
+
+from function_library import *
+
+class ParamCardWriter(object):
+    
+    header = \
+    """######################################################################\n""" + \
+    """## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO  #####################\n""" + \
+    """######################################################################\n"""   
+    
+    def __init__(self, filename, list_of_parameters=None, generic=False):
+        """write a valid param_card.dat"""
+        
+        if not list_of_parameters:
+            from parameters import all_parameters
+            list_of_parameters = [param for param in all_parameters if \
+                                                       param.nature=='external']
+        
+        self.generic_output = generic
+        if generic:
+            self.define_not_dep_param(list_of_parameters)
+
+        
+        self.fsock = open(filename, 'w')
+        self.fsock.write(self.header)
+        
+        self.write_card(list_of_parameters)
+    
+    def define_not_dep_param(self, list_of_parameters):
+        """define self.dep_mass and self.dep_width in case that they are 
+        requested in the param_card.dat"""
+        from particles import all_particles
+        
+        self.dep_mass = [(part, part.mass) for part in all_particles \
+                            if part.pdg_code > 0 and \
+                                            part.mass not in list_of_parameters]
+        self.dep_width = [(part, part.width) for part in all_particles\
+                             if part.pdg_code > 0 and \
+                                part.width not in list_of_parameters]        
+    
+    @staticmethod
+    def order_param(obj1, obj2):
+        """ order parameter of a given block """
+        
+        maxlen = min([len(obj1.lhacode), len(obj2.lhacode)])
+    
+        for i in range(maxlen):
+            if obj1.lhacode[i] < obj2.lhacode[i]:
+                return -1
+            elif obj1.lhacode[i] == obj2.lhacode[i]:
+                return 0
+            else:
+                return 1
+        #identical up to the first finish
+        if len(obj1.lhacode) > len(obj2.lhacode):
+            return 1
+        elif  len(obj1.lhacode) == len(obj2.lhacode):
+            return 0
+        else:
+            return -1
+        
+    def write_card(self, all_ext_param):
+        """ """
+        
+        # list all lhablock
+        all_lhablock = set([param.lhablock for param in all_ext_param])
+        
+        # ordonate lhablock alphabeticaly
+        all_lhablock = list(all_lhablock)
+        all_lhablock.sort()
+        # put at the beginning SMINPUT + MASS + DECAY
+        for name in ['DECAY', 'MASS','SMINPUTS']:
+            if name in all_lhablock:
+                all_lhablock.remove(name)
+                all_lhablock.insert(0, name)
+        
+        for lhablock in all_lhablock:
+            self.write_block(lhablock)
+            need_writing = [ param for param in all_ext_param if \
+                                                     param.lhablock == lhablock]
+            need_writing.sort(self.order_param)
+            [self.write_param(param, lhablock) for param in need_writing]
+            
+            if self.generic_output:
+                if lhablock in ['MASS', 'DECAY']:
+                    self.write_dep_param_block(lhablock)
+
+        if self.generic_output:
+            self.write_qnumber()
+                               
+    def write_block(self, name):
+        """ write a comment for a block"""
+        
+        self.fsock.writelines(
+        """\n###################################""" + \
+        """\n## INFORMATION FOR %s""" % name.upper() +\
+        """\n###################################\n"""
+         )
+        if name!='DECAY':
+            self.fsock.write("""Block %s \n""" % name)
+
+    def write_param(self, param, lhablock):
+        
+        lhacode=' '.join(['%3s' % key for key in param.lhacode])
+        if lhablock != 'DECAY':
+            text = """  %s %e # %s \n""" % (lhacode, complex(param.value).real, param.name ) 
+        else:
+            text = '''DECAY %s %e \n''' % (lhacode, complex(param.value).real)
+        self.fsock.write(text) 
+                    
+
+
+    
+    def write_dep_param_block(self, lhablock):
+        import cmath
+        from parameters import all_parameters
+        for parameter in all_parameters:
+            exec("%s = %s" % (parameter.name, parameter.value))
+        text = "##  Not dependent paramater.\n"
+        text += "## Those values should be edited following analytical the \n"
+        text += "## analytical expression. Some generator could simply ignore \n"
+        text += "## those values and use the analytical expression\n"
+        
+        if lhablock == 'MASS':
+            data = self.dep_mass
+            prefix = " "
+        else:
+            data = self.dep_width
+            prefix = "DECAY "
+        for part, param in data:
+            if isinstance(param.value, str):
+                value = complex(eval(param.value)).real
+            else:
+                value = param.value
+            
+            text += """%s %s %f # %s : %s \n""" %(prefix, part.pdg_code, 
+                        value, part.name, param.value)
+        self.fsock.write(text)    
+    
+    sm_pdg = [1,2,3,4,5,6,11,12,13,13,14,15,16,21,22,23,24,25]
+    data="""Block QNUMBERS %(pdg)d  # %(name)s 
+        1 %(charge)d  # 3 times electric charge
+        2 %(spin)d  # number of spin states (2S+1)
+        3 %(color)d  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 %(antipart)d  # Particle/Antiparticle distinction (0=own anti)\n"""
+    
+    def write_qnumber(self):
+        """ write qnumber """
+        from particles import all_particles
+        import particles
+        print particles.__file__
+        text="""#===========================================================\n"""
+        text += """# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)\n"""
+        text += """#===========================================================\n\n"""
+        
+        for part in all_particles:
+            if part.pdg_code in self.sm_pdg or part.pdg_code < 0:
+                continue
+            text += self.data % {'pdg': part.pdg_code,
+                                 'name': part.name,
+                                 'charge': 3 * part.charge,
+                                 'spin': 2 * part.spin + 1,
+                                 'color': part.color,
+                                 'antipart': part.name != part.antiname and 1 or 0}
+        
+        self.fsock.write(text)
+        
+        
+            
+            
+            
+            
+        
+            
+if '__main__' == __name__:
+    ParamCardWriter('./param_card.dat', generic=True)
+    print 'write ./param_card.dat'
+