In this work we have employed Bethe-Salpeter equation (BSE) under covariant instantaneous ansatz (CIA) to study electromagnetic decays of ground state equal mass vector mesons: $\rho$, $\omega$, $\phi$, $\psi$ and $Y$ through the process $V\rightarrow\gamma*\rightarrow e^+ + e^-$. We employ the generalized structure of hadron-quark vertex function $\Gamma$ which incorporates various Dirac structures from their complete set order-by-order in powers of inverse of meson mass. The electromagnetic decay constants for the above mesons are calculated using the leading order (LO) and the next-to-leading order (NLO) Dirac structures. The relevance of various Dirac structures in this calculation is studied.