No, this does not work this way as pointed out by another fellow in a different post of yours: the first \$RC\$ section is loaded by the Sallen-Key filter input impedance so you can't neglect it. You could if you would buffer the first network before driving the second filter. Another option would be to start on the left by the Sallen-Key structure then feeding the \$RC\$ filter. Considering a 0-\$\Omega\$ output impedance for the op amp then you could cascade (multiply) the transfer functions.
I've shown here without writing a single line of algebra how you can determine this transfer function using fast analytical circuits technique described here. Considering well-separated real poles, the transfer function can be approximated as

and if you plot this expression, you have

Edit: as correctly underlined by LvW, the above arrangement does not reflect a Bessel transfer function but rather a series of 3 cascaded poles. The transfer function describing this network actually obeys the following expression: \$H(s)=\frac{1}{1+b_1s+b_2s^2+b_3s^3}\$ in which
\$b_1=R_1C_1+(R_1+R_2+R_3)C_2\$
\$b_2=C_2(C_1R_1(R_2+R_3)+C_3R_3(R_1+R_2))\$
\$b_3=C_1C_2C_3R_1R_2R_3\$
The transfer function of a 3rd-order Bessel low-pass filter normalized at a 1-rad/s characteristic angular frequency is given by \$H(s)=\frac{1}{1+s+s^2\frac{6}{15}+\frac{s^3}{15}}\$. This expression can be refactored in a more friendly format, as its roots comprise 1 real pole and 2 complex poles: \$H(s)\approx\frac{1}{(1+0.43s)(1+0.57s+0.155s^2)}\$
Assume we would like to tune this filter to 1 kHz, then you have to scale the formula according to:
\$H(s)=\frac{1}{1+s\frac{1}{k}+s^2\frac{0.4}{k^2}+s^3\frac{0.067}{k^3}}\$
which is approximately equal to:
\$H(s)\approx\frac{1}{(1+\frac{0.43}{k}s)(1+\frac{0.57}{k}s+\frac{0.155}{k^2}s^2)}\$
with \$k=2\pi 1000=6.283\times 10^3\$
I can plot these two expressions and they give identical dynamic responses in magnitude and phase:

Now, you have to determine the component values so that the coefficients determined by the component combinations leads to:
\$b_1=159\;ms\$, \$b_2=1.013\times 10^4\;µs^2\$ and \$b_3=2.688\times 10^{-13}s^3\$
This is a system of 3 equations/3 unknown if you fix \$C_1\$, \$C_2\$ and \$C_3\$ for instance. If you use the tool here, then you have the following component values for a 1-kHz characteristic frequency:
\$R_1=9.1\;k\Omega\$ \$R_2=91\;k\Omega\$ \$R_3=36\;k\Omega\$ \$C_1=6.8\;nF\$ \$C_2=680\;pF\$ \$C_3=1.8\;nF\$
If now plot the expression that I determined here and featuring the above component values versus the response of the 3rd-order Bessel transfer function we have derived, the plots are identical:

H(w) =1/ (square root of 2) That sounds a bit confusing, we usually refer to the cutoff point as the -3 dB point.
That is the same though, -3 dB is half the power.
Let me explain: take your \$H(\omega) = \frac{1}{\sqrt2}\$
That means that at that \$\omega\$ the voltage is divided by \$\sqrt2\$, if this voltage is applied across a (load) resistor at the output of the filter then the current through that same load resistor will also be divided by \$\sqrt2\$.
What does that mean for the Power?
I means that the power is halved.
On a dB Power scale that means -3 dB
Using that "half of the power" as a reference point is useful because if we divide a wideband signal into a low frequency part and a high frequency part then we can do that using a low pass filter and a highpass filter with the same cutoff frequency. At that cutoff frequency half of the power ends up at the output of the lowpass filter and the other half ends up at the output of the highpass filter.
Best Answer
If you replace all components values in the circuit by 1 (\$1\;\Omega\$ and \$1\;F\$) and use the equation I derived here, then the below graphs show the dynamic response and I can extract the cutoff frequency to 53 mHz or 0.335 rad/s as you correctly found.