Skip to content
Snippets Groups Projects
Commit a0d92884 authored by Udo Eisenbarth's avatar Udo Eisenbarth :speech_balloon:
Browse files

Improve spectrum example with sub function

parent 8dbb5d97
No related branches found
No related tags found
No related merge requests found
use opossum::error::OpossumError; use opossum::error::OpossumError;
use opossum::spectrum::Spectrum; use opossum::spectrum::Spectrum;
use uom::si::{f64::Length, length::{nanometer,meter}}; use uom::si::{f64::Length, length::meter};
fn main() -> Result<(), OpossumError> { fn main() -> Result<(), OpossumError> {
let mut s = Spectrum::new( let mut s = Spectrum::new(
Length::new::<meter>(400.0)..Length::new::<meter>(450.0), Length::new::<meter>(400.0)..Length::new::<meter>(450.0),
Length::new::<meter>(0.1), Length::new::<meter>(0.1),
)?; )?;
//s.add_single_peak(Length::new::<meter>(405.0), 1.0)?; s.add_lorentzian_peak(Length::new::<meter>(415.0), Length::new::<meter>(3.2), 2.0)?;
s.add_lorentzian_peak(Length::new::<meter>(410.0), Length::new::<meter>(0.2), 1.0)?;
let mut s2=Spectrum::new( let mut s2=Spectrum::new(
Length::new::<meter>(400.0)..Length::new::<meter>(450.0), Length::new::<meter>(400.0)..Length::new::<meter>(450.0),
...@@ -16,6 +15,13 @@ fn main() -> Result<(), OpossumError> { ...@@ -16,6 +15,13 @@ fn main() -> Result<(), OpossumError> {
)?; )?;
s2.add_lorentzian_peak(Length::new::<meter>(430.0), Length::new::<meter>(1.2), 0.5)?; s2.add_lorentzian_peak(Length::new::<meter>(430.0), Length::new::<meter>(1.2), 0.5)?;
s.add(&s2); s.add(&s2);
let mut s3=Spectrum::new(
Length::new::<meter>(400.0)..Length::new::<meter>(450.0),
Length::new::<meter>(0.05),
)?;
s3.add_lorentzian_peak(Length::new::<meter>(420.0), Length::new::<meter>(0.3), 0.02)?;
s.sub(&s3);
s.to_plot("spectrum.svg"); s.to_plot("spectrum.svg");
Ok(()) Ok(())
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment