def test_ejas_bb():
""" unit test for blackbody photon count and other noises assuming JASMINE (2022 April)
"""
ejas = exocounts.InstClass()
ejas.lamb = 1.25 * u.micron #micron
ejas.dlam = 0.7 * u.micron #micron
ejas.dtel = 0.34 * u.m #telescope diameter m
ejas.dstel = 0.14 * u.m #secondary telescope diameter m or 12.4 (3 tels)
QE = 0.7
ejas.throughput = QE * 0.85 * 0.95
ejas.ndark = 15.5 / u.s #dark current
ejas.nread = 15.0 #nr
ejas.fullwell = 150000.
target = exocounts.TargetClass()
target.teff = 3000.0 * u.K #K
target.rstar = 0.2 * const.R_sun #Rsolar
target.d = 15.0 * u.pc #pc
obs = exocounts.ObsClass(ejas, target)
obs.texposure = 0.0833 * u.h #cadence [hour]
obs.tframe = 12.5 * u.s #time for one frame [sec]
obs.napix = 15 # number of the pixels in aperture
obs.mu = 1
S = 1.8 * 1.8 * np.pi #core size
obs.effnpix = S / 3.0 #3 is an approx. increment factor of PSF
obs.mu = 1
target.d = 16.0 * u.pc #change targets
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("H mag=", convmag.get_mag("H", obs.flux, magdict))
print("J mag=", convmag.get_mag("J", obs.flux, magdict))
print(
"Hw mag=", 0.7 * convmag.get_mag("J", obs.flux, magdict) +
0.3 * convmag.get_mag("H", obs.flux, magdict))
print("V mag=", convmag.get_mag("V", obs.flux, magdict))
print("=========================")
print("saturation?", obs.sat)
print("dark [ppm]=", obs.sigd)
print("readout [ppm]=", obs.sigr)
print("photon [ppm]=", obs.sign)
print("=========================")
print("photon relative=", obs.sign_relative)
assert not obs.sat
assert obs.sigd - 264.04942719119845 == 0.0
assert obs.sigr - 284.5480627240326 == 0.0
assert obs.sign - 3392.713732726309 == 0.0
assert obs.sign_relative - 294.74930064212117 == 0.0
from astropy import units as u import numpy as np ost = exocounts.InstClass() ost.lamb = 10.0 * u.micron #micron ost.dlam = 0.25 * u.micron #micron #ost.dtel = 9.1 #telescope diameter m ost.dtel = 5.9 * u.m #telescope diameter m ost.dstel = 0.00 * u.m #secondary telescope diameter m ost.throughput = 0.1 ost.ndark = 0.17 #dark current ost.nread = 14.0 #nr ost.fullwell = 80000. target = exocounts.TargetClass() target.name = "Trappist e" target.teff = 2559.0 * u.K #K target.rstar = 0.117 * const.R_sun #Rsolar target.d = 12.1 * u.pc #pc obs = exocounts.ObsClass(ost, target) obs.texposure = 30.0 * u.h #cadence [hour] # 30 x visits (1 hr=transit dur trappist e) obs.tframe = 7.1 * u.s #time for one frame [sec] obs.napix = 15 # number of the pixels in aperture obs.mu = 1 S = 1.8 * 1.8 * np.pi #core size obs.effnpix = S / 3.0 #3 is an approx. increment factor of PSF obs.mu = 1