Пример #1
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)
        dz.Axis.set_xlim(3500, 5250)

        #Save plot
        output_pickle = '{objFolder}{stepCode}_{objCode}_{ext}'.format(
            objFolder=ouput_folder,
            stepCode=dz.ScriptCode,
            objCode=objName,
            ext='StellarContinuum')
        dz.save_manager(output_pickle, save_pickle=True)

    except:
        dz.log_error(objName)

print 'All data treated', dz.display_errors()
        dz.data_plot(x_trendline, y_trendline, "Balmer jump regression", linestyle="--")

        # Format the graphs
        plotTitle = r"Object {} Nebular continuum calculation".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_Comparison"
        )
        dz.save_manager(output_pickle, save_pickle=True)

    except:
        dz.log_error(objName)

print "All data treated", dz.display_errors()

# from dazer_methods import Dazer
# from libraries.Astro_Libraries.Nebular_Continuum import NebularContinuumCalculator
#
# def nebular_continua_calculation_plots():
#
#     #Plotting the data
#     dz.data_plot(Wave,         Int,                'Reduced spectrum',                         dz.ColorVector[2][0])
#     dz.data_plot(NebWave,      NebularInt_Hbeta,    'Nebular continuum - Hbeta calibration',    dz.ColorVector[2][1])
#     dz.data_plot(NebWave,      NebularInt_Halpha,   'Nebular continuum - Halfa calibration',    dz.ColorVector[2][2])
#     dz.data_plot(x_trendline,  y_trendline,         'Balmer jump regression',                   dz.ColorVector[2][3],     linestyle='--')
#
#     #Format the graphs