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

[Olivier Mattelaer]

Add a check for charge conservation.
This checks all the charge define at UFO level as well as the electric charge.
This first check that all vertices preserve the charge conservation.

This provides a small optimization but also nicer error message.

[lp:~maddevelopers/mg5amcnlo/check_charge_conservation]

81. By Olivier Mattelaer on 2011-03-09

avoid double definition of InvalidCmd (only in __init__)

=== modified file 'madgraph/__init__.py'
--- madgraph/__init__.py	2010-10-19 14:51:20 +0000
+++ madgraph/__init__.py	2011-03-09 08:34:17 +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-09 08:34:17 +0000
@@ -87,12 +87,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
 
@@ -665,6 +665,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 +675,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
 

=== 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-09 08:34:17 +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
@@ -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/interface/cmd_interface.py'
--- madgraph/interface/cmd_interface.py	2011-01-27 12:30:18 +0000
+++ madgraph/interface/cmd_interface.py	2011-03-09 08:34:17 +0000
@@ -35,7 +35,7 @@
 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
@@ -592,9 +592,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 +602,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 +621,20 @@
         
         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")
+            raise InvalidCmd("No particle list currently active, please import a model first")
 
         if self._curr_model['particles'].find_name(args[0]):
             raise MadGraph5Error("label %s is a particle name in this model\n\
@@ -650,15 +647,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 +665,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 +681,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 +694,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 +712,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 +723,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 +736,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 +744,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,35 +763,35 @@
         
         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 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:
@@ -807,21 +804,21 @@
         
         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 +830,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,7 +856,7 @@
 
         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 get_default_path(self):
         """Set self._export_dir to the default (\'auto\') path"""
@@ -889,7 +886,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')    
             
         
@@ -1488,7 +1485,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)
@@ -1512,7 +1509,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 +1564,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 +1801,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:
@@ -2366,13 +2363,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"""

=== modified file 'madgraph/iolibs/files.py'
--- madgraph/iolibs/files.py	2010-11-06 18:34:32 +0000
+++ madgraph/iolibs/files.py	2011-03-09 08:34:17 +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=1299607013):
     """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/various/process_checks.py'
--- madgraph/various/process_checks.py	2011-01-26 13:23:54 +0000
+++ madgraph/various/process_checks.py	2011-03-09 08:34:17 +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-09 08:34:17 +0000
@@ -107,6 +107,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 +133,7 @@
         for interaction_info in self.ufomodel.all_vertices:
             self.add_interaction(interaction_info)
         
+        self.model.set('conserved_charge', self.conservecharge)
         return self.model
         
     
@@ -171,6 +173,10 @@
                     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 +187,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 +229,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*)\)')

=== modified file 'tests/acceptance_tests/test_cmd.py'
--- tests/acceptance_tests/test_cmd.py	2010-12-30 23:35:14 +0000
+++ tests/acceptance_tests/test_cmd.py	2011-03-09 08:34:17 +0000
@@ -50,7 +50,7 @@
         
     def test_generate(self):
         """command 'generate' works"""
-    
+        
         self.do('load model %s' % self.join_path(_pickle_path, 'sm.pkl'))
         self.cmd._curr_model.pass_particles_name_in_mg_default()
         self.do('generate e+ e- > e+ e-')
@@ -85,7 +85,7 @@
 
     def test_draw(self):
         """ command 'draw' works """
-        
+
         self.do('load processes %s' % self.join_path(_pickle_path,'e+e-_e+e-.pkl'))
         self.do('draw .')
         self.assertTrue(os.path.exists('./diagrams_0_epem_epem.eps'))
@@ -128,7 +128,7 @@
 
         if os.path.isdir(self.out_dir):
             shutil.rmdir(self.out_dir)
-
+            
         self.do('load processes %s' % self.join_path(_pickle_path,'e+e-_e+e-.pkl'))
         self.do('output %s -nojpeg' % self.out_dir)
         self.assertTrue(os.path.exists(self.out_dir))

=== modified file 'tests/input_files/e+e-_e+e-.pkl'
--- tests/input_files/e+e-_e+e-.pkl	2010-11-05 02:28:04 +0000
+++ tests/input_files/e+e-_e+e-.pkl	2011-03-09 08:34:17 +0000
@@ -8947,7 +8947,7 @@
 stp1750
 Rp1751
 ssg58
-S'sm_v4'
+S'sm'
 p1752
 sS'ref_dict_to0'
 p1753
@@ -10534,15 +10534,23 @@
 ag1751
 atp2357
 Rp2358
+sS'conserved_charge'
+p2359
+c__builtin__
+set
+p2360
+((lp2361
+tp2362
+Rp2363
 sg50
 g0
 (g51
 g2
-(lp2359
+(lp2364
 g0
 (g53
 g5
-(dp2360
+(dp2365
 g55
 g56
 sg71
@@ -10571,12 +10579,12 @@
 g56
 sg74
 I1
-stp2361
-Rp2362
+stp2366
+Rp2367
 ag0
 (g53
 g5
-(dp2363
+(dp2368
 g55
 g134
 sg71
@@ -10605,12 +10613,12 @@
 g134
 sg74
 I2
-stp2364
-Rp2365
+stp2369
+Rp2370
 ag0
 (g53
 g5
-(dp2366
+(dp2371
 g55
 g174
 sg71
@@ -10639,12 +10647,12 @@
 g174
 sg74
 I3
-stp2367
-Rp2368
+stp2372
+Rp2373
 ag0
 (g53
 g5
-(dp2369
+(dp2374
 g55
 g214
 sg71
@@ -10673,12 +10681,12 @@
 g214
 sg74
 I4
-stp2370
-Rp2371
+stp2375
+Rp2376
 ag0
 (g53
 g5
-(dp2372
+(dp2377
 g55
 g254
 sg71
@@ -10707,12 +10715,12 @@
 g254
 sg74
 I5
-stp2373
-Rp2374
+stp2378
+Rp2379
 ag0
 (g53
 g5
-(dp2375
+(dp2380
 g55
 g294
 sg71
@@ -10741,12 +10749,12 @@
 g294
 sg74
 I6
-stp2376
-Rp2377
+stp2381
+Rp2382
 ag0
 (g53
 g5
-(dp2378
+(dp2383
 g55
 g630
 sg71
@@ -10775,12 +10783,12 @@
 g630
 sg74
 I11
-stp2379
-Rp2380
+stp2384
+Rp2385
 ag0
 (g53
 g5
-(dp2381
+(dp2386
 g55
 g656
 sg71
@@ -10809,12 +10817,12 @@
 g656
 sg74
 I13
-stp2382
-Rp2383
+stp2387
+Rp2388
 ag0
 (g53
 g5
-(dp2384
+(dp2389
 g55
 g682
 sg71
@@ -10843,12 +10851,12 @@
 g682
 sg74
 I15
-stp2385
-Rp2386
+stp2390
+Rp2391
 ag0
 (g53
 g5
-(dp2387
+(dp2392
 g55
 g986
 sg71
@@ -10877,12 +10885,12 @@
 g986
 sg74
 I12
-stp2388
-Rp2389
+stp2393
+Rp2394
 ag0
 (g53
 g5
-(dp2390
+(dp2395
 g55
 g1012
 sg71
@@ -10911,12 +10919,12 @@
 g1012
 sg74
 I14
-stp2391
-Rp2392
+stp2396
+Rp2397
 ag0
 (g53
 g5
-(dp2393
+(dp2398
 g55
 g1038
 sg71
@@ -10945,12 +10953,12 @@
 g1038
 sg74
 I16
-stp2394
-Rp2395
+stp2399
+Rp2400
 ag0
 (g53
 g5
-(dp2396
+(dp2401
 g55
 g81
 sg71
@@ -10979,12 +10987,12 @@
 g81
 sg74
 I21
-stp2397
-Rp2398
+stp2402
+Rp2403
 ag0
 (g53
 g5
-(dp2399
+(dp2404
 g55
 g714
 sg71
@@ -11013,12 +11021,12 @@
 g714
 sg74
 I23
-stp2400
-Rp2401
+stp2405
+Rp2406
 ag0
 (g53
 g5
-(dp2402
+(dp2407
 g55
 g421
 sg71
@@ -11047,12 +11055,12 @@
 g421
 sg74
 I22
-stp2403
-Rp2404
+stp2408
+Rp2409
 ag0
 (g53
 g5
-(dp2405
+(dp2410
 g55
 g1070
 sg71
@@ -11081,86 +11089,86 @@
 g1070
 sg74
 I-24
-stp2406
-Rp2407
-ag0
-(g53
-g5
-(dp2408
-g55
-g1411
-sg71
-I01
-sg58
-g1412
-sg60
-I01
-sg61
-I1
-sg62
-g1413
-sg64
-F0.0
-sg65
-g1414
-sg67
-g1415
-sg69
-g1416
-sg57
-I01
-sg72
-I1
-sg73
-g1411
-sg74
-I25
-stp2409
-Rp2410
-atp2411
+stp2411
 Rp2412
+ag0
+(g53
+g5
+(dp2413
+g55
+g1411
+sg71
+I01
+sg58
+g1412
+sg60
+I01
+sg61
+I1
+sg62
+g1413
+sg64
+F0.0
+sg65
+g1414
+sg67
+g1415
+sg69
+g1416
+sg57
+I01
+sg72
+I1
+sg73
+g1411
+sg74
+I25
+stp2414
+Rp2415
+atp2416
+Rp2417
 sg91
 NsS'got_majoranas'
-p2413
+p2418
 I00
 sS'particle_dict'
-p2414
-(dp2415
+p2419
+(dp2420
 I1
-g2362
+g2367
 sI2
-g2365
+g2370
 sI3
-g2368
+g2373
 sI4
-g2371
+g2376
 sI5
-g2374
+g2379
 sI6
-g2377
+g2382
 sI11
-g2380
+g2385
 sI12
-g2389
+g2394
 sI13
-g2383
+g2388
 sI14
-g2392
+g2397
 sI15
-g2386
+g2391
 sI16
-g2395
+g2400
 sI21
-g2398
+g2403
 sI22
-g2404
+g2409
 sI23
-g2401
+g2406
 sI24
 g0
 (g53
 g5
-(dp2416
+(dp2421
 g55
 g1070
 sg57
@@ -11189,15 +11197,15 @@
 g1070
 sg74
 I-24
-stp2417
-Rp2418
+stp2422
+Rp2423
 sI25
-g2410
+g2415
 sI-1
 g0
 (g53
 g5
-(dp2419
+(dp2424
 g55
 g56
 sg57
@@ -11226,15 +11234,15 @@
 g56
 sg74
 I1
-stp2420
-Rp2421
+stp2425
+Rp2426
 sI-24
-g2407
+g2412
 sI-16
 g0
 (g53
 g5
-(dp2422
+(dp2427
 g55
 g1038
 sg57
@@ -11263,13 +11271,13 @@
 g1038
 sg74
 I16
-stp2423
-Rp2424
+stp2428
+Rp2429
 sI-15
 g0
 (g53
 g5
-(dp2425
+(dp2430
 g55
 g682
 sg57
@@ -11298,13 +11306,13 @@
 g682
 sg74
 I15
-stp2426
-Rp2427
+stp2431
+Rp2432
 sI-14
 g0
 (g53
 g5
-(dp2428
+(dp2433
 g55
 g1012
 sg57
@@ -11333,13 +11341,13 @@
 g1012
 sg74
 I14
-stp2429
-Rp2430
+stp2434
+Rp2435
 sI-13
 g0
 (g53
 g5
-(dp2431
+(dp2436
 g55
 g656
 sg57
@@ -11368,13 +11376,13 @@
 g656
 sg74
 I13
-stp2432
-Rp2433
+stp2437
+Rp2438
 sI-12
 g0
 (g53
 g5
-(dp2434
+(dp2439
 g55
 g986
 sg57
@@ -11403,13 +11411,13 @@
 g986
 sg74
 I12
-stp2435
-Rp2436
+stp2440
+Rp2441
 sI-11
 g0
 (g53
 g5
-(dp2437
+(dp2442
 g55
 g630
 sg57
@@ -11438,13 +11446,13 @@
 g630
 sg74
 I11
-stp2438
-Rp2439
+stp2443
+Rp2444
 sI-6
 g0
 (g53
 g5
-(dp2440
+(dp2445
 g55
 g294
 sg57
@@ -11473,13 +11481,13 @@
 g294
 sg74
 I6
-stp2441
-Rp2442
+stp2446
+Rp2447
 sI-5
 g0
 (g53
 g5
-(dp2443
+(dp2448
 g55
 g254
 sg57
@@ -11508,13 +11516,13 @@
 g254
 sg74
 I5
-stp2444
-Rp2445
+stp2449
+Rp2450
 sI-4
 g0
 (g53
 g5
-(dp2446
+(dp2451
 g55
 g214
 sg57
@@ -11543,13 +11551,13 @@
 g214
 sg74
 I4
-stp2447
-Rp2448
+stp2452
+Rp2453
 sI-3
 g0
 (g53
 g5
-(dp2449
+(dp2454
 g55
 g174
 sg57
@@ -11578,13 +11586,13 @@
 g174
 sg74
 I3
-stp2450
-Rp2451
+stp2455
+Rp2456
 sI-2
 g0
 (g53
 g5
-(dp2452
+(dp2457
 g55
 g134
 sg57
@@ -11613,81 +11621,81 @@
 g134
 sg74
 I2
-stp2453
-Rp2454
+stp2458
+Rp2459
 ssS'parameters'
-p2455
+p2460
 Nsg124
-Nstp2456
-Rp2457
+Nstp2461
+Rp2462
 sS'overall_orders'
-p2458
-(dp2459
+p2463
+(dp2464
 sg17
 I0
 sS'decay_chains'
-p2460
+p2465
 g0
 (cmadgraph.core.base_objects
 ProcessList
-p2461
+p2466
 g2
-(lp2462
-tp2463
-Rp2464
-stp2465
-Rp2466
+(lp2467
+tp2468
+Rp2469
+stp2470
+Rp2471
 sS'diagrams'
-p2467
+p2472
 g0
 (cmadgraph.core.base_objects
 DiagramList
-p2468
+p2473
 g2
-(lp2469
+(lp2474
 g0
 (cmadgraph.core.base_objects
 Diagram
-p2470
+p2475
 g5
-(dp2471
+(dp2476
 S'vertices'
-p2472
+p2477
 g0
 (cmadgraph.core.base_objects
 VertexList
-p2473
+p2478
 g2
-(lp2474
+(lp2479
 g0
 (cmadgraph.core.base_objects
 Vertex
-p2475
-g5
-(dp2476
-g11
-g0
-(g12
-g2
-(lp2477
-g0
-(g14
-g5
-(dp2478
-g16
-I00
-sg17
-I11
-sg18
-I00
-sg19
-I1
-stp2479
-Rp2480
-ag0
-(g14
+p2480
 g5
 (dp2481
+g11
+g0
+(g12
+g2
+(lp2482
+g0
+(g14
+g5
+(dp2483
+g16
+I00
+sg17
+I11
+sg18
+I00
+sg19
+I1
+stp2484
+Rp2485
+ag0
+(g14
+g5
+(dp2486
 g16
 I00
 sg17
@@ -11696,55 +11704,55 @@
 I00
 sg19
 I2
-stp2482
-Rp2483
-ag0
-(g14
-g5
-(dp2484
-g16
-I01
-sg17
-I22
-sg18
-I00
-sg19
-I1
-stp2485
-Rp2486
-atp2487
+stp2487
 Rp2488
-sg17
-I15
-stp2489
-Rp2490
 ag0
-(g2475
-g5
-(dp2491
-g11
-g0
-(g12
-g2
-(lp2492
-g0
 (g14
 g5
-(dp2493
+(dp2489
 g16
 I01
 sg17
-I-11
+I22
 sg18
 I00
 sg19
-I3
+I1
+stp2490
+Rp2491
+atp2492
+Rp2493
+sg17
+I15
 stp2494
 Rp2495
 ag0
-(g14
+(g2480
 g5
 (dp2496
+g11
+g0
+(g12
+g2
+(lp2497
+g0
+(g14
+g5
+(dp2498
+g16
+I01
+sg17
+I-11
+sg18
+I00
+sg19
+I3
+stp2499
+Rp2500
+ag0
+(g14
+g5
+(dp2501
 g16
 I01
 sg17
@@ -11753,77 +11761,77 @@
 I00
 sg19
 I4
-stp2497
-Rp2498
-ag0
-(g14
-g5
-(dp2499
-g16
-I01
-sg17
-I22
-sg18
-I00
-sg19
-I3
-stp2500
-Rp2501
-atp2502
+stp2502
 Rp2503
+ag0
+(g14
+g5
+(dp2504
+g16
+I01
+sg17
+I22
+sg18
+I00
+sg19
+I3
+stp2505
+Rp2506
+atp2507
+Rp2508
 sg17
 I15
-stp2504
-Rp2505
+stp2509
+Rp2510
 ag0
-(g2475
+(g2480
 g5
-(dp2506
+(dp2511
 g11
 g0
 (g12
 g2
-(lp2507
-g2486
-ag2501
-atp2508
-Rp2509
+(lp2512
+g2491
+ag2506
+atp2513
+Rp2514
 sg17
 I0
-stp2510
-Rp2511
-atp2512
-Rp2513
-sg35
-(dp2514
-g647
-I2
-sstp2515
+stp2515
 Rp2516
-ag0
-(g2470
-g5
-(dp2517
-g2472
-g0
-(g2473
-g2
-(lp2518
-g0
-(g2475
-g5
+atp2517
+Rp2518
+sg35
 (dp2519
+g647
+I2
+sstp2520
+Rp2521
+ag0
+(g2475
+g5
+(dp2522
+g2477
+g0
+(g2478
+g2
+(lp2523
+g0
+(g2480
+g5
+(dp2524
 g11
 g0
 (g12
 g2
-(lp2520
-g2480
-ag2483
+(lp2525
+g2485
+ag2488
 ag0
 (g14
 g5
-(dp2521
+(dp2526
 g16
 I01
 sg17
@@ -11832,94 +11840,94 @@
 I00
 sg19
 I1
-stp2522
-Rp2523
-atp2524
-Rp2525
+stp2527
+Rp2528
+atp2529
+Rp2530
 sg17
 I24
-stp2526
-Rp2527
-ag0
-(g2475
-g5
-(dp2528
-g11
-g0
-(g12
-g2
-(lp2529
-g2495
-ag2498
-ag0
-(g14
-g5
-(dp2530
-g16
-I01
-sg17
-I23
-sg18
-I00
-sg19
-I3
 stp2531
 Rp2532
-atp2533
-Rp2534
+ag0
+(g2480
+g5
+(dp2533
+g11
+g0
+(g12
+g2
+(lp2534
+g2500
+ag2503
+ag0
+(g14
+g5
+(dp2535
+g16
+I01
+sg17
+I23
+sg18
+I00
+sg19
+I3
+stp2536
+Rp2537
+atp2538
+Rp2539
 sg17
 I24
-stp2535
-Rp2536
+stp2540
+Rp2541
 ag0
-(g2475
+(g2480
 g5
-(dp2537
+(dp2542
 g11
 g0
 (g12
 g2
-(lp2538
-g2523
-ag2532
-atp2539
-Rp2540
+(lp2543
+g2528
+ag2537
+atp2544
+Rp2545
 sg17
 I0
-stp2541
-Rp2542
-atp2543
-Rp2544
-sg35
-(dp2545
-g935
-I2
-sstp2546
+stp2546
 Rp2547
-ag0
-(g2470
-g5
-(dp2548
-g2472
-g0
-(g2473
-g2
-(lp2549
-g0
-(g2475
-g5
+atp2548
+Rp2549
+sg35
 (dp2550
+g935
+I2
+sstp2551
+Rp2552
+ag0
+(g2475
+g5
+(dp2553
+g2477
+g0
+(g2478
+g2
+(lp2554
+g0
+(g2480
+g5
+(dp2555
 g11
 g0
 (g12
 g2
-(lp2551
-g2480
-ag2495
+(lp2556
+g2485
+ag2500
 ag0
 (g14
 g5
-(dp2552
+(dp2557
 g16
 I00
 sg17
@@ -11928,94 +11936,94 @@
 I00
 sg19
 I1
-stp2553
-Rp2554
-atp2555
-Rp2556
+stp2558
+Rp2559
+atp2560
+Rp2561
 sg17
 I15
-stp2557
-Rp2558
-ag0
-(g2475
-g5
-(dp2559
-g11
-g0
-(g12
-g2
-(lp2560
-g2483
-ag2498
-ag0
-(g14
-g5
-(dp2561
-g16
-I00
-sg17
-I22
-sg18
-I00
-sg19
-I2
 stp2562
 Rp2563
-atp2564
-Rp2565
+ag0
+(g2480
+g5
+(dp2564
+g11
+g0
+(g12
+g2
+(lp2565
+g2488
+ag2503
+ag0
+(g14
+g5
+(dp2566
+g16
+I00
+sg17
+I22
+sg18
+I00
+sg19
+I2
+stp2567
+Rp2568
+atp2569
+Rp2570
 sg17
 I15
-stp2566
-Rp2567
+stp2571
+Rp2572
 ag0
-(g2475
+(g2480
 g5
-(dp2568
+(dp2573
 g11
 g0
 (g12
 g2
-(lp2569
-g2554
-ag2563
-atp2570
-Rp2571
+(lp2574
+g2559
+ag2568
+atp2575
+Rp2576
 sg17
 I0
-stp2572
-Rp2573
-atp2574
-Rp2575
-sg35
-(dp2576
-g647
-I2
-sstp2577
+stp2577
 Rp2578
-ag0
-(g2470
-g5
-(dp2579
-g2472
-g0
-(g2473
-g2
-(lp2580
-g0
-(g2475
-g5
+atp2579
+Rp2580
+sg35
 (dp2581
+g647
+I2
+sstp2582
+Rp2583
+ag0
+(g2475
+g5
+(dp2584
+g2477
+g0
+(g2478
+g2
+(lp2585
+g0
+(g2480
+g5
+(dp2586
 g11
 g0
 (g12
 g2
-(lp2582
-g2480
-ag2495
+(lp2587
+g2485
+ag2500
 ag0
 (g14
 g5
-(dp2583
+(dp2588
 g16
 I00
 sg17
@@ -12024,74 +12032,74 @@
 I00
 sg19
 I1
-stp2584
-Rp2585
-atp2586
-Rp2587
+stp2589
+Rp2590
+atp2591
+Rp2592
 sg17
 I24
-stp2588
-Rp2589
-ag0
-(g2475
-g5
-(dp2590
-g11
-g0
-(g12
-g2
-(lp2591
-g2483
-ag2498
-ag0
-(g14
-g5
-(dp2592
-g16
-I00
-sg17
-I23
-sg18
-I00
-sg19
-I2
 stp2593
 Rp2594
-atp2595
-Rp2596
+ag0
+(g2480
+g5
+(dp2595
+g11
+g0
+(g12
+g2
+(lp2596
+g2488
+ag2503
+ag0
+(g14
+g5
+(dp2597
+g16
+I00
+sg17
+I23
+sg18
+I00
+sg19
+I2
+stp2598
+Rp2599
+atp2600
+Rp2601
 sg17
 I24
-stp2597
-Rp2598
+stp2602
+Rp2603
 ag0
-(g2475
+(g2480
 g5
-(dp2599
+(dp2604
 g11
 g0
 (g12
 g2
-(lp2600
-g2585
-ag2594
-atp2601
-Rp2602
+(lp2605
+g2590
+ag2599
+atp2606
+Rp2607
 sg17
 I0
-stp2603
-Rp2604
-atp2605
-Rp2606
-sg35
-(dp2607
-g935
-I2
-sstp2608
+stp2608
 Rp2609
 atp2610
 Rp2611
-stp2612
-Rp2613
-atp2614
-Rp2615
+sg35
+(dp2612
+g935
+I2
+sstp2613
+Rp2614
+atp2615
+Rp2616
+stp2617
+Rp2618
+atp2619
+Rp2620
 .
\ No newline at end of file

=== modified file 'tests/input_files/sm.pkl'
--- tests/input_files/sm.pkl	2010-09-24 07:14:59 +0000
+++ tests/input_files/sm.pkl	2011-03-09 08:34:17 +0000
@@ -127,2574 +127,2540 @@
 sbasS'interaction_dict'
 p54
 (dp55
-sg13
+sS'name'
+p56
 S'sm'
-p56
+p57
 sS'ref_dict_to0'
-p57
-(dp58
+p58
+(dp59
+sS'ref_dict_to1'
+p60
+(dp61
 sS'interactions'
-p59
+p62
 g0
 (cmadgraph.core.base_objects
 InteractionList
-p60
+p63
 c__builtin__
 list
-p61
-(lp62
+p64
+(lp65
 g0
 (cmadgraph.core.base_objects
 Interaction
-p63
+p66
 g2
-(dp64
+(dp67
 S'particles'
-p65
+p68
 g0
 (cmadgraph.core.base_objects
 ParticleList
-p66
-g61
-(lp67
+p69
+g64
+(lp70
 g0
 (cmadgraph.core.base_objects
 Particle
-p68
+p71
 g2
-(dp69
+(dp72
 S'texname'
-p70
+p73
 S'\\phi'
-p71
-sS'propagating'
-p72
-I01
-sg13
-S'h'
-p73
-sS'self_antipart'
 p74
-I01
-sS'color'
+sS'is_part'
 p75
-I1
-sS'width'
+I01
+sg56
+S'h'
 p76
-S'WH'
+sS'self_antipart'
 p77
-sS'charge'
+I01
+sS'color'
 p78
-F0.0
-sS'mass'
+I1
+sS'width'
 p79
-S'MH'
+S'WH'
 p80
-sS'antiname'
+sS'charge'
 p81
-S'h'
+F0.0
+sS'mass'
 p82
-sS'line'
+S'MH'
 p83
-S'dashed'
+sS'LeptonNumber'
 p84
-sS'is_part'
+I0
+sS'antiname'
 p85
+S'h'
+p86
+sS'line'
+p87
+S'dashed'
+p88
+sS'propagating'
+p89
 I01
 sS'spin'
-p86
+p90
 I1
 sS'antitexname'
-p87
-g71
+p91
+g74
 sS'pdg_code'
-p88
+p92
 I25
-stp89
-Rp90
-ag90
-ag90
-atp91
-Rp92
+stp93
+Rp94
+ag94
+ag94
+atp95
+Rp96
 sS'lorentz'
-p93
-(lp94
+p97
+(lp98
 S'SSS1'
-p95
-asg75
-(lp96
+p99
+asS'orders'
+p100
+(dp101
+S'QED'
+p102
+I1
+ssg78
+(lp103
 g0
 (cmadgraph.core.color_algebra
 ColorString
-p97
-g61
-(lp98
-tp99
-Rp100
-(dp101
+p104
+g64
+(lp105
+tp106
+Rp107
+(dp108
 S'coeff'
-p102
+p109
 cfractions
 Fraction
-p103
+p110
 (S'1'
-p104
-tp105
-Rp106
+p111
+tp112
+Rp113
 sS'is_imaginary'
-p107
+p114
 I00
 sS'Nc_power'
-p108
+p115
 I0
-sbasS'id'
-p109
-I1
-sS'orders'
-p110
-(dp111
-S'QED'
-p112
-I1
-ssS'couplings'
-p113
-(dp114
+sbasS'couplings'
+p116
+(dp117
 (I0
 I0
-tp115
+tp118
 S'GC_21'
-p116
-sstp117
-Rp118
+p119
+ssS'id'
+p120
+I1
+stp121
+Rp122
 ag0
-(g63
-g2
-(dp119
-g65
-g0
 (g66
-g61
-(lp120
-g0
-(g68
-g2
-(dp121
-g70
+g2
+(dp123
+g68
+g0
+(g69
+g64
+(lp124
+g0
+(g71
+g2
+(dp125
+g73
 S'G'
-p122
-sg72
+p126
+sg75
 I01
-sg13
+sg56
 S'g'
-p123
-sg74
+p127
+sg77
 I01
-sg75
+sg78
 I8
-sg76
+sg79
 S'ZERO'
-p124
-sg78
+p128
+sg81
 F0.0
-sg79
-g124
-sg81
+sg82
+g128
+sg84
+I0
+sg85
 S'g'
-p125
-sg83
+p129
+sg87
 S'curly'
-p126
-sg85
+p130
+sg89
 I01
-sg86
+sg90
 I3
-sg87
-g122
-sg88
+sg91
+g126
+sg92
 I21
-stp127
-Rp128
-ag128
-ag128
-atp129
-Rp130
-sg93
-(lp131
+stp131
+Rp132
+ag132
+ag132
+atp133
+Rp134
+sg97
+(lp135
 S'VVV1'
-p132
-asg75
-(lp133
+p136
+asg100
+(dp137
+S'QCD'
+p138
+I1
+ssg78
+(lp139
 g0
-(g97
-g61
-(lp134
+(g104
+g64
+(lp140
 cmadgraph.core.color_algebra
 f
-p135
+p141
 (I0
 I1
 I2
-tp136
-Rp137
-atp138
-Rp139
-(dp140
-g102
-g106
-sg107
+tp142
+Rp143
+atp144
+Rp145
+(dp146
+g109
+g110
+(g111
+tp147
+Rp148
+sg114
 I00
-sg108
+sg115
 I0
-sbasg109
-I2
-sg110
-(dp141
-S'QCD'
-p142
-I1
-ssg113
-(dp143
+sbasg116
+(dp149
 (I0
 I0
-tp144
+tp150
 S'GC_4'
-p145
-sstp146
-Rp147
+p151
+ssg120
+I2
+stp152
+Rp153
 ag0
-(g63
+(g66
 g2
-(dp148
-g65
+(dp154
+g68
 g0
-(g66
-g61
-(lp149
-g128
-ag128
-ag128
-ag128
-atp150
-Rp151
-sg93
-(lp152
+(g69
+g64
+(lp155
+g132
+ag132
+ag132
+ag132
+atp156
+Rp157
+sg97
+(lp158
 S'VVVV1'
-p153
+p159
 aS'VVVV3'
-p154
+p160
 aS'VVVV4'
-p155
-asg75
-(lp156
+p161
+asg100
+(dp162
+g138
+I2
+ssg78
+(lp163
 g0
-(g97
-g61
-(lp157
-g135
-(I1
-I2
-I-1
-tp158
-Rp159
-ag135
-(I-1
-I0
-I3
-tp160
-Rp161
-atp162
-Rp163
-(dp164
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbag0
-(g97
-g61
-(lp165
-g135
-(I1
-I3
-I-1
-tp166
-Rp167
-ag135
-(I-1
-I0
-I2
-tp168
-Rp169
-atp170
-Rp171
-(dp172
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbag0
-(g97
-g61
-(lp173
-g135
-(I2
-I3
-I-1
-tp174
-Rp175
-ag135
-(I-1
-I0
-I1
-tp176
-Rp177
-atp178
-Rp179
-(dp180
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I3
-sg110
+(g104
+g64
+(lp164
+g141
+(I1
+I2
+I-1
+tp165
+Rp166
+ag141
+(I-1
+I0
+I3
+tp167
+Rp168
+atp169
+Rp170
+(dp171
+g109
+g110
+(g111
+tp172
+Rp173
+sg114
+I00
+sg115
+I0
+sbag0
+(g104
+g64
+(lp174
+g141
+(I1
+I3
+I-1
+tp175
+Rp176
+ag141
+(I-1
+I0
+I2
+tp177
+Rp178
+atp179
+Rp180
 (dp181
-g142
-I2
-ssg113
-(dp182
+g109
+g110
+(g111
+tp182
+Rp183
+sg114
+I00
+sg115
+I0
+sbag0
+(g104
+g64
+(lp184
+g141
 (I2
-I0
-tp183
-S'GC_6'
-p184
-s(I1
-I1
+I3
+I-1
 tp185
-g184
-s(I0
-I2
-tp186
-g184
-sstp187
+Rp186
+ag141
+(I-1
+I0
+I1
+tp187
 Rp188
-ag0
-(g63
-g2
-(dp189
-g65
-g0
-(g66
-g61
-(lp190
-g0
-(g68
-g2
+atp189
+Rp190
 (dp191
-g70
+g109
+g110
+(g111
+tp192
+Rp193
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp194
+(I2
+I0
+tp195
+S'GC_6'
+p196
+s(I1
+I1
+tp197
+g196
+s(I0
+I2
+tp198
+g196
+ssg120
+I3
+stp199
+Rp200
+ag0
+(g66
+g2
+(dp201
+g68
+g0
+(g69
+g64
+(lp202
+g0
+(g71
+g2
+(dp203
+g73
 S'A'
-p192
-sg72
-I01
-sg13
-S'a'
-p193
-sg74
-I01
+p204
 sg75
+I01
+sg56
+S'a'
+p205
+sg77
+I01
+sg78
 I1
-sg76
-g124
-sg78
-F0.0
 sg79
-g124
+g128
 sg81
+F0.0
+sg82
+g128
+sg84
+I0
+sg85
 S'a'
-p194
-sg83
+p206
+sg87
 S'wavy'
-p195
-sg85
+p207
+sg89
 I01
-sg86
+sg90
 I3
-sg87
-g192
-sg88
+sg91
+g204
+sg92
 I22
-stp196
-Rp197
+stp208
+Rp209
 ag0
-(g68
+(g71
 g2
-(dp198
-g70
+(dp210
+g73
 S'W+'
-p199
-sg85
-I00
-sg13
+p211
+sg89
+I01
+sg56
 S'w+'
-p200
-sg74
+p212
+sg77
 I00
-sg75
+sg78
 I1
-sg76
+sg79
 S'WW'
-p201
-sg78
+p213
+sg81
 F1.0
-sg79
+sg82
 S'MW'
-p202
-sg81
+p214
+sg85
 S'w-'
-p203
-sg83
-g195
-sg72
-I01
-sg86
+p215
+sg84
+I0
+sg87
+g207
+sg75
+I00
+sg90
 I3
-sg87
-g199
-sg88
+sg91
+g211
+sg92
 I24
-stp204
-Rp205
+stp216
+Rp217
 ag0
-(g68
+(g71
 g2
-(dp206
-g70
-g199
-sg72
-I01
-sg13
-g200
-sg74
-I00
-sg75
-I1
-sg76
-g201
-sg78
-F1.0
-sg79
-g202
-sg81
-g203
-sg83
-g195
-sg85
-I01
-sg86
-I3
-sg87
-g199
-sg88
-I24
-stp207
-Rp208
-atp209
-Rp210
-sg93
-(lp211
-g132
-asg75
-(lp212
-g0
-(g97
-g61
-(lp213
-tp214
-Rp215
-(dp216
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I4
-sg110
-(dp217
-g112
-I1
-ssg113
 (dp218
+g73
+g211
+sg75
+I01
+sg56
+g212
+sg77
+I00
+sg78
+I1
+sg79
+g213
+sg81
+F1.0
+sg82
+g214
+sg84
+I0
+sg85
+g215
+sg87
+g207
+sg89
+I01
+sg90
+I3
+sg91
+g211
+sg92
+I24
+stp219
+Rp220
+atp221
+Rp222
+sg97
+(lp223
+g136
+asg100
+(dp224
+g102
+I1
+ssg78
+(lp225
+g0
+(g104
+g64
+(lp226
+tp227
+Rp228
+(dp229
+g109
+g110
+(g111
+tp230
+Rp231
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp232
 (I0
 I0
-tp219
+tp233
 S'GC_16'
-p220
-sstp221
-Rp222
+p234
+ssg120
+I4
+stp235
+Rp236
 ag0
-(g63
+(g66
 g2
-(dp223
-g65
+(dp237
+g68
 g0
-(g66
-g61
-(lp224
-g205
-ag208
-ag90
-ag90
-atp225
-Rp226
-sg93
-(lp227
+(g69
+g64
+(lp238
+g217
+ag220
+ag94
+ag94
+atp239
+Rp240
+sg97
+(lp241
 S'VVSS1'
-p228
-asg75
-(lp229
+p242
+asg100
+(dp243
+g102
+I2
+ssg78
+(lp244
 g0
-(g97
-g61
-(lp230
-tp231
-Rp232
-(dp233
-g102
-g106
-sg107
+(g104
+g64
+(lp245
+tp246
+Rp247
+(dp248
+g109
+g110
+(g111
+tp249
+Rp250
+sg114
 I00
-sg108
+sg115
 I0
-sbasg109
-I5
-sg110
-(dp234
-g112
-I2
-ssg113
-(dp235
+sbasg116
+(dp251
 (I0
 I0
-tp236
+tp252
 S'GC_10'
-p237
-sstp238
-Rp239
+p253
+ssg120
+I5
+stp254
+Rp255
 ag0
-(g63
+(g66
 g2
-(dp240
-g65
+(dp256
+g68
 g0
-(g66
-g61
-(lp241
-g205
-ag208
-ag90
-atp242
-Rp243
-sg93
-(lp244
+(g69
+g64
+(lp257
+g217
+ag220
+ag94
+atp258
+Rp259
+sg97
+(lp260
 S'VVS1'
-p245
-asg75
-(lp246
-g0
-(g97
-g61
-(lp247
-tp248
-Rp249
-(dp250
+p261
+asg100
+(dp262
 g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I6
-sg110
-(dp251
-g112
 I1
-ssg113
-(dp252
-(I0
-I0
-tp253
-S'GC_22'
-p254
-sstp255
-Rp256
-ag0
-(g63
-g2
-(dp257
-g65
-g0
-(g66
-g61
-(lp258
-g197
-ag197
-ag205
-ag208
-atp259
-Rp260
-sg93
-(lp261
-S'VVVV2'
-p262
-asg75
+ssg78
 (lp263
 g0
-(g97
-g61
+(g104
+g64
 (lp264
 tp265
 Rp266
 (dp267
+g109
+g110
+(g111
+tp268
+Rp269
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp270
+(I0
+I0
+tp271
+S'GC_22'
+p272
+ssg120
+I6
+stp273
+Rp274
+ag0
+(g66
+g2
+(dp275
+g68
+g0
+(g69
+g64
+(lp276
+g209
+ag209
+ag217
+ag220
+atp277
+Rp278
+sg97
+(lp279
+S'VVVV2'
+p280
+asg100
+(dp281
 g102
-g106
-sg107
+I2
+ssg78
+(lp282
+g0
+(g104
+g64
+(lp283
+tp284
+Rp285
+(dp286
+g109
+g110
+(g111
+tp287
+Rp288
+sg114
 I00
-sg108
-I0
-sbasg109
+sg115
+I0
+sbasg116
+(dp289
+(I0
+I0
+tp290
+S'GC_18'
+p291
+ssg120
 I7
-sg110
-(dp268
-g112
-I2
-ssg113
-(dp269
-(I0
-I0
-tp270
-S'GC_18'
-p271
-sstp272
-Rp273
+stp292
+Rp293
 ag0
-(g63
+(g66
 g2
-(dp274
-g65
+(dp294
+g68
 g0
-(g66
-g61
-(lp275
-g205
-ag208
+(g69
+g64
+(lp295
+g217
+ag220
 ag0
-(g68
+(g71
 g2
-(dp276
-g70
+(dp296
+g73
 S'Z'
-p277
-sg72
+p297
+sg75
 I01
-sg13
+sg56
 S'z'
-p278
-sg74
+p298
+sg77
 I01
-sg75
+sg78
 I1
-sg76
+sg79
 S'WZ'
-p279
-sg78
+p299
+sg81
 F0.0
-sg79
+sg82
 S'MZ'
-p280
-sg81
+p300
+sg84
+I0
+sg85
 S'z'
-p281
-sg83
-g195
-sg85
+p301
+sg87
+g207
+sg89
 I01
-sg86
+sg90
 I3
-sg87
-g277
-sg88
+sg91
+g297
+sg92
 I23
-stp282
-Rp283
-atp284
-Rp285
-sg93
-(lp286
-g132
-asg75
-(lp287
-g0
-(g97
-g61
-(lp288
-tp289
-Rp290
-(dp291
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I8
-sg110
-(dp292
-g112
-I1
-ssg113
-(dp293
-(I0
-I0
-tp294
-S'GC_7'
-p295
-sstp296
-Rp297
-ag0
-(g63
-g2
-(dp298
-g65
-g0
-(g66
-g61
-(lp299
-g205
-ag205
-ag208
-ag208
-atp300
-Rp301
-sg93
-(lp302
-g262
-asg75
-(lp303
-g0
-(g97
-g61
-(lp304
-tp305
-Rp306
+stp302
+Rp303
+atp304
+Rp305
+sg97
+(lp306
+g136
+asg100
 (dp307
 g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I9
-sg110
-(dp308
-g112
-I2
-ssg113
-(dp309
-(I0
-I0
+I1
+ssg78
+(lp308
+g0
+(g104
+g64
+(lp309
 tp310
-S'GC_8'
-p311
-sstp312
-Rp313
+Rp311
+(dp312
+g109
+g110
+(g111
+tp313
+Rp314
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp315
+(I0
+I0
+tp316
+S'GC_7'
+p317
+ssg120
+I8
+stp318
+Rp319
 ag0
-(g63
+(g66
 g2
-(dp314
-g65
-g0
-(g66
-g61
-(lp315
-g197
-ag205
-ag208
-ag283
-atp316
-Rp317
-sg93
-(lp318
-S'VVVV5'
-p319
-asg75
-(lp320
-g0
-(g97
-g61
+(dp320
+g68
+g0
+(g69
+g64
 (lp321
-tp322
+g217
+ag217
+ag220
+ag220
+atp322
 Rp323
-(dp324
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I10
-sg110
+sg97
+(lp324
+g280
+asg100
 (dp325
-g112
-I2
-ssg113
-(dp326
-(I0
-I0
-tp327
+g102
+I2
+ssg78
+(lp326
+g0
+(g104
+g64
+(lp327
+tp328
+Rp329
+(dp330
+g109
+g110
+(g111
+tp331
+Rp332
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp333
+(I0
+I0
+tp334
+S'GC_8'
+p335
+ssg120
+I9
+stp336
+Rp337
+ag0
+(g66
+g2
+(dp338
+g68
+g0
+(g69
+g64
+(lp339
+g209
+ag217
+ag220
+ag303
+atp340
+Rp341
+sg97
+(lp342
+S'VVVV5'
+p343
+asg100
+(dp344
+g102
+I2
+ssg78
+(lp345
+g0
+(g104
+g64
+(lp346
+tp347
+Rp348
+(dp349
+g109
+g110
+(g111
+tp350
+Rp351
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp352
+(I0
+I0
+tp353
 S'GC_17'
-p328
-sstp329
-Rp330
-ag0
-(g63
-g2
-(dp331
-g65
-g0
-(g66
-g61
-(lp332
-g283
-ag283
-ag90
-ag90
-atp333
-Rp334
-sg93
-(lp335
-g228
-asg75
-(lp336
-g0
-(g97
-g61
-(lp337
-tp338
-Rp339
-(dp340
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I11
-sg110
-(dp341
-g112
-I2
-ssg113
-(dp342
-(I0
-I0
-tp343
-S'GC_20'
-p344
-sstp345
-Rp346
-ag0
-(g63
-g2
-(dp347
-g65
-g0
-(g66
-g61
-(lp348
-g283
-ag283
-ag90
-atp349
-Rp350
-sg93
-(lp351
-g245
-asg75
-(lp352
-g0
-(g97
-g61
-(lp353
-tp354
-Rp355
-(dp356
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I12
-sg110
+p354
+ssg120
+I10
+stp355
+Rp356
+ag0
+(g66
+g2
 (dp357
-g112
-I1
-ssg113
-(dp358
-(I0
-I0
-tp359
-S'GC_23'
-p360
-sstp361
-Rp362
-ag0
-(g63
-g2
-(dp363
-g65
-g0
-(g66
-g61
+g68
+g0
+(g69
+g64
+(lp358
+g303
+ag303
+ag94
+ag94
+atp359
+Rp360
+sg97
+(lp361
+g242
+asg100
+(dp362
+g102
+I2
+ssg78
+(lp363
+g0
+(g104
+g64
 (lp364
-g205
-ag208
-ag283
-ag283
-atp365
+tp365
 Rp366
-sg93
-(lp367
-g262
-asg75
-(lp368
-g0
-(g97
-g61
-(lp369
-tp370
-Rp371
-(dp372
-g102
-g106
-sg107
+(dp367
+g109
+g110
+(g111
+tp368
+Rp369
+sg114
 I00
-sg108
+sg115
 I0
-sbasg109
-I13
-sg110
-(dp373
-g112
-I2
-ssg113
-(dp374
+sbasg116
+(dp370
 (I0
 I0
-tp375
-S'GC_9'
-p376
-sstp377
+tp371
+S'GC_20'
+p372
+ssg120
+I11
+stp373
+Rp374
+ag0
+(g66
+g2
+(dp375
+g68
+g0
+(g69
+g64
+(lp376
+g303
+ag303
+ag94
+atp377
 Rp378
-ag0
-(g63
-g2
-(dp379
-g65
-g0
-(g66
-g61
-(lp380
-g0
-(g68
-g2
-(dp381
-g70
-S'd'
-p382
-sg85
-I00
-sg13
-S'd'
-p383
-sg74
-I00
-sg75
-I3
-sg76
-g124
-sg78
-F-0.33333333333333331
-sg79
-g124
-sg81
-S'd~'
-p384
-sg83
-S'straight'
-p385
-sg72
-I01
-sg86
-I2
-sg87
-g382
-sg88
+sg97
+(lp379
+g261
+asg100
+(dp380
+g102
 I1
-stp386
+ssg78
+(lp381
+g0
+(g104
+g64
+(lp382
+tp383
+Rp384
+(dp385
+g109
+g110
+(g111
+tp386
 Rp387
-ag0
-(g68
-g2
+sg114
+I00
+sg115
+I0
+sbasg116
 (dp388
-g70
-g382
-sg72
-I01
-sg13
-g383
-sg74
-I00
-sg75
-I3
-sg76
-g124
-sg78
-F-0.33333333333333331
-sg79
-g124
-sg81
-g384
-sg83
-g385
-sg85
-I01
-sg86
-I2
-sg87
-g382
-sg88
-I1
-stp389
-Rp390
-ag197
-atp391
+(I0
+I0
+tp389
+S'GC_23'
+p390
+ssg120
+I12
+stp391
 Rp392
-sg93
-(lp393
-S'FFV1'
-p394
-asg75
-(lp395
+ag0
+(g66
+g2
+(dp393
+g68
 g0
-(g97
-g61
-(lp396
-cmadgraph.core.color_algebra
-T
-p397
-(I1
-I0
-tp398
-Rp399
-atp400
-Rp401
-(dp402
+(g69
+g64
+(lp394
+g217
+ag220
+ag303
+ag303
+atp395
+Rp396
+sg97
+(lp397
+g280
+asg100
+(dp398
 g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I14
-sg110
+I2
+ssg78
+(lp399
+g0
+(g104
+g64
+(lp400
+tp401
+Rp402
 (dp403
-g112
-I1
-ssg113
-(dp404
+g109
+g110
+(g111
+tp404
+Rp405
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp406
 (I0
 I0
-tp405
-S'GC_1'
-p406
-sstp407
-Rp408
+tp407
+S'GC_9'
+p408
+ssg120
+I13
+stp409
+Rp410
 ag0
-(g63
-g2
-(dp409
-g65
-g0
 (g66
-g61
-(lp410
-g0
-(g68
 g2
 (dp411
-g70
-S's'
-p412
-sg85
-I00
-sg13
-S's'
-p413
-sg74
-I00
-sg75
-I3
-sg76
-g124
-sg78
-F-0.33333333333333331
-sg79
-g124
-sg81
-S's~'
+g68
+g0
+(g69
+g64
+(lp412
+g0
+(g71
+g2
+(dp413
+g73
+S'd'
 p414
-sg83
-g385
-sg72
-I01
-sg86
-I2
-sg87
-g412
-sg88
-I3
-stp415
-Rp416
-ag0
-(g68
-g2
-(dp417
-g70
-g412
-sg72
-I01
-sg13
-g413
-sg74
+sg89
+I01
+sg56
+S'd'
+p415
+sg77
 I00
-sg75
-I3
-sg76
-g124
 sg78
-F-0.33333333333333331
+I3
 sg79
-g124
+g128
 sg81
+F-0.3333333333333333
+sg82
+g128
+sg85
+S'd~'
+p416
+sg84
+I0
+sg87
+S'straight'
+p417
+sg75
+I00
+sg90
+I2
+sg91
 g414
-sg83
-g385
-sg85
-I01
-sg86
-I2
-sg87
-g412
-sg88
-I3
+sg92
+I1
 stp418
 Rp419
-ag197
-atp420
-Rp421
-sg93
-(lp422
-g394
-asg75
-(lp423
+ag0
+(g71
+g2
+(dp420
+g73
+g414
+sg75
+I01
+sg56
+g415
+sg77
+I00
+sg78
+I3
+sg79
+g128
+sg81
+F-0.3333333333333333
+sg82
+g128
+sg84
+I0
+sg85
+g416
+sg87
+g417
+sg89
+I01
+sg90
+I2
+sg91
+g414
+sg92
+I1
+stp421
+Rp422
+ag209
+atp423
+Rp424
+sg97
+(lp425
+S'FFV1'
+p426
+asg100
+(dp427
+g102
+I1
+ssg78
+(lp428
 g0
-(g97
-g61
-(lp424
-g397
+(g104
+g64
+(lp429
+cmadgraph.core.color_algebra
+T
+p430
 (I1
 I0
-tp425
-Rp426
-atp427
-Rp428
-(dp429
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I15
-sg110
-(dp430
-g112
-I1
-ssg113
-(dp431
-(I0
-I0
-tp432
-g406
-sstp433
+tp431
+Rp432
+atp433
 Rp434
-ag0
-(g63
-g2
 (dp435
-g65
-g0
-(g66
-g61
-(lp436
-g0
-(g68
-g2
-(dp437
-g70
-S'b'
-p438
-sg85
-I00
-sg13
-S'b'
-p439
-sg74
-I00
-sg75
-I3
-sg76
-g124
-sg78
-F-0.33333333333333331
-sg79
-S'MB'
+g109
+g110
+(g111
+tp436
+Rp437
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp438
+(I0
+I0
+tp439
+S'GC_1'
 p440
-sg81
-S'b~'
-p441
-sg83
-g385
-sg72
-I01
-sg86
-I2
-sg87
-g438
-sg88
-I5
-stp442
-Rp443
-ag0
-(g68
-g2
-(dp444
-g70
-g438
-sg72
-I01
-sg13
-g439
-sg74
-I00
-sg75
-I3
-sg76
-g124
-sg78
-F-0.33333333333333331
-sg79
-g440
-sg81
-g441
-sg83
-g385
-sg85
-I01
-sg86
-I2
-sg87
-g438
-sg88
-I5
-stp445
-Rp446
-ag197
-atp447
-Rp448
-sg93
-(lp449
-g394
-asg75
-(lp450
-g0
-(g97
-g61
-(lp451
-g397
-(I1
-I0
-tp452
+ssg120
+I14
+stp441
+Rp442
+ag0
+(g66
+g2
+(dp443
+g68
+g0
+(g69
+g64
+(lp444
+g0
+(g71
+g2
+(dp445
+g73
+S's'
+p446
+sg89
+I01
+sg56
+S's'
+p447
+sg77
+I00
+sg78
+I3
+sg79
+g128
+sg81
+F-0.3333333333333333
+sg82
+g128
+sg85
+S's~'
+p448
+sg84
+I0
+sg87
+g417
+sg75
+I00
+sg90
+I2
+sg91
+g446
+sg92
+I3
+stp449
+Rp450
+ag0
+(g71
+g2
+(dp451
+g73
+g446
+sg75
+I01
+sg56
+g447
+sg77
+I00
+sg78
+I3
+sg79
+g128
+sg81
+F-0.3333333333333333
+sg82
+g128
+sg84
+I0
+sg85
+g448
+sg87
+g417
+sg89
+I01
+sg90
+I2
+sg91
+g446
+sg92
+I3
+stp452
 Rp453
+ag209
 atp454
 Rp455
-(dp456
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I16
-sg110
-(dp457
-g112
-I1
-ssg113
-(dp458
-(I0
+sg97
+(lp456
+g426
+asg100
+g427
+sg78
+(lp457
+g0
+(g104
+g64
+(lp458
+g430
+(I1
 I0
 tp459
-g406
-sstp460
-Rp461
-ag0
-(g63
-g2
-(dp462
-g65
-g0
-(g66
-g61
-(lp463
-g0
-(g68
-g2
-(dp464
-g70
-S'e-'
-p465
-sg85
-I00
-sg13
-S'e-'
-p466
-sg74
-I00
-sg75
-I1
-sg76
-g124
-sg78
-F-1.0
-sg79
-g124
-sg81
-S'e+'
-p467
-sg83
-g385
-sg72
-I01
-sg86
-I2
-sg87
-g465
-sg88
-I11
+Rp460
+atp461
+Rp462
+(dp463
+g109
+g110
+(g111
+tp464
+Rp465
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp466
+(I0
+I0
+tp467
+g440
+ssg120
+I15
 stp468
 Rp469
 ag0
-(g68
+(g66
 g2
 (dp470
-g70
-g465
-sg72
-I01
-sg13
-g466
-sg74
-I00
-sg75
-I1
-sg76
-g124
-sg78
-F-1.0
-sg79
-g124
-sg81
-g467
-sg83
-g385
-sg85
-I01
-sg86
-I2
-sg87
-g465
-sg88
-I11
-stp471
-Rp472
-ag197
-atp473
-Rp474
-sg93
-(lp475
-g394
-asg75
-(lp476
-g0
-(g97
-g61
-(lp477
-tp478
-Rp479
-(dp480
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I17
-sg110
-(dp481
-g112
-I1
-ssg113
-(dp482
+g68
+g0
+(g69
+g64
+(lp471
+g0
+(g71
+g2
+(dp472
+g73
+S'b'
+p473
+sg89
+I01
+sg56
+S'b'
+p474
+sg77
+I00
+sg78
+I3
+sg79
+g128
+sg81
+F-0.3333333333333333
+sg82
+S'MB'
+p475
+sg85
+S'b~'
+p476
+sg84
+I0
+sg87
+g417
+sg75
+I00
+sg90
+I2
+sg91
+g473
+sg92
+I5
+stp477
+Rp478
+ag0
+(g71
+g2
+(dp479
+g73
+g473
+sg75
+I01
+sg56
+g474
+sg77
+I00
+sg78
+I3
+sg79
+g128
+sg81
+F-0.3333333333333333
+sg82
+g475
+sg84
+I0
+sg85
+g476
+sg87
+g417
+sg89
+I01
+sg90
+I2
+sg91
+g473
+sg92
+I5
+stp480
+Rp481
+ag209
+atp482
+Rp483
+sg97
+(lp484
+g426
+asg100
+g427
+sg78
+(lp485
+g0
+(g104
+g64
+(lp486
+g430
+(I1
+I0
+tp487
+Rp488
+atp489
+Rp490
+(dp491
+g109
+g110
+(g111
+tp492
+Rp493
+sg114
+I00
+sg115
+I0
+sbasg116
+(dp494
 (I0
 I0
-tp483
-S'GC_3'
-p484
-sstp485
-Rp486
-ag0
-(g63
-g2
-(dp487
-g65
-g0
-(g66
-g61
-(lp488
-g0
-(g68
-g2
-(dp489
-g70
-S'm-'
-p490
-sg85
-I00
-sg13
-S'm-'
-p491
-sg74
-I00
-sg75
-I1
-sg76
-g124
-sg78
-F-1.0
-sg79
-g124
-sg81
-S'm+'
-p492
-sg83
-g385
-sg72
-I01
-sg86
-I2
-sg87
-g490
-sg88
-I13
-stp493
-Rp494
-ag0
-(g68
-g2
-(dp495
-g70
-g490
-sg72
-I01
-sg13
-g491
-sg74
-I00
-sg75
-I1
-sg76
-g124
-sg78
-F-1.0
-sg79
-g124
-sg81
-g492
-sg83
-g385
-sg85
-I01
-sg86
-I2
-sg87
-g490
-sg88
-I13
+tp495
+g440
+ssg120
+I16
 stp496
 Rp497
-ag197
-atp498
-Rp499
-sg93
-(lp500
-g394
-asg75
-(lp501
-g0
-(g97
-g61
-(lp502
-tp503
-Rp504
-(dp505
-g102
-g106
-sg107
-I00
-sg108
-I0
-sbasg109
-I18
-sg110
+ag0
+(g66
+g2
+(dp498
+g68
+g0
+(g69
+g64
+(lp499
+g0
+(g71
+g2
+(dp500
+g73
+S'e-'
+p501
+sg89
+I01
+sg56
+S'e-'
+p502
+sg77
+I00
+sg78
+I1
+sg79
+g128
+sg81
+F-1.0
+sg82
+g128
+sg85
+S'e+'
+p503
+sg84
+I1
+sg87
+g417
+sg75
+I00
+sg90
+I2
+sg91
+g501
+sg92
+I11
+stp504
+Rp505
+ag0
+(g71
+g2
 (dp506
-g112
-I1
-ssg113
-(dp507
-(I0
-I0
-tp508
-g484
-sstp509
+g73
+g501
+sg75
+I01
+sg56
+g502
+sg77
+I00
+sg78
+I1
+sg79
+g128
+sg81
+F-1.0
+sg82
+g128
+sg84
+I1
+sg85
+g503