esbo_ds/ESBO-ETC
esbo_ds/ESBO-ETC
11
0
Fork 1
Code Issues 1 Pull Requests Releases 1 Wiki Activity

Bugfix: divide by forward efficiency

Browse Source
This commit is contained in:
Lukas Klass 2020-10-12 14:22:12 +02:00
parent f205cc011e
commit 7fcb1850a6
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
esbo_etc/classes/sensor/Heterodyne.py
View File

@ -302,7 +302,7 @@ class Heterodyne(ASensor):
equivalencies=u.spectral_density( equivalencies=u.spectral_density(
background.wl))) background.wl)))
t_background = background * ( t_background = background * (
self.__main_beam_efficiency * background.wl ** 2 / (2 * k_B) * self.__eta_fss * u.sr) self.__main_beam_efficiency * background.wl ** 2 / (2 * k_B) / self.__eta_fss * u.sr)
t_background = SpectralQty(t_background.wl, t_background.qty.decompose()) t_background = SpectralQty(t_background.wl, t_background.qty.decompose())
# Calculate the incoming photon current of the target # Calculate the incoming photon current of the target
logger.info("Calculating the signal temperature.") logger.info("Calculating the signal temperature.")
@ -311,13 +311,13 @@ class Heterodyne(ASensor):
signal = SpectralQty(signal.wl, signal.qty.to(u.W / (u.m ** 2 * u.Hz), signal = SpectralQty(signal.wl, signal.qty.to(u.W / (u.m ** 2 * u.Hz),
equivalencies=u.spectral_density(signal.wl))) equivalencies=u.spectral_density(signal.wl)))
t_signal = signal * (self.__aperture_efficiency * self.__common_conf.d_aperture() ** 2 * t_signal = signal * (self.__aperture_efficiency * self.__common_conf.d_aperture() ** 2 *
np.pi / 4 / (2 * k_B) * self.__eta_fss) np.pi / 4 / (2 * k_B) / self.__eta_fss)
t_signal = SpectralQty(t_signal.wl, t_signal.qty.decompose()) t_signal = SpectralQty(t_signal.wl, t_signal.qty.decompose())
else: else:
signal = SpectralQty(signal.wl, signal.qty.to(u.W / (u.m ** 2 * u.Hz * u.sr), signal = SpectralQty(signal.wl, signal.qty.to(u.W / (u.m ** 2 * u.Hz * u.sr),
equivalencies=u.spectral_density(signal.wl))) equivalencies=u.spectral_density(signal.wl)))
t_signal = signal * (self.__main_beam_efficiency * signal.wl ** 2 / ( t_signal = signal * (
2 * k_B) * self.__eta_fss * u.sr) self.__main_beam_efficiency * signal.wl ** 2 / (2 * k_B) / self.__eta_fss * u.sr)
t_signal = SpectralQty(t_signal.wl, t_signal.qty.decompose()) t_signal = SpectralQty(t_signal.wl, t_signal.qty.decompose())
if self.__lambda_local_oscillator is None: if self.__lambda_local_oscillator is None: