select-points-in-polys.py

import arcpy


in_polys = "NcCountiesT.shp"
in_points = "NcNatRegPts.shp"

desc = arcpy.Describe(in_polys)
oid_fld = desc.OIDFieldName

# Make featurelayers for input into SelectByLocation
arcpy.MakeFeatureLayer_management(in_polys, "polys_lyr")
arcpy.MakeFeatureLayer_management(in_points, "points_lyr")

cursor = arcpy.SearchCursor("polys_lyr")
for row in cursor:
    # Get row id number
    row_id = row.getValue(oid_fld)
    # Select just the current polygon
    arcpy.SelectLayerByAttribute_management("polys_lyr", "NEW_SELECTION",
                                            '"{0}" = {1}'.format(oid_fld, row_id))
    # Run select by location on just the currently selected polygon
    pts = arcpy.SelectLayerByLocation_management("points_lyr", "intersect", "polys_lyr")
    matchcount = int(arcpy.GetCount_management('points_lyr').getOutput(0))
    print str(matchcount)