Ejemplos de InstClass en Python

Ejemplo n.º 1

import exocounts
import convmag
import numpy as np
from astropy import constants as const
from astropy import units as u

atec = exocounts.InstClass()
atec.lamb = 1.7 * u.micron  #micron
R = 100000.0
atec.dlam = atec.lamb / R  #micron
#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

Ejemplo n.º 2

import exocounts
import convmag

import numpy as np

ost = exocounts.InstClass()
ost.lamb = 2.3  #micron
ost.dlam = ost.lamb / 100000  #micron
#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

Ejemplo n.º 3

Archivo: spica.py Proyecto: HajimeKawahara/exocounts

import exocounts
import convmag

import numpy as np

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

Ejemplo n.º 4

import exocounts
import planet
import convmag
from astropy import constants as const
from astropy import units as u
import numpy as np

SUBARU=exocounts.InstClass()
SUBARU.lamb = 2.3*u.micron #micron
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

Ejemplo n.º 5

Archivo: tao.py Proyecto: HajimeKawahara/exocounts

import exocounts
import convmag
import diflimit
import nstar
import numpy as np
from astropy import constants as const
from astropy import units as u

michi = exocounts.InstClass()
michi.lamb = 10.0 * u.micron  #micron
michi.dlam = 1.0 * u.micron  #micron
#michi.dtel = 9.1 #telescope diameter m
michi.dtel = 6.5 * u.m  #telescope diameter m

michi.dstel = 0.00 * u.m  #secondary telescope diameter m
#michi.throughput = 0.8*0.9*0.3 #QE x Efficiency x Inst throughtput
michi.throughput = 0.5  #QE x Efficiency x Inst throughtput

michi.ndark = 0.0 / u.s  #dark current [1/s]
michi.nread = 0.0  #nr
michi.fullwell = 80000.

fgunit = u.s * u.m * u.m * u.arcsec * u.arcsec * u.micron
michi.fgtel = 0.05 * 1e10 / fgunit  #pt/s/m2/arcsec2/um
michi.fgatm = 0.05 * 1e10 / fgunit  #pt/s/m2/arcsec2/um

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

Ejemplo n.º 6

Archivo: ejas.py Proyecto: HajimeKawahara/exocounts

import exocounts
import csv
import convmag
from astropy import constants as const
from astropy import units as u
import numpy as np

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