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

[Pierre Artoisenet]

Nature of the problem:

The decomposition into color subamplitudes is incorrect for most processes
in the diquark model

Origin of the problem:

Although the current implementation in the trunk gives the correct value for
the matrix element squared, there is a problem with the decomposition into
color subamplitudes: in some cases there are too many subamplitudes, and
and they are not gauge invariant.

One example: for g g > g six six~, the correct decomposition involves
only 13 subamplitudes, whereas the current version in the trunk returns
104 jamp's.

The origin of the problem is that the algorithm fails to recognize
identical color strings in some cases, because the summation
indices are named in a different way and/or the symmetry propeties of
the Clebsch Gordan coefficients K6 and K6Bar are not taken into account.
As a result, each color subamplitude in the correct decomposition appears
in many different formats in the current version of the code, and this
gives rise to a very large number of (non gauge-invariant) Jamp's.

In this branch the problem has been solved by forcing a specific convention
for the summation indices and for the i,j indices in K6 and K6Bar.

I added a test in the file tests/unit_tests/various/test_diquark_models.py

ATTENTION:

  1. I had to change the previous tests in test_diquark_models.py
     because the color flow decomposition did not match exactly the
     reference results after I modified color_amp.py and color_algebra.py.
     The results before/after the modifications were consitent with each
     other, but up to a re-labeling of indices (which is conventional).
     So I had to use another convention for labeling the indices in the
     reference color flow strucures in test_diquark_models.py.
     Johan: can you confirm that it is ok I edited your own tests ?

  2. The tests/test_manager.py -p -A failed BEFORE I introduced any modification
     in the code. After the modifs, I got exactly the same error report.
     I copy it below

  3. The branch is originally a copy from the trunk, but it runs in the debbug mode

Pierre

[Johan Alwall (johan-alwall)]

Hello Pierre,

Many thanks for your excellent report, and for your great work with this revision.
To answer your questions:
1. It's not a problem if tests need to be modified due to corrections (or simply due to changes) in the code, as long as these modifications are carefully reviewed and fully understood. In your case, this is clearly the case.
2. Could you please copy the error report from tests/test_manager.py -p A ?
3. Yes this is normal.

Comments on your revision:
Awesome job! My only comment is that perhaps you might want to rename the function order_summation_for_K to something generic like order_summation? I guess this same issue might happen for other color objects in the future, in which case the function can be extended also for those objects, rather than adding new functions. As long as it's clearly explained in the comments and doc strings, it should still be clear.
Please also add some comments about how the ordering works. Right now, it's quite difficult to tell which ordering you have chosen for the indices. Better to over-comment than to undercomment.

Many thanks,
Johan

review: Needs Fixing

[Pierre Artoisenet (partois)]

Hi Johan,

Thank you for your prompt reaction on this.
Sorry, I forgot to copy the error report from tests/test_manager.py -p A
Now it can be read below.

As you pointed out, the new routine may be generalized in the future, so following
your suggestion I have renamed the routine into "order_summation" and put
also some comments about the order that is chosen for the summation indices.

I commited the modifications using

bzr commit -m "description of the modifs"

but they do not appear on the webpage. Should I use again

bzr push lp:~maddevelopers/madgraph5/fix_color_basis_sextet

to force these modifications on the Launchpad (public version) ?

Sorry for these basic questions, and thank you for your help!

Cheers,
Pierre

F.............FE.....E...
======================================================================
ERROR: test_output_standalone_directory (tests.acceptance_tests.test_cmd.TestCmdShell2)
Test command 'output' with path
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/acceptance_tests/test_cmd.py", line 316, in test_output_standalone_directory
    self.do('output standalone %s' % self.out_dir)
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/acceptance_tests/test_cmd.py", line 154, in do
    self.cmd.exec_cmd(line)
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/interface/extended_cmd.py", line 171, in exec_cmd
    stop = cmd.Cmd.onecmd(self, line)
  File "/localstore/theorie/partois/lib/python2.7/cmd.py", line 219, in onecmd
    return func(arg)
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/interface/cmd_interface.py", line 2740, in do_output
    self.finalize(nojpeg)
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/interface/cmd_interface.py", line 2983, in finalize
    online)
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/iolibs/export_v4.py", line 766, in finalize_v4_directory
    self.make()
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/iolibs/export_v4.py", line 754, in make
    misc.compile(arg=['../lib/libdhelas.a'], cwd=source_dir, mode='fortran')
  File "/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/iolibs/misc.py", line 157, in compile
    raise MadGraph5Error, error_text
MadGraph5Error: A compilation Error occurs when trying to compile /localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/AUTO_TEST_MG5/Source.
The compilation fails with the following output message:
  cd DHELAS; make
  make[1]: Entering directory `/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/AUTO_TEST_MG5/Source/DHELAS'
  g77 -O -w -ffixed-line-length-132 -c -o httsxx.o httsxx.F
  httsxx.F: In subroutine `httsxx':
  httsxx.F:24:
           include "dimension.inc"
           ^
  Unable to open INCLUDE file `dimension.inc' at (^)
  httsxx.F:25:
           double complex tc1(DIM),htts(DIM),sc(DIM),tc2...

[lp:~maddevelopers/mg5amcnlo/fix_color_basis_sextet]

175. By Pierre Artoisenet on 2011-10-10

Added first line of text

176. By mattelaer-olivier on 2011-10-10

fix the error in the acceptance test. A very stupid mistake

[Olivier Mattelaer (olivier-mattelaer)]

I've fix one of the error of the acceptance test.
The remaining three are fixed in the new_install_cmd branch.
So my proposal is to merge this branch in that one.
and merge both of them in the trunk.

[Johan Alwall (johan-alwall)]

Hello Pierre,

Excellent work. I approve the merge.

To perform the merge, first merge with the latest trunk and new_install_command branch
bzr merge lp:madgraph5
Fix the conflicts (there should be in VERSION and UpdateNotes), and do
bzr resolve
bzr commit -m "Merged with latest version of trunk, v. 1.3.20"
Next, merge with the new_install_command branch:
bzr merge lp:~maddevelopers/madgraph5/new_install_command
Same thing, fix the conflicts (note that since we do both these branches together, you only need one version number, v. 1.3.21, and don't forget to set today's date in both VERSION and UpdateNotes), resolve and commit. Then push, and let us take a look before you merge with the trunk.

For the actual merge, do the following:
bzr branch lp:madgraph5
cd madgraph5
bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
bzr commit -m "Merged up to version 1.3.21 (the fix_color_basis_sextet and new_install_command branches)"
Now create a release:
bin/create_release.py
Note that you need epydoc and gpg installed on your computer (and gpg set up) to do this.
Once this is done, and there are no errors, you can push to lp:madgraph5.

Now replace the download file on Launchpad (click on "All downloads"), add the new lines from UpdateNotes to the release info and the release 1.3.0 information (look through the page and you'll see where to make changes), and update both UIUC and UCL cluster. This step is crucial, since otherwise the download links won't work on the clusters.

Please let us know if you have any questions, and also let us know once it's all done so we can double-check.

Thanks!
Johan

review: Approve

[Olivier Mattelaer (olivier-mattelaer)]

> Note that you need epydoc and gpg installed on your computer (and
> gpg set up) to do this.

Concerning this point. gpg is not mandatory (I never succeed to make
it work due to compatibility issue with Mac10.5)
while epydoc is 100% require (you can install it from http://epydoc.sourceforge.net/installing.html)

Otherwise note that the script create_release test that
1) all your modification are committed.
2) that your version number has comment in the Update Note.
3) that the date is today.
And create a warning/crash if they are a problem.
As you will note, this script will also automatically test the full
set of unit test and acceptance test.
If a test fails, then the .tar.gz will not be created

Cheers,

Olivier

On 10-oct.-11, at 20:03, Johan Alwall wrote:

> Review: Approve
>
> Hello Pierre,
>
> Excellent work. I approve the merge.
>
> To perform the merge, first merge with the latest trunk and
> new_install_command branch
> bzr merge lp:madgraph5
> Fix the conflicts (there should be in VERSION and UpdateNotes), and do
> bzr resolve
> bzr commit -m "Merged with latest version of trunk, v. 1.3.20"
> Next, merge with the new_install_command branch:
> bzr merge lp:~maddevelopers/madgraph5/new_install_command
> Same thing, fix the conflicts (note that since we do both these
> branches together, you only need one version number, v. 1.3.21, and
> don't forget to set today's date in both VERSION and UpdateNotes),
> resolve and commit. Then push, and let us take a look before you
> merge with the trunk.
>
> For the actual merge, do the following:
> bzr branch lp:madgraph5
> cd madgraph5
> bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
> bzr commit -m "Merged up to version 1.3.21 (the
> fix_color_basis_sextet and new_install_command branches)"
> Now create a release:
> bin/create_release.py
> Note that you need epydoc and gpg installed on your computer (and
> gpg set up) to do this.
> Once this is done, and there are no errors, you can push to
> lp:madgraph5.
>
> Now replace the download file on Launchpad (click on "All
> downloads"), add the new lines from UpdateNotes to the release info
> and the release 1.3.0 information (look through the page and you'll
> see where to make changes), and update both UIUC and UCL cluster.
> This step is crucial, since otherwise the download links won't work
> on the clusters.
>
> Please let us know if you have any questions, and also let us know
> once it's all done so we can double-check.
>
> Thanks!
> Johan
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/
> fix_color_basis_sextet/+merge/78798
> You are subscribed to branch lp:madgraph5.

[Pierre Artoisenet (partois)]

Hi Johan, hi Olivier,

many thanks for the clear information below!

The merge with the version in the trunk (lp:madgraph5) went fine, I think,
but I have some problems for the merge with the new_install_command
branch:

The first problem is that I am not sure how to solve the conflict in
UpdateNotes, since there are 2 version 1.3.20 :

<<<<<<< TREE
1.3.21 (11/10/11) PA: fix the problem of degenerate color basis in the
                       diquark sextet model

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

The second problem is that "test_manager.py" is not working properly
after this second merge:

./tests/test_manager.py
Traceback (most recent call last):
   File "./tests/test_manager.py", line 377, in <module>
     package=options.path)
   File "./tests/test_manager.py", line 70, in run
     re_opt=re_opt):
   File "./tests/test_manager.py", line 117, in __iter__
     self._check_if_obj_build()
   File "./tests/test_manager.py", line 113, in _check_if_obj_build
     self.collect_dir(self.package, checking=True)
   File "./tests/test_manager.py", line 149, in collect_dir
     self.collect_dir(directory + '/' + name, local_check)
   File "./tests/test_manager.py", line 147, in collect_dir
     self.collect_file(directory + '/' + name, local_check)
   File "./tests/test_manager.py", line 158, in collect_file
     __import__(pyname)
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/unit_tests/interface/test_cmd.py",
line 23, in <module>
     class TestValidCmd(unittest.TestCase):
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/unit_tests/interface/test_cmd.py",
line 26, in TestValidCmd
     cmd = cmd.MadGraphCmdShell()
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/interface/cmd_interface.py",
line 1507, in __init__
     CmdExtended.__init__(self, *completekey, **stdin)
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/interface/cmd_interface.py",
line 111, in __init__
     info = misc.get_pkg_info()
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/iolibs/misc.py",
line 64, in get_pkg_info
     parse_info_str)
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/madgraph/iolibs/files.py",
line 40, in read_from_file
     except IOError, (errno, strerror):
ValueError: need more than 1 value to unpack
...

[Olivier Mattelaer (olivier-mattelaer)]

Hi Pierre,

Concerning the conflict:
The main idea is the same as cvs
in your file you have
tag <<<<<
=====
 >>>>>
delimitting conflict area.
In order to help you, you have three aditional file
NAME.this (your file before the merge)
NAM.other (the file with which you are doing the merge)
NAME.base (the common version of the file on which are done the
merging comparison)
Note that in 99% of the cases, you don't need those files.

The idea is to suppress from the file all tag
 >>>>>> ===== and <<<<<
(as in cvs) and have replace the area with the correct part off course.

When this is done, you can do
bzr resolve
and if you done forget any conflict this will work. (otherwise, it
will indicate which files still contains conflict)

For the test, manager, that's a problem linked to the conflict, when
those will be resolve, the test manager should pass.

Cheers,

Olivier

On 11-oct.-11, at 02:43, Pierre Artoisenet wrote:

>
> Hi Johan, hi Olivier,
>
> many thanks for the clear information below!
>
> The merge with the version in the trunk (lp:madgraph5) went fine, I
> think,
> but I have some problems for the merge with the new_install_command
> branch:
>
> The first problem is that I am not sure how to solve the conflict in
> UpdateNotes, since there are 2 version 1.3.20 :
>
> <<<<<<< TREE
> 1.3.21 (11/10/11) PA: fix the problem of degenerate color basis in the
> diquark sextet model
>
> 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
>
>
> The second problem is that "test_manager.py" is not working properly
> after this second merge:
>
> ./tests/test_manager.py
> Traceback (most recent call last):
> File "./tests/test_manager.py", line 377, in <module>
> package=options.path)
> File "./tests/test_manager.py", line 70, in run
> re_opt=re_opt):
> File "./tests/test_manager.py", line 117, in __iter__
> self._check_if_obj_build()
> File "./tests/test_manager.py", line 113, in _check_if_obj_build
> self.collect_dir(self.package, checking=True)
> File "./tests/test_manager.py", line 149, in collect_dir
> self.collect_dir(directory + '/' + name, local_check)
> File "./tests/test_manager.py", line 147, in collect_dir
> self.collect_file(directory + '/' + name, local_check)
> File "./tests/test_manager.py", line 158, in collect_file
> __import__(pyname)
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/unit_tests/interface/test_cmd....

[Pierre Artoisenet (partois)]

Hi Olivier,

Thanks for these info, now the tests are indeed running.

If I remove all tags >>>>>> ===== and <<<<<,
I will have two versions 1.3.20. Should I leave the UpdateNotes.txt
file with these two versions 1.3.20 ?

Also now I have got two error reports:

./tests/test_manager.py
..............................................................................................................................................................................................................................................................................................................................................F..........................................................................................................................................
======================================================================
FAIL: test_get_momenta
(tests.unit_tests.various.test_process_checks.TestMatrixElementChecker)
Test the get_momenta function
----------------------------------------------------------------------
Traceback (most recent call last):
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/unit_tests/various/test_process_checks.py",
line 76, in test_get_momenta
     self.assertAlmostEqual(mass, 0., 8)
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/unit_tests/__init__.py",
line 52, in assertAlmostEqual
     unittest.TestCase.assertAlmostEqual(self, a, b, *opt, **arg)
AssertionError: -6.83940015733e-09 != 0.0 within 8 digits

----------------------------------------------------------------------
Ran 473 tests in 37.643s

FAILED (failures=1)

./tests/test_manager.py -p -A
Traceback (most recent call last):
   File "./tests/test_manager.py", line 377, in <module>
     package=options.path)
   File "./tests/test_manager.py", line 70, in run
     re_opt=re_opt):
   File "./tests/test_manager.py", line 117, in __iter__
     self._check_if_obj_build()
   File "./tests/test_manager.py", line 113, in _check_if_obj_build
     self.collect_dir(self.package, checking=True)
   File "./tests/test_manager.py", line 134, in collect_dir
     for name in os.listdir(os.path.join(root_path,directory)):
OSError: [Errno 2] No such file or directory:
'/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/-A/'
[partois@clain fix_color_basis_sextet]$ ./tests/test_manager.py -p A
F........................
======================================================================
FAIL: test_config (tests.acceptance_tests.test_cmd.TestCmdShell1)
check that configuration file is at default value
----------------------------------------------------------------------
Traceback (most recent call last):
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/acceptance_tests/test_cmd.py",
line 127, in test_config
     self.assertEqual(launch_ext.open_file.web_browser, 'gv')
AssertionError: 'firefox' != 'gv'

----------------------------------------------------------------------
Ran 25 tests in 71.081s

FAILED (failures=1)

Cheers,
Pierre

On Tue, 11 Oct 2011, Olivier Mattelaer wrote:

> Hi Pierre,
>
> Concerning the conflict:
> The main idea is th...

[Johan Alwall (johan-alwall)]

Hello Pierre,

For the UpdateNotes: The only version which is released is (of course) the trunk. Any other version numbers are temporary and irrelevant. As I wrote in my explanation above, you should combine the two versions into one release:

1.3.21 (11/10/11) PA: fix the problem of degenerate color basis in the
                       diquark sextet model
                   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

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

Cheers,
Johan

[Pierre Artoisenet (partois)]

Hi Johan,

Ok I understand now... I made the modifications in the
UpdateNotes.txt file, and the test ./tests/test_manager.py ran
successfully (magic!)

Unfortunately I still have a problem with the second test:

./tests/test_manager.py -p A
F........................
======================================================================
FAIL: test_config (tests.acceptance_tests.test_cmd.TestCmdShell1)
check that configuration file is at default value
----------------------------------------------------------------------
Traceback (most recent call last):
   File
"/localstore/theorie/partois/madgraph5_branches/fix_color_basis_sextet/tests/acceptance_tests/test_cmd.py",
line 127, in test_config
     self.assertEqual(launch_ext.open_file.web_browser, 'gv')
AssertionError: 'firefox' != 'gv'

----------------------------------------------------------------------
Ran 25 tests in 73.986s

FAILED (failures=1)

Cheers,
Pierre

On Tue, 11 Oct 2011, Johan Alwall wrote:

> Hello Pierre,
>
> For the UpdateNotes: The only version which is released is (of course) the trunk. Any other version numbers are temporary and irrelevant. As I wrote in my explanation above, you should combine the two versions into one release:
>
> 1.3.21 (11/10/11) PA: fix the problem of degenerate color basis in the
> diquark sextet model
> 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
>
> 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.
> ...
>
> Cheers,
> Johan
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/fix_color_basis_sextet/+merge/78798
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/madgraph5/fix_color_basis_sextet.
>

[Olivier Mattelaer (olivier-mattelaer)]

Concerning this error:

> Traceback (most recent call last):
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/unit_tests/various/
> test_process_checks.py",
> line 76, in test_get_momenta
> self.assertAlmostEqual(mass, 0., 8)
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/unit_tests/__init__.py",
> line 52, in assertAlmostEqual
> unittest.TestCase.assertAlmostEqual(self, a, b, *opt, **arg)
> AssertionError: -6.83940015733e-09 != 0.0 within 8 digits

The test is in fact fine but probably just above the trigger of the
error. Since this test is base on a random number, this might be pass
the next time.

> ======================================================================
> FAIL: test_config (tests.acceptance_tests.test_cmd.TestCmdShell1)
> check that configuration file is at default value
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/acceptance_tests/test_cmd.py",
> line 127, in test_config
> self.assertEqual(launch_ext.open_file.web_browser, 'gv')
> AssertionError: 'firefox' != 'gv'

Did you pull my modifications on this test?
it looks not.

My advice is to do the following
a) bzr commit
b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
This one will probably fail, if so, you are then force to run
bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
and then see if the test pass.

Cheers,

olivier

On 11-oct.-11, at 03:38, Pierre Artoisenet wrote:

>
> Hi Olivier,
>
> Thanks for these info, now the tests are indeed running.
>
> If I remove all tags >>>>>> ===== and <<<<<,
> I will have two versions 1.3.20. Should I leave the UpdateNotes.txt
> file with these two versions 1.3.20 ?
>
>
> Also now I have got two error reports:
>
> ./tests/test_manager.py
> ..............................................................................................................................................................................................................................................................................................................................................F..........................................................................................................................................
> ======================================================================
> FAIL: test_get_momenta
> (tests
> .unit_tests.various.test_process_checks.TestMatrixElementChecker)
> Test the get_momenta function
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/unit_tests/various/
> test_process_checks.py",
> line 76, in test_get_momenta
> self.assertAlmostEqual(mass, 0., 8)
> File
> "/localstore/theorie/partois/madgraph5_branches/
> fix_color_basis_sextet/tests/unit_tests/__init__.py",
> line 52, in assertAlmostEqual
> unittest.TestCase.assertAlmostEqual(self, a, b, *opt, **arg)
> Asserti...

[Johan Alwall (johan-alwall)]

> My advice is to do the following
> a) bzr commit
> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
> This one will probably fail, if so, you are then force to run
> bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
> and then see if the test pass.

If you do a merge, you lose all commit messages between the commit and the merge, which looks quite ugly.

In fact, the correct way of doing this is
a) do NOT commit but save all files if you have them open
b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
if this fails, do
bzr uncommit
try to pull again - if it still doesn't work, do another uncommit, etc. until the pull works.
Finally, commit (and make sure you retain all commit messages that you have uncommitted, as well as describe additional changes).

Cheers,
Johan

[Pierre Artoisenet (partois)]

Ok, I eventually went through the merge procedure with the lattest trunk
and the new_install_command branch.

I am waiting for your green light to do the next step (update the version
in the trunk)

Cheers,
Pierre

On Tue, 11 Oct 2011, Johan Alwall wrote:

>> My advice is to do the following
>> a) bzr commit
>> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>> This one will probably fail, if so, you are then force to run
>> bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>> and then see if the test pass.
>
> If you do a merge, you lose all commit messages between the commit and the merge, which looks quite ugly.
>
> In fact, the correct way of doing this is
> a) do NOT commit but save all files if you have them open
> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
> if this fails, do
> bzr uncommit
> try to pull again - if it still doesn't work, do another uncommit, etc. until the pull works.
> Finally, commit (and make sure you retain all commit messages that you have uncommitted, as well as describe additional changes).
>
> Cheers,
> Johan
>
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/fix_color_basis_sextet/+merge/78798
> Your team MadTeam is subscribed to branch lp:madgraph5.
>

[Olivier Mattelaer (olivier-mattelaer)]

Ok for that we will need to see the last version.
So could you push it on the web.
(just make bzr push)

Cheers,

olivier

On Oct 11, 2011, at 4:37 AM, Pierre Artoisenet wrote:

>
> Ok, I eventually went through the merge procedure with the lattest
> trunk
> and the new_install_command branch.
>
> I am waiting for your green light to do the next step (update the
> version
> in the trunk)
>
> Cheers,
> Pierre
>
> On Tue, 11 Oct 2011, Johan Alwall wrote:
>
>>> My advice is to do the following
>>> a) bzr commit
>>> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>>> This one will probably fail, if so, you are then force to run
>>> bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>>> and then see if the test pass.
>>
>> If you do a merge, you lose all commit messages between the commit
>> and the merge, which looks quite ugly.
>>
>> In fact, the correct way of doing this is
>> a) do NOT commit but save all files if you have them open
>> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>> if this fails, do
>> bzr uncommit
>> try to pull again - if it still doesn't work, do another uncommit,
>> etc. until the pull works.
>> Finally, commit (and make sure you retain all commit messages that
>> you have uncommitted, as well as describe additional changes).
>>
>> Cheers,
>> Johan
>>
>> --
>> https://code.launchpad.net/~maddevelopers/madgraph5/
>> fix_color_basis_sextet/+merge/78798
>> Your team MadTeam is subscribed to branch lp:madgraph5.
>>
>
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/
> fix_color_basis_sextet/+merge/78798
> You are subscribed to branch lp:madgraph5.

[lp:~maddevelopers/mg5amcnlo/fix_color_basis_sextet]

177. By Pierre Artoisenet on 2011-10-11

Merged with the new_install_command branch

[Pierre Artoisenet (partois)]

ah yes, sorry, I have just done it

Cheers,
Pierre

On Tue, 11 Oct 2011, Olivier Mattelaer wrote:

> Ok for that we will need to see the last version.
> So could you push it on the web.
> (just make bzr push)
>
> Cheers,
>
> olivier
>
>
> On Oct 11, 2011, at 4:37 AM, Pierre Artoisenet wrote:
>
>>
>> Ok, I eventually went through the merge procedure with the lattest
>> trunk
>> and the new_install_command branch.
>>
>> I am waiting for your green light to do the next step (update the
>> version
>> in the trunk)
>>
>> Cheers,
>> Pierre
>>
>> On Tue, 11 Oct 2011, Johan Alwall wrote:
>>
>>>> My advice is to do the following
>>>> a) bzr commit
>>>> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>>>> This one will probably fail, if so, you are then force to run
>>>> bzr merge lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>>>> and then see if the test pass.
>>>
>>> If you do a merge, you lose all commit messages between the commit
>>> and the merge, which looks quite ugly.
>>>
>>> In fact, the correct way of doing this is
>>> a) do NOT commit but save all files if you have them open
>>> b) bzr pull lp:~maddevelopers/madgraph5/fix_color_basis_sextet
>>> if this fails, do
>>> bzr uncommit
>>> try to pull again - if it still doesn't work, do another uncommit,
>>> etc. until the pull works.
>>> Finally, commit (and make sure you retain all commit messages that
>>> you have uncommitted, as well as describe additional changes).
>>>
>>> Cheers,
>>> Johan
>>>
>>> --
>>> https://code.launchpad.net/~maddevelopers/madgraph5/
>>> fix_color_basis_sextet/+merge/78798
>>> Your team MadTeam is subscribed to branch lp:madgraph5.
>>>
>>
>> --
>> https://code.launchpad.net/~maddevelopers/madgraph5/
>> fix_color_basis_sextet/+merge/78798
>> You are subscribed to branch lp:madgraph5.
>
>
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/fix_color_basis_sextet/+merge/78798
> Your team MadDevelopers is subscribed to branch lp:~maddevelopers/madgraph5/fix_color_basis_sextet.
>

[Johan Alwall (johan-alwall)]

Hello Pierre,

Only one comment: Please put main modifs (such as new functionality) first in the UpdateNotes, and bug fixes at the end. The new install command is a pretty major thing, so it deserves to be first.

Cheers,
Johan

[lp:~maddevelopers/mg5amcnlo/fix_color_basis_sextet]

178. By Pierre Artoisenet on 2011-10-12

change the order of modifs in the UpdateNotes.txt

[Pierre Artoisenet (partois)]

Hi Johan,

This is done, I have just made a "bzr push"

Cheers,
Pierre

On Wed, 12 Oct 2011, Johan Alwall wrote:

> Hello Pierre,
>
> Only one comment: Please put main modifs (such as new functionality) first in the UpdateNotes, and bug fixes at the end. The new install command is a pretty major thing, so it deserves to be first.
>
> Cheers,
> Johan
> --
> https://code.launchpad.net/~maddevelopers/madgraph5/fix_color_basis_sextet/+merge/78798
> Your team MadTeam is subscribed to branch lp:madgraph5.
>

[lp:~maddevelopers/mg5amcnlo/fix_color_basis_sextet]

179. By Johan Alwall on 2011-10-12

Merged with fix_bug_from_Benj: Fix bug in cluster.f in case of incompatible t-channel processes.

=== modified file 'HELAS/hstlxx.F'
--- HELAS/hstlxx.F	2010-10-30 03:26:37 +0000
+++ HELAS/hstlxx.F	2011-10-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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/SubProcesses/cluster.f'
--- Template/SubProcesses/cluster.f	2011-09-29 11:31:12 +0000
+++ Template/SubProcesses/cluster.f	2011-10-12 10:39:26 +0000
@@ -174,6 +174,7 @@
       include 'nexternal.inc'
       include 'maxamps.inc'
       include 'cluster.inc'
+      include 'message.inc'
       integer ignum, idij, iproc, i
 
       if(idij.gt.n_max_cl) return
@@ -182,7 +183,8 @@
       enddo
       id_cl(iproc,idij,0)=id_cl(iproc,idij,0)+1
       id_cl(iproc,idij,id_cl(iproc,idij,0))=ignum
-c      print *,'Adding graph ',ignum,' to prop ',idij
+      if(btest(mlevel,5))
+     $     write(*,*)'Adding graph ',ignum,' to prop ',idij,' for proc ',iproc
       return
       end
 
@@ -226,6 +228,8 @@
       integer pow(-nexternal:0,lmaxconfigs)
       logical first_time
       save pmass,pwidth,pow
+      INTEGER CONFSUB(MAXSPROC,LMAXCONFIGS)
+      INCLUDE 'config_subproc_map.inc'
       data first_time /.true./
 
       integer combid
@@ -259,6 +263,7 @@
 c     Set pdg code for propagator
                   do l=1,2
                      do iproc=1,maxsproc
+                     if(confsub(iproc,ignum).eq.0) cycle
                      if(sprop(iproc,k,ignum).ne.0)then
                         ipdgcl(icmp(l),ignum,iproc)=sprop(iproc,k,ignum)
 c                       If this is radiation off heavy FS particle, set heavyrad to true
@@ -276,7 +281,6 @@
                         ipdgcl(icmp(l),ignum,iproc)=ipdgcl(2,ignum,iproc)
                      else
                         ipdgcl(icmp(l),ignum,iproc)=0
-                        cycle
                      endif
                      if(btest(mlevel,4))
      $                    write(*,*) 'add table entry for (',ipids(i,1,ipnum),
@@ -471,6 +475,7 @@
       include 'nexternal.inc'
       include 'maxamps.inc'
       include 'cluster.inc'
+      include 'message.inc'
 
       integer idij,nbw,ibwlist(nexternal),icgs(0:n_max_cg)
       logical foundbw
@@ -504,6 +509,9 @@
          if (icgs(0).gt.0)then
            findmt=.true.
          endif
+         if (btest(mlevel,5))
+     $        write (*,*)'findmt: ',findmt,' IPROC=',IPROC,' icgs(0)=',icgs(0),
+     $        ' icgs = ',(icgs(i),i=1,icgs(0))
          return
       else
 c     Check for common graphs
@@ -578,6 +586,8 @@
       enddo
 c     Check if any resonances are on the BW, store results in to_checkbw
       call checkbw(nbw,ibwlist,isbw)
+      if(btest(mlevel,4).and.nbw.gt.0)
+     $     write(*,*) 'Found BWs: ',(ibwlist(i),i=1,nbw)
 
 c     initialize index map
       do i=1,nexternal

=== modified file 'Template/bin/plot'
--- Template/bin/plot	2011-05-24 04:27:02 +0000
+++ Template/bin/plot	2011-10-12 10:39:26 +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-12 10:39:26 +0000
@@ -1,5 +1,18 @@
 Update notes for MadGraph 5 (in reverse time order)
 
+1.3.21 (12/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
+                  PA: fix the problem of degenerate color basis in the
+                      diquark sextet model
+                  JA: Fixed bug in cluster.f that created a hard stop,
+                      introduced in 1.3.18.
+
 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.

=== modified file 'aloha/aloha_writers.py'
--- aloha/aloha_writers.py	2011-09-30 13:55:43 +0000
+++ aloha/aloha_writers.py	2011-10-12 10:39:26 +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-12 10:39:26 +0000
@@ -1,3 +1,2 @@
-version = 1.3.20
-date = 2011-10-09
-
+version = 1.3.21
+date = 2011-10-12

=== modified file 'madgraph/core/color_algebra.py'
--- madgraph/core/color_algebra.py	2010-12-15 07:28:25 +0000
+++ madgraph/core/color_algebra.py	2011-10-12 10:39:26 +0000
@@ -757,6 +757,64 @@
                                             1, int(Nc ** abs(self.Nc_power))),
                     self.is_imaginary)
 
+    def order_summation(self, immutable=None):
+        """Force a specific order for the summation indices 
+           in case we have Clebsch Gordan coefficients K6's or K6Bar's
+           This is necessary to correctly recognize later on the equivalent 
+           color strings (otherwise the color basis is degenerate)
+           The new ordering is as follow:
+                1. put K and KBar Clebsch Gordan coefficients at the end of the list of color factors
+                   the other factors are re-arranged in the reversed order compared with immutable 
+                2. rename the summation indices so that they are increasing (starting from 10000)
+                   from left to right
+                3. finally, after the summation indices have been renamed, replace
+                   K6(a,i,j) by K6(a,j,i) and K6Bar(a,i,j) by K6Bar(a,j,i) IF j>i
+        """
+
+        if not immutable:
+            immutable = self.to_immutable()
+
+#       STEP 1: first scan to see whether there are some K's or KBar's,
+#       and put them at the en 
+        immutable_order2=[]
+        go_further=0
+        for  elem in immutable:
+          if elem[0]=="K6" or elem[0]=="K6Bar" :
+             immutable_order2.append(elem)
+             go_further=1
+          else: immutable_order2.insert(0,elem)
+
+        if go_further==0: return
+
+#       STEP 2: rename the summation indices so that they are increasing (starting from 10000)
+#               from left to right
+        replaced_indices = {}
+        curr_ind = 10000
+        return_list = []
+
+        for elem in immutable_order2:
+            can_elem = [elem[0], []]
+            for index in elem[1]:
+              if index>9999:  # consider only summation indices
+                try:
+                    new_index = replaced_indices[index]
+                except KeyError:
+                    new_index = curr_ind
+                    curr_ind += 1
+                    replaced_indices[index] = new_index
+              else: new_index=index
+              can_elem[1].append(new_index)
+#       STEP 3.  replace K6(a,i,j) by K6(a,j,i) and K6Bar(a,i,j) by K6Bar(a,j,i) IF j>i
+            if (can_elem[0]=="K6" or can_elem[0]=="K6Bar"):
+               if can_elem[1][2]>can_elem[1][1]: can_elem[1]=[can_elem[1][0], can_elem[1][2], can_elem[1][1] ]
+            return_list.append((can_elem[0], tuple(can_elem[1])))
+        return_list.sort()
+
+        self.from_immutable(return_list)
+        self.immutable=None   # don't use the information self.immutable later on in the code, 
+                              # since the summation indices have been modified 
+        return
+
     def to_canonical(self, immutable=None):
         """Returns the canonical representation of the immutable representation 
         (i.e., first index is 1, ...). This allow for an easy comparison of

=== modified file 'madgraph/core/color_amp.py'
--- madgraph/core/color_amp.py	2011-05-16 19:41:08 +0000
+++ madgraph/core/color_amp.py	2011-10-12 10:39:26 +0000
@@ -229,6 +229,10 @@
                 col_fact = color_algebra.ColorFactor([col_str])
                 col_fact = col_fact.full_simplify()
 
+                # Here we need to force a specific order for the summed indices
+                # in case we have K6 or K6bar Clebsch Gordan coefficients
+                for colstr in col_fact: colstr.order_summation()
+
                 # Save the result for further use
                 canonical_col_fact = col_fact.create_copy()
                 canonical_col_fact.replace_indices(rep_dict)
@@ -250,6 +254,10 @@
                 # Must simplify once to put traces in a canonical ordering
                 col_fact = col_fact.simplify()
 
+                # Here we need to force a specific order for the summed indices
+                # in case we have K6 or K6bar Clebsch Gordan coefficients
+                for colstr in col_fact: colstr.order_summation()
+
             # loop over color strings in the resulting color factor
             for col_str in col_fact:
                 immutable_col_str = col_str.to_immutable()

=== 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-12 10:39:26 +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 

=== modified file 'madgraph/iolibs/export_v4.py'
--- madgraph/iolibs/export_v4.py	2011-10-04 04:15:47 +0000
+++ madgraph/iolibs/export_v4.py	2011-10-12 10:39:26 +0000
@@ -1001,6 +1001,10 @@
             writers.FortranWriter(filename),
             matrix_element)
 
+        filename = 'config_subproc_map.inc'
+        self.write_config_subproc_map_file(writers.FortranWriter(filename),
+                                           s_and_t_channels)
+
         filename = 'coloramps.inc'
         self.write_coloramps_file(writers.FortranWriter(filename),
                              mapconfigs,
@@ -1640,6 +1644,23 @@
         return s_and_t_channels
 
     #===========================================================================
+    # write_config_subproc_map_file
+    #===========================================================================
+    def write_config_subproc_map_file(self, writer, s_and_t_channels):
+        """Write a dummy config_subproc.inc file for MadEvent"""
+
+        lines = []
+
+        for iconfig in range(len(s_and_t_channels)):
+            lines.append("DATA CONFSUB(1,%d)/1/" % \
+                         (iconfig + 1))
+
+        # Write the file
+        writer.writelines(lines)
+
+        return True
+
+    #===========================================================================
     # write_decayBW_file
     #===========================================================================
     def write_decayBW_file(self, writer, s_and_t_channels):

=== 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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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-12 10:39:26 +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'
+    

=== modified file 'tests/acceptance_tests/test_cmd.py'
--- tests/acceptance_tests/test_cmd.py	2011-09-18 02:38:48 +0000
+++ tests/acceptance_tests/test_cmd.py	2011-10-12 10:39:26 +0000
@@ -120,11 +120,12 @@
         if sys.platform == 'darwin':
             self.assertEqual(launch_ext.open_file.web_browser, None)
             self.assertEqual(launch_ext.open_file.text_editor, text_editor)
-            self.assertEqual(launch_ext.open_file.web_browser, None)
+            self.assertEqual(launch_ext.open_file.eps_viewer, None)
         else:
             self.assertEqual(launch_ext.open_file.web_browser, 'firefox')
             self.assertEqual(launch_ext.open_file.text_editor, text_editor)
-            self.assertEqual(launch_ext.open_file.web_browser, 'gv')            
+            self.assertEqual(launch_ext.open_file.eps_viewer, 'gv')
+                        
 class TestCmdShell2(unittest.TestCase,
                     test_file_writers.CheckFileCreate):
     """Test all command line related to MG_ME"""

=== modified file 'tests/unit_tests/iolibs/test_export_v4.py'
--- tests/unit_tests/iolibs/test_export_v4.py	2011-10-03 04:26:01 +0000
+++ tests/unit_tests/iolibs/test_export_v4.py	2011-10-12 10:39:26 +0000
@@ -3357,6 +3357,58 @@
       DATA MAPCONFIG(0)/42/
 """)
 
+        # Test dummy config_subproc_map.inc file
+        writer = writers.FortranWriter(self.give_pos('test'))
+        exporter.write_config_subproc_map_file(writer, s_and_t_channels)
+        writer.close()
+
+        self.assertFileContains('test',
+"""      DATA CONFSUB(1,1)/1/
+      DATA CONFSUB(1,2)/1/
+      DATA CONFSUB(1,3)/1/
+      DATA CONFSUB(1,4)/1/
+      DATA CONFSUB(1,5)/1/
+      DATA CONFSUB(1,6)/1/
+      DATA CONFSUB(1,7)/1/
+      DATA CONFSUB(1,8)/1/
+      DATA CONFSUB(1,9)/1/
+      DATA CONFSUB(1,10)/1/
+      DATA CONFSUB(1,11)/1/
+      DATA CONFSUB(1,12)/1/
+      DATA CONFSUB(1,13)/1/
+      DATA CONFSUB(1,14)/1/
+      DATA CONFSUB(1,15)/1/
+      DATA CONFSUB(1,16)/1/
+      DATA CONFSUB(1,17)/1/
+      DATA CONFSUB(1,18)/1/
+      DATA CONFSUB(1,19)/1/
+      DATA CONFSUB(1,20)/1/
+      DATA CONFSUB(1,21)/1/
+      DATA CONFSUB(1,22)/1/
+      DATA CONFSUB(1,23)/1/
+      DATA CONFSUB(1,24)/1/
+      DATA CONFSUB(1,25)/1/
+      DATA CONFSUB(1,26)/1/
+      DATA CONFSUB(1,27)/1/
+      DATA CONFSUB(1,28)/1/
+      DATA CONFSUB(1,29)/1/
+      DATA CONFSUB(1,30)/1/
+      DATA CONFSUB(1,31)/1/
+      DATA CONFSUB(1,32)/1/
+      DATA CONFSUB(1,33)/1/
+      DATA CONFSUB(1,34)/1/
+      DATA CONFSUB(1,35)/1/
+      DATA CONFSUB(1,36)/1/
+      DATA CONFSUB(1,37)/1/
+      DATA CONFSUB(1,38)/1/
+      DATA CONFSUB(1,39)/1/
+      DATA CONFSUB(1,40)/1/
+      DATA CONFSUB(1,41)/1/
+      DATA CONFSUB(1,42)/1/
+""")
+
+        #print open(self.give_pos('test')).read()
+
         # Test coloramps.inc output
         self.assertEqual("\n".join(\
                        exporter.get_icolamp_lines(mapconfigs,
@@ -3660,6 +3712,18 @@
       DATA MAPCONFIG(0)/3/
 """)
 
+        # Test dummy config_subproc_map.inc file
+        writer = writers.FortranWriter(self.give_pos('test'))
+        exporter.write_config_subproc_map_file(writer,
+                                               s_and_t_channels)
+        writer.close()
+        #print open(self.give_pos('test')).read()
+        self.assertFileContains('test',
+"""      DATA CONFSUB(1,1)/1/
+      DATA CONFSUB(1,2)/1/
+      DATA CONFSUB(1,3)/1/
+""")
+
     def test_generate_helas_diagrams_uu_susu(self):
         """Testing the helas diagram generation u u > su su with t-channel n1
         """

=== modified file 'tests/unit_tests/various/test_diquark_models.py'
--- tests/unit_tests/various/test_diquark_models.py	2011-04-15 04:15:48 +0000
+++ tests/unit_tests/various/test_diquark_models.py	2011-10-12 10:39:26 +0000
@@ -181,11 +181,11 @@
                                                                       2)
         self.assertEqual(len(col_flow), len(matrix_element.get('color_basis')))
         self.assertEqual(col_flow,
-                         [{1: [502, 0], 2: [0, 501],
-                           3: [502, -503], 4: [-503, 501]},
-                          {1: [502, 0], 2: [0, 501],
-                           3: [502, -503], 4: [-503, 501]},
-                          {1: [501, 0], 2: [0, 501],
+                         [{1: [503, 0], 2: [0, 501], 
+                           3: [502, -503], 4: [-501, 502]}, 
+                          {1: [502, 0], 2: [0, 501], 
+                           3: [502, -503], 4: [-503, 501]}, 
+                          {1: [501, 0], 2: [0, 501], 
                            3: [502, -503], 4: [-502, 503]}])
 
         # Test u u > six g
@@ -223,14 +223,14 @@
                                                                       2)
         self.assertEqual(len(col_flow), len(matrix_element.get('color_basis')))
         self.assertEqual(col_flow,
-                         [{1: [501, 0], 2: [502, 0],
-                           3: [501, -503], 4: [502, 503]},
-                          {1: [502, 0], 2: [503, 0],
-                           3: [501, -503], 4: [502, 501]},
-                          {1: [501, 0], 2: [502, 0],
-                           3: [501, -503], 4: [502, 503]},
-                          {1: [502, 0], 2: [501, 0],
-                           3: [501, -503], 4: [502, 503]}])
+                          [{1: [503, 0], 2: [502, 0], 
+                            3: [501, -503], 4: [502, 501]}, 
+                           {1: [502, 0], 2: [503, 0], 
+                            3: [501, -503], 4: [502, 501]}, 
+                           {1: [501, 0], 2: [502, 0], 
+                            3: [501, -503], 4: [502, 503]}, 
+                           {1: [502, 0], 2: [501, 0], 
+                            3: [501, -503], 4: [502, 503]}])
 
         # Test u u > six > u u g
 
@@ -267,14 +267,14 @@
                                                                       2)
         self.assertEqual(len(col_flow), len(matrix_element.get('color_basis')))
         self.assertEqual(col_flow,
-                         [{1: [501, 0], 2: [503, 0], 3: [501, 0],
-                           4: [502, 0], 5: [503, 502]},
-                          {1: [503, 0], 2: [501, 0], 3: [501, 0],
-                           4: [502, 0], 5: [503, 502]},
-                          {1: [502, 0], 2: [503, 0], 3: [501, 0],
-                           4: [502, 0], 5: [503, 501]},
-                          {1: [503, 0], 2: [502, 0], 3: [501, 0],
-                           4: [502, 0], 5: [503, 501]}])
+                          [{1: [501, 0], 2: [503, 0], 3: [501, 0], 
+                            4: [502, 0], 5: [503, 502]}, 
+                           {1: [503, 0], 2: [501, 0], 3: [501, 0], 
+                            4: [502, 0], 5: [503, 502]}, 
+                           {1: [502, 0], 2: [503, 0], 3: [501, 0], 
+                            4: [502, 0], 5: [503, 501]}, {1: [503, 0], 
+                            2: [502, 0], 3: [501, 0], 
+                            4: [502, 0], 5: [503, 501]}])
 
         # Test six > u u
 
@@ -304,7 +304,31 @@
                                                                       2)
         self.assertEqual(len(col_flow), len(matrix_element.get('color_basis')))
         self.assertEqual(col_flow,
-                         [{1: [502, -501], 2: [0, 501], 3: [502, 0]}])
+                          [{1: [501, -502], 2: [0, 501], 3: [502, 0]}])
+
+        # Test the size of the color basis in g g  > g six six 
+
+        myleglist = base_objects.LegList()
+
+
+        myleglist.append(base_objects.Leg({'id':21,
+                                           'state':False}))
+        myleglist.append(base_objects.Leg({'id':21,
+                                           'state':False}))
+        myleglist.append(base_objects.Leg({'id':21,
+                                           'state':True}))
+        myleglist.append(base_objects.Leg({'id':9000006,
+                                           'state':True}))
+        myleglist.append(base_objects.Leg({'id':-9000006,
+                                           'state':True}))
+
+        myproc = base_objects.Process({'legs':myleglist,
+                                       'model':self.base_model})
+
+        myamp = diagram_generation.Amplitude(myproc)
+        matrix_element = helas_objects.HelasMatrixElement(myamp)
+        self.assertEqual(13, len(matrix_element.get('color_basis')))
+
 
 #===============================================================================
 # TestColorTripletModel