Python extract_wifes_stdstar Examples

Example #1

File: reduce_blue_data.py Project: arizzuto/pipeline

def run_extract_stars(metadata, prev_suffix, curr_suffix, type='all',**args):
    # for each std, extract spectrum as desired
    std_obs_list = get_primary_std_obs_list(metadata, type=type)
    #print std_obs_list
    for fn in std_obs_list:
        in_fn  = '%s%s.p%s.fits' % (out_dir, fn, prev_suffix)
        out_fn = '%s%s.x%s.dat'  % (out_dir, fn, prev_suffix)
        print 'Extract %s standard star from %s' % (type, in_fn.split('/')[-1])
        pywifes.extract_wifes_stdstar(in_fn,
                                      save_fn=out_fn,
                                      save_mode='ascii',
                                      **args)
    return

Example #2

def run_extract_stars(metadata, prev_suffix, curr_suffix, type='all',**args):
    # for each std, extract spectrum as desired
    std_obs_list = get_primary_std_obs_list(metadata, type=type)
    #print std_obs_list
    for fn in std_obs_list:
        in_fn  = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix))
        out_fn = os.path.join(out_dir, '%s.x%s.dat'  % (fn, prev_suffix))
        print('Extract %s standard star from %s' % (type, in_fn.split('/')[-1]))
        pywifes.extract_wifes_stdstar(in_fn,
                                      save_fn=out_fn,
                                      save_mode='ascii',
                                      **args)
    return

Example #3

File: reduce_blue_data.py Project: waldo2590/pipeline

            ny_orig=76, offset_orig=2.0, **args)
        #print squirrel
    return

#------------------------------------------------------
# Standard star extraction
def run_extract_stars(metadata, prev_suffix, curr_suffix, type='all',**args):
    # for each std, extract spectrum as desired
    std_obs_list = get_primary_std_obs_list(metadata, type=type)
    #print std_obs_list
    for fn in std_obs_list:
        in_fn  = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix))
        out_fn = os.path.join(out_dir, '%s.x%s.dat'  % (fn, prev_suffix))
        print('Extract %s standard star from %s' % (type, in_fn.split('/')[-1]))
        pywifes.extract_wifes_stdstar(in_fn,
                                      save_fn=out_fn,
                                      save_mode='ascii',
                                      **args)
    return

# Sensitivity Function fit
def run_derive_calib(metadata, prev_suffix, curr_suffix, method = 'poly',**args):
    std_obs_list = get_primary_std_obs_list(metadata, type='flux')
    std_cube_list = [os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix))
                     for fn in std_obs_list]
    extract_list = [os.path.join(out_dir, '%s.x%s.dat' % (fn, prev_suffix))
                    for fn in std_obs_list]
    print('Deriving sensitivity function')
    best_calib = pywifes.derive_wifes_calibration(
        std_cube_list,
        calib_fn,
        extract_in_list=extract_list,method=method,