def make_bias(path=None):
"""
Make master bias from raw data files
"""
if not path:
path = set_path()
biases,bct = tp.get_files(dir=path, tag='Bias')
masterbias = tp.master_bias(biases,outdir='./')
return masterbias
pix0 = w.wcs_world2pix(world0,1)
x0 = pix0[0,0]
y0 = pix0[0,1]
plt.figure(1)
plt.plot([x0],[y0],'y+',markersize=30)
plt.xlim(0,2048)
plt.ylim(0,2048)
darks,dct = tp.get_files(tag='Dark')
biases,bct = tp.get_files(tag='Bias')
flats,fct = tp.get_files(tag='SkyFlat')
targetfiles,tct = tp.get_files(tag='KIC10935310')
bias = tp.master_bias(biases,outdir='./Photometry/',readnoise=False)
dark = tp.master_dark(darks,bias=bias,outdir='./Photometry/')
flat = tp.master_flat(flats,bias=bias,dark=dark,outdir='./Photometry/')
cal = (image - dark - bias)/flat
qi.display_image(cal,siglo=2,sighi=2)
plt.plot([x0],[y0],'y+',markersize=30,linewidth=1e6)
plt.xlim(0,2048)
plt.ylim(0,2048)
# define the aperture
radius,ynew,xnew,fwhm,aspect,snrmax,totflux,totap,chisq = \
tp.optimal_aperture(x0,y0,image,skyrad=[15,20])
import thacherphot as tp
import matplotlib as mpl
from matplotlib.dates import DayLocator, HourLocator, DateFormatter, drange
import matplotlib.gridspec as gridspec
import datetime
import jdutil
from astropysics.coords import AngularCoordinate as angcor
path = '/Users/jonswift/Astronomy/ThacherObservatory/data/22Dec2016/'
#path = '/Users/ONeill/astronomy/Data/04Dec2016/'
files,fct = tp.get_files(dir=path, prefix='WASP-50b',suffix='solved.fits')
bfiles,bct = tp.get_files(dir=path,prefix='WASP-50b',suffix='bias.fit')
master_bias = tp.master_bias(bfiles)
data = tp.batch_total_flux(files,object='WASP-50b',mag=11.6,SpT='g8v',bias=master_bias,
nearest=True,camera='iKON-L',filter='V',gain=3.0,network='swift')
#tau = data['tau']
#secz = data['secz']
#plt.clf()
#plt.ion()
#plt.figure()
#plt.plot(secz,tau,'o')
#plt.ylabel('tau')
#plt.xlabel('secz')
#plt.xlim(1.85,1.95)
#plt.ylim(0.64,0.67)
