Three and four point functions of stress energy tensors in D = 3 for the analysis of cosmological non-gaussianities