A textbook says that the geodesic equations (for both massive and massless) can be derived from the following action:
$$ S = -\frac{1}{2} \int d\tau \:\eta \: (\eta^{-2} \dot{x}^\mu \dot{x}^\nu g_{\mu\nu}(x) + m^2) $$ where $\eta$ is an auxiliary field. The signature convention is here $(+,-,-,-)$. The e.o.m. derived are $$ \eta^{-1}\dot{x}^\mu \dot{x}^\nu \partial_\sigma g_{\mu\nu}(x) - \frac{d}{d\tau}( 2 \eta^{-1}g_{\sigma \mu}(x) \dot{x}^\mu ) = 0, \quad \quad \eta^{-2} \dot{x}^\mu \dot{x}^\nu g_{\mu\nu}(x) = m^2.$$
In the massive case, by redefining $\tau$ to be proper time, $\eta^{-1} = m$ and the geodesic equations come out.
1) I'm confused about the massless case: It seems that $\eta$ is undetermined as we have $\eta^{-2} \:0 =0$.
The text then goes to say that one can think of $S$ as the action of 4 fields $x^{\mu}$ living on a 1D space with metric $\eta(\tau)^2$. I understand that $ds^2 = d\tau^2 \eta(\tau)^2$ and therefore $d\tau \eta$ is the invariant measure.
2) I don't really understand what would motivate the factor of $\eta^{-2}$ next to $\dot{x}^\mu \dot{x}^\nu g_{\mu\nu}(x)$ when writing a theory for these 4 fields though. I realise $$\eta^{-2}\dot{x}^\mu \dot{x}^\nu g_{\mu\nu}(x) = \frac{dx^\mu}{ds}\frac{dx^\nu}{ds} g_{\mu\nu}(x).$$ I guess I don't get why one would want to write derivatives w.r.t. $s$ instead of $\tau$ when writing a theory of 4 fields on 1D space with coordinate $\tau$.