-
Notifications
You must be signed in to change notification settings - Fork 1
Home
Mizukoshi Keita edited this page Nov 12, 2017
·
7 revisions
- ネットワークに接続された,Unix系のOS(virtual machine可)
- git
- python
- root
準備の詳細はこちら
NeuCBOT tutorial Mizukoshi Keita Osaka Univ. @Waseda Univ.
中性子の生成源
- 自発核分裂: Cf, U など
- 核破砕
- (α,n)反応 ←今回はこの計算方法
(α,n)反応の断面積は非常に小さいので,Geant4などの モンテカルロを使って計算することは困難
- SOURCES-4Cが業界標準だった
- エネルギー領域が狭い(0 - 5 MeV)
- データが更新されていない
- 高価
“NeuCBOT”をもちいて計算しよう!
- 最新の(α,n) 計算ツール
- 前述のSOURCESの欠点をすべて解決.
- Database: TALYS
- UI: python + GitHub ←手軽に使える！
>python --version
Python 2.7.10
>git --version
git version 2.13.4
#option
>root -v
------------------------------------------------------------
| Welcome to ROOT 6.10/08 http://root.cern.ch |
| (c) 1995-2017, The ROOT Team |
| Built for macosx64 |
| From tag v6-10-08, 16 October 2017 |
| Try '.help', '.demo', '.license', '.credits', '.quit'/'.q' |
————————————————————————————————————————————————————————————
全部入ってさえいればOK.
>git clone https://github.com/shawest/neucbot.git
Cloning into 'neucbot'...
remote: Counting objects: 338, done.
remote: Total 338 (delta 0), reused 0 (delta 0), pack-reused 338
Receiving objects: 100% (338/338), 233.25 KiB | 119.00 KiB/s, done.
Resolving deltas: 100% (121/121), done.
>cd neucbot
>ls
AlphaLists/ Chains/ Data/ License.txt
Materials/ README.md Scripts/ neucbot.py*
>./neucbot.py
No alpha list or chain specified
No target material specified
Usage: You must specify an alpha list or decay chain file and a target material file.
You may also specify a step size to for integrating the alphas as they slow down in MeV; the default value is 0.01 MeV
-l [alpha list file name]
-c [decay chain file name]
-m [material composition file name]
-s [alpha step size in MeV]
-t (to run TALYS for reactions not in libraries)
-d (download isotopic data for isotopes missing from database)
-o [output file name]
>./neucbot.py -l AlphaLists/Rn220Alphas.dat -m Materials/Acrylic.dat -o test1
load alpha list AlphaLists/Rn220Alphas.dat
load target material Materials/Acrylic.dat
Printing output to test1
Running alphas:
[============================================================] 100%
動くけど中身はからっぽ
-t and -d などの自動ダウンロードの オプションは環境によってはうまく動かない. 自分でダウンロードすれば解決. これもGitHub上に用意してある.
例:Fe data Download
>pwd
/hoge/hoge/neucbot
>cd Data/Isotopes
>git clone https://github.com/neucbot-datasets/fe.git
Cloning into 'fe'...
remote: Counting objects: 3, done.
remote: Compressing objects: 100% (2/2), done.
remote: Total 3 (delta 0), reused 3 (delta 0), pack-reused 0
Unpacking objects: 100% (3/3), done.
>tar xzf fe/Fe.tar.gz
>ls Fe
Fe54/ Fe56/ Fe57/ Fe58/
>cd /hoge/hoge/neucbot
>echo fe 0 100 > Materials/Iron.dat
>cat Materials/Iron.dat
fe 0 100
>./neucbot.py -c Chains/U238upperChain.dat -m Materials/Iron.dat -o iron_u.dat
load alpha chain Chains/U238upperChain.dat
load target material Materials/Iron.dat
Printing output to iron_u.dat
Running alphas:
[===========================================================]100%
Mizukoshi code download
>git clone https://github.com/mzks/neucbot_tutor.git
Cloning into 'neucbot_tutor'...
remote: Counting objects: 11, done.
remote: Compressing objects: 100% (8/8), done.
remote: Total 11 (delta 3), reused 6 (delta 2), pack-reused 0
Unpacking objects: 100% (11/11), done.
>root 'neucbot_tutor/show.C("iron_u.dat")'
root [0]
Processing neucbot_tutor/show.C("iron_u.dat")...
root [1]
この程度の入力で本当に計算できているのですか? 地下実験でよく来るα線はせいぜい数mmで 死ぬエネルギー領域 (数MeV) 確実に特定の物質中でα線が死ぬので, ジオメトリ,密度の情報は不要.
“NeuCBOT”が次の業界標準の ソフトウェアになる可能性が高い 正確,簡単,無料のNeuCBOTをつかおう! ナウでクールなヤングはGitHub+python.
Official repository https://github.com/shawest/neucbot see README.md Thesis Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment Volume 875, 11 December 2017, Pages 57-64 Radiogenic neutron yield calculations for low-background experiments - S. Westerdale *, P.D. Meyers
>cat Chains/U238Chain.dat U238 100
U234 100
Th230 100
Ra226 100
Rn222 100
Po218 99.98
At218 0.01998
Rn218 0.00002
Bi214 0.024
Po214 99.976
Pb210 0.0000019
Bi210 0.00013
Po210 99.99987
>cat Materials/RockKam.dat
h 0 0.6
c 0 0.3
o 0 60.7
na 0 3.9
mg 0 0.3
al 0 10.6
si 0 18.5
p 0 0.1
k 0 2.1
ca 0 1.8
ti 0 0.1
mn 0 0.0
fe 0 1.0