//define the line
gp_Lin line(origin, dir);

//declare the tool
IntCurvesFace_ShapeIntersector tool;

//loop through the shapes. Shapes are either SOLID or COMPOUND shapes

for( int i = 0; i < numShapes; ++i) {
    
    TopoDS_Shape shape = shapes[i];
    
    tool.Load(shape, tolerance);
    tool.PerformNearest(line, 0, RealLast());
    
    if( !tool.IsDone() || tool.NbPnt() == 0) {
          continue;
     }

    //process results
}