Initial commit

data/orca_footprint.dat

+1	-91.0	118.2
+2	-69.0	118.2
+3	-46.5	118.7
+4	-25.3	121.9
+5	-1.9	121.3
+6	18.8	118.6
+7	40.8	120.1
+8	66.3	118.2
+9	-79.4	100.4
+10	-58.7	100.1
+11	-35.1	101.0
+12	-14.2	101.3
+13	8.8	99.6
+14	30.9	100.5
+15	51.8	100.8
+16	73.4	98.8
+17	-90.9	76.4
+18	-67.0	78.7
+19	-46.8	78.4
+20	-24.0	79.6
+21	-3.9	77.8
+22	20.4	78.9
+23	41.8	78.7
+24	63.1	77.4
+25	83.1	78.5
+26	-102.8	56.4
+27	-79.6	59.5
+28	-56.8	56.1
+29	-35.5	57.0
+30	-15.4	56.6
+31	9.6	56.9
+32	30.5	58.4
+33	51.9	56.5
+34	74.8	57.2
+35	96.5	60.0
+36	-112.8	35.4
+37	-89.2	36.2
+38	-69.5	35.3
+39	-48.4	36.8
+40	-25.0	38.5
+41	-3.3	39.3
+42	19.6	37.7
+43	42.1	35.3
+44	63.2	37.9
+45	83.0	37.3
+46	105.7	38.3
+47	-100.8	15.8
+48	-79.9	14.2
+49	-59.8	15.7
+50	-35.9	15.4
+51	-14.8	16.8
+52	8.9	17.0
+53	30.8	16.1
+54	52.9	15.1
+55	75.0	15.2
+56	96.9	15.6
+57	117.6	16.1
+58	-113.4	-5.7
+59	-91.8	-6.5
+60	-70.7	-4.7
+61	-47.1	-4.3
+62	-25.5	-3.8
+63	-1.8	-6.2
+64	20.9	-5.0
+65	41.5	-2.5
+66	61.0	-3.5
+67	86.3	-5.7
+68	105.5	-3.5
+69	-101.0	-24.0
+70	-78.2	-25.9
+71	-56.2	-24.1
+72	-35.8	-23.8
+73	-14.5	-23.9
+74	7.5	-27.4
+75	28.0	-24.2
+76	50.7	-26.5
+77	72.2	-24.6
+78	97.0	-27.1
+79	117.8	-27.3
+80	-91.9	-46.4
+81	-70.1	-46.2
+82	-47.4	-47.4
+83	-24.3	-44.9
+84	-3.0	-46.4
+85	18.1	-48.0
+86	40.3	-47.9
+87	63.3	-48.1
+88	84.9	-44.6
+89	106.4	-48.2
+90	-103.1	-67.4
+91	-80.1	-67.4
+92	-58.0	-66.5
+93	-35.9	-68.3
+94	-12.9	-65.3
+95	8.4	-65.1
+96	30.7	-66.7
+97	52.8	-66.9
+98	75.9	-67.5
+99	96.7	-67.8
+100	-89.8	-89.3
+101	-67.3	-88.9
+102	-47.1	-88.3
+103	-23.3	-87.0
+104	-3.9	-88.7
+105	19.0	-86.5
+106	41.5	-87.4
+107	63.9	-85.9
+108	83.5	-87.4
+109	-80.0	-109.8
+110	-56.8	-108.1
+111	-35.3	-107.9
+112	-14.0	-106.2
+113	8.0	-106.8
+114	27.9	-107.9
+115	53.9	-108.2
+116	-45.0	-130.0
+117	-24.0	-128.7
+118	-1.0	-128.9
+119	21.2	-129.6
+120	41.0	-131.4
+

footprint.py

+#!/usr/bin/env python
+# coding=utf-8
+
+#from __future__ import division
+
+import math
+import numpy as np
+import matplotlib.pyplot as plt
+
+def load_footprint(filename):
+    i_xy = np.loadtxt(filename)
+    n = i_xy.shape[0]
+    fp = np.c_[ i_xy, np.zeros(n) ] 
+    return fp[:,[1,2,3]]
+
+def fluct(mu, sigma = 1.0):
+    return np.random.normal(mu, sigma)
+
+def radius(p):
+    return np.sqrt(p[:,0]**2 + p[:,1]**2)
+
+def subarea(footprint, rho):
+    return footprint[radius(footprint) < rho]
+
+def subroad(fp, rw):
+    return fp[(fp[:,1] < fp[:,0] + rw) & (fp[:,1] > fp[:,0] - rw)]
+
+def main():
+    fp0 = load_footprint("data/orca_footprint.dat")
+
+    fp1 = fluct(fp0)
+
+    fp3 = subarea(fp0, 30.0)
+
+    fp4 = subroad(fp0, 30.0)
+    
+    fig, ax = plt.subplots()
+
+    ax.set_title("KM3NeT/ORCA footprint")
+
+    ax.set_xlabel("x [m]")
+    
+    ax.set_ylabel("y [m]")
+
+    ax.scatter(fp0[:,0], fp0[:,1], marker='o')
+    #ax.scatter(fp1[:,0], fp1[:,1])
+    ax.scatter(fp3[:,0], fp3[:,1], marker='+')
+    ax.scatter(fp4[:,0], fp4[:,1], marker='x')
+    
+    plt.show()
+
+    
+if __name__ == '__main__':
+    main()