Skip to content

normalize spectrum #22

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
vanfossen merged 4 commits into from
Jan 15, 2023
Merged

normalize spectrum #22

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
9275655
add test notebook
vanfossen Jan 13, 2023
494b07e
update SerialSlabs(), process_spectrum return
vanfossen Jan 13, 2023
52db04e
normalize spectrum
vanfossen Jan 15, 2023
932374a
remove test notebook
vanfossen Jan 15, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
14 changes: 10 additions & 4 deletions app.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,11 +49,14 @@ def spectrum():
# --> detector response spectrum
processed_spectrum = __process_spectrum(params, spectrum, True)

# https://radis.readthedocs.io/en/latest/source/radis.spectrum.spectrum.html#radis.spectrum.spectrum.Spectrum.get
x_value, y_value = processed_spectrum.get("transmittance_noslit")

# convert dictionary values to strings and return as JSON
return {
"success": True,
"x": list(processed_spectrum.keys()),
"y": [str(flt) for flt in processed_spectrum.values()],
"x": list(x_value),
"y": list(map(str, y_value)),
}


Expand Down Expand Up @@ -88,11 +91,14 @@ def background():
# --> detector response spectrum
processed_spectrum = __process_spectrum(data, background_spectrum, True)

# https://radis.readthedocs.io/en/latest/source/radis.spectrum.spectrum.html#radis.spectrum.spectrum.Spectrum.get
x_value, y_value = processed_spectrum.get("transmittance_noslit")

# convert dictionary values to strings and return as JSON
return {
"success": True,
"x": list(processed_spectrum.keys()),
"y": [str(flt) for flt in processed_spectrum.values()],
"x": list(x_value),
"y": list(map(str, y_value)),
}


Expand Down
41 changes: 22 additions & 19 deletions functions.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -426,42 +426,48 @@ def __process_spectrum(params, raw_spectrum, find_peaks):
name="InSb",
)

# SerialSlabs() multiplies the transmittance values (y-values) of two provided spectra
# https://radis.readthedocs.io/en/latest/source/radis.los.slabs.html#radis.los.slabs.SerialSlabs
# list of spectra to multiply
slabs = []

# ----- a.) transmission spectrum of gas sample -----
slabs.append(raw_spectrum)

# ----- b.) blackbody spectrum of source -----
spectrum = SerialSlabs(raw_spectrum, spec_sPlanck)
slabs.append(spec_sPlanck)

# ----- c.) transmission spectrum of windows/beamsplitter -----
# ----- c.1) Beamsplitter -----
match params["beamsplitter"]:
case "AR_ZnSe":
spectrum = SerialSlabs(spectrum, spec_AR_ZnSe)
slabs.append(spec_AR_ZnSe)
case "AR_CaF2":
spectrum = SerialSlabs(spectrum, spec_AR_CaF2)
slabs.append(spec_AR_CaF2)

# ----- c.2) cell windows -----
match params["cellWindow"]:
case "CaF2":
spectrum = SerialSlabs(spectrum, spec_CaF2)
spectrum = SerialSlabs(spectrum, spec_CaF2)
slabs.extend([spec_CaF2, spec_CaF2])
case "ZnSe":
spectrum = SerialSlabs(spectrum, spec_ZnSe)
spectrum = SerialSlabs(spectrum, spec_ZnSe)
slabs.extend([spec_ZnSe, spec_ZnSe])

# ----- d.) detector response spectrum -----
match params["detector"]:
case "MCT":
spectrum = SerialSlabs(spectrum, spec_ZnSe)
spectrum = SerialSlabs(spectrum, spec_MCT)
slabs.extend([spec_ZnSe, spec_MCT])
case "InSb":
spectrum = SerialSlabs(spectrum, spec_sapphire)
spectrum = SerialSlabs(spectrum, spec_InSb)
slabs.extend([spec_sapphire, spec_InSb])

# SerialSlabs() multiplies the transmittance values (y-values) of the selected spectra
# https://radis.readthedocs.io/en/latest/source/radis.los.slabs.html#radis.los.slabs.SerialSlabs
spectrum = SerialSlabs(*slabs, modify_inputs="True")

spectrum.normalize(normalize_how="mean", inplace=True, force=True)

# add random noise to spectrum
# https://radis.readthedocs.io/en/latest/source/radis.spectrum.operations.html#radis.spectrum.operations.add_array
spectrum = add_array(
spectrum,
sum(np.random.normal(0, 800000000, (params["numScan"], len(w))))
sum(np.random.normal(0, 1, (params["numScan"], len(w))))
/ params["numScan"],
var="transmittance_noslit",
)
Expand All @@ -479,11 +485,8 @@ def __process_spectrum(params, raw_spectrum, find_peaks):
# print()
# print(lines)

# https://radis.readthedocs.io/en/latest/source/radis.spectrum.spectrum.html#radis.spectrum.spectrum.Spectrum.get
x_value, y_value = spectrum.get("transmittance_noslit")

# Return spectrum as a dictionary
return dict(zip(x_value, y_value))
# return processed spectrum
return spectrum


def __generate_spectrum(params):
Expand Down