eic-opticks/ana/sagitta.py

0001 #!/usr/bin/env python
0002 import numpy as np
0003 import matplotlib.pyplot as plt
0004
0005
0006 r = 17820
0007 c = plt.Circle((0, 0), r, fill=False, ec='red' )
0008
0009 a = 450
0010 s = r - np.sqrt( np.power(r,2) - np.power(a,2) )
0011
0012 b = plt.Rectangle( (-a,r-s), 2*a, 2*a, fill=False, ec='blue' )
0013
0014 bbox = plt.Rectangle( (-r,-r),    2*r, 2*r+2*a-s , fill=False, ec='blue' )
0015
0016
0017 fig, ax = plt.subplots()
0018
0019 ax.set_xlim( -r*1.1, r*1.1 )
0020 ax.set_ylim( -r*1.1, r*1.1 )
0021 ax.set_aspect('equal')
0022
0023
0024 ax.add_patch(b)
0025 ax.add_patch(c)
0026 ax.add_patch(bbox)
0027
0028 ax.plot( [0, 0], [0, r*1.1], linestyle="dashed" )
0029
0030
0031 fig.show()
0032
0033
0034
0035 #fig.savefig('/tmp/sagitta.png')
0036