It is not. You can 100% fix the gravitational constant using experiments. A simple way to see this is to pick a classical free-fall experiment. Here's the version with only some unnamed constant $\kappa$ where
$$G_{\mu\nu} = \kappa T_{\mu\nu}$$
If you consider the earth as a Kerr metric source, then, neglecting angular momentum (which is fairly small), we'll end up with the Schwarzschild metric
$$ds^2 = -(1 - \frac{r_S}{r}) c^2 dt^2 + (1 - \frac{r_s}{r})^{-1} dr^2 + r^2 d\Omega^2$$
$r_s$ is the Schwarzschild radius, some parameter of the metric. To compute it, we perform the analysis of the Komar mass. If we go 100% without any definite constants from previous theories, let's just write the Komar mass as proportional to some constant.
$$M = -\alpha \int_S dA n_\mu \sigma_\nu \nabla^\mu \xi^\nu$$
with $\xi$ the normalized Killing timelike vector (we'll just pick $ c\partial_t$), $\sigma$ the vector normal to the Cauchy surface, and $n$ the vector normal to the surface of integration. The constant $\alpha$ is related to the constant $\kappa$ since the Komar integral is alternatively written in term of the stress-energy tensor
$$M = -\alpha \int_\Sigma R_{\mu\nu} n^\mu \xi^\nu = -\alpha \kappa \int_\Sigma (T_{\mu\nu} + \frac 12 T g_{\mu\nu}) n^\mu \xi^\nu dV$$
We will really want $M$ to be of the same dimension as the mass. Since $T$ has the dimension of an energy density, we want $\alpha = \beta / c^2 \kappa $ for some constant $\beta$.
Then, picking the surface of constant $t$ and $r$ as $S$, we get
$$M = -\alpha \int_S r^2 \sin \theta d\theta d\varphi \nabla^r\xi^t$$
$\nabla^r\xi^t$ has for only non-zero term $c g^{rr} \Gamma^t_{rt}$, which is equal to $c g^{rr} g^{tt} g_{tt,r}/2 = -c r_S/2r^2 $, so
$$M = -\alpha \int_S \sin \theta d\theta d\varphi c r_s = 8\pi c \alpha r_S$$
We may then write $r_S = M / 8\pi c \alpha = M \kappa / 8\pi \beta c$.
The geodesic equation for a purely radial motion will then be
$$c^{-2} \dot r^2 + (1 - \frac{r_S}{r}) = E^2$$
Expanding the $r^{-1}$ term around $R_\oplus$, the radius of the earth, we get
$$\frac{1}{r} = \frac{1}{R_\oplus} - \frac{r - R_\oplus}{R_\oplus^2} + \mathcal O(R_\oplus^{-3})$$
This will give you the equation, cutting off the smallest term,
$$c^{-2} \dot r^2 + \frac{r_S}{R_\oplus^2} r = E^2 - 1 + \frac{2r_S}{R_\oplus} $$
If you take the proper time derivative,
$$c^{-3} \ddot r \dot r + 2 \frac{r_S}{cR_\oplus^2} \dot r= 0$$
or in other words,
$$\ddot r = - 2 c^2 \frac{r_S}{ R_\oplus^2}$$
The classic equation of free fall, which you can check experimentally. If you wish for less classical results, of course, you'll need to solve the geodesic equation properly, which I don't think can be done analytically (numerical simulations will do there).
You'll notice that this equation has the correct dimension and sign. We can then input back the expression of $r_S$
$$\ddot r = - 2 c^3 \frac{M \kappa}{8\pi \beta R_\oplus^2}$$
I think some $c$'s might have gotten lost in the shuffle, but you can see that globally, we get back roughly what we were looking for : the free fall of an object depends on some constants which, if you were to sort out a bit all the various constants, would turn out to be $G$.