In this paper we design and analyze a class of high order numerical methods to two dimensional Heaviside function integrals. Inspired by our high order numerical methods to two dimensional delta ...