I have read about convolution method in the Karlin & Taylor's book. I'm confused, why the domain of integral changed from $-\infty$ in to $z$?

Is it because $X$ and $Y$ non negative, so $Z=X+Y\geq 0$? So, the domain of integral changed?
On
If $\xi <0$ the $f_{X} (\xi)=0$ and if $\xi >z$ then $f_Y(z-\xi)=0$. So you only have integrate over the interval $0 <\xi <z$.
If $\Pr[X < 0] = 0$ and $\Pr[Y < 0] = 0$, then $f_X(x) = 0$ and $f_Y(y) = 0$ for $x < 0$ and $y < 0$. Consequently $f_Y(z - \xi) = 0$ if $\xi > z$, and $f_X(\xi) = 0$ if $\xi < 0$. Equation 2.13 then becomes $$f_Z(z) = \int_{\xi = -\infty}^\infty f_Y(z-\xi) f_X(\xi) \, d\xi = \int_{\xi = 0}^z f_Y(z - \xi) f_X(\xi) \, d\xi,$$ and if $z < 0$, this integral is $0$ since no $\xi \in \mathbb R$ satisfies the inequality $0 \le \xi \le z$.