\text{g}({\bf x})\big)$ is the Laplace functional of the received power. We consider the following general path-loss model $ \text{g}({\bf x}) = 1/(\alpha_1 + \alpha_2||{\bf x}||^\delta),$ where $\delta$ is the pathloss exponent, and $\alpha_1$ and $\alpha_2$ are control parameters. When $h$ follows Nakagami-$m$ fading, with the shaping and spread parameters of $m\in \mathbb Z^i$ and $\Omega>0$ respectively, the probability density function (PDF) of the power fading coefficient is given by: $$\text{p}_h(q)= \frac{1}{{\Gamma(m)}}(\frac{m}{\Omega})^m q^{m-1}\exp\big({-\frac{mq}{\Omega}}\big),\label{nm}$$ where $\Gamma$ is the Gamma function. Then, using Laplace table, $\mathcal L_{h}\big(sQ_j  P_k \text{g}({\bf x})\big)$ is derived as: $$\label{laph}L_h(s  Q_j P_k \text{g}({\bf x}))={\big(1+{\Omega}s P_k \text{g}({\bf x})/m\big)^{-m}}.$$ By inserting in and considering the fact that the received interferences from different devices are independent, we have: $$\mathcal L_{I_{\Psi}}(s)=  \prod\nolimits_{j,k}\mathbb E_{{\bf x},{\bf y}}\bigg[  \prod\nolimits_{{\bf y}\in \Theta_{k}}\big(  \prod\nolimits_{{\bf x}\in \theta_{\bf y}^{(j)}}u({\bf x},{\bf y})\big)\bigg],\nonumber$$ where $k\in\mathcal K$, $j\in\{1,2\}$, the set of parent points of type-$k$ is denoted by $\Theta_{k}$, and transmitting nodes which are daughter points of $y$ as $\theta_{\bf y}^{(j)}$. Also, $\mathbb E_x$ represents expectation over $x$, and $$u({\bf x},{\bf y})={\big(1 {+}{\Omega}s Q_j P_k \text{g}({\bf x} {-}{\bf y})/{m}\big)^{-m}}.$$ The received interference over the packet of interest can be decomposed into two parts: i) interference from transmitters belonging to the cluster of transmitter, 