How to write this tensor?

How can I write this tensor in Mathematica?

$ $ \mathcal{P}_{ijkl}(N) = \Big( \delta_{ik} – N_iN_k \Big) \Big( \delta_{jl} – N_jN_l \Big) – \dfrac{1}{2} \Big( \delta_{ij} – N_iN_j \Big)\Big( \delta_{kl} – N_kN_l \Big)$ $

I have tried with TensorProduct, but I realized that it is not so simple because of the position of the indices.