[guided]The purpose of this step is to compute the one quantity that controls the local power of the Wilcoxon signed-rank statistic: the derivative at the null of its limiting signed-rank mean.
The limiting signed-rank score is the map $g:\mathbb R\to[-1,1]$ given by
\begin{align*}
g(y)=\operatorname{sgn}(y)\,H(|y|),
\end{align*}
where $H:[0,\infty)\to[0,1]$ is given by
\begin{align*}
H(t)=2\int_0^t f(x)\,d\mathcal L^1(x)
\end{align*}
and is the distribution function of $|\varepsilon_1|$ under the null law. Under a location shift $\theta$, the density of $Y_1$ is $y\mapsto f(y-\theta)$, so the limiting signed-rank mean is
\begin{align*}
m(\theta)=\int_{\mathbb R} g(y)f(y-\theta)\,d\mathcal L^1(y).
\end{align*}
By the differentiability assumption in the theorem, we may differentiate this functional at $\theta=0$. The finite Fisher information hypothesis supplies the standard location-score integrability control behind this differentiation, while the theorem explicitly assumes the local signed-rank differentiation itself. Since differentiating $f(y-\theta)$ with respect to $\theta$ gives $-f'(y-\theta)$, we obtain
\begin{align*}
m'(0)
=
-\int_{\mathbb R} g(y)f'(y)\,d\mathcal L^1(y).
\end{align*}
Now we transfer the derivative from $f$ onto $g$. This is where the evenness of $f$ enters. For $y>0$, the score is $g(y)=H(y)$, hence
\begin{align*}
g'(y)=H'(y)=2f(y)
\quad\text{for }\mathcal L^1\text{-a.e. }y>0.
\end{align*}
For $y<0$, the score is $g(y)=-H(-y)$. Differentiating this expression gives
\begin{align*}
g'(y)=H'(-y)=2f(-y).
\end{align*}
Since $f$ is even, $f(-y)=f(y)$, and therefore
\begin{align*}
g'(y)=2f(y)
\quad\text{for }\mathcal L^1\text{-a.e. }y<0.
\end{align*}
The boundary terms in [integration by parts](/theorems/2098) vanish. At $y=0$, we have $H(0)=0$, so $g(0)=0$. At infinity, $|g(y)|\le1$ and $f(y)\to0$ as $|y|\to\infty$, so the endpoint contribution $g(y)f(y)$ tends to $0$ along both half-lines. To pass from local integration by parts to the full line, fix $a<0<b$. On $[a,0]$ and $[0,b]$, the functions are absolutely continuous, so integration by parts gives
\begin{align*}
-\int_a^b g(y)f'(y)\,d\mathcal L^1(y)= -g(b)f(b)+g(a)f(a)+\int_a^b g'(y)f(y)\,d\mathcal L^1(y),
\end{align*}
with no boundary contribution at $0$ because $g(0)=0$. Now let $a\to-\infty$ and $b\to\infty$. The endpoint terms tend to $0$ since $|g(y)|\le1$ and $f(y)\to0$ as $|y|\to\infty$. The integral containing $g'f$ has a finite full-line limit because $g'(y)f(y)=2f(y)^2$ for $\mathcal L^1$-a.e. $y\ne0$ and $f\in L^2(\mathbb R)$. Thus $|g'f|\in L^1(\mathbb R)$. By the [Dominated Convergence Theorem](/theorems/4), applied to $\mathbb{1}_{[a,b]}(y)g'(y)f(y)$ dominated by $|g'(y)f(y)|$, the finite-interval integrals of $g'f$ converge to the full-line integral as $a\to-\infty$ and $b\to\infty$. The left-hand side is interpreted through the assumed differentiability formula for $m'(0)$, whose score-integrability requirements are part of the finite-Fisher-information and local asymptotic framework in the statement. Hence the finite-interval identity passes to the full line and gives
\begin{align*}
-\int_{\mathbb R} g(y)f'(y)\,d\mathcal L^1(y)=\int_{\mathbb R} g'(y)f(y)\,d\mathcal L^1(y).
\end{align*}
Substituting $g'(y)=2f(y)$ for $\mathcal L^1$-a.e. $y\ne0$ yields
\begin{align*}
m'(0)
=
2\int_{\mathbb R} f(y)^2\,d\mathcal L^1(y).
\end{align*}
This is the analytic core of the proof: the local signed-rank drift is determined by the $L^2$ mass of the density.[/guided]