#-- correct the nr of independent frequencies for the frequency range # that is tested, but only if it is requested if correct_for_frange: nyqstat = kwargs.pop('nyqstat',np.min) nyquist = getNyquist(times,nyqstat=nyqstat) ni = int(freqs.ptp()/nyquist*ni) #p_value = 1. - (1.- (1-2*peak_value/nr_obs)**(nr_obs/2))**ni p_value = 1. - (1.- np.exp(-peak_value))**ni return p_value if __name__=="__main__": import sys import pylab as pl from ivs.aux import loggers logger = loggers.get_basic_logger() #-- run tests if '--test' in sys.argv[1] or '-t' in sys.argv[1]: import doctest doctest.testmod() pl.show() #-- command line interface else: method,args,kwargs = argkwargparser.parse() print "Running method %s with arguments %s and keyword arguments %s"%(method,args,kwargs) if '--help' in args or 'help' in args or 'help' in kwargs: sys.exit() times,signal = ascii.read2array(kwargs.pop('infile')).T[:2] freqs,ampls = globals()[method](times,signal, **kwargs)
hdulist[1] = table print("Axis:") for name in hdulist[1].columns.names: if name.islower() and not name == 'labs': ax = np.unique(hdulist[1].data.field(name)) print(name, len(ax), min(ax), max(ax)) teffs = hdulist[1].data.field('teff') loggs = hdulist[1].data.field('logg') ebvs = hdulist[1].data.field('ebv') #for logg in np.unique(loggs): #keep = loggs==logg #plt.figure() #plt.title(logg) #plt.plot(teffs[keep],ebvs[keep],'ko',ms=2) keep = hdulist[1].data.field('teff') > 0 logger.info('Removing {}/{} false entries'.format(sum(-keep), len(keep))) #print len(ff[1].data),sum(keep) hdulist[1].data = hdulist[1].data[keep] hdulist.close() #} if __name__ == "__main__": logger = loggers.get_basic_logger() imethod, iargs, ikwargs = argkwargparser.parse() output = globals()[imethod](*iargs, **ikwargs)
out, err = c.communicate() returncode = -1 else: # in the case the process completed normally out, err = c.communicate() returncode = c.poll() return out, err, returncode #} if __name__=="__main__": from ivs.aux import loggers logger = loggers.get_basic_logger(clevel='debug') filename = '/home/jorisv/sdB/Uli/hermesvr/reduced/2451576_HRF_OBJ_ext_CosmicsRemoved.fits' wvl_file = '/home/jorisv/sdB/Uli/hermesvr/reduced/2451577' unseq,output, error = run_hermesVR(filename, wvl_file=wvl_file, verbose=False, version='trunk', timeout=500) data = read_hermesVR_velocities(unseq=[unseq]) print data['vrad'], ' +- ', data['vraderr'] #import time #import sys #import doctest #import shutil #import pylab as pl