Changeset bcfc969 in sasview

Timestamp:

Aug 11, 2011 9:04:56 AM (13 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:

c0a30a24

Parents:

0f282a5

Message:

more accurate way to cal averaging sigmas

File:

• sanscalculator/src/sans/calculator/resolution_calculator.py (modified) (2 diffs)

sanscalculator/src/sans/calculator/resolution_calculator.py

@@ -166 +166 @@
             #    continue
             image = out
-            # set sigmas
-            sigma_1 += sig1_list[ind] * self.intensity
-            sigma_r += sigr_list[ind] * self.intensity
-            sigma_2 += sig2_list[ind] * self.intensity
-            sigma1d += sigma1d_list[ind] * self.intensity
+            # set variance as sigmas
+            sigma_1 += sig1_list[ind] *  sig1_list[ind] * self.intensity
+            sigma_r += sigr_list[ind] * sigr_list[ind] * self.intensity
+            sigma_2 += sig2_list[ind] * sig2_list[ind] * self.intensity
+            sigma1d += sigma1d_list[ind] * sigma1d_list[ind] * self.intensity
             total_intensity += self.intensity
+        
         if total_intensity != 0:
+            # average variance
             image_out = image / total_intensity
             sigma_1 = sigma_1 / total_intensity
@@ -179 +181 @@
             sigma1d = sigma1d / total_intensity
             # set sigmas
-            self.sigma_1 = sigma_1
-            self.sigma_lamd = sigma_r
-            self.sigma_2 = sigma_2
-            self.sigma_1d = sigma1d
+            self.sigma_1 = sqrt(sigma_1)
+            self.sigma_lamd = sqrt(sigma_r)
+            self.sigma_2 = sqrt(sigma_2)
+            self.sigma_1d = sqrt(sigma1d)
             # rescale 
             max_im_val = 1 #image_out.max()