cHbeta = catalogue_df.loc[objName, cHbeta_type]
flux_dered = dz.derreddening_spectrum(wave,
flux,
reddening_curve=red_curve,
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'])
