Ok great it worked!
I place the code here for future reference.
One thing that i would do is to search for symmetry along x-axis (i.e if y-coordinate of srf = 0 ) and split it because right now the outside of boundary symmetric surfaces are not removed.
import rhinoscriptsyntax as rs
filter = rs.filter.surface + rs.filter.polysurface
breps = rs.GetObjects("Select brep to split", filter)
cutter = rs.GetObjects("Select cutting brep", filter)
for brep in breps:
splits = rs.SplitBrep ( brep, cutter )
for obj in splits:
cog = rs.SurfaceAreaCentroid(obj)
if not rs.IsPointInSurface(cutter,cog):
Thank you all for the help.