h = 10; r = 2;

pts = RandomPoint[Sphere[{0, 0, h}, r],
   1000, {{-r, r}, {-r, r}, {h, h + r}}];

Graphics3D[{Cylinder[{{0, 0, 0}, {0, 0, h}}, r], Point[pts]}]