아침에 무턱대고 연구를 시작하는 것보다 1시간 성경을 읽고 시작하는 것이 훨씬 효과적이다. 성경을 읽고나면 반드시 정리된 생각과 새로운 아이디어를 주신다. 어제 아침에 사사기를 읽고나서 최근 문제가 되는 stylized model에서의 장기부피 계산방법이 떠올랐다.

보통은 이 작업을 위해서 Monte Carlo 적분 프로그램을 작성하게 되는데 여간 번거로운 작업이 아닐 수 없다. 특히 이곳에서 최근 업데이트한 모델의 경우 머리 속에만도 침샘, ET region등이 복잡하게 들어있어 equation을 풀어서 적분을 하자니 이 작업을 한번도 해보지 못한 나로서는 엄두가 나지 않았다.

하나님께서 주신 아이디어는 MCNP안에 분명히 부피계산기능이 들어있을꺼라는 생각이었다. 특정 장기를 source로 잡을 경우 일정 육면체나 구, 원통형으로 둘러싸게 되는데 이때 random하게 정해준 볼륨 안에서 starting point를 sampling한 뒤 그 point가 source로 정해준 cell에 들어있는지를 check하게 되고, 이 때 sampling efficiency가 떨어지면 에러가 난다. 이전에 이런 error를 접한 적이 많았고 그 때를 떠올리며 그렇다면 sampling efficiency가 바로 Monte Carlo적분을 말한다는 사실을 깨달았다.

특정 장기를 육면체로 둘러싸고 충분한 history를 줘서 source sampling을 한뒤 output을 보면 source efficiency가 나오는데 이것을 육면체의 부피에 곱하면 장기의 부피가 나온다. 육면체를 정할 때는 geometry plot을 해서 xy평면과 xz평면에 그려보고 range를 정하면 된다. 이 방법을 통해 최근 revised된 stylized model에서 잘못된 부피를 찾아낼 수 있었다.