assume the particle produces no light if it's outside the PSUP

1 files changed, 3 insertions, 3 deletions

diff --git a/src/likelihood.c b/src/likelihood.c
index b02aeef..90fdd82 100644
--- a/src/likelihood.c
+++ b/src/likelihood.c
@@ -1247,9 +1247,9 @@ double nll(event *ev, vertex *v, size_t n, double dx, int ns, const int fast, in
          * We calculate this here to make sure that we don't integrate outside
          * of the PSUP because otherwise there is no code to model the fact
          * that the PSUP is optically separated from the water outside. */
-        if (!intersect_sphere(v[j].pos,v[j].dir,PSUP_RADIUS,&distance_to_psup)) {
-            /* If the particle is outside the PSUP we just assume it produces
-             * no light. */
+        if (NORM(v[j].pos) > PSUP_RADIUS || !intersect_sphere(v[j].pos,v[j].dir,PSUP_RADIUS,&distance_to_psup)) {
+            /* If the particle is outside the PSUP or it doesn't intersect the
+             * PSUP we just assume it produces no light. */
             fprintf(stderr, "particle doesn't intersect PSUP!\n");
             continue;
         }