diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 10:24:59 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 10:24:59 -0500 |
| commit | a6b504dc3282f5775eb2ffaad42eff5545529f31 (patch) | |
| tree | 70cff0ba7a9b75286d7a17b97641c645b0af894d | |
| parent | b4be76539cbc0254df33e8080f965df3a4ce92a4 (diff) | |
| download | sddm-a6b504dc3282f5775eb2ffaad42eff5545529f31.tar.gz sddm-a6b504dc3282f5775eb2ffaad42eff5545529f31.tar.bz2 sddm-a6b504dc3282f5775eb2ffaad42eff5545529f31.zip | |
update likelihood calculation to use PMT_TTS macro
I accidentally hardcoded the single PE TTS to 1.5 ns in the likelihood
calculation.
| -rw-r--r-- | src/likelihood.c | 2 |
1 files changed, 1 insertions, 1 deletions
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)) { |