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))