Add user editable threshold for selecting states of the initial Hamiltonian

[mretegan]

At the moment the threshold is set to the machine epsilon for double precision, which is very conservative.