Changeset 0ac3db5 in sasmodels

Timestamp:

Feb 18, 2015 5:01:46 PM (10 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

Files:

• example/sesansfit.py (modified) (2 diffs)
• sasmodels/bumps_model.py (modified) (5 diffs)

example/sesansfit.py

@@ -6 +6 @@
 
 
-radius = 1000
-if False: # fix when data loader exists
-    from sasmodels.dataloader import load_sesans, load_sans
-    data = load_sesans('mydatfile.pz')
-    sans_data = load_sans('mysansfile.xml')
+if True: # fix when data loader exists
+#    from sas.dataloader.readers\
+    from sas.dataloader.loader import Loader
+    loader=Loader()
+    data=loader.load('testsasview1.ses')
+    
+#    data = load_sesans('mydatfile.pz')
+#    sans_data = load_sans('mysansfile.xml')
+
 else:
     SElength = np.linspace(0, 2400, 61) # [A]
@@ -16 +20 @@
     err_data = np.ones_like(SElength)*0.03
 
-    class SesansData:
+    class SESANSData1D:
         #q_zmax = 0.23 # [A^-1]
-        q_zmax = 0.1 # [A^-1]
-        SElength = np.linspace(0, 2400, 61) # [A]
-        Rmax = 3*radius # [A]
-        wavelength = 2e-10 # [m]
+        zacceptance = 0.1 # [A^-1]
+        lam = 2e-10 # [m]
         thickness = 0.2 # [cm]
-        SElength = SElength
-        data = data
-        err_data = err_data
+        x = SElength
+        y = data
+        dy = err_data
     data = SesansData()
+print dir(data)
+
+radius = 1000
+data.Rmax = 3*radius # [A]
 
 phi = Parameter(0.1, name="phi")

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")