#Plotting the data
dz.data_plot(wave, spectrum_dered, 'Reduced spectrum (without reddening)')
dz.data_plot(wave, NebularInt_Hbeta, 'Nebular flux')
dz.data_plot(wave, Int_dedNeb, 'Removed Nebular contribution')
#Format the graphs
PlotTitle = r'Object {} Nebular continuum substraction'.format(objName)
dz.FigWording(r'Wavelength $(\AA)$',
'Flux' + r'$(erg\,cm^{-2} s^{-1} \AA^{-1})$', PlotTitle)
mean_flux = spectrum_dered.mean()
dz.Axis.set_ylim(-0.05 * mean_flux, 15 * mean_flux)
dz.Axis.set_xlim(3500, 5250)
output_pickle = '{objFolder}{stepCode}_{objCode}_{ext}'.format(
objFolder=ouput_folder,
stepCode=dz.ScriptCode,
objCode=objName,
ext='NebularContinuum_substraction')
dz.save_manager(output_pickle, save_pickle=True)
#Export nebular continuum
dz.Data_2_Fits(ouput_folder,
objName + nebular_fits_exten,
header_0,
wave,
NebularInt_Hbeta,
NewKeyWord=['NEBUSPEC', 'zanstra_hbeta'])
print '\nAll data treated\n', dz.display_errors()
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")
#Set titles and legend
PlotTitle = 'Object ' + objName + ' emission and stellar and spectra'
dz.FigWording(r'Wavelength $(\AA)$',
'Flux' + r'$(erg\,cm^{-2} s^{-1} \AA^{-1})$', PlotTitle)
mean_flux = Input_Flux.mean()
dz.Axis.set_ylim(-0.05 * mean_flux, 15 * mean_flux)
'Flux' + r'$(erg\,cm^{-2} s^{-1} \AA^{-1})$', PlotTitle)
#Save data
output_pickle = '{objFolder}{stepCode}_{objCode}_{ext}'.format(
objFolder=ouput_folder,
stepCode=dz.ScriptCode,
objCode=objName,
ext='StellarContinuum_substraction')
dz.save_manager(output_pickle, save_pickle=True)
#Export fits
Int_E_redd = dz.reddening_spectrum(Wave_T,
Int_E,
reddening_curve=red_curve,
cHbeta=cHbeta.nominal_value,
R_v=R_v)
dz.Data_2_Fits(ouput_folder,
objName + emitting_ext,
ExtraData_T,
Wave_T,
Int_E_redd,
NewKeyWord=['EMISSPEC', 'only emission lines'])
catalogue_df.loc[objName,
'emission_fits'] = ouput_folder + objName + emitting_ext
#Store the dataframe
dz.save_excel_DF(
catalogue_df,
'/home/vital/Dropbox/Astrophysics/Data/WHT_observations/WHT_Galaxies_properties.xlsx',
df_sheet_format='catalogue_data')
