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
