Here is a example code i used to find the phase-shift between two different signals. I used cross-correlation to determine the phase shift.
In line 26 you can see the crosscorr() function which performs, well of course you know what that does ;)
Hope this helps. If you made the code any better share it with me plz.
Good luck ;)