Modelling the angular correlation function and its full covariance in photometric galaxy surveys