diff --git a/src/tally.F90 b/src/tally.F90
index 94143184ae..2c77119173 100644
--- a/src/tally.F90
+++ b/src/tally.F90
@@ -482,6 +482,12 @@ contains
               score = micro_xs(i_nuclide) % fission * nuclides(i_nuclide) % &
                    nu(E, EMISSION_PROMPT) * atom_density * flux
           else
+            ! make sure the correct energy is used
+            if (t % estimator == ESTIMATOR_TRACKLENGTH) then
+              E = p % E
+            else
+              E = p % last_E
+            end if
 
             score = ZERO
 
diff --git a/tests/test_mgxs_library_mesh/results_true.dat b/tests/test_mgxs_library_mesh/results_true.dat
index e3d7ca6473..03019cfd3d 100644
--- a/tests/test_mgxs_library_mesh/results_true.dat
+++ b/tests/test_mgxs_library_mesh/results_true.dat
@@ -126,7 +126,7 @@
 3      2  2  1        1   total  3.799163e-07  1.806470e-07
   mesh 1       group in nuclide      mean std. dev.
        x  y  z                                     
-0      1  1  1        1   total  0.025735  0.002840
-1      1  2  1        1   total  0.028773  0.006349
-2      2  1  1        1   total  0.022306  0.004010
-3      2  2  1        1   total  0.024549  0.009379
+0      1  1  1        1   total  0.025920  0.002893
+1      1  2  1        1   total  0.028922  0.006394
+2      2  1  1        1   total  0.022467  0.004039
+3      2  2  1        1   total  0.024923  0.009632