First some definitions:
For any $n > 0$, let $\mathcal{B}_n$ be the collection of covering of boxes of dimension $2 \times 2 \times n$.
For any $B_1 \in \mathcal{B}_p$, $B_2 \in \mathcal{B}_q$, we can "join" them together to obtain a $B_1B_2 \in \mathcal{B}_{p+q}$.
For any $B \in \mathcal{B}_n$, we call it irreducible if $B \ne B_1B_2$ for any $B_1 \in \mathcal{B}_p, B_2 \in \mathcal{B}_q$ where $0 < p, q \le n$.
Let $\mathcal{C}_n \subset \mathcal{B}_n$ be the collection of irreducible coverings of boxes of dimension $2 \times 2 \times n$.
It is easy to see there are three types of irreducible coverings.
there are two ways to place two dominoes to complete a single layer.
(counting from left to right, the $1^{st}$ and $2^{nd}$ figures in picture below)
there are four ways to place one domino into first layer.
After you do this, you place two domino vertically which stick out into second
layer. You can either use one more domino to complete the second layer
or place two more vertically which stick out into the third layer.
This process repeat, you obtain four irreducible coverings in each $\mathcal{B}_n$ for $n \ge 2$.
(the $3^{rd}$, $4^{th}$ and $5^{th}$ figures in picture)
- there are one extra way to place four domino vertically to complete
a double layer.
(the rightmost figure in picture)
$\hspace1in$ 
Combine this, we see
$$|\mathcal{C}_n| = \begin{cases}
2, & n = 1,\\
5, & n = 2,\\
4, & n > 2.
\end{cases}$$
Translate this to generating function, the OGF for irreducible covering is
$$\text{OGF}_{\mathcal{C}}(z) \stackrel{def}{=} \sum_{n=1}^\infty |\mathcal{C}_n| z^n = 2z + \frac{4z^2}{1-z} + z^2$$
Since we can obtain any covering in $\mathcal{B}_?$ by joining coverings in $\mathcal{C}_?$ in a unique manner, the OGF for all coverings is
$$\begin{align}
\text{OGF}_{\mathcal{B}}(z)
\stackrel{def}{=} 1 + \sum_{n=1}^\infty |\mathcal{B}_n| z^n
&= 1 + \sum_{k=1}^\infty \text{OGF}_{\mathcal{C}}(z)^k
= \frac{1}{1 - \text{OGF}_{\mathcal{C}}(z)}\\
&= \frac{1}{1 - \left(2z + 4\frac{z^2}{1-z} + z^2\right)} = \frac{1-z}{1-3z-3z^2+z^3}\\
\end{align}$$
Let $\alpha = 2+\sqrt{3}$ and $\beta = 2-\sqrt{3}$, we can simplify RHS as
$$\frac{1-z}{(1+z)(1-4z+z^2)} = \frac{1-z}{(1+z)(1-\alpha z)(1-\beta z)}
= \frac16\left(\frac{\alpha}{1-\alpha z} + \frac{\beta}{1-\beta z}\right) + \frac{1}{3(1+z)}$$
Expanding as a power series of $z$, we get
$$1 + \sum_{n=1}^\infty |\mathcal{B}_n| z^n = \text{OGF}_{\mathcal{B}}(z)
= 1 + \sum_{n=1}^\infty \left[\frac16(\alpha^{n+1}+\beta^{n+1}) + \frac13 (-1)^n\right] z^n$$
This leads to
$$|\mathcal{B}_n| = \frac16 (\alpha^{n+1}+\beta^{n+1}) + \frac13(-1)^n\quad\text{ for } n \ge 1$$
Notice $|\beta | < 1$ implies $\left|\frac16\beta^{n+1} + \frac13(-1)^n\right| < \frac12$, above formula reduces to
$$|\mathcal{B}_n| = \verb/NearestInteger/\left( \frac16 (2+\sqrt{3})^{n+1} \right)$$
Compute first few numbers of $|\mathcal{B}_n|$ explicitly, we have
$$|\mathcal{B}_1|, |\mathcal{B}|_2, \ldots
= 2, 9, 32, 121, 450, 1681, 6272, 23409, \ldots$$
An OEIS search will lead you to OEIS A006253.
Look at references there for more information.