  This paper is devoted to the numerical analysis of a piecewise constant
discontinuous Galerkin method for time fractional subdiffusion problems. The
regularity of weak solution is firstly established by using variational
approach and Mittag-Leffler function. Then several optimal error estimates are
derived with low regularity data. Finally, numerical experiments are conducted
to verify the theoretical results.
