Print total charges (within cut-off and all)
Add a comparison with all electron calculation