#! /usr/bin/env python 
#
# This program extracts the FWHM measurements from the .mag files using IRAF
#

from pyfits import getheader
from os  import getcwd, chdir, remove, path, environ, getenv, makedirs, listdir, system, unlink
import os.path
from pyraf import iraf
from pylab import *
import pyraf.iraf as iraf

from numpy import *
import numpy as numpy

def load_phot():
   iraf.noao(_doprint=0)
   iraf.digiphot(_doprint=0)
   iraf.daophot(_doprint=0)
   iraf.tables(_doprint=0)

workdir = getcwd()
magdir=workdir+"/mag/"

iraf.noao(_doprint=0)
#iraf.obsutil(_doprint=0)
#iraf.unlearn("psfmeasure")

chdir(magdir)

load_phot()
iraf.txdump(textfiles="*mag.1", fields="XCEN,YCEN,XSHIFT,YSHIFT,XERR,YERR,MSKY,STDEV,NSKY,NSREJ", parameters="no", headers="no", expr="yes", Stdout="info-all.list")
os.system('mv info-all.list ../.')


'''		FWHM
iraf.noao.obsutil.psfmeasure.images = "@images.list"
iraf.noao.obsutil.psfmeasure.coords = "mark1"
iraf.noao.obsutil.psfmeasure.display = "no"
iraf.noao.obsutil.psfmeasure.size = "FWHM"
iraf.noao.obsutil.psfmeasure.radius = "15"
iraf.noao.obsutil.psfmeasure.imagecu = "ds92.reg"
iraf.noao.obsutil.psfmeasure.logfile = "fwhm.log"


if os.path.isfile("fwhm.log"):
  os.system( 'rm fwhm.log' )

  #dt = [('filename',str)]
iraf.noao.obsutil.psfmeasure("*.fits")

F = loadtxt('fwhm.log',skiprows = 3, usecols=(0,4), dtype = str)
fwhma = F[:,1]
fwhma = fwhma[0:len(F)-1]

fwhm = []
for i in range(len(fwhma)):
  fwhm.append(fwhma[i])

fwhm = array(fwhm)
  
#ave_fwhm = numpy.mean(fwhm, dtype=np.float32)
#print ave_fwhm

f = open('fwhm.list', 'w+')
for j in range(len(fwhm)):
  print >> f, fwhm[j], 1
f.close()

os.system( 'mv fwhm.list ../data' )
'''