HDECAY 風波

最近因為要用到 Higgs 的 decay rate
從 Particle Data Group 的 review 一路找到這個程式 HDECAY
關於這個程式的使用方法可以在上 Arxiv
http://arxiv.org/abs/hep-ph/9704448

程式本身則可以上 Michael Spira 的個人網頁下載。

雖然說明文件表示這是一個使用簡單且快速的 fortran 程式，但我從下載到真正要可以使用上可是花了好大一番功夫。首先我不懂 fortran，連要用什麼編譯器都不知道。先是用在 window 裝了 g95，發現不知道 window 下要怎麼使用 makefile。之後捨棄改裝 MinGW，最後選擇了模擬 Linux 的 Cygwin 。雪上加霜的是對於基本指令相當不熟悉，另外開始時不知道 Cygwin 安裝時要點選需要的套件。一回兒 vim 沒有安裝，一回兒又不確定有沒有 g77，折騰了好多工夫。讓我不禁思考起投入 Linux 懷抱的可能性。

總算可以編譯 fortran 後，仍然跑不出預期的結果。懷疑了是否是 makefile 出問題、input 是否正確，找了許多幫手，最後甚至寄信給 Michael Spira。這位大叔很好心地幫我研究問題出在哪，最後發現是我的 input 檔案必須多加一行空白，程式才能判別文件結束。然而我還是沒跑出結果，但是 Michael 大叔說我的程式已經沒問題了呀。我只好回過頭來仔細閱讀說明文件，發現到預設 input 的參數並不會跑出結果，改過參數後才會開始 output 出數據。另外預設只會跑出一個數據，怪不得一瞬間程式就跑完了。無論如何，這起事件總算是落幕了。