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: ρ(1S, 2S), ψ(1S, 2S), Υ(1S, 2S) and ψt(1S, 2S). We also present predictions for the toponiuns masses ψt(1S, 2S) of m(1S) = 357 GeV and m(2S) = 374 GeV.