I am doing a difference-in-difference analysis of an event that affected several states in the US. I am interested to understand the effects of this event on state-level unemployment rates. I have state level data on demographics for several years before and after this event. My question is what covariates should I include if I am estimating equation of the following form:
$$ \mbox{unempRate}_{jt} = \gamma_j + \alpha_t + \beta D_{jt} + \delta X_{jt} + \epsilon_{jt} $$
where $\gamma_j$ are state-fixed effects; $\alpha_t$ are time fixed effects; $D_{jt}$ are dummies - 1 if state $j$ is affected at time $t$, and 0 otherwise; and $X_{jt}$ are time-varying state level covariates. Coefficient of interest is $\beta$.
So far, I have included population and average household income at the state level.
What more covariates can I include? How does one decide which covariates to include in a setting like this? What is the guiding philosophy?