import os
import sys
import numpy as np

dir = os.path.split(os.path.realpath(__file__))[0]
sys.path.append(dir + '/..')

from geometry import *
from solids import build_12inch_pmt, build_12inch_pmt_with_lc
from optics import *
from transform import rotate, make_rotation_matrix
from itertools import product
from make import cylinder

def build_lbne(radius, height, nstrings, pmts_per_string, endcap_spacing):
    pmt = build_12inch_pmt_with_lc()

    lbne = Geometry()

    # outer cylinder
    cylinder_mesh = cylinder(radius, radius, height+height/(pmts_per_string-1), theta=(2*np.pi/nstrings)/4)
    cylinder_mesh.vertices = rotate(cylinder_mesh.vertices, np.pi/2, (-1,0,0))
    lbne.add_solid(Solid(cylinder_mesh, water, vacuum, black_surface, 0xff))

    lbne.pmtids = []

    # construct the barrel
    for i in range(pmts_per_string):
        for j in range(nstrings):
            rotation = make_rotation_matrix(j*2*np.pi/nstrings, (0,0,1))
            displacement = rotate((0,-radius,-height/2+i*height/(pmts_per_string-1)), j*2*np.pi/nstrings, (0,0,1))
            lbne.pmtids.append(lbne.add_solid(pmt, rotation, displacement))

    # construct the top endcap
    for x, y in np.array(tuple(product(\
                np.arange(-radius+0.075, radius, endcap_spacing),
                np.arange(-radius+0.075, radius, endcap_spacing)))):
        if np.sqrt(x**2 + y**2) < radius - endcap_spacing/2:
            rotation = make_rotation_matrix(+np.pi/2, (1,0,0))
            displacement = (x,y,+height/2+height/(pmts_per_string-1)/2)
            lbne.pmtids.append(lbne.add_solid(pmt, rotation, displacement))

    # construct the bottom endcap
    for x, y in np.array(tuple(product(\
                np.arange(-radius+0.075, radius, endcap_spacing),
                np.arange(-radius+0.075, radius, endcap_spacing)))):
        if np.sqrt(x**2 + y**2) < radius - endcap_spacing/2:
            rotation = make_rotation_matrix(-np.pi/2, (1,0,0))
            displacement = (x,y,-height/2-height/(pmts_per_string-1)/2)
            lbne.pmtids.append(lbne.add_solid(pmt, rotation, displacement))

    return lbne