Background subtraction for high-fold gamma-ray coincidence data