From a6b504dc3282f5775eb2ffaad42eff5545529f31 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 25 Sep 2018 10:24:59 -0500 Subject: update likelihood calculation to use PMT_TTS macro I accidentally hardcoded the single PE TTS to 1.5 ns in the likelihood calculation. --- src/likelihood.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/likelihood.c b/src/likelihood.c index 43fe336..81331f9 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -495,7 +495,7 @@ double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, doubl if (ev->pmt_hits[i].hit) { for (j = 1; j < MAX_PE; j++) { - logp[j] = log(pq(ev->pmt_hits[i].qhs,j)) - mu[i] + j*log_mu - lnfact(j) + log_pt(ev->pmt_hits[i].t, j, mu_noise, mu_indirect, &mu_direct[i], 1, &ts[i], tmean, 1.5); + logp[j] = log(pq(ev->pmt_hits[i].qhs,j)) - mu[i] + j*log_mu - lnfact(j) + log_pt(ev->pmt_hits[i].t, j, mu_noise, mu_indirect, &mu_direct[i], 1, &ts[i], tmean, PMT_TTS); if (j == 1 || logp[j] > max_logp) max_logp = logp[j]; if (logp[j] - max_logp < MIN_RATIO*ln(10)) { -- cgit