sample.cfg

蒙特卡罗模拟光子成像C语言版,代码简洁专业

; From a semicolon to the end of the line is a comment

;; Number of photons (1e6, a small run) [MC only]
nphotons = 1000000
nphotons = 100000

;; Seed random number generator [MC only]
seed    = 393571066

;; Modulation frequence [MC only]
freq    = 70

;; Source position.  All lengths in mm.  Comma's and equal signs (=)
;;  are always optional.  White space is not significant.

source { pos = [ 50 51 1 ]
         dir = [  0  0 1 ]
         rad = 0.250       }

;; Times are in seconds.  This simulations covers a range of 10*0.5 = 5 ns

start_time = 0.00e-9
gate_width = 1.00e-9
ngate      = 20

;; segmentation file to load

segfile  demo_sphere.bin

;; Voxel size.  Square voxels only, for now (dz = dy = dx)

dx         2

;; Size of system (in voxels)

nxvox = 50
nyvox = 50
nzvox = 25

;; Region of interest within system volume (also in voxels)
;; If left out, 2pt file covers entire volume.

image_x 0 49
image_y 0 49
image_z 0 24

; The number of tissue types is computed on the fly.  Tissue types
; are numbered in the order they're declared.  Tissue type 0 is always
; implicitly air (n=1).  Any photons ending up in tissue type 0 is
; assumed to have left the medium, never to return again.  If you want
; "re-entrant air", define another tissue type and give it physically
; meaningful optical properties.
;
; There MUST be one tissue type for every tissue mentioned in the
; segmentation file.

tissue { mua = 0.004 mus = 0.800 g = 0.0 n = 1.4 }	; background
tissue { mua = 0.004 mus = 3.800 g = 0.0 n = 1.5 }	; sphere

; As with tissue types, there's no need to specify the number of
;  detectors, just their location

detector { pos = [ 10 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 10 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 20 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 30 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 40 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 50 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 60 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 70 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 80 90 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 10 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 20 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 30 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 40 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 50 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 60 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 70 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 80 48 ] dir = [ 0 0 -1 ] rad = 2 }
detector { pos = [ 90 90 48 ] dir = [ 0 0 -1 ] rad = 2 }