From: illysam <devnull@localhost>
Date: Sat, 6 Feb 2010 13:36:44 +0000 (+0000)
Subject: added plotmanipulator centerzero to procplots.py
X-Git-Url: http://git.tremily.us/?a=commitdiff_plain;h=65b49b7589e5c748d0921ec37a5cddddfa8b8ee7;p=hooke.git

added plotmanipulator centerzero to procplots.py
---

diff --git a/plugins/core.ini b/plugins/core.ini
index ac026ec..bcfeb5c 100644
--- a/plugins/core.ini
+++ b/plugins/core.ini
@@ -1,16 +1,21 @@
 [plotmanipulators]
     [[order]]
-        default = '"correct", "median", "absvalue", "flatten", "multiplier", "clamp", "threshold", "coincident", "showconvoluted"'
-        elements = correct, median, absvalue, flatten, multiplier, clamp, threshold, coincident, showconvoluted
-        #elements = correct, median, absvalue, flatten, multiplier, clamp, threshold, coincident, showconvoluted
+        default = '"correct", "median", "absvalue", "flatten", "multiplier", "clamp", "threshold", "coincident", "showconvoluted", "centerzero"'
+        elements = correct, median, absvalue, flatten, multiplier, clamp, threshold, coincident, showconvoluted, centerzero
+        #elements = correct, median, absvalue, flatten, multiplier, clamp, threshold, coincident, showconvoluted, centerzero
         type = arraystring
-        value = '"correct" "median" "absvalue" "flatten" "multiplier" "clamp" "threshold" "coincident" "showconvoluted"'
+        value = '"correct" "median" "absvalue" "flatten" "multiplier" "clamp" "threshold" "coincident" "showconvoluted" "centerzero"'
     
     [[absvalue]]
         default = False
         type = boolean
         value = False
     
+    [[centerzero]]
+        default = False
+        type = boolean
+        value = True
+    
     [[clamp]]
         default = False
         type = boolean
diff --git a/plugins/procplots.py b/plugins/procplots.py
index 7324b77..b7dd12c 100644
--- a/plugins/procplots.py
+++ b/plugins/procplots.py
@@ -16,7 +16,7 @@ import wxversion
 wxversion.select(lh.WX_GOOD)
 
 import copy
-from numpy import arange, diff, fft
+from numpy import arange, diff, fft, median
 from scipy.signal import medfilt
 
 from lib.peakspot import conv_dx
@@ -80,6 +80,7 @@ class procplotsCommands:
         '''
         column = self.GetIntFromConfig('procplots', 'derivplot', 'column')
         row = self.GetIntFromConfig('procplots', 'derivplot', 'row')
+        #TODO: what os select good for?
         select = self.GetBoolFromConfig('procplots', 'derivplot', 'select')
         whatset_str = self.GetStringFromConfig('procplots', 'derivplot', 'whatset')
         whatset = []
@@ -190,6 +191,31 @@ class procplotsCommands:
 
         return plot
 
+    def plotmanip_centerzero(self, plot, current, customvalue=False):
+        '''
+        Centers the force curve so the median (the free level) corresponds to 0 N
+        '''
+        #use only for force spectroscopy experiments!
+        if current.driver.experiment != 'smfs':
+            return plot
+    
+        if not customvalue:
+            customvalue = self.GetBoolFromConfig('procplots', 'centerzero')
+        if not customvalue:
+            return plot
+
+        levelapp = float(median(plot.curves[lh.EXTENSION].y))
+        levelret = float(median(plot.curves[lh.RETRACTION].y))
+    
+        level = (levelapp + levelret)/2    
+    
+        plot.curves[lh.EXTENSION].y = [i-level for i in plot.curves[lh.EXTENSION].y]
+        plot.curves[lh.RETRACTION].y = [i-level for i in plot.curves[lh.RETRACTION].y]
+        
+#        plot.vectors[0][1]=approach    
+#        plot.vectors[1][1]=retract    
+        return plot
+    
 #FFT---------------------------
     def fft_plot(self, curve, boundaries=[0, -1]):
         '''