Change negative infinity to min double in function LogPrior.
Introduce coefficient factor to posterior density of flux | theta because when we compute estimate of theta it is very important.