Abstract We report on the inclusive $$\text {J}/\psi $$ J / ψ production cross section measured at the CERN Large Hadron Collider in proton–proton collisions at a center-of-mass energy $$\sqrt{s}~=~13$$ s = 13  TeV. The $$\text {J}/\psi $$ J / ψ mesons are reconstructed in the $$\text {e}^{+}\text {e}^{-}$$ e + e - decay channel and the measurements are performed at midrapidity ( $$|y|<0.9$$ | y | < 0.9 ) in the transverse-momentum interval $$0<p_{\mathrm{T}} <40$$ 0 < p T < 40  GeV/ $$c$$ c , using a minimum-bias data sample corresponding to an integrated luminosity $$L_{\text {int}} = 32.2~\text {nb}^{-1}$$ L int = 32.2 nb - 1 and an Electromagnetic Calorimeter triggered data sample with $$L_{\text {int}} = 8.3~\mathrm {pb}^{-1}$$ L int = 8.3 pb - 1 . The $$p_{\mathrm{T}}$$ p T -integrated $$\text {J}/\psi $$ J / ψ production cross section at midrapidity, computed using the minimum-bias data sample, is $$\text {d}\sigma /\text {d}y|_{y=0} = 8.97\pm 0.24~(\text {stat})\pm 0.48~(\text {syst})\pm 0.15~(\text {lumi})~\mu \text {b}$$ d σ / d y | y = 0 = 8.97 ± 0.24 ( stat ) ± 0.48 ( syst ) ± 0.15 ( lumi ) μ b . An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and $$p_{\mathrm{T}}$$ p T -differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.