subroutine print_reference_notification(iStart, iEnd, title, legend)
implicit none
integer, intent(in) :: iStart, iEnd
character(*), intent(in) :: title
logical, optional :: legend
integer :: i
if (iProcIndex == root .and. .not. tSuppressSIOutput) then
! print out the given SIs sorted according to population
call sort(ilutRefAdi(0:NIfTot, iStart:iEnd), sign_gt, sign_lt)
write(stdout, *) title
if (present(legend)) write(stdout, "(4A25)") &
! TODO: Adapt legend for multiple runs
"Determinant (bitwise)", "Excitation level", "Coherence parameter", "Number of walkers"
do i = iStart, iEnd
call print_bit_rep(ilutRefAdi(:, i))
end do
end if
end subroutine print_reference_notification