Beispiel #1
0
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
Beispiel #7
0
'''
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
Beispiel #8
0
'''
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)
Beispiel #9
0
'''
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
Beispiel #10
0
'''
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)