Changeset 2684d45 in sasmodels

Timestamp:

Mar 21, 2016 5:05:46 AM (9 years ago)

Author:

awashington

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:

a6ebff8

Parents:

3f167d9

Message:

Force proper units in SESANS analysis

Files:

• example/sesans_parameters_css-hs.py (modified) (1 diff)
• sasmodels/sesans.py (modified) (4 diffs)

example/sesans_parameters_css-hs.py

@@ -22 +22 @@
 # Initial parameter values (if other than defaults)
 initial_vals = {
-    "sld_core" : 1.0592,
-    "sld_solvent" : 2.88,
-    "radius" : 890,
-    "thickness" : 130
+    "sld_core" : 1.83,
+    "sld_shell" : 1.83,
+    "sld_solvent" : 0,
+    "radius" : 730,
+    "thickness" : 20,
+    "volfraction" : phi,
+    "scale" : (1-phi)
 }
 

sasmodels/sesans.py

@@ -29 +29 @@
 
     q_min = dq = 0.1 * 2*pi / Rmax
-    return np.arange(q_min, Converter("1/A")(q_max[0], units=q_max[1]), dq)
+    return np.arange(q_min,
+                     Converter(q_max[1])(q_max[0],
+                                         units="1/A"),
+                     dq)
     
 def make_all_q(data):
@@ -75 +78 @@
 
 def call_hankel(data, q_calc, Iq_calc):
-    return hankel(data.x, data.lam * 1e-9,
-                  data.sample.thickness / 10,
+    return hankel((data.x, data.x_unit),
+                  (data.lam, data.lam_unit),
+                  (data.sample.thickness,
+                   data.sample.thickness_unit),
                   q_calc, Iq_calc)
   
@@ -182 +187 @@
     *I* [cm$^{-1}$] is the value of the SANS model at *q*
     """
+
+    from sas.sascalc.data_util.nxsunit import Converter
+    wavelength = Converter(wavelength[1])(wavelength[0],"A")
+    thickness = Converter(thickness[1])(thickness[0],"A")
+    Iq = Converter("1/cm")(Iq,"1/A") # All models default to inverse centimeters
+    SElength = Converter(SElength[1])(SElength[0],"A")
+
     G = np.zeros_like(SElength, 'd')
 #==============================================================================
@@ -192 +204 @@
 
     # [m^-1] step size in q, needed for integration
-    dq = (q[1]-q[0])*1e10
+    dq = (q[1]-q[0])
 
     # integration step, convert q into [m**-1] and 2 pi circle integration
-    G *= dq*1e10*2*pi
-    G0 *= dq*1e10*2*pi
+    G *= dq*2*pi
+    G0 = np.sum(Iq*q)*dq*2*np.pi
 
     P = exp(thickness*wavelength**2/(4*pi**2)*(G-G0))