montecarlomaplet.mw

A Monte-Carlo Maplet for the Study of the Optical Properties of Biological Tissues

PP:= pointplot({seq([RTd[n],temp[n]],n=1..Nr-1)}, labels=["radius (cm)", Transmittance], axes=boxed,colour=black):
#PPa:= pointplot({seq([RTd[n],tempa[n]],n=1..Nr-1)}, labels=["radius (cm)", Transmittance], axes=boxed,colour=red):
pprT:=[seq([RTd[n],temp[n]],n=1..Nr-1)];
#pprTa:=[seq([RTd[n],tempa[n]],n=1..Nr-1)];
save pprT,`Tra_rad_graph_points.m`:
RETURN(display[plots](PP));
<Font foreground="[153,0,153]">end:</Font><Font foreground="[0,0,0]">

</Font></Font></Text-field>
</Input>
</Group>
</Section>
</Section>
<Section collapsed="false" MultipleChoiceAnswerIndex="-1" MultipleChoiceRandomizeChoices="false" TrueFalseAnswerIndex="-1" EssayAnswerRows="5" EssayAnswerColumns="60"><Title>
<Text-field style="Heading 1" layout="Heading 1"><Font foreground="[51,51,255]">II - The Maplet</Font></Text-field></Title>
<Section collapsed="false" MultipleChoiceAnswerIndex="-1" MultipleChoiceRandomizeChoices="false" TrueFalseAnswerIndex="-1" EssayAnswerRows="5" EssayAnswerColumns="60"><Title>
<Text-field style="Heading 2" layout="Heading 2"><Font size="18">1. Construction of the Maplet</Font> </Text-field></Title>
<Group labelreference="L11" drawlabel="true">
<Input>
<Text-field prompt="&gt; " style="Maple Input" layout="Normal"><Equation executable="true" style="2D Input" input-equation="Typesetting:-mrow(Typesetting:-mi(&quot;&quot;))">LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2OVEhRicvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GJy8lJXNpemVHUSMxMEYnLyUlYm9sZEdRJmZhbHNlRicvJSdpdGFsaWNHUSV0cnVlRicvJSp1bmRlcmxpbmVHRjcvJSpzdWJzY3JpcHRHRjcvJSxzdXBlcnNjcmlwdEdGNy8lK2ZvcmVncm91bmRHUShbMCwwLDBdRicvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYnLyUnb3BhcXVlR0Y3LyUrZXhlY3V0YWJsZUdGOi8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjcvJSpjb252ZXJ0ZWRHRjcvJStpbXNlbGVjdGVkR0Y3LyUscGxhY2Vob2xkZXJHRjcvJTBmb250X3N0eWxlX25hbWVHUSkyRH5JbnB1dEYnLyUqbWF0aGNvbG9yR0ZDLyUvbWF0aGJhY2tncm91bmRHRkYvJStmb250ZmFtaWx5R0YxLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy8lKW1hdGhzaXplR0Y0</Equation><Font size="10" foreground="[0,0,255]">use Maplets[Elements] in
with(Maplets[Elements]):</Font><Font size="10">
</Font><Font size="9">    MonteCarloMaplet:= Maplet('onstartup' = RunWindow( 'W0' ),
 Font['F00']( 'family' = &quot;times&quot;,bold,'size'=12 ),
               Font['F0']( 'family' = &quot;arial&quot;,bold,'size'=12 ),
               Font['F1']( 'family' = &quot;arial&quot;,bold,'size'=14 ),
               Font['F2']( 'family' = &quot;times&quot;,bold, 'size'=20 ),
               Font['F3']( 'family' = &quot;times&quot;,bold, 'size'=14 ),

</Font><Font background="[255,255,153]" opaque="true" size="14">Window['W0']('title'= &quot;MAIN WINDOW&quot;,</Font><Font size="9">
BoxRow(
BoxColumn( #### BoxColumn 1 starts
'hscroll'='always',
'vscroll'='always','inset'=0,'spacing'=15,'valign'='top',
       Label(&quot;<Font background="[153,255,102]" opaque="true">INPUT PANEL</Font>      &quot;,'font'=Font(Arial,bold,16)),
[Button(&quot;Explanation of INPUT &quot;,'foreground'=black,'background'=white, 'font'='F1', RunWindow('Win'))],
[BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Number of photons </Font> &quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot;Enter a number &gt;500  &quot;), 
           GridCell('halign'='left',
              TextField['TF_2']('width'=6,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;   &quot;)
))) <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
</Font>,<Font size="9">
BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Entry Angle </Font> &quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot;Enter a number between 0 and 90  &quot;), 
           GridCell('halign'='left',
              TextField['TF_0']('width'=6,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;   &quot;)
)))] <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
</Font>,<Font size="9">
[BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Number of tissue layers</Font> &quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot;Enter an integer value  &quot;), 
           GridCell('halign'='left',
              TextField[TF_1]('width'=3,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;    &quot;)
))) <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; boxrow</Font>
,

BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Layers' refraction indices</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot; [n1, n2, ...]     &quot;), 
           GridCell('halign'='left',
              TextField['TF_3']('width'=5,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;   &quot;)
)))] <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; boxrow</Font>
,

[BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Layers' scattering coefficients</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot; [cs1, cs2, ...]      &quot;), 
           GridCell('halign'='left',
              TextField['TF_4']('width'=5,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;  1/cm &quot;)
))) <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
,
BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Layers' absorption coefficients </Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot; [ca1, ca2, ...]     &quot;), 
           GridCell('halign'='left',
              TextField['TF_5']('width'=5,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;  1/cm &quot;)
)))] <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
,
[BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Layers' anisotropy coefficients</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot;[g1, g2, ...]     &quot;), 
           GridCell('halign'='left',
              TextField['TF_6']('width'=5,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;   &quot;)
))) <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
,
BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Layers' depths</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot; [d1, d2, ...]     &quot;), 
           GridCell('halign'='left',
              TextField['TF_7']('width'=5,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;  cm &quot;)
)))] <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
),
[BoxColumn('hscroll'='always',
'vscroll'='always','inset'=0,'spacing'=15,'valign'='top',
       Label(&quot;   SELECT GEOMETRY TYPE</Font>   <Font size="9">  &quot;,'font'=Font(Arial,bold,16)),

BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Cylindrical Layers</Font>&quot; ,'font'='F1',     
GridRow('halign'='left',
           GridCell('halign'='left',&quot;Input maximum radius  &quot;), 
           GridCell('halign'='left',
              TextField['TF_r']('width'=3,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;  cm &quot;)
</Font>))), <Font size="9" foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font><Font size="9">

CheckBox['ChB1']('caption'=&quot;For cylindrical layers check here&quot;, 'value' = 'false' ),   
#[Button(&quot;Proceed &quot;,'background'=white,'foreground'=black, 'font'='F1', RunWindow('W1'))],
BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Prismatic Layers</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot;Input maximum x,y dimensions [x,y]    &quot;), 
           GridCell('halign'='left',
              TextField['TF_xy']('width'=6,'background'=grey, 'editable'='true', 'font'='F1')),  
GridCell('halign'='left',&quot;  cm &quot;)
</Font>)))<Font size="9">, <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>

CheckBox['ChB2']('caption'=&quot;For prismatic layers check here&quot;, 'value' = 'false' ),   

BoxRow(
GridLayout('border'='true', 'inset'=0, 
          'caption'=&quot;<Font background="[255,255,102]" opaque="true">Infinitely Long Layers</Font>&quot; ,'font'='F1',        
GridRow('halign'='left',
           GridCell('halign'='left',&quot; Infinite radius assumed   &quot;)#, 
           
</Font>))
)<Font size="9"> ,  <Font foreground="[0,0,0]">#  end gridrow ; gridlayout ; BoxRow</Font>
CheckBox['ChB3']('caption'=&quot;For inifnite layers check here&quot;, 'value' = 'false' ), 

MathMLViewer['ML']('foreground'=black,'background'=grey,'height'=2,
</Font>       'width'=1),<Font size="9"> 
Button(&quot;CLEAR Run Status BOX&quot;, 
Action(</Font>SetOption('ML' = &quot;&quot;)))<Font size="9">,
CheckBox['ChB4']('caption'=&quot;If only numerical results are required, check here&quot;, 'value' = 'false' ),   
Button(&quot;Start Simulation&quot;,'foreground'=black, 'background'=white,'font'='F00',Evaluate('ML' = 'SimulationMC()')) 
</Font> ,<Font size="9"> 
[Button(&quot;Explanation of OUTPUT&quot;,'background'=white, 'foreground'=black,'font'='F1',RunWindow('Wout'))],
[Button(&quot;Proceed to numerical results&quot;,'background'=white,'foreground'=black, 'font'='F1', RunWindow('W1')),
Button(&quot;Proceed to graphical results&quot;,'background'=white,'foreground'=black, 'font'='F1', RunWindow('W2'))]),
Button['B3']( &quot;Close Window&quot;, Shutdown() )]
)),   <Font foreground="[0,0,0]"># end of boxcolumn,boxrow,window</Font>

</Font><Font background="[255,255,153]" opaque="true" size="14">Window['Win']('title'= &quot;Explanation of Input&quot;,</Font><Font size="9">
BoxRow(
      BoxColumn( 
                 TextBox['TF_I']('editable' = 'false',16..60,  
 &quot; The INPUT WINDOW consists of two columns where data must be entered in order for 
   the Monte Carlo simulation to proceed. 
   On the left column  the entries are:
     number of photons (integer - &gt;500 recommended)
     entry angle (angle in degrees with respect to layer's surface) 
     number of tissue layers (integer)
     layers indices of refraction     (list) 
     layers absorption coefficients   (list)
     layers scattering coefficients   (list)
     layers anisotropy coefficient    (list)
     layers depths                    (list)
  Note that the lists have as many entries as the number of layers; the list format requires that 
  entries be separated by commas and enclosed in square brackets. Ranges of acceptable input values 
  are given in the literature (eg., </Font><Font size="10">L.Wang, S.L.Jacques and L.Zheng, Computer Methods and 
                                     Programs in Biomedicine, vol 47, (1995), 131-146)</Font></Text-field>
<Text-field prompt="&gt; " style="Maple Input" layout="Normal"><Font size="9">  
  On the right column, three options, for the layers geometry type, are presented: 
      cylindrical, 
      prismatic and 
      infinite layers.