diff --git a/tests/regression_tests/dagmc/legacy/test.py b/tests/regression_tests/dagmc/legacy/test.py
index 963e73226e..1551256f9a 100644
--- a/tests/regression_tests/dagmc/legacy/test.py
+++ b/tests/regression_tests/dagmc/legacy/test.py
@@ -8,42 +8,48 @@ pytestmark = pytest.mark.skipif(
     not openmc.lib._dagmc_enabled(),
     reason="DAGMC CAD geometry is not enabled.")
 
+class DAGMCLegacyTest(PyAPITestHarness):
+
+    def _build_inputs(self):
+        model = openmc.model.Model()
+
+        # settings
+        model.settings.batches = 5
+        model.settings.inactive = 0
+        model.settings.particles = 100
+
+        source_box = openmc.stats.Box([-4, -4, -4],
+                                      [ 4,  4,  4])
+        source = openmc.Source(space=source_box)
+
+        model.settings.source = source
+
+        model.settings.dagmc = True
+
+        # tally
+        tally = openmc.Tally()
+        tally.scores = ['total']
+        tally.filters = [openmc.CellFilter(1)]
+        model.tallies = [tally]
+
+        # materials
+        u235 = openmc.Material(name="fuel")
+        u235.add_nuclide('U235', 1.0, 'ao')
+        u235.set_density('g/cc', 11)
+        u235.id = 40
+
+        water = openmc.Material(name="water")
+        water.add_nuclide('H1', 2.0, 'ao')
+        water.add_nuclide('O16', 1.0, 'ao')
+        water.set_density('g/cc', 1.0)
+        water.add_s_alpha_beta('c_H_in_H2O')
+        water.id = 41
+
+        mats = openmc.Materials([u235, water])
+        model.materials = mats
+
+        model.export_to_xml()
+
 def test_dagmc():
-    model = openmc.model.Model()
-
-    # settings
-    model.settings.batches = 5
-    model.settings.inactive = 0
-    model.settings.particles = 100
-
-    source_box = openmc.stats.Box([-4, -4, -4],
-                                  [ 4,  4,  4])
-    source = openmc.Source(space=source_box)
-
-    model.settings.source = source
-
-    model.settings.dagmc = True
-
-    # tally
-    tally = openmc.Tally()
-    tally.scores = ['total']
-    tally.filters = [openmc.CellFilter(1)]
-    model.tallies = [tally]
-
-    # materials
-    u235 = openmc.Material(name="fuel")
-    u235.add_nuclide('U235', 1.0, 'ao')
-    u235.set_density('g/cc', 11)
-    u235.id = 40
-
-    water = openmc.Material(name="water")
-    water.add_nuclide('H1', 2.0, 'ao')
-    water.add_nuclide('O16', 1.0, 'ao')
-    water.set_density('g/cc', 1.0)
-    water.add_s_alpha_beta('c_H_in_H2O')
-    water.id = 41
-
-    mats = openmc.Materials([u235, water])
-    model.materials = mats
-
-    model.export_to_xml()
+    harness = DAGMCLegacyTest('statepoint.5.h5')
+    harness.main()