Changeset 2e83ff3 in sasview for DataLoader/test

Timestamp:

Sep 22, 2008 3:53:09 PM (16 years ago)

Author:

Mathieu Doucet <doucetm@…>

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:

1374174f

Parents:

55fd102

Message:

Added sector average (phi and q)

Location:

DataLoader/test

Files:

• sectorphi_testdata.txt (added)
• sectorq_testdata.txt (added)
• slabx_testdata.txt (added)
• slaby_testdata.txt (added)
• utest_averaging.py (modified) (1 diff)

DataLoader/test/utest_averaging.py

@@ -61 +61 @@
         self.assertAlmostEqual(ds, 0.011979881083557541, 4)
             
+    def test_slabX(self):
+        """
+            Test slab in X
+            The test data was not generated by IGOR.
+        """
+        from DataLoader.manipulations import SlabX
+        
+        r = SlabX(x_min=-.01, x_max=.01, y_min=-0.0002, y_max=0.0002, bin_width=0.0004)
+        r.fold = False
+        o = r(self.data)
+    
+        answer = Loader().load('slabx_testdata.txt')
+        for i in range(len(o.x)):
+            self.assertAlmostEqual(o.x[i], answer.x[i], 4)
+            self.assertAlmostEqual(o.y[i], answer.y[i], 4)
+            self.assertAlmostEqual(o.dy[i], answer.dy[i], 4)
+            
+    def test_slabY(self):
+        """
+            Test slab in Y
+            The test data was not generated by IGOR.
+        """
+        from DataLoader.manipulations import SlabY
+        
+        r = SlabY(x_min=.005, x_max=.01, y_min=-0.01, y_max=0.01, bin_width=0.0004)
+        r.fold = False
+        o = r(self.data)
+    
+        answer = Loader().load('slaby_testdata.txt')
+        for i in range(len(o.x)):
+            self.assertAlmostEqual(o.x[i], answer.x[i], 4)
+            self.assertAlmostEqual(o.y[i], answer.y[i], 4)
+            self.assertAlmostEqual(o.dy[i], answer.dy[i], 4)
+            
+    def test_sectorphi_full(self):
+        """
+            Test sector averaging I(phi)
+            When considering the whole azimuthal range (2pi), 
+            the answer should be the same as ring averaging.
+            The test data was not generated by IGOR.
+        """
+        from DataLoader.manipulations import SectorPhi
+        import math
+        
+        r = SectorPhi(r_min=.005, r_max=.01, phi_min=0, phi_max=math.pi*2.0)
+        r.nbins_phi = 20
+        o = r(self.data)
+    
+        answer = Loader().load('ring_testdata.txt')
+        for i in range(len(o.x)):
+            self.assertAlmostEqual(o.x[i], answer.x[i], 4)
+            self.assertAlmostEqual(o.y[i], answer.y[i], 4)
+            self.assertAlmostEqual(o.dy[i], answer.dy[i], 4)
+            
+    def test_sectorphi_quarter(self):
+        """
+            Test sector averaging I(phi)
+            The test data was not generated by IGOR.
+        """
+        from DataLoader.manipulations import SectorPhi
+        import math
+        
+        r = SectorPhi(r_min=.005, r_max=.01, phi_min=0, phi_max=math.pi/2.0)
+        r.nbins_phi = 20
+        o = r(self.data)
+    
+        answer = Loader().load('sectorphi_testdata.txt')
+        for i in range(len(o.x)):
+            self.assertAlmostEqual(o.x[i], answer.x[i], 4)
+            self.assertAlmostEqual(o.y[i], answer.y[i], 4)
+            self.assertAlmostEqual(o.dy[i], answer.dy[i], 4)
+            
+    def test_sectorq_full(self):
+        """
+            Test sector averaging I(q)
+            The test data was not generated by IGOR.
+        """
+        from DataLoader.manipulations import SectorQ
+        import math
+        
+        r = SectorQ(r_min=.005, r_max=.01, phi_min=0, phi_max=math.pi/2.0)
+        r.nbins_phi = 20
+        o = r(self.data)
+    
+        answer = Loader().load('sectorq_testdata.txt')
+        for i in range(len(o.x)):
+            self.assertAlmostEqual(o.x[i], answer.x[i], 4)
+            self.assertAlmostEqual(o.y[i], answer.y[i], 4)
+            self.assertAlmostEqual(o.dy[i], answer.dy[i], 4)
+            
 
 if __name__ == '__main__':