When you have an action with bosonic $X$ and fermionic $\psi$ (Majorana) fields and perform a SUSY transformation $\epsilon$ (the constant, infinitesimal parameter of transformation, a real, anticommuting spinor) can you do normal partial integration on the Dirac operator just like you're used doing with a derivative? For example: $$ \begin{align} S &= \int \mathrm d^{2}\sigma\; \bar{\epsilon}[\left(\gamma^{\alpha}\partial_{\alpha}\gamma^{\beta}\partial_{\beta}\right)X(\sigma)]\psi(\sigma)\\ &= \int \mathrm d^{2}\sigma\; \bar{\epsilon}[{\not}\partial{\not}\partial X(\sigma)]\psi(\sigma)\\ &= -\int \mathrm d^{2}\sigma\; \bar{\epsilon}[{\not}\partial X(\sigma)][{\not}\partial\psi(\sigma)] + `boundary` \end{align} $$
This would make calculations much easier because I'm lost with all matrix multiplications and properties of the Gamma matrices...I wanted to check because I couldn't find any information about this.