Changeset 997131a in sasview for sansview/perspectives/fitting/fitpage.py

Timestamp:

Apr 16, 2009 9:20:40 AM (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:

847091f

Parents:

d1e0473

Message:

plotting with smearer object fixed

File:

• sansview/perspectives/fitting/fitpage.py (modified) (9 diffs)

sansview/perspectives/fitting/fitpage.py

@@ -22 +22 @@
 from basepage import BasicPage
 from basepage import PageInfoEvent
-
+from DataLoader.qsmearing import smear_selection
 
 class FitPage(BasicPage):
@@ -48 +48 @@
         self._fill_range_sizer() 
         self._on_select_model(event=None)
-    
+        if self.data !=None:
+            self.smearer = smear_selection( self.data )
+            
         ## to update the panel according to the fit engine type selected
         self.Bind(EVT_FITTER_TYPE,self._on_engine_change)
@@ -452 +454 @@
         self.set_model_param_sizer(self.model)
         try:
-            self.compute_chisqr()
+            temp_smear= None
+            if self.enable_smearer.GetValue():
+                temp_smear= self.smearer
+            self.compute_chisqr(temp_smear)
         except:
             ## error occured on chisqr computation
@@ -603 +608 @@
             are compute when fitting
         """
-        smear =None
-        msg=""
+        
         if self.enable_smearer.GetValue():
-            from DataLoader.qsmearing import smear_selection
-            smear = smear_selection( self.data )
+            msg=""
+            temp_smearer= self.smearer
             if hasattr(self.data,"dxl"):
                 msg= ": Resolution smearing parameters"
             if hasattr(self.data,"dxw"):
                 msg= ": Slit smearing parameters"
-            if smear ==None:
+            if self.smearer ==None:
                 wx.PostEvent(self.manager.parent, StatusEvent(status=\
                             "Data contains no smearing information"))
@@ -618 +622 @@
                 wx.PostEvent(self.manager.parent, StatusEvent(status=\
                             "Data contains smearing information %s"%msg))
+        else:
+            temp_smearer = None
         ## set smearing value whether or not the data contain the smearing info
-        self.manager.set_smearer(smear, qmin= float(self.qmin_x),
+        self.manager.set_smearer(temp_smearer, qmin= float(self.qmin_x),
                                       qmax= float(self.qmax_x)) 
         ##Calculate chi2
-        self.compute_chisqr()  
+        self.compute_chisqr(smearer= temp_smearer)  
         ## save the state enable smearing
         self.save_current_state()
@@ -668 +674 @@
     
         
-    def compute_chisqr(self):
+    def compute_chisqr(self , smearer=None):
         """ 
             compute chi square given a model and data 1D and set the value
@@ -692 +698 @@
                         dy= numpy.asarray(self.data.dy)
                     dy[dy==0]=1
+                    
                     if self.qmin_x==None and self.qmax_x==None: 
                         fx =numpy.asarray([self.model.run(v) for v in x])
+                        if smearer!=None:
+                            fx= smearer(fx)
                         temp=(y - fx)/dy
                         res= temp*temp
@@ -699 +708 @@
                         idx = (x>= self.qmin_x) & (x <=self.qmax_x)
                         fx = numpy.asarray([self.model.run(item)for item in x[idx ]])
+                        if smearer!=None:
+                            fx= smearer(fx)
                         temp=(y[idx] - fx)/dy[idx]
                         res= temp*temp
@@ -708 +719 @@
                     self.tcChi.SetLabel(format_number(math.fabs(sum/ len(res))))
             except:
-                wx.PostEvent(self.parent.GrandParent, StatusEvent(status=\
-                            "Chisqr cannot be compute: %s"% sys.exc_value))
-                return 
+                raise
+                #wx.PostEvent(self.parent.GrandParent, StatusEvent(status=\
+                #            "Chisqr cannot be compute: %s"% sys.exc_value))
+                #return