color='black'
) #, y_error=log10(unumpy.std_devs(ratios_dict['blue_y'])))
#--Red points
if len(ratios_dict['red_x']) > 0:
dz.data_plot(
range_red,
log10(
unumpy.nominal_values(
ratios_dict['line_Flux_Red'])),
'Observed fluxes',
markerstyle='o',
color='black'
) #, y_error=log10(unumpy.std_devs(ratios_dict['red_y'])))
E_BV = dz.Ebv_from_cHbeta(cHbeta_in_MagEr, red_curve, R_v)
Xx_blue = dz.reddening_Xx(ratios_dict['line_wave_Blue'], red_curve,
R_v)
Xx_red = dz.reddening_Xx(ratios_dict['line_wave_Red'], red_curve,
R_v)
Blue_int = ratios_dict['line_Flux_Blue'] * unum_pow(
10, 0.4 * Xx_blue * E_BV)
Red_int = ratios_dict['line_Flux_Red'] * unum_pow(
10, 0.4 * Xx_red * E_BV)
label_regression = r'Line intensities {rd_curve}: $c(H\beta) = {Value} \pm {Error}$'.format(
rd_curve=red_curve,
Value=round(cHbeta_all_MagEr.nominal_value, 3),
Error=round(cHbeta_all_MagEr.std_dev, 3))
