Message: outp=p2.o  runtpe=p2.r  mctal=p2.m
      
Radial & Axial Dose Depth - Endovascular Radiation Therapy 32P beta source
c BEGINNING OF CELL CARDS
c
c NiTi Wire
   1   4   -6.45   -21   6 -11            imp:p,e 1      $ wire up to tip
   2   4   -6.45   -21  11 -10            imp:p,e 1      $ tip of the wire
c
c P-32 Source; 0.24 mm diameter
   3   5   -1.12     2  -3 -20            imp:p,e 1
c
c Tungsten Marker (19.3 g/cc)
   4   7   -19.3     3  -4 -20            imp:p,e 1
c
c NiTi Plug
   5   4   -6.45     4  -5 -20            imp:p,e 1
c
c NiTi Weld Hemispherical
   6   4   -6.45     5 -15                imp:p,e 1
c
c Symmetric Air Space
   7   6   -0.00120484     -5  2 -21 20   imp:p,e 1 
   8   6   -0.00120484     11 -2 -21 10   imp:p,e 1
c NiTi Tube
c     0.46 mm outside diameter; wall thickness 0.076 mm
   9   4   -6.45    21 -22 6 -5           imp:p,e 1       
c
c Blood within the Artery
c     (assumed to be water at 1.038 g/cc)
  10   1   -1.038    -22 1 -6             imp:p,e 1    
  11   1   -1.038    -22 5 15 -160        imp:p,e 1
c
c From the edge of the P-32 source to the end of the phantom y=0
c Rings within blood vessel wall.
 100   2   -1.04     23 -24   1 -101      imp:p,e 1
 101   2   -1.04     23 -24 101 -102      imp:p,e 1
 102   2   -1.04     23 -24 102 -103      imp:p,e 1
 103   2   -1.04     23 -24 103 -104      imp:p,e 1
 104   2   -1.04     23 -24 104 -105      imp:p,e 1
 105   2   -1.04     23 -24 105 -106      imp:p,e 1
 106   2   -1.04     23 -24 106 -107      imp:p,e 1
 107   2   -1.04     23 -24 107 -108      imp:p,e 1
 108   2   -1.04     23 -24 108 -109      imp:p,e 1
 109   2   -1.04     23 -24 109 -110      imp:p,e 1
 110   2   -1.04     23 -24 110 -111      imp:p,e 1
 111   2   -1.04     23 -24 111 -112      imp:p,e 1
 112   2   -1.04     23 -24 112 -113      imp:p,e 1
 113   2   -1.04     23 -24 113 -114      imp:p,e 1
 114   2   -1.04     23 -24 114 -115      imp:p,e 1
 115   2   -1.04     23 -24 115 -116      imp:p,e 1
 116   2   -1.04     23 -24 116 -117      imp:p,e 1
 117   2   -1.04     23 -24 117 -118      imp:p,e 1
 118   2   -1.04     23 -24 118 -119      imp:p,e 1
 119   2   -1.04     23 -24 119 -120      imp:p,e 1
 120   2   -1.04     23 -24 120 -121      imp:p,e 1
 121   2   -1.04     23 -24 121 -122      imp:p,e 1
 122   2   -1.04     23 -24 122 -123      imp:p,e 1
 123   2   -1.04     23 -24 123 -124      imp:p,e 1
 124   2   -1.04     23 -24 124 -125      imp:p,e 1
 125   2   -1.04     23 -24 125 -126      imp:p,e 1
 126   2   -1.04     23 -24 126 -127      imp:p,e 1
 127   2   -1.04     23 -24 127 -128      imp:p,e 1
 128   2   -1.04     23 -24 128 -129      imp:p,e 1
 129   2   -1.04     23 -24 129 -130      imp:p,e 1
 130   2   -1.04     23 -24 130 -131      imp:p,e 1
 131   2   -1.04     23 -24 131 -132      imp:p,e 1
 132   2   -1.04     23 -24 132 -133      imp:p,e 1
 133   2   -1.04     23 -24 133 -134      imp:p,e 1
 134   2   -1.04     23 -24 134 -135      imp:p,e 1
 135   2   -1.04     23 -24 135 -136      imp:p,e 1
 136   2   -1.04     23 -24 136 -137      imp:p,e 1
 137   2   -1.04     23 -24 137 -138      imp:p,e 1
 138   2   -1.04     23 -24 138 -139      imp:p,e 1
 139   2   -1.04     23 -24 139 -140      imp:p,e 1
 140   2   -1.04     23 -24 140 -141      imp:p,e 1
 141   2   -1.04     23 -24 141 -142      imp:p,e 1
 142   2   -1.04     23 -24 142 -143      imp:p,e 1
 143   2   -1.04     23 -24 143 -144      imp:p,e 1
 144   2   -1.04     23 -24 144 -145      imp:p,e 1
 145   2   -1.04     23 -24 145 -146      imp:p,e 1
 146   2   -1.04     23 -24 146 -147      imp:p,e 1
 147   2   -1.04     23 -24 147 -148      imp:p,e 1
 148   2   -1.04     23 -24 148 -149      imp:p,e 1
 149   2   -1.04     23 -24 149 -150      imp:p,e 1
 150   2   -1.04     23 -24 150 -151      imp:p,e 1
 151   2   -1.04     23 -24 151 -152      imp:p,e 1
 152   2   -1.04     23 -24 152 -153      imp:p,e 1
 153   2   -1.04     23 -24 153 -154      imp:p,e 1
 154   2   -1.04     23 -24 154 -155      imp:p,e 1
 155   2   -1.04     23 -24 155 -156      imp:p,e 1
 156   2   -1.04     23 -24 156 -157      imp:p,e 1
 157   2   -1.04     23 -24 157 -158      imp:p,e 1
 158   2   -1.04     23 -24 158 -159      imp:p,e 1
 159   2   -1.04     23 -24 159 -160      imp:p,e 1
c
c Radial Dose Depth - concentric cylinders
c   within the blood and this also defines the blood
 500   1   -1.038   (22 -200 1 -6):(22 -200 5 15 -160):(22 -200 6 -5)
                                          imp:p,e 1
 501   1   -1.038   200  -23 1 -160       imp:p,e 1
c   1 mm into the arterial wall will be defined in the axial dose above
c   remaining arterial wall
 502   2   -1.04     24  -25 1 -160       imp:p,e 1
c   within the water
 503   1   -1.00     25 -201 1 -160       imp:p,e 1
 504   1   -1.00    201 -202 1 -160       imp:p,e 1
 505   1   -1.00    202 -203 1 -160       imp:p,e 1
 506   1   -1.00    203 -204 1 -160       imp:p,e 1
 507   1   -1.00    204 -205 1 -160       imp:p,e 1
c
c Outside World
 999   0           205:-1:160             imp:p,e 0
c
c END OF CELL CARDS

c BEGINNING OF SURFACE CARDS
c
c Planes marking ends of catheter innards and water phantom
c      in the axial or y direction (units are cm)
  1  py    0.0          $ beginning of water phantom
  2  py   13.65         $ beginning of P-32 source cell
c                         The P-32 source is axially centered at 15.0 cm
  3  py   16.35         $ end of P-32 source cell, begin W marker
  4  py   16.45         $ end of W marker, begin NiTi Plug
  5  py   16.55         $ end of NiTi Plug
  6  py    5.7          $ end of NiTi guide wire 
c
c Sphere approximating the end of the NiTi wire
 10  sy   13.6346   0.0154 
 11  py   13.6346       $ plane cutting the sphere in half to approximate the
c                         end of the NiTi wire 
c
c Sphere approximating the end of the catheter
 15  sy   16.55     0.023
c
c Cylinders marking ends of catheter innards and water phantom
 20  cy    0.012        $ radius of P-32, W marker, NiTi wire
c                         and NiTi Plug
 21  cy    0.0154       $ inside radius of NiTi Tube
 22  cy    0.023        $ outside radius of NiTi Tube
 23  cy    0.15         $ inside radius of arterial wall
 24  cy    0.25         $ 1 mm into arterial wall
 25  cy    0.35         $ outside radius of arterial wall
c                         and the beginning of water phantom
c                         cylinder 1068 is the end of the water phantom
c
c Planes to be used in axial depth of water phantom (units are cm)
101  py    1.0
102  py    2.0
103  py    3.0
104  py    4.0
105  py    5.0
106  py    5.5
107  py    6.0
108  py    6.5
109  py    7.0
110  py    7.5
111  py    8.0
112  py    8.5
113  py    9.0
114  py    9.5
115  py   10.0
116  py   10.5
117  py   11.0
118  py   11.5
119  py   12.0
120  py   12.5
121  py   13.0
122  py   13.25
123  py   13.5
124  py   13.75
125  py   14.0
126  py   14.25
127  py   14.5
128  py   14.7
129  py   14.9
130  py   15.0
131  py   15.1
132  py   15.3
133  py   15.5
134  py   15.75
135  py   16.0
136  py   16.25
137  py   16.5
138  py   16.75
139  py   17.0
140  py   17.5
141  py   18.0
142  py   18.5
143  py   19.0
144  py   19.5
145  py   20.0
146  py   20.5
147  py   21.0
148  py   21.5
149  py   22.0
150  py   22.5
151  py   23.0
152  py   23.5
153  py   24.0
154  py   24.5
155  py   25.0
156  py   26.0
157  py   27.0
158  py   28.0
159  py   29.0
160  py   30.0
c
c Concentric cylinders including catheter walls which represent
c      radial depth into water phantom (units are cm)
200  cy    0.075     $ within blood
201  cy    0.5       $ within water
202  cy    1.0       $ within water
203  cy    2.0       $ within water
204  cy    5.0       $ within water
205  cy   15.0       $ to contain water and define outside of problem
c
c END OF SURFACE CARDS

c BEGINNING OF DATA CARDS
c
mode P E
c
print
c prdmp 2j 1
c
nps 10000000
c
c Source Definition
c Isotropic cylindrical volume source
sdef  pos=0 15 0 axs=0 1 0 rad=d1 ext=d2 erg=d3 par 3
si1 0.012
sp1 -21 1
sb1 -21 4
si2 -1.35 1.35
sp2 -21 0
sb2 -21 4
si3 0 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2
      0.225 0.25 0.275 0.3 0.325 0.35 0.375 0.4 0.425
      0.45 0.475 0.5 0.525 0.55 0.575 0.6 0.625 0.65
      0.675 0.7 0.725 0.75 0.775 0.8 0.825 0.85 0.875
      0.9 0.925 0.95 0.975 1
c
sp3 0 0.04275 0.0855 0.12825 0.171 0.21375 0.2565 0.29925
      0.342 0.38475 0.4275 0.47025 0.513 0.55575 0.5985
      0.64125 0.684 0.72675 0.7695 0.81225 0.855 0.89775
      0.9405 0.98325 1.026 1.06875 1.1115 1.15425 1.197
      1.23975 1.2825 1.32525 1.368 1.41075 1.4535 1.49625
      1.539 1.58175 1.6245 1.66725 1.71
c
c Weight Window Generator and Mesh Cards
wwg 6 0 0
c
c I think this is a cylindrical mesh with axis from -0.5 to 30.6 along
c  the y-axis and radii including 3.8 and 16 cm.  I'm not exactly sure
c  what the kmesh is doing unless it is cutting up the rings into wedges
mesh  geom=cyl  ref=0 15 0  origin=0.1 -0.5 0  axs=0 1 0  vec=1 0 0
      imesh  3.8 16    iints 9.5 4
      jmesh  10 21 31  jints 4 40 4
      kmesh  1         kints 2 
c
c Cell Importances
c
c imp:e 1 231r 1 0
c imp:p 1 231r 1 0
c
c Tallies
c
 fc6 Energy deposition (MeV/g) per source electron only from
      secondary photons (electrons produced from these photons
      deposit energy locally)
      f6 tallies are the cylindrical/radial dose depth
 f6:p  500 501 (100 58i 159) 502 503 504 505 506 507
 fm6   1.6e-10      $ convert MeV/g to Gy
c
 f16c f16 tallies are for the axial dose depth
 f16:p  100 58i 159
 fm16   1.6e-10     $ convert MeV/g to Gy
c
 fc26 Energy deposition (MeV/g) to midway point 1 mm into artery wall 
      - Question iv) $ IE. FIND THE MAX DOSE RATE TO THE ARTERY WALL
 f26:p  (127 128)
 fm26   1.6e-10     $ convert MeV/g to Gy
c
 fc8 Total energy deposition (MeV) in each of the cells per
      particle - Question iii)
      f8 tallies are the cylindrical/radial dose depth
 *f8:p 500 501 (100 58i 159) 502 503 504 505 506 507
c
 fc18 f18 tallies are for the axial dose depth - Question ii)
 *f18:p   100 58i 159
c
 fc28 f28 tallies are for Dose to midway point 1 mm into artery wall
 *f28:p   (127 128)
c
c Beginning of Material Data Cards
c
c Water
c density = 1.00 g/cc
 m1  001000  -0.112 
     008000  -0.888 
c
c Soft tissue
c density = 1.04 g/cc
 m2  001000  -0.101
     006000  -0.111
     007000  -0.026
     008000  -0.762
c
c Plaque
c density = 1.45 g/cc
 m3  001000  -0.06141
     006000  -0.40302
     007000  -0.01919
     008000  -0.29884
     015000  -0.10254
     020000  -0.115
c
c NiTi - alloy
c density = 6.45 g/cc
 m4  024000  -0.0025
     028000  -0.558
     022000  -0.4395
c
c P-32 Source in a polyethylene matrix
c density = 1.12 g/cc
 m5  001000  -0.09677
     006000  -0.38710
     008000  -0.51613
c
c Air
c density = 0.00120484
 m6  007000  -0.790
     008000  -0.210
c
c Tungsten
c density = 19.3 g/cc
 m7  074000  1
c
c END OF DATA CARDS
                         

