I'm trying to calculate the interference of two non-entangled photons,
That's like saying you are trying to compute the wave function of a unicorn.
The individual particles have the wave function $\psi_i\in\mathcal{H}_i$ where $i\in\{1,2\}$.
That's not quite right. There is a space $\mathcal{H}_i$ of single particle states. And the space of states of a multiparticle space is a subset of the tensor product of the single particle states. And the tensor product is $\mathcal{H}_1\otimes\mathcal{H}_2$ and what subset it is depends on whether the single particles are fermions or bosons and whether the two particles are identical or different.
Now why did I object so early? Because a general state of the multiple system is not of the form $\psi_1\otimes \psi_2$ with $\psi_i\in\mathcal{H}_i$ even when they are distinct. A state like $\psi_1\otimes \psi_2$ is a special multiparticle state, one where there is a single particle state for each particle (namely $\psi_1$ and $\psi_2$). And in general this doesn't happen. In general these states are like the x axis and y axis in $\mathbb R^2$ they are possible, but not general. The general states are linear combinations of these basis states and hardly ever are these states.
The tensor product of their wave function is
$$
\Psi(1,2) = c_1\Psi_1 + c_2\Psi_2 = c_1(\psi_1\otimes 1_2) + c_2(1_1\otimes\psi_2)
$$
The state you just wrote out as $c_1(\psi_1\otimes 1_2) + c_2(1_1\otimes\psi_2)$ is not the tensor product of $\psi_1$ and $\psi_2$ the tensor product of $\psi_1$ and $\psi_2$ is $\psi_1\otimes\psi_2$ which is different. And it is, for instance, factorizable. Whereas the one you have is in general not.
Which is another problem, you said you wanted a non entangled state. But that means you want a factorizable state. And in general that isn't possible. Sure you could have bosons in a state that is identical i.e. if $\psi_1=\psi_2$ then $\psi_1\otimes\psi_2$ is possible but in general only states like $c\psi_1\otimes\psi_2+c\psi_2\otimes\psi_1$ Are possible which means you have to learn to distinguish between states on the left of the tensor symbol and states on the right rather than letting the numbers do the distinguishing for you. I mean you could let $a=\psi_1$ and let $b=\psi_2$ and then talk about a state like $a\otimes b +b\otimes a $ but guess what? That state is entangled. Entanglement isn't just common, like being not on the x axis or the y axis is common entanglement is almost required.
And when it isn't required the states are the same, so it is hard to have a state aimed at just one slit and then not be entangled.
Hence you are asking for the wave function of a unicorn. Except evolution could have produced a unicorn and scientists could try to someday make one. Whereas these symmetry rules come from the statistics required to observationally get the thermodynamics we observe. So there are really good reasons to think that it is literally impossible to have what you want.
So tour whole goal of a non entangled state with each particle aimed at just one slit is misguided and impossible. And you didn't write down the tensor product correctly, and you wrongly assumed a multiparticle state allows there to be single particle states instead of being made out of linear combinations of products of single particle states.
But what if it were? Then, if you compute $|\psi_1(x_1,y_1,z_1)\otimes\psi_2(x_2,y_2,z_2)|^2$ you get $|\psi_1(x_1,y_1,z_1)|^2|\psi_2(x_2,y_2,z_2)|^2$. And check my work questions are off topic which is why I spent so much time addressing more fundamental conceptual errors. And I can't answer your question about how to proceed towards showing interference for non entangled photons where each one is heading towards just one slit since that doesn't sound physically possible to have such a thing even before it went through the slits.