def check_sensitivity():
    fileList = []
    for arg in sys.argv[1:]:
        if 1:
            fileList.append(arg)

    badregions = Search_slow.badregions
    interp_regions = Search_slow.interp_regions
    trimsize = Search_slow.trimsize
    prim_vsini = StarData.get_vsini(fileList)

    # Remove anything without a vsini
    new_file_list = []
    new_prim_vsini = []
    for vsini, fname in zip(prim_vsini, fileList):
        if vsini is not None:
            new_file_list.append(fname)
            new_prim_vsini.append(vsini)

    Sensitivity.Analyze(new_file_list, new_prim_vsini,
                        hdf5_file='/media/ExtraSpace/PhoenixGrid/IGRINS_Grid.hdf5',
                        extensions=True,
                        resolution=None,
                        trimsize=trimsize,
                        badregions=badregions, interp_regions=interp_regions,
                        metal_values=(0.0,),
                        vsini_values=(75, 100, 125, 150, 175, 200, 225, 250),
                        Tvalues=range(3000, 3300, 100),
                        debug=False,
                        addmode='all',
                        output_mode='hdf5',
                        output_file='Sensitivity.hdf5')
def check_sensitivity():
    fileList = []
    for arg in sys.argv[1:]:
        if 1:
            fileList.append(arg)

    badregions = Search_slow.badregions
    interp_regions = Search_slow.interp_regions
    trimsize = Search_slow.trimsize
    prim_vsini = StarData.get_vsini(fileList)

    Sensitivity.Analyze(
        fileList,
        prim_vsini,
        hdf5_file="/media/ExtraSpace/PhoenixGrid/TS23_Grid.hdf5",
        extensions=True,
        resolution=None,
        trimsize=trimsize,
        badregions=badregions,
        interp_regions=interp_regions,
        metal_values=(0.0,),
        vsini_values=(100, 150, 200, 250),
        Tvalues=range(7000, 12100, 1000),
        debug=False,
        addmode="all",
        output_mode="hdf5",
        output_file="Sensitivity.hdf5",
    )
def check_sensitivity():
    fileList = []
    for arg in sys.argv[1:]:
        if 1:
            fileList.append(arg)

    badregions = Search_slow.badregions
    interp_regions = Search_slow.interp_regions
    trimsize = Search_slow.trimsize
    prim_vsini = StarData.get_vsini(fileList, vsini_filename='../Useful_Datafiles/Vsini.csv')

    Sensitivity.Analyze(fileList, prim_vsini,
                        hdf5_file='/media/ExtraSpace/PhoenixGrid/CHIRON_Grid.hdf5',
			extensions=True,
                        resolution=None,
                        trimsize=trimsize,
                        badregions=badregions, interp_regions=interp_regions,
                        metal_values=(0.0,),
                        #vsini_values=(0, 10, 20, 30, 40, 50),
                        #Tvalues=range(4000, 6000, 100),
                        vsini_values=(40,),
                        Tvalues=(3600,),
                        debug=False,
                        addmode='all',
                        output_mode='hdf5')
def check_sensitivity():
    fileList = []
    for arg in sys.argv[1:]:
        if 1:
            fileList.append(arg)

    badregions = Search_slow.badregions
    interp_regions = Search_slow.interp_regions
    trimsize = Search_slow.trimsize
    prim_vsini = StarData.get_vsini(fileList)

    Sensitivity.Analyze(fileList, prim_vsini,
                        hdf5_file=hdf5_filename,
                        extensions=True,
                        resolution=None,
                        trimsize=trimsize,
                        badregions=badregions, interp_regions=interp_regions,
                        metal_values=(0.0,),
                        vsini_values=(0, 10, 20, 30, 40, 50),
                        Tvalues=range(5200, 7000, 100),
                        debug=False,
                        addmode='all',
                        output_mode='hdf5')
Ejemplo n.º 5
0
    modeldir = "/media/FreeAgent_Drive/SyntheticSpectra/Sorted/Stellar/Vband/"
    hdf5_filename = "/media/ExtraSpace/PhoenixGrid/HRS_Grid.hdf5"
else:
    modeldir = raw_input("sys.platform not recognized. Please enter model directory below: ")
    if not modeldir.endswith("/"):
        modeldir = modeldir + "/"

if __name__ == "__main__":
    # Parse command line arguments:
    fileList = []
    for arg in sys.argv[1:]:
        if 1:
            fileList.append(arg)

    # Get the primary star vsini values
    prim_vsini = StarData.get_vsini(fileList)

    GenericSearch.slow_companion_search(
        fileList,
        prim_vsini,
        hdf5_file=hdf5_filename,
        extensions=True,
        resolution=None,
        trimsize=trimsize,
        modeldir=modeldir,
        badregions=badregions,
        metal_values=(0.0, -0.5, 0.5),
        vsini_values=(1, 5.0, 10.0, 20.0, 30.0),
        Tvalues=range(3000, 9000, 100),
        observatory="McDonald",
        debug=False,