Changeset b6422c7 in sasmodels for explore/beta/sasfit_compare_new.py

Timestamp:

Mar 25, 2019 1:06:41 PM (6 years ago)

Author:

richardh

Branches:

master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

0507e09

Parents:

2fe39d1

Message:

more edits to sasfit_compare_new, vesicle and hollow cylinder improving

File:

• explore/beta/sasfit_compare_new.py (modified) (14 diffs)

explore/beta/sasfit_compare_new.py

@@ -21 +21 @@
 #Used to calculate F(q) for the cylinder, sphere, ellipsoid models
 # RKH adding vesicle and hollow_cylinder to test sasview special cases of ER and VR
-# BEWARE there are issues here if you run this from python3 (i.e. qt5 version of sasview), so do "source activate sasview"
+# There were issues here from python3 (i.e. qt5 version of sasview),fixed by Paul K (else do "source activate sasview")
 def sas_sinx_x(x):
     with np.errstate(all='ignore'):
@@ -275 +275 @@
     if radius_effective is None:
         radius_effective = radius
-        print("radius_effective  ",radius_effective)
         average_volume = total_volume/total_weight
     if norm == 'sasfit':
@@ -285 +284 @@
         F2 *= 1e-12
         IQD = F2/average_volume*1e8*volfraction
+    print("in sphere_r : radius_effective  ",radius_effective,"  volfraction",volfraction)
     SQ = hardsphere_simple(q, radius_effective, volfraction)
     beta = F1**2/F2
@@ -332 +332 @@
     if radius_effective is None:
         radius_effective = radius_eff/total_weight
-        print("radius_effective  ",radius_effective)
+        print("in vesicle_pe : radius_effective  ",radius_effective)
     if norm == 'sasfit':
         IQD = F2
@@ -348 +348 @@
         F2 *= 1e-12
         IQD = F2/average_volume*1e8*volfraction
-    # RKH SORT THIS OUT WHEN HAVE NEW MODEL
-    vfff=volfraction*(30./20.)**3
-    print("radius_effective, fudged volfaction for S(Q) ",radius_effective,vfff)
+    # RKH SORT THIS OUT WHEN HAVE NEW MODEL  Vtot/Vshell = (R+T)^3/ ( (R+T)^3 -R^3 )
+    vfff=volfraction*( 1.0/(1.0 - ( (radius/(radius+thickness))**3 )))
+    print("in vesicle_pe : radius_effective, fudged volfraction for S(Q) ",radius_effective,vfff)
     SQ = hardsphere_simple(q, radius_effective, vfff)
     beta = F1**2/F2
@@ -359 +359 @@
 #
 #polydispersity for hollow_cylinder
-def hollow_cylinder_pe(q,radius=20, thickness=10, sld=4, sld_solvent=1, volfraction=0.15,
-        radius_pd=0.1, thickness_pd=0.2, length_pd=0.05, radius_pd_type="gaussian", thickness_pd_type="gaussian",
-        radius_effective=None, background=0, scale=1,
+def hollow_cylinder_pe(q,radius=20, thickness=10, length=80, sld=4, sld_solvent=1, volfraction=0.15,
+        radius_pd=0.1, thickness_pd=0.2, length_pd=0.05, radius_pd_type="gaussian", length_pd_type="gaussian", 
+        thickness_pd_type="gaussian", radius_effective=None, background=0, scale=1,
                  norm='sasview'):
     if norm not in ['sasview', 'sasfit', 'yun']:
@@ -374 +374 @@
         T_val, T_prob = schulz_distribution(thickness, thickness_pd*thickness, 0, inf)
     if length_pd_type == 'gaussian':
-        T_val, T_prob = gaussian_distribution(length, length_pd*length, 0, inf)
+        L_val, L_prob = gaussian_distribution(length, length_pd*length, 0, inf)
     elif length_pd_type == 'schulz':
         L_val, L_prob = schulz_distribution(length, length_pd*length, 0, inf)
@@ -402 +402 @@
     if radius_effective is None:
         radius_effective = radius_eff/total_weight
-        print("radius_effective  ",radius_effective)
+        print("in hollow_cylinder : radius_effective  ",radius_effective)
     if norm == 'sasfit':
         IQD = F2
@@ -420 +420 @@
 # RKH SORT THIS OUT WHEN HAVE NEW MODEL
     vfff=volfraction
-    print("radius_effective, volfaction for S(Q) ",radius_effective,vfff)
+    print("in hollow_cylinder : radius_effective, volfaction for S(Q) ",radius_effective,vfff)
     SQ = hardsphere_simple(q, radius_effective, vfff)
     beta = F1**2/F2
@@ -548 +548 @@
     #Iq = None#= Sq = None
     r = dict(I._kernel.results())
-    #print(r)
+    print(r)
     return Theory(Q=q, F1=None, F2=None, P=Pq, S=None, I=None, Seff=r["S_eff(Q)"], Ibeta=Iq)
 
@@ -591 +591 @@
         radius_pd=.1, radius_pd_type=pd_type,
         volfraction=0.15,
-        radius_effective=12.59921049894873,  # equivalent average sphere radius
+        radius_effective=20  # equivalent average sphere radius
+#       NOTE sasview computes its own radius_effective in "target" (the print(r) at end of sasmodels_theory will reveal its value),
+#       this one is only used locally for "actual"
         )
     target = sasmodels_theory(q, model, effective_radius_mode=0, structure_factor_mode=1, **pars)
@@ -604 +606 @@
     q = np.logspace(-5, 0, 250)
     model = 'vesicle'
-    print("F2 seems OK, but big issues with S(Q), wo work in progress")
-# NOTE parameers for vesicle model are soon to be changed to remove "volfraction"
+    print("F2 seems OK, but big issues with S(Q), so work in progress")
+# NOTE parameters for vesicle model are soon to be changed to remove "volfraction" (though still there in 5.0)
     pars = dict(
         radius=20, thickness=10, sld=4, sld_solvent=1, volfraction=0.03,
         background=0,
-        radius_pd=0.1, thickness_pd=0.0, radius_pd_type=pd_type, thickness_pd_type=pd_type
-        #,radius_effective=12.59921049894873,  # equivalent average sphere radius
-        )
+        radius_pd=0.01, thickness_pd=0.01, radius_pd_type=pd_type, thickness_pd_type=pd_type,
+        radius_effective=30. ) 
+        # equivalent average sphere radius for local "actual" to match what sasview uses, use None to compute average outer radius here,
+
     target = sasmodels_theory(q, model, effective_radius_mode=0, structure_factor_mode=1, **pars)
     actual = vesicle_pe(q, norm='sasview', **pars)
@@ -629 +632 @@
         radius=20, thickness=10, length=80, sld=4, sld_solvent=1,
         background=0,
-        radius_pd=0.1, thickness_pd=0.0, length_pd=0.0, radius_pd_type=pd_type, thickness_pd_type=pd_type, length_pd_type=pd_type
-        #,radius_effective=12.59921049894873,  # equivalent average sphere radius
-        )
+        radius_pd=0.1, thickness_pd=0.0, length_pd=0.0, radius_pd_type=pd_type, thickness_pd_type=pd_type, length_pd_type=pd_type,
+        radius_effective=40.687)  
+        # equivalent average sphere radius for local "actual" to match what sasview uses
     target = sasmodels_theory(q, model, effective_radius_mode=0, structure_factor_mode=1, **pars)
     actual = hollow_cylinder_pe(q, norm='sasview', **pars)
@@ -651 +654 @@
         volfraction=0.15,
         radius_effective=270.7543927018,
-        #radius_effective=12.59921049894873,
+# if change radius_effective to some other value, the S(Q) from sasview does not agree
         )
     target = sasmodels_theory(q, model, effective_radius_mode=0, structure_factor_mode=1, **pars)