Combining source direction vectors with wavefield decomposition to calculate angle gathers from reverse time migration