docs: def find_collisions(Signal, tolerance=50):
      """
      Finds collision events in the signal from the shift in phase of the signal.
  
      Parameters
      ----------
      Signal : array_like
          Array containing the values of the signal of interest containing a single frequency.
      tolerance : float
          Percentage tolerance, if the value of the FM Discriminator varies from the mean by this
          percentage it is counted as being during a collision event (or the aftermath of an event).
  
      Returns
      -------
      Collisions : ndarray
          Array of booleans, true if during a collision event, false otherwise.
      """
      fmd = fm_discriminator(Signal)
      mean_fmd = _np.mean(fmd)
  
      Collisions = [_is_this_a_collision(
          [value, mean_fmd, tolerance]) for value in fmd]
  
      return Collisions