I have something wrong in my understanding of the theory behind confidence intervals, as an example I have an exercise that states:
A sample of size $n = 100$ produced the sample mean of $\hat{X} = 16$. Assuming the population standard deviation is $\sigma = 3$ , compute a 95% confidence interval for the population mean $\mu$.
Now my reasoning is the following:
I know from the central limit theorem that $ \frac{\hat{X} - \mu}{\sigma / \sqrt{n}} $ is going to approximately normally distributed. I am looking for an $a$ s.t. $P(-a < \mu < a) = 0.95$ (since I know the distribution is symmetrical) but
$$P(-a < \mu < a) = P(-a < -\mu < a) = P(\frac{-a + \hat{X}}{\sigma / \sqrt{n}} < -\mu < \frac{+a + \hat{X}}{\sigma / \sqrt{n}})$$
So now I just need to look up in a $z$-table the value of $\frac{+a + \hat{X}}{\sigma / \sqrt{n}}$ and divide it by two since the distribution is symmetric.
I realized while writing this that probably saying $P(-a < \mu < a) = 0.95$ makes no sense since the mean of the population is not a random quantity but a fixed constant. I will leave my previous reasoning so an answerer can better understand where I am getting confused.
Could someone explain to me the theoretical passages that solve this exercise and clear up my confusion?