Given the Hilbert action: $$ S_{H}=\int \sqrt{|g|}R d^{n}x $$ and the metric written in terms of Minkowski and perturbed metric: $$ g_{\mu \nu}=\eta_{\mu \nu}+h_{\mu \nu}. $$
I am able to derived the suitable Lagrangian for Linearised theory is: $$ \mathcal{L}=\frac{1}{2}[ (\partial_{\mu}h^{\mu \nu})(\partial_{\nu}h)-(\partial_{\mu}h^{\rho \sigma})(\partial_{\rho}h_{\sigma}^{\mu})+\frac{1}{2}\eta^{\mu \nu}(\partial_{\mu}h^{\rho \sigma})(\partial_{\nu}h_{\rho \sigma})-\frac{1}{2}\eta^{\mu \nu} (\partial_{\mu}h)(\partial_{\nu}h) ] $$ which using the Euler-Lagrange equation should give rise to the linearised Einstein Tensor, nonetheless, I am not able to produce the result. So I have the following questions:
The Lagrangian I obtained is it correct?
If it's correct, can anyone show the remaining steps lead to linearised Einstein Tensor.