diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 10:26:23 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 10:26:23 -0500 |
| commit | f23450aa8b225fc38a3fdbf9d9db489f4aba7158 (patch) | |
| tree | f7d7db8c224e48ef79e6085e1f00697aa5fba078 /src | |
| parent | a6b504dc3282f5775eb2ffaad42eff5545529f31 (diff) | |
| download | sddm-f23450aa8b225fc38a3fdbf9d9db489f4aba7158.tar.gz sddm-f23450aa8b225fc38a3fdbf9d9db489f4aba7158.tar.bz2 sddm-f23450aa8b225fc38a3fdbf9d9db489f4aba7158.zip | |
update indirect scattering PDF start time
Currently the PDF for scattered light is modelled as a flat distribution
starting at some time t. Previously I was using the mean hit time for all PMTs,
however this should really be a flat distribution in the time *residual* after
the main peak. Therefore, the PDF now starts at the estimated time for direct
photons.
Diffstat (limited to 'src')
| -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 81331f9..ef6622d 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, PMT_TTS); + 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], ts[i], PMT_TTS); if (j == 1 || logp[j] > max_logp) max_logp = logp[j]; if (logp[j] - max_logp < MIN_RATIO*ln(10)) { |