← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to MadSpin/decay.py

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

Show diffs side-by-side

added added

removed removed

Lines of Context:
MadSpin/decay.py
3776
3776
                assert production['total_br'] - min(total_br) < 1e-4
3777
3777
        
3778
3778
        self.branching_ratio = max(total_br) * eff
3779
 
        
3780
 
        
 
3779
 
3781
3780
        #self.banner['madspin'] += ms_banner
3782
3781
        # Update cross-section in the banner
3783
3782
        if 'mggenerationinfo' in self.banner:
3785
3784
            for i,line in enumerate(mg_info):
3786
3785
                if 'Events' in line:
3787
3786
                    if eff == 1:
 
3787
                        self.err_branching_ratio = 0
3788
3788
                        continue
3789
 
                    nb_event =  int(int(mg_info[i].split()[-1]) * eff) 
 
3789
                    initial_event = int(mg_info[i].split()[-1])
 
3790
                    nb_event =  int(initial_event * eff) 
3790
3791
                    mg_info[i] = '#  Number of Events        :       %i' % nb_event
 
3792
                    if eff >0.5:
 
3793
                        self.err_branching_ratio = max(total_br) * math.sqrt(initial_event - eff * initial_event)/initial_event
 
3794
                    else:
 
3795
                        self.err_branching_ratio = max(total_br) * math.sqrt(eff * initial_event)/initial_event
3791
3796
                    continue
3792
3797
                if ':' not in line:
3793
3798
                    continue
3801
3806
                else:
3802
3807
                    mg_info[i] = '%s : %s' % (info, value * self.branching_ratio)
3803
3808
                self.banner['mggenerationinfo'] = '\n'.join(mg_info)
 
3809
                
 
3810
            
 
3811
        
3804
3812
        if 'init' in self.banner:
3805
3813
            new_init =''
3806
3814
            for line in self.banner['init'].split('\n'):