def grid_getmatch_423to404(self, ratio, eratio): match,indbest,chi2r = grid_fitter.grid_getmatch(ratio, eratio, self.modelratio_423_404) return chi2r
def grid_getmatch_322to321(self, ratio, eratio): match,indbest,chi2r = grid_fitter.grid_getmatch(ratio, eratio, self.modelratio2) return chi2r
def grid_getmatch_line(self, tline, etline, transition): match,indbest,chi2r = grid_fitter.grid_getmatch(tline, etline, self.tline[transition]) return chi2r
def grid_getmatch_404to303(self, ratio, eratio): match,indbest,chi2r = grid_fitter.grid_getmatch(ratio, eratio, self.modelratio_404_303) return chi2r
def grid_getmatch_322to321(self, ratio, eratio, chi2_thresh=1): match, indbest, chi2r = grid_fitter.grid_getmatch( ratio, eratio, self.modelratio2, chi2_thresh=chi2_thresh) return chi2r
def grid_getmatch_322(self, taline322, etaline322, chi2_thresh=1): match, indbest, chi2r = grid_fitter.grid_getmatch( taline322, etaline322, self.tline322, chi2_thresh=chi2_thresh) return chi2r
def grid_getmatch_322to321(self, ratio, eratio, chi2_thresh=1): match,indbest,chi2r = grid_fitter.grid_getmatch(ratio, eratio, self.modelratio2, chi2_thresh=chi2_thresh) return chi2r
def grid_getmatch_322to321(self, ratio, eratio): match, indbest, chi2r = grid_fitter.grid_getmatch( ratio, eratio, self.modelratio2) return chi2r
def grid_getmatch_423to404(self, ratio, eratio): match, indbest, chi2r = grid_fitter.grid_getmatch( ratio, eratio, self.modelratio_423_404) return chi2r
def grid_getmatch_404to303(self, ratio, eratio): match, indbest, chi2r = grid_fitter.grid_getmatch( ratio, eratio, self.modelratio_404_303) return chi2r
xx,yy,zz = (90-1, 110-1, 110-1) # low dens, high *column* std but low *dens* std cont1 = cont11[yy,xx] cont2 = cont22[yy,xx] tline1 = cube1.cube[zz,yy,xx] tline2 = cube2.cube[zz,yy,xx] etline1 = cube1.errorcube[zz,yy,xx] etline2 = cube2.errorcube[zz,yy,xx] pargrid1 = (cont1*np.exp(-taugrid1) + (1-np.exp(-taugrid1))*texgrid1) pargrid2 = (cont2*np.exp(-taugrid2) + (1-np.exp(-taugrid2))*texgrid2) #spec = (1.0-np.exp(-np.array(tau_nu_cumul)))*(Tex-Tbackground) match, indbest, chi2 = grid_2p_getmatch(tline1+cont1, etline1, pargrid1, tline2+cont2, etline2, pargrid2) match11, indbest11, chi2_11 = grid_getmatch(tline1+cont1, etline1, pargrid1) match22, indbest22, chi2_22 = grid_getmatch(tline2+cont2, etline2, pargrid2) bw,bz,by,bx = np.unravel_index(indbest, pargrid1.shape) pl.jet() pl.clf() extent = [densityarr.flat[0],densityarr.flat[-1], columnarr.flat[0],columnarr.flat[-1]] pl.imshow(pargrid1[bw,bz,:,:], alpha=0.5, extent=extent) pl.contour(densityarr[bw,bz,:,:],columnarr[bw,bz,:,:], pargrid2[bw,bz,:,:], levels=np.linspace(tline2-20*etline2+cont2, tline2+20*etline2+cont2, 10))