I know there are many questions about this topic and also various answers, but it's never stated explicitly, why there is a certain sign before the mass term in the Dirac Lagrangian. I'm also confused, since the text I'm following, it is stated without proof that
The relative sign between the two Lorentz-Scalars can be deduced, by implying that the equation of motion must satisfy the Klein-Gordon-Equation for every components.
But when I tried to get to this result by keeping the $\pm$ explicitly through all calculations, I saw that it drops out anyways
$$\mathcal{L} = i\bar{\Psi}(\partial_\mu \gamma^\mu \pm m)\Psi$$
Euler-Lagrange equation yields
$$\frac{\partial\mathcal{L}}{\partial\bar{\Psi}} - \partial_\mu\frac{\partial\mathcal{L}}{\partial(\partial\bar{\Psi})} = (i\partial_\mu\gamma^\mu \pm m)\Psi = 0. $$
Multiplying the hermitian conjugate of the operator from the left
$$ \begin{align} 0 &= (-i\partial_\nu\gamma^\nu \pm m)(i\partial_\mu\gamma^\mu \pm m)\Psi \\ &=(\partial_\mu\gamma^\mu\partial_\nu\gamma^\nu \pm im\partial_\mu\gamma^\mu\mp im\partial_\nu\gamma^\nu+m^2)\Psi \\ &= (\frac{1}{2} \left[\partial_\mu\gamma^\mu\partial_\nu\gamma^\nu + \partial_\mu\gamma^\mu\partial_\nu\gamma^\nu\right] + m^2)\Psi \\ &= (\frac{1}{2}\{\gamma^\mu,\gamma^\nu\}\partial_\mu\partial_\nu + m^2)\Psi \\ &= (\partial_\mu\partial^\mu + m^2)\Psi = 0, \end{align}$$
where the identity matrix has been suppressed. So this seems indeed to be the Klein-Gordon-Equation for each component of the spinor, but the $\pm$ drops out in the beginning.
Hence, my question does it matter which sign we choose (at the level of the Lagrangian), and is there a deeper reason of why to choose one or the other?