diff options
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index 97db51c..0855d08 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -266,10 +266,12 @@ double F(double t, double mu_noise, double mu_indirect, double *mu_direct, doubl mu_total = mu_noise + mu_indirect; for (i = 0; i < n; i++) { + if (mu_direct[i] == 0.0) continue; p += mu_direct[i]*norm_cdf(t,ts[i],sigma); mu_total += mu_direct[i]; } for (i = 0; i < n; i++) { + if (mu_shower[i] == 0.0) continue; p += mu_shower[i]*norm_cdf(t,ts_shower[i],ts_sigma[i]); mu_total += mu_shower[i]; } @@ -300,10 +302,12 @@ double f(double t, double mu_noise, double mu_indirect, double *mu_direct, doubl mu_total = mu_noise + mu_indirect; for (i = 0; i < n; i++) { + if (mu_direct[i] == 0.0) continue; p += mu_direct[i]*norm(t,ts[i],sigma); mu_total += mu_direct[i]; } for (i = 0; i < n; i++) { + if (mu_shower[i] == 0.0) continue; p += mu_shower[i]*norm(t,ts_shower[i],ts_sigma[i]); mu_total += mu_shower[i]; } |