import numpy as np
import pylab as pl
import math
import random

num_events = 1000

# this uses a forloop...
'''
xrad = np.array([random.uniform(-np.pi,np.pi) for z in range(num_events)])
y = np.array([random.uniform(-1,1) for z in range(num_events)])
yrad = np.arcsin(y)
'''
xrad = np.array([np.random.uniform(-np.pi,np.pi,size=num_events)])
y = np.array([np.random.uniform(-1,1,size=num_events)])
yrad = np.arcsin(y)


fig2 = pl.figure(2)
ax1 = fig2.add_subplot(111, projection="mollweide")

ax1.scatter(xrad,yrad)

ax1.plot([-math.pi, math.pi], [0,0],'r-')
ax1.plot([0,0],[-math.pi, math.pi], 'r-')


pl.title("Random Black Hole Merger events")

pl.grid(True)


pl.show()
