We use the method of double pole QCD sum rule which is basically a fit with two exponentials of the correlation function, where we can extract the masses and decay constants of mesons as a function of the Borel mass. We apply this method to study the mesons: $rho(1S,2S)$, $psi(1S,2S)$, $Upsilon(1S,2S)$ and $psi_t(1S,2S)$. We also present predictions for the toponiuns masses $psi_t(1S,2S)$ of m(1S)=357 GeV and m(2S)=374 GeV.