How to get face normal's x or y axis vector (or matrix)?
obj = bpy.context.active_object
obj.data.calc_tangents()

print(obj.data.loops[obj.data.polygons[0].loop_indices[0]].tangent)
print(obj.data.loops[obj.data.polygons[0].loop_indices[0]].bitangent)

I could find a solution by myself.
Thanks!