if len(bn835) != len(old_combinations):
print "len not matching: ", len(bn835), len(old_combinations)
sys.exit(1)
for (c1,c2),(v1,e1,v2,e2) in zip(old_combinations, bn835):
n1 = all_channels[c1]
n2 = all_channels[c2]
print "%10s:%10s: %.4f %.4f" %(n1, n2, v1, v2)
old_eff_svd1[c1,c2] = v1
old_eff_svd2[c1,c2] = v2
old_err_svd1[c1,c2] = e1
old_err_svd2[c1,c2] = e2
b0 = pr.chain(filenames,"B0")
b0g = pr.chain(filenames, "B0gen")
h_gen_svd1 = b0g.draw("channelP:channelM",cut="svdVs==0", option="goff", range=(4,6,10,4,6,10))
h_gen_svd2 = b0g.draw("channelP:channelM",cut="svdVs==1", option="goff", range=(4,6,10,4,6,10))
h_rec_svd1 = b0.draw("bestLHsig.channelP:bestLHsig.channelM",
cut="svdVs==0 && bestLHsig.mcInfo&1 && bestLHsig.flag==0 && bestLHsig.tag.flavour!=0", option="goff", range=(4,6,10,4,6,10))
h_rec_svd2 = b0.draw("bestLHsig.channelP:bestLHsig.channelM",
cut="svdVs==1 && bestLHsig.mcInfo&1 && bestLHsig.flag==0 && bestLHsig.tag.flavour!=0", option="goff", range=(4,6,10,4,6,10))
gen_svd1, dummy = r2mpl.getHistData(h_gen_svd1)
gen_svd2, dummy = r2mpl.getHistData(h_gen_svd2)
rec_svd1, dummy = r2mpl.getHistData(h_rec_svd1)
rec_svd2, dummy = r2mpl.getHistData(h_rec_svd2)
import sys
import utils
import pyroot as pr
b0 = pr.chain(sys.argv[1:],"B0")
b0.draw("(bestLHsig.tag.z - mcVtxTag)*1e4", range=(200,-1e3,1e3), cut=utils.cuts)
sys.stdin.readline()
