We consider the point process \begin{align*} \frac{1}{Z_{n}}\prod_{1 \leq j < k \leq n} |e^{iθ_{j}}-e^{-iθ_{k}}|^β\prod_{j=1}^{n} dθ_{j}, \qquad θ_{1},\ldots,θ_{n} \in (-π,π], \quad β> 0, \end{align*} where $Z_{n}$ is the normalization constant. The feature of this process is that the points $e^{iθ_{1}},\ldots,e^{iθ_{n}}$ interact with the mirror points reflected over the real line $e^{-iθ_{1}},\ldots,e^{-iθ_{n}}$. We study smooth linear statistics of the form $\sum_{j=1}^{n}g(θ_{j})$ as $n \to \infty$, where $g$ is $2π$-periodic. We prove that a wide range of asymptotic scenarios can occur: depending on $g$, the leading order fluctuations around the mean can (i) be of order $n$ and purely Bernoulli, (ii) be of order $1$ and purely Gaussian, (iii) be of order $1$ and purely Bernoulli, or (iv) be of order $1$ and of the form $BN_{1}+(1-B)N_{2}$, where $N_{1},N_{2}$ are two independent Gaussians and $B$ is a Bernoulli that is independent of $N_{1}$ and $N_{2}$. The above list is not exhaustive: the fluctuations can be of order $n$, of order $1$ or $o(1)$, and other random variables can also emerge in the limit. We also obtain large $n$ asymptotics for $Z_{n}$ (and some generalizations), up to and including the term of order $1$. Our proof is inspired by a method developed by McKay and Wormald [12] to estimate related $n$-fold integrals.