Skip to content
Snippets Groups Projects
Commit 59fc3da2 authored by Kevin Kuns's avatar Kevin Kuns
Browse files

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
parent 1906253b
No related branches found
No related tags found
No related merge requests found
......@@ -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 \
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment