For any $f\in L^1(\mathbb{R}^n)$ and $1\le j\le n$, recall that the Riesz transform $\mathcal{R}_jf\in L^{1,\infty}(\mathbb{R}^n)$ is given by $$ \mathcal{R}_jf:=c_n\lim_{\epsilon\to 0}\left(\frac{x_j}{|x|^{n+1}}1_{\mathbb{R}^n\setminus B_\epsilon(0)}(x)\right)*f, $$ for some irrelevant dimensional constant $c_n$. It is well-known that the above limit holds in the $L^{1,\infty}$-topology.
Assume now that, for some $f\in L^1(\mathbb{R}^n)$, we have $\mathcal{R}_j f\in L^1(\mathbb{R}^n)$. Can we conclude that the usual formula $\widehat{\mathcal{R}_j f}=-i\frac{\xi_j}{|\xi|}\widehat{f}(\xi)$ holds for this $f$?
This question is not trivial since it is not clear whether the above limit holds in the $L^1$-topology.
Moreover, given $f\in L^1(\mathbb{R}^n)$, is the following equivalence $$\mathcal{R}_j f\in L^1(\mathbb{R}^n) \ \Longleftrightarrow\ \mathcal{F}^{-1}\left(-i\frac{\xi_j}{|\xi|}\widehat{f}(\xi)\right)\in L^1(\mathbb{R}^n)$$ true? (The second formula is meant in the sense of tempered distributions)