In "Solution Property Reconstruction for Finite Volume scheme: a BVD+MOOD framework", Int. J. Numer. Methods Fluids, 2020, we have designed a novel solution property preserving reconstruction, so-called multi-stage BVD-MOOD scheme. The scheme is able to maintain a high accuracy in smooth profile, eliminate the oscillations in the vicinity of discontinuity, capture sharply discontinuity and preserve some physical properties like the positivity of density and pressure for the Euler equations of compressible gas dynamics. In this paper, we present an extension of this approach for the compressible Euler equations supplemented with source terms (e.g., gravity, chemical reaction). One of the main challenges when simulating these models is the occurrence of negative density or pressure during the time evolution, which leads to a blow-up of the computation. General compressible Euler equations with different type of source terms are considered as models for physical situations such as detonation waves. Then, we illustrate the performance of the proposed scheme via a numerical test suite including genuinely demanding numerical tests. We observe that the present scheme is able to preserve the physical properties of the numerical solution still ensuring robustness and accuracy when and where appropriate.Read less <