Changeset a1b2471 in sasview for sansmodels/src/sans/models/CoreMultiShellModel.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/CoreMultiShellModel.py (modified) (1 diff)

sansmodels/src/sans/models/CoreMultiShellModel.py

@@ -129 +129 @@
                         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.
+        """
+        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_shell%s'% str(n)])"
+
+            # Right side of each shells
+            exec "r0 += self.params['thick_shell%s'% str(n)]"
+            r.append(r0)
+            exec "beta.append(self.params['sld_shell%s'% str(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):
         """