EROSITAwiki: attachment:call_get_param.py of VariabilitySearch

   1 import get_param
   2 import numpy as np
   3 lcfil="flare_am15.0_du03.9_040_LightCurve_source_1.fits"
   4 
   5 chi2_lin,chi2_flare,chi2_eclipse,snr_flare,snr_eclipse=get_param.call_char_var(lcfil)
   6 #params for linear fit
   7 lich2ts = float(chi2_lin[0]) #chi2
   8 lich2pv = float(chi2_lin[1]) #p-value
   9 lich2nd = int(chi2_lin[2])   #d.o.f
  10 lich2p1 = float(chi2_lin[3]) #param1
  11 lich2e1 = float(chi2_lin[4]) #err_param1
  12 lich2p2 = float(chi2_lin[5]) #param2
  13 lich2e2 = float(chi2_lin[6]) #err_param2
  14 #params for flare fit
  15 flch2ts = float(chi2_flare[0]) #chi2
  16 flch2pv = float(chi2_flare[1]) #p-value
  17 flch2nd = int(chi2_flare[2])   #d.o.f 
  18 flch2p1 = float(chi2_flare[3]) #param1
  19 flch2e1 = float(chi2_flare[4]) #err_param1
  20 flch2p2 = float(chi2_flare[5]) #param2
  21 flch2e2 = float(chi2_flare[6]) #err_param2
  22 flch2p3 = float(chi2_flare[7]) #param3
  23 flch2e3 = float(chi2_flare[8]) #err_param3
  24 flch2p4 = float(chi2_flare[9]) #param4
  25 flch2e4 = float(chi2_flare[10])#err_param4 
  26 #params for eclipse fit
  27 ecch2ts = float(chi2_eclipse[0]) #chi2
  28 ecch2pv = float(chi2_eclipse[1]) #p-value
  29 ecch2nd = int(chi2_eclipse[2])   #d.o.f  
  30 ecch2p1 = float(chi2_eclipse[3]) #param1
  31 ecch2e1 = float(chi2_eclipse[4]) #err_param1
  32 ecch2p2 = float(chi2_eclipse[5]) #param2
  33 ecch2e2 = float(chi2_eclipse[6]) #err_param2
  34 ecch2p3 = float(chi2_eclipse[7]) #param3
  35 ecch2e3 = float(chi2_eclipse[8]) #err_param3
  36 ecch2p4 = float(chi2_eclipse[9]) #param4
  37 ecch2e4 = float(chi2_eclipse[10])#err_param4 
  38 
  39 print "p-val flare over linear (log):",np.log10(flch2pv/lich2pv)
  40 print "p-val eclipse over linear (log):",np.log10(ecch2pv/lich2pv)
  41 
  42 print "SNR flare:",snr_flare
  43 print "SNR eclipse:",snr_eclipse