of the order of $O(1/Q^2)$ if $x_B > 0.1$, not $O(A^{1/3}/Q^2)$.

On the other hand, the numerator in Eq. (\[LTavg\]) depends on two physically measured hadrons: the pion and the nucleus, and therefore, the OPE alone cannot be used to expand the numerator. However, the generalized factorization theorem introduced in Ref. [@QS_fac] can be used to factorize the numerator up to the $1/Q^2$ power corrections. Similarly to the Eq. (1) of Ref. [@LQS2], we expand the numerator in Eq. (\[LTavg\]) as $$\begin{aligned}
\int d\ell_{T}^2 \, \ell_{T}^2 \, 
\frac{d\sigma_{eA\rightarrow e\pi X}}{dx_B  dQ^2 d\ell_{T}^2} 
&=&
\sum_{a,c} \phi_{a/A}(x) \otimes 
   C^{(0)}_{ea\rightarrow ecX}\left(x_B/x,z=\ell/p_c,Q^2\right)
           \otimes D_{c\rightarrow\pi}(z)
\nonumber \\
&+&\frac{1}{Q^2}
\sum_{a,c} \left[
           T_{a/A}(x) \otimes 
   C^{(2)}_{ea\rightarrow ecX}\left(x_B/x,z=\ell/p_c,Q^2\right)
           \otimes D_{c\rightarrow\pi}(z) \right.
\nonumber \\ 
&& {\hskip 0.4in} \left.
          +\phi_{a/A}(x) \otimes 
\bar{C}^{(2)}_{ea\rightarrow ecX}\left(x_B/x,z=\ell/p_c,Q^2\right)
           \otimes d_{c\rightarrow\pi}(z) \right]
\nonumber \\
&+& ... 
\nonumber \\
&\equiv &
H^{(0)}_A + H^{(2)}_A + \bar{H}^{(2)}_A 
+ ...\, ,
\label{Hexp}\end{aligned}$$ where $...$ represents terms further suppressed in $1/Q^2$, $\otimes $ represents the convolution over partons’ momentum fractions, and explicit dependence on the factorization and/or renormalization scale is suppressed for simplicity. In Eq. (\[Hexp\]), $C^{(0)}$, $C^{(2)}$, and $\bar{C}^{(2)}$ are perturbatively calculable short-distance hard parts [@QS_fac], and $D_{c\rightarrow \pi}$ and $d_{c\rightarrow \pi}$ are twist-2 and twist-4 parton-to-pion fragmentation functions, respectively. 