Changeset e806077 in sasmodels for example

Timestamp:

Feb 18, 2015 9:25:56 AM (10 years ago)

Author:

pkienzle

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:

3a45c2c, 5428233, 0ac3db5

Parents:

c97724e

Message:

add example of sans/sesans simultaneous system

File:

• example/sesansfit.py (modified) (2 diffs)

example/sesansfit.py

@@ -5 +5 @@
 kernel = sas.load_model('sphere', dtype='single')
 
+
 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')
+else:
+    SElength = np.linspace(0, 2400, 61) # [A]
+    data = np.ones_like(SElength)
+    err_data = np.ones_like(SElength)*0.03
 
-SElength = np.linspace(0, 2400, 61) # [A]
-data = np.ones_like(SElength)
-err_data = np.ones_like(SElength)*0.03
-
-class SesansData:
-    #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]
-    thickness = 0.2 # [cm]
-    SElength = SElength
-    data = data
-    err_data = err_data
-data = SesansData()
+    class SesansData:
+        #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]
+        thickness = 0.2 # [cm]
+        SElength = SElength
+        data = data
+        err_data = err_data
+    data = SesansData()
 
 phi = Parameter(0.1, name="phi")
@@ -31 +36 @@
 model.background.range(0,5)
 
-problem = FitProblem(model)
 
+if False: # have sans data
+    sansmodel = sas.BumpsModel(sans_data, kernel, **model.parameters())
+    problem = FitProblem([model, sansmodel])
+else:
+    problem = FitProblem(model)
+