def plot_chi2_scan(self, sizeX=100, sizeY=100, edge=None, toyMC=False, bootstrap=False, simulation=None): ### Read chi2 scan bestFit = self.read_chi2_scan(sizeX, sizeY) if edge is not None: self._chi2_scan_edge = edge if toyMC: bestFit = numpy.append(bestFit, self.read_toyMC(), axis=1) if bootstrap: bestFit = numpy.append(bestFit, self.read_bootStrap(), axis=1) if simulation is not None: bestFit = numpy.append(bestFit, simulation, axis=1) ### Plot the data myTools.plot_chi2_scan( self._chi2_scan, self._chi2_scan_edge, True, bestFit, "\\alpha_{\\perp}", "\\alpha_{\\parallel}", "\\Delta \\chi^{2} = \\chi^{2}-\\chi^{2}_{best \\, fit}", ) return
def plot_chi2_scan(self, sizeX=100, sizeY=100, edge=None, bestFitData=None): self._chi2_scan = numpy.zeros( shape=(sizeX,sizeY) ) nb = numpy.zeros( shape=(sizeX,sizeY) ) if (bestFitData is None): bestFit = numpy.zeros( shape=(2,1) ) else: bestFit = bestFitData for el in self._listFit: bestFit = numpy.append(bestFit,el.get_best_fit(),axis=1) el.read_chi2_scan(sizeX, sizeY) cut = numpy.isfinite(el._chi2_scan) self._chi2_scan[cut] += el._chi2_scan[cut] nb[cut] += 1. self._chi2_scan /= nb ### Plot the data myTools.plot_chi2_scan(self._chi2_scan,self._listFit[0]._chi2_scan_edge,True, bestFit,'\\alpha_{\\perp}','\\alpha_{\\parallel}','\\Delta \\chi^{2} = \\chi^{2}-\\chi^{2}_{best \\, fit}') return