# -*- coding: utf-8 -*-
"""
Created on Wed Apr 29 19:07:38 2015

@author: slacour
"""

from numpy import *   
from matplotlib.pyplot import *    
import pyfits as fits

dir = "../DATA/data_2015-05-26_Julien/"
sof = "p2vm_data.sof"

text=open(dir+sof,'r').readlines()
files={}
flat=0
for line in text:
    if line.split()[1] == 'DARK_RAW':
        files[line.split()[1]]=dir+line.split()[0]

dark=files['DARK_RAW']

dark_sc = fits.getdata(dark,"IMAGING_DATA_SC")
dark_ft = fits.getdata(dark,"IMAGING_DATA_FT").field("PIX")


# FOR DARK CALCULATION 

median_dark_sc=median(dark_sc,axis=0)
median_dark_ft=median(dark_ft,axis=0)
var_dark_sc=mean(((dark_sc-median_dark_sc)**2),axis=0)
var_dark_ft=mean(((dark_ft-median_dark_ft)**2),axis=0)

# FOR DARK QC PARAMETERS

QC_med_sc=median(median_dark_sc)
QC_rms_sc=sqrt(median(var_dark_sc))
print("QC Median SC = "+str(QC_med_sc))
print("QC Median RMS SC = "+str(QC_rms_sc))


QC_mea_ft=mean(median_dark_ft)
QC_rms_ft=sqrt(mean(var_dark_ft))
print("QC Mean FT = "+str(QC_mea_ft))
print("QC Mean RMS FT = "+str(QC_rms_ft))