Merge lp:~maddevelopers/mg5amcnlo/2.2.0_tutofix into lp:~maddevelopers/mg5amcnlo/2.2.0

[Olivier Mattelaer]

A couple of small change that I would like to have release for tommorow (sorry for that).
the change are
- allow madwidth for NLO model/generation
- have automatic plot for HW6 (with MA4)
- adding the possibility to script the use of showerkt merging
- fixing a bug for DJR plot with root6

[Valentin Hirschi (valentin-hirschi)]

Reviewed via skype.

review: Approve

=== modified file 'Template/LO/bin/internal/plot_tree.C'
--- Template/LO/bin/internal/plot_tree.C	2011-06-29 14:15:21 +0000
+++ Template/LO/bin/internal/plot_tree.C	2014-09-24 19:46:10 +0000
@@ -32,11 +32,11 @@
   Float_t Xsecfact;
   leaf_Xsec->SetAddress(&Xsecfact);
   xsecs->GetEntry(0);
-
-  for(int i=0;i<maxjets && events->GetEntries()>0;i++){
+  if (events->GetEntries()>0) {
+  for(int i=0;i<maxjets;i++){
     events->SetLineWidth(2);
-    events->SetLineColor(color[i]);
-    events->SetLineStyle(style[i]);
+    events->SetLineColor(i+2);
+    events->SetLineStyle(i+2);
     
     if(log) 
       sprintf(tmp1,"log10(%s)>>%s%i%s",quantity,quantity,i,plotdim);
@@ -54,7 +54,7 @@
       return false;
     }
   }
-
+  }
 
   TH1F *hsum = (TH1F*)hists[0]->Clone();
   sprintf(tmp3,"%ssum",quantity);
 
=== modified file 'UpdateNotes.txt'
--- UpdateNotes.txt	2014-09-17 11:50:05 +0000
+++ UpdateNotes.txt	2014-09-24 19:46:10 +0000
@@ -36,10 +36,16 @@
                     set lhc 14    # configure for LHC 14TeV
                     set ilc 1000  # configure for ilc 1TeV
                     set fixed_scale 100 # set all scale to fixed and at 100GeV
+		    set showerkt T # set showerkt on T in the shower card
+ 		    set qcut 20    # set the qctu to 20 in the shower card 
                 OM: Fix a bug in the card edition mode which was sometimes returning to default value
                     which were edited by hand and not via the set command.
                 Seoyoung Kim (by OM): Implementation of the htcaas (super-)cluster support.
 		Juan Rojo (by RF): extended the 3 internal NNPDF sets for scales relevant for a 100TeV collider.
+                OM: Fix a problem with the creation of DJR plot with root 6
+                OM: allow the set the width to Auto in NLO computation (width computated at LO accuracy)
+                OM: Adding the possibility to have automatic plot after the parton shower for Herwig6/Pythia6.
+                    This require MadAnalysis and the pythia-pgs package. 
 
 2.1.2(03/07/14) OM: Fix a bug in ALOHA in presence of customized propagator (Thanks Saurabh)
                 OM: Fixing some compilation issue with MadWeight (Thanks A. Pin)
 
=== modified file 'madgraph/interface/amcatnlo_run_interface.py'
--- madgraph/interface/amcatnlo_run_interface.py	2014-09-17 07:02:56 +0000
+++ madgraph/interface/amcatnlo_run_interface.py	2014-09-24 19:46:10 +0000
@@ -404,8 +404,6 @@
         """Check the argument for the plot command
         plot run_name modes"""
 
-        if options['force']:
-            self.force = True
 
         madir = self.options['madanalysis_path']
         td = self.options['td_path']
@@ -966,7 +964,8 @@
         self.check_plot(args)
         logger.info('plot for run %s' % self.run_name)
         
-        self.ask_edit_cards([], args, plot=True)
+        if not self.force:
+            self.ask_edit_cards([], args, plot=True)
                 
         if any([arg in ['parton'] for arg in args]):
             filename = pjoin(self.me_dir, 'Events', self.run_name, 'events.lhe')
@@ -1024,14 +1023,17 @@
                 misc.gunzip(filename, keep=True, stdout=pjoin(self.me_dir, 'Events','pythia_events.hep'))
                 
                 if not os.path.exists(pjoin(self.me_dir, 'Cards', 'pythia_card.dat')):
-                    files.cp(pjoin(self.me_dir, 'Cards', 'pythia_card_default.dat'),
+                    if aMCatNLO and not self.options['mg5_path']:
+                        raise "plotting NLO HEP file needs MG5 utilities"
+                    
+                    files.cp(pjoin(self.options['mg5_path'], 'Template','LO', 'Cards', 'pythia_card_default.dat'),
                              pjoin(self.me_dir, 'Cards', 'pythia_card.dat'))
                 self.run_hep2lhe()
             else:
                 filename = filenames[0]
                 misc.gunzip(filename, keep=True, stdout=pjoin(self.me_dir, 'Events','pythia_events.hep'))
 
-            self.create_plot('Pythia')
+            self.create_plot('shower')
             lhe_file_name = filename.replace('.hep.gz', '.lhe')
             shutil.move(pjoin(self.me_dir, 'Events','pythia_events.lhe'), 
                         lhe_file_name)
@@ -1125,6 +1127,8 @@
     ############################################################################
     def do_treatcards(self, line, amcatnlo=True):
         """Advanced commands: this is for creating the correct run_card.inc from the nlo format"""
+                #check if no 'Auto' are present in the file
+        self.check_param_card(pjoin(self.me_dir, 'Cards','param_card.dat'))
         return super(aMCatNLOCmd,self).do_treatcards(line, amcatnlo)
     
     ############################################################################
@@ -2545,6 +2549,13 @@
                 raise aMCatNLOError('No file has been generated, an error occurred.'+\
              ' More information in %s' % pjoin(os.getcwd(), 'amcatnlo_run.log'))
 
+            # run the plot creation in a secure way
+            try:
+                self.do_plot('%s -f' % self.run_name)
+            except Exception, error:
+                logger.info("Fail to make the plot. Continue...")
+                pass
+
         elif out_id == 'TOP':
             #copy the topdrawer file(s) back in events
             topfiles = []
 
=== modified file 'madgraph/interface/common_run_interface.py'
--- madgraph/interface/common_run_interface.py	2014-09-17 07:02:56 +0000
+++ madgraph/interface/common_run_interface.py	2014-09-24 19:46:10 +0000
@@ -111,6 +111,14 @@
         logger.info("     will be performed automaticaly during the event generation.")
         logger.info("   -f options: answer all question by default.")
 
+    def help_compute_widths(self):
+        logger.info("syntax: compute_widths Particle [Particles] [--precision=] [--path=Param_card] [--output=PATH]")
+        logger.info("-- Compute the widths for the particles specified.")
+        logger.info("   By default, this takes the current param_card and overwrites it.") 
+        logger.info("   Precision allows to define when to include three/four/... body decays (LO).")
+        logger.info("   If this number is an integer then all N-body decay will be included.")      
+
+
     def help_pythia(self):
         logger.info("syntax: pythia [RUN] [--run_options]")
         logger.info("-- run pythia on RUN (current one by default)")
@@ -168,6 +176,111 @@
             if not args[1].isdigit():
                 raise self.InvalidCmd('timeout values should be a integer')
 
+    def check_compute_widths(self, args):
+        """check that the model is loadable and check that the format is of the
+        type: PART PATH --output=PATH -f --precision=N
+        return the model.
+        """
+        
+        # Check that MG5 directory is present .
+        if MADEVENT and not self.options['mg5_path']:
+            raise self.InvalidCmd, '''The automatic computations of widths requires that MG5 is installed on the system.
+            You can install it and set his path in ./Cards/me5_configuration.txt'''
+        elif MADEVENT:
+            sys.path.append(self.options['mg5_path'])
+        try:
+            import models.model_reader as model_reader
+            import models.import_ufo as import_ufo
+        except ImportError:
+            raise self.ConfigurationError, '''Can\'t load MG5.
+            The variable mg5_path should not be correctly configure.'''
+        
+        ufo_path = pjoin(self.me_dir,'bin','internal', 'ufomodel')
+        # Import model
+        if not MADEVENT:
+            modelname = self.find_model_name()
+            #restrict_file = None
+            #if os.path.exists(pjoin(ufo_path, 'restrict_default.dat')):
+            #    restrict_file = pjoin(ufo_path, 'restrict_default.dat')
+            model = import_ufo.import_model(modelname, decay=True, 
+                        restrict=True)
+            if self.mother and self.mother.options['complex_mass_scheme']:
+                model.change_mass_to_complex_scheme()
+        else:
+            model = import_ufo.import_model(pjoin(self.me_dir,'bin','internal', 'ufomodel'),
+                                        decay=True)
+            #pattern for checking complex mass scheme.
+            has_cms = re.compile(r'''set\s+complex_mass_scheme\s*(True|T|1|true|$|;)''')
+            if has_cms.search(open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')\
+                                   ).read()):
+                model.change_mass_to_complex_scheme()
+   
+            
+#        if not hasattr(model.get('particles')[0], 'partial_widths'):
+#            raise self.InvalidCmd, 'The UFO model does not include partial widths information. Impossible to compute widths automatically'
+            
+        # check if the name are passed to default MG5
+        if '-modelname' in open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')).read():
+            model.pass_particles_name_in_mg_default()        
+        model = model_reader.ModelReader(model)
+        particles_name = dict([(p.get('name'), p.get('pdg_code'))
+                                               for p in model.get('particles')])
+        particles_name.update(dict([(p.get('antiname'), p.get('pdg_code'))
+                                               for p in model.get('particles')]))        
+        
+        output = {'model': model, 'force': False, 'output': None, 
+                  'path':None, 'particles': set(), 'body_decay':4.0025,
+                  'min_br':None, 'precision_channel':0.01}
+        for arg in args:
+            if arg.startswith('--output='):
+                output_path = arg.split('=',1)[1]
+                if not os.path.exists(output_path):
+                    raise self.InvalidCmd, 'Invalid Path for the output. Please retry.'
+                if not os.path.isfile(output_path):
+                    output_path = pjoin(output_path, 'param_card.dat')
+                output['output'] = output_path       
+            elif arg == '-f':
+                output['force'] = True
+            elif os.path.isfile(arg):
+                ftype = self.detect_card_type(arg)
+                if ftype != 'param_card.dat':
+                    raise self.InvalidCmd , '%s is not a valid param_card.' % arg
+                output['path'] = arg
+            elif arg.startswith('--path='):
+                arg = arg.split('=',1)[1]
+                ftype = self.detect_card_type(arg)
+                if ftype != 'param_card.dat':
+                    raise self.InvalidCmd , '%s is not a valid param_card.' % arg
+                output['path'] = arg
+            elif arg.startswith('--'):
+                name, value = arg.split('=',1)
+                try:
+                    value = float(value)
+                except Exception:
+                    raise self.InvalidCmd, '--%s requires integer or a float' % name
+                output[name[2:]] = float(value)                
+            elif arg in particles_name:
+                # should be a particles
+                output['particles'].add(particles_name[arg])
+            elif arg.isdigit() and int(arg) in particles_name.values():
+                output['particles'].add(eval(arg))
+            elif arg == 'all':
+                output['particles'] = set(['all'])
+            else:
+                self.help_compute_widths()
+                raise self.InvalidCmd, '%s is not a valid argument for compute_widths' % arg
+        if self.force:
+            output['force'] = True
+
+        if not output['particles']:
+            raise self.InvalidCmd, '''This routines requires at least one particle in order to compute
+            the related width'''
+            
+        if output['output'] is None:
+            output['output'] = output['path']
+
+        return output
+
     def check_open(self, args):
         """ check the validity of the line """
 
@@ -739,6 +852,12 @@
                 event_path = pjoin(self.me_dir, 'Events', self.run_name,'%s_delphes_events.lhco' % tag)
                 output = pjoin(self.me_dir, 'HTML',self.run_name,
                               'plots_delphes_%s.html' % tag)
+            elif mode == "shower":
+                event_path = pjoin(self.me_dir, 'Events','pythia_events.lhe')
+                output = pjoin(self.me_dir, 'HTML',self.run_name,
+                              'plots_shower_%s.html' % tag)
+                if not self.options['pythia-pgs_path']:
+                    return
             else:
                 raise self.InvalidCmd, 'Invalid mode %s' % mode
         elif mode == 'reweight' and not output:
@@ -815,7 +934,7 @@
 
         # Creating LHE file
         if misc.is_executable(pjoin(pydir, 'hep2lhe')):
-            self.update_status('Creating Pythia LHE File', level='pythia')
+            self.update_status('Creating shower LHE File (for plot)', level='pythia')
             # Write the banner to the LHE file
             out = open(pjoin(self.me_dir,'Events','pythia_events.lhe'), 'w')
             #out.writelines('<LesHouchesEvents version=\"1.0\">\n')
@@ -829,9 +948,10 @@
 
             self.cluster.launch_and_wait(self.dirbin+'/run_hep2lhe',
                                          argument= [pydir],
-                                        cwd=pjoin(self.me_dir,'Events'))
+                                        cwd=pjoin(self.me_dir,'Events'),
+                                        stdout=open(os.devnull,'w'))
 
-            logger.info('Warning! Never use this pythia lhe file for detector studies!')
+            logger.info('Warning! Never use this lhe file for detector studies!')
             # Creating ROOT file
             if eradir and misc.is_executable(pjoin(eradir, 'ExRootLHEFConverter')):
                 self.update_status('Creating Pythia LHE Root File', level='pythia')
@@ -960,6 +1080,33 @@
 
         self.update_status('finish', level='pgs', makehtml=False)
 
+    ############################################################################                                                                                                           
+    def do_compute_widths(self, line):
+        """Require MG5 directory: Compute automatically the widths of a set 
+        of particles"""
+
+        args = self.split_arg(line)
+        opts = self.check_compute_widths(args)
+        
+        
+        from madgraph.interface.master_interface import MasterCmd
+        cmd = MasterCmd()
+        self.define_child_cmd_interface(cmd, interface=False)
+        cmd.exec_cmd('set automatic_html_opening False --no_save')
+        if not opts['path']:
+            opts['path'] = pjoin(self.me_dir, 'Cards', 'param_card.dat')
+            if not opts['force'] :
+                self.ask_edit_cards(['param_card'],[], plot=False)
+        
+        
+        line = 'compute_widths %s %s' % \
+                (' '.join([str(i) for i in opts['particles']]),
+                 ' '.join('--%s=%s' % (key,value) for (key,value) in opts.items()
+                        if key not in ['model', 'force', 'particles'] and value))
+        
+        cmd.exec_cmd(line, model=opts['model'])
+        self.child = None
+        del cmd
 
     ############################################################################ 
     def do_print_results(self, line):
@@ -1313,6 +1460,24 @@
             cluster_name = opt['cluster_type']
             self.cluster = cluster.from_name[cluster_name](**opt)
 
+    def check_param_card(self, path, run=True):
+        """Check that all the width are define in the param_card.
+        If some width are set on 'Auto', call the computation tools."""
+        
+        pattern = re.compile(r'''decay\s+(\+?\-?\d+)\s+auto''',re.I)
+        text = open(path).read()
+        pdg = pattern.findall(text)
+        if pdg:
+            if run:
+                logger.info('Computing the width set on auto in the param_card.dat')
+                self.do_compute_widths('%s %s' % (' '.join(pdg), path))
+            else:
+                logger.info('''Some width are on Auto in the card. 
+    Those will be computed as soon as you have finish the edition of the cards.
+    If you want to force the computation right now and being able to re-edit
+    the cards afterwards, you can type \"compute_wdiths\".''')
+
+
     def add_error_log_in_html(self, errortype=None):
         """If a ME run is currently running add a link in the html output"""
 
@@ -1696,6 +1861,27 @@
                                         os.path.join('.',*[a for a in args \
                                                     if a.endswith(os.path.sep)]))
 
+    def complete_compute_widths(self, text, line, begidx, endidx):
+        "Complete the compute_widths command"
+
+        args = self.split_arg(line[0:begidx])
+        
+        if args[-1] in  ['--path=', '--output=']:
+            completion = {'path': self.path_completion(text)}
+        elif line[begidx-1] == os.path.sep:
+            current_dir = pjoin(*[a for a in args if a.endswith(os.path.sep)])
+            if current_dir.startswith('--path='):
+                current_dir = current_dir[7:]
+            if current_dir.startswith('--output='):
+                current_dir = current_dir[9:]                
+            completion = {'path': self.path_completion(text, current_dir)}
+        else:
+            completion = {}            
+            completion['options'] = self.list_completion(text, 
+                            ['--path=', '--output=', '--min_br=0.\$'
+                             '--precision_channel=0.\$', '--body_decay='])            
+        
+        return self.deal_multiple_categories(completion)
         
 
 # lhapdf-related functions
@@ -2080,7 +2266,7 @@
                           self.list_completion(text, categories)
         
         if 'shortcut' in allowed.keys():
-            possibilities['special values'] = self.list_completion(text, self.special_shortcut.keys()+['qcut'])
+            possibilities['special values'] = self.list_completion(text, self.special_shortcut.keys()+['qcut', 'showerkt'])
 
         if 'run_card' in allowed.keys():
             opts = self.run_set
@@ -2219,6 +2405,19 @@
                     p_card = '%s \n QCUT= %s' % (p_card, args[1])
                 open(pythia_path, 'w').write(p_card)
                 return
+        # Special case for the showerkt value
+        if args[0].lower() == 'showerkt':
+            pythia_path = pjoin(self.me_dir, 'Cards','pythia_card.dat')
+            if os.path.exists(pythia_path):
+                logger.info('add line SHOWERKT = %s in pythia_card.dat' % args[1].upper())
+                p_card = open(pythia_path,'r').read()
+                p_card, n = re.subn('''^\s*SHOWERKT\s*=\s*[default\de\+\-\.]*\s*$''',
+                                    ''' SHOWERKT = %s ''' % args[1].upper(), \
+                                    p_card, flags=(re.M+re.I))
+                if n==0:
+                    p_card = '%s \n SHOWERKT= %s' % (p_card, args[1].upper())
+                open(pythia_path, 'w').write(p_card)
+                return
             
 
         card = '' #store which card need to be modify (for name conflict)
 
=== modified file 'madgraph/interface/madevent_interface.py'
--- madgraph/interface/madevent_interface.py	2014-09-20 07:40:10 +0000
+++ madgraph/interface/madevent_interface.py	2014-09-24 19:46:10 +0000
@@ -418,13 +418,6 @@
         logger.info("   asked in the run_card.")
         self.run_options_help([])
         
-    def help_compute_widths(self):
-        logger.info("syntax: compute_widths Particle [Particles] [--precision=] [--path=Param_card] [--output=PATH]")
-        logger.info("-- Compute the widths for the particles specified.")
-        logger.info("   By default, this takes the current param_card and overwrites it.") 
-        logger.info("   Precision allows to define when to include three/four/... body decays.")
-        logger.info("   If this number is an integer then all N-body decay will be included.")      
-
     def help_store_events(self):
         """ """
         logger.info("syntax: store_events [--run_options]")
@@ -867,112 +860,7 @@
                 raise self.InvalidCmd('refine arguments are suppose to be number')
             
         return True
-    
-    def check_compute_widths(self, args):
-        """check that the model is loadable and check that the format is of the
-        type: PART PATH --output=PATH -f --precision=N
-        return the model.
-        """
-        
-        # Check that MG5 directory is present .
-        if MADEVENT and not self.options['mg5_path']:
-            raise self.InvalidCmd, '''The automatic computations of widths requires that MG5 is installed on the system.
-            You can install it and set his path in ./Cards/me5_configuration.txt'''
-        elif MADEVENT:
-            sys.path.append(self.options['mg5_path'])
-        try:
-            import models.model_reader as model_reader
-            import models.import_ufo as import_ufo
-        except ImportError:
-            raise self.ConfigurationError, '''Can\'t load MG5.
-            The variable mg5_path should not be correctly configure.'''
-        
-        ufo_path = pjoin(self.me_dir,'bin','internal', 'ufomodel')
-        # Import model
-        if not MADEVENT:
-            modelname = self.find_model_name()
-            #restrict_file = None
-            #if os.path.exists(pjoin(ufo_path, 'restrict_default.dat')):
-            #    restrict_file = pjoin(ufo_path, 'restrict_default.dat')
-            model = import_ufo.import_model(modelname, decay=True, 
-                        restrict=True)
-            if self.mother and self.mother.options['complex_mass_scheme']:
-                model.change_mass_to_complex_scheme()
-        else:
-            model = import_ufo.import_model(pjoin(self.me_dir,'bin','internal', 'ufomodel'),
-                                        decay=True)
-            #pattern for checking complex mass scheme.
-            has_cms = re.compile(r'''set\s+complex_mass_scheme\s*(True|T|1|true|$|;)''')
-            if has_cms.search(open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')\
-                                   ).read()):
-                model.change_mass_to_complex_scheme()
-   
-            
-#        if not hasattr(model.get('particles')[0], 'partial_widths'):
-#            raise self.InvalidCmd, 'The UFO model does not include partial widths information. Impossible to compute widths automatically'
-            
-        # check if the name are passed to default MG5
-        if '-modelname' in open(pjoin(self.me_dir,'Cards','proc_card_mg5.dat')).read():
-            model.pass_particles_name_in_mg_default()        
-        model = model_reader.ModelReader(model)
-        particles_name = dict([(p.get('name'), p.get('pdg_code'))
-                                               for p in model.get('particles')])
-        particles_name.update(dict([(p.get('antiname'), p.get('pdg_code'))
-                                               for p in model.get('particles')]))        
-        
-        output = {'model': model, 'force': False, 'output': None, 
-                  'path':None, 'particles': set(), 'body_decay':4.0025,
-                  'min_br':None, 'precision_channel':0.01}
-        for arg in args:
-            if arg.startswith('--output='):
-                output_path = arg.split('=',1)[1]
-                if not os.path.exists(output_path):
-                    raise self.InvalidCmd, 'Invalid Path for the output. Please retry.'
-                if not os.path.isfile(output_path):
-                    output_path = pjoin(output_path, 'param_card.dat')
-                output['output'] = output_path       
-            elif arg == '-f':
-                output['force'] = True
-            elif os.path.isfile(arg):
-                type = self.detect_card_type(arg)
-                if type != 'param_card.dat':
-                    raise self.InvalidCmd , '%s is not a valid param_card.' % arg
-                output['path'] = arg
-            elif arg.startswith('--path='):
-                arg = arg.split('=',1)[1]
-                type = self.detect_card_type(arg)
-                if type != 'param_card.dat':
-                    raise self.InvalidCmd , '%s is not a valid param_card.' % arg
-                output['path'] = arg
-            elif arg.startswith('--'):
-                name, value = arg.split('=',1)
-                try:
-                    value = float(value)
-                except Exception:
-                    raise self.InvalidCmd, '--%s requires integer or a float' % name
-                output[name[2:]] = float(value)                
-            elif arg in particles_name:
-                # should be a particles
-                output['particles'].add(particles_name[arg])
-            elif arg.isdigit() and int(arg) in particles_name.values():
-                output['particles'].add(eval(arg))
-            elif arg == 'all':
-                output['particles'] = set(['all'])
-            else:
-                self.help_compute_widths()
-                raise self.InvalidCmd, '%s is not a valid argument for compute_widths' % arg
-        if self.force:
-            output['force'] = True
-
-        if not output['particles']:
-            raise self.InvalidCmd, '''This routines requires at least one particle in order to compute
-            the related width'''
-            
-        if output['output'] is None:
-            output['output'] = output['path']
-
-        return output
-    
+        
     def check_combine_events(self, arg):
         """ Check the argument for the combine events command """
         
@@ -1506,28 +1394,6 @@
         else:
             return self.complete_generate_events(*args, **opts)
 
-    def complete_compute_widths(self, text, line, begidx, endidx):
-        "Complete the compute_widths command"
-
-        args = self.split_arg(line[0:begidx])
-        
-        if args[-1] in  ['--path=', '--output=']:
-            completion = {'path': self.path_completion(text)}
-        elif line[begidx-1] == os.path.sep:
-            current_dir = pjoin(*[a for a in args if a.endswith(os.path.sep)])
-            if current_dir.startswith('--path='):
-                current_dir = current_dir[7:]
-            if current_dir.startswith('--output='):
-                current_dir = current_dir[9:]                
-            completion = {'path': self.path_completion(text, current_dir)}
-        else:
-            completion = {}            
-            completion['options'] = self.list_completion(text, 
-                            ['--path=', '--output=', '--min_br=0.\$'
-                             '--precision_channel=0.\$', '--body_decay='])            
-        
-        return self.deal_multiple_categories(completion)
-
     def complete_calculate_decay_widths(self, text, line, begidx, endidx):
         """ Complete the calculate_decay_widths command"""
         
@@ -2694,33 +2560,7 @@
                 self.create_root_file(output='%s/unweighted_events.root' % \
                                                                   self.run_name)
     
-    ############################################################################                                                                                                           
-    def do_compute_widths(self, line):
-        """Require MG5 directory: Compute automatically the widths of a set 
-        of particles"""
 
-        args = self.split_arg(line)
-        opts = self.check_compute_widths(args)
-        
-        
-        from madgraph.interface.master_interface import MasterCmd
-        cmd = MasterCmd()
-        self.define_child_cmd_interface(cmd, interface=False)
-        cmd.exec_cmd('set automatic_html_opening False --no_save')
-        if not opts['path']:
-            opts['path'] = pjoin(self.me_dir, 'Cards', 'param_card.dat')
-            if not opts['force'] :
-                self.ask_edit_cards(['param_card'],[], plot=False)
-        
-        
-        line = 'compute_widths %s %s' % \
-                (' '.join([str(i) for i in opts['particles']]),
-                 ' '.join('--%s=%s' % (key,value) for (key,value) in opts.items()
-                        if key not in ['model', 'force', 'particles'] and value))
-        
-        cmd.exec_cmd(line, model=opts['model'])
-        self.child = None
-        del cmd
 
 
     
@@ -4048,22 +3888,8 @@
                 
   
             
-    def check_param_card(self, path, run=True):
-        """Check that all the width are define in the param_card.
-        If some width are set on 'Auto', call the computation tools."""
-        
-        pattern = re.compile(r'''decay\s+(\+?\-?\d+)\s+auto''',re.I)
-        text = open(path).read()
-        pdg = pattern.findall(text)
-        if pdg:
-            if run:
-                logger.info('Computing the width set on auto in the param_card.dat')
-                self.do_compute_widths('%s %s' % (' '.join(pdg), path))
-            else:
-                logger.info('''Some width are on Auto in the card. 
-    Those will be computed as soon as you have finish the edition of the cards.
-    If you want to force the computation right now and being able to re-edit
-    the cards afterwards, you can type \"compute_wdiths\".''')
+
+                
 #===============================================================================
 # MadEventCmd
 #===============================================================================
 
=== modified file 'madgraph/various/gen_crossxhtml.py'
--- madgraph/various/gen_crossxhtml.py	2014-09-22 07:17:43 +0000
+++ madgraph/various/gen_crossxhtml.py	2014-09-24 19:46:10 +0000
@@ -791,6 +791,10 @@
                glob.glob(pjoin(path,"*.hep.gz")):
                 self.shower.append('hep')
 
+            if 'plot' not in self.shower and \
+                          exists(pjoin(html_path,"plots_shower_%s.html" % tag)):
+                self.shower.append('plot')                
+
             if glob.glob(pjoin(path,"*.hepmc")) + \
                glob.glob(pjoin(path,"*.hepmc.gz")):
                 self.shower.append('hepmc')
@@ -988,6 +992,8 @@
                       glob.glob(pjoin(self.me_dir, 'Events', self['run_name'], '*.' + kind)) + \
                       glob.glob(pjoin(self.me_dir, 'Events', self['run_name'], '*.' + kind + '.gz')):
                         out += " <a href=\"%s\">%s</a> " % (f, kind.upper())
+            if 'plot' in self.shower:
+                out += """ <a href="./HTML/%(run_name)s/plots_shower_%(tag)s.html">plots</a>"""
             
             return out % self
                 
 
=== modified file 'tests/test_manager.py'
--- tests/test_manager.py	2014-05-15 00:02:06 +0000
+++ tests/test_manager.py	2014-09-24 19:46:10 +0000
@@ -115,6 +115,38 @@
             self.stream.writeln(" ".join(self.bypassed))
         return result 
 
+    def run_border(self, test):
+        "Run the given test case or test suite."
+        MyTextTestRunner.stream = self.stream
+        result = self._makeResult()
+        startTime = time.time()
+        test(result)
+        stopTime = time.time()
+        timeTaken = float(stopTime - startTime)
+        result.printErrors()
+        self.stream.writeln(result.separator2)
+        run = result.testsRun
+        #self.stream.writeln("Ran %d test%s in %.3fs" %
+        #                    (run, run != 1 and "s" or "", timeTaken))
+        #self.stream.writeln()
+        if not result.wasSuccessful():
+            self.stream.write("FAILED (")
+            failed, errored = map(len, (result.failures, result.errors))
+            if failed:
+                self.stream.write("failures=%d" % failed)
+            if errored:
+                if failed: self.stream.write(", ")
+                self.stream.write("errors=%d" % errored)
+            self.stream.writeln(")")
+            sys.exit(0)
+        #else:
+        #    self.stream.writeln("OK")
+        #if self.bypassed:
+        #    self.stream.writeln("Bypassed %s:" % len(self.bypassed))
+        #    self.stream.writeln(" ".join(self.bypassed))
+        return result 
+
+
             
 #===============================================================================
 # run
@@ -153,6 +185,46 @@
     #return out
 
 #===============================================================================
+# run
+#===============================================================================
+def run_border_search(to_crash='',expression='', re_opt=0, package='./tests/unit_tests', verbosity=1,
+        timelimit=0):
+    """ running the test associated to expression one by one. and follow them by the to_crash one
+        up to the time that to_crash is actually crashing. Then the run stops and print the list of the 
+        routine tested. Then the code re-run itself(via a fork) to restrict the list. 
+        The code stops when the list is of order 1. The order of the test is randomize at each level!
+    """
+    #init a test suite
+    collect = unittest.TestLoader()
+    TestSuiteModified.time_limit =  float(timelimit)
+    all_test = TestFinder(package=package, expression=expression, re_opt=re_opt)
+    import random
+    random.shuffle(all_test)
+    print "to_crash"
+    to_crash = TestFinder(package=package, expression=to_crash, re_opt=re_opt)
+    to_crash.collect_dir(package, checking=True)
+    print dir(to_crash)
+
+    for test_fct in all_test:
+        testsuite = unittest.TestSuite()
+        data = collect.loadTestsFromName(test_fct)        
+        assert(isinstance(data,unittest.TestSuite))        
+        data.__class__ = TestSuiteModified
+        testsuite.addTest(data)
+        data = collect.loadTestsFromName(to_crash[0])        
+        assert(isinstance(data,unittest.TestSuite))        
+        data.__class__ = TestSuiteModified
+        testsuite.addTest(data)
+        # Running it
+        print "run it for %s" % test_fct
+        output =  MyTextTestRunner(verbosity=verbosity).run_border(testsuite)
+    
+    return output
+    #import tests
+    #print 'runned %s checks' % tests.NBTEST
+    #return out
+
+#===============================================================================
 # listIOTests
 #===============================================================================
 
@@ -262,8 +334,8 @@
     IOTestManager.testFolders_filter = arg.split('/')[0].split('&')
     IOTestManager.testNames_filter = arg.split('/')[1].split('&')
     IOTestManager.filesChecked_filter = '/'.join(arg.split('/')[2:]).split('&')
-    #print "INFO:: Using folders %s"%str(IOTestManager.testFolders_filter)
-    #print "INFO:: Using test names %s"%str(IOTestManager.testNames_filter)
+    #print "INFO:: Using folders %s"%str(IOTestManager.testFolders_filter)    
+    #print "INFO:: Using test names %s"%str(IOTestManager.testNames_filter)         
     #print "INFO:: Using file paths %s"%str(IOTestManager.filesChecked_filter)
     
     # Initiate all the IOTests from all the setUp()
@@ -820,7 +892,10 @@
                                       "content of the folder IOTestsComparison")
     parser.add_option("-t", "--timed", default="Auto",
           help="limit the duration of each test. Negative number re-writes the information file.")    
-    
+
+    parser.add_option("", "--border_effect", default=None,
+          help="Define the test which are sensitive to a border effect, the test will find which test creates this border effect")        
+
     (options, args) = parser.parse_args()
 
     if options.IOTestsUpdate:
@@ -897,9 +972,13 @@
           " MadGraph5_aMCatNLO is configured not to compress the references files."
     
     if options.IOTests=='No' and not options.synchronize:
-        #logging.basicConfig(level=vars(logging)[options.logging])
-        run(args, re_opt=options.reopt, verbosity=options.verbose, \
-            package=options.path, timelimit=options.timed)
+        if not options.border_effect:
+            #logging.basicConfig(level=vars(logging)[options.logging])
+            run(args, re_opt=options.reopt, verbosity=options.verbose, \
+                package=options.path, timelimit=options.timed)
+        else:
+            run_border_search(options.border_effect, args, re_opt=options.reopt, verbosity=options.verbose, \
+                package=options.path, timelimit=options.timed)
     else:
         if options.IOTests=='L':
             print "Listing all tests defined in the reference files ..."