Regarding computing DOS for graphene see Appendix in Phys. Rev. B 87, 075414 – Published 7 February 2013. As Tue suggested also try increasing the number of k-points significantly. Did you also use tetrahedron method with significantly increased k-point density? It would also be helpful to see your script.