From 34ffbd9b554a2634dd4a207cba220815ede12083 Mon Sep 17 00:00:00 2001
From: Lee McCuller <Lee.McCuller@gmail.com>
Date: Fri, 24 Aug 2018 11:21:02 -0400
Subject: [PATCH] added the function for optimal squeeze angle. Works for A+

---
 gwinc/noise/quantum.py | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/gwinc/noise/quantum.py b/gwinc/noise/quantum.py
index 73f0112d..e974f048 100644
--- a/gwinc/noise/quantum.py
+++ b/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
-- 
GitLab