The so-called ρπ puzzle of J/ψ and ψ(2S) decays is examined using the experimental data available to date. Two different approaches were taken to estimate the ratio of J/ψ and ψ(2S) hadronic decay rates. While one of the estimates could not yield the exact ratio of ψ(2S) to J/ψ inclusive hadronic decay rates, the other, based on a computation of the inclusive ggg decay rate for ψ(2S) (J/ψ) by subtracting other decay rates from the total decay rate, differs by two standard deviations from the naive prediction of perturbative QCD, even though its central value is nearly twice as large as what was naively expected. A comparison between this ratio, upon making corrections for specific exclusive two-body decay modes, and the corresponding experimental data confirms the puzzles in J/ψ and ψ(2S) decays. We find from our analysis that the exclusively reconstructed hadronic decays of the ψ(2S) account for only a small fraction of its total decays, and a ratio exceeding the above estimate should be expected to occur for a considerable number of the remaining decay channels. We also show that the recent new results from the BES experiment provide crucial tests of various theoretical models proposed to explain the puzzle.PACS numbers: 13.25. Gv, 12.38.Qk One of the outstanding problems in heavy quarkonium physics is the strong suppression of the ψ(2S) decays to vector plus pseudoscalar-meson (VP) final states, ρπ and K + K * − + c.c., which is referred to as the ρπ puzzle [1]. Following the first observation of this anomaly [2], meagre experimental progress was made over the years, and theoretical analysis based on limited data often lead to unsatisfactory, sometimes premature, inferences. The situation has been changed dramatically in the last few years. A wealth of interesting new information, which extended the puzzle considerably, has emerged from intense studies of ψ(2S) hadronic decays at the BES experiment, using a large sample of 3.79 million ψ(2S) decays [3]. It is hoped for that new concerted efforts on both the theoretical and experimental side would eventually lead to a solution of this long-standing conundrum.In this paper we seek to examine the ρπ puzzle based purely on existing experimental data. We begin with an analysis for estimating the ratio of hadronic decay rates of J/ψ and ψ(2S), which we shall denote by Q, by using the data compiled by the Particle Data Group [4], in an attempt to avoid as many theoretical ambiguities as possible in the analysis. Two different approaches to this estimate are performed. First we compare the results between themselves, and then the naive prediction of perturbative QCD (PQCD) is used. Subsequently, possible corrections to Q are discussed, as they associate with specific exclusive decay modes, and the corrected values of Q are used as standards to compare with the corresponding experimental data. Comments on the issue of J/ψ and ψ(2S) decays to multi-hadron final states and on the potential similarity of the η c -η c (2S) decays to the ρπ puzzle are profusely added. F...