Changes in / [49da079:5430333] in sasmodels

Files:

• example/sesans_parameters_css-hs.py (modified) (1 diff)
• sasmodels/models/core_shell_sphere.py (modified) (1 diff)
• sasmodels/models/img/spherical_sld_1d.jpg (added)
• sasmodels/models/img/spherical_sld_default_profile.jpg (added)
• sasmodels/models/img/spherical_sld_profile.gif (added)
• 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_core" : 1.05,
+    "sld_shell" : 2.88*pen-0.05*(1-pen),
     "sld_solvent" : 2.88,
-    "radius" : 890,
-    "thickness" : 130
+    "radius" : 730,
+    "thickness" : 20,
+    "volfraction" : phi,
+    "scale" : (1-phi)
 }
 

sasmodels/models/core_shell_sphere.py

@@ -91 +91 @@
         @param thickness: shell thickness
     """
+    return (1,1)
     whole = 4.0 * pi / 3.0 * pow((radius + thickness), 3)
     core = 4.0 * pi / 3.0 * radius * radius * radius

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