\section{Kinky string decay simulation.}

\hspace{1.0em}
For kinky string decay simulation we have assumed the two steps process:
\begin{enumerate}
\item split gluon $g\rightarrow q_1 \bar{q}_1$ and create two longitudinal 
strings;
\item decay longitudunal strings $q\bar{q}_1 \rightarrow h$ and
 $\bar{q}q_1 \rightarrow h$ into hadrons $h$.
\end{enumerate}
The production of $q_1\bar{q}_1$ is considered similar (the same sampling of 
quark flavors and the same $p_t$-distribution for the quarks) as a production 
of $q\bar{q}$-pairs during a longitudinal string decay. 
We use the $g\rightarrow q_1\bar{q}_1$ splitting function \cite{APKW80}:
\begin{equation}
\label{KSD1}f^{q}_{g}(z)=z^2+(1-z)^2,
\end{equation}
where $z=\frac{E_q+p^z_q}{E_g+p^z_g}$, derived by Altarelli and 
Parisi \cite{AP77}.