In this paper, we work with 16 different single scalar particle extensions of the Standard Model. For the first time, we present the exhaustive sets of dimension-6 effective operators and the associated Wilson coefficients as functions of model parameters after integrating out the heavy scalars up to 1-loop, including the light-heavy mixing, for each such scenario. Using the correspondence between the effective operators and the observables at electroweak scale, and employing Bayesian statistics, we compute the allowed ranges of new physics parameters that are further translated and depicted in 2-dimensional Wilson coefficient space in the light of the latest CMS and ATLAS data up to 137 fb −1 and 139 fb −1 , respectively. We also adjudge the status of those new physics extensions that offer similar sets of relevant effective operators. In addition, we provide a model-independent fit of 23 Standard Model effective field theory Wilson coefficients using electroweak precision observables, single and di-Higgs data as well as kinematic distributions of di-boson production.