def reportCache(self) : likes = self.likelihoods() fileName = '/'.join(self.outputFileName.split('/')[:-1]+[self.name]) optstat = r.gStyle.GetOptStat() r.gStyle.SetOptStat(0) c = r.TCanvas() utils.tCanvasPrintPdf(c, fileName, False, '[') with open(fileName+'.txt', "w") as file : print >> file, '\n'.join(str(d) for d in self.populations) print >> file for key in sorted((set(likes[0]) & set(likes[1]) & set(likes[2])), key = lambda k: -1 if self.name==k else k) : if not all(item[key] for item in likes) : continue dilutions = [utils.dilution(bins1,bins2) for bins1,bins2 in itertools.combinations([utils.binValues(item[key]) for item in likes],2)] print >> file, "\t".join([key.ljust(25)]+["%.4f"%d for d in dilutions]) h = 1.1*max(item[key].GetMaximum() for item in likes) l = r.TLegend(0.75,0.85,0.98,0.98) l.SetHeader("Dilutions : %.3f %.3f %.3f"%tuple(dilutions)) for i,(item,color,pop) in enumerate(zip(likes,[r.kBlack,r.kRed,r.kBlue],self.populations)) : item[key].SetLineColor(color) item[key].SetMaximum(h) item[key].SetMinimum(0) item[key].GetYaxis().SetTitle("pdf") item[key].Draw("hist" + ('same' if i else '')) l.AddEntry(item[key], pop['pre'],'l') l.Draw() utils.tCanvasPrintPdf(c, fileName, False) c.Clear() c.Divide(3,2) covs = [self.fromCache([item['pre']], ["_cov_%s_%s"%pair for pair in itertools.combinations(self.dists.keys(),2)], tag = item['tag'])[item['pre']] for item in self.populations] for key in (set(covs[0]) & set(covs[1]) & set(covs[2])) : if not all(item[key] for item in covs) : continue hists = [item[key] for item in covs] deps = [utils.dependence(h,limit=1,inSigma=False) for h in hists] for i,(h,d,pop) in enumerate(zip(hists,deps,self.populations)) : d.SetTitle(pop['pre']) c.cd(i+1); h.Draw('colz') c.cd(i+4); d.Draw('colz') utils.tCanvasPrintPdf(c, fileName, False) c.Clear() for key in sorted((set(likes[0]) & set(likes[1]) & set(likes[2])), key = lambda k: -1 if self.name==k else k) : if not all(item[key] for item in likes) : continue for i,j in itertools.combinations(range(3),2) : k = next(k for k in range(3) if k not in [i,j]) templates = [ utils.binValues(likes[ii][key]) for ii in [i,j] ] observed = utils.binValues(likes[k][key]) try: cs = utils.fractions.componentSolver(observed,templates) except: cs = None if not cs : continue notK = likes[k][key].Clone('bf_not_'+likes[k][key].GetName()) notK.Reset() notK.Add(likes[i][key],likes[j][key],cs.fractions[0],cs.fractions[1]) dil = utils.dilution(utils.binValues(likes[k][key]),utils.binValues(notK)) l = r.TLegend(0.75,0.85,0.98,0.98) l.SetHeader("Dilution : %.3f"%dil) l.AddEntry(likes[k][key], self.populations[k]['pre'], 'l') l.AddEntry(notK, 'not '+self.populations[k]['pre'], 'l') notK.SetLineStyle(2) likes[k][key].Draw('hist') notK.Draw('hist same') l.Draw() utils.tCanvasPrintPdf(c, fileName, False) utils.tCanvasPrintPdf(c, fileName, True, ']') r.gStyle.SetOptStat(optstat)