Let
$T_{1}$,
$T_{2}$ be two Calderón–Zygmund operators and
$T_{1,b}$ be the commutator of
$T_{1}$ with symbol
$b\in \text{BMO}(\mathbb{R}^{n})$. In this paper, by establishing new bilinear sparse dominations and a new weighted estimate for bilinear sparse operators, we prove that the composite operator
$T_{1}T_{2}$ satisfies the following estimate: for
$\unicode[STIX]{x1D706}>0$ and weight
$w\in A_{1}(\mathbb{R}^{n})$,
$$\begin{eqnarray}\displaystyle & & \displaystyle w(\{x\in \mathbb{R}^{n}:\,|T_{1}T_{2}f(x)|>\unicode[STIX]{x1D706}\})\nonumber\\ \displaystyle & & \displaystyle \qquad \lesssim [w]_{A_{1}}[w]_{A_{\infty }}\log (\text{e}+[w]_{A_{\infty }})\int _{\mathbb{R}^{n}}\frac{|f(x)|}{\unicode[STIX]{x1D706}}\log \bigg(\text{e}+\frac{|f(x)|}{\unicode[STIX]{x1D706}}\bigg)w(x)\,dx,\nonumber\end{eqnarray}$$ while the composite operator
$T_{1,b}T_{2}$ satisfies
$$\begin{eqnarray}\displaystyle & & \displaystyle w(\{x\in \mathbb{R}^{n}:\,|T_{1,b}T_{2}f(x)|>\unicode[STIX]{x1D706}\})\nonumber\\ \displaystyle & & \displaystyle \qquad \lesssim [w]_{A_{1}}[w]_{A_{\infty }}^{2}\log (\text{e}+[w]_{A_{\infty }})\int _{\mathbb{R}^{n}}\frac{|f(x)|}{\unicode[STIX]{x1D706}}\log ^{2}\bigg(\text{e}+\frac{|f(x)|}{\unicode[STIX]{x1D706}}\bigg)w(x)\,dx.\nonumber\end{eqnarray}$$