We present results on transverse momentum ([Formula: see text]) and rapidity ([Formula: see text]) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive [Formula: see text] and [Formula: see text] at forward rapidity ([Formula: see text]) as well as [Formula: see text]-to-[Formula: see text] cross section ratios. These quantities are measured in pp collisions at center of mass energies [Formula: see text] and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at [Formula: see text], 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full [Formula: see text] range, provided that both contributions are summed. In particular, it is found that for [Formula: see text] GeV/c the non-prompt contribution reaches up to 50% of the total charmonium yield.