From fb28a85dc84830e0d729aafbe584abddb80a13be Mon Sep 17 00:00:00 2001
From: naotoY <naotino@gmail.com>
Date: Thu, 11 Jan 2024 14:03:18 +0900
Subject: [PATCH] [fix] IIR function
---
mbtrack2/tracking/rf.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/mbtrack2/tracking/rf.py b/mbtrack2/tracking/rf.py
index 10711aa..3d37de9 100644
--- a/mbtrack2/tracking/rf.py
+++ b/mbtrack2/tracking/rf.py
@@ -1492,10 +1492,11 @@ class ProportionalIntegralLoop():
T = self.ring.T1 * self.every
alpha = np.cos(omega*T)-1
tmp = alpha * alpha - 2 * alpha
+ print(tmp)
if tmp > 0:
self.IIRcoef = alpha + np.sqrt(tmp)
else:
- self.IIRcoef = self.every * self.ring.T1 * cutoff * 2 * np.pi
+ self.IIRcoef = T * cutoff * 2 * np.pi
self.IIRout = self.cav_res.Vc
def IIR(self,input):
--
GitLab