Operational earthquake forecasting protocols commonly use statistical models for their recognized ease of implementation and robustness in describing the short-term spatiotemporal patterns of triggered seismicity. However, recent advances on physics-based aftershock forecasting reveal comparable performance to the standard statistical counterparts with significantly improved predictive skills when fault and stress field heterogeneities are considered. Here, we perform a pseudo-prospective forecasting experiment during the first month of the 2019 Ridgecrest (California) earthquake sequence. We develop seven Coulomb rate-and-state models that couple static stress change estimates with continuum mechanics expressed by the rate-and-state friction laws. Our model parametrization supports a gradually increasing complexity; we start from a preliminary model implementation with simplified slip distributions and spatially homogeneous receiver faults to reach an enhanced one featuring optimized fault constitutive parameters, finite-fault slip models, secondary triggering effects, and spatially heterogenous planes informed by pre-existing ruptures. The data-rich environment of Southern California allows us to test whether incorporating data collected in near real-time during an unfolding earthquake sequence boosts our predictive power. We assess the absolute and relative performance of the forecasts by means of statistical tests used within the Collaboratory for the Study of Earthquake Predictability (CSEP) and compare their skills against a standard benchmark ETAS model for the short (24 hours after the two Ridgecrest mainshocks) and intermediate-term (one month). Stress-based forecasts expect heightened rates along the whole near-fault region and increased expected seismicity rates in Central Garlock Fault. Our comparative model evaluation supports that faulting heterogeneities coupled with secondary triggering effects are the most critical success components behind physics-based forecasts, but also underlines the importance of model updates incorporating near real-time available aftershock data reaching better performance than ETAS models. We explore the physical basis behind our results by investigating the localized shut down of pre-existing normal faults in the Ridgecrest near-source area.