##^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^Lance Simms, Stanford University 2009
#! /usr/bin/env python

##This script will run through a given fits file trying various values of 
##C_threshold using DAOFIND

import sys,os
from pyraf import iraf 
import pyfits
from numarray import *
from pylab import *

#Import the good packages from IRAF
iraf.digiphot(_doprint=0)
iraf.daophot(_doprint=0)

#DATAPARS
datapars=iraf.datapars.getParList()
findpars=iraf.findpars.getParList()
centerpars=iraf.centerpars.getParList()

print datapars[0]
iraf.centerpars.setParam('cbox','10.0')
iraf.datapars.setParam('fwhmpsf','4.0')
iraf.centerpars.saveParList(filename='~/iraf/uparm/datcentes.par')
iraf.datapars.saveParList(filename='~/iraf/uparm/datdataps.par')

#Set up daofind to go without prompting for input
iraf.daofind.setParam('image','hr5132a_001_raw.fits')	#Set filename
iraf.daofind.setParam('output','')
iraf.daofind.setParam('verify','no')			#Don't verify
iraf.daofind.saveParList(filename='daofind.par')	#Save values

daostring=iraf.daofind(mode='h',Stdout=1)		#Run hidden
lendaostring=len(daostring)
print lendaostring
first_string=daostring[0]
second_string=daostring[1]
third_string=daostring[2]

splitdao1=daostring[1].split()		#[0] is empty [1] contains header
splitdao2=daostring[3].split()		#[2] is empty [3]-end contain data

print first_string
print second_string
print third_string
print splitdao1
print splitdao2