#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')