[guided]We apply the previously established [Hörmander weighted $L^2$ existence theorem for $\bar{\partial}$](/page/Hormander%20Weighted%20L2%20Existence%20Theorem). This is a prerequisite theorem, not the statement currently being proved. It requires a pseudoconvex domain, a plurisubharmonic $C^2$ weight, and a $\bar{\partial}$-closed form with finite weighted $L^2$ norm. These are exactly the hypotheses in the present situation, except that we should explicitly check plurisubharmonicity of the weight. The inequality $i\partial\bar{\partial}\varphi \ge c\omega$ with $c>0$ says that the complex Hessian of $\varphi$ is positive semidefinite, indeed uniformly positive, so $\varphi$ is plurisubharmonic.
Hörmander's theorem therefore gives a function $u \in L^2(\Omega,e^{-\varphi})$ such that $\bar{\partial}u=f$ in the distributional sense and
\begin{align*}
\int_\Omega |u|^2 e^{-\varphi}\,d\mathcal L^{2n}
&\le
\int_\Omega \langle (i\partial\bar{\partial}\varphi)^{-1}f,f\rangle_\omega e^{-\varphi}\,d\mathcal L^{2n}.
\end{align*}
The remaining point is to compute the constant. For each point $z \in \Omega$, let $T_z^{0,1*}\Omega$ denote the complex [vector space](/page/Vector%20Space) of $(0,1)$-covectors at $z$, with Hermitian inner product induced by the Euclidean Kähler metric $\omega$. Define
\begin{align*}
A_z: T_z^{0,1*}\Omega &\to T_z^{0,1*}\Omega
\end{align*}
to be the positive Hermitian endomorphism representing the Hermitian form $i\partial\bar{\partial}\varphi$ on $(0,1)$-covectors with respect to this $\omega$-inner product. The inequality $i\partial\bar{\partial}\varphi \ge c\omega$ says that every eigenvalue of $A_z$ is at least $c$. Taking inverses reverses the order of positive Hermitian endomorphisms, so every eigenvalue of $A_z^{-1}$ is at most $1/c$. Therefore, for the coefficient vector of the $(0,1)$-form $f$,
\begin{align*}
\langle (i\partial\bar{\partial}\varphi)^{-1}f,f\rangle_\omega
\le \frac{1}{c}|f|_\omega^2.
\end{align*}
Multiplying this pointwise inequality by the non-negative weight $e^{-\varphi}$ and integrating with respect to $\mathcal L^{2n}$ gives
\begin{align*}
\int_\Omega |u|^2 e^{-\varphi}\,d\mathcal L^{2n}
\le
\frac{1}{c}
\int_\Omega |f|_\omega^2 e^{-\varphi}\,d\mathcal L^{2n}.
\end{align*}
This is precisely the asserted estimate.[/guided]