#Perform the reddening correction
cHbeta = catalogue_df.iloc[i][cHbeta_type]
Obs_spectrum_dered = dz.derreddening_spectrum(Wave_T,
Int_T,
reddening_curve=red_curve,
cHbeta=cHbeta.nominal_value,
R_v=R_v)
spectrum_dered = Obs_spectrum_dered - Int_N
#Generating the starlight files
FileName = basename(fits_file) + stellar_files_exten
Grid_FileName, Sl_OutputFile, Sl_OutputFolder, X_1Angs, Y_1Angs = dz.GenerateStarlightFiles(
ouput_folder,
FileName,
objName,
catalogue_df.iloc[i],
Wave_T,
spectrum_dered,
ext_loop=stellar_files_exten)
print '--Output file ', Sl_OutputFile,
#Plot the data
dz.data_plot(Wave_T, spectrum_dered, "Reduced spectrum")
dz.data_plot(X_1Angs, Y_1Angs, "Resampled spectrum", linestyle='--')
# Set titles and legend
PlotTitle = 'Object {} Resampled spectrum for starlight'.format(objName)
dz.FigWording(r'Wavelength $(\AA)$',
'Flux' + r'$(erg\,cm^{-2} s^{-1} \AA^{-1})$', PlotTitle)
