The reaction N N → N N π offers a good testing ground for chiral effective field theory at intermediate energies. It challenges our understanding of the first inelastic channel in nucleon-nucleon scattering and of the charge-symmetry breaking pattern in hadronic reactions. In our previous studies, we presented a complete calculation of the pion-production operator for s-wave pions upto-and-including next-to-next-to-leading order (NNLO) in the formulation of chiral effective field theory, which includes pions, nucleons and ∆(1232) degrees of freedom. In this paper we calculate the near threshold cross section for the pp → dπ + reaction by performing the convolution of the obtained operators with nuclear wave functions based on modern phenomenological and chiral potentials. The available chiral N N wave functions are constructed with a cutoff comparable with the momentum transfer scale inherent in pion production reactions. Hence, a significant portion of the dynamical intermediate-range physics is thereby cut off by them. On the other hand, the NNLO amplitudes evaluated with phenomenological wave functions appear to be largely independent of the N N model used and give corrections to the dominant leading order contributions as expected from dimensional analysis. The result gives support to the counting scheme used to classify the pion production operators, which is a precondition for a reliable investigation of the chirally suppressed neutral pion production. The explicit inclusion of the ∆(1232) is found to be important but smaller than expected due to cancellations.