The two sides of the equation are equal up to a factor of $3$. That is, the left side is $\, q + q^4 + 2q^7 + \dots, \,$ the generating function of OEIS sequence A033687, and the right side is the generating function of OEIS sequence A005882 which is $3$ times that. It is also the right side of equation $(63)$ on page $111$ of Conway and Sloane "Sphere Packings, Lattices and Groups". On page $103$ is equation $(11)$ with the definition $\, \psi_k(z) = e^{\pi i/ z^2} \, \theta_3(\pi z/k|z) = \sum_{m=-\infty}^\infty q^{(m+1/k)^2}. \,$
The left side is equal to
$$\, \frac{q}3(2\, \psi(q^6)\, f(q^6, q^{12}) + \phi(q^3)\, f(q^3, q^{15}))$$ where $\phi(), \psi()$ are Ramanujan theta functions and $f(, )$ is Ramanujan's general theta function. The right side can be written as
$$ \theta_2(0, q^3)\, q^{1/4} f(q^6, q^{12}) + \theta_3(0, q^3)\, qf(q^3, q^{15}). $$ This can be shown by splitting the infinite sum into two infinite sums depending on the parity of $\,a\,$ and then identifying each sum with one of the two products of a theta null function and a Ramanujan general theta function.
In more detail, the infinite sum exponent of $q$ function
$$ e(a,b) := 3 (a^2 + b^2 + a b + a + b) + 1 $$
regarded as a matrix $\, \{e(i,j)\} \,$ has rank $3$ while the matrices $\, \{e(2i,j-i)\} \,$ and $\, \{e(2i+1,j-i)\} \,$
have rank $2$. Note that the matrix $\, \{q^{e(i,j)}\} \,$ has infinite rank, while the two matrices $\, \{q^{e(2i,j-i)}\} \,$ and $\, \{q^{e(2i+1,j-i)}\} \,$ have rank $1$, and thus the two infinite sums of their entries factor into products of two Ramanujan theta functions each.