diff --git a/src/tallies/tally_header.F90 b/src/tallies/tally_header.F90
index 474e77d1c..9604d2382 100644
--- a/src/tallies/tally_header.F90
+++ b/src/tallies/tally_header.F90
@@ -732,8 +732,10 @@ contains
     integer :: MT
     character(C_CHAR), pointer :: string(:)
     character(len=:, kind=C_CHAR), allocatable :: score_
+    logical :: depletion_rx
 
     err = E_UNASSIGNED
+    depletion_rx = .false.
     if (index >= 1 .and. index <= size(tallies)) then
       associate (t => tallies(index) % obj)
         if (allocated(t % score_bins)) deallocate(t % score_bins)
@@ -757,10 +759,13 @@ contains
             t % score_bins(i) = SCORE_NU_SCATTER
           case ('(n,2n)')
             t % score_bins(i) = N_2N
+            depletion_rx = .true.
           case ('(n,3n)')
             t % score_bins(i) = N_3N
+            depletion_rx = .true.
           case ('(n,4n)')
             t % score_bins(i) = N_4N
+            depletion_rx = .true.
           case ('absorption')
             t % score_bins(i) = SCORE_ABSORPTION
           case ('fission', '18')
@@ -829,8 +834,10 @@ contains
             t % score_bins(i) = N_NC
           case ('(n,gamma)')
             t % score_bins(i) = N_GAMMA
+            depletion_rx = .true.
           case ('(n,p)')
             t % score_bins(i) = N_P
+            depletion_rx = .true.
           case ('(n,d)')
             t % score_bins(i) = N_D
           case ('(n,t)')
@@ -839,6 +846,7 @@ contains
             t % score_bins(i) = N_3HE
           case ('(n,a)')
             t % score_bins(i) = N_A
+            depletion_rx = .true.
           case ('(n,2a)')
             t % score_bins(i) = N_2A
           case ('(n,3a)')
@@ -879,6 +887,7 @@ contains
         end do
 
         err = 0
+        t % depletion_rx = depletion_rx
       end associate
     else
       err = E_OUT_OF_BOUNDS