In the paper we consider the following quasilinear Schrödinger–Poisson system in the whole space \({\mathbb {R}}^{3}\) \(\begin{aligned} {\left\{ \begin{array}{ll} - \varepsilon ^2 \Delta u + ({V + \phi }) u = u \left| {u}\right| ^{p - 1} \\ - \Delta \phi - \beta \Delta _4 \phi = u^2, \end{array}\right. } \end{aligned}\) where \(1< p < 5, \beta > 0,V:{\mathbb {R}}^{3}\rightarrow ]0, \infty [\) , and look for solutions \(u,\phi :{\mathbb {R}}^{3}\rightarrow {\mathbb {R}}\) in the semiclassical regime, namely when \(\varepsilon \rightarrow 0.\) By means of the Lyapunov–Schmidt method we estimate the number of solutions by the cup-length of the critical manifold of the external potential V.