New covariant Feynman rules for effective field theories

We provide a new and completely general formalism to compute the effective field theory matching contributions from integrating out massive fields in a manifestly gauge covariant way, at any desired loop order. The formalism is based on old ideas such as the background field method and the heat kernel, however we add some crucial new ingredients that greatly improve the simplicity and general applicability of the approach. We formulate our method in terms of Feynman rules, the resulting effective action is expressed in terms of local heat kernel coefficients. We also provide as supplementary material a mathematica code that facilitates the computation of these coefficients.