||We present a lattice calculation of the masses and decay constants of D-(s)* and B-(s)* mesons using the gauge configurations produced by the European Twisted Mass Collaboration (ETMC) with N-f=2+1+1 dynamical quarks at three values of the lattice spacing a similar to(0.06-0.09) fm. Pion masses are simulated in the range M-pi similar or equal to(210-450) MeV, while the strange and charm sea-quark masses are close to their physical values. We compute the ratios of vector to pseudoscalar masses and decay constants for various values of the heavy-quark mass mh in the range 0.7m(c)(phys) less than or similar to m(h) less than or similar to 3m(c)(phys). In order to reach the physical b-quark mass, we exploit the Heavy Quark Effective Theory prediction that, in the static limit of infinite heavy-quark mass, the considered ratios are equal to one. At the physical point our results are: M-D*/M-D=1.0769(79), M-D*(s)/M-Ds=1.0751(56), f(D)*/f(D)=1.078(36), f(D)*s/f(Ds)=1.087(20), M-B*/M-B=1.0078(15), M-B*(s)/M-Bs=1.0083(10), f(B)*/f(B)=0.958(22) and f(B)*s/f(Bs)=0.974(10). Combining them with the experimental values of the pseudoscalar meson masses (used as input to fix the quark masses) and the values of the pseudoscalar decay constants calculated by ETMC, we get: M-D*=2013(14) MeV, M-D*(s)=2116(11) MeV, f(D)*=223.5(8.4) MeV, f(D)*(s)=268.8(6.6) MeV, M-B*=5320.5(7.6) MeV, M-B*(s)=5411.36(5.3) MeV, f(B)*=185.9(7.2) MeV and f(B)*(s)=223.1(5.4) MeV.