In Subsection 33.2.2. of Schwartz's Quantum Field Theory and the Standard Model, he starts to derive the Euler-Heisenberg effective Lagrangian by "replacing" the field which is being integrated out with a corresponding classical expectation value.
Specifically, suppose that $|A\rangle$ is a state corresponding to some fixed background configuration of the potential $A^\mu$. Schwartz starts with the standard QED Lagrangian $$\mathcal{L} = -\frac{1}{4}F_{\mu\nu}^2 + \overline{\psi}(i\gamma^\mu\partial_\mu - m)\psi - eA_\mu\overline{\psi}\gamma^\mu\psi,\tag{33.31}$$ and then simply declares that the result of integrating out the fields $(\psi,\overline{\psi})$ is given by $$\mathcal{L}_{\mathrm{eff}} = -\frac{1}{4}F_{\mu\nu}^2 - eA_\mu J_A^\mu,\tag{33.32}$$ where $$J_A^\mu = \langle A|\overline{\psi}\gamma^\mu\psi|A\rangle\tag{33.33}.$$ What is the justification for this step? What happened to quadratic term $\overline{\psi}(i\gamma^\mu\partial_\mu - m)\psi$?
Note that Schwartz has, throughout his book, a rather idiosyncratic way of integrating out fields by solving their equations of motion (see for example the discussion here). This does not seem to be what's happening here however.
Also, just to be clear, I am not asking about the usual process of integrating out fields by performing the marginal path integrals. Nor am I really interested in the derivation of the Euler-Heisenberg Lagrangian itself. I am primarily just curious about whether the "shortcut" performed by Schwartz here to integrate out the $\psi$ field has any real justification. Something akin to Adam's answer (point 1) in the linked question would be a perfect example of the sort of thing I'm looking for.