diff --git a/lisainstrument/noises.py b/lisainstrument/noises.py index 5507df9e27dd974def3386a075dc94b57410d26d..f990775785395e80fa824aa809f64fecedcd8362 100644 --- a/lisainstrument/noises.py +++ b/lisainstrument/noises.py @@ -134,7 +134,7 @@ def infrared(fs, size, asd): return np.cumsum(red_noise) * (2 * pi / fs) -def laser(fs, size, asd=30, fknee=2E-3, shape='white+infrared'): +def laser(fs, size, asd, shape): """Generate laser noise [Hz]. This is a white noise with an infrared relaxation towards low frequencies, @@ -153,6 +153,7 @@ def laser(fs, size, asd=30, fknee=2E-3, shape='white+infrared'): fknee: cutoff frequency [Hz] shape: spectral shape, either 'white' or 'white+infrared' """ + fknee = 2E-3 logger.debug( "Generating laser noise (fs=%s Hz, size=%s, asd=%s " "Hz/sqrt(Hz), fknee=%s Hz, shape=%s)", @@ -165,7 +166,7 @@ def laser(fs, size, asd=30, fknee=2E-3, shape='white+infrared'): raise ValueError(f"invalid laser noise spectral shape '{shape}'") -def clock(fs, size, asd=6.32E-14): +def clock(fs, size, asd): """Generate clock noise fluctuations [ffd]. The power spectral density in fractional frequency deviations is a pink noise, @@ -179,7 +180,7 @@ def clock(fs, size, asd=6.32E-14): return pink(fs, size, asd) -def modulation(fs, size, asd=5.2E-14): +def modulation(fs, size, asd): """Generate modulation noise [ffd]. The power spectral density as fractional frequency deviations reads @@ -195,7 +196,7 @@ def modulation(fs, size, asd=5.2E-14): return powerlaw(fs, size, asd, 1/3) -def backlink(fs, size, asd=3E-12, fknee=2E-3): +def backlink(fs, size, asd, fknee): """Generate backlink noise as fractional frequency deviation [ffd]. The power spectral density in displacement is given by @@ -220,7 +221,7 @@ def backlink(fs, size, asd=3E-12, fknee=2E-3): + red(fs, size, 2 * pi * asd * fknee**2 / c) -def ranging(fs, size, asd=3E-9): +def ranging(fs, size, asd): """Generate stochastic ranging noise [s]. This is a white noise as a timing jitter, @@ -234,7 +235,7 @@ def ranging(fs, size, asd=3E-9): return white(fs, size, asd) -def testmass(fs, size, asd=2.4E-15, fknee=0.4E-3): +def testmass(fs, size, asd, fknee): """Generate test-mass acceleration noise [m/s]. Expressed in acceleration, the noise power spectrum reads @@ -299,7 +300,7 @@ def jitter(fs, size, asd): logger.debug("Generating jitter (fs=%s Hz, size=%s, asd=%s rad/sqrt(Hz))", fs, size, asd) return violet(fs, size, 2 * pi * asd) -def dws(fs, size, asd=7E-8/335): +def dws(fs, size, asd): """Generate DWS measurement noise. The power spectral density in angle is given by