mailr19381 - /branches/relax_disp/specific_analyses/relax_disp.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on April 05, 2013 - 09:57:
Author: bugman
Date: Fri Apr  5 09:57:42 2013
New Revision: 19381

URL: http://svn.gna.org/viewcvs/relax?rev=19381&view=rev
Log:
Created a custom base_data_loop() method for the relaxation dispersion 
analysis.

This defines the base data as the peak intensities of a single exponential 
curve and yields the spin
container and exponential curve key identifying the individual curves.


Modified:
    branches/relax_disp/specific_analyses/relax_disp.py

Modified: branches/relax_disp/specific_analyses/relax_disp.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/specific_analyses/relax_disp.py?rev=19381&r1=19380&r2=19381&view=diff
==============================================================================
--- branches/relax_disp/specific_analyses/relax_disp.py (original)
+++ branches/relax_disp/specific_analyses/relax_disp.py Fri Apr  5 09:57:42 
2013
@@ -53,7 +53,6 @@
         super(Relax_disp, self).__init__()
 
         # Place methods into the API.
-        self.base_data_loop = self._base_data_loop_spin
         self.data_init = self._data_init_spin
         self.model_loop = self._model_loop_spin
         self.return_conversion_factor = self._return_no_conversion_factor
@@ -970,6 +969,31 @@
 
         # Printout.
         print("Setting the '%s' spectrum spin-lock field strength to %s 
kHz." % (spectrum_id, cdp.spin_lock_nu1[spectrum_id]/1000.0))
+
+
+    def base_data_loop(self):
+        """Custom generator method for looping over spins and exponential 
curves.
+
+        The base data is hereby identified as the peak intensity data 
defining a single exponential curve.
+
+
+        @return:    The tuple of the spin container and the exponential 
curve identifying key (the CPMG frequency or R1rho spin-lock field strength).
+        @rtype:     tuple of SpinContainer instance and float
+        """
+
+        # Loop over the sequence.
+        for spin in spin_loop():
+            # Skip deselected spins.
+            if not spin.select:
+                continue
+
+            # Skip spins with no peak intensity data.
+            if not hasattr(spin, 'intensities'):
+                continue
+
+            # Loop over each exponential curve.
+            for exp_index, key in self._exp_curve_loop():
+                yield spin, key
 
 
     def create_mc_data(self, data_id):




Related Messages


Powered by MHonArc, Updated Fri Apr 05 12:20:02 2013