Changeset a1b2471 in sasview for sansmodels/src/sans/models/SphereExpShellModel.py

Timestamp:

Sep 21, 2010 5:44:24 PM (14 years ago)

Author:

Jae Cho <jhjcho@…>

Branches:

master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

db08737

Parents:

60c320b

Message:

added sld plot for onion model and etc…

File:

• sansmodels/src/sans/models/SphereExpShellModel.py (modified) (2 diffs)

sansmodels/src/sans/models/SphereExpShellModel.py

@@ -42 +42 @@
         for item in self.model.orientation_params:
             self.orientation_params.append(item)
-                
+        self.getProfile()        
     def _clone(self, obj):
         """
@@ -133 +133 @@
                             self.model.setParam(key, value)
                         except: pass
-
-   
+    def getProfile(self):
+        """
+        Get SLD profile 
+        
+        : return: (r, beta) where r is a list of radius of the transition points
+                beta is a list of the corresponding SLD values 
+        : Note: This works only for func_shell# = 2.
+        """
+        # max_pts for each shells
+        max_pts = 10
+        r = []
+        beta = []
+        # for core at r=0
+        r.append(0)
+        beta.append(self.params['sld_core'])
+        # for core at r=rad_core
+        r.append(self.params['rad_core'])
+        beta.append(self.params['sld_core'])
+        
+        # for shells
+        for n in range(1,self.n_shells+1):
+            # Left side of each shells
+            r0 = r[len(r)-1]            
+            r.append(r0)
+            exec "beta.append(self.params['sld_in_shell%s'% str(n)])"
+            
+            exec "A = self.params['A_shell%s'% str(n)]"
+            if A ==0:
+                # Right side of each shells
+                exec "r0 += self.params['thick_shell%s'% str(n)]"
+                r.append(r0)
+                exec "beta.append(self.params['sld_in_shell%s'% str(n)])"
+            else:
+                from math import exp
+                rn = r0
+                for n_sub in range(0,max_pts):
+                    # Right side of each sub_shells
+                    exec "rn += self.params['thick_shell%s'% str(n)]/10.0"
+                    r.append(rn)
+                    exec "slope = (self.params['sld_out_shell%s'% str(n)] \
+                                        -self.params['sld_in_shell%s'% str(n)]) \
+                                        /(exp(self.params['A_shell%s'% str(n)])-1)"
+                    exec "const = (self.params['sld_in_shell%s'% str(n)]-slope)"
+                    exec "beta_n = slope*exp((self.params['A_shell%s'% str(n)]*(rn-r0) \
+                                        /self.params['thick_shell%s'% str(n)])) + const"
+                    beta.append(beta_n)
+            
+        # for solvent
+        r0 = r[len(r)-1]            
+        r.append(r0)
+        beta.append(self.params['sld_solv'])
+        r_solv = 5*r0/4
+        r.append(r_solv)
+        beta.append(self.params['sld_solv'])
+        
+        return r, beta
+    
     def setParam(self, name, value):
         """