--- /dev/null
+And here is my analysis of data on the same cantilever that I took:
+
+Variable (units) : mean +/- std. dev. (relative error)
+Cantilever k (N/m) : 0.0686691 +/- 0.00134223 (0.0195464)
+photoSensitivity**2 (V/nm)**2 : 0.000333038 +/- 1.88694e-06 (0.00566586)
+T (K) : 295.15 +/- 5.68434e-14 (1.92592e-16)
+1/Vphoto**2 (1/V)**2 : 50598.9 +/- 946.564 (0.0187072)
+
# read the data
data = T_load(path)
Ts.extend(data)
- T_analyze(Ts, convert_to_K=K_to_K)
+ T_analyze(Ts, convert_to_K=convert_to_K)
return numpy.array(Ts, dtype=numpy.float)
# commandline interface functions
)
parser = OptionParser(usage=usage_string, version='%prog 0.1')
parser.add_option('-C', '--celsius', dest='celsius',
- help='Use Celsius input temperatures instead of Kelvin (defaul %default)\n',
+ help='Use Celsius input temperatures instead of Kelvin (default %default)\n',
action='store_true', default=False)
parser.add_option('-o', '--output-file', dest='ofilename',
help='write output to FILE (default stdout)',
if options.celsius :
convert_to_K = C_to_K
else :
- convert_to_K = lambda T : T # no-change function
-
+ convert_to_K = K_to_K
+
if options.tweakmode == False :
data = read_data(ifilename)
- Ts = T_analyze(T, convert_to_K)
+ Ts = T_analyze(data, convert_to_K)
else : # tweak file mode
Ts = T_load_analyze_tweaked(ifilename, convert_to_K, textVerboseFile=vfile)