import os, sys, glob
reddir = '[PATHTOZSALT]/zSALT/zsalt/'
sys.path.append(reddir)
bpmfile = '[PATHTOBPM]/bpm.fits'
propid = '[PROPID]'
from imred import imred
from specred import specred
ddir = sys.argv[1]
os.chdir(ddir)
if not os.path.isdir('sci'): os.mkdir('sci')
os.chdir('sci')
#basic image reuctions
infile_list = glob.glob('../raw/P*fits')
imred(infile_list, './', bpmfile, cleanup=True)
#spectroscopic reductions
infile_list = glob.glob('m*fits')
specred(infile_list, propid, inter=True)
#UNEDIT IF YOU WANT TO REDUCE MOS DATA
#mosxml = [MOSXMLFILE]
#mosred = infile_list, propid, mosxml, inter=True)
from specpolextract_sc import specpolextract_sc
from specpolrawstokes import specpolrawstokes
from specpolfinalstokes import specpolfinalstokes
obsdate = sys.argv[1]
os.chdir(obsdate)
if not os.path.isdir('sci'): os.mkdir('sci')
shutil.copy(scrdir + 'script.py', 'sci')
os.chdir('sci')
#basic image reductions
infilelist = sorted(glob.glob('../raw/P*fits'))
imred(infilelist,
'./',
datadir + 'bpm_rss_11.fits',
crthresh=False,
cleanup=True)
#basic polarimetric reductions
logfile = 'specpol' + obsdate + '.log'
#wavelength map
infilelist = sorted(glob.glob('m*fits'))
linelistlib = ""
specpolwavmap(infilelist, linelistlib=linelistlib, logfile=logfile)
#background subtraction and extraction
infilelist = sorted(glob.glob('wm*fits'))
dyasec = 5. # star +/-5, bkg= +/-(25-35)arcsec: 2nd order is 9-20 arcsec away
yoffo = 0. # optional offset of target (bins) from brightest in O (bottom) image
os.chdir('sci')
#basic image reuctions
infile_list = glob.glob('../raw/P*fits')
for img in infile_list:
hdu = fits.open(img, 'update')
hdu[1].header['XTALK'] = 1474
hdu[2].header['XTALK'] = 1474
hdu[3].header['XTALK'] = 1166
hdu[4].header['XTALK'] = 1111
hdu[5].header['XTALK'] = 1377
hdu[6].header['XTALK'] = 1377
hdu.close()
imred(infile_list, './', bpmfile, cleanup=True)
#spectroscopic reductions
infile_list = glob.glob('m*fits')
specred(infile_list, propid, inter=True)
#UNEDIT IF YOU WANT TO REDUCE MOS DATA
#mosxml = [MOSXMLFILE]
#mosred = infile_list, propid, mosxml, inter=True)
from specpolwavmap import specpolwavmap
from specpolextract import specpolextract
from specpolrawstokes import specpolrawstokes
from specpolfinalstokes import specpolfinalstokes
obsdate = sys.argv[1]
os.chdir(obsdate)
if not os.path.isdir('sci'): os.mkdir('sci')
os.chdir('sci')
#basic image reductions
infile_list = glob.glob('../raw/P*fits')
#imred(infile_list, './', datadir+'bpm_rss_11.fits', cleanup=False)
imred(infile_list, './', datadir+'bpm_rss_11.fits', cleanup=True)
#basic polarimetric reductions
# debug=True
debug=False
logfile='specpol'+obsdate+'.log'
#target and wavelength map
infile_list = sorted(glob.glob('m*fits'))
linelistlib=""
specpolwavmap(infile_list, linelistlib=linelistlib, logfile=logfile)
#background subtraction and extraction
infile_list = sorted(glob.glob('wm*fits'))
specpolextract(infile_list, logfile=logfile, debug=debug)
args.slitmask = os.path.abspath(args.slitmask)
ddir = args.ddir
os.chdir(ddir)
if not os.path.isdir('sci'): os.mkdir('sci')
os.chdir('sci')
#basic image reuctions
infile_list = glob.glob('../raw/P*fits')
# remove alignment images
spec_list = []
for infile in infile_list:
hdu = fits.open(infile)
if hdu[0].header['OBSMODE'].strip() != 'IMAGING':
spec_list.append(infile)
infile_list = spec_list
print spec_list
if args.basic_red: imred.imred(infile_list, './', bpmfile, cleanup=True)
#spectroscopic reductions
propid=None
infile_list = glob.glob('m*fits')
mosred(infile_list, args.slitmask, dy=args.dy, preprocess=args.preprocess)
#extract the data
from specpolwavmap import specpolwavmap
from specpolextract import specpolextract
from specpolrawstokes import specpolrawstokes
from specpolfinalstokes import specpolfinalstokes
obsdate = sys.argv[1]
os.chdir(obsdate)
if not os.path.isdir('sci'): os.mkdir('sci')
os.chdir('sci')
#basic image reductions
infile_list = glob.glob('../raw/P*fits')
#imred(infile_list, './', datadir+'bpm_rss_11.fits', cleanup=False)
imred(infile_list, './', datadir+'bpm_rss_11.fits', cleanup=True)
#basic polarimetric reductions
# debug=True
debug=False
logfile='specpol'+obsdate+'.log'
#target and wavelength map
infile_list = sorted(glob.glob('m*fits'))
linelistlib=""
specpolwavmap(infile_list, linelistlib=linelistlib, inter=True, logfile=logfile)
#background subtraction and extraction
#infile_list = sorted(glob.glob('wm*099.fits')+glob.glob('wm*10[0-6].fits')) # optional subselection
from specpolwavmap import specpolwavmap
from specpolextract_sc import specpolextract_sc
from specpolrawstokes import specpolrawstokes
from specpolfinalstokes import specpolfinalstokes
print sys.argv
obsdate = sys.argv[1]
print obsdate
os.chdir(obsdate)
if not os.path.isdir('sci'): os.mkdir('sci')
shutil.copy(scrdir+'script.py','sci')
os.chdir('sci')
#basic image reductions
infilelist = sorted(glob.glob('../raw/P*fits'))
imred(infilelist, './', datadir+'bpm_rss_11.fits', crthresh=False, cleanup=True)
#basic polarimetric reductions
logfile='specpol'+obsdate+'.log'
#wavelength map
infilelist = sorted(glob.glob('m*fits'))
linelistlib=""
specpolwavmap(infilelist, linelistlib=linelistlib, logfile=logfile)
#background subtraction and extraction
infilelist = sorted(glob.glob('wm*fits'))
extract = 10. # star +/-5, bkg= +/-(25-35)arcsec: 2nd order is 9-20 arcsec away
locate = (-120.,120.) # science target is brightest target in whole slit
#locate = (-20.,20.)
