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