In this article, we investigate the Kirchhoff-Schrödinger-Poisson type systems on the Heisenberg group of the following form:

$$\{\begin{array}{ll}-(a+b{\int }_{\mathrm{\Omega }}|{\mathrm{\nabla }}_{H}u{|}^{p}d\xi ){\mathrm{\Delta }}_{H,p}u-\mu \varphi |u{|}^{p-2}u=\lambda |u{|}^{q-2}u+|u{|}^{Q^{\ast }-2}u& \text{in}\mathrm{\Omega },\\ -{\mathrm{\Delta }}_H\varphi =|u{|}^p& \text{in}\mathrm{\Omega },\\ u=\varphi =0& \text{on}\mathrm{\partial }\mathrm{\Omega },\end{array}$$

where $a,b$ are positive real numbers, $\mathrm{\Omega }\subset {\mathbb{H}}^N$ is a bounded region with smooth boundary, $1<p<Q$, $Q=2N+2$ is the homogeneous dimension of the Heisenberg group ${\mathbb{H}}^N$, $Q^{\ast }=\frac{pQ}{Q-p}$, $q\in (2p,Q^{\ast })$ and ${\mathrm{\Delta }}_{H,p}u=\text{div}(|{\mathrm{\nabla }}_{H}u{|}^{p-2}{\mathrm{\nabla }}_{H}u)$ is the $p$-horizontal Laplacian. Under some appropriate conditions for the parameters $\mu$ and $\lambda$, we establish existence and multiplicity results for the system above. To some extent, we generalize the results of An and Liu (Israel J. Math., 2020) and Liu et al. (Adv. Nonlinear Anal., 2022).