~maddevelopers/mg5amcnlo/2.7.4_array

                    if name not in self.list_branches or len(self.list_branches[name]) == 0:

                        continue

                    raise self.InvalidCmd("The onshell mode of MadSpin does not support event files where events do not *all* share the same set of final state particles to be decayed.")

        self.branching_ratio = br

        self.efficiency = 1

        self.cross, self.error = self.banner.get_cross(witherror=True)

        #4. determine the maxwgt

        maxwgt = self.get_maxwgt_for_onshell(orig_lhe, evt_decayfile)

        

        #5. generate the decay 

        orig_lhe.seek(0)

        output_lhe = lhe_parser.EventFile(orig_lhe.name.replace('.lhe', '_decayed.lhe'), 'w')

        # event_decay is a dict pdg -> list of event file (contain the decay)

        

        

        nevents = self.options['Nevents_for_max_weigth']

        if nevents == 0 :

            nevents = 75

        logger.info("*****************************")

        logger.info("Probing the first %s events with %s phase space points" % (nevents, self.options['max_weight_ps_point']))

 

        self.efficiency = 1. / self.options['max_weight_ps_point']

        start = time.time()

        for i in range(nevents):

                base_event = base_event[0]

            for j in range(self.options['max_weight_ps_point']):

                decays = self.get_decay_from_file(base_event, evt_decayfile, nevents-i)   

                _, wgt = self.get_onshell_evt_and_wgt(base_event, decays)

                maxwgt = max(wgt, maxwgt)

            all_maxwgt.append(maxwgt)

                

            if all_maxwgt[1] > base_max_weight:

                base_max_weight = 1.05 * all_maxwgt[1]

        return base_max_weight

            

            

            raise

        import copy

        

        if hasattr(production, 'me_wgt'):

            production_me = production.me_wgt

        else:

            production_me = self.calculate_matrix_element(production)

            production.me_wgt = production_me

        decay_me = 1.0

        for pdg in decays:

            for dec in decays[pdg]:

            final = tuple(-i for i in final)

            tag = (init, final)

            orig_order = self.all_me[tag]['order']

        pdir = self.all_me[tag]['pdir']

        if pdir in self.all_f2py:

            p = event.get_momenta(orig_order)

            

            mymod = __import__("%s.matrix2py" % (pdir))

            reload(mymod)

            with misc.chdir(pjoin(self.path_me, 'madspin_me', 'SubProcesses', pdir)):

                with misc.stdchannel_redirected(sys.stdout, os.devnull):

            return self.calculate_matrix_element(event)