anchor.diff:    real(kind=fPrec), allocatable, save :: sumimpact(:) !(50)
     real(kind=fPrec), allocatable, save :: sqsumimpact(:) !(50)
   
  -  character(len=:), allocatable, save :: ename(:) !(mNameLen)(nblz)
     integer, allocatable, save :: nampl(:) !(nblz)
     real(kind=fPrec), allocatable, save :: sum_ax(:) !(nblz)
     real(kind=fPrec), allocatable, save :: sqsum_ax(:) !(nblz)
  
positive.code:         end if
  
          sampl(ie) = totals
          ename(ie) = bez(myix)
        end do
      end if
    end if
  
negative.code:   use hdf5_tracks2
  #endif
  
    integer i, ix
  
    ! Just call it here since samples are no longer supported
    call collimate_end_sample(1)
  