We evaluate the entanglement entropy and entropic function of massive two dimensional QED (Schwinger model) at finite temperature, density, and $\theta$-angle. In the strong coupling regime, the entropic function is dominated by the boson mass for large spatial intervals, and reduces to the CFT result for small spatial intervals. We also discuss the entanglement spectrum at finite temperature and a finite $\theta$-angle.
Comment: v2: matches version published in PRD; 10+3 pages, 7+1 figures