cHbeta=cHbeta.nominal_value,
R_v=R_v)
#Define the maximum sigma for the fitting from the Hbeta line
Sigma = ImportDispersionVelocity(lineslog_frame)
UpperDispersionVelocity_Limit = str(round(
Sigma, 1)) + ' [vd_upp (km/s)] = upper allowed vd\n'
dz.replace_line(GridFileAddress, 21, UpperDispersionVelocity_Limit)
#Launch starlight
print '--Initiating starlight for ', fits_name, Sigma
dz.Starlight_Launcher(Grid_FileName, dz.RootFolder)
print '--Starlight finished succesfully ended:', Sl_OutputFile
#Get stellar spectrum from starlight file
Input_Wavelength, Input_Flux, Output_Flux, MaskPixels, ClippedPixels, FlagPixels, Parameters = dz.File_to_data(
Sl_OutputFolder, Sl_OutputFile)
#Export data to fits file
stellar_cont_fits = objName + '_StellarContinuum.fits'
dz.Data_2_Fits(ouput_folder,
stellar_cont_fits,
header_0,
Input_Wavelength,
Output_Flux,
NewKeyWord=['STALIGHT', 'Basic Treatment'])
#Plot the data
dz.data_plot(wave, flux_dered, "obs de-red")
dz.data_plot(Input_Wavelength, Input_Flux, "Input Spectra")
dz.data_plot(Input_Wavelength, Output_Flux, "Stellar absorption")
from dazer_methods import Dazer
#Declare object
dz = Dazer()
dz.FigConf(n_colors=5)
CodeName1 = '70'
CodeName2 = 'SDSS2'
extension = '_WHT.fits'
Folder_1 = '/home/vital/Dropbox/Astrophysics/Data/WHT_Catalogue_SulfurRegression/Objects/' + CodeName1 + '/'
Folder_2 = '/home/vital/Dropbox/Astrophysics/Data/WHT_Catalogue_SulfurRegression/Objects/' + CodeName2 + '/'
Wave1, Flux1, ExtraData1 = dz.File_to_data(Folder_1,
'obj' + CodeName1 + extension)
Wave2, Flux2, ExtraData2 = dz.File_to_data(Folder_2,
'obj' + CodeName2 + extension)
dz.Axis.set_yscale('log')
print Flux1
print Flux2
dz.data_plot(Wave1, Flux1, label=CodeName1, color=dz.ColorVector[2][0])
dz.data_plot(Wave2, Flux2, label=CodeName2, color=dz.ColorVector[2][1])
dz.FigWording(
r'Wavelength $(\AA)$', 'Flux ' + r'$(erg\,cm^{-2} s^{-1} \AA^{-1})$',
'Objects {obj1} and {obj2} comparison'.format(obj1=CodeName1,
obj2=CodeName2))
dz.display_fig()
