方法

• 目的の物体の中を通っているときのみ各ステップが終わった毎に　w = density×rxy²×steplength　 なるウエイトをつけて加える。rxyはstepの中心（出口や入口でもいいだろう）の座標をxy平面に投影したときの原点からの半径。 densityは通過物質の密度。steplengthは最後のstepの長さ。
• geantinoを(ηrange)×(φrange)の領域に一様にNray個飛ばしたあと total_weight=(ηrange)×(φrange)×Σw /Nray と計算すれば全重量が出る。

サンプルプログラム（実証のために作った（linux上 g++でcompileする））

• box_weight.ccと計算結果
• box_weight_c++.ccと計算結果（object orient化したら結果は同じで2.8倍速かった!!!）
• cylinder_weight.ccと計算結果

理由
geantinoが走る各ステップがカバーする小体積は図示すると右図のように計算できる。この体積ΔVに通過中の物質の密度ρをかければこのステップが カバーした物質の重量になる。 ここで重要なのはrapidityとθの微分関係

今後
現存のradiation length分布をとるプログラムを少し変更するだけで重量が求められるはず→ 中村浩のcodingと 結果。