metrics = [ fulloutput[0]['throughput'], fulloutput[0]['noise'], fulloutput[0]['sensitivity (Student)'] ] metrics = np.array(metrics) return metrics, fulloutput[3] sp.save_obs = False sp.show_cube = False sp.return_cube = True sp.show_wframe = False ap.companion = False sp.num_processes = 20 iop.date = '180627' #'180407/' iop.update(iop.date) theta = 45 lod = 8 # ************************ H2RG **************************** tp.detector = 'H2RG' tp.NCPA_type = 'Static' #'Wave' # ++++++++ H2RG - Ref_PSF ++++++++ iop.hyperFile = iop.datadir + '/nocomp_BinH2RG_noCoron_hyper.pkl' tp.occulter_type = 'None' num_exp = 10 ap.exposure_time = 0.01 #0.001 ap.numframes = int(num_exp * ap.exposure_time / cp.frame_time) tp.nwsamp = 1 hp.use_readnoise = False
# tp.CPA_type = 'Static' # tp.aber_params['OOPP'] = [8,4] tp.aber_params = { 'CPA': True, 'NCPA': True, 'QuasiStatic': False, # or Static 'Phase': True, 'Amp': False, 'n_surfs': 4, 'OOPP': [16, 8, 8, 4] } #False}# mp.date = '180416mkids/' cp.date = '1804171hr8m/' import os cp.atmosdir = os.path.join(cp.rootdir, cp.data, cp.date) iop.update(mp.date) sp.num_processes = 48 # tp.occulter_type = '8th_Order' tp.occulter_type = 'Vortex' # num_exp = 2000#500#1000#50#50#1000 # ap.exposure_time = 0.001 # 0.001 # cp.frame_time = 0.001 num_exp = 300 #500#1000#50#50#1000 ap.exposure_time = 0.01 # 0.001 cp.frame_time = 0.01 ap.numframes = int(num_exp * ap.exposure_time / cp.frame_time) # ap.startframe=ap.numframes tp.piston_error = True xlocs = range(0, 128) # range(0,128)#65 ylocs = range(0, 128) # range(0,128)#85 tp.band = np.array([800, 1500])
import medis.Detector.get_photon_data as gpd import Examples.SSD_example as SSD import Analysis.stats as stats import traceback import itertools # os.system("taskset -p 0xfffff %d" % os.getpid()) ap.star_photons = int(1e5) ap.companion = True ap.contrast = [0.01] #[0.1,0.1] ap.lods = [[-2.5, 2.5]] ap.numframes = 500 tp.detector = 'MKIDs' sp.save_obs = True mp.date = '180406b/' iop.update('180406b/') sp.show_cube = False sp.return_cube = False cp.vary_r0 = True tp.use_atmos = True tp.ao_act = 29 #odd numbers appear to remove the cross pattern tp.use_ao = True # True # tp.occulter_type = 'None' # tp.occulter_type = '8TH_ORDER' # tp.NCPA_type = 'Static' tp.CPA_type = 'Static' tp.nwsamp = 1 tp.satelite_speck = False tp.speck_peakIs = [0.05] #[0.01, 0.0025] tp.speck_phases = [np.pi / 2] #[np.pi / 2.,np.pi / 2.] tp.speck_locs = [[64, 30]] #[[50, 50], [100,100]]