posts:comparing_velocity_clamp_experiments: add expected unfolding force line.

For comparing crunched results against previously published numbers.

diff --git a/posts/Comparing_velocity_clamp_experiments/crunch.py b/posts/Comparing_velocity_clamp_experiments/crunch.py
index 99fc50ab7ff67adadb3c534f5067171e58271e67..5b261c8c598cffa26eabfbf4e51b31e31de2a48f 100755 (executable)
--- a/posts/Comparing_velocity_clamp_experiments/crunch.py
+++ b/posts/Comparing_velocity_clamp_experiments/crunch.py
@@ -98,6 +98,8 @@ pull_parameters = {
     #'contour_space_parameters': None,  # use defaults
     'contour_space_parameters': {'b': _b, 'gamma': _gamma},
     'contour_space_lines': _numpy.arange(75e-9, 350e-9, 29e-9),
+    # Carrion-Vazquez, 1999, doi: 10.1073/pnas.96.7.3694
+    'expected_unfolding_force': 200e-12,
     }
 
 ## readers for your data
@@ -594,7 +596,7 @@ class Pull (Experiment):
             'vibration_temperature', 'drift_slope', 'alpha_f_offset',
             'smooth_steps', 'start_sawteeth', 'stop_sawteeth',
             'contour_space_model', 'contour_space_parameters',
-            'contour_space_lines']
+            'contour_space_lines', 'expected_unfolding_force']
         kwargs['output_parameters'] = [
             'pull_alpha_d', 'pull_velocity', 'pull_temperature']
         kwargs['output_parameter_units'] = {
@@ -653,7 +655,8 @@ class Pull (Experiment):
               xlabel='z piezo (m)', ylabel='deflection (N)',
               bounds=None, log=None, slider=None,
               contour_space_model=None, contour_space_parameters=None,
-              contour_space_lines=None, **kwargs):
+              contour_space_lines=None, expected_unfolding_force=None,
+              **kwargs):
         """Generate pull-plots with three possible abscissa.
 
         The contour-space plots are following Puchner et al., 2008.
@@ -709,6 +712,8 @@ class Pull (Experiment):
         if xlabel == 'contour-space (m)' and contour_space_lines is not None:
             for zi in contour_space_lines:
                 axes.axvline(zi, color='g')
+        if ylabel == 'deflection (N)' and expected_unfolding_force is not None:
+            axes.axhline(expected_unfolding_force, color='g')
         if bounds:
             axes.axis(bounds)
         axes.set_title('velocity clamp unfolding')

diff --git a/posts/Comparing_velocity_clamp_experiments/pull-piezo.png b/posts/Comparing_velocity_clamp_experiments/pull-piezo.png
index aa0af7896e2c5cd8135f25755c20db2633fefe33..eccd7fb92f06614ba868b5d65f1e8fd2c473bfa9 100644 (file)
Binary files a/posts/Comparing_velocity_clamp_experiments/pull-piezo.png and b/posts/Comparing_velocity_clamp_experiments/pull-piezo.png differ

diff --git a/posts/Comparing_velocity_clamp_experiments/pull-protein.png b/posts/Comparing_velocity_clamp_experiments/pull-protein.png
index af0cf52a0854226099f8bec6f606de0b64000c1b..e2b8cccee829bfc7ff333ab8e693e47b5cb7e3c7 100644 (file)
Binary files a/posts/Comparing_velocity_clamp_experiments/pull-protein.png and b/posts/Comparing_velocity_clamp_experiments/pull-protein.png differ