Warning: the following answer is (a) maybe a bit careless in a nineteenth-century sort of way, (b) missing its final step (which may be obvious to others, and/or amount to looking things up in the right table)
We are looking for an even function $g:\mathbb{R}\to \mathbb{R}$ with support on $[-1,1]$; we will define $f$ to be its Fourier transform. The condition $\int_{\mathbb{R}} f(x) dx=1$ becomes $g(0)=1$. What we need to minimize is the quantity $$I = \int_{\mathbb{R}} |x| \widehat{g}(x) dx.$$ Here are two attempts to express $I$ more directly in terms of $g$.
The Fourier transform of $\frac{g'(t)}{2\pi i}$ equals $x \widehat{g}(x)$. We can write $|x| \widehat{g}(x) = \mathrm{sgn}(x) x \widehat{g}(x)$. The Fourier transform of $\mathrm{sgn}(x)$ is $\frac{1}{i \pi t}$ (in some sense). Hence, the Fourier transform of $|x| \widehat{g}(x)$ should be the convolution of $\frac{1}{i \pi t}$ and $\frac{g'(-t)}{2\pi i} = \frac{g'(t)}{2\pi i}$. In particular, $I$ should equal the value of the Fourier transform of $|x| \widehat{g}(x)$ at $0$, i.e., $$I = \int_{\mathbb{R}} \frac{1}{i \pi t} \frac{g'(t)}{2\pi i} dt = - \frac{1}{2 \pi^2} \int_{\mathbb{R}} \frac{g'(t)}{t} dt.$$ We can assume $g'(0)=0$, so the integral above should make sense.
Since $I = - 2 \int_{-\infty}^0 x \widehat{g}(x) dx$, we can write $I = - 2\int_{-\infty}^0 \int_{-\infty}^y \widehat{g}(x) dx dy$. Now, the second antiderivative of $\widehat{g}(x)$ should have Fourier transform $\frac{g(t)}{(2\pi i)^2 t^2}$. Its value at $0$ equals both $-I/2$ (by definition) and $$\int_\mathbb{R} \frac{g(t)}{(2\pi i)^2 t^2} dt = -\frac{1}{4\pi^2} \int_{\mathbb{R}} \frac{g(t)}{t^2} dt;$$ hence, $$I = \frac{1}{2\pi^2} \int_{\mathbb{R}} \frac{g(t)}{t^2} dt.$$ Of course this diverges.
At the same time, by integration by parts, $-\int_\mathbb{R} \frac{g'(t)}{t} dt$ equals $\int_\mathbb{R} \frac{g(t)-1}{t^2} dt$, which converges. So, it looks like $$I =\frac{1}{2\pi^2} \int_\mathbb{R} \frac{g(t)-1}{t^2} dt.$$
We recall we must also fulfill the constraint that $f$ take non-negative values. This will certainly be true if $g$ is defined as a convolution $h\ast h$, with $h$ symmetric and real-valued, as then $\widehat{h}$ will be real valued, and $\widehat{h\ast h} = \widehat{h}^2$. We can require that $h$ also be an even function, and that the support of $h$ be contained in $[-1/2,1/2]$. I think these are all necessary conditions, so I am not making the search space for my optimum any smaller, but I'd be delighted if others can double-check and confirm.
So, we've reduced our problem to: find a symmetric function $h:[-1/2,1/2]\to \mathbb{R}$ with $(h\ast h)(0) = |h|_2^2 = 1$ such that $$\int_{\mathbb{R}} \frac{(h\ast h)(t)-1}{t^2} dt$$ is minimal.
A bit of calculus variations seems to show that the optimal $h(t)$ has to have $$\frac{1}{h(t_0)} \int \frac{\frac{1}{2} (h(t+t_0)+h(t-t_0)) - h(t_0)}{t^2} dt$$ equal to a constant independent of $t_0$, for $t_0\in (-1/2,1/2)$. Again by integration by parts, this is just $$\frac{1}{h(t_0)} \int \frac{\frac{1}{2} (h'(t+t_0)+h'(t-t_0))}{t} dt,$$ which equals $\frac{1}{h(t_0)} H(h')(t_0)$, where $H$ is our old friend the Hilbert transform.
In other words: we are to find a (continuous) function $h:\mathbb{R}\to \mathbb{R}$, supported on $[-1/2,1/2]$, with $|h|_2=1$, such that $H(h')(t) = \lambda h(t)$ for all $t\in (-1/2,1/2)$ and some $\lambda$.
Surely such a function must be known (if it exists)?