def _init_instrument(self):
"""Set up the instrument and contamination model."""
self.instrument = Instrument('example',
[sdss_g, sdss_r, sdss_i, sdss_z])
self.cm = SMContamination(self.instrument, "i'")
self.lnpriors.append(lambda pv: where(pv[:, 4] < pv[:, 5], 0, -inf))
self.lnpriors.append(lambda pv: where(pv[:, 8] < pv[:, 5], 0, -inf))
def _init_instrument(self):
"""Set up the instrument and contamination model."""
qe = TabulatedFilter('MockQE',
[300, 350, 500, 550, 700, 800, 1000, 1050],
[0.10, 0.20, 0.90, 0.96, 0.90, 0.75, 0.11, 0.05])
self.instrument = Instrument('MockInstrument',
[sdss_g, sdss_r, sdss_i, sdss_z],
(qe, qe, qe, qe))
self.cm = SMContamination(self.instrument, "i'")
self.lnpriors.append(lambda pv: where(pv[:, 4] < pv[:, 5], 0, -inf))
def __init__(self, setup: SimulationSetup, **kwargs):
self.setup = self.s = s = setup
self.t_exposure_d = Qty(kwargs.get('exptime', 60), 's').rescale('d')
self.t_baseline_d = Qty(s.t_baseline, 'h').rescale('d')
self.ldcs = s.ldcs
self.tm = QuadraticModel(klims=(0.01, 0.99), nk=512)
self.filters = "g' r' i' z'".split()
self.npb = len(self.filters)
self.k_apparent, self.p, self.a, self.b, self.i = s.orbital_parameters
self.duration_d = Qty(
duration_eccentric(self.p, self.k_apparent, self.a, self.i, 0, 0,
1), 'd')
# Contamination
# -------------
qe_be = TabulatedFilter('1024B_eXcelon', [
300, 325, 350, 400, 450, 500, 700, 800, 850, 900, 950, 1050, 1150
], [
0.0, 0.1, 0.25, 0.60, 0.85, 0.92, 0.96, 0.85, 0.70, 0.50, 0.30,
0.05, 0.0
])
qe_b = TabulatedFilter(
'2014B', [300, 350, 500, 550, 700, 800, 1000, 1050],
[0.10, 0.20, 0.90, 0.96, 0.90, 0.75, 0.11, 0.05])
qes = qe_be, qe_b, qe_be, qe_be
self.instrument = instrument = Instrument(
'MuSCAT2', (sdss_g, sdss_r, sdss_i, sdss_z), qes)
self.contaminator = SMContamination(instrument, "i'")
self.hteff = setup.hteff
self.cteff = setup.cteff
self.i_contamination = setup.c
self.k_true = setup.k_apparent / sqrt(1 - self.i_contamination)
self.contamination = self.contaminator.contamination(
self.i_contamination, self.hteff, self.cteff)
def _init_instrument(self):
self.instrument = Instrument('MuSCAT2',
[sdss_g, sdss_r, sdss_i, sdss_z])
self.cm = SMContamination(self.instrument, "i'")
def _init_instrument(self):
self.instrument = Instrument('example', [sdss_g, sdss_r, sdss_i])
self.cm = SMContamination(self.instrument, "i'")
self.lnpriors.append(lambda pv: where(pv[:, 4] < pv[:, 5], 0, -inf))