--- /dev/null
+The output of analyze.py was mislabeled as pN/nm, when it should be N/m.
+
+Variable (units) : mean +/- std. dev. (relative error)
+Cantilever k (pN/nm) : 2.98914 +/- 0.333093 (0.111434)
+photoSensitivity**2 (V/nm)**2 : 0.000859293 +/- 2.5841e-05 (0.0300723)
+T (K) : 295.15 +/- 2.27374e-13 (7.70367e-16)
+1/Vphoto**2 (1/V)**2 : 853648 +/- 91596.3 (0.1073)
+
--- /dev/null
+
+
+
+Content-type=text/plain
+
+
+
+
+
+
+Date=Mon, 17 Nov 2008 14:59:48 +0000
+
+
+
+
+
+
+From=wking
+
+
+
--- /dev/null
+
+
+
+creator=wking
+
+
+
+
+
+
+severity=serious
+
+
+
+
+
+
+status=fixed
+
+
+
+
+
+
+summary=analyze.py unit error for spring constant.
+
+
+
+
+
+
+time=Mon, 17 Nov 2008 14:56:57 +0000
+
+
+
--- /dev/null
+
+
+
+rcs_name=git
+
+
+
--- /dev/null
+Bugs Everywhere Tree 1 0
Inputs (all are arrays of recorded data) :
bumps measured (V_photodiode / nm_tip) proportionality constant
Ts measured temperature (K)
- vibs measured V_photodiode variance in free solution
+ vibs measured V_photodiode variance in free solution (V**2)
Outputs :
k cantilever spring constant (in N/m, or equivalently nN/nm)
k_s standard deviation in our estimate of k
T_r = T_s / T_m
one_o_Vphoto2_r = one_o_Vphoto2_s / one_o_Vphoto2_m
string = "Variable (units) : mean +/- std. dev. (relative error)\n"
- string += "Cantilever k (pN/nm) : %g +/- %g (%g)\n" \
+ string += "Cantilever k (N/m) : %g +/- %g (%g)\n" \
% (k_m, k_s, k_r)
string += "photoSensitivity**2 (V/nm)**2 : %g +/- %g (%g)\n" \
% (photoSensitivity2_m, photoSensitivity2_s, photoSensitivity2_r)
vib_plot(deflection_bits, freq_axis, power, A, B, C, plotVerbose=plotVerbose)
- # Our A is in uV**2, so convert back to Volts
+ # Our A is in uV**2, so convert back to Volts**2
return lorentzian_area(A,B,C) * 1e-12
def fit_psd(freq_axis, psd_data, minFreq=500, maxFreq=7000) :