def hierarchical_tile(masterfile,tilefile): """ Create Hierarchical tile from Master prior :param masterfile: Master prior file :param tilefile: File containing Tiling scheme """ try: taskid = np.int(os.environ['SGE_TASK_ID']) task_first=np.int(os.environ['SGE_TASK_FIRST']) task_last=np.int(os.environ['SGE_TASK_LAST']) except KeyError: print("Error: could not read SGE_TASK_ID from environment") taskid = int(input("Please enter task id: ")) print("you entered", taskid) with open(tilefile, 'rb') as f: obj = pickle.load(f) tiles = obj['tiles'] order = obj['order'] tiles_large = obj['tiles_large'] order_large = obj['order_large'] obj=xidplus.io.pickle_load(masterfile) priors = obj['priors'] moc = moc_routines.get_fitting_region(order_large, tiles_large[taskid - 1]) for p in priors: p.moc = moc p.cut_down_prior() outfile = 'Tile_'+ str(tiles_large[taskid - 1]) + '_' + str(order_large) + '.pkl' with open(outfile, 'wb') as f: pickle.dump({'priors':priors, 'version':xidplus.io.git_version()}, f)
def hierarchical_tile(masterfile, tilefile): """ Create Hierarchical tile from Master prior :param masterfile: Master prior file :param tilefile: File containing Tiling scheme """ try: taskid = np.int(os.environ['SGE_TASK_ID']) task_first = np.int(os.environ['SGE_TASK_FIRST']) task_last = np.int(os.environ['SGE_TASK_LAST']) except KeyError: print("Error: could not read SGE_TASK_ID from environment") taskid = int(input("Please enter task id: ")) print("you entered", taskid) with open(tilefile, 'rb') as f: obj = pickle.load(f) tiles = obj['tiles'] order = obj['order'] tiles_large = obj['tiles_large'] order_large = obj['order_large'] obj = xidplus.io.pickle_load(masterfile) priors = obj['priors'] moc = moc_routines.get_fitting_region(order_large, tiles_large[taskid - 1]) for p in priors: p.moc = moc p.cut_down_prior() outfile = 'Tile_' + str( tiles_large[taskid - 1]) + '_' + str(order_large) + '.pkl' with open(outfile, 'wb') as f: pickle.dump({'priors': priors, 'version': xidplus.io.git_version()}, f)
############################# #outfile=output_folder+'failed_tiles.pkl' #with open(outfile, 'rb') as f: # obj=pickle.load(f) tiles=obj['tiles'] order=obj['order'] order_large=obj['order_large'] tile_large=moc_routines.tile_in_tile(order,tiles[taskid-1],order_large) outfile=output_folder+'Tile_'+str(tile_large)+'_'+str(order_large)+'.pkl' with open(outfile, 'rb') as f: obj=pickle.load(f) priors=obj['priors'] moc=moc_routines.get_fitting_region(order,tiles[taskid-1]) priors[0].moc=moc priors[0].cut_down_prior() print('fitting '+ str(priors[0].nsrc)+' sources \n') print('there are '+ str(priors[0].snpix)+' pixels') print('fitting tile:'+str(tiles[taskid-1])+' order: '+str(order)+' of '+str(len(tiles))+' tiles') from xidplus.stan_fit import MIPS priors[0].prior_bkg(0.0,1) priors[0].get_pointing_matrix() #priors[0].upper_lim_map() #priors[0].prior_flux_upper=(priors[0].prior_flux_upper-10.0+0.02)/np.max(priors[0].prf) fit=MIPS.MIPS_24(priors[0],iter=1000)