Multi-Purpose Particle and Heavy Ion Transport code System

Slides:



Advertisements
Similar presentations
P HI T S 陽子ビームで雪だるまを溶かそう Multi-Purpose Particle and Heavy Ion Transport code System title 年 3 月改訂.
Advertisements

P HI T S α 線, β 線, γ 線,中性子線を止めるに は? Multi-Purpose Particle and Heavy Ion Transport code System title 年 5 月改訂.
P HI T S PSFC4PHITS の使い方 Multi-Purpose Particle and Heavy Ion Transport code System Title 年 2 月改訂.
P HI T S スクリプト言語を用いた PHITS の連続 実行 Multi-Purpose Particle and Heavy Ion Transport code System title 年 2 月改訂.
第5章 JMPのインストールと基本操作 廣野元久
初年次セミナー 第13回 2次元グラフィックス(1).
ボクセルファントムを用いた X線治療のシミュレーション
情報・知能工学系 山本一公 プログラミング演習Ⅱ 第3回 配列(1) 情報・知能工学系 山本一公
PHITS 講習会 基礎実習(III): 計算条件の設定
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
PHITS講習会 基礎実習(I): 体系及び線源の定義
PHITS講習会 基礎実習(II): Tally(タリー)の定義
強度変調回転放射線治療(VMAT)の シミュレーション
Multi-Purpose Particle and Heavy Ion Transport code System
REIMEI EISA Viewerの使い方
実習用サーバの利用開始手順 (Windowsの場合) TeraTerm Proをインストール 公開鍵をメールで送付
Multi-Purpose Particle and Heavy Ion Transport code System
DICOM医療画像を使用したPHITSシミュレーション
実習問題の解き方 JIS3-2の場合.
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
情報基礎A 第13週 VBAプログラミング VBAの基本文法7・実際のデータ処理
PHITSによるX線治療シミュレーション 基本操作復習編
第13回授業(7/10-1) での学習内容 6月19日に宿題にした平均値の差の検定結果、及び7月3日の授業で実習した同検定結果のウェブ上での検算のやり方を学習する。 この検算の宿題は、春期定期試験の時に、今日渡す2枚目の出欠表の裏に授業中の手計算による結果と比較して、手計算の結果が正しかったかを報告する。
ParaViewを用いたPHITS 計算結果の3次元表示
EGSに対応した粒子軌跡と 計算体系の3次元表示ソフト - CGVIEW -
第5回授業(5/08) の目標 第1章のウェブ上での宿題のやり方の続き(ここまでで、終了)。
軌跡とジオメトリー表示プログラム CGVIEW(Ver2.2)の改良
Outlook で送信したメールの 添付ファイルが消える
情報コミュニケーション入門 2016年度版 情報コミュニケーション入門m 総合実習(1).
スクリプト言語を用いたPHITSの連続実行
PHITS 講習会 基礎実習(III): 計算条件の設定
Multi-Purpose Particle and Heavy Ion Transport code System
誘導放射能の評価方法: PHITSとDCHAIN-SPの接続計算
α線,β線,γ線,中性子線を止めるには?
情報工学科 3年生対象 専門科目 システムプログラミング 第5回、第6回 ヒアドキュメント レポート課題 情報工学科 篠埜 功.
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
平山 英夫、波戸 芳仁 KEK, 高エネルギー加速器研究機構
Multi-Purpose Particle and Heavy Ion Transport code System
応用実習用資料 Environmental radioactivity
Gmsh, Elmer, Garfield++のInstall
Multi-Purpose Particle and Heavy Ion Transport code System
MeV internal meeting Oct. 2, 2015
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 A
PHITS 講習会 基礎実習(III): 計算条件の設定
応用実習用資料 Neutron target
IAEA phase space fileを用いた X線治療シミュレーション
単語登録(1) ◎MS-IMEの「単語登録」に、単語、語句、記号など自分がよく使うものを登録しておくと、便利である。
実習問題の解き方 JIS3-3の場合.
Processing使用の準備.
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
EGSに対応した粒子軌跡と 計算体系の3次元表示ソフト - CGVIEW -
プログラミング言語Ⅰ(実習を含む。), 計算機言語Ⅰ・計算機言語演習Ⅰ, 情報処理言語Ⅰ(実習を含む。)
ファイルのアップロード HTMLファイルをWebサーバにアップロード 名商大のWebサーバ(opinion.nucba.ac.jp)
コンピュータ プレゼンテーション.
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
Multi-Purpose Particle and Heavy Ion Transport code System
プログラミング入門2 第13回、14回 総合演習 情報工学科 篠埜 功.
α線,β線,γ線,中性子線を止めるには?
Multi-Purpose Particle and Heavy Ion Transport code System
平山 英夫、波戸 芳仁 KEK, 高エネルギー加速器研究機構
BLJ2013 BentleyArchitecture
計算効率を上げるためのvariance reduction(粒子のウェイト)の利用 B
目的とするユーザーコードを 作成するために
単語登録(1) ◎MS-IMEの「単語登録」に、単語、語句、記号など自分がよく使うものを登録しておくと、便利である。
Presentation transcript:

Multi-Purpose Particle and Heavy Ion Transport code System PHITS Multi-Purpose Particle and Heavy Ion Transport code System PSFC4PHITSの使い方 PHITS講習会 入門実習 2016年5月改訂 Title 1

Phase Space Fileとは 医療用加速器や60Co線源により生成されたビーム情報をまとめたデータファイル What is PSF メーカーが公開していない場合が多い加速器部分(上流)のシミュレーションを省略し、ビーム照射部分(下流)のみのシミュレーションが可能 IAEAのNuclear Data Servicesから入手可能 IAEA形式により作成されたヘッダーファイルとバイナリデータファイルの組 https://www-nds.iaea.org/phsp/phsp.htmlx What is PSF 2

Phase Space Fileとは 医療用加速器や60Co線源により生成されたビーム情報をまとめたデータファイル What is PSF メーカーが公開していない場合が多い加速器部分(上流)のシミュレーションを省略し、ビーム照射部分(下流)のみのシミュレーションが可能 IAEAのNuclear Data Servicesから入手可能 IAEA形式により作成されたヘッダーファイルとバイナリデータファイルの組 本資料では下のURLからダウンロードできる Varian_Clinac_600C_6MV_1x1 を例にPSFC4PHITSの使用方法についてご紹介します ヘッダーファイルURL: https://www-nds.iaea.org/phsp/photon1/Varian_Clinac_600C_6MV_1x1.IAEAheader バイナリデータファイルURL: https://www-nds.iaea.org/phsp/photon/Varian_Clinac_600C_6MV/Varian_Clinac_600C_6MV_1x1.IAEAphsp What is PSF 3

Phase Space Fileとは What is PSF 4 ヘッダーファイル(Varian_Clinac_600C_6MV_1x1.IAEAheader) テキストエディタで開くことができます $IAEA_INDEX: 600 // IAEA website: www-nds.iaea.org/phsp/photon/ $TITLE: PHASESPACE in IAEA format for Varian Clinac 600C, 6MV photon ・ ・ ・ ・ ・ ・ $BYTE_ORDER: 1234 $ORIG_HISTORIES: 69184770 $PARTICLES: 9414373 規格化のために必要となるオリジナル計算のヒストリー数($ORIG_HISTORIES)と粒子データ数($PARTICLES)の情報も書かれています。 What is PSF 4

変換プログラム(psfc4phits) バイナリデータであるphase space fileをPHITS形式のdump dataに変換するプログラム IAEAにより提供されている変換プログラム*をPHITS用に改変*https://www-nds.iaea.org/reports-new/indc-reports/indc-nds/indc-nds-0484.pdfを参照 ①入力ファイルを作成(psfc4phits.inp)⇒詳しくはREADME参照 ②入力ファイルがあるフォルダにphase space fileのheaderファイルとデータファイルをおく ③実行 (Windows) psfc4phits_win.batにpsfc4phits.inpをドラッグ&ドロップ    (Mac) psfc4phits_mac.commnadをダブルクリック、現れる窓にpsfc4phits.inpと入力 Phase space fileの名前(拡張子を除いた部分) Varian_Clinac_600C_6MV_1x1 $ Phase space filename dmp-PHITS.out $ Output filename 1 $ Starting line 100000 $ Ending line -1 $ iopt: -1 (ASCII), 1 (binary[unformatted]) PHITS用dump data file名 変換する範囲の指定(変換できる数値はheaderファイルにある$PARTICLESの値が上限) iopt: 変換後のデータ形式 PSFC4PHITS HowTo 5

PSFC4PHITS実行結果 入力ファイルで指定したPHITS用dump data fileが生成されます dmp-PHITS.out (iopt=-1の場合) 左から粒子のkf-code, x,y,z座標, 速度ベクトルのx,y,z成分, エネルギー, ウエイト値の9個がアスキー形式で出力されます PSFC4PHITS outputs 6

setting to use dump data PHITSインプットファイル 変換したPHITS形式のdump dataを線源として利用する場合は下のように[source]セクションを設定します PHITSの計算結果をオリジナルの線源粒子あたりの量にする場合は規格化が必要 phits.inp [ S o u r c e ] set:c1[69184770] set:c2[9414373] totfact = c2/c1 s-type = 17 file = dmp-PHITS.out dump = -9 1 2 3 4 5 6 7 8 9 C1: Phase space file作成時のオリジナルのヒストリー数(headerファイルにある$ORIG_HISTORIES) C2: Phase space fileに含まれる粒子データ数(headerファイルにある$PARTICLESを確認) Dump定義文 9個のデータが出力されるので、アスキー形式の場合(iopt=-1)はdump=-9, バイナリの場合(iopt=1)はdump=9 次行には出力されるデータのID番号を順番通りに並べる setting to use dump data 7

Result of t-track tally PHITSの計算結果 track_xz.eps Z=66.8cmの面からz軸の正の方向に光子ビームが飛んでいる Result of t-track tally 8

注意点 本変換プログラムはIAEAのINDC (International Nuclear Data Committee)のもと Roberto Capote Noy氏(IAEA) Iwan Kawrakow氏(National Research Council of Canada Ottawa) により開発されたプログラムをPHITS用に改変したものです。 本プログラムでPHITS形式に変換できるのはIAEA形式で作成されたデータファイルのみです。 変換したデータを用いてPHITSを実行する際、総ヒストリー数(maxcas * maxbch)を変換したデータの個数より多く設定した場合、重複してデータが使用されることになります。この場合、統計的に偏った結果を与えることになります。 Important notice 9