#atec.dtel = 9.1 #telescope diameter m
atec.dtel = 2.0 * u.m #telescope diameter m
atec.dstel = 0.00 * u.m #secondary telescope diameter m
atec.throughput = 0.05
atec.ndark = 0.17 / u.s #dark current
atec.nread = 14.0 #nr
atec.fullwell = 80000.
target = exocounts.TargetClass()
target.name = "HR8799bcde"
target.teff = 1000.0 * u.K #K
target.rstar = 0.1 * const.R_sun #Rsolar
target.d = 39 * u.pc #pc
obs = exocounts.ObsClass(atec, target)
obs.texposure = 1.0 * 24.0 * 1.0 * u.h #1 d
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
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("photon noise [ppm]=", obs.sign_relative)
print("photon N for exp=", obs.nphoton_exposure)
#ost.dtel = 9.1 #telescope diameter m
ost.dtel = 10.0 #telescope diameter m
ost.dstel = 0.00 #secondary telescope diameter m
ost.throughput = 0.1
ost.ndark = 0.0 #dark current
ost.nread = 0.0 #nr
ost.fullwell = 80000.
target = exocounts.TargetClass()
target.name = "self luminous"
target.teff = 1300.0 #K
target.rstar = 0.1 #Rsolar
target.dpc = 40 #pc
obs = exocounts.ObsClass(ost, target)
obs.texposure = 5.0 #cadence [hour] # 30 x visits (1 hr=transit dur trappist e)
obs.tframe = 7.1 #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
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("photon noise [ppm]=", obs.sign_relative)
print("photon N for exp=", obs.nphoton_exposure)
spica = exocounts.InstClass()
spica.lamb = 15.0 #micron
spica.dlam = 10.0 #micron
spica.dtel = 2.5 #telescope diameter m
spica.dstel = 0.00 #secondary telescope diameter m
spica.throughput = 0.2
spica.ndark = 0.17 #dark current
spica.nread = 14.0 #nr
spica.fullwell = 80000.
target = exocounts.TargetClass()
target.teff = 5800.0 #K
target.rstar = 1.0 #Rsolar
target.dpc = 10.0 #pc
obs = exocounts.ObsClass(spica, target)
obs.texposure = 1.0 #cadence [hour]
obs.tframe = 7.1 #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
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("photon noise [ppm]=", obs.sign_relative)
print("photon N for exp=", obs.nphoton_exposure)
R=100000.0
SUBARU.dlam = 2.3/R*u.micron #micron
SUBARU.dtel = 8.2*u.m #telescope diameter m
SUBARU.dstel = 0.00*u.m #secondary telescope diameter m
SUBARU.throughput = 0.05
SUBARU.ndark = 0.0/u.s #dark current [e-/pix/s]
SUBARU.nread = 0.00 #nr [e-/pix/read]
SUBARU.fullwell = 1.e7
target=exocounts.TargetClass()
target.teff = 1300.0*u.K #K
target.rstar = 0.1*const.R_sun #Rsolar
target.d = 40.0*u.pc #pc
obs=exocounts.ObsClass(SUBARU,target)
obs.texposure = 1.0*u.h #= [hour]
obs.tframe = 10.0*u.s #time for one frame [sec]
obs.napix = 1 # number of the pixels in aperture
obs.mu = 1
obs.effnpix = 0.5
obs.target = target
obs.update()
magdict=convmag.get_magdict()
print("magnitude=",convmag.get_mag("V",obs.flux,magdict))
print(obs.nphoton_brightest/SUBARU.fullwell)
print(obs.sat)
print("photon/pix/frame=",obs.nphoton_frame)
target = exocounts.TargetClass()
target.name = "Alpha Cen A"
target.teff = 255.0 * u.K #K ## optimistic (no cloud)
target.rstar = 1.0 * const.R_earth
#tau_ceti_teff = 5344.0*u.K #K ## optimistic (no cloud)
#tau_ceti_rstar = 0.79*const.R_sun
#c=(((target.rstar)**2*nstar.Blunitless(target.teff,michi.lamb))/(nstar.Blunitless(tau_ceti_teff,michi.lamb)*(tau_ceti_rstar)**2)).to(1)
#print("contrast=",c)
#target.d = 3.65*u.pc #pc
target.d = 1.3 * u.pc #pc
obs = exocounts.ObsClass(michi, target)
obs.texposure = 100.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
obs.fgaparture = ((diflimit.ld(michi.lamb, michi.dtel)) / 2.0)**2 * np.pi
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("photon count for foreground", obs.nphoton_foreground)
ejas = exocounts.InstClass()
ejas.lamb = 1.4 * u.micron #micron
ejas.dlam = 0.6 * u.micron #micron
ejas.dtel = 0.31 * u.m #telescope diameter m
ejas.dstel = 0.09 * u.m #secondary telescope diameter m
ejas.throughput = 0.7
ejas.ndark = 60.0 / u.s #dark current [1/s]
ejas.nread = 30.0 #nr
ejas.fullwell = 80000.
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 = 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
target.dpc = 16.0 #change targets
obs.target = target
obs.update()
magdict = convmag.get_magdict()
print("H mag=", convmag.get_mag("H", obs.flux, magdict))
