#!/usr/bin/python
import os,sys,fpformat
from ROOT import *
fn = 'sum84-9-1_sum584-1_Tam.root'        #data-sim comparison
#fn = 'sum583-1_sum542-1-1_Tam.root'       #sim with and without DS beam package
f = TFile(fn)

def integrate_histo(h):
    print 'Integrating',h.GetName()
    nx = h.GetNbinsX()
    l = []
    for ix in range(1,nx+1):
        bc = h.GetBinCenter(ix)
        n = h.GetBinContent(ix)
#        print ' ',ix,fpformat.fix(bc,3).ljust(7),str(n).ljust(20)
        l.append(n)
    return sum(l)

for hs in ['dataminusmcUS','dataminusmc']:
    print integrate_histo(f.Get(hs))

#print integrate_histo(f.Get('Data_wintype8plus15_Norm'))
#print integrate_histo(f.Get('MC_wintype8plus15_Norm'))

#Data_wintype8plus15_Norm.Draw()
au = raw_input('>')

##c = TCanvas()
##c.Divide(2,1)
##c.cd(1)
##US.Draw()
##c.cd(2)
##DS.Draw()
##
##US_sum = []
##DS_sum = []
##
##nx = dataminusmc.GetNbinsX()
##for ix in range(nx):
##    bc = US.GetBinCenter(ix)
##    US_n = US.GetBinContent(ix)
##    DS_n = DS.GetBinContent(ix)
##    print ix,str(bc).ljust(20),str(US_n).ljust(20),str(DS_n).ljust(20)
##    US_sum.append(US_n)
##    DS_sum.append(DS_n)
##
##print 'US_sum=',sum(US_sum)
##print 'DS_sum=',sum(DS_sum)
##
##au = raw_input('>')
##