diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-09-20 14:35:26 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-09-20 14:35:26 -0500 |
| commit | 433d6f3a23bb2a987428a125683dc01c20062a42 (patch) | |
| tree | 33929d4b43d81963c56677b739cdcbe75abf896d /src/likelihood.c | |
| parent | 378ecbc7af093382ed6b4faa4ceb645de1f00f1d (diff) | |
| download | sddm-433d6f3a23bb2a987428a125683dc01c20062a42.tar.gz sddm-433d6f3a23bb2a987428a125683dc01c20062a42.tar.bz2 sddm-433d6f3a23bb2a987428a125683dc01c20062a42.zip | |
don't include the OWL PMTs in the likelihood calculation
For some reason the OWL tubes have 9999.00 for the x, y, and z coordinates of
the normal vector in the PMT file. For now, I'm just going to remove them from
the likelihood calculation.
Diffstat (limited to 'src/likelihood.c')
| -rw-r--r-- | src/likelihood.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/likelihood.c b/src/likelihood.c index 33bf095..db59040 100644 --- a/src/likelihood.c +++ b/src/likelihood.c @@ -368,7 +368,7 @@ double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, doubl total_charge = 0.0; for (i = 0; i < MAX_PMTS; i++) { - if (ev->pmt_hits[i].flags || (pmts[i].pmt_type != PMT_NORMAL && pmts[i].pmt_type != PMT_OWL)) continue; + if (ev->pmt_hits[i].flags || pmts[i].pmt_type != PMT_NORMAL) continue; params.i = i; @@ -406,13 +406,13 @@ double nll_muon(event *ev, double T0, double *pos, double *dir, double t0, doubl mu_indirect = total_charge*CHARGE_FRACTION/10000.0; for (i = 0; i < MAX_PMTS; i++) { - if (ev->pmt_hits[i].flags || (pmts[i].pmt_type != PMT_NORMAL && pmts[i].pmt_type != PMT_OWL)) continue; + if (ev->pmt_hits[i].flags || pmts[i].pmt_type != PMT_NORMAL) continue; mu[i] = mu_direct[i] + mu_indirect + mu_noise; } nhit = 0; for (i = 0; i < MAX_PMTS; i++) { - if (ev->pmt_hits[i].flags || (pmts[i].pmt_type != PMT_NORMAL && pmts[i].pmt_type != PMT_OWL)) continue; + if (ev->pmt_hits[i].flags || pmts[i].pmt_type != PMT_NORMAL) continue; log_mu = log(mu[i]); |