Frequency domain phenomenological models have been actively utilised for the description of gravitational waves from binary black hole systems (BBH). Finale state quantities as the total energy radiated and the final spin are implicitly used to characterise the merger-ringdown part of the phenomenological waveforms by means of two parameter fits on the asymmetric mass-ratio and some effective spin term which accounts for the contribution of the two individual spins. However, subdominant effects as the 'unequal spin difference' terms might play an important role to disentangle the well known spin degeneracy present in our effective spin when running parameter estimation codes. In this work we describe a bottom-up hierarchical methodology for building up 3 dimensional fits with an application to final spin and energy radiated including the extreme mass ratio limit and we investigate its effect in the near future phenomenological waveform descriptions.