import get_param
import numpy as np
lcfil="flare_am15.0_du03.9_040_LightCurve_source_1.fits"

chi2_lin,chi2_flare,chi2_eclipse,snr_flare,snr_eclipse=get_param.call_char_var(lcfil)
#params for linear fit
lich2ts = float(chi2_lin[0]) #chi2
lich2pv = float(chi2_lin[1]) #p-value
lich2nd = int(chi2_lin[2])   #d.o.f
lich2p1 = float(chi2_lin[3]) #param1
lich2e1 = float(chi2_lin[4]) #err_param1
lich2p2 = float(chi2_lin[5]) #param2
lich2e2 = float(chi2_lin[6]) #err_param2
#params for flare fit
flch2ts = float(chi2_flare[0]) #chi2
flch2pv = float(chi2_flare[1]) #p-value
flch2nd = int(chi2_flare[2])   #d.o.f 
flch2p1 = float(chi2_flare[3]) #param1
flch2e1 = float(chi2_flare[4]) #err_param1
flch2p2 = float(chi2_flare[5]) #param2
flch2e2 = float(chi2_flare[6]) #err_param2
flch2p3 = float(chi2_flare[7]) #param3
flch2e3 = float(chi2_flare[8]) #err_param3
flch2p4 = float(chi2_flare[9]) #param4
flch2e4 = float(chi2_flare[10])#err_param4 
#params for eclipse fit
ecch2ts = float(chi2_eclipse[0]) #chi2
ecch2pv = float(chi2_eclipse[1]) #p-value
ecch2nd = int(chi2_eclipse[2])   #d.o.f  
ecch2p1 = float(chi2_eclipse[3]) #param1
ecch2e1 = float(chi2_eclipse[4]) #err_param1
ecch2p2 = float(chi2_eclipse[5]) #param2
ecch2e2 = float(chi2_eclipse[6]) #err_param2
ecch2p3 = float(chi2_eclipse[7]) #param3
ecch2e3 = float(chi2_eclipse[8]) #err_param3
ecch2p4 = float(chi2_eclipse[9]) #param4
ecch2e4 = float(chi2_eclipse[10])#err_param4 

print "p-val flare over linear (log):",np.log10(flch2pv/lich2pv)
print "p-val eclipse over linear (log):",np.log10(ecch2pv/lich2pv)

print "SNR flare:",snr_flare
print "SNR eclipse:",snr_eclipse