On P&S's QFT book, chapter 9.5, the book discussed how to derive two point correlation function for dirac field using generating functional.
Start with $$ Z[\bar{\eta}, \eta]=\int \mathcal{D} \bar{\psi} \mathcal{D} \psi \exp \left[i \int d^4 x[\bar{\psi}(i \not \partial-m) \psi+\bar{\eta} \psi+\bar{\psi} \eta]\right] \tag{9.73} $$ we can shift the field as: $$\begin{aligned}\psi &\leadsto \psi-i \int d^4 y S_F(x-y) \eta(y)\\ \bar{\psi} &\leadsto \bar{\psi}-i \int d^4 y \bar{\eta}(y) S_F(x-y) \end{aligned}$$
I met following term in the derivation: $$\int d^4 x d^4 y \bar{\eta}(y)S_F(x-y)(i\not\partial-m)\psi^{\prime}(x) $$
So how to handle this term? I know that $(i\not\partial-m)S_F(x-y)=i\delta^4(x-y)$, but this time the differential operator appear on the right side of $S_F(x-y)$.
In my attempt of derivation:
(1) Should I write them in components form? (i.e. $S_F(x-y)_{ab}$)
(2) If not, in the Integration by part, how do we handle $i\not\partial-m$? Since there is gamma matrices.
Another point is that I fail to shift the field by directly using $\psi^{\prime}\equiv \psi+(i\not\partial-m)^{-1}\eta$, I summarized this question in this post.
--
EDIT: $\stackrel{\star}{=}$ expansion in J.G.'s excellent answer! $$\begin{aligned}&\int d^4y \overline{\eta(y)}S_F(x-y)(i\not\partial_x-m)\psi(x) \\ &=\text{Surface}-\int d^4y \overline{\eta(y)}(i\not\partial_x+m)S_F(x-y)\psi(x)\\ &=\int d^4y \overline{\eta(y)}(i\not\partial_y-m)S_F(x-y)\psi(x) \\ &=\int d^4y \overline{\eta(y)}(i\not\partial_y-m)S_F(y-x)\psi(x)\end{aligned}$$