############################################################## 結果の合成(sumtally) ############################################################## * モンテカルロで別々に実行した計算結果は加算できる. ========================================================= 設定 ========================================================= * 実行時のファイルはそのまま流用すればよい. * 他の条件はすべて同じにする必要あり. * [parameters] セクションにて、icntrl = 13 を設定 * sumtallyしたいtally設定ファイルをコピー (sumtally_phits.inpなどへ) * sumtallyしたいtallyセクション中に以下を追記. :: sumtally start isumtally = 2 nfile = 2 sum/file_1.out 1.0 sum/file_2.out 1.0 sfile = sum/file_s.out sumfactor = 1.0 sumtally end * 通常通り、phitsを実行 ========================================================= 実行例 ========================================================= * 同じファイルを2個 1:1で読み込ませたとき. - 平均値に変化はない - 相対誤差が 0.707倍 (1/root(2) ) 倍になる. - 統計数が2倍 ⇒ 不確かさが減少.