[step:Realize the exceptional diagrams inside the standard $E_8$ root system]
Let $(e_1,\dots,e_8)$ be the standard orthonormal basis of $\mathbb{R}^8$, and define
\begin{align*}
\Phi(E_8)
&=
\{\pm e_i \pm e_j : 1 \le i<j \le 8\}\\
&\qquad \cup
\left\{
\frac12\sum_{i=1}^8 \varepsilon_i e_i :
\varepsilon_i \in \{\pm 1\},\ \prod_{i=1}^8 \varepsilon_i=1
\right\}.
\end{align*}
This is the standard finite reduced crystallographic root system of type $E_8$. Define roots
\begin{align*}
\alpha_1&=\frac12(e_1-e_2-e_3-e_4-e_5-e_6-e_7+e_8),\\
\alpha_2&=e_1+e_2,\\
\alpha_3&=e_2-e_1,\\
\alpha_4&=e_3-e_2,\\
\alpha_5&=e_4-e_3,\\
\alpha_6&=e_5-e_4,\\
\alpha_7&=e_6-e_5,\\
\alpha_8&=e_7-e_6.
\end{align*}
Each $\alpha_i$ belongs to $\Phi(E_8)$ and has squared length $2$. Direct computation gives the nonzero off-diagonal inner products
\begin{align*}
(\alpha_1,\alpha_3)&=-1,\\
(\alpha_2,\alpha_4)&=-1,\\
(\alpha_3,\alpha_4)&=-1,\\
(\alpha_4,\alpha_5)&=-1,\\
(\alpha_5,\alpha_6)&=-1,\\
(\alpha_6,\alpha_7)&=-1,\\
(\alpha_7,\alpha_8)&=-1,
\end{align*}
and all other off-diagonal inner products are $0$. Therefore the Dynkin diagram of $\alpha_1,\dots,\alpha_8$ is $E_8$.
The subsets $\{\alpha_1,\dots,\alpha_7\}$ and $\{\alpha_1,\dots,\alpha_6\}$ have the induced diagrams $E_7$ and $E_6$, respectively. The root subsystems generated by these subsets inside $\Phi(E_8)$ are finite, reduced, and crystallographic because they are closed subsystems of the finite reduced crystallographic root system $\Phi(E_8)$. Hence $E_6$, $E_7$, and $E_8$ are all realized by finite reduced crystallographic root systems.
Together with the coordinate realizations of $A_\ell$ and $D_\ell$, this proves the converse and completes the classification.
[/step]