Changeset 0ac3db5 in sasmodels for sasmodels/bumps_model.py

Timestamp:

Feb 18, 2015 3:01:46 PM (9 years ago)

Author:

jhbakker

Branches:

master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

a006b63

Parents:

e806077

Message:

SESANSfit now uses sesans_reader dataloader

File:

• sasmodels/bumps_model.py (modified) (5 diffs)

sasmodels/bumps_model.py

@@ -217 +217 @@
 def _plot_sesans(data, theory, view):
     import matplotlib.pyplot as plt
-    resid = (theory - data.data)/data.err_data
+    resid = (theory - data.y)/data.dy
     plt.subplot(121)
-    plt.errorbar(data.SElength, data.data, yerr=data.err_data)
-    plt.plot(data.SElength, theory, '-', hold=True)
+    plt.errorbar(data.x, data.y, yerr=data.dy)
+    plt.plot(data.x, theory, '-', hold=True)
     plt.xlabel('spin echo length (A)')
     plt.ylabel('polarization')
     plt.subplot(122)
-    plt.plot(data.SElength, resid, 'x')
+    plt.plot(data.x, resid, 'x')
     plt.xlabel('spin echo length (A)')
     plt.ylabel('residuals')
@@ -267 +267 @@
         self.model = model
         self.cutoff = cutoff
-        if hasattr(data, 'SElength'):
+# TODO       if  isinstance(data,SESANSData1D)        
+        if hasattr(data, 'lam'):
             self.data_type = 'sesans'
         elif hasattr(data, 'qx_data'):
@@ -278 +279 @@
         # interpret data
         if self.data_type == 'sesans':
-            q = sesans.make_q(data.q_zmax, data.Rmax)
+            q = sesans.make_q(data.sample.zacceptance, data.Rmax)
             self.index = slice(None,None)
-            self.iq = data.data
-            self.diq = data.err_data
+            self.iq = data.y
+            self.diq = data.dy
             self._theory = np.zeros_like(q)
             q_vectors = [q]
@@ -353 +354 @@
             #self._theory[:] = self._fn.eval(pars, pd_pars)
             if self.data_type == 'sesans':
-                P = sesans.hankel(self.data.SElength, self.data.wavelength,
-                                  self.data.thickness, self._fn_inputs[0],
+                P = sesans.hankel(self.data.x, self.data.lam,
+                                  self.data.sample.thickness/10, self._fn_inputs[0],
                                   self._theory)
                 self._cache['theory'] = P
@@ -401 +402 @@
             self.data.data[self.index] = y
         elif self.data_type == 'sesans':
-            self.data.data[self.index] = y
+            self.data.y[self.index] = y
         else:
             raise ValueError("Unknown model")