as opposed to Abe's technical appendix, we do draw tau also for case z=1; hence mu can be drawn according to standard update formula, as gamma is conjugate prior of exponential distribution
avoid numeric overflow