Changeset 1d67243 in sasview for sansview/perspectives/fitting/sldPanel.py

Timestamp:

Aug 20, 2009 12:37:45 PM (15 years ago)

Author:

Gervaise Alina <gervyh@…>

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:

6c5239c

Parents:

8e36cdd

Message:

sld calculator redesigned

File:

• sansview/perspectives/fitting/sldPanel.py (modified) (11 diffs)

sansview/perspectives/fitting/sldPanel.py

@@ -1 +1 @@
 """
-    This module intends to compute the neutron scattering length density of molecule
+    This module provide GUI for the neutron scattering length density calculator
     @author: Gervaise B. Alina
 """
@@ -7 +7 @@
 import sys
 
-import periodictable
-from periodictable import formula
-from periodictable.xsf import xray_energy, xray_sld_from_atoms
-from periodictable.constants import avogadro_number
-import periodictable.nsf
-neutron_sld_from_atoms= periodictable.nsf.neutron_sld_from_atoms 
-
 from sans.guiframe.utils import format_number, check_float
 from sans.guicomm.events import StatusEvent  
-
+from sldCalculator import SldCalculator
 
 _BOX_WIDTH = 76
 _STATICBOX_WIDTH = 350
 _SCALE = 1e-6
-_DEFAULT_WAVELENGTH = 1.798
-
 
 class SldPanel(wx.Panel):
@@ -32 +23 @@
         # Object that receive status event
         self.base= base
+        self.calculator = SldCalculator()
+        self.wavelength = self.calculator.wavelength
+        
         self._do_layout()
         self.SetAutoLayout(True)
@@ -58 +52 @@
         wavelength_txt = wx.StaticText(self, -1, 'Wavelength (A)')
         self.wavelength_ctl = wx.TextCtrl(self, -1, size=(_BOX_WIDTH,-1))
+        self.wavelength_ctl.SetValue(str(self.wavelength))
         iy = 0
         ix = 0
@@ -231 +226 @@
         #---------layout----------------
         vbox  = wx.BoxSizer(wx.VERTICAL)
-
         vbox.Add(sizer1)
         vbox.Add(sizer2)
@@ -237 +231 @@
         vbox.Fit(self) 
         self.SetSizer(vbox)
+        
         
     def check_inputs(self):
@@ -250 +245 @@
         self.wavelength= self.wavelength_ctl.GetValue()
         if self.wavelength.lstrip().rstrip()=="":
-            self.wavelength = _DEFAULT_WAVELENGTH
+            self.wavelength = self.calculator.wavelength
         else:
             if check_float(self.wavelength_ctl):
@@ -269 +264 @@
         return flag 
         
+        
     def onHelp(self, event):
         """
@@ -294 +290 @@
                 #get ready to compute
                 try:
-                    self.new_formula = formula(self.formulata_text, density= self.density)
-                    atom = self.new_formula.atoms
+                    self.calculator.setValue(self.formulata_text,self.density,self.wavelength)
                 except:
                     if self.base !=None:
@@ -303 +298 @@
                         raise
                     return
+                
                 # Compute the Cu SLD
-                Cu_reel, Cu_im = self.calculateXRaySld( "Cu", density= self.density,
-                                                        user_formula= self.new_formula)
+                Cu_reel, Cu_im = self.calculator.calculateXRaySld( "Cu")
                 self.cu_ka_sld_reel_ctl.SetValue(format_number(Cu_reel*_SCALE))
                 self.cu_ka_sld_im_ctl.SetValue(format_number(Cu_im*_SCALE))
+                
                 # Compute the Mo SLD
-                Mo_reel, Mo_im = self.calculateXRaySld( "Mo", density= self.density,
-                                                        user_formula= self.new_formula)
+                Mo_reel, Mo_im = self.calculator.calculateXRaySld( "Mo")
                 self.mo_ka_sld_reel_ctl.SetValue(format_number(Mo_reel*_SCALE))
                 self.mo_ka_sld_im_ctl.SetValue(format_number(Mo_im*_SCALE))
                
-                coh,absorp,inc= self.calculateNSld(self.density, wavelength= self.wavelength,
-                                                     user_formula= self.new_formula)
-                #Don't know if value is return in cm or  cm^(-1).assume return in cm
-                # to match result of neutron inc of Alan calculator
-                inc= inc*1/10
-                #Doesn't match result of Alan calculator for absorption factor of 2
-                #multiplication of 100 is going around
-                absorp= absorp *2*100
-                volume= (self.new_formula.mass /self.density)/avogadro_number*1.0e24
-                #im: imaginary part of neutron SLD
-                im=0
-                for el, count in atom.iteritems():
-                    if el.neutron.b_c_i !=None:
-                        im += el.neutron.b_c_i*count 
-                im = im/volume
-                
+                coh,absorp,inc= self.calculator.calculateNSld()
+                im = self.calculator.absorptionIm()
+                length = self.calculator.computeLength()
+                # Neutron SLD
                 self.neutron_sld_reel_ctl.SetValue(format_number(coh*_SCALE))
                 self.neutron_sld_im_ctl.SetValue(format_number(im*_SCALE))
                 self.neutron_inc_ctl.SetValue(format_number(inc ))
                 self.neutron_abs_ctl.SetValue(format_number(absorp))
-                #Don't know if value is return in cm or  cm^(-1).assume return in cm
-                # to match result of neutron inc of Alan calculator
-                length= (coh+ absorp+ inc)/volume
+                # Neutron length
                 self.neutron_length_ctl.SetValue(format_number(length))
+                # display wavelength
+                self.wavelength_ctl.SetValue(str(self.wavelength))
         except:
             if self.base !=None:
@@ -346 +329 @@
             return   
 
-    def calculateXRaySld(self, element, density,user_formula):
-        """
-            Get an element and compute the corresponding SLD for a given formula
-            @param element:  elementis a string of existing atom
-            @param formula: molecule enters by the user
-        """
-        try:
-            myformula = formula(str (element))
-            if len(myformula.atoms)!=1:
-                return 
-            element= myformula.atoms.keys()[0] 
-            energy = xray_energy(element.K_alpha)
-            atom = user_formula.atoms
-            atom_reel, atom_im = xray_sld_from_atoms( atom,
-                                                  density=density,
-                                                  energy= energy )
-            return atom_reel, atom_im
-        except:
-            if self.base !=None:
-                  msg= "SLD Calculator: %s" % (sys.exc_value)
-                  wx.PostEvent(self.base, StatusEvent(status= msg ))
-            else:
-                raise
-            return   
-        
-    def calculateNSld(self,density,wavelength,user_formula ):
-        """
-            Compute the neutron SLD for a given molecule
-            @return absorp: absorption
-            @return coh: coherence cross section
-            @return inc: incoherence cross section
-          
-        """
-        if density ==0:
-            raise ZeroDivisionError,"integer division or modulo by zero for density"
-            return 
-        atom = user_formula.atoms
-        coh,absorp,inc = neutron_sld_from_atoms(atom,density,wavelength)
-
-        return coh,absorp,inc
+   
+        
+