Based on the data of BES and Belle, the production of DD in the e + e − → DD scattering process is studied in this paper. We analyze the continuum and resonant contributions in the energy region from 3.7 to 4.4 GeV. In the χ 2 fit to data, we obtain the resonance parameters of ψ(3770), the branching ratio of ψ(3770) → DD decay by confronting the data to the theoretical formula where both the contributions of the resonances, continuum and interference effects are included. We obtain the branching ratio of ψ(3770) → DD decay is 97.2% ± 8.9%, as well as the branching ratio of ψ(4040), ψ(4160) → DD decays. PACS numbers: 13.25.Gv, 13.66.Bc, 14.40.Gx Physics of e + e − annihilation at the energy region of 3 − 5 GeV is interesting. It attracts the focus of both experimental and theoretical studies. The physics at this energy region involves several well-established vector resonances, J/ψ, ψ(2S), ψ(3770), ψ(4040), ψ(4160), ··· , which are bound states of quark and antiquark pair cc. Studying the production and decays of these resonances can deepen our knowledge about dynamics of interactions between quarks. The resonance ψ(3770) has a mass just slightly above the threshold of DD pair production. Its decays evade the suppression of the Okubo-Zweig-Iizuka (OZI) rule [1]. This is consistent with the fact that the width of ψ(3770) is about 2 orders larger than those of J/ψ and ψ(2S) [2]. The J/ψ and ψ(2S) can only decay into non-DD final states, which is suppressed according to the OZI rule. It is believed that ψ(3770) decays dominantly into DD pair. All the well-established non-DD decay modes of ψ(3770) only show up with the branching ratios at the order of 10 −3 − 10 −4 . The measurement of many other decay modes of ψ(3770) only gives upper bounds, which shows that decay rates of these non-DD decay modes should be smaller than 10 −3 or 10 −4 [2,3]. The sum of all these well-measured decay rates is at most at the order of several percent.In the e + e − collider the properties of the resonance ψ(3770) are measured through the scattering process e + e − → ψ(3770) → f , where f can be any final states like DD or any other hadrons. The branching ratio of ψ(3770) → DD and ψ(3770) →non-DD can be derived from the measured scattering cross section of e + e − → DD and e + e − → hadrons. Both BES and CLEO-c Collaborations measured the cross section of e + e − → DD at the center-of-mass energy E c.m. = 3773 MeV several years ago [4, 5], and their results are in good agreement with each other.The CLEO-c Collaboration also measured the cross section of e + e − → ψ(3770) → hadrons at E c.m. = 3773 MeV [6]. The difference between this and the cross section of e + e − → ψ(3770) → DD is found to be (−0.01 ± 0.08 +0.41 −0.30 ) nb, which indicates that the decay rate of ψ(3770) to non-DD is tiny. However, the measurement of the BES Collaboration gives that the branching ratio of ψ(3770) → DD is (85.5 ± 1.7 ± 5.8)% or (83.6 ± 7.3 ± 4.2)%, and the decay branching ratio of ψ(3770) to non-DD is (14.5 ± 1.7 ± 5.8)% or (16.4 ± 7.3 ± 4.2)% [7,8],...