Bare bones bi-photon density matrix simulation