Objective Although HIFU has been successfully applied in various clinical applications in the past two decades for the ablation of many types of tumors, one bottleneck in its wider applications is the lack of a reliable and affordable strategy to guide the therapy. This study aims at estimating the therapeutic beam path at the pre-treatment stage to guide the therapeutic procedure. Methods An incident beam mapping technique using passive beamforming was proposed based on a clinical HIFU system and an ultrasound imaging research system. An optimization model was created to map the cross-like beam pattern by maximizing the total energy within the mapped area. This beam mapping technique was validated by comparing the estimated focal region with the HIFU-induced actual focal region (damaged region) through simulation, in-vitro, ex-vivo and in-vivo experiments. Results The results of this study showed that the proposed technique was, to a large extent, tolerant of sound speed inhomogeneities, being able to estimate the focal location with errors of 0.15 mm and 0.93 mm under in-vitro and ex-vivo situations respectively, and slightly over 1 mm under the in-vivo situation. It should be noted that the corresponding errors were 6.8 mm, 3.2 mm, and 9.9 mm respectively when the conventional geometrical method was used. Conclusion This beam mapping technique can be very helpful in guiding the HIFU therapy and can be easily applied in clinical environments with an ultrasound-guided HIFU system. Significance The technique is non-invasive and can potentially be adapted to other ultrasound-related beam manipulating applications.