Merge remote-tracking branch 'origin/feature/7452_ILL_SANS_Data_reduc…

…tion'
mantidproject · Apr 24, 2014 · 863b0a4 · 863b0a4
2 parents d008ce7 + 2d43339
commit 863b0a4
Show file tree

Hide file tree

Showing 8 changed files with 169 additions and 0 deletions.
diff --git a/Data/ILL/001420.nxs b/Data/ILL/001420.nxs
diff --git a/Data/ILL/001422.nxs b/Data/ILL/001422.nxs
diff --git a/Data/ILL/001425.nxs b/Data/ILL/001425.nxs
diff --git a/Data/ILL/001427.nxs b/Data/ILL/001427.nxs
diff --git a/Data/ILL/001428.nxs b/Data/ILL/001428.nxs
diff --git a/Data/ILL/001431.nxs b/Data/ILL/001431.nxs
diff --git a/Data/ILL/readme.txt b/Data/ILL/readme.txt
@@ -0,0 +1,54 @@
+1425 sample Latex
+1431 sample transmission
+1427 direct beam and transmission
+1420 B4C
+1422 empty cell
+1428 empty cell transmission
+
+
+SPREAD SHEET WRITTEN FOR D33 BY LAMP on Mon Mar 18 14:10:59 2013
+
+Lambda: 2.09 Angstrom.   Delta_L/Lambda:0.10
+
+Dead time correction is disabled.   Solid angle correction is disabled.
+MONITORS are used for normalisations. 
+
+Reduction set for 2*Pi radial integration.
+
+**************
+*WATER  RUNS.  Notes-> 
+**************
+     Water run   Cd/B4C run  W. cell run Thickness Cm Tr.Water run Tr.Wcell run
+  ------------ ------------ ------------ ------------ ------------ ------------
+
+
+Tr.Wbeam run : 
+<->Tr.Wvalue : 0.000
+<->Tr.Cvalue : 0.000
+Calibration: d33.cal
+
+
+**************
+*LOW    ANGLE.[Mask=425,437,429,440] [CxCy=430.400,431.740 dTh:0.0001]  Notes-> 
+*MIDDLE ANGLE.[CxCy=428.850,430.050 dTh:0.0001] [Mask=421,438,423,439]  Notes-> 
+*HIGH   ANGLE. Notes-> 
+**************
+               LOW    Q run MIDDLE Q run HIGH   Q run  Transm. run <-> Tr.value
+               ------------ ------------ ------------ ------------ ------------
+     Beam Tr.         1031         1427                      1427            -
+      Cd/B4C                       1420                         -            -
+     Cell Bg.                      1422                      1428        0.840
+
+
+**************
+*SAMPLE inputs.
+**************
+  LOW    Q run MIDDLE Q run HIGH   Q run  Transm. run Thickness Cm <-> Tr.value
+  ------------ ------------ ------------ ------------ ------------ ------------
+ ID1: J
+         1030         1425                      1431           .1        0.993
+ ID2: D2O
+                                                               .1             
+ ID3: S3
+
+
diff --git a/SystemTests/AnalysisTests/ILLD33Test.py b/SystemTests/AnalysisTests/ILLD33Test.py
@@ -0,0 +1,115 @@
+import stresstesting
+
+from mantid.api import mtd
+from mantid.simpleapi import SetupILLD33Reduction, SANSReduction,Rebin,SANSAzimuthalAverage1D
+
+import unittest
+
+class ILLD33SANSTest(unittest.TestCase):
+
+    prefix = "D33"
+
+    def tearDown(self):
+        for wsName in mtd.getObjectNames():
+            if wsName.startswith(self.prefix):
+                mtd.remove(wsName)
+
+    def test_all(self):
+
+        SetupILLD33Reduction(
+            # Beam center shouldn't work
+            #BeamCenterMethod="None",
+            MaskedDetectorList=[14709,14710,14711,14712,14713,14714,14715,14716,14717,14718,14719,
+                                14720,14721,14722,14723,14724,14725,14726,14727,14728,14729,14730,
+                                14731,14732,14733,14734,14735,14965,14966,14967,14968,14969,14970,
+                                14971,14972,14973,14974,14975,14976,14977,14978,14979,14980,14981,
+                                14982,14983,14984,14985,14986,14987,14988,14989,14990,14991,15221,
+                                15222,15223,15224,15225,15226,15227,15228,15229,15230,15231,15232,
+                                15233,15234,15235,15236,15237,15238,15239,15240,15241,15242,15243,
+                                15244,15245,15246,15247,15477,15478,15479,15480,15481,15482,15483,
+                                15484,15485,15486,15487,15488,15489,15490,15491,15492,15493,15494,
+                                15495,15496,15497,15498,15499,15500,15501,15502,15503,15733,15734,
+                                15735,15736,15737,15738,15739,15740,15741,15742,15743,15744,15745,
+                                15746,15747,15748,15749,15750,15751,15752,15753,15754,15755,15756,
+                                15757,15758,15759,15989,15990,15991,15992,15993,15994,15995,15996,
+                                15997,15998,15999,16000,16001,16002,16003,16004,16005,16006,16007,
+                                16008,16009,16010,16011,16012,16013,16014,16015,16245,16246,16247,
+                                16248,16249,16250,16251,16252,16253,16254,16255,16256,16257,16258,
+                                16259,16260,16261,16262,16263,16264,16265,16266,16267,16268,16269,
+                                16270,16271,16501,16502,16503,16504,16505,16506,16507,16508,16509,
+                                16510,16511,16512,16513,16514,16515,16516,16517,16518,16519,16520,
+                                16521,16522,16523,16524,16525,16526,16527,16757,16758,16759,16760,
+                                16761,16762,16763,16764,16765,16766,16767,16768,16769,16770,16771,
+                                16772,16773,16774,16775,16776,16777,16778,16779,16780,16781,16782,
+                                16783,17013,17014,17015,17016,17017,17018,17019,17020,17021,17022,
+                                17023,17024,17025,17026,17027,17028,17029,17030,17031,17032,17033,
+                                17034,17035,17036,17037,17038,17039,17269,17270,17271,17272,17273,
+                                17274,17275,17276,17277,17278,17279,17280,17281,17282,17283,17284,
+                                17285,17286,17287,17288,17289,17290,17291,17292,17293,17294,17295,
+                                17525,17526,17527,17528,17529,17530,17531,17532,17533,17534,17535,
+                                17536,17537,17538,17539,17540,17541,17542,17543,17544,17545,17546,
+                                17547,17548,17549,17550,17551],
+            BeamCenterMethod="DirectBeam",
+            BeamCenterFile='ILL/001427.nxs',
+            Normalisation="Timer",
+            DarkCurrentFile= 'ILL/001420.nxs',
+            TransmissionMethod="DirectBeam",
+            TransmissionSampleDataFile= 'ILL/001431.nxs',
+            TransmissionEmptyDataFile= 'ILL/001427.nxs',
+            BckTransmissionEmptyDataFile= 'ILL/001427.nxs',
+            TransmissionBeamRadius = 3,
+            TransmissionUseSampleDC=False,
+            BackgroundFiles='ILL/001422.nxs',
+            BckTransmissionSampleDataFile='ILL/001428.nxs',
+            DoAzimuthalAverage=False,    
+            Do2DReduction=False,
+            ComputeResolution=True,
+            ReductionProperties=self.prefix + "props")
+
+        output=SANSReduction(Filename='ILL/001425.nxs', ReductionProperties=self.prefix + "props",
+                             OutputWorkspace=self.prefix + "out")
+        Rebin(InputWorkspace=self.prefix + 'out',OutputWorkspace=self.prefix + 'out_rebin',
+              Params='4,0.1,15')
+        SANSAzimuthalAverage1D(InputWorkspace=self.prefix + 'out_rebin',Binning='0.001,0.0002,0.03',
+                               OutputWorkspace=self.prefix + 'final')        
+
+        # Check some data
+        wsOut = mtd[self.prefix + 'out']
+        self.assertEqual(wsOut.getNumberHistograms(), 65538)
+        wsOut = mtd[self.prefix + 'out_rebin']
+        self.assertEqual(wsOut.getNumberHistograms(), 65538)
+        wsOut = mtd[self.prefix + 'final']
+        self.assertEqual(wsOut.getNumberHistograms(), 1)
+
+
+
+    #================== Failure cases ================================
+
+    # TODO
+
+
+
+
+#====================================================================================
+
+class ILLD33Test(stresstesting.MantidStressTest):
+
+    def requiredMemoryMB(self):
+        """Set a limit of 2.5Gb to avoid 32-bit environment"""
+        return 2500
+
+    def runTest(self):
+        self._success = False
+        # Custom code to create and run this single test suite
+        suite = unittest.TestSuite()
+        suite.addTest( unittest.makeSuite(ILLD33SANSTest, "test") )
+        runner = unittest.TextTestRunner()
+        # Run using either runner
+        res = runner.run(suite)
+        if res.wasSuccessful():
+            self._success = True 
+        else:
+            self._success = False
+
+    def validate(self):
+        return self._success