Added analysis function

stitch_project/api/Stitch/Stitchtility.py

+import numpy as np
+
 def frequencyHzToDphi(frequency):
     dphi = int(frequency/0.18626451492)
     if dphi <= 1:
@@ -24,3 +26,25 @@ def numPeriodsToNringup(numPeriods, frequencyHz):
 
 def numPeriodsToNmeas(numPeriods, frequencyHz):
     return int(numPeriods*100.0e6/frequencyHz)
+
+#this spits out nm and rad
+def dataToAmpsPhase(data):
+    amps = []
+    numdata = [int.from_bytes(x[0], byteorder='little') for x in data[0]]
+    cosdata = [int.from_bytes(x[0], byteorder='little') for x in data[1]]
+    sindata = [int.from_bytes(x[0], byteorder='little') for x in data[2]]
+
+    for i in range(len(cosdata)):
+        if cosdata[i]&(2**63) != 0:
+            cosdata[i] = cosdata[i]-2**64
+    for i in range(len(sindata)):
+        if sindata[i]&(2**63) != 0:
+            sindata[i] = sindata[i]-2**64
+
+    amps = [2**-15*np.sqrt(2.*cosdata[i]**2+2.*sindata[i]**2)/numdata[i] for i in range(len(numdata))] #probably pm
+    angles = [2.*np.arctan(sindata[i]/cosdata[i]) for i in range(len(numdata))]
+    angles = np.angle(np.array(cosdata)+1j*np.array(sindata))
+    angles_unwr = np.unwrap(angles)
+    cosportion = [2**-15*np.abs(np.sqrt(2)*cosdata[i])/numdata[i] for i in range(len(amps))]
+    sinportion = [2**-15*np.abs(np.sqrt(2)*sindata[i])/numdata[i] for i in range(len(amps))]
+    return [amps,cosportion,sinportion], angles_unwr
\ No newline at end of file