We present a new tensor network algorithm for calculating the partition function of interacting quantum field theories in 2 dimensions. It is based on the Tensor Renormalization Group (TRG) protocol, adapted to operate entirely at the level of fields. This strategy was applied in Ref. to the much simpler case of a free boson, obtaining an excellent performance. Here we include an arbitrary self-interaction and treat it in the context of perturbation theory. A real space analogue of the Wilsonian effective action and its expansion in Feynman graphs is proposed. Using a λφ4 theory for benchmark, we evaluate the order λ correction to the free energy. The results show a fast convergence with the bond dimension, implying that our algorithm captures well the effect of interaction on entanglement.