added the function for optimal squeeze angle. Works for A+

gwinc/noise/quantum.py

@@ -5,6 +5,12 @@ import scipy.constants
 import logging
 from ..struct import Struct
 
+def sqzOptimalSqueezeAngle(Mifo, eta):
+    vHD = np.array([[sin(eta), cos(eta)]])
+    H = getProdTF(vHD, Mifo)[0]
+    alpha = arctan(abs(H[1]), abs(H[0]))
+    return alpha
+
 
 def shotrad(f, ifo):
     """Quantum noise