[step:Verify the $b$-weak boundedness property by direct kernel estimation]
Let $R > 0$ and let $\varphi^x_R, \psi^x_R \in C_c^\infty(\mathbb{R})$ be bump functions adapted at scale $R$ to a point $x \in \mathbb{R}$: $\operatorname{supp}\varphi^x_R, \operatorname{supp}\psi^x_R \subseteq [x - R, x + R]$, with the normalisations $\|\varphi^x_R\|_\infty, \|\psi^x_R\|_\infty \le 1$ and $\|(\varphi^x_R)'\|_\infty, \|(\psi^x_R)'\|_\infty \le R^{-1}$. The $b$-weak boundedness property requires
\begin{align*}
|\langle C_A(b\varphi^x_R), b\psi^x_R\rangle| \le C(M)\, R \quad \text{(uniformly in } x \text{ and } R\text{).}
\end{align*}
We verify this by direct estimation of the Cauchy kernel against the bump support, **avoiding any appeal to a Schur-test argument on $|x-y|^{-1}$** (which would diverge logarithmically — this is the same obstruction that makes WBP a non-trivial property in the first place). The standard verification pairs the antisymmetric structure of the Cauchy kernel with the smoothness of the bumps.
Write $z(y) := y + iA(y)$, so the kernel is $K(x,y) = \frac{1 + iA'(y)}{z(x) - z(y)}$ and $b(y)\,d\mathcal{L}^1(y) = dz(y)$. Expanding the pairing,
\begin{align*}
\langle C_A(b\varphi^x_R), b\psi^x_R\rangle = \mathrm{p.v.}\int_{\mathbb{R}}\!\int_{\mathbb{R}} \frac{b(u)\,\psi^x_R(u)\,b(y)\,\varphi^x_R(y)}{z(u) - z(y)}\,d\mathcal{L}^1(y)\,d\mathcal{L}^1(u),
\end{align*}
and, since $b(u)\,d\mathcal{L}^1(u) = dz(u)$ and similarly for $y$,
\begin{align*}
\langle C_A(b\varphi^x_R), b\psi^x_R\rangle = \mathrm{p.v.}\int_{\mathbb{R}}\!\int_{\mathbb{R}} \frac{\varphi^x_R(y)\,\psi^x_R(u)}{z(u) - z(y)}\,dz(y)\,dz(u).
\end{align*}
Antisymmetrising the integrand in $(y, u)$ — exploiting $1/(z(u) - z(y)) = -1/(z(y) - z(u))$ — gives
\begin{align*}
\langle C_A(b\varphi^x_R), b\psi^x_R\rangle = \tfrac{1}{2}\,\mathrm{p.v.}\int\!\!\int \frac{\varphi^x_R(y)\,\psi^x_R(u) - \varphi^x_R(u)\,\psi^x_R(y)}{z(u) - z(y)}\,dz(y)\,dz(u).
\end{align*}
The numerator vanishes on the diagonal $y = u$, and by the Lipschitz property of $\varphi^x_R, \psi^x_R$ (both with derivatives $\le R^{-1}$ in $L^\infty$),
\begin{align*}
|\varphi^x_R(y)\,\psi^x_R(u) - \varphi^x_R(u)\,\psi^x_R(y)| \le 2\,R^{-1}\,|y - u|.
\end{align*}
Combined with $|z(u) - z(y)| \ge |y - u|/\sqrt{1 + M^2}$, the antisymmetrised integrand-modulus (with respect to *real* Lebesgue measure $d\mathcal{L}^1(y)\,d\mathcal{L}^1(u)$) satisfies
\begin{align*}
\Big|\frac{\varphi^x_R(y)\,\psi^x_R(u) - \varphi^x_R(u)\,\psi^x_R(y)}{z(u) - z(y)}\,\frac{dz(y)\,dz(u)}{d\mathcal{L}^1(y)\,d\mathcal{L}^1(u)}\Big| \le 2\,R^{-1}\sqrt{1 + M^2}\cdot |b(y)|\,|b(u)|,
\end{align*}
where we used $|dz/d\mathcal{L}^1| = |b| \le 1 + M$ pointwise (i.e. $|b(y)| = |1 + iA'(y)| \le 1 + M$ a.e., not the global $\|b\|_\infty$ as a multiplicative constant, but pointwise bounded by the same number $1 + M$). The integrand is supported in $[x - R, x + R]^2$, which has $\mathcal{L}^2$-measure $4 R^2$. Therefore
\begin{align*}
|\langle C_A(b\varphi^x_R), b\psi^x_R\rangle| \le \tfrac{1}{2}\cdot 2\,R^{-1}\sqrt{1 + M^2}\cdot (1 + M)^2\cdot 4 R^2 = 4\,R\,\sqrt{1 + M^2}\,(1 + M)^2 =: C(M)\,R.
\end{align*}
The factor $(1 + M)^2$ — pointwise upper bound for $|b(y)\,b(u)|$ — replaces the spurious $\|b\|_\infty^2$ multiplicative-constant claim; the substitution $dz = b\,d\mathcal{L}^1$ is a *pointwise* substitution that introduces $b$ into the integrand, and the integrand bound uses $|b(y)| \le 1 + M$ pointwise. This is the WBP at scale $R$ around $x$ with constant $C(M) = 4\sqrt{1 + M^2}\,(1 + M)^2$, depending only on $M$.
[/step]