def add_flux (gui, graph_obj, label, atmos_site, galactic_file, mirror_temp,
mirror_constant, zodiac_file, cib, cmb):
temp_tot = temp.total(gui, atmos_site.file, galactic_file.file, mirror_temp,
mirror_constant, zodiac_file.file, cib, cmb)
flux_tot = flux_convert(gui, cal.intensity(gui.interp.freq_array, temp_tot))
data_set = new_dataset("Total Flux ("+label+")", gui.energy_form,
"Flux", flux_units(gui), graph_list(gui, flux_tot))
graph_obj.dataset_list.append(data_set)
def add_radiance (gui, graph_obj, site_file, spec_res):
data_name = "Atmos Radiance ("+site_file.name+")"
if gui.noise_what == 0: # BLING
bling_list = graph_list(gui, bling.radiance(gui, site_file.file, spec_res))
data_set = new_dataset(data_name, gui.energy_form,
"BLING", bling_units(gui), bling_list)
elif gui.noise_what == 1: # flux
flux_list = graph_list(gui,
flux_convert(gui, cal.intensity(gui.interp.freq_array,
temp.radiance(gui, site_file.file))))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", flux_units(gui), flux_list)
elif gui.noise_what == 2: # temperature
temp_list = graph_list(gui, temp.radiance(gui, site_file.file))
data_set = new_dataset(data_name, gui.energy_form,
"Temperature", "K", temp_list)
else: # diffraction-limited flux
intensity = cal.intensity(gui.interp.freq_array,
temp.radiance(gui, site_file.file))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", signal_units(gui), flux_dl(gui, intensity, spec_res))
graph_obj.dataset_list.append(data_set)
def add_cmb (gui, graph_obj, spec_res):
data_name = "Cosmic Microwave Bkgd"
if gui.noise_what == 0: # BLING
bling_list = graph_list(gui, bling.cmb(gui, spec_res))
data_set = new_dataset(data_name, gui.energy_form,
"BLING", bling_units(gui), bling_list)
elif gui.noise_what == 1: # flux
flux_list = graph_list(gui,
flux_convert(gui,
cal.intensity(gui.interp.freq_array, temp.cmb(gui))))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", flux_units(gui), flux_list)
elif gui.noise_what == 2: # temperature
temp_list = graph_list(gui, temp.cmb(gui))
data_set = new_dataset(data_name, gui.energy_form,
"Temperature", "K", temp_list)
else: # diffraction limited flux
intensity = cal.intensity(gui.interp.freq_array, temp.cmb(gui))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", signal_units(gui), flux_dl(gui, intensity, spec_res))
graph_obj.dataset_list.append(data_set)
def add_cib (gui, graph_obj, spec_res):
data_name = "Cosmic Infrared Bkgd"
# TODO: convert to equation fit
if gui.noise_what == 0: # BLING
bling_list = graph_list(gui, bling.generic(gui,
"data/Backgrounds/CIB/cib.xlsx", spec_res))
data_set = new_dataset(data_name, gui.energy_form,
"BLING", bling_units(gui), bling_list)
elif gui.noise_what == 1: # flux
flux_list = graph_list(gui,
flux_convert(cal.intensity(gui.interp.freq_array,
temp.generic(gui, "data/Backgrounds/CIB/cib.xlsx"))))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", flux_units(gui), flux_list)
elif gui.noise_what == 2: # temperature
temp_list = graph_list(gui,
temp.generic(gui, "data/Backgrounds/CIB/cib.xlsx"))
data_set = new_dataset(data_name, gui.energy_form,
"Temperature", "K", temp_list)
else: # diffraction limited flux
intensity = cal.intensity(gui.interp.freq_array,
temp.generic(gui, "data/Backgrounds/CIB/cib.xlsx"))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", signal_units(gui), flux_dl(gui, intensity, spec_res))
graph_obj.dataset_list.append(data_set)
def add_mirror (gui, graph_obj, metal_name, mirror_temp, constant, spec_res):
data_name = "Thermal Mirror ("+metal_name+", "+str(mirror_temp)+" K)"
if gui.noise_what == 0: # BLING
bling_list = graph_list(gui, bling.mirror(gui,
mirror_temp, constant, spec_res))
data_set = new_dataset(data_name,
gui.energy_form, "BLING", bling_units(gui), bling_list)
elif gui.noise_what == 1: # flux
flux_list = graph_list(gui,
flux_convert(gui, cal.intensity(gui.interp.freq_array,
temp.mirror(gui, mirror_temp, constant))))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", flux_units(gui), flux_list)
elif gui.noise_what == 2: # temperature
temp_list = graph_list(gui, temp.mirror(gui, mirror_temp, constant))
data_set = new_dataset(data_name,
gui.energy_form, "Temperature", "K", temp_list)
else: # diffraction limited flux
intensity = cal.intensity(gui.interp.freq_array,
temp.mirror(gui, mirror_temp, constant))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", signal_units(gui), flux_dl(gui, intensity, spec_res))
graph_obj.dataset_list.append(data_set)
def add_galactic (gui, graph_obj, galactic_file, spec_res):
data_name = "Galactic Emission ("+galactic_file.name+")"
if gui.noise_what == 0: # BLING
bling_list = graph_list(gui, bling.generic(gui, galactic_file.file, spec_res))
data_set = new_dataset(data_name,
gui.energy_form, "BLING", bling_units(gui), bling_list)
elif gui.noise_what == 1: # flux
flux_list = graph_list(gui,
flux_convert(gui, cal.intensity(gui.interp.freq_array,
temp.generic(gui, galactic_file.file))))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", flux_units(gui), flux_list)
elif gui.noise_what == 2: # temperature
temp_list = graph_list(gui, temp.generic(gui, galactic_file.file))
data_set = new_dataset(data_name,
gui.energy_form, "Temperature", "K", temp_list)
else:
intensity = cal.intensity(gui.interp.freq_array,
temp.generic(gui, galactic_file.file))
data_set = new_dataset(data_name, gui.energy_form,
"Flux", signal_units(gui), flux_dl(gui, intensity, spec_res))
graph_obj.dataset_list.append(data_set)
