On the numerical characterization of the reachability cone for an essentially nonnegative matrix