From 59fc3da29b3659e0f3ccce29ba0de8db1abba1ac Mon Sep 17 00:00:00 2001
From: Kevin Kuns <kevin.a.kuns@gmail.com>
Date: Sat, 28 Mar 2020 14:44:28 -0400
Subject: [PATCH] Fix length to strain conversion when only high frequencies
are computed
Removed normalization of sinc_sqr in dhdl because it is one at DC. Previously the normalization was set by the first frequency point, which is wrong if only high frequencies are considered
---
gwinc/ifo/noises.py | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/gwinc/ifo/noises.py b/gwinc/ifo/noises.py
index 344969e7..555c01b1 100644
--- a/gwinc/ifo/noises.py
+++ b/gwinc/ifo/noises.py
@@ -112,13 +112,12 @@ def dhdl(f, armlen):
omega_arm_b = (1 + sin(nu_small)) * pi * f * armlen / c
sinc_sqr = 4 / abs(sin(omega_arm_f) * exp(-1j * omega_arm) / omega_arm_f
+ sin(omega_arm_b) * exp(1j * omega_arm) / omega_arm_b)**2
- # keep DC value equal to 1
- sinc_sqr /= sinc_sqr[0]
dhdl_sqr = sinc_sqr / armlen**2
return dhdl_sqr, sinc_sqr
##################################################
+
def precomp_mirror(f, ifo):
ifo.Materials.MirrorVolume = \
pi * ifo.Materials.MassRadius**2 \
--
GitLab