Functional methods provide an alternate, fast and automatizable theoretical framework for model matching using an effective action formula. When evaluated for a BSM, this formula provides the Wilson coefficients of the effective operators as a function of model parameters. In this seminar, I shall discuss the integration-out of heavy fields from BSM theories and the computation of Wilson coefficients automatically. We examine the compatibility of multiple BSMs with the currently available data using their SMEFT Wilson coefficients. I shall also discuss the Renormalisation Group running of the SMEFT bosonic sector at 1-loop level up to order 1/Λ^4, where Λ is the SMEFT cut-off scale. Throughout my presentation, I will highlight several key points, including: (a) notable deviations in certain ADM elements from the expectation based on naive dimensional analysis and the largest mixing induced across different classes, (b) the mixing of tree-level generated operators into loop-level generated operators, which differs from what is known at dimension six, and (c) the contributions from LNVs to ADM elements, with comments on non-trivial vanishing entries.