We suggest the modified matching procedure for TMD PDF to the integrated PDF aimed to increase the amount of perturbative information in the TMD PDF expression. The procedure consists in the selection and usage of the non-minimal operator basis, whic
h restricts the expansion to desired general behavior. The implication of OPE allows to systematic account of the higher order corrections. In the case of TMD PDF we assume the Gaussian behavior, which suggests Laguerre polynomial basis as the best for the convergence of OPE. We present the leading and next-to-leading expression of TMD PDF in this basis. The obtained perturbative expression for the TMD PDF is valid in the wide region of $b_T$ (we estimate this region as $b_Tlesssim 2-3$ GeV$^{-1}$ depending on $x$).
