diff options
| -rw-r--r-- | src/likelihood.c | 6 | ||||
| -rw-r--r-- | src/pmt.h | 5 |
2 files changed, 8 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]); @@ -1,6 +1,11 @@ #ifndef PMT_H #define PMT_H +/* FIXME: The OWL PMTs normal vector isn't defined in pmt.txt! + * + * I'm not sure why, but for now we shouldn't use the OWL tubes in the + * likelihood calculation. */ + #define MAX_PMTS 10000 /* Radius of the PMT concentrator (cm). */ #define PMT_RADIUS 13.5 |