Generate random unit vector in d dimensional space in Python

I was wondering what is the best way to:

Generate t = 160 unit vectors in Rd for d = 100.

Plot of cdf of their pairwise dot products aka calculate nC2 dot products.