ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/OSUT3Analysis/Configuration/scripts/makePlots.py
(Generate patch)

Comparing UserCode/OSUT3Analysis/Configuration/scripts/makePlots.py (file contents):
Revision 1.53 by wulsin, Mon Jul 15 11:20:00 2013 UTC vs.
Revision 1.54 by lantonel, Tue Jul 16 14:38:50 2013 UTC

# Line 131 | Line 131 | def ratioHistogram( dataHist, mcHist, re
131      def groupErr(group):
132          Data,MC = [float(sum(hist.GetBinContent(i) for i in group)) for hist in [dataHist,mcHist]]
133          dataErr2,mcErr2 = [sum(hist.GetBinError(i)**2 for i in group) for hist in [dataHist,mcHist]]
134 <        return abs(math.sqrt( (dataErr2+mcErr2)/(Data-MC)**2 + mcErr2/MC**2 ) * (Data-MC)/MC) if Data and MC else 0
134 >        if Data > 0 and MC > 0 and Data != MC:
135 >            return abs(math.sqrt( (dataErr2+mcErr2)/(Data-MC)**2 + mcErr2/MC**2 ) * (Data-MC)/MC)
136 >        else:
137 >            return 0
138  
139      def regroup(groups):
140          err,iG = max( (groupErr(g),groups.index(g)) for g in groups )

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines