# pretty-print the result
from B2DXFitters.FitResult import getDsHBlindFitResult
result = getDsHBlindFitResult(fitconfig['IsData'], fitconfig['Blinding'],
rawfitresult)
print 'RESULT',result
rawfitresult.floatParsFinal().Print()
ds.get().Print();
'''
for i in range(mistagresultList.GetSize()):
print "mistag =",mistagresultList.At(i).getValV(),"+-",mistagresultList.At(i).getError(), ", avg. eta =",etaAvgList.At(i).getValV(),"+-",etaAvgList.At(i).getError()
'''
print 'GENETALIST'
genEtaList.Print();
print 'ETAAVG'
etaAvg.Print();
print 'ETAAVGLIST'
etaAvgList.Print();
from ROOT import TFile
g = TFile('fitresultlist123a/fitresultlist123a_%04d.root' % SEED, 'recreate')
#g.WriteTObject(p0End, 'fitresultlist123/fitresultlist123003_%04d' % SEED)
#g.WriteTObject(p1End, 'fitresultlist123/fitresultlist123003_%04d' % SEED)
g.WriteTObject(genEtaList,'fitresultlist123a/fitresultlist123a003_%04d' % SEED)
g.WriteTObject(etaAvg, 'fitresultlist123a/fitresultlist123a003_%04d' % SEED)
g.WriteTObject(etaAvgList, 'fitresultlist123a/fitresultlist123a003_%04d' % SEED)
g.Close()
#parse etaHist directory for .root files
for i in fileList[:1]:
if (i[-5:] != '.root'):
continue
print i
inFile = TFile(i)
inFile.GetListOfKeys().Print()
inFile.ls()
#s = raw_input("Press Enter to continue");
theTH1DHist = TH1F()
inFile.GetObject("etaHist", theTH1DHist)
mainResultList.AddLast(
divideEtaTH1Ds(theTH1DHist, numTagCat, i, None, True))
from ROOT import TFile
import time
os.chdir(os.environ['B2DXFITTERSROOT'] + '/tutorial')
if (not (os.path.isdir('fits'))):
os.mkdir('fits')
os.chdir('fits')
mainResultList.Print("", 5)
f = TFile('tagRegionFitList_%f.root' % time.time(), 'recreate')
f.WriteTObject(mainResultList, 'tagRegionFitList_%f' % time.time())
f.Close()
del f
