# 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