Added header comments. Renamed phd-of-rho-and-z-two-layer-compute-fluence.py to phd-of-rho-and-z-compute-fluence.py because homogeneous solution not two layer.

Author: hayakawa16

forward-solvers/fluence-of-rho-and-z-two-layer.py

@@ -1,5 +1,9 @@
 # This is an example of python code using VTS to Compute fluence for a two-layer 
-# medium as a function of radial extent and depth at a given set of optical properties 
+# medium as a function of radial extent and depth at a given set of optical properties.
+# A two-layer SDA forward solver is used to compute the fluence with optical properties
+# defined in opRegions[0] opRegions[1], and top layer thickness defined in topLayerThickness [mm].
+# The fluence as a function of rho and z is determined and when displayed, it is mirrored
+# to show full fluence.
 #
 # Import PythonNet
 from pythonnet import load
@@ -92,4 +96,5 @@ def heatmap(values, x, y, x_label="", y_label="", title=""):
     return fig
 
 fluenceChart = heatmap(allFluenceRowsToPlot.tolist(), allRhos.tolist(), list(zs), "ρ [mm]", "z [mm]", "log(Φ(ρ, z) [mm-2])")
-fluenceChart.show(renderer="browser")
+fluenceChart.add_hline(y=topLayerThickness, line_dash="dash", line_color="white", line_width=2)
+fluenceChart.show(renderer="browser")

…f-rho-and-z-two-layer-compute-fluence.py …vers/phd-of-rho-and-z-compute-fluence.pyforward-solvers/phd-of-rho-and-z-two-layer-compute-fluence.py renamed to forward-solvers/phd-of-rho-and-z-compute-fluence.py forward-solvers/phd-of-rho-and-z-two-layer-compute-fluence.py renamed to forward-solvers/phd-of-rho-and-z-compute-fluence.py

@@ -1,5 +1,5 @@
-# This is an example of python code using VTS to Compute photon hitting density for a two-layer 
-# medium at a given set of optical properties 
+# This is an example of python code using VTS to Compute photon hitting density for homogeneous 
+# medium at a given set of optical properties opRegions[0].
 # This sample uses ComputeFluence in place of calling FluenceOfRhoAndZ on the forward solver object
 # and it uses a distributed point source SDA Forward Solver
 #
@@ -35,9 +35,8 @@
 solver = DistributedPointSourceSDAForwardSolver()
 
 topLayerThickness = 5
-opRegions = Array.CreateInstance(IOpticalPropertyRegion, 2)
+opRegions = Array.CreateInstance(IOpticalPropertyRegion, 1)
 opRegions[0] = LayerOpticalPropertyRegion(DoubleRange(0, topLayerThickness, 2), OpticalProperties(0.1, 1, 0.8, 1.4))
-opRegions[1] = LayerOpticalPropertyRegion(DoubleRange(topLayerThickness, Double.PositiveInfinity, 2), OpticalProperties(0.01, 1, 0.8, 1.4))
 # Create the DoubleRange instance
 rhos_range = DoubleRange(0.1, 19.9, 100) # range of s-d separations in mm
 
@@ -105,4 +104,4 @@ def heatmap(values, x, y, x_label="", y_label="", title=""):
     return fig
 
 fluenceChart = heatmap(phdRowsToPlot.tolist(), allRhos.tolist(), list(zs), "ρ [mm]", "z [mm]", "log(phd(ρ, z) [mm-2])")
-fluenceChart.show(renderer="browser")
+fluenceChart.show(renderer="browser")

forward-solvers/phd-of-rho-and-z-two-layer.py

@@ -1,5 +1,7 @@
 # This is an example of python code using VTS to Compute photon hitting density for a two-layer 
-# medium at a given set of optical properties 
+# medium at a given set of optical properties.  A two-layer SDA forward solver is used to
+# compute the photon hitting density with optical properties defined in opRegions[0] and
+# opRegions[1], and top layer thickness defined in topLayerThickness [mm].
 #
 # Import PythonNet
 from pythonnet import load