diff --git a/tests/regression_tests/photon_production/inputs_true.dat b/tests/regression_tests/photon_production/inputs_true.dat
index 1d43c89137..f82b4a6a78 100644
--- a/tests/regression_tests/photon_production/inputs_true.dat
+++ b/tests/regression_tests/photon_production/inputs_true.dat
@@ -49,17 +49,20 @@
   </tally>
   <tally id="2">
     <filters>2</filters>
-    <scores>total</scores>
+    <nuclides>Al27 total</nuclides>
+    <scores>total heating</scores>
     <estimator>tracklength</estimator>
   </tally>
   <tally id="3">
     <filters>2</filters>
-    <scores>total</scores>
+    <nuclides>Al27 total</nuclides>
+    <scores>total heating</scores>
     <estimator>collision</estimator>
   </tally>
   <tally id="4">
     <filters>2</filters>
-    <scores>total</scores>
+    <nuclides>Al27 total</nuclides>
+    <scores>total heating</scores>
     <estimator>analog</estimator>
   </tally>
 </tallies>
diff --git a/tests/regression_tests/photon_production/results_true.dat b/tests/regression_tests/photon_production/results_true.dat
index 7df7e5dbcf..a7725576c8 100644
--- a/tests/regression_tests/photon_production/results_true.dat
+++ b/tests/regression_tests/photon_production/results_true.dat
@@ -2,11 +2,29 @@ tally 1:
 9.403000E-01
 8.841641E-01
 tally 2:
+1.819886E-01
+3.311985E-02
+1.960159E+05
+3.842224E+10
 8.281718E-01
 6.858685E-01
+1.960159E+05
+3.842224E+10
 tally 3:
+1.799069E-01
+3.236650E-02
+1.945225E+05
+3.783901E+10
 8.242000E-01
 6.793056E-01
+1.945225E+05
+3.783901E+10
 tally 4:
+2.308000E-01
+5.326864E-02
+1.988563E+05
+3.954384E+10
 8.242000E-01
 6.793056E-01
+1.988612E+05
+3.954578E+10
diff --git a/tests/regression_tests/photon_production/test.py b/tests/regression_tests/photon_production/test.py
index 47215fff30..c6a1fc3908 100644
--- a/tests/regression_tests/photon_production/test.py
+++ b/tests/regression_tests/photon_production/test.py
@@ -51,20 +51,23 @@ class SourceTestHarness(PyAPITestHarness):
         current_tally = openmc.Tally()
         current_tally.filters = [surface_filter, particle_filter]
         current_tally.scores = ['current']
-        total_tally_tracklength = openmc.Tally()
-        total_tally_tracklength.filters = [particle_filter]
-        total_tally_tracklength.scores = ['total']
-        total_tally_tracklength.estimator = 'tracklength'
-        total_tally_collision = openmc.Tally()
-        total_tally_collision.filters = [particle_filter]
-        total_tally_collision.scores = ['total']
-        total_tally_collision.estimator = 'collision'
-        total_tally_analog = openmc.Tally()
-        total_tally_analog.filters = [particle_filter]
-        total_tally_analog.scores = ['total']
-        total_tally_analog.estimator = 'analog'
-        tallies = openmc.Tallies([current_tally, total_tally_tracklength,
-                                  total_tally_collision, total_tally_analog])
+        tally_tracklength = openmc.Tally()
+        tally_tracklength.filters = [particle_filter]
+        tally_tracklength.scores = ['total', 'heating']
+        tally_tracklength.nuclides = ['Al27', 'total']
+        tally_tracklength.estimator = 'tracklength'
+        tally_collision = openmc.Tally()
+        tally_collision.filters = [particle_filter]
+        tally_collision.scores = ['total', 'heating']
+        tally_collision.nuclides = ['Al27', 'total']
+        tally_collision.estimator = 'collision'
+        tally_analog = openmc.Tally()
+        tally_analog.filters = [particle_filter]
+        tally_analog.scores = ['total', 'heating']
+        tally_analog.nuclides = ['Al27', 'total']
+        tally_analog.estimator = 'analog'
+        tallies = openmc.Tallies([current_tally, tally_tracklength,
+                                  tally_collision, tally_analog])
         tallies.export_to_xml()
 
     def _get_results(self):