diff --git a/Code/Mantid/Testing/Data/SystemTest/REF_L_119690_event.nxs.md5 b/Code/Mantid/Testing/Data/SystemTest/REF_L_119690_event.nxs.md5
new file mode 100644
index 0000000000000000000000000000000000000000..b6a45fd7261ef11bd82a7a9940b308d54ad4c0a8
--- /dev/null
+++ b/Code/Mantid/Testing/Data/SystemTest/REF_L_119690_event.nxs.md5
@@ -0,0 +1 @@
+79694adb3d0d8a7573e0d9b0ff210edc
diff --git a/Code/Mantid/Testing/Data/SystemTest/REF_L_119814_event.nxs.md5 b/Code/Mantid/Testing/Data/SystemTest/REF_L_119814_event.nxs.md5
new file mode 100644
index 0000000000000000000000000000000000000000..e25141c4c3c67967dcd18a50e4b17f0d0e85ae14
--- /dev/null
+++ b/Code/Mantid/Testing/Data/SystemTest/REF_L_119814_event.nxs.md5
@@ -0,0 +1 @@
+43d80c791f44419e242ec2e5ff34ec46
diff --git a/Code/Mantid/Testing/Data/SystemTest/directBeamDatabaseFall2014_IPTS_11601_2.cfg.md5 b/Code/Mantid/Testing/Data/SystemTest/directBeamDatabaseFall2014_IPTS_11601_2.cfg.md5
new file mode 100644
index 0000000000000000000000000000000000000000..5c1e010c0bbeae3ea7efc981bbad3e2e81d6ccb2
--- /dev/null
+++ b/Code/Mantid/Testing/Data/SystemTest/directBeamDatabaseFall2014_IPTS_11601_2.cfg.md5
@@ -0,0 +1 @@
+0f14c83d41bad464038778d2dd02424d
diff --git a/Code/Mantid/Testing/SystemTests/tests/analysis/REFLReduction.py b/Code/Mantid/Testing/SystemTests/tests/analysis/REFLReduction.py
new file mode 100644
index 0000000000000000000000000000000000000000..adae2d6a95035c91847b6f0e01f0ddf5745427d7
--- /dev/null
+++ b/Code/Mantid/Testing/SystemTests/tests/analysis/REFLReduction.py
@@ -0,0 +1,45 @@
+#pylint: disable=no-init
+import stresstesting
+from mantid import *
+
+from mantid.simpleapi import *
+
+class REFLReduction(stresstesting.MantidStressTest):
+    def runTest(self):
+        #TODO: The reduction algorithm should not require an absolute path
+        scaling_factor_file = FileFinder.getFullPath("directBeamDatabaseFall2014_IPTS_11601_2.cfg")
+        
+        RefLReduction(RunNumbers=[119814],
+                      NormalizationRunNumber=119690,
+                      SignalPeakPixelRange=[154, 166],
+                      SubtractSignalBackground=True,
+                      SignalBackgroundPixelRange=[151, 169],
+                      NormFlag=True,
+                      NormPeakPixelRange=[154, 160],
+                      NormBackgroundPixelRange=[151, 163],
+                      SubtractNormBackground=True,
+                      LowResDataAxisPixelRangeFlag=True,
+                      LowResDataAxisPixelRange=[99, 158],
+                      LowResNormAxisPixelRangeFlag=True,
+                      LowResNormAxisPixelRange=[98, 158],
+                      TOFRange=[29623.0, 42438.0],
+                      IncidentMediumSelected='2InDiamSi',
+                      GeometryCorrectionFlag=False,
+                      QMin=0.005,
+                      QStep=0.01,
+                      AngleOffset=0.009,
+                      AngleOffsetError=0.001,
+                      ScalingFactorFile=scaling_factor_file,
+                      SlitsWidthFlag=True,
+                      OutputWorkspace='reflectivity_119814')
+
+    def validate(self):
+        # Be more tolerant with the output, mainly because of the errors.
+        # The following tolerance check the errors up to the third digit.
+        self.tolerance = 0.25
+        self.disableChecking.append('Instrument')
+        self.disableChecking.append('Sample')
+        self.disableChecking.append('SpectraMap')
+        self.disableChecking.append('Axes')
+        return "reflectivity_119814", 'REFL_119814_combined_data.nxs'
+
diff --git a/Code/Mantid/Testing/SystemTests/tests/analysis/reference/REFL_119814_combined_data.nxs.md5 b/Code/Mantid/Testing/SystemTests/tests/analysis/reference/REFL_119814_combined_data.nxs.md5
new file mode 100644
index 0000000000000000000000000000000000000000..ad207d1de237a2322bc8cd16095cc6704ccb3f6f
--- /dev/null
+++ b/Code/Mantid/Testing/SystemTests/tests/analysis/reference/REFL_119814_combined_data.nxs.md5
@@ -0,0 +1 @@
+f7949eee903277da4a6c1fa957500546