def bootstrapHD103060(): global pup HD10360corr = [] HD10360 = [] for k, f in enumerate(files_hd10360): print k+1, '/', len(files_hd10360) HD10360corr.append(astromNew.bootstrapListOfFiles(f, data_directory, firstGuess={'SEP':8., 'PA':135, 'M0':-0.03, 'AT4 a_2':pup['BEST']['AT4 a_2'], 'AT4 phi_2':pup['BEST']['AT4 phi_2']}, doNotFit=['AT4 phi_2', 'AT4 a_2'], N=100)) HD10360.append(astromNew.bootstrapListOfFiles(f, data_directory, N=100, firstGuess={'SEP':8., 'PA':135, 'M0':-0.03})) pyplot.figure(10) pyplot.clf() pyplot.axes().set_aspect('equal', 'datalim') th = np.linspace(0, 2*np.pi, 100) colors=['r', 'g', 'b', 'y', 'm', 'c'] for k,d in enumerate(HD10360corr): ### observation: pyplot.plot(d['RA_mas']/1000., d['DEC_mas']/1000., '+'+colors[k], markersize=12) pyplot.plot(d['RA_mas']/1000.+ d['errs_uas'][0]*d['PCA'][0][1]*np.cos(th)/1e6+ d['errs_uas'][1]*d['PCA'][0][0]*np.sin(th)/1e6, d['DEC_mas']/1000.+ -d['errs_uas'][1]*d['PCA'][0][1]*np.sin(th)/1e6+ d['errs_uas'][0]*d['PCA'][0][0]*np.cos(th)/1e6, color=colors[k], linewidth=3, alpha=0.5) for k,d in enumerate(HD10360): ### observation: pyplot.plot(d['RA_mas']/1000., d['DEC_mas']/1000., '.'+colors[k], markersize=8) pyplot.plot(d['RA_mas']/1000.+ d['errs_uas'][0]*d['PCA'][0][1]*np.cos(th)/1e6+ d['errs_uas'][1]*d['PCA'][0][0]*np.sin(th)/1e6, d['DEC_mas']/1000.+ -d['errs_uas'][1]*d['PCA'][0][1]*np.sin(th)/1e6+ d['errs_uas'][0]*d['PCA'][0][0]*np.cos(th)/1e6, color=colors[k], linewidth=3, alpha=0.5) return
def bootstrapHD10360(N=100, ATPUP_scale=1.0, plot=False): """ bootstraping. ATPUP_scale=0.0 sets the pupil correction to 0. N is the number of bootstrapping iteration. """ t0 = time.time() fgNACO['ATPUP scale'] = ATPUP_scale b = astromNew.bootstrapListOfFiles(files_hd10360_NACO, directory=data_directory, firstGuess=fgNACO, doNotFit=['AT3', 'AT4', 'ATPUP scale'], plot=plot, N=N) print 'done in', round(time.time()-t0,2), \ 's (', round((time.time()-t0)/float(N), 2),\ 's per iteration)' if plot: t = np.linspace(0, 2*np.pi, 100) pyplot.plot(naco['Delta RA cos(DEC)'], naco['Delta DEC'], '*r', label='NACO position') pyplot.plot(naco['Delta RA cos(DEC)']+np.cos(t)*13e-3, naco['Delta DEC']+np.sin(t)*13e-3, 'r', label='13mas circle') pyplot.legend(numpoints=1) return b