Suppose that two unit-norm vectors $\boldsymbol{a}\in \mathbb{R}^m$ and $\boldsymbol{b}\in\mathbb{R}^n$ are given with $m\leq n$. Furthermore, let $\boldsymbol{F}_{m,n}$ denote the first $m$ rows of the $n$-point DFT matrix. To bound the tail-probability of some random process, I need to upper bound $N\left(\mathcal{S},\left\Vert\cdot\right\Vert,\epsilon\right)$, the $\epsilon$-covering number of the set \begin{align*} \mathcal{S} =\left\{\boldsymbol{F}_{m,n}\odot\left(\boldsymbol{a}\boldsymbol{x}^\mathrm{T}+\boldsymbol{y}\boldsymbol{b}^\mathrm{T}\right)\mid \boldsymbol{x}\in\mathbb{R}^n,\boldsymbol{y}\in\mathbb{R}^m,\left\Vert \boldsymbol{a}\boldsymbol{x}^\mathrm{T}+\boldsymbol{y}\boldsymbol{b}^\mathrm{T}\right\Vert_{\mathrm{F}}\leq 1\right\}, \end{align*} where $\odot$ denotes the Hadamard product.

The upper bound that I managed to achieve is $\left(1+\frac{2}{\epsilon}\right)^{m+n-1}$. Is there any way that I can verify whether my naive bound is sharp or not?

To achieve the above bound I simply bounded the spectral norm of $\boldsymbol{F}_{m,n}\odot\boldsymbol{X}$ by the Frobenius-norm and reduced the covering to covering of a unit ball in $\mathbb{R}^{m+n-1}$. I asked here if there is any better bound for the spectral norm I mentioned, but I didn't get any response that is obviously better than the Frobenius norm.