import os
import csv
from make_assoc import *
from astropy.io import fits






f_psf_data = open('psf_fwhm_data_list_A.txt', 'w')
with open('list_fits_files.csv','r') as csvfile:
    rows = csv.reader(csvfile, delimiter=';')
    i=0
    for row in rows:
      print row
      os.chdir(row[2])
      hdulist = fits.open('res_catalog_A.psf')
      hdulist2 = fits.open('res_gal.fits')
      hdr = hdulist[1].header
      dat = hdulist2[2].data
      if len(dat)>0:   
        fwhm_gal = dat[0]["FWHM_IMAGE"]
        fwhm = hdr["PSF_FWHM"]
        ratio=fwhm_gal/fwhm
        f_psf_data.write(str(row[2])+';'+str(fwhm)+';'+str(fwhm_gal)+';'+str(ratio)+'\n')
      else:
        print 'ALARM', row[2]
      os.chdir('../../../')
    f_psf_data.close()