時間カット、ウエイトカット、ウエイトウインドウ -------------------------------------------------------------------------------- .. rst-class:: no-caption-number .. list-table:: **tmax(i)** :header-rows: 0 * - 値 - 説明 * - | (D=1.e+9) - | 粒子番号i(i=1-20)の粒子の計算打切時間[nsec]。粒子番号は :numref:`tbl-partt` を参照。 .. rst-class:: no-caption-number .. list-table:: **wc1(i)** :header-rows: 0 * - 値 - 説明 * - | (D=-0.5) - | i-th 粒子のminimum weight .. rst-class:: no-caption-number .. list-table:: **wc2(i)** :header-rows: 0 * - 値 - 説明 * - | (D=wc1/2) - | i-th 粒子のcutoff weight .. rst-class:: no-caption-number .. list-table:: **swtm(i)** :header-rows: 0 * - 値 - 説明 * - | (D=1.0) - | i-th 粒子のminimum source weight .. rst-class:: no-caption-number .. list-table:: **wupn** :header-rows: 0 * - 値 - 説明 * - | (D=5) - | ウェイトウィンドウの上限値 | = [Weight Window]セクションで指定した下限値 :math:`\times` wupn. (wupn :math:`\ge` 2) .. rst-class:: no-caption-number .. list-table:: **wsurvn** :header-rows: 0 * - 値 - 説明 * - | (0.6\*wupn) - | サバイバルウエイト値(1 :math:`<` wsurvn :math:`<` wupn) .. rst-class:: no-caption-number .. list-table:: **mxspln** :header-rows: 0 * - 値 - 説明 * - | (D=5) - | スプリット数の最大値、サバイバルの最大倍数(mxspln :math:`>` 1) .. rst-class:: no-caption-number .. list-table:: **mwhere** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | ウエイトウインドウのアクションの場所 | -1: 核反応時、0: 両方、1: 境界横断時 .. rst-class:: no-caption-number .. list-table:: **iwwbias** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | [ww bias]に関するオプション。 * - | =0 - | [ww bias]の値を用いない。 * - | =1 - | [ww bias]で与えたバイアスの値の逆数が[weight window]で指定したウエイトの下限値にかかる。 .. rst-class:: no-caption-number .. list-table:: **nfcseg** :header-rows: 0 * - 値 - 説明 * - | (D=50) - | 荷電粒子に対してforced collisionを適用する際の対象領域分割数。この値を小さくすると計算時間が短くなりますが、厚いターゲットの場合は、ターゲット内での詳細なエネルギー変化を考慮できなくなるため計算精度が落ちる可能性があります。 .. rst-class:: no-caption-number .. list-table:: **istdcut** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | 個別タリー打ち切り機能の切り替え。 * - | =0 - | 統計誤差がstdcut以下になっても計算終了時まで全てのタリー計算を続行します。なお,istdcutでタリーを打ち切った後に再開始計算を行うと,乱数の重複が起きてしまう可能性があります。再開始計算を行う場合は,前回のシミュレーションで最後まで実行されたタリーをインプットファイルの最初に書くようにしてください。 * - | =1 - | 統計誤差がstdcut以下になったタリーから計算を打ち切ります。 .. rst-class:: no-caption-number .. list-table:: **istdbat** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | 計算初期にstdcutを無効にするパラメータ。バッチ番号がistdbat以下の場合,stdcutによる計算打ち切りを行いません。OpenMPを利用した場合,バッチ分散から統計誤差を推定するため(istdev = 1),計算初期に統計誤差が不安定になります。そのような場合は,誤差が安定してくるバッチ数をistdbatに指定してください。 .. rst-class:: no-caption-number .. list-table:: **ctmin(1-3)** :header-rows: 0 * - 値 - 説明 * - (省略可, D=-9999) - カウンター1-3の下限値。各カウンター値がこの値以下になると粒子輸送が打ち切られます。粒子毎に異なる値は設定できません。 .. rst-class:: no-caption-number .. list-table:: **ctmax(1-3)** :header-rows: 0 * - 値 - 説明 * - (省略可, D=9999) - カウンター1-3の上限値。各カウンター値がこの値以上になると粒子輸送が打ち切られます。粒子毎に異なる値は設定できません。 時間のカットオフは粒子毎に tmax(i)= で指定します。単位はnsecです。カットオフ時間を上回った粒子は輸送を終了し、殺されます。高エネルギーの場合ほとんど意味がありませんが、低エネルギーの輸送で利用できます。 ウエイトカットオフは、インポータンス、強制衝突、implicit capture、ウエイトウインドウを利用し、粒子のウエイトが変化する場合に、指定のウエイトカットオフ値になった粒子にロシアンルーレットを施すものです。ウエイトウインドウを指定した粒子については作用しません。 粒子のウエイトが、 :math:`WC2` とソースの発生場所のインポータンスと現在のインポータンスの比 :math:`R` の積、 :math:`WC2 \times R` より小さくなった時、現在の粒子ウエイト :math:`WGT` の関数、 :math:`WGT/(WC1 \times R)` の確率で存続させ、その時粒子のウエイトを :math:`WGT = WC1 \times R` とします。それ以外は、その粒子を殺します。 :math:`WC1, WC2` の入力が負の数で与えられた場合は、 :math:`|WC1| \times SWTM` また、 :math:`|WC2| \times SWTM` が :math:`WC1, WC2` として設定されます。インポータンスが指定されていない粒子や、領域は、インポータンスは全て1にセットされます。 ウエイトウインドウのパラメータについては、通常はデフォルト値で省略できます。 iwwbias=1とした場合、input echoの[weight window]には[ww bias]で定義した値を考慮した結果が書き出され、[ww bias]にはoffが付きます。もし、[ww bias]が定義されていない場合は、input echoにおいて全ての値が1となった[ww bias]が書き出されます。