If it is a photon, then you know that $x^{a}x_{a}$ is always zero. You have measured the spatial direction of the photon with your mirror apparatus, so you know the values of the $x^{i}(0)$ at some time${}^{1}$, which we will call zero. Furthermore, we can infer the value of $x^{0}$ at this time from the fact that $x^{a}x_{a} = 0$.
Then, all we need is to know the metric tensor, and we can simply solve the set of ODEs given by,
$$\frac{\mathrm{d}^2 x^a (\lambda)}{\mathrm{d}\lambda^2} + \Gamma^a_{bc} \, \frac{\mathrm{d} x^b(\lambda)}{\mathrm{d}\lambda} \frac{\mathrm{d}x^c (\lambda)}{\mathrm{d}\lambda} = 0$$
for all subsequent times${}^{2}$. Typically, this will have to be done numerically, but there are some known exact solutions for simplified cases. $\lambda$ is the null affine parameter, and the $\Gamma$'s are the Christoffel symbols, which are defined from the metric tensor. For the familiar example of a spherically symmeteric gravitating body AND for weak gravity in nonrelativistic motion, the only relevant component of the metric tensor is $g_{tt} = -\left(1-\frac{2M}{r}\right)$, and the only relevant Christoffel symbol is $\Gamma_{tt}{}^{r} = \frac{M}{r^{2}}$, and our theory becomes Newtonian gravity.
${}^{1}$Here, I take letters from the beginning of the alphabet to be spacetime indices, while letters from the middle of the alphabet are sums over only the spatial components.
${}^{2}$ note that "time" (given by $\lambda$) here doesn't really mean time in the proper sense. Null geodesics don't experience any clock motion. Null geodesics, however, DO trace out "null affine parameter", which you can interpret as how many perpendicular null geodesics they cross as the photon moves.