Let $C\left[0,\infty\right)$, the space of all continuous, real-valued functions on $\left[0,\infty\right)$ with metric
$\rho\left(\omega_{1},\omega_{2}\right)=\sum_{n=1}^{\infty}\frac{1}{2^{n}}\max_{0\le t\le n}\left(\min\left\{ 1,\left|\omega_{1}\left(t\right)-\omega_{2}\left(t\right)\right|\right\} \right)$
I have proved that Under $\rho$, $C\left[0,\infty\right)$ is a complete, separable metric space
Let $\mathscr{C}$ be the collection of finite-dimensional cylinder sets i.e.,
$C=\left\{ \omega\in C\left[0,\infty\right):\left(\omega\left(t_{1}\right),\cdots,\omega\left(t_{n}\right)\right)\in A\right\} ,n\geq1,A\in\mathscr{B}\left(\mathbb{R}^{n}\right) $
where, for all $i=1,2,\cdots,n, t_{i}\in\left[0,\infty\right)$. Denote by $\mathscr{G}$ the smallest $\sigma$-field containing $\mathscr{C}$. Now we need to show that $\mathscr{B}\left(C\left[0,\infty\right)\right)=\mathscr{G}$.
I have proved that $\mathscr{G}\subset\mathscr{B}\left(C\left[0,\infty\right)\right)$. To prove the converse statement. By Karatzas&Shreve Brownian Motion and Stochastic Calculus, I just need to show that $B\left(\omega_{0},\epsilon\right)=\left\{ \omega\in C\left[0,\infty\right):\rho\left(\omega_{0},\omega\right)<\epsilon\right\} \in \mathscr{G}$. Let $Q$ be the set of rationals in $\left[0,\infty\right)$. We have
$ B\left(\omega_{0},\epsilon\right) =\left\{ \omega:\sum_{n=1}^{\infty}\frac{1}{2^{n}}\sup_{0\le t\le n,t\in Q}\left(\min\left\{ 1,\left|\omega\left(t\right)-\omega_{0}\left(t\right)\right|\right\} \right)<\epsilon\right\}$
But how can I show the set on the right is in $\mathscr{G}$? Can anyone give a hint? Many thanks in advance.

Statements like this are typically proved by inspecting all the "building blocks" of the object being defined, often from the inside out, and checking that they all preserve measurability.
For each $t$, the map from $C[0,\infty)$ to $\mathbb{R}$ defined by $\omega \mapsto \omega(t)$ is measurable with respect to $\mathcal{C}$.
The function from $\mathbb{R}$ to $\mathbb{R}$ defined by $x \mapsto \min\{1, |x-\omega_0(t)|\}$ is Borel, so its composition with measurable functions is measurable. Thus $\omega \mapsto \min\{1, |\omega(t) - \omega_0(t)|$ is measurable with respect to $\mathcal{C}$.
The pointwise supremum of countably many measurable functions is measurable. So $\omega \mapsto \sup_{0 \le t \le n, t \in \mathbb{Q}} \min\{1, |\omega(t) - \omega_0(t)|$ is measurable with respect to $\mathcal{C}$.
Scalar multiples and infinite sums of measurable functions are measurable. So the map $\omega \mapsto \sum_{n=1}^{\infty}\frac{1}{2^{n}}\sup_{0\le t\le n,t\in Q}\min\left\{ 1,\left|\omega\left(t\right)-\omega_{0}\left(t\right)\right|\right\} $ is measurable.
Your set $B(\omega_0, \epsilon)$ is the preimage of the Borel set $[0,\epsilon)$ under that measurable function.