diff --git a/src/black_holes/EAGLE/black_holes.h b/src/black_holes/EAGLE/black_holes.h
index 0cb18d6ed181f0458d9ad59ea87c22b28e3875ed..5311b1bb85f7dbc3128c1621d759d86f4b62a844 100644
--- a/src/black_holes/EAGLE/black_holes.h
+++ b/src/black_holes/EAGLE/black_holes.h
@@ -30,6 +30,7 @@
 #include "minmax.h"
 #include "physical_constants.h"
 #include "random.h"
+#include "rays.h"
 
 /* Standard includes */
 #include <float.h>
@@ -180,6 +181,9 @@ __attribute__((always_inline)) INLINE static void black_holes_init_bpart(
   bp->f_visc = FLT_MAX;
   bp->accretion_boost_factor = -FLT_MAX;
   bp->mass_at_start_of_step = bp->mass; /* bp->mass may grow in nibbling mode */
+
+  /* Reset the rays carried by this BH */
+  ray_init(bp->rays, eagle_blackhole_number_of_rays);
 }
 
 /**