Sự khác nhau về số đếm diện tích đỉnh giữa thực nghiệm và mô phỏng được trình bày ở bảng 3.5
Bảng 3.5. Độ sai biệt về số đếm diện tích đỉnh giữa thực nghiệm và mô phỏng (%)
Nguồn E (keV) 20cm 15m 10cm 5cm 0cm 60 Co 1173 4,34 6,02 1,09 12,84 10,55 1332 1,82 5,21 6,11 11,85 12,68 54Mn 835 1,08 17,55 21,78 19,15 11,28 22 Na 1275 1,76 12,65 3,09 12,17 15,83 511 3,74 6,70 5,90 11,14 11,07 133 Ba 81 2,69 0,87 0,29 5,75 2,49 356 1,51 1,37 3,79 0,05 3,83 57Co 122 13,79 2,67 1,57 1,85 2,63 137 Cs 662 2,75 4,52 6,08 7,73 8,52 109 Cd 88 1,70 0,70 17,81 2,72 5,66 65 Zn 1115 9,85 10,77 9,46 10,49 6,43
32
Có sự sai biệt khá lớn ở 22Na, 60Co ở vị trí 5cm và 0cm do hiện tượng trùng phùng trong detector. 54Mn ở các vị trí 0cm, 5cm, 10cm, 15cm và 109Cd ở vị trí 10cm cũng có sự sai biệt khá lớn, có thể do sai sót trong quá trình đo hoặc cấu hình nguồn trong mô phỏng có khác biệt so với cấu hình thực tế.
Các phổ so sánh giữa mô phỏng và thực nghiệm của 60Co, 137Cs 133Ba ở vị trí 10cm và 22Na ở vị trí 5cm.
Hình 3.3. Phổ so sánh thực nghiệm và mô phỏng của 60Co ở vị trí 10cm
33
Hình 3.5. Phổ so sánh thực nghiệm và mô phỏng của 133Ba ở vị trí 10cm
34
Qua các phổ so sánh giữa thực nghiệm và mô phỏng, ta thấy các đỉnh mô phỏng phù hợp với thực nghiệm, các vùng tán xạ cũng phù hợp ở các nguồn đơn năng và đa năng. Chỉ có vùng năng lượng cao do sự xuất hiện của phông đến từ môi trường nên có sự khác biệt giữa mô phỏng và thực nghiệm. Đặc biệt, ở phổ 22Na trong thực nghiệm xuất hiện đỉnh do hiện tượng trùng phùng, nhưng trong mô phỏng không có.
Hiệu suất ghi của detector được mô tả qua hình 3.7
Hình 3.7. So sánh hiệu suất ghi nhận của detector giữa thực nghiệm và mô phỏng của 8 nguồn ở các vị trí khác nhau
Về mặt hiệu suất ghi nhận của detector, giữa thực nghiệm và mô phỏng có sự phù hợp khá tốt,điều đó thể hiện qua bảng 3.6.
Có sự sai biệt khá lớn 60Co và 22Na ở vị trí 0cm và 5cm do trong thực nghiệm xảy ra hiệu ứng trùng phùng trong quá trình ghi nhận.
35
Nguồn 54Mn cũng có sự sai biệt khá lớn ở vị trí 5cm và 10cm, có thể do thao tác trong quá trình đo bị lỗi ( đặt vị trí nguồn bị lệch) hoặc do cấu hình nguồn mô phỏng không giống hoàn toàn với thực nghiệm.
Bảng 3.6. Độ sai biệt về hiệu suất ghi nhận giữa thực nghiệm và mô phỏng (%)
Nguồn E (keV) 20cm 15cm 10cm 5cm 0cm 60 Co 1173 4,38 6,24 0,95 14,56 11,62 1332 1,83 5,48 6,49 13,42 14,50 54 Mn 835 4,31 8,70 14,58 10,86 1,02 22 Na 1275 11,10 13,50 2,30 12,89 17,79 511 4,75 8,08 7,16 13,48 13,39 133Ba 81 0,95 0,51 0,65 5,71 6,03 356 0,27 0,42 1,93 1,84 6,91 57Co 122 11,38 3,60 0,72 2,73 3,56 137Cs 662 3,26 5,17 6,91 8,83 7,47 109 Cd 88 3,32 2,28 0,51 4,41 7,66 65 Zn 1115 5,26 6,35 4,82 6,02 1,41 3.4. Nhận xét chương 3
Trong chương này đã trình bày kết quả về số đếm diện tích đỉnh và hiệu suất ghi của detector trong thực nghiệm và mô phỏng.
36
KẾT LUẬN VÀ KIẾN NGHỊ
Kết luận
Trong khóa luận đã xây dựng được hệ đo thực nghiệm và mô phỏng và thu được kết quả sau:
Số đếm diện tích đỉnh giữa thực nghiệm và mô phỏng có sự phù hợp khá tốt qua số liệu thu được từ việc đo và mô phỏng 8 nguồn: 60Co, 54Mn, 22Na, 133Ba, 57Co, 137Cs, 109Cd, 65Zn với 5 khoảng cách: 20cm, 15cm, 10cm, 5cm và 0cm. Chỉ có 22Na, 60Co ở vị trí 5cm và 0cm, 54Mn ở các vị trí 0cm, 5cm, 10cm, 15cm và 109
Cd ở vị trí 10cm có sự sai biệt khá lớn.
Hiệu suất ghi nhận của detector giữa thực nghiệm và mô phỏng tốt, chỉ có 60Co và 22Na ở vị trí 0cm và 5cm, 54Mn ở vị trí 5cm và 10cm có sự sai biệt khá lớn.
Nguyên nhân của các sự sai biệt trên là do: -Hiệu ứng trùng phùng ở nguồn 60Co và 22Na.
-Thao tác khi đo không chính xác (vị trí đặt nguồn bị lệch)
-Cấu hình nguồn thực nghiệm với cấu hình nhà sản xuất đưa ra không hoàn toàn giống nhau dẫn đến việc mô phỏng không hoàn toàn chính xác.
Mặc dù vậy chương trình PENELOPE vẫn mô phỏng khá tốt và đầy đủ các hiện tượng vật lý.
Kiến nghị
Nếu có thời gian hoặc các khóa luận tiếp theo có thể: Kiểm tra lại những kết quả có độ sai biệt lớn.
Mô phỏng bằng các chương trình mô phỏng khác nhằm đánh giá sự khác nhau giữa các chương trình.
37
TÀI LIỆU THAM KHẢO Tiếng Việt
[1]. Bùi Hải Âu (2009), Xây dựng chương trình mô phỏng hiệu suất đỉnh năng lượng toàn phần của detector”, Khóa luận tốt nghiệp, Trường Đại Học KHTN, Đại học Quốc Gia Tp.HCM.
[2]. Trần Phong Dũng, Châu Văn Tạo, Nguyễn Hải Dương (2005), Phương pháp ghi bức xạ ion hóa”, Trường Đại Học KHTN, Đại học Quốc Gia Tp.HCM. [3]. Ngô uang Huy (2006), Cơ sở vật lý hạt nhân”, Nhà xuất bản khoa học và kỹ
thuật.
[4]. Mai Văn Nhơn (2001), Vật lý hạt nhân đại cương”, Nhà xuất bản Đại học quốc gia Tp.HCM.
[5]. Trần Thiện Thanh (2013), Hiệu chỉnh phổ gamma bằng phương pháp Monte Carlo , Luận án Tiến sĩ, Trường Đại Học KHTN, Đại học Quốc Gia Tp.HCM. [6]. Trần Thiện Thanh, Châu Văn Tạo, Trương Thị Hồng Loan, Mai Văn Nhơn
(2012), “Nghiên cứu đường cong hiệu suất đỉnh năng lượng toàn phần sử dụng chương trình DETEFF”, Tạp chí Khoa Học, Trường Đại Học Sư Phạm Tp.HCM. Số 10 năm 2012.
Tiếng nước ngoài
[7]. F.Salvat, J.M.Fernandez-Varea, J.Sempau (2008), PENELOPE-2008: A code system for Monte Carlo simulation of electron and photon transport, Publisher OECD, Spain.
[8]. Hu-Xia Shi, Bo-Xian Chen, Ti-Zhu Li, Di Yun (2002), Precise Monte Carlo Simulation of gamma-ray response function for an NaI(Tl) detector, Applied Radiation and Isotopes; 517-524.
[9]. S.Ashrafi, S.Anvarian, S.Sobhanian (2005), Monte-Carlo Modeling of a NaI(Tl) scintillator, Journal of Radioanalytical and Nuclear Chemistry, Vol.269, 95-98.
38
[10]. U.Akar Tarim, E.N.Ozmutlu, O.Gurler, S.Yalcin (2012), The Effect of the housing material on the NaI(Tl) detector response function, Journal of Radioanalytical and Nuclear Chemistry, Vol.293, 425-429.
Các trang web
39
PHỤ LỤC
Phụ lục a- cấu hình mô phỏng ở vị trí 0cm
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx university of science, vnu
detector model 802-3x3
support and point source (0cm)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 1) plane z=0cm indices=( 0, 0, 0, 1, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 2) plane z=7.62cm indices=( 0, 0, 0, 1,-1) z-scale=( 7.620000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 3) plane z=7.78cm indices=( 0, 0, 0, 1,-1) z-scale=( 7.780000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 4) plane z=7.89cm indices=( 0, 0, 0, 1,-1) z-scale=( 7.890000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 5) plane z=7.94cm indices=( 0, 0, 0, 1,-1) z-scale=( 7.940000000000000e+00, 0) 00000000000000000000000000000000000000000000000000000000000000 00
40 surface ( 6) plane z=-0.3cm indices=( 0, 0, 0, 1, 1) z-scale=( 0.300000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 7) plane z=-3.3cm indices=( 0, 0, 0, 1, 1) z-scale=( 3.300000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 8) plane z=8.14cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 8.140000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 9) plane z=28.14cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 2.814000000000000e+01, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 10) plane z=28.29cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 2.829000000000000e+01, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 11) plane z=28.44cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 2.844000000000000e+01, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 12) plane z=7.14cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 7.140000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 13) plane z=5.14cm #support
41 indices=( 0, 0, 0, 1,-1)
z-scale=( 5.140000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 14) plane x=4.9cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 4.900000000000000e+00, 0)
theta=( 90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 15) plane x=-4.9cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 4.900000000000000e+00, 0)
theta=(-90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 16) plane x=-5.08cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 5.080000000000000e+00, 0)
theta=(-90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 17) plane y=4.475cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 4.475000000000000e+00, 0)
theta=( 90.00000000000000e+00, 0) deg phi=( 90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 18) plane y=-4.475cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 4.475000000000000e+00, 0)
theta=( 90.00000000000000e+00, 0) deg phi=(-90.00000000000000e+00, 0) deg
42
0000000000000000000000000000000000000000000000000000000000000000 surface ( 19) plane y=5.3cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 5.300000000000000e+00, 0)
theta=( 90.00000000000000e+00, 0) deg phi=( 90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 20) plane y=-5.3cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 5.300000000000000e+00, 0)
theta=( 90.00000000000000e+00, 0) deg phi=(-90.00000000000000e+00, 0) deg
0000000000000000000000000000000000000000000000000000000000000000 surface ( 22) plane z=8.39cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 8.390000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 23) plane z=8.44cm #support
indices=( 0, 0, 0, 1,-1)
z-scale=( 8.440000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 24) plane z=8.467cm #source
indices=( 0, 0, 0, 1,-1)
z-scale=( 8.467000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 25) plane z=8.475cm #source
indices=( 0, 0, 0, 1,-1)
z-scale=( 8.475000000000000e+00, 0)
43
surface ( 26) plane z=8.785cm #source indices=( 0, 0, 0, 1,-1)
z-scale=( 8.785000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 27) plane z=8.825cm #source
indices=( 0, 0, 0, 1,-1) z-scale=( 8.825000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 101) cylinder r=3.81cm indices=( 1, 1, 0, 0,-1) x-scale=( 3.810000000000000e+00, 0) y-scale=( 3.810000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 102) cylinder r=4.0cm indices=( 1, 1, 0, 0,-1) x-scale=( 4.000000000000000e+00, 0) y-scale=( 4.000000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 103) cylinder r=4.05cm indices=( 1, 1, 0, 0,-1) x-scale=( 4.050000000000000e+00, 0) y-scale=( 4.050000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 104) cylinder r=4.13cm indices=( 1, 1, 0, 0,-1) x-scale=( 4.130000000000000e+00, 0) y-scale=( 4.130000000000000e+00, 0) 0000000000000000000000000000000000000000000000000000000000000000 surface ( 105) cylinder r=1.15cm #support
44 indices=( 1, 1, 0, 0,-1)
x-scale=( 1.150000000000000e+00, 0) y-scale=( 1.150000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 106) cylinder r=1.3cm #support
indices=( 1, 1, 0, 0,-1)
x-scale=( 1.300000000000000e+00, 0) y-scale=( 1.300000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 107) cylinder r=3.4cm #support
indices=( 1, 1, 0, 0,-1)
x-scale=( 3.400000000000000e+00, 0) y-scale=( 3.400000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 108) cylinder r=4.1cm #support
indices=( 1, 1, 0, 0,-1)
x-scale=( 4.100000000000000e+00, 0) y-scale=( 4.100000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 109) cylinder r=0.25cm #source
indices=( 1, 1, 0, 0,-1)
x-scale=( 0.250000000000000e+00, 0) y-scale=( 0.250000000000000e+00, 0)
0000000000000000000000000000000000000000000000000000000000000000 surface ( 110) cylinder r=1.27cm #source
indices=( 1, 1, 0, 0,-1)
x-scale=( 1.270000000000000e+00, 0) y-scale=( 1.270000000000000e+00, 0)
45 surface ( 201) sphere r=100 indices=( 1, 1, 1, 0,-1) x-scale=( 1.000000000000000e+02, 0) y-scale=( 1.000000000000000e+02, 0) z-scale=( 1.000000000000000e+02, 0) 0000000000000000000000000000000000000000000000000000000000000000 body ( 1) material( 1)
surface ( 1), side pointer=( 1) surface ( 2), side pointer=(-1) surface ( 101), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 2)
material( 2)
surface ( 1), side pointer=( 1) surface ( 3), side pointer=(-1) surface ( 102), side pointer=(-1) body ( 1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 3)
material( 3)
surface ( 3), side pointer=( 1) surface ( 4), side pointer=(-1) surface ( 102), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 4)
material( 5)
surface ( 6), side pointer=( 1) surface ( 1), side pointer=(-1)
46 surface ( 102), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 5)
material( 4)
surface ( 6), side pointer=( 1) surface ( 5), side pointer=(-1) surface ( 103), side pointer=(-1) body ( 2) body ( 3) body ( 4) 0000000000000000000000000000000000000000000000000000000000000000 body ( 6) material( 4)
surface ( 7), side pointer=( 1) surface ( 6), side pointer=(-1) surface ( 104), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 101) #source
material( 9)
surface ( 24), side pointer=( 1) surface ( 25), side pointer=(-1) surface ( 109), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 102) #source
material( 8)
surface ( 25), side pointer=( 1) surface ( 26), side pointer=(-1) surface ( 109), side pointer=(-1)
47 body ( 103) #source
material( 7)
surface ( 22), side pointer=( 1) surface ( 27), side pointer=(-1) surface ( 110), side pointer=(-1) body ( 101)
body ( 102)
0000000000000000000000000000000000000000000000000000000000000000 body ( 201) #support
material( 6)
surface ( 12), side pointer=(-1) surface ( 13), side pointer=( 1) surface ( 14), side pointer=(-1) surface ( 15), side pointer=(-1) surface ( 17), side pointer=(-1) surface ( 18), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 202) #support
material( 6)
surface ( 5), side pointer=(-1) surface ( 12), side pointer=( 1) surface ( 108), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 203) #support
material( 6)
surface ( 8), side pointer=(-1) surface ( 5), side pointer=( 1) surface ( 107), side pointer=(-1)
48 body ( 205) #support
material( 6)
surface ( 22), side pointer=(-1) surface ( 8), side pointer=( 1) surface ( 105), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 206) #support
material( 6)
surface ( 23), side pointer=(-1) surface ( 22), side pointer=( 1) surface ( 106), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 207) #support
material( 6)
surface ( 9), side pointer=(-1) surface ( 23), side pointer=( 1) surface ( 14), side pointer=(-1) surface ( 15), side pointer=(-1) surface ( 17), side pointer=(-1) surface ( 18), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 208) #support
material( 6)
surface ( 10), side pointer=(-1) surface ( 9), side pointer=( 1) surface ( 105), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 209) #support
49 surface ( 11), side pointer=(-1)
surface ( 10), side pointer=( 1) surface ( 106), side pointer=(-1)
0000000000000000000000000000000000000000000000000000000000000000 body ( 210) #support
material( 7)
surface ( 11), side pointer=(-1) surface ( 13), side pointer=( 1) surface ( 14), side pointer=(-1) surface ( 16), side pointer=(-1) surface ( 19), side pointer=(-1) surface ( 20), side pointer=(-1) body ( 103) body ( 201) body ( 202) body ( 203) body ( 205) body ( 206) body ( 207) body ( 208) body ( 209) 000000000000000000000000000000000000000000000000000000000000000 body (1000) material( 6)
surface ( 201), side pointer=(-1) body ( 5)
body ( 6) body ( 210) end
50
Phụ lục b- file input của 133Ba ở vị trí 0cm: title gamma spectrometry with nai(tl) model 802-3x3. >>>>>>>> source definition.
skpar 2 [primary particles: 1=electron, 2=photon, 3=positron] 133ba gamma-ray spectrum [lara database].
spectr 3.0625e4 1.3e-1 [e bin: lower-end and total probability] spectr 3.0625e4 1.0e-35 [e bin: lower-end and total probability] spectr 3.0973e4 2.5e-1 [e bin: lower-end and total probability] spectr 3.0973e4 1.0e-35 [e bin: lower-end and total probability] spectr 3.5053e4 7.1e-2 [e bin: lower-end and total probability] spectr 3.5053e4 1.0e-35 [e bin: lower-end and total probability] spectr 3.5900e4 1.8e-2 [e bin: lower-end and total probability] spectr 3.5900e4 1.0e-35 [e bin: lower-end and total probability] spectr 5.3162e4 8.4e-3 [e bin: lower-end and total probability] spectr 5.3162e4 1.0e-35 [e bin: lower-end and total probability] spectr 7.9614e4 1.0e-2 [e bin: lower-end and total probability] spectr 7.9614e4 1.0e-35 [e bin: lower-end and total probability] spectr 8.0998e4 1.3e-1 [e bin: lower-end and total probability] spectr 8.0998e4 1.0e-35 [e bin: lower-end and total probability] spectr 1.6061e5 2.5e-3 [e bin: lower-end and total probability] spectr 1.6061e5 1.0e-35 [e bin: lower-end and total probability] spectr 2.2324e5 1.8e-3 [e bin: lower-end and total probability] spectr 2.2324e5 1.0e-35 [e bin: lower-end and total probability] spectr 2.7640e5 2.8e-2 [e bin: lower-end and total probability] spectr 2.7640e5 1.0e-35 [e bin: lower-end and total probability] spectr 3.0285e5 7.2e-2 [e bin: lower-end and total probability] spectr 3.0285e5 1.0e-35 [e bin: lower-end and total probability] spectr 3.5601e5 2.4e-1 [e bin: lower-end and total probability] spectr 3.5601e5 1.0e-35 [e bin: lower-end and total probability]
51
spectr 3.8385e5 3.5e-2 [e bin: lower-end and total probability] spectr 3.8385e5 1.0e-35 [e bin: lower-end and total probability] sposit 0 0 0.8571e1 [coordinates of the source] sbody 7
scone 0 0 180 [conical beam; angles in deg] .
>>>>>>>> material data and simulation parameters.
mfname nai.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname al2o3.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname si.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname al.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname sio2.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname air.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname polymethyl.mat [material file, up to 20 chars] msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname epoxy.mat [material file, up to 20 chars] msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] mfname ba.mat [material file, up to 20 chars]
msimpa 1.0e3 1.0e3 1.0e3 0.1 0.1 1.0e3 1.0e3 [eabs(1:3),c1,c2,wcc,wcr] .
>>>>>>>> geometry definition file.
geomfn 802-3x3_0.geo [geometry file, up to 20 chars] .
52
>>>>>>>> energy deposition detectors (up to 25).
endetc 0.5e3 2.0801e6 2014 [energy window and number of bins] edbody 1 [active body; one line for each body]
.
>>>>>>>> job properties
resume dump.dmp [resume from this dump file, 20 chars] dumpto dump.dmp [generate this dump file, 20 chars] dumpp 3600 [dumping period, in sec] .
nsimsh 2e7 [desired number of simulated showers] time 1e9 [allotted simulation time, in sec] .