Note
Go to the end to download the full example code.
Cylinder: Qsca vs wavelength std#
Importing the package dependencies: numpy, PyMieSim
import numpy as np
from TypedUnit import ureg
from PyMieSim.experiment.scatterer import Cylinder
from PyMieSim.experiment.source import Gaussian
from PyMieSim.experiment import Setup
from PyOptik import Material
source = Gaussian(
wavelength=np.linspace(200, 1800, 300) * ureg.nanometer,
polarization=0 * ureg.degree,
optical_power=1 * ureg.watt,
NA=0.2 * ureg.AU,
)
scatterer = Cylinder(
diameter=np.linspace(400, 1400, 10) * ureg.nanometer,
property=Material.silver,
medium_property=1 * ureg.RIU,
source=source,
)
experiment = Setup(scatterer=scatterer, source=source)
dataframe = experiment.get("Qsca")
dataframe.plot(x="source:wavelength", std="scatterer:diameter")

<Figure size 800x500 with 1 Axes>
Total running time of the script: (0 minutes 0.510 seconds)