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

[Olivier Mattelaer]

The created param_card follow the SLHA convention:
- presence of QNUMBERS
- all particles have mass/width defined

=== modified file 'bin/mg5'
--- bin/mg5	2011-01-19 12:10:40 +0000
+++ bin/mg5	2011-03-21 08:51:32 +0000
@@ -43,7 +43,7 @@
  
     # charge history file
     try:
-        history_file = os.path.join(os.environ['HOME'], '.pythonhistory')
+        history_file = os.path.join(os.environ['HOME'], '.mg5history')
         readline.read_history_file(history_file)
     except:
         pass
@@ -123,6 +123,7 @@
         cmd_interface.MadGraphCmdShell(mgme_dir = options.mgme_dir).cmdloop()
 
 try: 
+    readline.set_history_length(100)
     readline.write_history_file(history_file)
 except:
     pass

=== modified file 'madgraph/VERSION'
--- madgraph/VERSION	2011-03-09 20:12:43 +0000
+++ madgraph/VERSION	2011-03-21 08:51:32 +0000
@@ -1,4 +1,4 @@
-version = 0.6.2.3
-date = 2011-03-09
+version = 0.7.0-beta
+date = 2011-03-18
 
 

=== modified file 'madgraph/__init__.py'
--- madgraph/__init__.py	2010-10-19 14:51:20 +0000
+++ madgraph/__init__.py	2011-03-21 08:51:32 +0000
@@ -16,6 +16,9 @@
     """Exception raised if an exception is find 
     Those Types of error will stop nicely in the cmd interface"""
 
+class InvalidCmd(MadGraph5Error):
+    """a class for the invalid syntax call"""
+
 import os
 
 #Look for basic file position MG5DIR and MG4DIR

=== modified file 'madgraph/core/base_objects.py'
--- madgraph/core/base_objects.py	2011-01-26 13:23:54 +0000
+++ madgraph/core/base_objects.py	2011-03-21 08:51:32 +0000
@@ -22,6 +22,7 @@
 import numbers
 import os
 import re
+import StringIO
 import madgraph.core.color_algebra as color
 from madgraph import MadGraph5Error, MG5DIR
 
@@ -87,12 +88,12 @@
 
         return self[name]
 
-    def set(self, name, value):
+    def set(self, name, value, force=False):
         """Set the value of the property name. First check if value
         is a valid value for the considered property. Return True if the
         value has been correctly set, False otherwise."""
 
-        if not __debug__:
+        if not __debug__ or force:
             self[name] = value
             return True
 
@@ -263,7 +264,7 @@
             if not isinstance(value, str):
                 raise self.PhysicsObjectError, \
                     "Line type %s is not a string" % repr(value)
-            if value not in ['dashed', 'straight', 'wavy', 'curly', 'double']:
+            if value not in ['dashed', 'straight', 'wavy', 'curly', 'double','swavy','scurly']:
                 raise self.PhysicsObjectError, \
                    "Line type %s is unknown" % value
 
@@ -665,6 +666,7 @@
         self['ref_dict_to0'] = {}
         self['ref_dict_to1'] = {}
         self['got_majoranas'] = None
+        self['conserved_charge'] = set()
 
     def filter(self, name, value):
         """Filter for model property values"""
@@ -674,43 +676,45 @@
                 raise self.PhysicsObjectError, \
                     "Object of type %s is not a string" % \
                                                             type(value)
-        if name == 'particles':
+        elif name == 'particles':
             if not isinstance(value, ParticleList):
                 raise self.PhysicsObjectError, \
                     "Object of type %s is not a ParticleList object" % \
                                                             type(value)
-        if name == 'interactions':
+        elif name == 'interactions':
             if not isinstance(value, InteractionList):
                 raise self.PhysicsObjectError, \
                     "Object of type %s is not a InteractionList object" % \
                                                             type(value)
-        if name == 'particle_dict':
-            if not isinstance(value, dict):
-                raise self.PhysicsObjectError, \
-                    "Object of type %s is not a dictionary" % \
-                                                        type(value)
-        if name == 'interaction_dict':
-            if not isinstance(value, dict):
-                raise self.PhysicsObjectError, \
-                    "Object of type %s is not a dictionary" % \
-                                                        type(value)
-
-        if name == 'ref_dict_to0':
-            if not isinstance(value, dict):
-                raise self.PhysicsObjectError, \
-                    "Object of type %s is not a dictionary" % \
-                                                        type(value)
-        if name == 'ref_dict_to1':
-            if not isinstance(value, dict):
-                raise self.PhysicsObjectError, \
-                    "Object of type %s is not a dictionary" % \
-                                                        type(value)
-
-        if name == 'got_majoranas':
+        elif name == 'particle_dict':
+            if not isinstance(value, dict):
+                raise self.PhysicsObjectError, \
+                    "Object of type %s is not a dictionary" % \
+                                                        type(value)
+        elif name == 'interaction_dict':
+            if not isinstance(value, dict):
+                raise self.PhysicsObjectError, \
+                    "Object of type %s is not a dictionary" % type(value)
+
+        elif name == 'ref_dict_to0':
+            if not isinstance(value, dict):
+                raise self.PhysicsObjectError, \
+                    "Object of type %s is not a dictionary" % type(value)
+                    
+        elif name == 'ref_dict_to1':
+            if not isinstance(value, dict):
+                raise self.PhysicsObjectError, \
+                    "Object of type %s is not a dictionary" % type(value)
+
+        elif name == 'got_majoranas':
             if not (isinstance(value, bool) or value == None):
                 raise self.PhysicsObjectError, \
-                    "Object of type %s is not a boolean" % \
-                                                        type(value)
+                    "Object of type %s is not a boolean" % type(value)
+
+        elif name == 'conserved_charge':
+            if not (isinstance(value, set)):
+                raise self.PhysicsObjectError, \
+                    "Object of type %s is not a set" % type(value)
 
         return True
 
@@ -867,45 +871,19 @@
                 
     def write_param_card(self):
         """Write out the param_card, and return as string."""
-
-        def write_param(param, lhablock):
-
-            lhacode=' '.join(['%3s' % key for key in param.lhacode])
-            if lhablock == 'DECAY':
-                return 'DECAY %s %e # %s' % (lhacode, param.value, param.name)
-            else:
-                return "  %s %e # %s" % (lhacode, param.value, param.name ) 
-
-        if not self.get('parameters'):
-            raise self.PhysicsObjectError,\
-                  "Attempt to write param_card from non-UFO model"
-
-        external_params = self.get('parameters')[('external',)]
-
-        # list all lhablock
-        all_lhablock = set([param.lhablock for param in external_params])
-        
-        # sort lhablock alphabeticaly
-        all_lhablock = sorted(list(all_lhablock))
-        # place DECAY blocks last
-        all_lhablock.remove('DECAY')
-        all_lhablock.append('DECAY')
-
-        ret_list = ["# SLHA param_card for %s written by MadGraph 5" % \
-                    self.get('name')]
-        
-        for lhablock in all_lhablock:
-            if lhablock != 'DECAY':
-                ret_list.append("BLOCK %s" % lhablock)
-            ret_list.extend(sorted([write_param(param, lhablock) \
-                                    for param in external_params if \
-                                    param.lhablock == lhablock]))
-        return "\n".join(ret_list) + "\n"
+        
+        import models.write_param_card as writter
+        out = StringIO.StringIO() # it's suppose to be written in a file
+        param = writter.ParamCardWriter(self)
+        param.define_output_file(out)
+        param.write_card()
+        return out.getvalue()
         
     @ staticmethod
     def load_default_name():
         """ load the default for name convention """
         
+        logger.info('Change particles name to pass to MG5 convention')    
         default = {}
         for line in open(os.path.join(MG5DIR, 'input', \
                                                  'particles_name_default.txt')):
@@ -1385,6 +1363,7 @@
         self['model'] = Model()
         # Optional number to identify the process
         self['id'] = 0
+        self['uid'] = 0 # should be a uniq id number
         # Required s-channels are given as a list of id lists. Only
         # diagrams with all s-channels in any of the lists are
         # allowed. This enables generating e.g. Z/gamma as s-channel
@@ -1412,10 +1391,10 @@
             if not isinstance(value, Model):
                 raise self.PhysicsObjectError, \
                         "%s is not a valid Model object" % str(value)
-        if name == 'id':
+        if name in ['id', 'uid']:
             if not isinstance(value, int):
                 raise self.PhysicsObjectError, \
-                    "Process id %s is not an integer" % repr(value)
+                    "Process %s %s is not an integer" % (name, repr(value))
 
         if name == 'required_s_channels':
             if not isinstance(value, list):
@@ -1637,7 +1616,7 @@
         # Remove last space
         return mystr[:-1]
 
-    def shell_string(self):
+    def shell_string(self, schannel=True, forbid=True, main=True):
         """Returns process as string with '~' -> 'x', '>' -> '_',
         '+' -> 'p' and '-' -> 'm', including process number,
         intermediate s-channels and forbidden particles"""
@@ -1655,7 +1634,7 @@
                 mystr = mystr + '_'
                 # Add required s-channels
                 if self['required_s_channels'] and \
-                       self['required_s_channels'][0]:
+                       self['required_s_channels'][0] and schannel:
                     mystr += "_or_".join(["".join([self['model'].\
                                        get('particle_dict')[req_id].get_name() \
                                                 for req_id in id_list]) \
@@ -1668,7 +1647,7 @@
             prevleg = leg
 
         # Check for forbidden particles
-        if self['forbidden_particles']:
+        if self['forbidden_particles'] and forbid:
             mystr = mystr + '_no_'
             for forb_id in self['forbidden_particles']:
                 forbpart = self['model'].get('particle_dict')[forb_id]
@@ -1684,7 +1663,19 @@
         mystr = mystr.replace(' ', '')
 
         for decay in self.get('decay_chains'):
-            mystr = mystr + "_" + decay.shell_string()
+            mystr = mystr + "_" + decay.shell_string(schannel,forbid, main=False)
+
+        # Too long name are problematic so restrict them to a maximal of 70 char
+        if len(mystr) > 64 and main:
+            if schannel and forbid:
+                return self.shell_string(True, False, False)+ '-%s' % self['uid']
+            elif schannel:
+                return self.shell_string(False, False, False)+'-%s' % self['uid']
+            else:
+                return mystr[:64]+'-%s' % self['uid']
+            
+            
+            
 
         return mystr
 

=== modified file 'madgraph/core/diagram_generation.py'
--- madgraph/core/diagram_generation.py	2010-12-19 03:57:07 +0000
+++ madgraph/core/diagram_generation.py	2011-03-21 08:51:32 +0000
@@ -24,7 +24,7 @@
 
 import madgraph.core.base_objects as base_objects
 
-from madgraph import MadGraph5Error
+from madgraph import MadGraph5Error, InvalidCmd
 logger = logging.getLogger('madgraph.diagram_generation')
 
 #===============================================================================
@@ -165,11 +165,11 @@
                   
 
         res = base_objects.DiagramList()
-
         # First check that the number of fermions is even
         if len(filter(lambda leg: model.get('particle_dict')[\
                         leg.get('id')].is_fermion(), legs)) % 2 == 1:
             self['diagrams'] = res
+            raise InvalidCmd, 'The number of fermion is odd' 
             return res
 
         # Then check same number of incoming and outgoing fermions (if
@@ -178,8 +178,30 @@
            len(filter(lambda leg: leg.is_incoming_fermion(model), legs)) != \
            len(filter(lambda leg: leg.is_outgoing_fermion(model), legs)):
             self['diagrams'] = res
+            raise InvalidCmd, 'The number of of incoming/outcoming fermions are different' 
             return res
+        
+        # Finally check that charge (conserve by all interactions) of the process
+        #is globally conserve for this process.
+        for charge in model.get('conserved_charge'):
+            total = 0
+            for leg in legs:
+                part = model.get('particle_dict')[leg.get('id')]
+                try:
+                    value = part.get(charge)
+                except AttributeError:
+                    value = 0
+                    
+                if (leg.get('id') != part['pdg_code']) != leg['state']:
+                    total -= value
+                else:
+                    total += value
 
+            if abs(total) > 1e-10:
+                self['diagrams'] = res
+                raise InvalidCmd, 'No %s conservation for this process ' % charge
+                return res
+                    
         logger.info("Trying %s " % process.nice_string().replace('Process', 'process'))
 
         # Give numbers to legs in process
@@ -623,7 +645,7 @@
         """Test if object obj is a valid Amplitude for the list."""
 
         return isinstance(obj, Amplitude)
-
+    
 #===============================================================================
 # DecayChainAmplitude
 #===============================================================================
@@ -920,15 +942,23 @@
                     continue
 
                 amplitude = Amplitude({"process": process})
-                if not amplitude.generate_diagrams():
-                    # Add process to failed_procs
+                try:
+                    result = amplitude.generate_diagrams()
+                except InvalidCmd as error:
                     failed_procs.append(tuple(sorted_legs))
-                if amplitude.get('diagrams'):
-                    amplitudes.append(amplitude)
+                else:
+                    if amplitude.get('diagrams'):
+                        amplitudes.append(amplitude)
+                    elif not result:
+                        failed_procs.append(tuple(sorted_legs))
+ 
 
         # Raise exception if there are no amplitudes for this process
         if not amplitudes:
-            raise MadGraph5Error, \
+            if len(failed_procs) == 1 and 'error' in locals():
+                raise error
+            else:
+                raise MadGraph5Error, \
             "No amplitudes generated from process %s. Please enter a valid process" % \
                   process_definition.nice_string()
         
@@ -1092,14 +1122,18 @@
                         continue
 
                     amplitude = Amplitude({"process": process})
-                    if not amplitude.generate_diagrams():
-                        # Add process to failed_procs
+                    try:
+                        amplitude.generate_diagrams()
+                    except InvalidCmd:
                         failed_procs.append(tuple(sorted_legs))
+                    else:
+                        if amplitude.get('diagrams'):
+                            # We found a valid amplitude. Return this order number
+                            logger.setLevel(oldloglevel)
+                            return {coupling: max_order_now}
+                        else:
+                            failed_procs.append(tuple(sorted_legs))
 
-                    if amplitude.get('diagrams'):
-                        # We found a valid amplitude. Return this order number
-                        logger.setLevel(oldloglevel)
-                        return {coupling: max_order_now}
             # No processes found, increase max_order_now
             max_order_now += 1
             logger.setLevel(oldloglevel)

=== modified file 'madgraph/core/drawing.py'
--- madgraph/core/drawing.py	2010-12-01 09:33:43 +0000
+++ madgraph/core/drawing.py	2011-03-21 08:51:32 +0000
@@ -1809,8 +1809,11 @@
         # Find the type line of the particle [straight, wavy, ...]
         line_type = line.get_info('line')
         # Call the routine associate to this type [self.draw_straight, ...]
-        getattr(self, 'draw_' + line_type)(line)
-
+        if hasattr(self, 'draw_' + line_type):
+            getattr(self, 'draw_' + line_type)(line)
+        else:
+            self.draw_straight(line)
+            
         # Finalize the line representation with adding the name of the particle
         name = line.get_name()
         self.associate_name(line, name)

=== modified file 'madgraph/interface/cmd_interface.py'
--- madgraph/interface/cmd_interface.py	2011-01-27 12:30:18 +0000
+++ madgraph/interface/cmd_interface.py	2011-03-21 08:51:32 +0000
@@ -17,25 +17,18 @@
 """
 
 import atexit
-import cmd
 import logging
 import optparse
 import os
 import pydoc
 import re
-import signal
 import subprocess
 import sys
 import traceback
 import time
 
-# Optional Library (not present on all platform)
-try:
-    import readline
-except:
-    readline = None
-    
-from madgraph import MG4DIR, MG5DIR, MadGraph5Error
+
+from madgraph import MG4DIR, MG5DIR, MadGraph5Error, InvalidCmd
 
 import madgraph.core.base_objects as base_objects
 import madgraph.core.diagram_generation as diagram_generation
@@ -52,8 +45,9 @@
 import madgraph.iolibs.misc as misc
 import madgraph.iolibs.save_load_object as save_load_object
 
+import madgraph.interface.extended_cmd as cmd
 import madgraph.interface.tutorial_text as tutorial_text
-
+import madgraph.interface.launch_ext_program as launch_ext
 import madgraph.various.process_checks as process_checks
 
 import models as ufomodels
@@ -69,13 +63,11 @@
 # CmdExtended
 #===============================================================================
 class CmdExtended(cmd.Cmd):
-    """Extension of the cmd.Cmd command line.
-    This extensions supports line breaking, history, comments,
-    internal call to cmdline,..."""
+    """Particularisation of the cmd command for MG5"""
 
     #suggested list of command
     next_possibility = {
-        'mg5_start': ['import model ModelName', 'import command PATH',
+        'start': ['import model ModelName', 'import command PATH',
                       'import proc_v4 PATH', 'tutorial'],
         'import model' : ['generate PROCESS','define MULTIPART PART1 PART2 ...', 
                                    'display particles', 'display interactions'],
@@ -83,20 +75,29 @@
                                                     'display multiparticles'],
         'generate': ['add process PROCESS','output [OUTPUT_TYPE] [PATH]','draw .'],
         'add process':['output [OUTPUT_TYPE] [PATH]', 'display processes'],
-        'output':['history PATH', 'exit'],
+        'output':['launch','history PATH', 'exit'],
         'display': ['generate PROCESS', 'add process PROCESS', 'output [OUTPUT_TYPE] [PATH]'],
         'draw': ['shell CMD'],
-        'export':['history PATH', 'exit'],
-        'import proc_v4' : ['exit'],
+        'import proc_v4' : ['launch','exit'],
         'tutorial': ['generate PROCESS', 'import model MODEL', 'help TOPIC']
     }
-        
+    
+    debug_output = 'MG5_debug'
+    error_debug = 'Please report this bug on https://bugs.launchpad.net/madgraph5\n'
+    error_debug += 'More information is found in \'%s\'.\n' 
+    error_debug += 'Please attach this file to your report.'
+    
+    keyboard_stop_msg = """stopping all operation
+            in order to quit mg5 please enter exit"""
+    
+    class InvalidCmd(MadGraph5Error):
+        pass
+    
     def __init__(self, *arg, **opt):
         """Init history and line continuation"""
         
         # If possible, build an info line with current version number 
         # and date, from the VERSION text file
-
         info = misc.get_pkg_info()
         info_line = ""
 
@@ -162,128 +163,7 @@
         "*                                                          *\n" + \
         "************************************************************")
 
-        self.log = True
-        self.history = []
-        self.save_line = ''
         cmd.Cmd.__init__(self, *arg, **opt)
-        self.__initpos = os.path.abspath(os.getcwd())
-        
-    def precmd(self, line):
-        """ A suite of additional function needed for in the cmd
-        this implement history, line breaking, comment treatment,...
-        """
-        
-        if not line:
-            return line
-        line = line.lstrip()
-
-        # Update the history of this suite of command,
-        # except for useless commands (empty history and help calls)
-        if line != "history" and \
-            not line.startswith('help') and \
-            not line.startswith('#*'):
-            self.history.append(line)
-
-        # Check if we are continuing a line:
-        if self.save_line:
-            line = self.save_line + line 
-            self.save_line = ''
-        
-        # Check if the line is complete
-        if line.endswith('\\'):
-            self.save_line = line[:-1]
-            return '' # do nothing   
-        
-        # Remove comment
-        if '#' in line:
-            line = line.split('#')[0]
-
-        # Deal with line splitting
-        if ';' in line and not (line.startswith('!') or line.startswith('shell')):
-            for subline in line.split(';'):
-                stop = self.onecmd(subline)
-                stop = self.postcmd(stop, subline)
-            return ''
-        
-        # execute the line command
-        return line
-
-    def nice_error_handling(self, error, line):
-        """ """ 
-        # Make sure that we are at the initial position
-        os.chdir(self.__initpos)
-        # Create the debug files
-        self.log = False
-        cmd.Cmd.onecmd(self, 'history MG5_debug')
-        debug_file = open('MG5_debug', 'a')
-        traceback.print_exc(file=debug_file)
-        # Create a nice error output
-        if self.history and line == self.history[-1]:
-            error_text = 'Command \"%s\" interrupted with error:\n' % line
-        elif self.history:
-            error_text = 'Command \"%s\" interrupted in sub-command:\n' %line
-            error_text += '\"%s\" with error:\n' % self.history[-1]
-        else:
-            error_text = ''
-        error_text += '%s : %s\n' % (error.__class__.__name__, str(error).replace('\n','\n\t'))
-        error_text += 'Please report this bug on https://bugs.launchpad.net/madgraph5\n'
-        error_text += 'More information is found in \'%s\'.\n' % \
-                        os.path.realpath("MG5_debug")
-        error_text += 'Please attach this file to your report.'
-        logger_stderr.critical(error_text)
-        #stop the execution if on a non interactive mode
-        if self.use_rawinput == False:
-            sys.exit('Exit on error')
-        return False
-
-    def nice_user_error(self, error, line):
-        # Make sure that we are at the initial position
-        os.chdir(self.__initpos)
-        if line == self.history[-1]:
-            error_text = 'Command \"%s\" interrupted with error:\n' % line
-        else:
-            error_text = 'Command \"%s\" interrupted in sub-command:\n' %line
-            error_text += '\"%s\" with error:\n' % self.history[-1] 
-        error_text += '%s : %s' % (error.__class__.__name__, str(error).replace('\n','\n\t'))
-        logger_stderr.error(error_text)
-        #stop the execution if on a non interactive mode
-        if self.use_rawinput == False:
-            sys.exit()
-        # Remove failed command from history
-        self.history.pop()
-        return False
-
-    
-    def onecmd(self, line):
-        """catch all error and stop properly command accordingly"""
-        
-        try:
-            cmd.Cmd.onecmd(self, line)
-        except MadGraph5Error as error:
-            if __debug__:
-                self.nice_error_handling(error, line)
-            else:
-                self.nice_user_error(error, line)
-        except Exception as error:
-            self.nice_error_handling(error, line)
-            
-    def exec_cmd(self, line):
-        """for third party call, call the line with pre and postfix treatment"""
-
-        logger.info(line)
-        line = self.precmd(line)
-        stop = cmd.Cmd.onecmd(self, line)
-        stop = self.postcmd(stop, line)
-        return stop      
-
-    def run_cmd(self, line):
-        """for third party call, call the line with pre and postfix treatment"""
-        
-        logger.info(line)
-        line = self.precmd(line)
-        stop = self.onecmd(line)
-        stop = self.postcmd(stop, line)
-        return stop 
     
     def postcmd(self,stop, line):
         """ finishing a command
@@ -313,96 +193,17 @@
                 logger_tuto.info(getattr(tutorial_text, args[0]).replace('\n','\n\t'))
             except:
                 pass
-            
-
-
-    def emptyline(self):
-        """If empty line, do nothing. Default is repeat previous command."""
-        pass
-    
-    def default(self, line):
-        """Default action if line is not recognized"""
-
-        # Faulty command
-        logger.warning("Command \"%s\" not recognized, please try again" % \
-                                                                line.split()[0])
-    # Quit
-    def do_quit(self, line):
-        sys.exit(1)
-        
-    do_exit = do_quit
-
-    # Aliases
-    do_EOF = do_quit
-    do_exit = do_quit
-
-    def do_help(self, line):
-        """ propose some usefull possible action """
-        
-        cmd.Cmd.do_help(self,line)
-        
-        # if not basic help -> simple call is enough    
-        if line:
-            return
-
-        if len(self.history) == 0:
-            last_action_2 = last_action = 'mg5_start'
-        else:
-            last_action_2 = last_action = 'none'
-        
-        pos = 0
-        authorize = self.next_possibility.keys() 
-        while last_action_2  not in authorize and last_action not in authorize:
-            pos += 1
-            if pos > len(self.history):
-                last_action_2 = last_action = 'mg5_start'
-                break
-            
-            args = self.history[-1 * pos].split()
-            last_action = args[0]
-            if len(args)>1: 
-                last_action_2 = '%s %s' % (last_action, args[1])
-            else: 
-                last_action_2 = 'none'
-        
-        print 'Contextual Help'
-        print '==============='
-        if last_action_2 in authorize:
-            options = self.next_possibility[last_action_2]
-        elif last_action in authorize:
-            options = self.next_possibility[last_action]
-        
-        text = 'The following command(s) may be useful in order to continue.\n'
-        for option in options:
-            text+='\t %s \n' % option      
-        print text
-
+        
+        return stop
 
     def timed_input(self, question, default, timeout=None):
         """ a question with a maximal time to answer take default otherwise"""
-        class TimeOutError(Exception):
-            """Class for run-time error"""
-            pass
-        def handle_alarm(signum, frame): 
-            raise TimeOutError
-        if timeout is None:
+        
+        if not timeout:
             timeout = self.timeout
-        signal.signal(signal.SIGALRM, handle_alarm)
-        signal.alarm(timeout)
-        if timeout:
-            question += '[%ss to answer] ' % timeout
-        try:
-            result = raw_input(question)
-        except TimeOutError:
-            print default
-            return default
-        finally:
-            signal.alarm(0)
-        return result
-
-
-
-
+        
+        return misc.timed_input(question, default, timeout) 
+    
 
 #===============================================================================
 # Helper function
@@ -474,6 +275,10 @@
         logger.info("   example display particles e+.")
         logger.info("   For \"checks\", can specify only to see failed checks.")
 
+    def help_launch(self):
+        """help for launch command"""
+        _launch_parser.print_help()
+
     def help_tutorial(self):
         logger.info("syntax: tutorial [" + "|".join(self._tutorial_opts) + "]")
         logger.info("-- start/stop the tutorial mode")
@@ -562,6 +367,13 @@
         logger.info("-- define a multiparticle")
         logger.info("   Example: define p = g u u~ c c~ d d~ s s~ b b~")
 
+    def help_restrict(self):
+        logger.info("syntax: restrict [model] param_card")
+        logger.info('   Suppress in the model all the interactions with zero')
+        logger.info('   couplings according to the param_card given in parameter.')
+        logger.info('   All zero parameter of the param_card are also suppress of')
+        logger.info('   the model.')
+        
     def help_history(self):
         logger.info("syntax: history [FILEPATH|clean|.] ")
         logger.info("   If FILEPATH is \'.\' and \'output\' is done,")
@@ -592,9 +404,6 @@
 class CheckValidForCmd(object):
     """ The Series of help routine for the MadGraphCmd"""
     
-    class InvalidCmd(MadGraph5Error):
-        """a class for the invalid syntax call"""
-    
     class RWError(MadGraph5Error):
         """a class for read/write errors"""
     
@@ -605,10 +414,10 @@
     
         if len(args) < 2:
             self.help_add()
-            raise self.InvalidCmd('\"add\" requires two arguments')
+            raise InvalidCmd('\"add\" requires two arguments')
         
         if args[0] != 'process':
-            raise self.InvalidCmd('\"add\" requires the argument \"process\"')
+            raise InvalidCmd('\"add\" requires the argument \"process\"')
 
         if not self._curr_model:
             raise MadGraph5Error, \
@@ -624,20 +433,21 @@
         
         if len(args) < 2:
             self.help_define()
-            raise self.InvalidCmd('\"define\" command requires at least two arguments')
+            raise InvalidCmd('\"define\" command requires at least two arguments')
 
         if args[1] == '=':
             del args[1]
             if len(args) < 2:
                 self.help_define()
-                raise self.InvalidCmd('\"define\" command requires at least one particles name after \"=\"')
+                raise InvalidCmd('\"define\" command requires at least one particles name after \"=\"')
         
         if '=' in args:
             self.help_define()
-            raise self.InvalidCmd('\"define\" command requires symbols \"=\" at the second position')
+            raise InvalidCmd('\"define\" command requires symbols \"=\" at the second position')
         
         if not self._curr_model:
-            raise self.InvalidCmd("No particle list currently active, please import a model first")
+            logger.info('No model currently active. Try with the Standard Model')
+            self.do_import('model sm')
 
         if self._curr_model['particles'].find_name(args[0]):
             raise MadGraph5Error("label %s is a particle name in this model\n\
@@ -650,15 +460,15 @@
             
         if len(args) < 1 or args[0] not in self._display_opts:
             self.help_display()
-            raise self.InvalidCmd
+            raise InvalidCmd
 
         if not self._curr_model:
-            raise self.InvalidCmd("No model currently active, please import a model!")
+            raise InvalidCmd("No model currently active, please import a model!")
 
         if args[0] in ['processes', 'diagrams'] and not self._curr_amps:
-            raise self.InvalidCmd("No process generated, please generate a process!")
+            raise InvalidCmd("No process generated, please generate a process!")
         if args[0] == 'checks' and not self._comparisons:
-            raise self.InvalidCmd("No check results to display.")
+            raise InvalidCmd("No check results to display.")
 
 
     def check_draw(self, args):
@@ -668,13 +478,13 @@
         
         if len(args) < 1:
             self.help_draw()
-            raise self.InvalidCmd('\"draw\" command requires a directory path')
+            raise InvalidCmd('\"draw\" command requires a directory path')
         
         if not self._curr_amps:
-            raise self.InvalidCmd("No process generated, please generate a process!")
+            raise InvalidCmd("No process generated, please generate a process!")
             
         if not os.path.isdir(args[0]):
-            raise self.InvalidCmd( "%s is not a valid directory for export file" % args[0])
+            raise InvalidCmd( "%s is not a valid directory for export file" % args[0])
             
     def check_check(self, args):
         """check the validity of args"""
@@ -684,12 +494,12 @@
             self.do_import('model sm')
 
         if self._model_v4_path:
-            raise self.InvalidCmd(\
+            raise InvalidCmd(\
                 "\"check\" not possible for v4 models")
 
         if len(args) < 2:
             self.help_check()
-            raise self.InvalidCmd("\"check\" requires an argument and a process.")
+            raise InvalidCmd("\"check\" requires an argument and a process.")
 
         param_card = None
         if os.path.isfile(args[1]):
@@ -697,7 +507,7 @@
 
         if args[0] not in self._check_opts:
             self.help_check()
-            raise self.InvalidCmd("\"check\" called with wrong argument")
+            raise InvalidCmd("\"check\" called with wrong argument")
         
         if any([',' in elem for elem in args]):
             raise MadGraph5Error('Decay chains not allowed in check')
@@ -715,7 +525,7 @@
 
         if len(args) < 1:
             self.help_generate()
-            raise self.InvalidCmd("\"generate\" requires a process.")
+            raise InvalidCmd("\"generate\" requires a process.")
 
         self.check_process_format(" ".join(args))
  
@@ -726,7 +536,7 @@
         
         #check balance of paranthesis
         if process.count('(') != process.count(')'):
-            raise self.InvalidCmd('Invalid Format, no balance between open and close parenthesis')
+            raise InvalidCmd('Invalid Format, no balance between open and close parenthesis')
         #remove parenthesis for fututre introspection
         process = process.replace('(',' ').replace(')',' ')
         
@@ -739,7 +549,7 @@
         
         # request that we have one or two > in the process
         if process.count('>') not in [1,2]:
-            raise self.InvalidCmd(
+            raise InvalidCmd(
                'wrong format for \"%s\" this part requires one or two symbols \'>\', %s found' 
                % (process, process.count('>')))
         
@@ -747,16 +557,16 @@
         particles_parts = process.split('>')
         for particles in particles_parts:
             if re.match(r'^\s*$', particles):
-                raise self.InvalidCmd(
+                raise InvalidCmd(
                 '\"%s\" is a wrong process format. Please try again' % process)  
         
         # '/' and '$' sould be used only after the process definition
         for particles in particles_parts[:-1]:
             if re.search('\D/', particles):
-                raise self.InvalidCmd(
+                raise InvalidCmd(
                 'wrong process format: restriction should be place after the final states')
             if re.search('\D\$', particles):
-                raise self.InvalidCmd(
+                raise InvalidCmd(
                 'wrong process format: restriction should be place after the final states')
         
     
@@ -766,62 +576,119 @@
         
         if len(args) > 1:
             self.help_history()
-            raise self.InvalidCmd('\"history\" command takes at most one argument')
+            raise InvalidCmd('\"history\" command takes at most one argument')
         
         if not len(args):
             return
         
         if args[0] =='.':
             if not self._export_dir:
-                raise self.InvalidCmd("No default directory is defined for \'.\' option")
+                raise InvalidCmd("No default directory is defined for \'.\' option")
         elif args[0] != 'clean':
                 dirpath = os.path.dirname(args[0])
                 if dirpath and not os.path.exists(dirpath) or \
                        os.path.isdir(args[0]):
-                    raise self.InvalidCmd("invalid path %s " % dirpath)
+                    raise InvalidCmd("invalid path %s " % dirpath)
     
     def check_import(self, args):
         """check the validity of line"""
+ 
+        if '-modelname' in args:
+            if args[-1] != '-modelname':
+                args.remove('-modelname')
+                args.append('-modelname') 
         
         if not args or args[0] not in self._import_formats:
             self.help_import()
-            raise self.InvalidCmd('wrong \"import\" format')
+            raise InvalidCmd('wrong \"import\" format')
         
         if args[0].startswith('model') and len(args) != 2:
             if not (len(args) == 3 and args[-1] == '-modelname'):
                 self.help_import()
-                raise self.InvalidCmd('incorrect number of arguments')
+                raise InvalidCmd('incorrect number of arguments')
         
         if args[0] == 'proc_v4' and len(args) != 2 and not self._export_dir:
             self.help_import()
-            raise self.InvalidCmd('PATH is mandatory in the current context\n' + \
+            raise InvalidCmd('PATH is mandatory in the current context\n' + \
                                   'Did you forget to run the \"output\" command')
                         
-        if '-modelname' in args:
-            if args[-1] != '-modelname':
-                args.remove('-modelname')
-                args.append('-modelname')
+                
+    def check_launch(self, args, options):
+        """check the validity of the line"""
+        # modify args in order to be MODE DIR
+        # mode being either standalone or madevent
+        if not( 0 <= int(options.cluster) <= 2):
+            return self.InvalidCmd, 'cluster mode should be between 0 and 2'
+        
+        if not args:
+            if self._done_export:
+                args.append(self._done_export[1])
+                args.append(self._done_export[0])
+                return
+            else:
+                self.help_launch()
+                raise self.InvalidCmd, \
+                       'Impossible to use default location: No output command runned'
+        
+        if len(args) != 1:
+            self.help_launch()
+            return self.InvalidCmd, 'Invalid Syntax: Too many argument'
+        
+        # search for a valid path
+        if os.path.sep in args[0] and os.path.isdir(args[0]):
+            path = args[0]
+        elif os.path.isdir(os.path.join(MG5DIR,args[0])):
+            path = os.path.join(MG5DIR,args[0])
+        elif  MG4DIR and os.path.isdir(os.path.join(MG4DIR,args[0])):
+            path = os.path.join(MG4DIR,args[0])
+        elif os.path.isdir(args[0]):
+            path = args[0]
+        else:    
+            raise self.InvalidCmd, '%s is not a valid directory' % args[0]
+            
+        mode = self.find_output_type(path)
+        args[0] = mode
+        args.append(path)
+        
+    
+    def find_output_type(self, path):
+        """ identify the type of output of a given directory:
+        valid output: madevent/standalone/standalone_cpp"""
+        
+        card_path = os.path.join(path,'Cards')
+        bin_path = os.path.join(path,'bin')
+        subproc_path = os.path.join(path,'SubProcesses')
+        
+        if not os.path.isdir(card_path) or not os.path.isdir(subproc_path):
+            raise self.InvalidCmd, '%s : Not a valid directory' % path
+
+        if not os.path.isdir(bin_path):
+            return 'standalone_cpp'
+        elif os.path.isfile(os.path.join(bin_path,'generate_events')):
+            return 'madevent'
+        else:
+            return 'standalone'
         
     def check_load(self, args):
         """ check the validity of the line"""
         
         if len(args) != 2 or args[0] not in self._save_opts:
             self.help_load()
-            raise self.InvalidCmd('wrong \"load\" format')
+            raise InvalidCmd('wrong \"load\" format')
             
         
     def check_save(self, args):
         """ check the validity of the line"""
         if len(args) != 2 or args[0] not in self._save_opts:
             self.help_save()
-            raise self.InvalidCmd('wrong \"save\" format')
+            raise InvalidCmd('wrong \"save\" format')
     
     def check_output(self, args):
         """ check the validity of the line"""
         
         if not self._curr_model:
             text = 'No model found. Please import a model first and then retry.'
-            raise self.InvalidCmd(text)
+            raise InvalidCmd(text)
 
         if args and args[0] in self._export_formats:
             self._export_format = args.pop(0)
@@ -833,7 +700,7 @@
             text += " output for " + args[0] + ", you have to use a UFO model.\n"
             text += " Those model can be imported with mg5> import model NAME."
             logger.warning(text)
-            raise self.InvalidCmd('')
+            raise InvalidCmd('')
 
         if args and args[0][0] != '-':
             # This is a path
@@ -859,8 +726,34 @@
 
         if not self._curr_amps and self._export_format != "pythia8_model":
             text = 'No processes generated. Please generate a process first.'
-            raise self.InvalidCmd(text)
-
+            raise InvalidCmd(text)
+
+    def check_restrict(self,args):
+        """ check the format: restrict [model] param_card.dat"""
+        
+        if len(args) > 2:
+            self.help_restrict()
+            raise self.InvalidCmd, 'Wrong restrict format'
+        
+        if len(args) == 2:
+            if  args[0] != "model":
+                self.help_restrict()
+                raise self.InvalidCmd, 'Wrong restrict format'
+            else:
+                del args[0]
+        
+        if self._model_v4_path:
+            raise self.InvalidCmd, 'Operation not possible with v4 model. ' + \
+                            'Please use a UFO model as a starting point'
+                            
+        if self._restrict_file:
+            raise MadGraph5Error, 'This model is already restricted to the ' + \
+                        'card %s. In order to always keep track ' % self._restrict_file + \
+                        'of model modifications. We forbids multiple restrictions files.'
+
+        if not os.path.isfile(args[0]):
+            raise self.InvalidCmd, 'path \"%s\" is not a file' % args[0]
+            
     def get_default_path(self):
         """Set self._export_dir to the default (\'auto\') path"""
         
@@ -889,7 +782,7 @@
                 self._export_dir = auto_path(i) 
                 break
         if not self._export_dir:
-            raise self.InvalidCmd('Can\'t use auto path,' + \
+            raise InvalidCmd('Can\'t use auto path,' + \
                                   'more than 500 dirs already')    
             
         
@@ -908,7 +801,12 @@
         syntax: draw FILEPATH [option=value]
         """
         raise self.WebRestriction('direct call to draw is forbidden on the web')
-      
+    
+    def check_check(self, args):
+        """ Not authorize for the Web"""
+        
+        raise self.WebRestriction('Check call is forbidden on the web')
+    
     def check_history(self, args):
         """check the validity of line
         No Path authorize for the Web"""
@@ -977,41 +875,15 @@
                             for f in list
                             if f.startswith(text)
                             ]
-        return completions
-
-    def path_completion(self, text, base_dir = None, only_dirs = False):
-        """Propose completions of text to compose a valid path"""
-
-        if base_dir is None:
-            base_dir = os.getcwd()
-
-        if only_dirs:
-            completion = [f
-                          for f in os.listdir(base_dir)
-                          if f.startswith(text) and \
-                          os.path.isdir(os.path.join(base_dir, f)) and \
-                          (not f.startswith('.') or text.startswith('.'))
-                          ]
-        else:
-            completion = [f
-                          for f in os.listdir(base_dir)
-                          if f.startswith(text) and \
-                          os.path.isfile(os.path.join(base_dir, f)) and \
-                          (not f.startswith('.') or text.startswith('.'))
-                          ]
-
-            completion = completion + \
-                         [f + os.path.sep
-                          for f in os.listdir(base_dir)
-                          if f.startswith(text) and \
-                          os.path.isdir(os.path.join(base_dir, f)) and \
-                          (not f.startswith('.') or text.startswith('.'))
-                          ]
-
-        completion += [f for f in ['.'+os.path.sep, '..'+os.path.sep] if \
-                       f.startswith(text)]
-
-        return completion
+        
+        def put_space(name): 
+            if name.endswith(' '): 
+                return name
+            else:
+                return '%s ' % name 
+            
+        return [put_space(name) for name in completions] 
+
 
     def model_completion(self, text, process):
         """ complete the line with model information """
@@ -1089,7 +961,21 @@
             couplings = [c + "=" for c in self._couplings] + ['@','$','/','>']
         return self.list_completion(text, self._particle_names + \
                                     self._multiparticles.keys() + couplings)
-        
+    
+    def complete_restrict(self, text, line, begidx, endidx):
+        "Complete the restrict command"
+    
+        args = split_arg(line[0:begidx])
+        
+        if len(args) == 1 and text == 'model'[:len(text)]:
+            return ['model ']
+        
+        # Directory continuation
+        return self.path_completion(text,
+                                        os.path.join('.',*[a for a in args if a.endswith(os.path.sep)]),
+                                        only_dirs = False)
+            
+          
     def complete_check(self, text, line, begidx, endidx):
         "Complete the add command"
 
@@ -1163,6 +1049,33 @@
                                 'non_propagating', '--']
             return self.list_completion(text, opt)
 
+    def complete_launch(self, text, line, begidx, endidx):
+        """ complete the launch command"""
+
+        args = split_arg(line[0:begidx])
+
+        # Directory continuation
+        if args[-1].endswith(os.path.sep):
+            return self.path_completion(text,
+                                        os.path.join('.',*[a for a in args if a.endswith(os.path.sep)]),
+                                        only_dirs = True)
+        # Format
+        if len(args) == 1:
+            complete = self.path_completion(text, '.', only_dirs = True)
+            if MG5DIR != os.path.realpath('.'):
+                complete += self.path_completion(text, MG5DIR, only_dirs = True, 
+                                                                 relative=False)
+            if MG4DIR and MG4DIR != os.path.realpath('.'):
+                complete += self.path_completion(text, MG4DIR, only_dirs = True,
+                                                                 relative=False)
+            return complete
+
+        #option
+        if len(args) >= 2:
+            opt = ['--cluster=', '--name=', '-f']
+            return self.list_completion(text, opt)
+
+
     def complete_load(self, text, line, begidx, endidx):
         "Complete the load command"
 
@@ -1259,22 +1172,48 @@
         "Complete the import command"
 
         args=split_arg(line[0:begidx])
-
+        
         # Format
         if len(args) == 1:
             return self.list_completion(text, self._import_formats)
 
         # Directory continuation
-        if args[-1].endswith(os.path.sep):
+        if os.path.sep in args[-1] + text:
             if args[1].startswith('model'):
-                return self.path_completion(text,
-                                    os.path.join('.',*[a for a in args if \
-                                                      a.endswith(os.path.sep)]),
-                                    only_dirs = True)
+                # Directory continuation
+                return self.path_completion(text, os.path.join('.',*[a for a in args \
+                                                   if a.endswith(os.path.sep)]),
+                                                only_dirs = True)
             else:
                 return self.path_completion(text,
                                     os.path.join('.',*[a for a in args if \
                                                       a.endswith(os.path.sep)]))
+
+        # restriction continuation (for UFO)
+        if args[1] == 'model' and ('-' in args[-1] + text):
+            # deal with - in 2.7 as in 2.6
+            if sys.version_info[1] == 7:
+                prefix = '-'.join([part for part in text.split('-')[:-1]])+'-'
+                args.append(prefix)
+                text = text.split('-')[-1]
+            #model name
+            path = args[-1][:-1] # remove the final - for the model name
+            # find the different possibilities
+            all_name = self.find_restrict_card(path, no_restrict=False)
+            all_name += self.find_restrict_card(path, no_restrict=False,
+                                        base_dir=os.path.join(MG5DIR,'models'))
+
+            # select the possibility according to the current line            
+            all_name = [name.split('-')[-1] for name in  all_name ]
+            all_name = [name+' ' for name in  all_name if name.startswith(text)
+                                                       and name.strip() != text]
+            # adapt output for python2.7 (due to different splitting)
+            if sys.version_info[1] == 7:
+                all_name = [prefix + name for name in  all_name ]
+                
+            if all_name:
+                return all_name                  
+               
         # Model directory name if directory is not given
         if len(split_arg(line[0:begidx])) == 2:
             if args[1] == 'model':
@@ -1287,22 +1226,64 @@
             else:
                 return []
                 
-            return [mod_name(name) for name in \
-                    self.path_completion(text,
+            model_list = [mod_name(name) for name in \
+                                            self.path_completion(text,
                                             os.path.join(MG5DIR,'models'),
                                             only_dirs = True) \
-                       if file_cond(name)]
-                                
+                                            if file_cond(name)]
+            
+            if args[1] == 'model_v4':
+                return model_list
+            else:
+                # need to update the  list with the possible restriction
+                all_name = []
+                for model_name in model_list:
+                    all_name += self.find_restrict_card(model_name, 
+                                        base_dir=os.path.join(MG5DIR,'models'))
+                return all_name                
 
         # Options
         if len(args) > 2 and args[1].startswith('model') and args[-1][0] != '-':
                 return ['-modelname']
+            
         if len(args) > 3 and args[1].startswith('model') and args[-1][0] == '-':
-                return ['modelname']
-
-  
-
-
+                if sys.version_info[1] == 6:
+                    return ['modelname']
+                else: 
+                    return ['-modelname']
+    
+    def find_restrict_card(self, model_name, base_dir='./', no_restrict=True):
+        """find the restriction file associate to a given model"""
+
+        # check if the model_name should be keeped as a possibility
+        if no_restrict:
+            output = [model_name]
+        else:
+            output = []
+        
+        # check that the model is a valid model
+        if not os.path.exists(os.path.join(base_dir, model_name, 'couplings.py')):
+            # not valid UFO model
+            return output
+        
+        if model_name.endswith(os.path.sep):
+            model_name = model_name[:-1]
+        
+        # look for _default and treat this case
+        if os.path.exists(os.path.join(base_dir, model_name, 'restrict_default.dat')):
+            output.append('%s-full' % model_name)
+        
+        # look for other restrict_file
+        for name in os.listdir(os.path.join(base_dir, model_name)):
+            if name.startswith('restrict_') and not name.endswith('default.dat') \
+                and name.endswith('.dat'):
+                tag = name[9:-4] #remove restrict and .dat
+                while model_name.endswith(os.path.sep):
+                    model_name = model_name[:-1]
+                output.append('%s-%s' % (model_name, tag))
+
+        # return
+        return output
     
 #===============================================================================
 # MadGraphCmd
@@ -1342,6 +1323,7 @@
     _export_format = 'madevent'
     _mgme_dir = MG4DIR
     _comparisons = None
+    _restrict_file = None
     
     def __init__(self, mgme_dir = '', *completekey, **stdin):
         """ add a tracker of the history """
@@ -1422,6 +1404,7 @@
                                                     amp.nice_string_processes()
                         logger.warning(warning)
 
+
                 # Reset _done_export, since we have new process
                 self._done_export = False
 
@@ -1488,7 +1471,7 @@
                 else:
                     particle = self._curr_model['particles'].find_name(arg)
                 if not particle:
-                    raise self.InvalidCmd, 'no particle %s in current model' % arg
+                    raise InvalidCmd, 'no particle %s in current model' % arg
 
                 print "Particle %s has the following properties:" % particle.get_name()
                 print str(particle)
@@ -1496,8 +1479,8 @@
         elif args[0] == 'interactions' and len(args) == 1:
             text = "Current model contains %i interactions\n" % \
                     len(self._curr_model['interactions'])
-            for inter in self._curr_model['interactions']:
-                text += str(inter['id']) + ':'
+            for i, inter in enumerate(self._curr_model['interactions']):
+                text += str(i+1) + ':'
                 for part in inter['particles']:
                     if part['is_part']:
                         text += part['name']
@@ -1512,7 +1495,7 @@
         elif args[0] == 'interactions':
             for arg in args[1:]:
                 if int(arg) > len(self._curr_model['interactions']):
-                    raise self.InvalidCmd, 'no interaction %s in current model' % arg
+                    raise InvalidCmd, 'no interaction %s in current model' % arg
                 if int(arg) == 0:
                     print 'Special interactions which identify two particles'
                 else:
@@ -1567,21 +1550,21 @@
                 ufomodel = ufomodels.load_model(self._curr_model.get('name'))
                 print eval('ufomodel.couplings.%s.nice_string()'%args[1])
             except:
-                raise self.InvalidCmd, 'no couplings %s in current model' % args[1]
+                raise InvalidCmd, 'no couplings %s in current model' % args[1]
         
         elif args[0] == 'lorentz':
             if self._model_v4_path:
                 print 'No lorentz information available in V4 model'
                 return
             elif len(args) == 1: 
-                raise self.InvalidCmd,\
+                raise InvalidCmd,\
                      'display lorentz require an argument: the name of the lorentz structure.'
                 return
             try:
                 ufomodel = ufomodels.load_model(self._curr_model.get('name'))
                 print eval('ufomodel.lorentz.%s.nice_string()'%args[1])
             except:
-                raise self.InvalidCmd, 'no lorentz %s in current model' % args[1]
+                raise InvalidCmd, 'no lorentz %s in current model' % args[1]
             
         elif args[0] == 'checks':
             comparisons = self._comparisons[0]
@@ -1804,7 +1787,7 @@
         if not line.count('>') in [1,2]:
             self.do_help('generate')
             print
-            raise self.InvalidCmd('Wrong use of \">\" special character.')
+            raise InvalidCmd('Wrong use of \">\" special character.')
         
 
         # Perform sanity modifications on the lines:
@@ -2112,20 +2095,16 @@
         
         if args[0].startswith('model'):
             self._model_v4_path = None
+            self._restrict_file = None
             if args[0].endswith('_v4'):
                 self._curr_model, self._model_v4_path = \
                                  import_v4.import_model(args[1], self._mgme_dir)
                 self._curr_fortran_model = \
                       helas_call_writers.FortranHelasCallWriter(\
-                                                             self._curr_model)
+                                                             self._curr_model)                
             else:
-                self._curr_model = import_ufo.import_model(args[1])
-                self._curr_fortran_model = \
-                      helas_call_writers.FortranUFOHelasCallWriter(\
-                                                             self._curr_model)
-                self._curr_cpp_model = \
-                      helas_call_writers.CPPUFOHelasCallWriter(\
-                                                             self._curr_model)
+                self.import_ufo_model(args[1])
+
             if '-modelname' not in args:
                 self._curr_model.pass_particles_name_in_mg_default()
 
@@ -2141,7 +2120,7 @@
             if not os.path.isfile(args[1]):
                 raise MadGraph5Error("Path %s is not a valid pathname" % args[1])
             else:
-                # Check the status of export and try to use file position is no
+                # Check the status of export and try to use file position if no
                 #self._export dir are define
                 self.check_for_export_dir(args[1])
                 # Execute the card
@@ -2166,7 +2145,16 @@
  
             #convert and excecute the card
             self.import_mg4_proc_card(proc_card)
-                                     
+    
+    def import_ufo_model(self, model_name):
+        """ import the UFO model """
+        
+        self._curr_model = import_ufo.import_model(model_name)
+        self._curr_fortran_model = \
+                helas_call_writers.FortranUFOHelasCallWriter(self._curr_model)
+        self._curr_cpp_model = \
+                helas_call_writers.CPPUFOHelasCallWriter(self._curr_model)
+                
     def process_model(self):
         """Set variables _particle_names and _couplings for tab
         completion, defined multiparticles"""
@@ -2220,7 +2208,7 @@
         self.timeout, old_time_out = 20, self.timeout
         
         # Read the lines of the file and execute them
-        for line in CmdFile(filepath):
+        for line in cmd.CmdFile(filepath):
             #remove pointless spaces and \n
             line = line.replace('\n', '').strip()
             # execute the line
@@ -2287,6 +2275,32 @@
             self._export_dir = os.path.sep.join(path_split[:-2])
                 
     
+    def do_launch(self, line):
+        """Ask for editing the parameter and then 
+        Execute the code (madevent/standalone/...)
+        """
+        
+        args = split_arg(line)
+        # check argument validity and normalise argument
+        (options, args) = _launch_parser.parse_args(args)
+        self.check_launch(args, options)
+        options = options.__dict__
+        # args is now MODE PATH
+        
+        if args[0].startswith('standalone'):
+            ext_program = launch_ext.SALauncher(args[1], self.timeout, **options)
+        elif args[0] == 'madevent':
+            ext_program = launch_ext.MELauncher(args[1], self.timeout, **options)            
+        else:
+            raise self.InvalidCmd , '%s cannot be run from MG5 interface' % args[0]
+        
+        
+        
+        ext_program.run()
+        
+        
+        
+    
     def do_load(self, line):
         """Load information from file"""
 
@@ -2366,13 +2380,13 @@
                 if save_load_object.save_to_file(args[1], self._curr_model):
                     logger.info('Saved model to file %s' % args[1])
             else:
-                raise self.InvalidCmd('No model to save!')
+                raise InvalidCmd('No model to save!')
         elif args[0] == 'processes':
             if self._curr_amps:
                 if save_load_object.save_to_file(args[1], self._curr_amps):
                     logger.info('Saved processes to file %s' % args[1])
             else:
-                raise self.InvalidCmd('No processes to save!')
+                raise InvalidCmd('No processes to save!')
             
     def do_output(self, line):
         """Initialize a new Template or reinitialize one"""
@@ -2403,9 +2417,11 @@
         if self._export_format == 'madevent':
             export_v4.copy_v4template(self._mgme_dir, self._export_dir,
                                       not noclean)
+            export_v4.cp_model_restriction(self._restrict_file, self._export_dir)
         elif self._export_format == 'standalone':
             export_v4.copy_v4standalone(self._mgme_dir, self._export_dir,
                                         not noclean)
+            export_v4.cp_model_restriction(self._restrict_file, self._export_dir)
         elif self._export_format == 'standalone_cpp':
             export_cpp.setup_cpp_standalone_dir(self._export_dir, self._curr_model)
         elif not os.path.isdir(self._export_dir):
@@ -2419,7 +2435,7 @@
         if nojpeg:
             options = '-nojpeg'
 
-        self.export(options)
+        self.export(options)        
 
     # Export a matrix element
     def export(self, line):
@@ -2473,10 +2489,14 @@
         path = self._export_dir
         if self._export_format in ['madevent', 'standalone', 'standalone_cpp']:
             path = os.path.join(path, 'SubProcesses')
-            
+        
         if self._export_format == 'madevent':
+            uid = 0 # assign a unique id number to all process
             for ime, me in \
                 enumerate(self._curr_matrix_elements.get('matrix_elements')):
+                for proc in me.get('processes'):
+                    uid += 1 # update the identification number
+                    proc.set('uid', uid) # assign it to the process
                 calls = calls + \
                         export_v4.generate_subprocess_directory_v4_madevent(\
                             me, self._curr_fortran_model, ime, path)
@@ -2607,6 +2627,21 @@
         if self._export_format == 'madevent':
             logger.info('Please see ' + self._export_dir + '/README')
             logger.info('for information about how to generate events from this process.')
+            logger.info('You can also use the launch command.')
+
+    def do_restrict(self, line):
+        """ from a param_card.dat remove all zero interactions 
+            and all zero external parameter."""
+        
+        args = split_arg(line)
+        # Check args validity
+        self.check_restrict(args)
+        
+        
+        self._curr_model = import_ufo.RestrictModel(self._curr_model)
+        self._curr_model.restrict_model(args[0])
+        self._restrict_file = args[0]
+        
 
     def do_help(self, line):
         """ propose some usefull possible action """
@@ -2661,18 +2696,6 @@
 
         self.prompt = 'mg5>'
         
-        if readline:
-            readline.parse_and_bind("tab: complete")
-
-        # initialize command history if HOME exists
-        if os.environ.has_key('HOME') and readline:
-            history_file = os.path.join(os.environ['HOME'], '.mg5history')
-            try:
-                readline.read_history_file(history_file)
-            except IOError:
-                pass
-            atexit.register(readline.write_history_file, history_file)
-
         # By default, load the UFO Standard Model
         logger.info("Loading default model: sm")
         self.do_import('model sm')
@@ -2688,38 +2711,7 @@
             logging.info("running shell command: " + line)
             subprocess.call(line, shell=True)
 
-#===============================================================================
-# 
-#===============================================================================
-class CmdFile(file):
-    """ a class for command input file -in order to debug cmd \n problem"""
-    
-    def __init__(self, name, opt='rU'):
-        
-        file.__init__(self, name, opt)
-        self.text = file.read(self)
-        self.close()
-        self.lines = self.text.split('\n')
-    
-    def readline(self, *arg, **opt):
-        """readline method treating correctly a line whithout \n at the end
-           (add it)
-        """
-        if self.lines:
-            line = self.lines.pop(0)
-        else:
-            return ''
-        
-        if line.endswith('\n'):
-            return line
-        else:
-            return line + '\n'
-    
-    def __next__(self):
-        return self.lines.__next__()    
-    def __iter__(self):
-        return self.lines.__iter__()
-  
+
 #===============================================================================
 # Command Parser
 #=============================================================================== 
@@ -2743,7 +2735,18 @@
 _draw_parser.add_option("", "--add_gap", default=0, type='float', \
                           help="set the x-distance between external particles")  
 
-    
+# LAUNCH PROGRAM
+_launch_usage = "launch [DIRPATH] [options]\n" + \
+         "-- execute the madevent/standalone output present in DIRPATH\n" + \
+         "   By default DIRPATH is the latest created directory \n" + \
+         "   Example: launch PROC_SM_1 --name=run2 \n"
+_launch_parser = optparse.OptionParser(usage=_launch_usage)
+_launch_parser.add_option("-f", "--force", default=False, action='store_true',
+                                help="Use the card present in the directory in order to launch the different program")
+_launch_parser.add_option("-n", "--name", default='', type='str',
+                                help="Provide a name to the run (for madevent run)")
+_launch_parser.add_option("-c", "--cluster", default='0', type='str',
+                                help="Choose the cluster mode (0: single machine, 1: qsub cluster, 2: multi-core) (for madevent run)")
     
     
 #===============================================================================
@@ -2752,6 +2755,7 @@
 
 if __name__ == '__main__':
     
+    print 'pass here'
     run_option = sys.argv
     if len(run_option) > 1:
         # The first argument of sys.argv is the name of the program

=== added file 'madgraph/interface/extended_cmd.py'
--- madgraph/interface/extended_cmd.py	1970-01-01 00:00:00 +0000
+++ madgraph/interface/extended_cmd.py	2011-03-21 08:51:32 +0000
@@ -0,0 +1,381 @@
+################################################################################
+#
+# Copyright (c) 2011 The MadGraph Development team and Contributors
+#
+# This file is a part of the MadGraph 5 project, an application which 
+# automatically generates Feynman diagrams and matrix elements for arbitrary
+# high-energy processes in the Standard Model and beyond.
+#
+# It is subject to the MadGraph license which should accompany this 
+# distribution.
+#
+# For more information, please visit: http://madgraph.phys.ucl.ac.be
+#
+################################################################################
+"""  A file containing different extension of the cmd basic python library"""
+
+
+import cmd
+import logging
+import os
+import signal
+import traceback
+logger = logging.getLogger('cmdprint') # for stdout
+logger_stderr = logging.getLogger('fatalerror') # for stderr
+
+#===============================================================================
+# CmdExtended
+#===============================================================================
+class Cmd(cmd.Cmd):
+    """Extension of the cmd.Cmd command line.
+    This extensions supports line breaking, history, comments,
+    internal call to cmdline, path completion,...
+    this class should be MG5 independent"""
+
+    #suggested list of command
+    next_possibility = {} # command : [list of suggested command]
+    
+    class InvalidCmd(Exception):
+        """expected error for wrong command"""
+        pass    
+ 
+        debug_output = 'debug'
+        error_debug = """Please report this bug to developers\n
+           More information is found in '%s'.\n
+           Please attach this file to your report."""
+           
+        keyboard_stop_msg = """stopping all current operation
+            in order to quit the program please enter exit"""
+ 
+    
+    def __init__(self, *arg, **opt):
+        """Init history and line continuation"""
+        
+        self.log = True
+        self.history = []
+        self.save_line = ''
+        cmd.Cmd.__init__(self, *arg, **opt)
+        self.__initpos = os.path.abspath(os.getcwd())
+        
+
+        
+        
+    def precmd(self, line):
+        """ A suite of additional function needed for in the cmd
+        this implement history, line breaking, comment treatment,...
+        """
+        
+        if not line:
+            return line
+        line = line.lstrip()
+
+        # Update the history of this suite of command,
+        # except for useless commands (empty history and help calls)
+        if line != "history" and \
+            not line.startswith('help') and \
+            not line.startswith('#*'):
+            self.history.append(line)
+
+        # Check if we are continuing a line:
+        if self.save_line:
+            line = self.save_line + line 
+            self.save_line = ''
+        
+        # Check if the line is complete
+        if line.endswith('\\'):
+            self.save_line = line[:-1]
+            return '' # do nothing   
+        
+        # Remove comment
+        if '#' in line:
+            line = line.split('#')[0]
+
+        # Deal with line splitting
+        if ';' in line and not (line.startswith('!') or line.startswith('shell')):
+            for subline in line.split(';'):
+                stop = self.onecmd(subline)
+                stop = self.postcmd(stop, subline)
+            return ''
+        
+        # execute the line command
+        return line
+
+    def nice_error_handling(self, error, line):
+        """ """ 
+        # Make sure that we are at the initial position
+        os.chdir(self.__initpos)
+        # Create the debug files
+        self.log = False
+        cmd.Cmd.onecmd(self, 'history %s' % self.debug_output)
+        debug_file = open(self.debug_output, 'a')
+        traceback.print_exc(file=debug_file)
+        # Create a nice error output
+        if self.history and line == self.history[-1]:
+            error_text = 'Command \"%s\" interrupted with error:\n' % line
+        elif self.history:
+            error_text = 'Command \"%s\" interrupted in sub-command:\n' %line
+            error_text += '\"%s\" with error:\n' % self.history[-1]
+        else:
+            error_text = ''
+        error_text += '%s : %s\n' % (error.__class__.__name__, 
+                                                str(error).replace('\n','\n\t'))
+        error_text += self.error_debug % self.debug_output
+        logger_stderr.critical(error_text)
+        #stop the execution if on a non interactive mode
+        if self.use_rawinput == False:
+            return True 
+        return False
+
+    def nice_user_error(self, error, line):
+        # Make sure that we are at the initial position
+        os.chdir(self.__initpos)
+        if line == self.history[-1]:
+            error_text = 'Command \"%s\" interrupted with error:\n' % line
+        else:
+            error_text = 'Command \"%s\" interrupted in sub-command:\n' %line
+            error_text += '\"%s\" with error:\n' % self.history[-1] 
+        error_text += '%s : %s' % (error.__class__.__name__, 
+                                                str(error).replace('\n','\n\t'))
+        logger_stderr.error(error_text)
+        #stop the execution if on a non interactive mode
+        if self.use_rawinput == False:
+            return True
+        # Remove failed command from history
+        self.history.pop()
+        return False
+
+    def onecmd(self, line):
+        """catch all error and stop properly command accordingly"""
+        
+        try:
+            return cmd.Cmd.onecmd(self, line)
+        except self.InvalidCmd as error:
+            if __debug__:
+                self.nice_error_handling(error, line)
+            else:
+                self.nice_user_error(error, line)
+        except Exception as error:
+            self.nice_error_handling(error, line)
+        except KeyboardInterrupt:
+            print self.keyboard_stop_msg
+            
+    def exec_cmd(self, line, errorhandling=False):
+        """for third party call, call the line with pre and postfix treatment
+        without global error handling """
+
+        logger.info(line)
+        line = self.precmd(line)
+        if errorhandling:
+            stop = self.onecmd(line)
+        else:
+            stop = cmd.Cmd.onecmd(self, line)
+        stop = self.postcmd(stop, line)
+        return stop      
+
+    def run_cmd(self, line):
+        """for third party call, call the line with pre and postfix treatment
+        with global error handling"""
+        
+        return self.exec_cmd(line, errorhandling=True)
+    
+    def emptyline(self):
+        """If empty line, do nothing. Default is repeat previous command."""
+        pass
+    
+    def default(self, line):
+        """Default action if line is not recognized"""
+
+        # Faulty command
+        logger.warning("Command \"%s\" not recognized, please try again" % \
+                                                                line.split()[0])
+    # Quit
+    def do_quit(self, line):
+        """ exit the mainloop() """
+        print
+        return True
+ 
+    # Aliases
+    do_EOF = do_quit
+    do_exit = do_quit
+
+    def do_help(self, line):
+        """ propose some usefull possible action """
+        
+        cmd.Cmd.do_help(self,line)
+        
+        # if not basic help -> simple call is enough    
+        if line:
+            return
+
+        if len(self.history) == 0:
+            last_action_2 = last_action = 'start'
+        else:
+            last_action_2 = last_action = 'none'
+        
+        pos = 0
+        authorize = self.next_possibility.keys() 
+        while last_action_2  not in authorize and last_action not in authorize:
+            pos += 1
+            if pos > len(self.history):
+                last_action_2 = last_action = 'start'
+                break
+            
+            args = self.history[-1 * pos].split()
+            last_action = args[0]
+            if len(args)>1: 
+                last_action_2 = '%s %s' % (last_action, args[1])
+            else: 
+                last_action_2 = 'none'
+        
+        print 'Contextual Help'
+        print '==============='
+        if last_action_2 in authorize:
+            options = self.next_possibility[last_action_2]
+        elif last_action in authorize:
+            options = self.next_possibility[last_action]
+        
+        text = 'The following command(s) may be useful in order to continue.\n'
+        for option in options:
+            text+='\t %s \n' % option      
+        print text
+
+    def list_completion(self, text, list):
+        """Propose completions of text in list"""
+        if not text:
+            completions = list
+        else:
+            completions = [ f
+                            for f in list
+                            if f.startswith(text)
+                            ]
+        return completions
+
+    @staticmethod
+    def path_completion(text, base_dir = None, only_dirs = False, 
+                                                                 relative=True):
+        """Propose completions of text to compose a valid path"""
+        
+        if base_dir is None:
+            base_dir = os.getcwd()
+            
+        prefix, text = os.path.split(text)
+        base_dir = os.path.join(base_dir, prefix)
+        if prefix:
+            prefix += os.path.sep
+        
+        
+
+        if only_dirs:
+            completion = [prefix + f
+                          for f in os.listdir(base_dir)
+                          if f.startswith(text) and \
+                          os.path.isdir(os.path.join(base_dir, f)) and \
+                          (not f.startswith('.') or text.startswith('.'))
+                          ]
+        else:
+            completion = [ prefix + f
+                          for f in os.listdir(base_dir)
+                          if f.startswith(text) and \
+                          os.path.isfile(os.path.join(base_dir, f)) and \
+                          (not f.startswith('.') or text.startswith('.'))
+                          ]
+
+            completion = completion + \
+                         [prefix + f + os.path.sep
+                          for f in os.listdir(base_dir)
+                          if f.startswith(text) and \
+                          os.path.isdir(os.path.join(base_dir, f)) and \
+                          (not f.startswith('.') or text.startswith('.'))
+                          ]
+
+        if relative:
+            completion += [prefix + f for f in ['.'+os.path.sep, '..'+os.path.sep] if \
+                       f.startswith(text) and not prefix.startswith('.')]
+
+        return completion
+
+
+
+
+
+#===============================================================================
+# Question in order to return a path with auto-completion
+#===============================================================================
+class OneLinePathCompletion(cmd.Cmd):
+    """ a class for answering a question with the path autocompletion"""
+
+    def preloop(self):
+        """Initializing before starting the main loop"""
+        self.prompt = ''
+        self.value = None
+
+    def __init__(self,  allow_arg=[],*arg, **opt):
+        self.allow_arg = allow_arg
+        self.history_header = ''
+        cmd.Cmd.__init__(self, *arg, **opt)
+
+    def completenames(self, text, *ignored):
+        signal.alarm(0) # avoid timer if any
+        try:
+            return Cmd.path_completion(text,'.', only_dirs = False)
+        except Exception, error:
+            print error
+            
+    def default(self, line):
+        """Default action if line is not recognized"""
+        self.value = line
+
+    def postcmd(self, stop, line):
+        
+        try:    
+            if self.value in self.allow_arg or os.path.isfile(self.value):
+                return True
+            else:
+                raise Exception
+        except Exception:
+            print """not valid argument. Valid argument are file path or value in (%s).""" \
+                          % ','.join(self.allow_arg)
+            print 'please retry'
+            return False
+            
+    def cmdloop(self, intro=None):
+        cmd.Cmd.cmdloop(self, intro)
+        return self.value
+    
+# a function helper
+def raw_path_input(input_text, allow_arg=[]):
+    print input_text
+    obj = OneLinePathCompletion(allow_arg=allow_arg)
+    return os.path.relpath(obj.cmdloop())
+    
+#===============================================================================
+# 
+#===============================================================================
+class CmdFile(file):
+    """ a class for command input file -in order to debug cmd \n problem"""
+    
+    def __init__(self, name, opt='rU'):
+        
+        file.__init__(self, name, opt)
+        self.text = file.read(self)
+        self.close()
+        self.lines = self.text.split('\n')
+    
+    def readline(self, *arg, **opt):
+        """readline method treating correctly a line whithout \n at the end
+           (add it)
+        """
+        if self.lines:
+            line = self.lines.pop(0)
+        else:
+            return ''
+        
+        if line.endswith('\n'):
+            return line
+        else:
+            return line + '\n'
+    
+    def __next__(self):
+        return self.lines.__next__()    
+    def __iter__(self):
+        return self.lines.__iter__()

=== added file 'madgraph/interface/launch_ext_program.py'
--- madgraph/interface/launch_ext_program.py	1970-01-01 00:00:00 +0000
+++ madgraph/interface/launch_ext_program.py	2011-03-21 08:51:32 +0000
@@ -0,0 +1,296 @@
+################################################################################
+#
+# Copyright (c) 2009 The MadGraph Development team and Contributors
+#
+# This file is a part of the MadGraph 5 project, an application which 
+# automatically generates Feynman diagrams and matrix elements for arbitrary
+# high-energy processes in the Standard Model and beyond.
+#
+# It is subject to the MadGraph license which should accompany this 
+# distribution.
+#
+# For more information, please visit: http://madgraph.phys.ucl.ac.be
+#
+################################################################################
+
+import logging
+import os
+import subprocess
+
+
+import madgraph.iolibs.files as files
+import madgraph.iolibs.misc as misc
+import madgraph.interface.extended_cmd as cmd
+from madgraph import MG4DIR, MG5DIR, MadGraph5Error
+from madgraph.iolibs.files import cp
+
+logger = logging.getLogger('cmdprint.ext_program')
+
+class ExtLauncher(object):
+    """ Generic Class for executing external program """
+    
+    program_dir = ''
+    executable = ''  # path from program_dir
+    
+    cards = [] # files can be modified (path from self.card_dir)
+    force = False
+    
+    def __init__(self, running_dir, card_dir='', timeout=None, **options):
+        """ initialize an object """
+        
+        self.running_dir = running_dir
+        self.card_dir = os.path.join(self.running_dir, card_dir)
+        self.timeout = timeout
+        self.found_editor()
+        
+        #include/overwrite options
+        for key,value in options.items():
+            setattr(self, key, value)
+            
+            
+    def found_editor(self):
+        """ found a (shell) program for editing file """
+        
+        # let first try to use the prefer editor (if EDITOR is define)
+        # if not define use the first define in a pre-define list
+            
+        if os.environ.has_key('EDITOR') and misc.which(os.environ['EDITOR']):
+            self.editor =  os.environ['EDITOR']
+            return
+        
+        logger.info('INFO: You can choose your prefered editor by defining the shell variable EDITOR')
+        
+        possibility = ['vi', 'emacs', 'vim', 'gedit', 'nano']
+        for editor in possibility:
+            if misc.which(editor):
+                self.editor = editor
+                return
+        
+        logger.warning(
+         'No valid editor detected. Please configure the shell variable EDITOR')
+            
+    def run(self):
+        """ execute the main code """
+
+        self.prepare_run()        
+        for card in self.cards:
+            self.treat_input_file(card, default = 'n')
+
+        self.launch_program()
+
+        
+    def prepare_run(self):
+        """ aditional way to prepare the run"""
+        pass
+    
+    def launch_program(self):
+        """launch the main program"""
+        subprocess.call([self.executable], cwd=self.running_dir)
+    
+    def edit_file(self, path):
+        """edit a file"""
+
+        path = os.path.realpath(path)
+        subprocess.call([self.editor, path], cwd=os.getcwd())
+        
+    def ask(self, question, default):
+        """ ask a question """
+        
+        if not self.force:
+            try:
+                out =  misc.timed_input(question, default, timeout=self.timeout,
+                                        noerror=False)
+            except misc.TimeOutError:
+                # avoid to always wait a given time for the next answer
+                self.force = True
+            else:
+                self.timeout=None # answer at least one question so wait...
+                return out
+        else:
+            return default
+        
+    def treat_input_file(self, filename, default=None, msg=''):
+        """ask to edit a file"""
+        
+        if msg == '' and filename == 'param_card.dat':
+            msg = \
+            """WARNING: If you edit this file don\'t forget to modify 
+            consistently the different parameters, especially 
+            the width of all particles.""" 
+        
+        fct = lambda q: cmd.raw_path_input(q, allow_arg=['y','n'])     
+                                    
+        if not self.force:
+            if msg:  print msg
+            question = 'Do you want to edit file: %(card)s? [y/n/path of the new %(card)s]' 
+            question = question % {'card':filename}
+            try:
+                ans =  misc.timed_input(question, default, timeout=self.timeout,
+                                        noerror=False, fct=fct)
+            except misc.TimeOutError:
+                # avoid to always wait a given time for the next answer
+                self.force = True
+            else:
+                self.timeout=None # answer at least one question so wait...
+        else:
+            ans = default
+ 
+        if ans == 'y':
+            path = os.path.join(self.card_dir, filename)
+            self.edit_file(path)
+        elif ans == 'n':
+            return
+        else:
+            path = os.path.join(self.card_dir, filename)
+            files.cp(ans, path)
+        
+        
+        
+                    
+class SALauncher(ExtLauncher):
+    """ A class to launch a simple Standalone test """
+    
+    cards = ['param_card.dat']
+
+    def __init__(self, running_dir, timeout, **options):
+        """ initialize the StandAlone Version"""
+        
+        ExtLauncher.__init__(self, running_dir, './Cards', timeout, **options)
+    
+    def launch_program(self):
+        """launch the main program"""
+        sub_path = os.path.join(self.running_dir, 'SubProcesses')
+        for path in os.listdir(sub_path):
+            if path.startswith('P') and \
+                                   os.path.isdir(os.path.join(sub_path, path)):
+                cur_path =  os.path.join(sub_path, path)
+                # make
+                subprocess.call(['make'], cwd=cur_path)
+                # check
+                subprocess.call(['./check'], cwd=cur_path)
+                
+        
+class MELauncher(ExtLauncher):
+    """A class to launch MadEvent run"""
+    
+    cards = ['param_card.dat', 'run_card.dat']
+
+    def __init__(self, running_dir, timeout, **option):
+        """ initialize the StandAlone Version"""
+        
+        
+        ExtLauncher.__init__(self, running_dir, './Cards', timeout, **option)
+        self.executable = os.path.join('.', 'bin','generate_events')
+
+        assert hasattr(self, 'cluster')
+        assert hasattr(self, 'name')
+         
+        # Check for pythia-pgs directory
+        if os.path.isdir(os.path.join(MG5DIR,'pythia-pgs')):
+            self.pythia = os.path.join(MG5DIR,'pythia-pgs')
+        elif MG4DIR and os.path.isdir(os.path.join(MG4DIR,'pythia-pgs')):
+            self.pythia = os.path.join(MG4DIR,'pythia-pgs')
+        else:
+            self.pythia = None
+            
+        # Check for DELPHES directory
+        if os.path.isdir(os.path.join(MG5DIR,'Delphes')):
+            self.delphes = os.path.join(MG5DIR,'Delphes')
+        elif MG4DIR and os.path.isdir(os.path.join(MG4DIR,'Delphes')):
+            self.delphes = os.path.join(MG4DIR,'Delphes')
+        else:
+            self.delphes = None
+        
+        
+        # Assign a valid run name if not put in options
+        if self.name == '':
+            for i in range(1000):
+                path = os.path.join(self.running_dir, 'Events','run_%02i_banner.txt' % i)
+                if not os.path.exists(path):
+                    self.name = 'run%s' % i
+                    break
+        
+        if self.name == '':
+            raise MadGraph5Error, 'too much run in this directory'
+    
+    
+    def copy_default_card(self, name):
+
+        dico = {'dir': self.card_dir, 'name': name }
+
+        if not os.path.exists('%(dir)s/%(name)s_card.dat' % dico):
+            cp('%(dir)s/%(name)s_card_default.dat' % dico,
+                '%(dir)s/%(name)s_card.dat' % dico)
+    
+          
+    def prepare_run(self):
+        """ ask for pythia-pgs/delphes run """
+        
+        # Check If we Need to run pythia 
+        if not self.pythia or self.force:
+            return
+        
+        answer = self.ask('Do you want to run pythia? [y/n]','n')
+        if answer == 'y':
+            self.copy_default_card('pythia')
+            self.cards.append('pythia_card.dat')
+        else:
+            path = os.path.join(self.card_dir, 'pythia_card.dat')
+            try: os.remove(path)
+            except OSError: pass
+            return # no Need to ask for PGS
+        
+        answer = self.ask('Do you want to run PGS? [y/n]','n')
+        if answer == 'y':
+            self.copy_default_card('pgs')
+            self.cards.append('pgs_card.dat')
+            return # No Need to ask for Delphes
+        else:
+            path = os.path.join(self.card_dir, 'pgs_card.dat')
+            try: os.remove(path)
+            except OSError: pass
+    
+        if not self.delphes:
+            return
+        
+        answer = self.ask('Do you want to run Delphes? [y/n]','n')
+        if answer == 'y':
+            self.copy_default_card('delphes')
+            self.cards.append('delphes_card.dat')
+        else:
+            path = os.path.join(self.card_dir, 'delphes_card.dat')
+            try: os.remove(path)
+            except OSError: pass        
+    
+    def launch_program(self):
+        """launch the main program"""
+        
+        mode = str(self.cluster)
+        if mode == "0":
+            subprocess.call([self.executable, mode, self.name], 
+                                                           cwd=self.running_dir)
+        elif mode == "1":
+            subprocess.call([self.executable, mode, self.name, self.name], 
+                                                           cwd=self.running_dir)
+        elif mode == "2":
+            nb_node = self.ask('How many core do you want to use?', '2')
+            subprocess.call([self.executable, mode, nb_node, self.name], 
+                                                           cwd=self.running_dir)
+        
+        # Display the cross-section to the screen
+        path = os.path.join(self.running_dir, 'SubProcesses', '%s_results.dat' 
+                                                                    % self.name) 
+        fsock = open(path)
+        line = fsock.readline()
+        cross, error = line.split()[0:2]
+        
+        logger.info('The total cross-section is %s +- %s pb' % (cross, error))
+        logger.info('more information in %s' 
+                                 % os.path.join(self.running_dir, 'index.html'))
+        
+        
+        
+        
+        
+                    
+            

=== modified file 'madgraph/interface/tutorial_text.py'
--- madgraph/interface/tutorial_text.py	2010-10-07 00:08:46 +0000
+++ madgraph/interface/tutorial_text.py	2011-03-21 08:51:32 +0000
@@ -26,6 +26,7 @@
 produce the output for MadEvent. In this part we will learn
 a) How to generate a process
 b) How to create output for MadEvent
+c) How to run the MadEvent output
 
 Let's start with the first point, how to generate a process:
 mg5>generate p p > t t~
@@ -77,9 +78,15 @@
 output = """
 If you are following the tutorial, a directory MY_FIRST_MG5_RUN has
 been created which can be used in order to run MadEvent exactly as if
-it was coming from MG4, see MY_FIRST_MG5_RUN/README.
-
-This step ends the tutorial of the basic commands of MG5. You can
+it was coming from MG4.
+
+Additionally to the MG4 command (see MY_FIRST_MG5_RUN/README), you can also 
+generate your events/compute the cross-section from this interface:
+Please Enter:
+mg5> launch MY_FIRST_MG5_RUN
+"""
+
+launch = """This step ends the tutorial of the basic commands of MG5. You can
 always use the help to see the options available for different
 commands. For example, if you want to know all valid output formats,
 you can enter
@@ -92,11 +99,11 @@
 
 But you can also continue the tutorial to learn some other useful
 commands:
-c) How to load a model
-d) How to define a multi-particle label 
-e) How to store a history of the commands in a session
-f) How to call shell commands from MG5
-g) How to draw the diagrams for your processes without generating
+d) How to load a model
+e) How to define a multi-particle label 
+f) How to store a history of the commands in a session
+g) How to call shell commands from MG5
+h) How to draw the diagrams for your processes without generating
    MadEvent output
 
 To import a model, write:

=== modified file 'madgraph/iolibs/drawing_eps.py'
--- madgraph/iolibs/drawing_eps.py	2010-12-09 15:12:32 +0000
+++ madgraph/iolibs/drawing_eps.py	2011-03-21 08:51:32 +0000
@@ -40,7 +40,7 @@
 from __future__ import division
 
 import os
-
+import math
 import madgraph.core.drawing as draw
 import madgraph.core.base_objects as base_objects
 
@@ -166,15 +166,15 @@
                          line.end.pos_x, line.end.pos_y, 'Fhiggs')
 
 
-    def draw_wavy(self, line):
+    def draw_wavy(self, line, opt=0, type='d'):
         """ADD the EPS code for this photon line."""
 
         #add the code in the correct format
         self.text += self.line_format(line.start.pos_x, line.start.pos_y,
-                         line.end.pos_x, line.end.pos_y, '0 Fphotond')
-
-
-    def draw_curly(self, line):
+                         line.end.pos_x, line.end.pos_y, '%d Fphoton%s' % (opt,type))
+
+
+    def draw_curly(self, line, type=''):
         """ADD the EPS code for this gluon line."""
 
         # Due to the asymmetry in the way to draw the gluon (everything is draw
@@ -185,13 +185,45 @@
                                 line.start.pos_y > line.end.pos_y):
             self.text += self.line_format(line.start.pos_x,
                         line.start.pos_y, line.end.pos_x,
-                        line.end.pos_y, '0 Fgluon')
+                        line.end.pos_y, '0 Fgluon%s' % type)
         else:
             self.text += self.line_format(line.end.pos_x,
                         line.end.pos_y, line.start.pos_x,
-                        line.start.pos_y, '0 Fgluon')
-
-
+                        line.start.pos_y, '0 Fgluon%s' % type)
+    
+    def draw_scurly(self, line):
+        """ADD the EPS code for this gluino line."""
+        self.draw_curly(line, type='r'
+                        )
+        self.draw_straight(line)
+        
+    def draw_swavy(self, line):
+        """ADD the EPS code for this neutralino line."""
+        self.draw_wavy(line, type='r')
+        self.draw_straight(line)  
+             
+    def draw_double(self, line, type='r'):
+        """ADD the EPS code for this neutralino line."""
+        
+        
+        length = math.sqrt((line.end.pos_y - line.start.pos_y)**2 + (line.end.pos_x - line.start.pos_x) **2)
+        c1 = (line.end.pos_x - line.start.pos_x)/length
+        c2 = (line.end.pos_y - line.start.pos_y)/length
+        
+        gap = 0.013
+        start2_x = line.start.pos_x + gap * c1  
+        start2_y = line.start.pos_y + gap * c2
+        stop1_x = line.end.pos_x - gap * c1
+        stop1_y = line.end.pos_y - gap * c2
+        
+        
+        self.text += self.line_format(line.start.pos_x, line.start.pos_y,
+                         stop1_x, stop1_y, '0 Fphoton%s' % (type))
+        #add the code in the correct format
+        self.text += self.line_format(start2_x, start2_y,
+                         line.end.pos_x, line.end.pos_y, '0 Fphoton%s' % (type))
+        
+        
     def put_diagram_number(self, number=0):
         """ADD the comment 'diagram [number]' just below the diagram."""
 

=== modified file 'madgraph/iolibs/export_cpp.py'
--- madgraph/iolibs/export_cpp.py	2010-11-24 20:10:07 +0000
+++ madgraph/iolibs/export_cpp.py	2011-03-21 08:51:32 +0000
@@ -1469,6 +1469,8 @@
                                                  p.depend))
             else:
                 for p in self.model['parameters'][key]:
+                    if p.name == 'ZERO':
+                        continue
                     self.params_indep.append(base_objects.ModelVariable(p.name,
                                                  self.p_to_cpp.parse(p.expr),
                                                  p.type,
@@ -1479,15 +1481,16 @@
             param = params_ext.pop(0)
             # Read value from the slha variable
             expression = ""
+            assert param.value.imag == 0
             if len(param.lhacode) == 1:
                 expression = "slha.get_block_entry(\"%s\", %d, %e);" % \
                              (param.lhablock.lower(), param.lhacode[0],
-                              param.value)
+                              param.value.real)
             elif len(param.lhacode) == 2:
                 expression = "indices[0] = %d;\nindices[1] = %d;\n" % \
                              (param.lhacode[0], param.lhacode[1])
                 expression += "%s=slha.get_block_entry(\"%s\", indices, %e);" \
-                              % (param.name, param.lhablock.lower(), param.value)
+                              % (param.name, param.lhablock.lower(), param.value.real)
             else:
                 raise MadGraph5Error("Only support for SLHA blocks with 1 or 2 indices")
             self.params_indep.insert(0,

=== modified file 'madgraph/iolibs/export_v4.py'
--- madgraph/iolibs/export_v4.py	2011-03-09 20:12:43 +0000
+++ madgraph/iolibs/export_v4.py	2011-03-21 08:51:32 +0000
@@ -32,8 +32,7 @@
 import madgraph.iolibs.ufo_expression_parsers as parsers
 
 import aloha.create_aloha as create_aloha
-
-import models.sm.write_param_card as write_param_card
+import models.write_param_card as param_writer
 from madgraph import MadGraph5Error, MG5DIR
 from madgraph.iolibs.files import cp, ln, mv
 _file_path = os.path.split(os.path.dirname(os.path.realpath(__file__)))[0] + '/'
@@ -125,7 +124,34 @@
         open(os.path.join(dir_path, 'SubProcesses', 'MGVersion.txt'), 'w').write(
                                                           MG_version['version'])
 
-      
+#===============================================================================
+# Copy the model restriction in the Model Directory
+#===============================================================================
+def cp_model_restriction(file_path, dir_path):
+    """Copy the model restriction in the Model Directory."""
+    
+    if not file_path:
+        return
+
+    assert os.path.isfile(file_path)
+    assert os.path.isdir(os.path.join(dir_path,'Source','MODEL'))
+    
+    output_path = os.path.join(dir_path,'Source','MODEL','restrict_model.dat')
+    
+    header="""#*********************************************************************
+#  THIS FILE WAS USED TO RESTRICT THE ORIGINAL MODEL
+#  PLEASE DON'T EDIT THIS FILE. HE IS IMPORTANT IN ORDER TO BE ABLE 
+#  TO REPRODUCE THE RESULT IN THE FUTURE.
+#*********************************************************************\n"""    
+    ff = open(output_path,'w')
+    ff.writelines(header)
+    ff.writelines(open(file_path).read())
+    ff.close()
+    
+                         
+                         
+    
+    
 #===============================================================================
 # Make the Helas and Model directories for Standalone directory
 #===============================================================================
@@ -1507,7 +1533,6 @@
        
         self.model = model
         self.model_name = model['name']
-        
         self.dir_path = output_path
         
         self.coups_dep = []    # (name, expression, type)
@@ -1588,14 +1613,15 @@
         index = self.params_dep.index('sqrt__aS')
         self.params_indep.insert(0, self.params_dep.pop(index))
         
-    def build(self, wanted_couplings = []):
+    def build(self, wanted_couplings = [], full=True):
         """modify the couplings to fit with MG4 convention and creates all the 
         different files"""
         
         self.pass_parameter_to_case_insensitive()
         self.refactorize(wanted_couplings)
         # write the files
-        self.write_all()
+        if full:
+            self.write_all()
 
     def open(self, name, comment='c', format='default'):
         """ Open the file name in the correct directory and with a valid
@@ -1641,6 +1667,7 @@
 
         # All the standard files
         self.copy_standard_file()
+
     ############################################################################
     ##  ROUTINE CREATING THE FILES  ############################################
     ############################################################################
@@ -1760,6 +1787,8 @@
         fsock.writelines("if(readlha) then\n")
         
         for param in self.params_indep:
+            if param.name == 'ZERO':
+                continue
             fsock.writelines("%s = %s\n" % (param.name,
                                             self.p_to_f.parse(param.expr)))
         
@@ -1885,7 +1914,8 @@
         fsock.writelines("""write(*,*)  ' Internal Params'
                             write(*,*)  ' ---------------------------------'
                             write(*,*)  ' '""")        
-        lines = [format(data.name) for data in self.params_indep]
+        lines = [format(data.name) for data in self.params_indep 
+                                                         if data.name != 'ZERO']
         fsock.writelines('\n'.join(lines))
         fsock.writelines("""write(*,*)  ' Internal Params evaluated point by point'
                             write(*,*)  ' ----------------------------------------'
@@ -1901,7 +1931,7 @@
     
         def format(parameter):
             """return the line for the ident_card corresponding to this parameter"""
-            colum = [parameter.lhablock] + \
+            colum = [parameter.lhablock.lower()] + \
                     [str(value) for value in parameter.lhacode] + \
                     [parameter.name]
             return ' '.join(colum)+'\n'
@@ -1918,11 +1948,10 @@
         def format_line(parameter):
             """return the line for the ident_card corresponding to this 
             parameter"""
-
             template = \
             """ call LHA_get_real(npara,param,value,'%(name)s',%(name)s,%(value)s)""" \
                 % {'name': parameter.name,
-                               'value': self.p_to_f.parse(str(parameter.value))}
+                   'value': self.p_to_f.parse(str(parameter.value.real))}
         
             return template        
     
@@ -1944,19 +1973,10 @@
 
     def create_param_card(self):
         """ create the param_card.dat """
-        
-        write_param_card.ParamCardWriter(
-                os.path.join(self.dir_path, 'param_card.dat'),
-                self.params_ext)
-
-#    def search_type(self, expr):
-#        """return the type associate to the expression"""
-#        
-#        for param in self.model.all_parameters:
-#            if param.name == expr:
-#                return param.type
-#        
-#        return CompactifyExpression.search_type(self, expr)
+
+        out_path = os.path.join(self.dir_path, 'param_card.dat')
+        param_writer.ParamCardWriter(self.model, out_path)
+
   
  
 

=== modified file 'madgraph/iolibs/files.py'
--- madgraph/iolibs/files.py	2010-11-06 18:34:32 +0000
+++ madgraph/iolibs/files.py	2011-03-21 08:51:32 +0000
@@ -88,7 +88,7 @@
 #===============================================================================
 # check piclke validity
 #===============================================================================
-def is_uptodate(picklefile, path_list=None, min_time=1289300000):
+def is_uptodate(picklefile, path_list=None, min_time=1300120445):
     """Check if the pickle files is uptodate compare to a list of files. 
     If no files are given, the pickle files is checked against it\' current 
     directory"""

=== modified file 'madgraph/iolibs/misc.py'
--- madgraph/iolibs/misc.py	2010-06-16 08:51:57 +0000
+++ madgraph/iolibs/misc.py	2011-03-21 08:51:32 +0000
@@ -17,6 +17,7 @@
 
 import os
 import re
+import signal
 import StringIO
 import time
 
@@ -85,6 +86,9 @@
     def is_exe(fpath):
         return os.path.exists(fpath) and os.access(fpath, os.X_OK)
 
+    if not program:
+        return None
+
     fpath, fname = os.path.split(program)
     if fpath:
         if is_exe(program):
@@ -96,3 +100,37 @@
                 return exe_file
 
     return None
+
+#===============================================================================
+# Ask a question with a maximum amount of time to answer
+#===============================================================================
+class TimeOutError(Exception):
+    """Class for run-time error"""
+         
+def timed_input(question, default, timeout=None, noerror=True, fct=None):
+    """ a question with a maximal time to answer take default otherwise"""
+    
+    def handle_alarm(signum, frame): 
+            raise TimeOutError
+        
+    signal.signal(signal.SIGALRM, handle_alarm)
+    
+    if fct is None:
+        fct = raw_input
+        
+    if timeout:
+        signal.alarm(timeout)
+        question += '[%ss to answer [%s]] ' % (timeout, str(default))    
+    try:
+        result = fct(question)
+    except TimeOutError:
+        if noerror:
+            print '\nuse %s' % default
+            return default
+        else:
+            signal.alarm(0)
+            raise
+    finally:
+        signal.alarm(0)
+    return result
+

=== modified file 'madgraph/iolibs/template_files/drawing_eps_header.inc'
--- madgraph/iolibs/template_files/drawing_eps_header.inc	2010-08-06 09:18:24 +0000
+++ madgraph/iolibs/template_files/drawing_eps_header.inc	2011-03-21 08:51:32 +0000
@@ -30,7 +30,7 @@
 def /Fx2 exch def /Fy1 exch def /Fx1 exch def
 Fx2 Fx1 sub Fy2 Fy1 sub Fbasis /Fttype 1
 Ftype 2 mod 2 mul abs sub def Ftype 0 ge
-{/Fddist Fdist Fr div def /Fn Fddist round
+{/Fddist Fdist Fr div 2 div def /Fn Fddist round
 def Ftype 1 gt {/Fn Fn Fddist Fn sub Fsign
 add def} if} {/Fn Fdist Fr div 2 div round 2
 mul def} ifelse Fx1 Fy1 moveto 0 1 Fnopoints
@@ -70,18 +70,25 @@
 
 /Fphoton
 {{ Fx1 Fx2 Fx1 sub Fi mul Fnopoints div Fn
+div Fxt Fi 360 mul Fnopoints div Frmod sin
+mul Fttype mul 2 div add add Fy1 Fy2 Fy1 sub
+Fi mul Fnopoints div Fn div Fyt Fi 360 mul
+Fnopoints div Frmod sin mul Fttype mul 2 div
+add add } Fstraight } def
+
+/Fphotonr
+{{ Fx1 Fx2 Fx1 sub Fi mul Fnopoints div Fn
 div Fxt Fi 180 mul Fnopoints div Frmod sin
-mul Fttype mul 2 div add add Fy1 Fy2 Fy1 sub
+mul Fttype mul 1 div add add Fy1 Fy2 Fy1 sub
 Fi mul Fnopoints div Fn div Fyt Fi 180 mul
-Fnopoints div Frmod sin mul Fttype mul 2 div
+Fnopoints div Frmod sin mul Fttype mul 1 div
 add add } Fstraight } def
 
-
 /Fphotond
 {{ Fx1 Fx2 Fx1 sub Fi mul Fnopoints div Fn
-div Fxt Fi 180 mul Fnopoints div Frmod sin
+div Fxt Fi 360 mul Fnopoints div Frmod sin
 mul Fttype mul 2 div add add Fy1 Fy2 Fy1 sub
-Fi mul Fnopoints div Fn div Fyt Fi 180 mul
+Fi mul Fnopoints div Fn div Fyt Fi 360 mul
 Fnopoints div Frmod sin mul Fttype mul 2 div
 add add } Fstraight Fx1
 Fx2 add 2 div Fy1 Fy2 add 2 div Farrow} def
@@ -108,6 +115,17 @@
 180 mul Fnopoints div sin mul add add add }
 Fstraight } def
 
+/Fgluonr
+{2 sub { Fx1 Fx2 Fx1 sub Fi mul Fnopoints div
+Fn div Fxt 0 Fi 120 mul Fnopoints div cos sub
+mul Fttype mul Fxl Fi 120 mul Fnopoints div
+sin mul add add add Fy1 Fy2 Fy1 sub Fi mul
+Fnopoints div Fn div Fyt 0 Fi 120 mul
+Fnopoints div cos sub mul Fttype mul Fyl Fi
+120 mul Fnopoints div sin mul add add add }
+Fstraight } def
+
+
 /Fgluonl
 {exch 2 sub exch { Fxc Fth cos Frr mul Fxt 1
 Fi 180 mul Fnopoints div cos sub mul Fttype

=== modified file 'madgraph/iolibs/ufo_expression_parsers.py'
--- madgraph/iolibs/ufo_expression_parsers.py	2010-10-26 20:13:20 +0000
+++ madgraph/iolibs/ufo_expression_parsers.py	2011-03-21 08:51:32 +0000
@@ -182,6 +182,7 @@
     def p_error(self, p):
         if p:
             logger.error("Syntax error at '%s'" % p.value)
+            raise MadGraph5Error
         else:
             logger.error("Syntax error at EOF")
         self.parsed_string = "Error"

=== modified file 'madgraph/various/process_checks.py'
--- madgraph/various/process_checks.py	2011-01-26 13:23:54 +0000
+++ madgraph/various/process_checks.py	2011-03-21 08:51:32 +0000
@@ -41,7 +41,7 @@
 
 import madgraph.various.rambo as rambo
 
-from madgraph import MG5DIR, MadGraph5Error
+from madgraph import MG5DIR, MadGraph5Error, InvalidCmd
 
 import models.model_reader as model_reader
 import aloha.template_files.wavefunctions as wavefunctions
@@ -191,7 +191,9 @@
                               multiprocess.get('is_decay_chain'),
                 'overall_orders': \
                               multiprocess.get('overall_orders')})
+            
             result = function(process, stored_quantities, *args)
+                        
             if result:
                 results.append(result)
                 
@@ -478,8 +480,14 @@
         newproc.set('legs', legs)
 
         # Generate the amplitude for this process
-        amplitude = diagram_generation.Amplitude(newproc)
-        if not amplitude.get('diagrams'):
+        try:
+            amplitude = diagram_generation.Amplitude(newproc)
+        except InvalidCmd:
+            result=False
+        else:
+            result = amplitude.get('diagrams')
+
+        if not result:
             # This process has no diagrams; go to next process
             logging.info("No diagrams for %s" % \
                          process.nice_string().replace('Process', 'process'))
@@ -905,7 +913,13 @@
     legs = process.get('legs')
     # Generate a process with these legs
     # Generate the amplitude for this process
-    amplitude = diagram_generation.Amplitude(process)
+    try:
+        amplitude = diagram_generation.Amplitude(process)
+    except InvalidCmd:
+        logging.info("No diagrams for %s" % \
+                         process.nice_string().replace('Process', 'process'))
+        return None
+    
     if not amplitude.get('diagrams'):
         # This process has no diagrams; go to next process
         logging.info("No diagrams for %s" % \

=== modified file 'models/import_ufo.py'
--- models/import_ufo.py	2011-03-01 12:51:49 +0000
+++ models/import_ufo.py	2011-03-21 08:51:32 +0000
@@ -28,17 +28,22 @@
 import madgraph.iolibs.save_load_object as save_load_object
 from madgraph.core.color_algebra import *
 
+
 import aloha.create_aloha as create_aloha
 
 import models as ufomodels
-
+import models.model_reader as model_reader
 logger = logging.getLogger('models.import_ufo')
-
-class UFOFormatError(Exception):
-    pass
-
-def import_model(model_name):
-    """ a practical and efficient way to import one of those models """
+logger_mod = logging.getLogger('madgraph.model')
+
+root_path = os.path.dirname(os.path.realpath( __file__ ))
+sys.path.append(root_path)
+
+class UFOImportError(MadGraph5Error):
+    """ a error class for wrong import of UFO model""" 
+
+def find_ufo_path(model_name):
+    """ find the path to a model """
 
     # Check for a valid directory
     if model_name.startswith('./') and os.path.isdir(model_name):
@@ -48,7 +53,58 @@
     elif os.path.isdir(model_name):
         model_path = model_name
     else:
-        raise MadGraph5Error("Path %s is not a valid pathname" % model_name)
+        raise UFOImportError("Path %s is not a valid pathname" % model_name)
+
+    return model_path
+
+def import_model(model_name):
+    """ a practical and efficient way to import a model"""
+        
+    # check if this is a valid path or if this include restriction file       
+    try:
+        model_path = find_ufo_path(model_name)
+    except UFOImportError:
+        if '-' not in model_name:
+            raise
+        split = model_name.split('-')
+        model_name = '-'.join([text for text in split[:-1]])
+        model_path = find_ufo_path(model_name)
+            
+        restrict_file = os.path.join(model_path, 'restrict_%s.dat'% split[-1])
+        #if restriction is full, then we by pass restriction (avoid default)
+        if split[-1] == 'full':
+            restrict_file = None
+    else:
+        # Check if by default we need some restrictions
+        if os.path.exists(os.path.join(model_path,'restrict_default.dat')):
+            restrict_file = os.path.join(model_path,'restrict_default.dat')
+        else:
+            restrict_file = None
+            
+    #import the FULL model
+    model = import_full_model(model_path) 
+
+    #restrict it if needed       
+    if restrict_file:
+        # but doing this in silence
+        old_level = logger_mod.level
+        if old_level < 30:
+            logger_mod.setLevel(30) # WARNING
+        
+        # Modify the mother class of the object in order to allow restriction
+        model = RestrictModel(model)
+        model.restrict_model(restrict_file)
+        
+        # put logger in normal mode
+        logger_mod.setLevel(old_level) 
+    
+    return model
+
+def import_full_model(model_path):
+    """ a practical and efficient way to import one of those models 
+        (no restriction file use)"""
+
+    assert model_path == find_ufo_path(model_path)
             
     # Check the validity of the model
     files_list_prov = ['couplings.py','lorentz.py','parameters.py',
@@ -57,7 +113,7 @@
     for filename in files_list_prov:
         filepath = os.path.join(model_path, filename)
         if not os.path.isfile(filepath):
-            raise MadGraph5Error,  "%s directory is not a valid UFO model: \n %s is missing" % \
+            raise UFOImportError,  "%s directory is not a valid UFO model: \n %s is missing" % \
                                                          (model_path, filename)
         files_list.append(filepath)
         
@@ -66,22 +122,20 @@
         try:
             model = save_load_object.load_from_file( \
                                           os.path.join(model_path, 'model.pkl'))
-        except Exception:
+        except Exception, error:
+            print error
             logger.info('failed to load model from pickle file. Try importing UFO from File')
         else:
             return model
 
     # Load basic information
-    ufo_model = ufomodels.load_model(model_name)
+    ufo_model = ufomodels.load_model(model_path)
     ufo2mg5_converter = UFOMG5Converter(ufo_model)
     model = ufo2mg5_converter.load_model()
-    model.set('name', os.path.split(model_name)[-1])
+    
+    if model_path[-1] == '/': model_path = model_path[:-1] #avoid empty name
+    model.set('name', os.path.split(model_path)[-1])
  
-    # Load Abstract Helas routine from Aloha
-    #abstract_model = create_aloha.AbstractALOHAModel(model_name)
-    #abstract_model.compute_all(save=False)
-    #model.set('lorentz', abstract_model)
-    
     # Load the Parameter/Coupling in a convinient format.
     parameters, couplings = OrganizeModelExpression(ufo_model).main()
     model.set('parameters', parameters)
@@ -91,6 +145,10 @@
     # save in a pickle files to fasten future usage
     save_load_object.save_to_file(os.path.join(model_path, 'model.pkl'), model) 
  
+    #if default and os.path.exists(os.path.join(model_path, 'restrict_default.dat')):
+    #    restrict_file = os.path.join(model_path, 'restrict_default.dat') 
+    #    model = import_ufo.RestrictModel(model)
+    #    model.restrict_model(restrict_file)
     return model
     
 
@@ -107,6 +165,7 @@
         self.model = base_objects.Model()
         self.model.set('particles', self.particles)
         self.model.set('interactions', self.interactions)
+        self.conservecharge = set(['charge'])
         
         self.ufomodel = model
         
@@ -132,6 +191,7 @@
         for interaction_info in self.ufomodel.all_vertices:
             self.add_interaction(interaction_info)
         
+        self.model.set('conserved_charge', self.conservecharge)
         return self.model
         
     
@@ -146,6 +206,7 @@
         # MG5 doesn't use ghost (use unitary gauges)
         if particle_info.spin < 0:
             return 
+        
         # MG5 doesn't use goldstone boson 
         if hasattr(particle_info, 'GoldstoneBoson'):
             if particle_info.GoldstoneBoson:
@@ -167,10 +228,14 @@
                         particle.set(key, value)
                 elif key == 'charge':
                     particle.set(key, float(value))
-                elif key in ['charge','width']:
+                elif key in ['mass','width']:
                     particle.set(key, str(value))
                 else:
                     particle.set(key, value)
+            elif key not in ('GhostNumber','selfconjugate','goldstoneboson'):
+                # add charge -we will check later if those are conserve 
+                self.conservecharge.add(key)
+                particle.set(key,value, force=True)
             
         assert(12 == nb_property) #basic check that all the information is there         
         
@@ -181,7 +246,6 @@
         # Add the particles to the list
         self.particles.append(particle)
 
-
     def add_interaction(self, interaction_info):
         """add an interaction in the MG5 model. interaction_info is the 
         UFO vertices information."""
@@ -224,6 +288,17 @@
                     # add to the interactions
                     self.interactions.append(interaction)
 
+        # check if this interaction conserve the charge defined
+        for charge in list(self.conservecharge): #duplicate to allow modification
+            total = 0
+            for part in interaction_info.particles:
+                try:
+                    total += getattr(part, charge)
+                except AttributeError:
+                    pass
+            if abs(total) > 1e-12:
+                logger.info('The model has interaction violating the charge: %s' % charge)
+                self.conservecharge.discard(charge)
     
     _pat_T = re.compile(r'T\((?P<first>\d*),(?P<second>\d*)\)')
     _pat_id = re.compile(r'Identity\((?P<first>\d*),(?P<second>\d*)\)')
@@ -501,7 +576,219 @@
         except:
             return 'complex'
             
-            
+class RestrictModel(model_reader.ModelReader):
+    """ A class for restricting a model for a given param_card.
+    Two rules apply:
+     - Vertex with zero couplings are throw away
+     - external parameter with zero input are changed into internal parameter."""
+     
+    def restrict_model(self, param_card):
+        """apply the model restriction following param_card"""
+        
+        # compute the value of all parameters
+        self.set_parameters_and_couplings(param_card)
+        
+        # deal with couplings
+        zero_couplings = self.detect_zero_couplings()
+        self.remove_couplings(zero_couplings)
+        
+        # deal with parameters
+        parameters = self.detect_special_parameters()
+        self.fix_parameter_values(*parameters)
+        
+        # deal with identical parameters
+        iden_parameters = self.detect_identical_parameters()
+        for iden_param in iden_parameters:
+            self.merge_identical_parameters(iden_param)
+        
+
+    def detect_zero_couplings(self):
+        """return a list with the name of all vanishing couplings"""
+        
+        zero_coupling = []
+        
+        for name, value in self['coupling_dict'].items():
+            if value == 0:
+                zero_coupling.append(name)
+        return zero_coupling
+    
+    
+    def detect_special_parameters(self):
+        """ return the list of (name of) parameter which are zero """
+        
+        null_parameters = []
+        one_parameters = []
+        for name, value in self['parameter_dict'].items():
+            if value == 0 and name != 'ZERO':
+                null_parameters.append(name)
+            elif value == 1:
+                one_parameters.append(name)
+                
+        return null_parameters, one_parameters
+    
+    def detect_identical_parameters(self):
+        """ return the list of tuple of name of parameter with the same 
+        input value """
+
+        # Extract external parameters
+        external_parameters = self['parameters'][('external',)]
+        
+        # define usefull variable to detect identical input
+        block_value_to_var={} #(lhablok, value): list_of_var
+        mult_param = set([])       # key of the previous dict with more than one
+                              #parameter.
+                              
+        #detect identical parameter and remove the duplicate parameter
+        for param in external_parameters[:]:
+            value = self['parameter_dict'][param.name]
+            if value == 0:
+                continue
+            key = (param.lhablock, value) 
+            if key in block_value_to_var:
+                block_value_to_var[key].append(param)
+                mult_param.add(key)
+                #remove the duplicate parameter
+                #external_parameters.remove(param)
+            else: 
+                block_value_to_var[key] = [param]        
+        
+        output=[]  
+        for key in mult_param:
+            output.append(block_value_to_var[key])
+            
+        return output
+            
+    def merge_identical_parameters(self, parameters):
+        """ merge the identical parameters given in argument """
+            
+        logger_mod.info('Parameters set to identical values: %s '% \
+                        ', '.join([obj.name for obj in parameters]))
+        
+        # Extract external parameters
+        external_parameters = self['parameters'][('external',)]
+        for i, obj in enumerate(parameters):
+            # Keeped intact the first one and store information
+            if i == 0:
+                obj.info = 'set of param :' + \
+                                     ', '.join([param.name for param in parameters])
+                expr = obj.name
+                continue
+            # delete the old parameters
+            external_parameters.remove(obj)
+            # replace by the new one pointing of the first obj of the class
+            new_param = base_objects.ModelVariable(obj.name, expr, 'real')
+            self['parameters'][()].insert(0, new_param)
+        
+    
+    def remove_couplings(self, zero_couplings):
+        """ remove the interactions associated to couplings"""
+        
+        # clean the interactions
+        for interaction in self['interactions'][:]:
+            modified = False
+            for key, coupling in interaction['couplings'].items()[:]:
+                if coupling in zero_couplings:
+                    modified = True
+                    del interaction['couplings'][key]
+                    
+            if modified: 
+                part_name = [part['name'] for part in interaction['particles']]
+                orders = ['%s=%s' % (order,value) 
+                               for order,value in interaction['orders'].items()]                    
+            if not interaction['couplings']:
+                logger_mod.info('remove interactions: %s at order: %s' % \
+                                (' '.join(part_name),', '.join(orders)))
+                self['interactions'].remove(interaction)
+            elif modified:
+                logger_mod.info('modify interactions: %s at order: %s' % \
+                                (' '.join(part_name),', '.join(orders)))                
+        #clean the coupling list:
+        for name, data in self['couplings'].items():
+            for coupling in data[:]:
+                if coupling.name in zero_couplings:
+                    data.remove(coupling)
+        
+    def fix_parameter_values(self, zero_parameters, one_parameters):
+        """ Remove all instance of the parameters in the model and replace it by 
+        zero when needed."""
+
+        special_parameters = zero_parameters + one_parameters
+        
+        # treat specific cases for masses and width
+        for particle in self['particles']:
+            if particle['mass'] in zero_parameters:
+                particle['mass'] = 'ZERO'
+            if particle['width'] in zero_parameters:
+                particle['width'] = 'ZERO'
+        for pdg, particle in self['particle_dict'].items():
+            if particle['mass'] in zero_parameters:
+                particle['mass'] = 'ZERO'
+            if particle['width'] in zero_parameters:
+                particle['width'] = 'ZERO'            
+        
+        # check if the parameters is still usefull:
+        re_str = '|'.join(special_parameters)
+        re_pat = re.compile(r'''\b(%s)\b''' % re_str)
+        used = set()
+        # check in coupling
+        for name, coupling_list in self['couplings'].items():
+            for coupling in coupling_list:
+                for use in  re_pat.findall(coupling.expr):
+                    used.add(use)  
+        
+        # simplify the regular expression
+        re_str = '|'.join([param for param in special_parameters 
+                                                          if param not in used])
+        re_pat = re.compile(r'''\b(%s)\b''' % re_str)
+        
+        param_info = {}
+        # check in parameters
+        for dep, param_list in self['parameters'].items():
+            for tag, parameter in enumerate(param_list):
+                # update information concerning zero/one parameters
+                if parameter.name in special_parameters:
+                    param_info[parameter.name]= {'dep': dep, 'tag': tag, 
+                                                               'obj': parameter}
+                    continue
+                                    
+                # Bypass all external parameter
+                if isinstance(parameter, base_objects.ParamCardVariable):
+                    continue
+
+                # check the presence of zero/one parameter
+                for use in  re_pat.findall(parameter.expr):
+                    used.add(use)
+
+        # modify the object for those which are still used
+        for param in used:
+            data = self['parameters'][param_info[param]['dep']]
+            data.remove(param_info[param]['obj'])
+            tag = param_info[param]['tag']
+            data = self['parameters'][()]
+            if param in zero_parameters:
+                data.insert(0, base_objects.ModelVariable(param, '0.0', 'real'))
+            else:
+                data.insert(0, base_objects.ModelVariable(param, '1.0', 'real'))
+                
+        # remove completely useless parameters
+        for param in special_parameters:
+            #by pass parameter still in use
+            if param in used:
+                logger_mod.info('fix parameter value: %s' % param)
+                continue 
+            logger_mod.info('remove parameters: %s' % param)
+            data = self['parameters'][param_info[param]['dep']]
+            data.remove(param_info[param]['obj'])
+                  
+                
+                
+        
+        
+        
+         
+      
+    
+      
         
         
     

=== modified file 'models/model_reader.py'
--- models/model_reader.py	2010-12-02 15:45:16 +0000
+++ models/model_reader.py	2011-03-21 08:51:32 +0000
@@ -104,7 +104,8 @@
                               value))
                         parameter_dict['decay'][(pid,)].value = complex(value)
                     except KeyError:
-                        logger.warning('No decay parameter found for %d' % pid)
+                        pass
+                        #logger.warning('No decay parameter found for %d' % pid)
                     continue
                 # Look for blocks
                 block_match = re_block.match(line)
@@ -121,9 +122,8 @@
                         exec("locals()[\'%s\'] = %s" % (parameter_dict[block][(i1,i2)].name,
                                           value))
                         parameter_dict[block][(i1,i2)].value = float(value)
-
                     except KeyError:
-                        logger.warning('No parameter found for block %s index %d %d' %\
+                            logger.warning('No parameter found for block %s index %d %d' %\
                                        (block, i1, i2))
                     continue
                 # Look for single indices
@@ -136,7 +136,8 @@
                                           value))
                         parameter_dict[block][(i1,)].value = complex(value)
                     except KeyError:
-                        logger.warning('No parameter found for block %s index %d' %\
+                        if block not in  ['qnumbers','mass']:
+                            logger.warning('No parameter found for block %s index %d' %\
                                        (block, i1))
                     continue
         else:
@@ -193,7 +194,7 @@
                                          derived_parameters]))
 
         # Add "zero"
-        self.get('parameter_dict')['ZERO'] = 0.
+        self.get('parameter_dict')['ZERO'] = complex(0.)
 
         self.set('coupling_dict', dict([(coup.name, coup.value) \
                                         for coup in couplings]))

=== modified file 'models/mssm/MSSM_UFO.log'
--- models/mssm/MSSM_UFO.log	2010-12-05 05:38:08 +0000
+++ models/mssm/MSSM_UFO.log	2011-03-21 08:51:32 +0000
@@ -1,13 +1,13 @@
-# This file was automatically created by FeynRules $Revision: 364 $
-# Mathematica version: 7.0 for Linux x86 (32-bit) (February 18, 2009)
-# Date: Sat 4 Dec 2010 11:43:16
+# This file was automatically created by FeynRules $Revision: 510 $
+# Mathematica version: 7.0 for Linux x86 (64-bit) (February 18, 2009)
+# Date: Wed 2 Mar 2011 15:02:57
 
 
 # 
 # This is the logfile for the model MSSM
 
-# Authors: Benjamin Fuks, Neil D. Christensen
-# Model version: 1.1.9
+# Authors: Benjamin Fuks
+# Model version: 1.3.0
 # Checking the Lagrangians
    * All Lagrangians are ok.
 #
@@ -33,34 +33,36 @@
 # Parameter definitions
 #
 
-   * All parameters are ok.
+   * Parameter gp has no value. Default value 1 assigned.
+   * Parameter gw has no value. Default value 1 assigned.
 
 
 # Vertices
-   * Calling FeynmanRules for 10 Lagrangians.
-   * Number of classes vertices: 268
-   * Number of flavored vertices: 931
+   * Calling FeynmanRules for 1 Lagrangians.
+   * Number of classes vertices: 174
+   * Number of flavored vertices: 822
    * Saved vertices in InterfaceRun[ 1 ].
    * Checked QNumber conservation.
       - Quantum number GhostNumber conserved in all vertices.
       - Quantum number Q conserved in all vertices.
+      - Quantum number Y conserved in all vertices.
    * particles.py written.
    * parameters.py written.
 #
 # Vertex definitions
 #
 
-   * 931 vertices written.
+   * 822 vertices written.
    * vertices.py written.
 #
 # Lorentz structure definitions
 #
 
-   * 19 lorentz structures written.
+   * 17 lorentz structures written.
    * lorentz.py written.
 #
 # Coupling definitions
 #
 
-   * 1063 couplings written.
+   * 945 couplings written.
    * couplings.py written.

=== modified file 'models/mssm/__init__.py'
--- models/mssm/__init__.py	2010-08-19 09:06:28 +0000
+++ models/mssm/__init__.py	2011-03-21 08:51:32 +0000
@@ -13,3 +13,8 @@
 all_lorentz = lorentz.all_lorentz
 all_parameters = parameters.all_parameters
 all_functions = function_library.all_functions
+
+
+__author__ = "Benjamin Fuks"
+__version__ = "1.3.0"
+__email__ = "[e-mail]"

=== modified file 'models/mssm/couplings.py'
--- models/mssm/couplings.py	2010-12-05 05:38:08 +0000
+++ models/mssm/couplings.py	2011-03-21 08:51:32 +0000
@@ -1,6 +1,6 @@
-# This file was automatically created by FeynRules $Revision: 364 $
-# Mathematica version: 7.0 for Linux x86 (32-bit) (February 18, 2009)
-# Date: Sat 4 Dec 2010 12:26:08
+# This file was automatically created by FeynRules $Revision: 510 $
+# Mathematica version: 7.0 for Linux x86 (64-bit) (February 18, 2009)
+# Date: Wed 2 Mar 2011 15:09:56
 
 
 from object_library import all_couplings, Coupling
@@ -46,4218 +46,3746 @@
                 order = {'QCD':2})
 
 GC_10 = Coupling(name = 'GC_10',
-                 value = '-(complex(0,1)*G*I1044)',
-                 order = {'QCD':1})
+                 value = '(2*ee**2*complex(0,1)*I1211)/9.',
+                 order = {'QED':2})
 
 GC_11 = Coupling(name = 'GC_11',
-                 value = 'complex(0,1)*G**2*I1044',
-                 order = {'QCD':2})
+                 value = '(-2*ee*complex(0,1)*G*I1211)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_12 = Coupling(name = 'GC_12',
-                 value = '-(complex(0,1)*G*I1055)',
-                 order = {'QCD':1})
+                 value = 'complex(0,1)*G**2*I1211',
+                 order = {'QCD':2})
 
 GC_13 = Coupling(name = 'GC_13',
-                 value = 'complex(0,1)*G**2*I1055',
-                 order = {'QCD':2})
+                 value = '(2*ee**2*complex(0,1)*I1222)/9.',
+                 order = {'QED':2})
 
 GC_14 = Coupling(name = 'GC_14',
-                 value = '(ee*complex(0,1)*I111)/3.',
-                 order = {'QED':1})
+                 value = '(-2*ee*complex(0,1)*G*I1222)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_15 = Coupling(name = 'GC_15',
-                 value = '(-2*ee*complex(0,1)*G*I111)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'complex(0,1)*G**2*I1222',
+                 order = {'QCD':2})
 
 GC_16 = Coupling(name = 'GC_16',
-                 value = '(-2*ee*complex(0,1)*I1111)/3.',
-                 order = {'QED':1})
+                 value = '(2*ee**2*complex(0,1)*I1233)/9. + (2*ee**2*complex(0,1)*I1333)/9.',
+                 order = {'QED':2})
 
 GC_17 = Coupling(name = 'GC_17',
-                 value = '-(complex(0,1)*G*I1111)',
-                 order = {'QCD':1})
+                 value = '(-2*ee*complex(0,1)*G*I1233)/3. - (2*ee*complex(0,1)*G*I1333)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_18 = Coupling(name = 'GC_18',
-                 value = '(-2*ee*complex(0,1)*I1122)/3.',
-                 order = {'QED':1})
+                 value = 'complex(0,1)*G**2*I1233 + complex(0,1)*G**2*I1333',
+                 order = {'QCD':2})
 
 GC_19 = Coupling(name = 'GC_19',
-                 value = '-(complex(0,1)*G*I1122)',
-                 order = {'QCD':1})
+                 value = '(2*ee**2*complex(0,1)*I1236)/9. + (2*ee**2*complex(0,1)*I1336)/9.',
+                 order = {'QED':2})
 
 GC_20 = Coupling(name = 'GC_20',
-                 value = '(ee*complex(0,1)*I122)/3.',
-                 order = {'QED':1})
+                 value = '(-2*ee*complex(0,1)*G*I1236)/3. - (2*ee*complex(0,1)*G*I1336)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_21 = Coupling(name = 'GC_21',
-                 value = '(-2*ee*complex(0,1)*G*I122)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'complex(0,1)*G**2*I1236 + complex(0,1)*G**2*I1336',
+                 order = {'QCD':2})
 
 GC_22 = Coupling(name = 'GC_22',
-                 value = '(-2*ee*complex(0,1)*I1133)/3. - (2*ee*complex(0,1)*I1233)/3.',
-                 order = {'QED':1})
+                 value = '(2*ee**2*complex(0,1)*I1344)/9.',
+                 order = {'QED':2})
 
 GC_23 = Coupling(name = 'GC_23',
-                 value = '-(complex(0,1)*G*I1133) - complex(0,1)*G*I1233',
-                 order = {'QCD':1})
+                 value = '(-2*ee*complex(0,1)*G*I1344)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_24 = Coupling(name = 'GC_24',
-                 value = '(-2*ee*complex(0,1)*I1136)/3. - (2*ee*complex(0,1)*I1236)/3.',
-                 order = {'QED':1})
+                 value = 'complex(0,1)*G**2*I1344',
+                 order = {'QCD':2})
 
 GC_25 = Coupling(name = 'GC_25',
-                 value = '-(complex(0,1)*G*I1136) - complex(0,1)*G*I1236',
-                 order = {'QCD':1})
+                 value = '(2*ee**2*complex(0,1)*I1355)/9.',
+                 order = {'QED':2})
 
 GC_26 = Coupling(name = 'GC_26',
-                 value = '(-2*ee*complex(0,1)*I1244)/3.',
-                 order = {'QED':1})
+                 value = '(-2*ee*complex(0,1)*G*I1355)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_27 = Coupling(name = 'GC_27',
-                 value = '-(complex(0,1)*G*I1244)',
-                 order = {'QCD':1})
+                 value = 'complex(0,1)*G**2*I1355',
+                 order = {'QCD':2})
 
 GC_28 = Coupling(name = 'GC_28',
-                 value = '(-2*ee*complex(0,1)*I1255)/3.',
-                 order = {'QED':1})
+                 value = '(2*ee**2*complex(0,1)*I1263)/9. + (2*ee**2*complex(0,1)*I1363)/9.',
+                 order = {'QED':2})
 
 GC_29 = Coupling(name = 'GC_29',
-                 value = '-(complex(0,1)*G*I1255)',
-                 order = {'QCD':1})
+                 value = '(-2*ee*complex(0,1)*G*I1263)/3. - (2*ee*complex(0,1)*G*I1363)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_30 = Coupling(name = 'GC_30',
-                 value = '(2*ee*complex(0,1)*I1163)/3. + (2*ee*complex(0,1)*I1263)/3.',
-                 order = {'QED':1})
+                 value = 'complex(0,1)*G**2*I1263 + complex(0,1)*G**2*I1363',
+                 order = {'QCD':2})
 
 GC_31 = Coupling(name = 'GC_31',
-                 value = 'complex(0,1)*G*I1163 + complex(0,1)*G*I1263',
-                 order = {'QCD':1})
+                 value = '(2*ee**2*complex(0,1)*I1266)/9. + (2*ee**2*complex(0,1)*I1366)/9.',
+                 order = {'QED':2})
 
 GC_32 = Coupling(name = 'GC_32',
-                 value = '(-2*ee*complex(0,1)*I1166)/3. - (2*ee*complex(0,1)*I1266)/3.',
-                 order = {'QED':1})
+                 value = '(-2*ee*complex(0,1)*G*I1266)/3. - (2*ee*complex(0,1)*G*I1366)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_33 = Coupling(name = 'GC_33',
-                 value = '-(complex(0,1)*G*I1166) - complex(0,1)*G*I1266',
-                 order = {'QCD':1})
+                 value = 'complex(0,1)*G**2*I1266 + complex(0,1)*G**2*I1366',
+                 order = {'QCD':2})
 
 GC_34 = Coupling(name = 'GC_34',
-                 value = '(8*ee**2*complex(0,1)*I1311)/9.',
-                 order = {'QED':2})
+                 value = 'ee*complex(0,1)*I2511',
+                 order = {'QED':1})
 
 GC_35 = Coupling(name = 'GC_35',
-                 value = '(4*ee*complex(0,1)*G*I1311)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'ee*complex(0,1)*I2522',
+                 order = {'QED':1})
 
 GC_36 = Coupling(name = 'GC_36',
-                 value = 'complex(0,1)*G**2*I1311',
-                 order = {'QCD':2})
+                 value = 'ee*complex(0,1)*I2533 + ee*complex(0,1)*I2633',
+                 order = {'QED':1})
 
 GC_37 = Coupling(name = 'GC_37',
-                 value = '(8*ee**2*complex(0,1)*I1322)/9.',
-                 order = {'QED':2})
+                 value = 'ee*complex(0,1)*I2536 + ee*complex(0,1)*I2636',
+                 order = {'QED':1})
 
 GC_38 = Coupling(name = 'GC_38',
-                 value = '(4*ee*complex(0,1)*G*I1322)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'ee*complex(0,1)*I2644',
+                 order = {'QED':1})
 
 GC_39 = Coupling(name = 'GC_39',
-                 value = 'complex(0,1)*G**2*I1322',
-                 order = {'QCD':2})
+                 value = 'ee*complex(0,1)*I2655',
+                 order = {'QED':1})
 
 GC_40 = Coupling(name = 'GC_40',
-                 value = '(8*ee**2*complex(0,1)*I1333)/9. + (8*ee**2*complex(0,1)*I1433)/9.',
-                 order = {'QED':2})
+                 value = '-(ee*complex(0,1)*I2563) - ee*complex(0,1)*I2663',
+                 order = {'QED':1})
 
 GC_41 = Coupling(name = 'GC_41',
-                 value = '(4*ee*complex(0,1)*G*I1333)/3. + (4*ee*complex(0,1)*G*I1433)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'ee*complex(0,1)*I2566 + ee*complex(0,1)*I2666',
+                 order = {'QED':1})
 
 GC_42 = Coupling(name = 'GC_42',
-                 value = 'complex(0,1)*G**2*I1333 + complex(0,1)*G**2*I1433',
-                 order = {'QCD':2})
+                 value = '2*ee**2*complex(0,1)*I3111',
+                 order = {'QED':2})
 
 GC_43 = Coupling(name = 'GC_43',
-                 value = '(8*ee**2*complex(0,1)*I1336)/9. + (8*ee**2*complex(0,1)*I1436)/9.',
+                 value = '2*ee**2*complex(0,1)*I3122',
                  order = {'QED':2})
 
 GC_44 = Coupling(name = 'GC_44',
-                 value = '(4*ee*complex(0,1)*G*I1336)/3. + (4*ee*complex(0,1)*G*I1436)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '2*ee**2*complex(0,1)*I3133 + 2*ee**2*complex(0,1)*I3233',
+                 order = {'QED':2})
 
 GC_45 = Coupling(name = 'GC_45',
-                 value = 'complex(0,1)*G**2*I1336 + complex(0,1)*G**2*I1436',
-                 order = {'QCD':2})
+                 value = '2*ee**2*complex(0,1)*I3136 + 2*ee**2*complex(0,1)*I3236',
+                 order = {'QED':2})
 
 GC_46 = Coupling(name = 'GC_46',
-                 value = '(8*ee**2*complex(0,1)*I1444)/9.',
+                 value = '2*ee**2*complex(0,1)*I3244',
                  order = {'QED':2})
 
 GC_47 = Coupling(name = 'GC_47',
-                 value = '(4*ee*complex(0,1)*G*I1444)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '2*ee**2*complex(0,1)*I3255',
+                 order = {'QED':2})
 
 GC_48 = Coupling(name = 'GC_48',
-                 value = 'complex(0,1)*G**2*I1444',
-                 order = {'QCD':2})
+                 value = '2*ee**2*complex(0,1)*I3163 + 2*ee**2*complex(0,1)*I3263',
+                 order = {'QED':2})
 
 GC_49 = Coupling(name = 'GC_49',
-                 value = '(8*ee**2*complex(0,1)*I1455)/9.',
+                 value = '2*ee**2*complex(0,1)*I3166 + 2*ee**2*complex(0,1)*I3266',
                  order = {'QED':2})
 
 GC_50 = Coupling(name = 'GC_50',
-                 value = '(4*ee*complex(0,1)*G*I1455)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '(ee*complex(0,1)*I511)/3.',
+                 order = {'QED':1})
 
 GC_51 = Coupling(name = 'GC_51',
-                 value = 'complex(0,1)*G**2*I1455',
-                 order = {'QCD':2})
+                 value = '-(complex(0,1)*G*I511)',
+                 order = {'QCD':1})
 
 GC_52 = Coupling(name = 'GC_52',
-                 value = '(8*ee**2*complex(0,1)*I1363)/9. + (8*ee**2*complex(0,1)*I1463)/9.',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5111)/3.',
+                 order = {'QED':1})
 
 GC_53 = Coupling(name = 'GC_53',
-                 value = '(4*ee*complex(0,1)*G*I1363)/3. + (4*ee*complex(0,1)*G*I1463)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '-(complex(0,1)*G*I5111)',
+                 order = {'QCD':1})
 
 GC_54 = Coupling(name = 'GC_54',
-                 value = 'complex(0,1)*G**2*I1363 + complex(0,1)*G**2*I1463',
-                 order = {'QCD':2})
+                 value = '(-2*ee*complex(0,1)*I5122)/3.',
+                 order = {'QED':1})
 
 GC_55 = Coupling(name = 'GC_55',
-                 value = '(8*ee**2*complex(0,1)*I1366)/9. + (8*ee**2*complex(0,1)*I1466)/9.',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5122)',
+                 order = {'QCD':1})
 
 GC_56 = Coupling(name = 'GC_56',
-                 value = '(4*ee*complex(0,1)*G*I1366)/3. + (4*ee*complex(0,1)*G*I1466)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '(ee*complex(0,1)*I522)/3.',
+                 order = {'QED':1})
 
 GC_57 = Coupling(name = 'GC_57',
-                 value = 'complex(0,1)*G**2*I1366 + complex(0,1)*G**2*I1466',
-                 order = {'QCD':2})
+                 value = '-(complex(0,1)*G*I522)',
+                 order = {'QCD':1})
 
 GC_58 = Coupling(name = 'GC_58',
-                 value = '-(ee**2*complex(0,1)*I1511)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5133)/3. - (2*ee*complex(0,1)*I5233)/3.',
+                 order = {'QED':1})
 
 GC_59 = Coupling(name = 'GC_59',
-                 value = '-(ee**2*complex(0,1)*I1522)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5133) - complex(0,1)*G*I5233',
+                 order = {'QCD':1})
 
 GC_60 = Coupling(name = 'GC_60',
-                 value = '-(ee**2*complex(0,1)*I1533)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5136)/3. - (2*ee*complex(0,1)*I5236)/3.',
+                 order = {'QED':1})
 
 GC_61 = Coupling(name = 'GC_61',
-                 value = '-(ee**2*complex(0,1)*I1536)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5136) - complex(0,1)*G*I5236',
+                 order = {'QCD':1})
 
 GC_62 = Coupling(name = 'GC_62',
-                 value = '-(ee**2*complex(0,1)*I1563)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5244)/3.',
+                 order = {'QED':1})
 
 GC_63 = Coupling(name = 'GC_63',
-                 value = '-(ee**2*complex(0,1)*I1566)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5244)',
+                 order = {'QCD':1})
 
 GC_64 = Coupling(name = 'GC_64',
-                 value = '-(ee**2*complex(0,1)*I1611)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5255)/3.',
+                 order = {'QED':1})
 
 GC_65 = Coupling(name = 'GC_65',
-                 value = '-(ee**2*complex(0,1)*I1622)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5255)',
+                 order = {'QCD':1})
 
 GC_66 = Coupling(name = 'GC_66',
-                 value = '-(ee**2*complex(0,1)*I1633)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(2*ee*complex(0,1)*I5163)/3. + (2*ee*complex(0,1)*I5263)/3.',
+                 order = {'QED':1})
 
 GC_67 = Coupling(name = 'GC_67',
-                 value = '-(ee**2*complex(0,1)*I1636)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = 'complex(0,1)*G*I5163 + complex(0,1)*G*I5263',
+                 order = {'QCD':1})
 
 GC_68 = Coupling(name = 'GC_68',
-                 value = '-(ee**2*complex(0,1)*I1663)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(-2*ee*complex(0,1)*I5166)/3. - (2*ee*complex(0,1)*I5266)/3.',
+                 order = {'QED':1})
 
 GC_69 = Coupling(name = 'GC_69',
-                 value = '-(ee**2*complex(0,1)*I1666)/(3.*cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I5166) - complex(0,1)*G*I5266',
+                 order = {'QCD':1})
 
 GC_70 = Coupling(name = 'GC_70',
-                 value = '(ee**2*complex(0,1)*I1711)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(ee*complex(0,1)*I533)/3. + (ee*complex(0,1)*I633)/3.',
+                 order = {'QED':1})
 
 GC_71 = Coupling(name = 'GC_71',
-                 value = '(ee**2*complex(0,1)*I1722)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I533) - complex(0,1)*G*I633',
+                 order = {'QCD':1})
 
 GC_72 = Coupling(name = 'GC_72',
-                 value = '(ee**2*complex(0,1)*I1733)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(ee*complex(0,1)*I536)/3. + (ee*complex(0,1)*I636)/3.',
+                 order = {'QED':1})
 
 GC_73 = Coupling(name = 'GC_73',
-                 value = '(ee**2*complex(0,1)*I1736)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I536) - complex(0,1)*G*I636',
+                 order = {'QCD':1})
 
 GC_74 = Coupling(name = 'GC_74',
-                 value = '(ee**2*complex(0,1)*I1811)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(ee*complex(0,1)*I644)/3.',
+                 order = {'QED':1})
 
 GC_75 = Coupling(name = 'GC_75',
-                 value = '(ee**2*complex(0,1)*I1822)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I644)',
+                 order = {'QCD':1})
 
 GC_76 = Coupling(name = 'GC_76',
-                 value = '(ee**2*complex(0,1)*I1833)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '(ee*complex(0,1)*I655)/3.',
+                 order = {'QED':1})
 
 GC_77 = Coupling(name = 'GC_77',
-                 value = '(ee**2*complex(0,1)*I1836)/(cw*cmath.sqrt(2))',
-                 order = {'QED':2})
+                 value = '-(complex(0,1)*G*I655)',
+                 order = {'QCD':1})
 
 GC_78 = Coupling(name = 'GC_78',
-                 value = '(ee*complex(0,1)*I133)/3. + (ee*complex(0,1)*I233)/3.',
+                 value = '-(ee*complex(0,1)*I563)/3. - (ee*complex(0,1)*I663)/3.',
                  order = {'QED':1})
 
 GC_79 = Coupling(name = 'GC_79',
-                 value = '(-2*ee*complex(0,1)*G*I133)/3. - (2*ee*complex(0,1)*G*I233)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'complex(0,1)*G*I563 + complex(0,1)*G*I663',
+                 order = {'QCD':1})
 
 GC_80 = Coupling(name = 'GC_80',
-                 value = '(ee*complex(0,1)*I136)/3. + (ee*complex(0,1)*I236)/3.',
+                 value = '(ee*complex(0,1)*I566)/3. + (ee*complex(0,1)*I666)/3.',
                  order = {'QED':1})
 
 GC_81 = Coupling(name = 'GC_81',
-                 value = '(-2*ee*complex(0,1)*G*I136)/3. - (2*ee*complex(0,1)*G*I236)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '-(complex(0,1)*G*I566) - complex(0,1)*G*I666',
+                 order = {'QCD':1})
 
 GC_82 = Coupling(name = 'GC_82',
-                 value = '(ee*complex(0,1)*I244)/3.',
-                 order = {'QED':1})
+                 value = '(8*ee**2*complex(0,1)*I7411)/9.',
+                 order = {'QED':2})
 
 GC_83 = Coupling(name = 'GC_83',
-                 value = '(-2*ee*complex(0,1)*G*I244)/3.',
+                 value = '(4*ee*complex(0,1)*G*I7411)/3.',
                  order = {'QCD':1,'QED':1})
 
 GC_84 = Coupling(name = 'GC_84',
-                 value = '(ee*complex(0,1)*I255)/3.',
-                 order = {'QED':1})
+                 value = 'complex(0,1)*G**2*I7411',
+                 order = {'QCD':2})
 
 GC_85 = Coupling(name = 'GC_85',
-                 value = '(-2*ee*complex(0,1)*G*I255)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = '(8*ee**2*complex(0,1)*I7422)/9.',
+                 order = {'QED':2})
 
 GC_86 = Coupling(name = 'GC_86',
-                 value = '-(ee*complex(0,1)*I163)/3. - (ee*complex(0,1)*I263)/3.',
-                 order = {'QED':1})
+                 value = '(4*ee*complex(0,1)*G*I7422)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_87 = Coupling(name = 'GC_87',
-                 value = '(-2*ee*complex(0,1)*G*I163)/3. - (2*ee*complex(0,1)*G*I263)/3.',
-                 order = {'QCD':1,'QED':1})
+                 value = 'complex(0,1)*G**2*I7422',
+                 order = {'QCD':2})
 
 GC_88 = Coupling(name = 'GC_88',
-                 value = '(ee*complex(0,1)*I166)/3. + (ee*complex(0,1)*I266)/3.',
-                 order = {'QED':1})
+                 value = '(8*ee**2*complex(0,1)*I7433)/9. + (8*ee**2*complex(0,1)*I7533)/9.',
+                 order = {'QED':2})
 
 GC_89 = Coupling(name = 'GC_89',
-                 value = '(-2*ee*complex(0,1)*G*I166)/3. - (2*ee*complex(0,1)*G*I266)/3.',
+                 value = '(4*ee*complex(0,1)*G*I7433)/3. + (4*ee*complex(0,1)*G*I7533)/3.',
                  order = {'QCD':1,'QED':1})
 
 GC_90 = Coupling(name = 'GC_90',
-                 value = '(2*ee**2*complex(0,1)*I311)/9.',
-                 order = {'QED':2})
+                 value = 'complex(0,1)*G**2*I7433 + complex(0,1)*G**2*I7533',
+                 order = {'QCD':2})
 
 GC_91 = Coupling(name = 'GC_91',
-                 value = '(2*ee**2*complex(0,1)*I322)/9.',
+                 value = '(8*ee**2*complex(0,1)*I7436)/9. + (8*ee**2*complex(0,1)*I7536)/9.',
                  order = {'QED':2})
 
 GC_92 = Coupling(name = 'GC_92',
-                 value = '(2*ee**2*complex(0,1)*I333)/9. + (2*ee**2*complex(0,1)*I433)/9.',
-                 order = {'QED':2})
+                 value = '(4*ee*complex(0,1)*G*I7436)/3. + (4*ee*complex(0,1)*G*I7536)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_93 = Coupling(name = 'GC_93',
-                 value = '(2*ee**2*complex(0,1)*I336)/9. + (2*ee**2*complex(0,1)*I436)/9.',
-                 order = {'QED':2})
+                 value = 'complex(0,1)*G**2*I7436 + complex(0,1)*G**2*I7536',
+                 order = {'QCD':2})
 
 GC_94 = Coupling(name = 'GC_94',
-                 value = '(2*ee**2*complex(0,1)*I444)/9.',
+                 value = '(8*ee**2*complex(0,1)*I7544)/9.',
                  order = {'QED':2})
 
 GC_95 = Coupling(name = 'GC_95',
-                 value = '(2*ee**2*complex(0,1)*I455)/9.',
-                 order = {'QED':2})
+                 value = '(4*ee*complex(0,1)*G*I7544)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_96 = Coupling(name = 'GC_96',
-                 value = '(2*ee**2*complex(0,1)*I363)/9. + (2*ee**2*complex(0,1)*I463)/9.',
-                 order = {'QED':2})
+                 value = 'complex(0,1)*G**2*I7544',
+                 order = {'QCD':2})
 
 GC_97 = Coupling(name = 'GC_97',
-                 value = '(2*ee**2*complex(0,1)*I366)/9. + (2*ee**2*complex(0,1)*I466)/9.',
+                 value = '(8*ee**2*complex(0,1)*I7555)/9.',
                  order = {'QED':2})
 
 GC_98 = Coupling(name = 'GC_98',
-                 value = 'ee*complex(0,1)*I511',
-                 order = {'QED':1})
+                 value = '(4*ee*complex(0,1)*G*I7555)/3.',
+                 order = {'QCD':1,'QED':1})
 
 GC_99 = Coupling(name = 'GC_99',
-                 value = 'ee*complex(0,1)*I522',
-                 order = {'QED':1})
+                 value = 'complex(0,1)*G**2*I7555',
+                 order = {'QCD':2})
 
 GC_100 = Coupling(name = 'GC_100',
-                  value = 'ee*complex(0,1)*I533 + ee*complex(0,1)*I633',
-                  order = {'QED':1})
+                  value = '(8*ee**2*complex(0,1)*I7463)/9. + (8*ee**2*complex(0,1)*I7563)/9.',
+                  order = {'QED':2})
 
 GC_101 = Coupling(name = 'GC_101',
-                  value = 'ee*complex(0,1)*I536 + ee*complex(0,1)*I636',
-                  order = {'QED':1})
+                  value = '(4*ee*complex(0,1)*G*I7463)/3. + (4*ee*complex(0,1)*G*I7563)/3.',
+                  order = {'QCD':1,'QED':1})
 
 GC_102 = Coupling(name = 'GC_102',
-                  value = 'ee*complex(0,1)*I644',
-                  order = {'QED':1})
+                  value = 'complex(0,1)*G**2*I7463 + complex(0,1)*G**2*I7563',
+                  order = {'QCD':2})
 
 GC_103 = Coupling(name = 'GC_103',
-                  value = 'ee*complex(0,1)*I655',
-                  order = {'QED':1})
+                  value = '(8*ee**2*complex(0,1)*I7466)/9. + (8*ee**2*complex(0,1)*I7566)/9.',
+                  order = {'QED':2})
 
 GC_104 = Coupling(name = 'GC_104',
-                  value = '-(ee*complex(0,1)*I563) - ee*complex(0,1)*I663',
-                  order = {'QED':1})
+                  value = '(4*ee*complex(0,1)*G*I7466)/3. + (4*ee*complex(0,1)*G*I7566)/3.',
+                  order = {'QCD':1,'QED':1})
 
 GC_105 = Coupling(name = 'GC_105',
-                  value = 'ee*complex(0,1)*I566 + ee*complex(0,1)*I666',
-                  order = {'QED':1})
+                  value = 'complex(0,1)*G**2*I7466 + complex(0,1)*G**2*I7566',
+                  order = {'QCD':2})
 
 GC_106 = Coupling(name = 'GC_106',
-                  value = '2*ee**2*complex(0,1)*I711',
-                  order = {'QED':2})
+                  value = '-(complex(0,1)*G*Rd11*cmath.sqrt(2))',
+                  order = {'QCD':1})
 
 GC_107 = Coupling(name = 'GC_107',
-                  value = '2*ee**2*complex(0,1)*I722',
-                  order = {'QED':2})
+                  value = '-(complex(0,1)*G*Rd22*cmath.sqrt(2))',
+                  order = {'QCD':1})
 
 GC_108 = Coupling(name = 'GC_108',
-                  value = '2*ee**2*complex(0,1)*I733 + 2*ee**2*complex(0,1)*I833',
-                  order = {'QED':2})
+                  value = '-(complex(0,1)*G*Rd33*cmath.sqrt(2))',
+                  order = {'QCD':1})
 
 GC_109 = Coupling(name = 'GC_109',
-                  value = '2*ee**2*complex(0,1)*I736 + 2*ee**2*complex(0,1)*I836',
-                  order = {'QED':2})
+                  value = 'complex(0,1)*G*Rd36*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_110 = Coupling(name = 'GC_110',
-                  value = '2*ee**2*complex(0,1)*I844',
-                  order = {'QED':2})
+                  value = 'complex(0,1)*G*Rd44*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_111 = Coupling(name = 'GC_111',
-                  value = '2*ee**2*complex(0,1)*I855',
-                  order = {'QED':2})
+                  value = 'complex(0,1)*G*Rd55*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_112 = Coupling(name = 'GC_112',
-                  value = '2*ee**2*complex(0,1)*I763 + 2*ee**2*complex(0,1)*I863',
-                  order = {'QED':2})
+                  value = '-(complex(0,1)*G*Rd63*cmath.sqrt(2))',
+                  order = {'QCD':1})
 
 GC_113 = Coupling(name = 'GC_113',
-                  value = '2*ee**2*complex(0,1)*I766 + 2*ee**2*complex(0,1)*I866',
-                  order = {'QED':2})
+                  value = 'complex(0,1)*G*Rd66*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_114 = Coupling(name = 'GC_114',
-                  value = '-(complex(0,1)*G*I911)',
+                  value = '-(complex(0,1)*G*Ru11*cmath.sqrt(2))',
                   order = {'QCD':1})
 
 GC_115 = Coupling(name = 'GC_115',
-                  value = 'complex(0,1)*G**2*I911',
-                  order = {'QCD':2})
+                  value = '-(complex(0,1)*G*Ru22*cmath.sqrt(2))',
+                  order = {'QCD':1})
 
 GC_116 = Coupling(name = 'GC_116',
-                  value = '-(complex(0,1)*G*I922)',
+                  value = '-(complex(0,1)*G*Ru33*cmath.sqrt(2))',
                   order = {'QCD':1})
 
 GC_117 = Coupling(name = 'GC_117',
-                  value = 'complex(0,1)*G**2*I922',
-                  order = {'QCD':2})
+                  value = 'complex(0,1)*G*Ru36*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_118 = Coupling(name = 'GC_118',
-                  value = '-(complex(0,1)*G*I1033) - complex(0,1)*G*I933',
+                  value = 'complex(0,1)*G*Ru44*cmath.sqrt(2)',
                   order = {'QCD':1})
 
 GC_119 = Coupling(name = 'GC_119',
-                  value = 'complex(0,1)*G**2*I1033 + complex(0,1)*G**2*I933',
-                  order = {'QCD':2})
+                  value = 'complex(0,1)*G*Ru55*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_120 = Coupling(name = 'GC_120',
-                  value = '-(complex(0,1)*G*I1036) - complex(0,1)*G*I936',
+                  value = '-(complex(0,1)*G*Ru63*cmath.sqrt(2))',
                   order = {'QCD':1})
 
 GC_121 = Coupling(name = 'GC_121',
-                  value = 'complex(0,1)*G**2*I1036 + complex(0,1)*G**2*I936',
-                  order = {'QCD':2})
+                  value = 'complex(0,1)*G*Ru66*cmath.sqrt(2)',
+                  order = {'QCD':1})
 
 GC_122 = Coupling(name = 'GC_122',
-                  value = 'complex(0,1)*G*I1063 + complex(0,1)*G*I963',
-                  order = {'QCD':1})
+                  value = '-(ee**2*complex(0,1)) + (ee**2*complex(0,1))/sw**2',
+                  order = {'QED':2})
 
 GC_123 = Coupling(name = 'GC_123',
-                  value = 'complex(0,1)*G**2*I1063 + complex(0,1)*G**2*I963',
-                  order = {'QCD':2})
+                  value = '(ee**2*complex(0,1))/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_124 = Coupling(name = 'GC_124',
-                  value = '-(complex(0,1)*G*I1066) - complex(0,1)*G*I966',
-                  order = {'QCD':1})
+                  value = '-((ee**2*complex(0,1))/sw**2)',
+                  order = {'QED':2})
 
 GC_125 = Coupling(name = 'GC_125',
-                  value = 'complex(0,1)*G**2*I1066 + complex(0,1)*G**2*I966',
-                  order = {'QCD':2})
+                  value = '(ee**2*complex(0,1)*I9611)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_126 = Coupling(name = 'GC_126',
-                  value = '-(complex(0,1)*G*Rd11*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9622)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_127 = Coupling(name = 'GC_127',
-                  value = '-(complex(0,1)*G*Rd22*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9633)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_128 = Coupling(name = 'GC_128',
-                  value = '-(complex(0,1)*G*Rd33*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9636)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_129 = Coupling(name = 'GC_129',
-                  value = 'complex(0,1)*G*Rd36*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9663)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_130 = Coupling(name = 'GC_130',
-                  value = '-(complex(0,1)*I3133*NN13) - (ee*complex(0,1)*NN11*Rd36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9666)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_131 = Coupling(name = 'GC_131',
-                  value = '-(complex(0,1)*I3133*NN23) - (ee*complex(0,1)*NN21*Rd36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9711)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_132 = Coupling(name = 'GC_132',
-                  value = '-(complex(0,1)*I3133*NN33) - (ee*complex(0,1)*NN31*Rd36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9722)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_133 = Coupling(name = 'GC_133',
-                  value = '-(complex(0,1)*I3133*NN43) - (ee*complex(0,1)*NN41*Rd36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9733)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_134 = Coupling(name = 'GC_134',
-                  value = 'complex(0,1)*G*Rd44*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9736)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_135 = Coupling(name = 'GC_135',
-                  value = '-(ee*complex(0,1)*NN11*Rd44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9763)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_136 = Coupling(name = 'GC_136',
-                  value = '-(ee*complex(0,1)*NN21*Rd44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9766)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_137 = Coupling(name = 'GC_137',
-                  value = '-(ee*complex(0,1)*NN31*Rd44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9811)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_138 = Coupling(name = 'GC_138',
-                  value = '-(ee*complex(0,1)*NN41*Rd44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9822)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_139 = Coupling(name = 'GC_139',
-                  value = 'complex(0,1)*G*Rd55*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I9833)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_140 = Coupling(name = 'GC_140',
-                  value = '-(ee*complex(0,1)*NN11*Rd55*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9836)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_141 = Coupling(name = 'GC_141',
-                  value = '-(ee*complex(0,1)*NN21*Rd55*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9863)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_142 = Coupling(name = 'GC_142',
-                  value = '-(ee*complex(0,1)*NN31*Rd55*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I9866)/(2.*sw**2)',
+                  order = {'QED':2})
 
 GC_143 = Coupling(name = 'GC_143',
-                  value = '-(ee*complex(0,1)*NN41*Rd55*cmath.sqrt(2))/(3.*cw)',
+                  value = '(ee*complex(0,1))/(sw*cmath.sqrt(2))',
                   order = {'QED':1})
 
 GC_144 = Coupling(name = 'GC_144',
-                  value = '-(complex(0,1)*G*Rd63*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(cw*ee*complex(0,1))/sw',
+                  order = {'QED':1})
 
 GC_145 = Coupling(name = 'GC_145',
-                  value = 'complex(0,1)*G*Rd66*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(-2*cw*ee**2*complex(0,1))/sw',
+                  order = {'QED':2})
 
 GC_146 = Coupling(name = 'GC_146',
-                  value = '-(complex(0,1)*I3136*NN13) - (ee*complex(0,1)*NN11*Rd66*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I3911)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_147 = Coupling(name = 'GC_147',
-                  value = '-(complex(0,1)*I3136*NN23) - (ee*complex(0,1)*NN21*Rd66*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I3922)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_148 = Coupling(name = 'GC_148',
-                  value = '-(complex(0,1)*I3136*NN33) - (ee*complex(0,1)*NN31*Rd66*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I3933)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_149 = Coupling(name = 'GC_149',
-                  value = '-(complex(0,1)*I3136*NN43) - (ee*complex(0,1)*NN41*Rd66*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I3936)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_150 = Coupling(name = 'GC_150',
-                  value = '-(complex(0,1)*I3833*NN13) - (ee*complex(0,1)*NN11*Rl36*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I4511)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_151 = Coupling(name = 'GC_151',
-                  value = '-(complex(0,1)*I3833*NN23) - (ee*complex(0,1)*NN21*Rl36*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I4522)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_152 = Coupling(name = 'GC_152',
-                  value = '-(complex(0,1)*I3833*NN33) - (ee*complex(0,1)*NN31*Rl36*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I4533)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_153 = Coupling(name = 'GC_153',
-                  value = '-(complex(0,1)*I3833*NN43) - (ee*complex(0,1)*NN41*Rl36*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '-((ee**2*complex(0,1)*I4536)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_154 = Coupling(name = 'GC_154',
-                  value = '-((ee*complex(0,1)*NN11*Rl44*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5311)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_155 = Coupling(name = 'GC_155',
-                  value = '-((ee*complex(0,1)*NN21*Rl44*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5311*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_156 = Coupling(name = 'GC_156',
-                  value = '-((ee*complex(0,1)*NN31*Rl44*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5322)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_157 = Coupling(name = 'GC_157',
-                  value = '-((ee*complex(0,1)*NN41*Rl44*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5322*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_158 = Coupling(name = 'GC_158',
-                  value = '-((ee*complex(0,1)*NN11*Rl55*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5333)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_159 = Coupling(name = 'GC_159',
-                  value = '-((ee*complex(0,1)*NN21*Rl55*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5333*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_160 = Coupling(name = 'GC_160',
-                  value = '-((ee*complex(0,1)*NN31*Rl55*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5336)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_161 = Coupling(name = 'GC_161',
-                  value = '-((ee*complex(0,1)*NN41*Rl55*cmath.sqrt(2))/cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5336*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_162 = Coupling(name = 'GC_162',
-                  value = '-(complex(0,1)*I3836*NN13) - (ee*complex(0,1)*NN11*Rl66*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5363)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_163 = Coupling(name = 'GC_163',
-                  value = '-(complex(0,1)*I3836*NN23) - (ee*complex(0,1)*NN21*Rl66*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5363*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_164 = Coupling(name = 'GC_164',
-                  value = '-(complex(0,1)*I3836*NN33) - (ee*complex(0,1)*NN31*Rl66*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I5366)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_165 = Coupling(name = 'GC_165',
-                  value = '-(complex(0,1)*I3836*NN43) - (ee*complex(0,1)*NN41*Rl66*cmath.sqrt(2))/cw',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I5366*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_166 = Coupling(name = 'GC_166',
-                  value = '-(complex(0,1)*G*Ru11*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I6611)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_167 = Coupling(name = 'GC_167',
-                  value = '-(complex(0,1)*G*Ru22*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee*complex(0,1)*G*I6611*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_168 = Coupling(name = 'GC_168',
-                  value = '-(complex(0,1)*G*Ru33*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I6622)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_169 = Coupling(name = 'GC_169',
-                  value = 'complex(0,1)*G*Ru36*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(ee*complex(0,1)*G*I6622*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_170 = Coupling(name = 'GC_170',
-                  value = '-(complex(0,1)*I4933*NN14) + (2*ee*complex(0,1)*NN11*Ru36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I6633)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_171 = Coupling(name = 'GC_171',
-                  value = '-(complex(0,1)*I4933*NN24) + (2*ee*complex(0,1)*NN21*Ru36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I6633*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_172 = Coupling(name = 'GC_172',
-                  value = '-(complex(0,1)*I4933*NN34) + (2*ee*complex(0,1)*NN31*Ru36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I6636)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_173 = Coupling(name = 'GC_173',
-                  value = '-(complex(0,1)*I4933*NN44) + (2*ee*complex(0,1)*NN41*Ru36*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I6636*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_174 = Coupling(name = 'GC_174',
-                  value = 'complex(0,1)*G*Ru44*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '(ee**2*complex(0,1)*I6663)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_175 = Coupling(name = 'GC_175',
-                  value = '(2*ee*complex(0,1)*NN11*Ru44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I6663*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_176 = Coupling(name = 'GC_176',
-                  value = '(2*ee*complex(0,1)*NN21*Ru44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee**2*complex(0,1)*I6666)/(3.*sw*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_177 = Coupling(name = 'GC_177',
-                  value = '(2*ee*complex(0,1)*NN31*Ru44*cmath.sqrt(2))/(3.*cw)',
-                  order = {'QED':1})
+                  value = '(ee*complex(0,1)*G*I6666*cmath.sqrt(2))/sw',
+                  order = {'QCD':1,'QED':1})
 
 GC_178 = Coupling(name = 'GC_178',
-                  value = '(2*ee*complex(0,1)*NN41*Ru44*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9211)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_179 = Coupling(name = 'GC_179',
-                  value = 'complex(0,1)*G*Ru55*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '-((ee*complex(0,1)*I9222)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':1})
 
 GC_180 = Coupling(name = 'GC_180',
-                  value = '(2*ee*complex(0,1)*NN11*Ru55*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9233)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_181 = Coupling(name = 'GC_181',
-                  value = '(2*ee*complex(0,1)*NN21*Ru55*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9236)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_182 = Coupling(name = 'GC_182',
-                  value = '(2*ee*complex(0,1)*NN31*Ru55*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9263)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_183 = Coupling(name = 'GC_183',
-                  value = '(2*ee*complex(0,1)*NN41*Ru55*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9266)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_184 = Coupling(name = 'GC_184',
-                  value = '-(complex(0,1)*G*Ru63*cmath.sqrt(2))',
-                  order = {'QCD':1})
+                  value = '-((ee*complex(0,1)*I9311)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':1})
 
 GC_185 = Coupling(name = 'GC_185',
-                  value = 'complex(0,1)*G*Ru66*cmath.sqrt(2)',
-                  order = {'QCD':1})
+                  value = '-((ee*complex(0,1)*I9322)/(sw*cmath.sqrt(2)))',
+                  order = {'QED':1})
 
 GC_186 = Coupling(name = 'GC_186',
-                  value = '-(complex(0,1)*I4936*NN14) + (2*ee*complex(0,1)*NN11*Ru66*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9333)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_187 = Coupling(name = 'GC_187',
-                  value = '-(complex(0,1)*I4936*NN24) + (2*ee*complex(0,1)*NN21*Ru66*cmath.sqrt(2))/(3.*cw)',
+                  value = '-((ee*complex(0,1)*I9336)/(sw*cmath.sqrt(2)))',
                   order = {'QED':1})
 
 GC_188 = Coupling(name = 'GC_188',
-                  value = '-(complex(0,1)*I4936*NN34) + (2*ee*complex(0,1)*NN31*Ru66*cmath.sqrt(2))/(3.*cw)',
+                  value = '(ee*complex(0,1)*I9411)/(sw*cmath.sqrt(2))',
                   order = {'QED':1})
 
 GC_189 = Coupling(name = 'GC_189',
-                  value = '-(complex(0,1)*I4936*NN44) + (2*ee*complex(0,1)*NN41*Ru66*cmath.sqrt(2))/(3.*cw)',
+                  value = '(ee*complex(0,1)*I9422)/(sw*cmath.sqrt(2))',
                   order = {'QED':1})
 
 GC_190 = Coupling(name = 'GC_190',
-                  value = '(ee**2*complex(0,1))/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9433)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_191 = Coupling(name = 'GC_191',
-                  value = '-((ee**2*complex(0,1))/sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9436)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_192 = Coupling(name = 'GC_192',
-                  value = '(cw**2*ee**2*complex(0,1))/sw**2',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9463)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_193 = Coupling(name = 'GC_193',
-                  value = '(ee**2*complex(0,1)*I2511)/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9466)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_194 = Coupling(name = 'GC_194',
-                  value = '(ee**2*complex(0,1)*I2522)/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9511)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_195 = Coupling(name = 'GC_195',
-                  value = '(ee**2*complex(0,1)*I2533)/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9522)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_196 = Coupling(name = 'GC_196',
-                  value = '(ee**2*complex(0,1)*I2536)/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9533)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_197 = Coupling(name = 'GC_197',
-                  value = '(ee**2*complex(0,1)*I2563)/(2.*sw**2)',
-                  order = {'QED':2})
+                  value = '(ee*complex(0,1)*I9536)/(sw*cmath.sqrt(2))',
+                  order = {'QED':1})
 
 GC_198 = Coupling(name = 'GC_198',
-                  value = '(ee**2*complex(0,1)*I2566)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9211)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_199 = Coupling(name = 'GC_199',
-                  value = '(ee**2*complex(0,1)*I2611)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9222)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_200 = Coupling(name = 'GC_200',
-                  value = '(ee**2*complex(0,1)*I2622)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9233)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_201 = Coupling(name = 'GC_201',
-                  value = '(ee**2*complex(0,1)*I2633)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9236)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_202 = Coupling(name = 'GC_202',
-                  value = '(ee**2*complex(0,1)*I2636)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9263)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_203 = Coupling(name = 'GC_203',
-                  value = '(ee**2*complex(0,1)*I2663)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9266)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_204 = Coupling(name = 'GC_204',
-                  value = '(ee**2*complex(0,1)*I2666)/(2.*sw**2)',
+                  value = '-((cw*ee**2*complex(0,1)*I9311)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
                   order = {'QED':2})
 
 GC_205 = Coupling(name = 'GC_205',
-                  value = '(ee**2*complex(0,1)*I911)/(2.*sw**2)',
+                  value = '-((cw*ee**2*complex(0,1)*I9322)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
                   order = {'QED':2})
 
 GC_206 = Coupling(name = 'GC_206',
-                  value = '(ee**2*complex(0,1)*I922)/(2.*sw**2)',
+                  value = '-((cw*ee**2*complex(0,1)*I9333)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
                   order = {'QED':2})
 
 GC_207 = Coupling(name = 'GC_207',
-                  value = '(ee**2*complex(0,1)*I933)/(2.*sw**2)',
+                  value = '-((cw*ee**2*complex(0,1)*I9336)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
                   order = {'QED':2})
 
 GC_208 = Coupling(name = 'GC_208',
-                  value = '(ee**2*complex(0,1)*I936)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9411)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_209 = Coupling(name = 'GC_209',
-                  value = '(ee**2*complex(0,1)*I963)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9422)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_210 = Coupling(name = 'GC_210',
-                  value = '(ee**2*complex(0,1)*I966)/(2.*sw**2)',
+                  value = '(cw*ee**2*complex(0,1)*I9433)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
                   order = {'QED':2})
 
 GC_211 = Coupling(name = 'GC_211',
-                  value = '(ee*complex(0,1))/(sw*cmath.sqrt(2))',
-                  order = {'QED':1})
+                  value = '(cw*ee**2*complex(0,1)*I9436)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_212 = Coupling(name = 'GC_212',
-                  value = '-(cw*ee*complex(0,1))/(2.*sw)',
-                  order = {'QED':1})
+                  value = '(cw*ee**2*complex(0,1)*I9463)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_213 = Coupling(name = 'GC_213',
-                  value = '(cw*ee*complex(0,1))/(2.*sw)',
-                  order = {'QED':1})
+                  value = '(cw*ee**2*complex(0,1)*I9466)/(3.*(-1 + sw)*(1 + sw)*cmath.sqrt(2))',
+                  order = {'QED':2})
 
 GC_214 = Coupling(name = 'GC_214',
-                  value = '(cw*ee*complex(0,1))/sw',
-                  order = {'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I9511)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_215 = Coupling(name = 'GC_215',
-                  value = '(-2*cw*ee**2*complex(0,1))/sw',
+                  value = '-((cw*ee**2*complex(0,1)*I9522)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
                   order = {'QED':2})
 
 GC_216 = Coupling(name = 'GC_216',
-                  value = '(ee*complex(0,1)*I1511)/(sw*cmath.sqrt(2))',
-                  order = {'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I9533)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_217 = Coupling(name = 'GC_217',
-                  value = '(ee*complex(0,1)*I1522)/(sw*cmath.sqrt(2))',
-                  order = {'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I9536)/((-1 + sw)*(1 + sw)*cmath.sqrt(2)))',
+                  order = {'QED':2})
 
 GC_218 = Coupling(name = 'GC_218',
-                  value = '(ee*complex(0,1)*I1533)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN11*Rd44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_219 = Coupling(name = 'GC_219',
-                  value = '(ee*complex(0,1)*I1536)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN21*Rd44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_220 = Coupling(name = 'GC_220',
-                  value = '(ee*complex(0,1)*I1563)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN31*Rd44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_221 = Coupling(name = 'GC_221',
-                  value = '(ee*complex(0,1)*I1566)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN41*Rd44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_222 = Coupling(name = 'GC_222',
-                  value = '-((ee*complex(0,1)*I1611)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN11*Rd55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_223 = Coupling(name = 'GC_223',
-                  value = '-((ee*complex(0,1)*I1622)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN21*Rd55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_224 = Coupling(name = 'GC_224',
-                  value = '-((ee*complex(0,1)*I1633)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN31*Rd55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_225 = Coupling(name = 'GC_225',
-                  value = '-((ee*complex(0,1)*I1636)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN41*Rd55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_226 = Coupling(name = 'GC_226',
-                  value = '-((ee*complex(0,1)*I1663)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN11*Rl44*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_227 = Coupling(name = 'GC_227',
-                  value = '-((ee*complex(0,1)*I1666)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN21*Rl44*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_228 = Coupling(name = 'GC_228',
-                  value = '-((ee*complex(0,1)*I1711)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN31*Rl44*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_229 = Coupling(name = 'GC_229',
-                  value = '-((ee*complex(0,1)*I1722)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN41*Rl44*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_230 = Coupling(name = 'GC_230',
-                  value = '-((ee*complex(0,1)*I1733)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN11*Rl55*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_231 = Coupling(name = 'GC_231',
-                  value = '-((ee*complex(0,1)*I1736)/(sw*cmath.sqrt(2)))',
+                  value = '(cw*ee*complex(0,1)*NN21*Rl55*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_232 = Coupling(name = 'GC_232',
-                  value = '(ee*complex(0,1)*I1811)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN31*Rl55*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_233 = Coupling(name = 'GC_233',
-                  value = '(ee*complex(0,1)*I1822)/(sw*cmath.sqrt(2))',
+                  value = '(cw*ee*complex(0,1)*NN41*Rl55*cmath.sqrt(2))/((-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_234 = Coupling(name = 'GC_234',
-                  value = '(ee*complex(0,1)*I1833)/(sw*cmath.sqrt(2))',
+                  value = '(-2*cw*ee*complex(0,1)*NN11*Ru44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_235 = Coupling(name = 'GC_235',
-                  value = '(ee*complex(0,1)*I1836)/(sw*cmath.sqrt(2))',
+                  value = '(-2*cw*ee*complex(0,1)*NN21*Ru44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QED':1})
 
 GC_236 = Coupling(name = 'GC_236',
-                  value = '-((ee**2*complex(0,1)*I1911)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN31*Ru44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_237 = Coupling(name = 'GC_237',
-                  value = '-((ee**2*complex(0,1)*I1922)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN41*Ru44*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_238 = Coupling(name = 'GC_238',
-                  value = '-((ee**2*complex(0,1)*I1933)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN11*Ru55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_239 = Coupling(name = 'GC_239',
-                  value = '-((ee**2*complex(0,1)*I1936)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN21*Ru55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_240 = Coupling(name = 'GC_240',
-                  value = '-((ee**2*complex(0,1)*I2011)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN31*Ru55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_241 = Coupling(name = 'GC_241',
-                  value = '-((ee**2*complex(0,1)*I2022)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*NN41*Ru55*cmath.sqrt(2))/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_242 = Coupling(name = 'GC_242',
-                  value = '-((ee**2*complex(0,1)*I2033)/(sw*cmath.sqrt(2)))',
+                  value = '-(ee**2*complex(0,1))/(2.*(-1 + sw)*sw**2*(1 + sw))',
                   order = {'QED':2})
 
 GC_243 = Coupling(name = 'GC_243',
-                  value = '-((ee**2*complex(0,1)*I2036)/(sw*cmath.sqrt(2)))',
-                  order = {'QED':2})
+                  value = '-(cw*ee*complex(0,1))/(2.*(-1 + sw)*sw*(1 + sw))',
+                  order = {'QED':1})
 
 GC_244 = Coupling(name = 'GC_244',
-                  value = '(ee*complex(0,1)*G*I2111*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(cw*ee*complex(0,1))/(2.*(-1 + sw)*sw*(1 + sw))',
+                  order = {'QED':1})
 
 GC_245 = Coupling(name = 'GC_245',
-                  value = '(ee*complex(0,1)*G*I2122*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '-(cw*ee*complex(0,1)*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_246 = Coupling(name = 'GC_246',
-                  value = '(ee*complex(0,1)*G*I2133*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(2*cw*ee*complex(0,1)*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_247 = Coupling(name = 'GC_247',
-                  value = '(ee*complex(0,1)*G*I2136*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '-((cw*ee*complex(0,1)*sw)/((-1 + sw)*(1 + sw)))',
+                  order = {'QED':1})
 
 GC_248 = Coupling(name = 'GC_248',
-                  value = '(ee*complex(0,1)*G*I2163*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(cw*ee*complex(0,1)*I10044*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_249 = Coupling(name = 'GC_249',
-                  value = '(ee*complex(0,1)*G*I2166*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(cw*ee*complex(0,1)*I10055*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_250 = Coupling(name = 'GC_250',
-                  value = '(ee**2*complex(0,1)*I2211)/(3.*sw*cmath.sqrt(2))',
-                  order = {'QED':2})
+                  value = '(cw*ee*complex(0,1)*I10144*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_251 = Coupling(name = 'GC_251',
-                  value = '(ee**2*complex(0,1)*I2222)/(3.*sw*cmath.sqrt(2))',
-                  order = {'QED':2})
+                  value = '(cw*ee*complex(0,1)*I10155*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_252 = Coupling(name = 'GC_252',
-                  value = '(ee**2*complex(0,1)*I2233)/(3.*sw*cmath.sqrt(2))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*I10244*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_253 = Coupling(name = 'GC_253',
-                  value = '(ee**2*complex(0,1)*I2236)/(3.*sw*cmath.sqrt(2))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*I10255*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':1})
 
 GC_254 = Coupling(name = 'GC_254',
-                  value = '(ee**2*complex(0,1)*I2263)/(3.*sw*cmath.sqrt(2))',
+                  value = '(2*cw*ee**2*complex(0,1)*I2644*sw)/((-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_255 = Coupling(name = 'GC_255',
-                  value = '(ee**2*complex(0,1)*I2266)/(3.*sw*cmath.sqrt(2))',
+                  value = '(2*cw*ee**2*complex(0,1)*I2655*sw)/((-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_256 = Coupling(name = 'GC_256',
-                  value = '(ee*complex(0,1)*G*I2311*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(8*cw*ee**2*complex(0,1)*I5244*sw)/(9.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_257 = Coupling(name = 'GC_257',
-                  value = '(ee*complex(0,1)*G*I2322*cmath.sqrt(2))/sw',
+                  value = '(4*cw*ee*complex(0,1)*G*I5244*sw)/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QCD':1,'QED':1})
 
 GC_258 = Coupling(name = 'GC_258',
-                  value = '(ee*complex(0,1)*G*I2333*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(8*cw*ee**2*complex(0,1)*I5255*sw)/(9.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_259 = Coupling(name = 'GC_259',
-                  value = '(ee*complex(0,1)*G*I2336*cmath.sqrt(2))/sw',
+                  value = '(4*cw*ee*complex(0,1)*G*I5255*sw)/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QCD':1,'QED':1})
 
 GC_260 = Coupling(name = 'GC_260',
-                  value = '(ee*complex(0,1)*G*I2363*cmath.sqrt(2))/sw',
-                  order = {'QCD':1,'QED':1})
+                  value = '(2*cw*ee**2*complex(0,1)*I644*sw)/(9.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_261 = Coupling(name = 'GC_261',
-                  value = '(ee*complex(0,1)*G*I2366*cmath.sqrt(2))/sw',
+                  value = '(-2*cw*ee*complex(0,1)*G*I644*sw)/(3.*(-1 + sw)*(1 + sw))',
                   order = {'QCD':1,'QED':1})
 
 GC_262 = Coupling(name = 'GC_262',
-                  value = '(ee**2*complex(0,1)*I2411)/(3.*sw*cmath.sqrt(2))',
+                  value = '(2*cw*ee**2*complex(0,1)*I655*sw)/(9.*(-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_263 = Coupling(name = 'GC_263',
-                  value = '(ee**2*complex(0,1)*I2422)/(3.*sw*cmath.sqrt(2))',
-                  order = {'QED':2})
+                  value = '(-2*cw*ee*complex(0,1)*G*I655*sw)/(3.*(-1 + sw)*(1 + sw))',
+                  order = {'QCD':1,'QED':1})
 
 GC_264 = Coupling(name = 'GC_264',
-                  value = '(ee**2*complex(0,1)*I2433)/(3.*sw*cmath.sqrt(2))',
+                  value = '(-2*ee**2*complex(0,1)*I10044*sw**2)/(9.*(-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_265 = Coupling(name = 'GC_265',
-                  value = '(ee**2*complex(0,1)*I2436)/(3.*sw*cmath.sqrt(2))',
+                  value = '(-2*ee**2*complex(0,1)*I10055*sw**2)/(9.*(-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_266 = Coupling(name = 'GC_266',
-                  value = '(ee**2*complex(0,1)*I2463)/(3.*sw*cmath.sqrt(2))',
+                  value = '(-2*ee**2*complex(0,1)*I10144*sw**2)/((-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_267 = Coupling(name = 'GC_267',
-                  value = '(ee**2*complex(0,1)*I2466)/(3.*sw*cmath.sqrt(2))',
+                  value = '(-2*ee**2*complex(0,1)*I10155*sw**2)/((-1 + sw)*(1 + sw))',
                   order = {'QED':2})
 
 GC_268 = Coupling(name = 'GC_268',
-                  value = '-(ee*complex(0,1)*sw)/(6.*cw)',
-                  order = {'QED':1})
+                  value = '(-8*ee**2*complex(0,1)*I10244*sw**2)/(9.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_269 = Coupling(name = 'GC_269',
-                  value = '(ee*complex(0,1)*sw)/(2.*cw)',
-                  order = {'QED':1})
+                  value = '(-8*ee**2*complex(0,1)*I10255*sw**2)/(9.*(-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_270 = Coupling(name = 'GC_270',
-                  value = '-(ee*complex(0,1)*I1044*sw)/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I2511)/((-1 + sw)*sw*(1 + sw))) + (2*cw*ee**2*complex(0,1)*I2511*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_271 = Coupling(name = 'GC_271',
-                  value = '(2*ee*complex(0,1)*G*I1044*sw)/(3.*cw)',
-                  order = {'QCD':1,'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I2522)/((-1 + sw)*sw*(1 + sw))) + (2*cw*ee**2*complex(0,1)*I2522*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_272 = Coupling(name = 'GC_272',
-                  value = '-(ee*complex(0,1)*I1055*sw)/(3.*cw)',
-                  order = {'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I2533)/((-1 + sw)*sw*(1 + sw))) + (2*cw*ee**2*complex(0,1)*I2533*sw)/((-1 + sw)*(1 + sw)) + (2*cw*ee**2*complex(0,1)*I2633*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_273 = Coupling(name = 'GC_273',
-                  value = '(2*ee*complex(0,1)*G*I1055*sw)/(3.*cw)',
-                  order = {'QCD':1,'QED':1})
+                  value = '-((cw*ee**2*complex(0,1)*I2536)/((-1 + sw)*sw*(1 + sw))) + (2*cw*ee**2*complex(0,1)*I2536*sw)/((-1 + sw)*(1 + sw)) + (2*cw*ee**2*complex(0,1)*I2636*sw)/((-1 + sw)*(1 + sw))',
+                  order = {'QED':2})
 
 GC_274 = Coupling(name = 'GC_274',
-                  value = '(-8*ee**2*complex(0,1)*I1244*sw)/(9.*cw)',
+                  value = '-((cw*ee**2*complex(0,1)*I2563)/((-1 + sw)*sw*(1 + sw))) + (2*cw*ee**2*complex(0,1)*I2563*sw)/((-1 + sw)*(1 + sw)) + (2*cw*ee**2*complex(0,1)*I2663*sw)/((-1 + sw)*(1 + sw))',
                   order = {'QED':2})