This example runs well. Two concerns to me:
1) The reconstructed maps look noisy and pretty bad with respect to what Wu2008 shows. Do we have better data then that to show?

2) Both the rtop signal and pdf have a colorbar with mininum value below zero. Are there negative rtop values reconstructed? This is impossible, right? RTOP is the sum of the signal... 

Please check these two issues.