NRQCD factorization for twist-2 light-cone wave-functions of charmonia