This python routine should use obspy f-k analysis functionality to locate two seismic events relative to each other using cross-correlation trace waveforms with "coordinates" set as the slowness vectors describing how the waves leave the source region to the station.