add the part for calculating the density of state at all sites
the results seem OK