import numpy
from splittable_kwargs import splittableKwargsFunction, \
make_splittable_kwargs_function
+import data_logger
+
import common # common module for the calibcant package
import config # config module for the calibcant package
import T_analyze # T_analyze module for the calibcant package
if log_dir != None :
log = data_logger.data_log(log_dir, noclobber_logsubdir=False,
log_name="calib_analysis_text")
- log.write_binary((
- "k (N/m) : %g +/- %g\n" % (k, k_s)
- + "photoSensitivity**2 (V/nm)**2 : %g +/- %g\n" % \
- (photoSensitivity2_m, photoSensitivity2_s)
- + "T (K) : %g +/- %g\n" % (T_m, T_s)
- + "1/Vp**2 (1/V)**2 : %g +/- %g\n" % \
- (one_o_Vphoto2_m, one_o_Vphoto2_s)
- ))
+ log.write_binary(string_errors(k, k_s,
+ photoSensitivity2_m, photoSensitivity2_s,
+ T_m, T_s,
+ one_o_Vphoto2_m, one_o_Vphoto2_s)
+ +'\n')
@splittableKwargsFunction()
def calib_plot(bumps, Ts, vibs, plotVerbose=False) :
if plotVerbose or config.PYLAB_VERBOSE :
- common.import_pylab()
+ common._import_pylab()
common._pylab.figure(config.BASE_FIGNUM+4)
common._pylab.subplot(311)
- common._pylab.plot(bumps, 'g.')
+ common._pylab.plot(bumps, 'g.-')
common._pylab.title('Photodiode sensitivity (V/nm)')
common._pylab.subplot(312)
- common._pylab.plot(Ts, 'r.')
+ common._pylab.plot(Ts, 'r.-')
common._pylab.title('Temperature (K)')
common._pylab.subplot(313)
- common._pylab.plot(vibs, 'b.')
+ common._pylab.plot(vibs, 'b.-')
common._pylab.title('Thermal deflection variance (Volts**2)')
common._flush_plot()
options,args = parser.parse_args()
parser.destroy()
-
+
+ config.TEXT_VERBOSE = options.verbose
+ config.PYLAB_INTERACTIVE = False
+ config.PYLAB_VERBOSE = options.plot
+
bumps = get_array(options.bump_string, options.bump_file, "bump")
temps = get_array(options.temp_string, options.temp_file, "temp")
vibs = get_array(options.vib_string, options.vib_file, "vib")
- if options.plot :
- calib_plot(bumps, temps, vibs)
+ #if options.plot :
+ # calib_plot(bumps, temps, vibs)
if options.celsius :
for i in range(len(temps)) :
temps[i] = T_analyze.C_to_K(temps[i])
- km,ks,ps2m,ps2s,Tm,Ts,ooVp2m,ooVp2s = calib_analyze(bumps, temps, vibs)
+ km,ks,ps2m,ps2s,Tm,Ts,ooVp2m,ooVp2s = \
+ calib_analyze(bumps, temps, vibs, plotVerbose=options.plot)
if options.verbose :
print >> sys.stderr, \