from CodeTools.PlottingManager import myPickle from Plotting_Libraries.dazer_plotter import Plot_Conf from numpy import array, power, savetxt, transpose, unique, where, zeros, ones from Math_Libraries import sigfig from Scientific_Lib.IrafMethods import Pyraf_Workflow from collections import OrderedDict from os import mkdir #Declare Classes pv = myPickle() dz = Plot_Conf() py_w = Pyraf_Workflow('WHT') # #-----------------------------------------STARBURST EQUIVALENT WIDTH EVOLUTION---------------------------------------- FilesFolder = '/home/vital/Dropbox/Astrophysics/Data/Starburst_Spectra_z0.004/' FilesPattern = '_txt_LinesLog_v3.txt' #Locate files on hard drive FilesList = pv.Folder_Explorer(FilesPattern, FilesFolder, CheckComputer=False) # #Define figure format dz.FigConf() #Lines to plot # H_Lines = ['H1_3970A','H1_4102A','H1_4340A', 'H1_4861A', 'H1_6563A'] H_Lines = [ 'He1_3188A', 'He1_4026A', 'He1_4471A', 'He2_4686A', 'He1_5016A', 'He1_5876A', 'He1_6678A' ] #Define dictionary to store the data
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from PipeLineMethods.ManageFlow import DataToTreat from Scientific_Lib.IrafMethods import Pyraf_Workflow pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') # input_spec = 'obj70_s_c_cr_f_a_bg0001.fits[0]' # output_spec = 'obj70_s_c_cr_f_a_t_bg_fxW_TESTPYRAF.fits' # Folder = '/home/vital/Desktop/SDSS_Galaxies/Flux_Calibration_Steps/Night1_Blue/' # sens_func = 'a_std_wideCombAll_se.fits[0]' Pattern = 'obj70_s_c_cr_f_a_bg0001.fits' CatalogueFolder = '/home/vital/Desktop/SDSS_Galaxies/Flux_Calibration_Steps/Night1_Blue/' #Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) #Generate plot frame and colors pv.FigFormat_One(ColorConf='Night1') #Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j])
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from os.path import isfile #WHILE YOU USE THE SPLOT TASK USE THE X KEY TO ABSORPTION LINES FROM THE STAR #USE KEY "i" TO SAVE SPECTRA #SAVE THE SPECTRA AS std + "STANDARDCODE" + '_clean.fits' pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') CatalogueFolder, PropertiesDict = py_w.DataFormat(Reduction_Operation = 'Std_Clean') Pattern = PropertiesDict['Suffix_FluxCalibrated_Trimmed'] #Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) #Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): #WARNING NOT THE BEST APPROACH... CLEAN: if ('std' in FilesList[i][j]) and ('Red' in FilesList[i][j]): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j]) StarCode = py_w.StarnameFromFileName(FileName)
""" Created on May 18, 2015 @author: vital """ import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from CodeTools.File_Managing_Tools import File_Manager # Define the library pv = plotMan.myPickle() py_w = Pyraf_Workflow("WHT") Fm = File_Manager() # Establish the location of files CatalogueFolder, PropertiesDict = py_w.DataFormat(Reduction_Operation="Obs_TrimCalibrated") Pattern = PropertiesDict["Suffix_FluxCalibrated"] # THIS SHOULR GO INSIDE THE CLASS Catalogue_Instructions = CatalogueFolder + py_w.SpectraTreating + Fm.Extensions_dictionary["Reduction Instructions"] Fm.select_Table(Catalogue_Instructions, 1, loadheaders_check=True) Objects, Trimming_wavelengths = Fm.get_ColumnData( ["CodeName", "Wavelength_Trim"], 1, StringIndexes=True, datatype=str, unpack_check=True ) # Find and organize files from terminal command or .py file FilesList = Fm.Folder_Explorer(Pattern, CatalogueFolder, CheckComputer=False) # Generate plot frame and colors pv.FigFormat_One(ColorConf="Night1")
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from numpy import median, loadtxt from Scientific_Lib.AstroMethods import Fnu_to_Flam, mAb_to_Fnu pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') Pattern, CatalogueFolder = py_w.DataFormat(Reduction_Operation = 'Std_FluxCalibration') #Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) #Generate plot frame and colors pv.FigFormat_One(ColorConf='Night1') #Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j]) #From the fits file determine which is the star being treated StarName = py_w.StarnameFromFileName(FileName) #Star sensitivity curve SensitivityCurve = 'a_sen_' + StarName
""" Created on May 18, 2015 @author: vital """ import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from numpy import median pv = plotMan.myPickle() py_w = Pyraf_Workflow("WHT") Pattern, CatalogueFolder = py_w.DataFormat(Reduction_Operation="Obs_FluxCalibration") # Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) # Generate plot frame and colors pv.FigFormat_One(ColorConf="Night1") # Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j]) CodeName = FileName[3 : FileName.find("_")] # From the fits file determine which is the star being treated StarName = py_w.StarnameFromFileName(FileName) Sensitivity_CurveCombined = py_w.GlobalSensitiviyCurve # Get Airmass and Exptime for the object
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from os.path import isfile from CodeTools.File_Managing_Tools import File_Manager #Define the library pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') Fm = File_Manager() CatalogueFolder, PropertiesDict = py_w.DataFormat(Reduction_Operation = 'dopcor') PatternBlue = PropertiesDict['Suffix_FluxCalibrated_Trimmed'] PatternRed = PropertiesDict['Suffix_FluxCalibrated_Trimmed_Telluric'] #Find and organize files from terminal command or .py file FilesList_Blue = Fm.Folder_Explorer(PatternBlue, CatalogueFolder + 'Blue/', CheckComputer=False) FilesList_Red = Fm.Folder_Explorer(PatternRed, CatalogueFolder + 'Red/', CheckComputer=False) FilesList = FilesList_Blue + FilesList_Red #THIS SHOULR GO INSIDE THE CLASS Catalogue_Instructions = CatalogueFolder + py_w.SpectraTreating + Fm.Extensions_dictionary['Reduction Instructions'] Fm.select_Table(Catalogue_Instructions, 1, loadheaders_check = True) Objects, Redshifts = Fm.get_ColumnData(['CodeName', 'Dopcor'], 1, StringIndexes = True, datatype = str, unpack_check = True) #Generate plot frame and colors
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from numpy import loadtxt pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') Pattern, CatalogueFolder = py_w.DataFormat() #Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) #Generate plot frame and colors pv.FigFormat_One(ColorConf='Night1') #Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j]) #Get Star name StarName = py_w.StarnameFromFileName(FileName) Star_StandardFile = 'a_std_' + StarName.lower() #Open standard task file Wavelength, Flux_nu_I = loadtxt(FileFolder + Star_StandardFile, skiprows=1, usecols = (0, 1), unpack = True)
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow from CodeTools.File_Managing_Tools import File_Manager #USE A SMALL ORDER 5? pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') Fm = File_Manager() CatalogueFolder, PropertiesDict = py_w.DataFormat(Reduction_Operation = 'Std_Clean') Pattern = '_fxW_t.fits' #THIS SHOULR GO INSIDE THE CLASS Catalogue_Instructions = CatalogueFolder + py_w.SpectraTreating + Fm.Extensions_dictionary['Reduction Instructions'] Fm.select_Table(Catalogue_Instructions, 1, loadheaders_check = True) Objects, TelluricStar_List = Fm.get_ColumnData(['CodeName', 'Telluric_Star'], 1, StringIndexes = True, datatype = str, unpack_check = True) #Find and organize files from terminal command or .py file FilesList = Fm.Folder_Explorer(Pattern, CatalogueFolder, CheckComputer = False) #Generate plot frame and colors pv.FigFormat_One(ColorConf='Night1') #Loop through files
''' Created on May 18, 2015 @author: vital ''' import CodeTools.PlottingManager as plotMan from Scientific_Lib.IrafMethods import Pyraf_Workflow pv = plotMan.myPickle() py_w = Pyraf_Workflow('WHT') Pattern, CatalogueFolder = py_w.DataFormat(Reduction_Operation = 'Std_FluxCalibration') #Find and organize files from terminal command or .py file FilesList = pv.FindAndOrganize(Pattern, CatalogueFolder, CheckComputer=False) GlobalStandard_File = py_w.Airmass_Key #Loop through files for i in range(len(FilesList)): for j in range(len(FilesList[i])): CodeName, FileName, FileFolder = pv.FileAnalyzer(FilesList[i][j]) #Red Airmass_Value = py_w.get_HeaderValue(py_w.Airmass_Key, FileFolder + FileName, Digitalize=True) ExpTime_Value = py_w.get_HeaderValue(py_w.Exptime_Key, FileFolder + FileName, Digitalize=True) #Run the task OutputFile_StandardFile = py_w.StandardTask(FileName, OutputFile = None, FitsFolder =FileFolder, airmass_value = Airmass_Value, exptime_value = ExpTime_Value)