  In this article, we investigate the contribution of the high twist Feynman
diagrams to the large $p_{T}$ single pseudoscalar and vector mesons inclusive
production cross section in two-photon collisions and we present the general
formulae for the high and leading twist differential cross sections. The pion
wave function where two non-trivial Gegenbauer coefficients $a_2$ and $a_4$
have been extracted from the CLEO data, Braun-Filyanov pion wave function, the
asymptotic and the Chernyak-Zhitnitsky wave functions are used in the
calculations. For $\rho$-meson we used Ball-Braun wave function. The results of
all the calculations reveal that the high twist cross sections, the ratio R,
the dependence transverse momentum $p_{T}$ and the rapidity $y$ of meson in
$\Phi_{CLEO}(x,Q^2)$ wave function case is very close to the $\Phi_{asy}(x)$
asymptotic wave function case. It is shown that the high twist contribution to
the cross section depends on the choice of the meson wave functions.
