Charmonia production in p+p collisions under NRQCD formalism